Source author record

Jose H. Groh

Jose H. Groh appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2022arXiv

Grids of stellar models with rotation VII: Models from 0.8 to 300 M$_\odot$ at super-solar metallicity (Z = 0.020)

We present a grid of stellar models at super-solar metallicity (Z = 0.020) extending the previous grids of Geneva models at solar and sub-solar metallicities. A metallicity of Z = 0.020 was chosen to match that of the inner Galactic disk. A modest increase of 43% (=0.02/0.014) in metallicity compared to solar models means that the models evolve similarly to solar models but with slightly larger mass loss. Mass loss limits the final total masses of the super-solar models to 35 M$_\odot$ even for stars with initial masses much larger than 100 M$_\odot$. Mass loss is strong enough in stars above 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) to remove the entire hydrogen-rich envelope. Our models thus predict SNII below 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) and SNIb (possibly SNIc) above that. We computed both isochrones and synthetic clusters to compare our super-solar models to the Westerlund 1 (Wd1) massive young cluster. A synthetic cluster combining rotating and non-rotating models with an age spread between log10 (age/yr) = 6.7 and 7.0 is able to reproduce qualitatively the observed populations of WR, RSG and YSG stars in Wd1, in particular their simultaneous presence at log10(L/L$_\odot$) = 5-5.5. The quantitative agreement is imperfect and we discuss the likely causes: synthetic cluster parameters, binary interactions, mass loss and their related uncertainties. In particular, mass loss in the cool part of the HRD plays a key role.

preprint2021arXiv

Gaia EDR3 confirms that Westerlund 1 is closer and older than previously thought

Using {\it Gaia} Early Data Release 3 (EDR3) parallaxes and Bayesian inference, we infer a parallax of the Westerlund 1 (Wd1) cluster. We find a parallax of $0.34\pm{0.05}$ mas corresponding to a distance of $2.8^{+0.7}_{-0.6}$ kpc. The new {\it Gaia} EDR3 distance is consistent with our previous result using {\it Gaia} DR2 parallaxes. This confirms that Wd1 is less massive and older than previously assumed. Compared to DR2, the EDR3 individual parallax uncertainties for each star decreased by 30\%. However, the aggregate parallax uncertainty for the cluster remained the same. This suggests that the uncertainty is dominated by systematics, which is possibly due to crowding, motions within the cluster, or motions due to binary orbits.

preprint2021arXiv

Is GW190521 the merger of black holes from the first stellar generations?

GW190521 challenges our understanding of the late-stage evolution of massive stars and the effects of the pair-instability in particular. We discuss the possibility that stars at low or zero metallicity could retain most of their hydrogen envelope until the pre-supernova stage, avoid the pulsational pair-instability regime and produce a black hole with a mass in the mass gap by fallback. We present a series of new stellar evolution models at zero and low metallicity computed with the Geneva and MESA stellar evolution codes and compare to existing grids of models. Models with a metallicity in the range 0-0.0004 have three properties which favour higher BH masses as compared to higher metallicity models. These are (i) lower mass-loss rates during the post-MS phase, (ii) a more compact star disfavouring binary interaction and (iii) possible H-He shell interactions which lower the CO core mass. We conclude that it is possible that GW190521 may be the merger of black holes produced directly by massive stars from the first stellar generations. Our models indicate BH masses up to 70-75 Msun. Uncertainties related to convective mixing, mass loss, H-He shell interactions and pair-instability pulsations may increase this limit to ~85 Msun.

preprint2020arXiv

CO, Water, and Possible Methanol in Eta Carinae Approaching Periastron

In circumstellar gas, the complex organic molecule methanol has been found almost exclusively around young stellar objects, and is thus regarded as a signpost of recent star formation. Here we report the first probable detection of methanol around an evolved high-mass star, in the complex circumstellar environment around the Luminous Blue Variable $η$ Carinae, while using ALMA to investigate molecular cloud conditions traced by CO (2-1) in an orbit phase of the massive binary preceding the 2020 periastron. Favoring methanol over a $^{13}$CS alternative, the emission originates from hot ($T_{\rm{gas}} \simeq$ 700 K) material, $\sim$2$''$ (0.02 pc) across, centered on the dust-obscured binary in contrast to the CO which traces inner layers of the extended massive equatorial torus, and is accompanied by prominent absorption in a cooler ($T_{\rm{gas}} \simeq$ 110 K) layer of gas. We also report detections of water in $Herschel$/HIFI observations at 557 GHz and 988 GHz. The methanol abundance is several to 50 times higher than observed towards several lower mass stars, while water abundances are similar to those observed in cool, dense molecular clouds. The very high methanol:water abundance ratio in the core of $η$ Carinae may suggest methanol formation processes similar to Fischer-Tropsch-type catalytic reactions on dust grains. These observations prove that complex molecule formation can occur in the chemically evolved environments around massive stars in the end stages of their evolution, given sufficient gas densities and shielding conditions as may occur in material around massive interacting companions and merger remnants.

preprint2020arXiv

Inferring the parallax of Westerlund 1 from Gaia DR2

Westerlund 1 (Wd1) is potentially the largest star cluster in the Galaxy. That designation critically depends upon the distance to the cluster, yet the cluster is highly obscured, making luminosity-based distance estimates difficult. Using {\it Gaia} Data Release 2 (DR2) parallaxes and Bayesian inference, we infer a parallax of $0.35^{+0.07}_{-0.06}$ mas corresponding to a distance of $2.6^{+0.6}_{-0.4}$ kpc. To leverage the combined statistics of all stars in the direction of Wd1, we derive the Bayesian model for a cluster of stars hidden among Galactic field stars; this model includes the parallax zero-point. Previous estimates for the distance to Wd1 ranged from 1.0 to 5.5 kpc, although values around 5 kpc have usually been adopted. The {\it Gaia} DR2 parallaxes reduce the uncertainty from a factor of 3 to 18\% and rules out the most often quoted value of 5 kpc with 99\% confidence. This new distance allows for more accurate mass and age determinations for the stars in Wd1. For example, the previously inferred initial mass at the main-sequence turn-off was around 40 M$_{\odot}$; the new {\it Gaia} DR2 distance shifts this down to about 22 M$_{\odot}$. This has important implications for our understanding of the late stages of stellar evolution, including the initial mass of the magnetar and the LBV in Wd1. Similarly, the new distance suggests that the total cluster mass is about four times lower than previously calculated.

preprint2020arXiv

Massive black holes regulated by luminous blue variable mass loss and magnetic fields

We investigate the effects of mass loss during the main-sequence (MS) and post-MS phases of massive star evolution on black hole (BH) birth masses. We compute solar metallicity Geneva stellar evolution models of an 85 $M_{\odot}$ star with mass-loss rate ($\dot{M}$) prescriptions for MS and post-MS phases and analyze under which conditions such models could lead to very massive BHs. Based on the observational constraints for $\dot{M}$ of luminous stars, we discuss two possible scenarios that could produce massive BHs at high metallicity. First, if a massive BH progenitor evolves from the observed population of massive MS stars known as WNh stars, we show that its average post-MS mass-loss rate has to be less than $1\,\times10^{-5}\,M_{\odot}$/yr. However, this is lower than the typical observed mass-loss rates of luminous blue variables (LBV). Second, a massive BH progenitor could evolve from a yet undetected population of $80-85$ $M_{\odot}$ stars with strong surface magnetic fields, which could quench mass loss during the evolution. In this case, the average mass-loss rate during the post-MS LBV phase has to be less than $5\,\times10^{-5}\,M_{\odot}$/yr to produce 70 $M_{\odot}$ BHs. We suggest that LBVs that explode as SNe have large envelopes and small cores that could be prone to explosion, possibly evolving from binary interaction (either mergers or mass gainers that do not fully mix). Conversely, LBVs that directly collapse to BHs could have evolve from massive single stars or binary-star mergers that fully mix, possessing large cores that would favor BH formation.

preprint2020arXiv

Progenitors of early-time interacting supernovae

We compute an extensive set of early-time spectra of supernovae interacting with circumstellar material using the radiative transfer code CMFGEN. Our models are applicable to events observed from 1 to a few days after explosion. Using these models, we constrain the progenitor and explosion properties of a sample of 17 observed interacting supernovae at early-times. Because massive stars have strong mass loss, these spectra provide valuable information about supernova progenitors, such as mass-loss rates, wind velocities, and surface abundances. We show that these events span a wide range of explosion and progenitor properties, exhibiting supernova luminosities in the 1e8 to 1e12 Lsun range, temperatures from 10000 to 60000 K, progenitor mass-loss rates from a few 1e-4 up to 1 Msun/yr, wind velocities from 100 to 800 km/s, and surface abundances from solar-like to H-depleted. Our results suggest that many progenitors of supernovae interacting with circumstellar material have significantly increased mass-loss before explosion compared to what massive stars show during the rest of their lifetimes. We also infer a lack of correlation between surface abundances and mass-loss rates. This may point to the pre-explosion mass-loss mechanism being independent of stellar mass. We find that the majority of these events have CNO-processed surface abundances. In the single star scenario this points to a preference towards high-mass RSGs as progenitors of interacting SNe, while binary evolution could impact this conclusion. Our models are publicly available and readily applicable to analyze results from ongoing and future large scale surveys such as the Zwicky Transient Factory.

preprint2016arXiv

The fossil wind structures of Eta Carinae: changes across on 5.54-year cycle

Eta Carinae, the closest, active, massive binary containing a highly unstable LBV, exhibits expanding, compressed wind shells, seen in emission, that are spatially and spectrally resolved by \hst/\stis. Starting in June 2009, these structures were mapped across its 5.54-year, highly elliptical, binary orbit to follow temporal changes in the light of [Fe~III] 4659Å and [Fe~II] 4815Å. The emissions trace portions of fossil wind-shells, that were formed by wind-wind interactions across each cycle. Over the high-ionization state, dense arcs, photo-ionized by far ultraviolet radiation from the hot secondary, are seen in [Fe~III]. Other arcs, ionized by mid-ultraviolet radiation from the primary star, are seen in [Fe II]. The [Fe III] structures tend to be interior to [Fe II] structures that trace extensive, less disturbed primary wind. During the brief periastron passage when the secondary plunges deep into the primary's extremely dense wind, on the far side of primary star, high-ionization [Fe III] structures fade and reappear in [Fe II]. Multiple fossil wind-structures were traced across the 5.7-year monitoring interval. The strong similarity of the expanding [Fe II] shells suggests that the wind and photo-ionization properties of the massive binary have not changed substantially from one orbit to the next over the past several orbital cycles. These observations trace structures that can be used to test three-dimensional hydrodynamical and radiative-transfer models of massive, interacting winds. They also provide a baseline for following future changes in Eta Carinae, especially of its winds and photoionization properties.

preprint2015arXiv

Early Emission from the Type IIn Supernova 1998S at High Resolution

The well-studied Type IIn supernova (SN) 1998S is often dubbed the prototypical SN IIn, and it provides a unique opportunity to study its progenitor star from within as the SN lights up dense circumstellar material (CSM) launched from the progenitor. Here we present a Keck HIRES spectrum of SN 1998S taken within a few days after core collapse - both the earliest high-resolution ($Δλ< 1.0\,Å$) spectrum published of a SN IIn and the earliest published spectrum of SN 1998S. Modern SN studies achieve impressively short turn-around times between SN detection and the first observed spectrum, but high-resolution spectra of very young supernovae are rare; the unique spectrum presented here provides a useful case study for observations of other young SN systems including SN 2013cu, which displayed a remarkably similar spectrum when very young. We examine the fully resolved emission-line profiles of SN 1998S, finding evidence for extreme mass loss from the progenitor at velocities much less than those characteristic of Wolf-Rayet (WR) stars. We model our high-resolution SN 1998S spectrum using the radiative-transfer code CMFGEN and explore the composition, density, and velocity gradients within the SN system. We find a mass-loss rate of $6.0 \times 10^{-3} M_{\odot}$ per year during the $\sim$15 yr before core collapse, while other studies indicate a much lower rate at earlier times (>15yr before core collapse). A comparison with a spectrum of SN 2013cu indicates many similarities, though SN 2013cu was of Type IIb - indicating that very different supernovae can arise from progenitors with extreme mass loss in the last few years before explosion.

preprint2014arXiv

Early-time spectra of supernovae and their precursor winds: the luminous blue variable/yellow hypergiant progenitor of SN 2013cu

We present the first quantitative spectroscopic modeling of an early-time supernova that interacts with its progenitor wind. Using the radiative transfer code CMFGEN, we investigate the recently-reported 15.5 h post-explosion spectrum of the type IIb SN 2013cu. For the first time, we are able to directly measure the chemical abundances of a SN progenitor and find a relatively H-rich wind, with H and He abundances (by mass) of X=0.46 +- 0.2 and Y=0.52 +- 0.2, respectively. The wind is enhanced in N and depleted in C relative to solar values (mass fractions of 8.2e-3 and 1e-5). We obtain that a dense wind/circumstellar medium, with a mass-loss rate of Mdot= 3e-3 Msun/yr and wind velocity vwind=100 km/s, surrounds the star at the pre-SN stage. These values are lower than previous analytical estimates, although we find Mdot/vinf consistent with previous work. We also compute a CMFGEN model to constrain the progenitor spectral type and find that the high Mdot and low vwind imply that the star had an effective temperature of ~8000 K immediately before the SN explosion. Our models suggest that the progenitor was either an unstable luminous blue variable or a yellow hypergiant undergoing an eruptive phase, and rule out a WR star. We classify the post-explosion spectra at 15.5 h as XWN5(h) and advocate for the use of the prefix `X' (eXplosion) to avoid confusion between post-explosion, non-stellar spectra with those of massive stars. We show that the progenitor spectral type is significantly different than the early post-explosion spectral type owing to the huge differences in the ionization structure before and after the SN event. We find the following temporal evolution: LBV/YHG -> XWN5(h) -> SN IIb. Future early-time spectroscopy in the UV will give access to additional spectroscopic diagnostics and further constrain the properties of SN precursors, such as their metallicities.

preprint2013arXiv

Constraints on decreases in Eta Carinae's mass loss from 3D hydrodynamic simulations of its binary colliding winds

Recent work suggests that the mass-loss rate of the primary star (Eta A) in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results from large- (r=1545au) and small- (r=155au) domain, 3D smoothed particle hydrodynamic (SPH) simulations of Eta Car's colliding winds for 3 Eta A mass-loss rates (2.4, 4.8, and 8.5 x 10^-4 M_sun/yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling, and radiative forces. We find that Eta A's mass-loss rate greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star (Eta B) switches from the adiabatic to the radiative-cooling regime during periastron passage. The SPH simulations together with 1D radiative transfer models of Eta A's spectra reveal that a factor of 2 or more drop in Eta A's mass-loss rate should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in Eta A's mass-loss rate was likely by a factor < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line-of-sight to the system lies close to the dense walls of the WWC zone. A modest decrease in Eta A's mass-loss rate may be responsible, but changes in the wind/stellar parameters of Eta B cannot yet be fully ruled out. We suggest observations during Eta Car's next periastron in 2014 to further test for decreases in Eta A's mass-loss rate. If Eta A's mass-loss rate is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

preprint2013arXiv

Episodic modulations in supernova radio light curves from luminous blue variable supernova progenitor models

<Context> Ideally, one would like to know which type of core-collapse SNe is produced by different progenitors and the channels of stellar evolution leading to these progenitors. These links have to be very well known to use the observed frequency of different types of SN events for probing the star formation rate and massive star evolution in different types of galaxies. <Aims> We investigate the link between LBV as SN progenitors and the appearance of episodic radio light curve modulations of the SN event. <Methods> We use the 20Msun and 25Msun models with rotation at solar metallicity, part of an extended grid of stellar models computed by the Geneva team. At their pre-SN stage, these two models have recently been shown to have spectra similar to those of LBV stars and possibly explode as Type IIb SNe. Based on the wind properties before the explosion, we derive the density structure of their circumstellar medium. This structure is used as input for computing the SN radio light curve. <Results> We find that the 20Msun model shows radio light curves with episodic luminosity modulations, similar to those observed in some Type IIb SNe. This occurs because the evolution of the 20Msun model terminates in a region of the HR diagram where radiative stellar winds present strong density variations, caused by the bistability limit. The 25Msun model, ending its evolution in a zone of the HR diagram where no change of the mass-loss rates is expected, presents no such modulations in its radio SN light curve. <Conclusions> Our results reinforce the link between SN progenitors and LBV stars. We also confirm the existence of a physical mechanism for a single star to have episodic radio light curve modulations. In the case of the 25Msun progenitors, we do not obtain modulations in the radio light curve, but our models may miss some outbursting behavior in the late stages of massive stars.

preprint2013arXiv

Fundamental properties of core-collapse Supernova and GRB progenitors: predicting the look of massive stars before death

We investigate the fundamental properties of core-collapse Supernova (SN) progenitors from single stars at solar metallicity. We combine Geneva stellar evolutionary models with initial masses of Mini=20-120 Msun with atmospheric/wind models using CMFGEN. We provide synthetic photometry and high-resolution spectra of hot stars at the pre-SN stage. For Mini=9-20 Msun, we supplement our analysis using publicly available MARCS model atmospheres of RSGs. We employ observational criteria of spectroscopic classification and find that massive stars, depending on Mini and rotation, end their lives as red supergiants (RSG), yellow hypergiants (YHG), luminous blue variables (LBV), and Wolf-Rayet (WR) stars of the WN and WO spectral types. For rotating models, we obtain the following types of SN progenitors: WO1-3 (Mini <= 32 Msun), WN10-11 (25 < Mini < 32 Msun), LBV (20 <= Mini < 25 Msun), G1 Ia+ (18 < Mini < 20 Msun), and RSGs (9 <= Mini <= 18 Msun). For non-rotating models, we find spectral types WO1-3 (Mini>40 Msun), WN7-8 (25 < Mini <= 40 Msun), WN11h/LBV (20 < Mini <= 25 Msun), and RSGs (9 <= Mini <= 20 Msun). Our rotating models indicate that SN IIP progenitors are all RSG, SN IIL/b progenitors are 56% LBVs and 44% YHGs, SN Ib progenitors are 96% WN10-11 and 4% WOs, and SN Ic progenitors are all WO stars. We find that not necessarily the most massive and luminous SN progenitors are the brighter ones in a given filter. We show that SN IIP progenitors (RSGs) are bright in the RIJHK_S filters and faint in the UB filters. SN IIL/b progenitors (LBVs and YHGs), and SN Ib progenitors (WNs) are relatively bright in optical/IR filters, while SN Ic progenitors (WOs) are faint in all optical filters. We argue that SN Ib and Ic progenitors from single stars should be undetectable in the available pre-explosion images with the current magnitude limits, in agreement with observational results.[abridged]

preprint2013arXiv

Massive star evolution: Luminous Blue Variables as unexpected Supernova progenitors

Stars more massive than about 8 Msun end their lives as a Supernova (SN), an event of fundamental importance Universe-wide. Theoretically, these stars have been expected to be either at the red supergiant, blue supergiant, or Wolf-Rayet stage before the explosion. We performed coupled stellar evolution and atmospheric modeling of stars with initial masses between 20 Msun and 120 Msun. We found that the 20 Msun and 25 Msun rotating models, before exploding as SN, have spectra that do not resemble any of the aforementioned classes of massive stars. Rather, they have remarkable similarities with rare, unstable massive stars known as Luminous Blue Variables (LBV). While observations show that some SNe seem to have had LBVs as progenitors, no theoretical model had yet predicted that a star could explode at this stage. Our models provide theoretical support for relatively low-luminosity LBVs exploding as SN in the framework of single stellar evolution. This is a significant shift in paradigm, meaning that a fraction of LBVs could be the end stage of massive star evolution, rather than a transitory evolutionary phase. We suggest that type IIb SN could have LBV as progenitors, and a prime example could be SN 2008ax.

preprint2013arXiv

Progenitors of supernova Ibc: a single Wolf-Rayet star as the possible progenitor of the SN Ib iPTF13bvn

Core-collapse supernova (SN) explosions mark the end of the tumultuous life of massive stars. Determining the nature of their progenitors is a crucial step towards understanding the properties of SNe. Until recently, no progenitor has been directly detected for SN of type Ibc, which are believed to come from massive stars that lose their hydrogen envelope through stellar winds and from binary systems where the companion has stripped the H envelope from the primary. Here we analyze recently reported observations of iPTF13bvn, which could possibly be the first detection of a SN Ib progenitor based on pre-explosion images. Very interestingly, the recently published Geneva models of single stars can reproduce the observed photometry of the progenitor candidate and its mass-loss rate, confirming a recently proposed scenario. We find that a single WR star with initial mass in the range 31-35 Msun fits the observed photometry of the progenitor of iPTF13bvn. The progenitor likely has a luminosity of log (L/Lsun)~5.55, surface temperature ~45000 K, and mass of ~10.9 Msun at the time of explosion. Our non-rotating 32 Msun model overestimates the derived radius of the progenitor, although this could likely be reconciled with a fine-tuned model of a more massive (between 40 and 50 Msun), hotter, and luminous progenitor. Our models indicate a very uncertain ejecta mass of ~8 Msun, which is higher than the average of the SN Ib ejecta mass that is derived from the lightcurve (2-4 Msun). This possibly high ejecta mass could produce detectable effects in the iPTF13bvn lightcurve and spectrum. If the candidate is indeed confirmed to be the progenitor, our results suggest that stars with relatively high initial masses (>30 Msun) can produce visible SN explosions at their deaths and do not collapse directly to a black hole.

preprint2012arXiv

A companion as the cause of latitude-dependent effects in the wind of Eta Carinae

We analyze spatially resolved spectroscopic observations of the Eta Carinae binary system obtained with HST/STIS. Eta Car is enshrouded by the dusty Homunculus nebula, which scatters light emitted by the central binary and provides a unique opportunity to study a massive binary system from different vantage points. We investigate the latitudinal and azimuthal dependence of H$α$ line profiles caused by the presence of a wind-wind collision (WWC) cavity created by the companion star. Using two-dimensional radiative transfer models, we find that the wind cavity can qualitatively explain the observed line profiles around apastron. Regions of the Homunculus which scatter light that propagated through the WWC cavity show weaker or no H alpha absorption. Regions scattering light that propagated through a significant portion of the primary wind show stronger P Cygni absorption. Our models overestimate the H alpha absorption formed in the primary wind, which we attribute to photoionization by the companion, not presently included in the models. We can qualitatively explain the latitudinal changes that occur during periastron, shedding light on the nature of Eta Car's spectroscopic events. Our models support the idea that during the brief period of time around periastron when the primary wind flows unimpeded toward the observer, H alpha absorption occurs in directions toward the central object and Homunculus SE pole, but not toward equatorial regions close to the Weigelt blobs. We suggest that observed latitudinal and azimuthal variations are dominated by the companion star via the WWC cavity, rather than by rapid rotation of the primary star.

preprint2012arXiv

A Lighthouse Effect in Eta Carinae

We present a new model for the behavior of scattered time-dependent, asymmetric near-UV emission from the nearby ejecta of η Car. Using a 3-D hydrodynamical simulation of η Car's binary colliding winds, we show that the 3-D binary orientation derived by Madura et al. (2012) is capable of explaining the asymmetric near-UV variability observed in the Hubble Space Telescope Advanced Camera for Surveys/High Resolution Camera (HST ACS/HRC) F220W images of Smith et al. (2004b). Models assuming a binary orientation with i ~ 130 to 145 degrees, ω ~ 230 to 315 degrees, PAz ~ 302 to 327 degrees are consistent with the observed F220W near-UV images. We find that the hot binary companion does not significantly contribute to the near-UV excess observed in the F220W images. Rather, we suggest that a bore-hole effect and the reduction of Fe II optical depths inside the wind-wind collision cavity carved in the extended photosphere of the primary star lead to the time-dependent directional illumination of circum-binary material as the companion moves about in its highly elliptical orbit.

preprint2012arXiv

On the influence of the companion star in Eta Carinae: 2D radiative transfer modeling of the ultraviolet and optical spectra

We present 2D radiative transfer modeling of the Eta Carinae binary system accounting for the presence of a wind-wind collision (WWC) cavity carved in the optically-thick wind of the primary star. By comparing synthetic line profiles with HST/STIS spectra obtained near apastron, we show that the WWC cavity has a strong influence on multi-wavelength diagnostics. This influence is regulated by the modification of the optical depth in the continuum and spectral lines. We find that H-alpha, H-beta, and Fe II lines are the most affected by the WWC cavity, since they form over a large volume of the primary wind. These spectral lines depend on latitude and azimuth since, according to the orientation of the cavity, different velocity regions of a spectral line are affected. For 2D models with orientation corresponding to orbital inclination angle 110deg < i < 140deg and longitude of periastron 210deg < omega < 330deg, the blueshifted and zero-velocity regions of the line profiles are the most affected. These orbital orientations are required to simultaneously fit the UV and optical spectrum of Eta Car, for a half-opening angle of the cavity in the range 50-70deg. We find that the excess P-Cygni absorption seen in H-alpha, H-beta and optical Fe II lines in spherical models becomes much weaker or absent in the 2D models, in agreement with the observations. The observed UV spectrum of Eta Car, dominated by Fe II absorption lines, is superbly reproduced by our 2D cavity models. Small discrepancies still remain, as H-gamma and H-delta absorptions are overestimated by our models. We suggest that photoionization of the wind of the primary by the hot companion star is responsible for the weak absorption seen in these lines. Our CMFGEN models indicate that the primary star has a mass-loss rate of 8.5x10e-4 Msun/yr and wind terminal velocity of 420 km/s around the 2000 apastron.

preprint2012arXiv

Summary of IAU GA SpS5 - III. Matter ejection and feedback

The last part of SpS5 dealt with the circumstellar environment. Structures are indeed found around several types of massive stars, such as blue and red supergiants, as well as WRs and LBVs. As shown in the last years, the potential of IR for their study is twofold: first, IR can help discover many previously unknown nebulae, leading to the identification of new massive stars as their progenitors; second, IR can help characterize the nebular features. Current and new IR facilities thus pave the way to a better understanding of the feedback from massive stars.

preprint2011arXiv

Constraining the Absolute Orientation of Eta Carinae's Binary Orbit: A 3-D Dynamical Model for the Broad [Fe III] Emission

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA θ that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38 degrees, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i \approx 130 to 145 degrees, θ \approx -15 to +30 degrees, and an orbital axis projected on the sky at a PA \approx 302 to 327 degrees east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta B, thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.

preprint2011arXiv

Imaging with HST the time evolution of Eta Carinae's colliding winds

We report new HST/STIS observations that map the high-ionization forbidden line emission in the inner arcsecond of Eta Car, the first that fully image the extended wind-wind interaction region of the massive colliding wind binary. These observations were obtained after the 2009.0 periastron at orbital phases 0.084, 0.163, and 0.323 of the 5.54-year spectroscopic cycle. We analyze the variations in brightness and morphology of the emission, and find that blue-shifted emission (-400 to -200 km s-1) is symmetric and elongated along the northeast-southwest axis, while the red-shifted emission (+100 to +200 km s-1) is asymmetric and extends to the north- northwest. Comparison to synthetic images generated from a 3-D dynamical model strengthens the 3-D orbital orientation found by Madura et al. (2011), with an inclination i \approx 138\degree, argument of periapsis ω \approx 270\degree, and an orbital axis that is aligned at the same PA on the sky as the symmetry axis of the Homunculus, 312\degree. We discuss the potential that these and future mappings have for constraining the stellar parameters of the companion star and the long-term variability of the system.

preprint2011arXiv

Multi-Wavelength Implications of the Companion Star in Eta Carinae

Eta Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e \sim 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of Eta Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of Eta Car's observables at multiple wavelengths.

preprint2011arXiv

On the nature of the prototype LBV AG Carinae II. Witnessing a massive star evolving close to the Eddington and bistability limits

We show that the significantly different effective temperatures (Teff) achieved by the luminous blue variable AG Carinae during the consecutive visual minima of 1985-1990 (Teff= 22,800 K) and 2000-2001 (Teff=17,000 K) place the star on different sides of the bistability limit, which occurs in line-driven stellar winds around Teff ~ 21,000 K. Decisive evidence is provided by huge changes in the optical depth of the Lyman continuum in the inner wind as Teff changes during the S Dor cycle. These changes cause different Fe ionization structures in the inner wind. The bistability mechanism is also related to the different wind parameters during visual minima: the wind terminal velocity was 2-3 times higher and the mass-loss rate roughly two times smaller in 1985-1990 than in 2000-2003. We obtain a projected rotational velocity of 220 +- 50 km/s during 1985-1990 which, combined with the high luminosity (L=1.5x10^6 Lsun), puts AG Car extremely close to the Eddington limit modified by rotation (Omega-Gamma limit): for an inclination angle of 90 deg, Gamma_Omega > 1.0 for M< 60 Msun. Based on evolutionary models and mass budget, we obtain an initial mass of ~ 100 Msun and a current mass of 60-70 Msun for AG Car. Therefore, AG Car is close to, if not at, the Omega-Gamma limit during visual minimum. Assuming M=70 Msun, we find that Gamma_Omega decreases from 0.93 to 0.72 as AG Car expands toward visual maximum, suggesting that the star is not above the Eddington limit during maximum phases.

preprint2011arXiv

The bi-stability jump as the origin for multiple P-Cygni absorption components in Luminous Blue Variables

Luminous Blue Variables (LBVs) oftentimes show double-troughed absorption in their strong Halpha lines, which are as yet not understood. Intriguingly, the feature has also been seen in the interacting supernova SN 2005gj, which was for this reason suggested to have an LBV progenitor. Our aims are to understand the double-troughed absorption feature in LBVs and investigate whether this phenomenon is related to wind variability. To this purpose, we perform time-dependent radiative transfer modeling using CMFGEN. We find that abrupt changes in the wind-terminal velocity - as expected from the bi-stability jump - are required to explain the double-troughed absorption profiles in LBVs. This strengthens scenarios that discuss the link between LBVs and SNe utilizing the progenitor's wind variability resulting from the bi-stability jump. We also discuss why the presence of double-troughed P-Cygni components may become an efficient tool to detect extra-galactic LBVs and how to analyze their mass-loss history on the basis of just one single epoch of spectral observations.

preprint2010arXiv

Detection of high-velocity material from the wind-wind collision zone of Eta Carinae across the 2009.0 periastron passage

We report near-IR spectroscopic observations of the Eta Carinae massive binary system during 2008-2009 using VLT/CRIRES. We detect a strong, broad absorption wing in He I 10833 extending up to -1900 km/s across the 2009.0 spectroscopic event. Archival HST/STIS ultraviolet and optical data shows a similar high-velocity absorption (up to -2100 km/s) in the UV resonance lines of Si IV 1394, 1403 across the 2003.5 event. UV lines from low-ionization species, such as Si II 1527, 1533 and C II 1334, 1335, show absorption up to -1200 km/s, indicating that the absorption with v from -1200 to -2100 km/s originates in a region markedly faster and more ionized than the nominal wind of the primary star. Observations obtained at the OPD/LNA during the last 4 spectroscopic cycles (1989-2009) also display high-velocity absorption in He I 10833 during periastron. Based on the OPD/LNA dataset, we determine that material with v < -900 km/s is present in the phase range 0.976 < phi < 1.023 of the spectroscopic cycle, but absent in spectra taken at phi < 0.947 and phi > 1.049. Therefore, we constrain the duration of the high-velocity absorption to be 95 to 206 days (or 0.047 to 0.102 in phase). We suggest that the high-velocity absorption originates from shocked gas in the wind-wind collision zone, at distances of 15 to 45 AU in the line-of-sight to the primary star. Using 3-D hydrodynamical simulations of the wind-wind collision zone, we find that the dense high-velocity gas is in the line-of-sight to the primary star only if the binary system is oriented in the sky so that the companion is behind the primary star during periastron, corresponding to a longitude of periastron of omega ~ 240 to 270 degrees. We study a possible tilt of the orbital plane relative to the Homunculus equatorial plane and conclude that our data are broadly consistent with orbital inclinations in the range i=40 to 60 degrees.

preprint2010arXiv

Is Eta Carinae a fast rotator, and how much does the companion influence the inner wind structure?

We analyze interferometric measurements of the Luminous Blue Variable Eta Carinae with the goal of constraining the rotational velocity of the primary star and probing the influence of the companion. Using 2-D radiative transfer models of latitude-dependent stellar winds, we find that prolate wind models with a ratio of the rotational velocity (vrot) to the critical velocity (vcrit) of W=0.77-0.92, inclination angle of i=60-90 degrees, and position angle PA=108-142 degrees reproduce simultaneously K-band continuum visibilities from VLTI/VINCI and closure phase measurements from VLTI/AMBER. Interestingly, oblate models with W=0.73-0.90 and i=80-90 degrees produce similar fits to the interferometric data, but require PA=210-230 degrees. Therefore, both prolate and oblate models suggest that the rotation axis of the primary star is not aligned with the Homunculus polar axis. We also compute radiative transfer models of the primary star allowing for the presence of a cavity and dense wind-wind interaction region created by the companion star. We find that the wind-wind interaction has a significant effect on the K-band image mainly via free-free emission from the compressed walls and, for reasonable model parameters, can reproduce the VLTI/VINCI visibilities taken at phase 0.92-0.93. We conclude that the density structure of the primary wind can be sufficiently disturbed by the companion, thus mimicking the effects of fast rotation in the interferometric observables. Therefore, fast rotation may not be the only explanation for the interferometric observations. Intense temporal monitoring and 3-D modeling are needed to resolve these issues.

preprint2009arXiv

Bona-fide, strong-variable galactic Luminous Blue Variable stars are fast rotators: detection of a high rotational velocity in HR Carinae

We report optical observations of the Luminous Blue Variable (LBV) HR Carinae which show that the star has reached a visual minimum phase in 2009. More importantly, we detected absorptions due to Si IV 4088-4116 Angstroms. To match their observed line profiles from 2009 May, a high rotational velocity of vrot=150 +- 20 km/s is needed (assuming an inclination angle of 30 degrees), implying that HR Car rotates at ~0.88 +- 0.2 of its critical velocity for break-up (vcrit). Our results suggest that fast rotation is typical in all strong-variable, bona-fide galactic LBVs, which present S Dor-type variability. Strong-variable LBVs are located in a well-defined region of the HR diagram during visual minimum (the "LBV minimum instability strip"). We suggest this region corresponds to where vcrit is reached. To the left of this strip, a forbidden zone with vrot/vcrit>1 is present, explaining why no LBVs are detected in this zone. Since dormant/ex LBVs like P Cygni and HD 168625 have low vrot, we propose that LBVs can be separated in two groups: fast-rotating, strong-variable stars showing S-Dor cycles (such as AG Car and HR Car) and slow-rotating stars with much less variability (such as P Cygni and HD 168625). We speculate that SN progenitors which had S-Dor cycles before exploding (such as in SN 2001ig, SN 2003bg, and SN 2005gj) could have been fast rotators. We suggest that the potential difficulty of fast-rotating Galactic LBVs to lose angular momentum is an additional evidence that such stars could explode during the LBV phase.

preprint2009arXiv

On the nature of the prototype LBV AG Carinae I. Fundamental parameters during visual minimum phases and changes in the bolometric luminosity during the S-Dor cycle

We present a detailed spectroscopic analysis of the luminous blue variable AG Carinae during the last two visual minimum phases of its S-Dor cycle (1985-1990 and 2000-2003). The analysis reveals an overabundance of He, N, and Na, and a depletion of H, C, and O, on the surface of AG Car, indicating the presence of CNO-processed material. Furthermore, the ratio N/O is higher on the stellar surface than in the nebula. We found that the minimum phases of AG Car are not equal to each other, since we derived a noticeable difference between the maximum effective temperature achieved during 1985-1990 (22,800 K) and 2000-2001 (17,000 K). While the wind terminal velocity was 300 km/s in 1985-1990, it was as low as 105 km/s in 2001. The mass-loss rate, however, was lower from 1985-1990 (1.5 x 10^(-5) Msun/yr) than from 2000-2001 (3.7 x 10^(-5) Msun/yr). We found that the wind of AG Car is significantly clumped (f=0.10 - 0.25) and that clumps must be formed deep in the wind. We derived a bolometric luminosity of 1.5 x 10^6 Lsun during both minimum phases which, contrary to the common assumption, decreases to 1.0 x 10^6 Lsun as the star moves towards maximum flux in the V band. Assuming that the decrease in the bolometric luminosity of AG Car is due to the energy used to expand the outer layers of the star (Lamers 1995), we found that the expanding layers contain roughly 0.6 - 2 Msun. Such an amount of mass is an order of magnitude lower than the nebular mass around AG Car, but is comparable to the nebular mass found around lower-luminosity LBVs and to that of the Little Homunculus of Eta Car. If such a large amount of mass is indeed involved in the S Dor-type variability, we speculate that such instability could be a failed Giant Eruption, with several solar masses never becoming unbound from the star.(abridged)