The Atomic Physics Underlying the Spectroscopic Analysis of Massive Stars and Supernovae

We have developed a radiative transfer code, CMFGEN, which allows us to model the spectra of massive stars and supernovae. Using CMFGEN we can derive fundamental parameters such as effective temperatures and surface gravities, derive abundances, and place constraints on stellar wind properties. The last of these is important since all massive stars are losing mass via a stellar wind that is driven from the star by radiation pressure, and this mass loss can substantially influence the spectral appearance and evolution of the star. Recently we have extended CMFGEN to allow us to undertake time-dependent radiative transfer calculations of supernovae. Such calculations will be used to place constraints on the supernova progenitor, to place constraints on the supernova explosion and nucleosynthesis, and to derive distances using a physical approach called the "Expanding Photosphere Method". We describe the assumptions underlying the code and the atomic processes involved. A crucial ingredient in the code is the atomic data. For the modeling we require accurate transition wavelengths, oscillator strengths, photoionization cross-sections, collision strengths, autoionization rates, and charge exchange rates for virtually all species up to, and including, cobalt. Presently, the available atomic data varies substantially in both quantity and quality.Comment: 8 pages, 2 figures, Accepted for publication in Astrophysics & Space Scienc

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.Comment: 5 pages, 4 figures, To be published in the proceedings of the meeting 'Four Decades of Research on Massive Stars' in honor of Tony Moffat, 11-15 July 2011, Saint-Michel-des-Saints, Quebe

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.Comment: 20 pages, 14 figures, accepted for publication in MNRA

Broad Line Emission in Low-Metallicity Blue Compact Dwarf Galaxies: Evidence for Stellar Wind, Supernova and Possible AGN Activity

We present spectra of a large sample of low-metallicity blue compact dwarf galaxies which exhibit broad components in their strong emission lines, mainly in Hbeta, [O III]4959, 5007 and Halpha. Twenty-three spectra have been obtained with the MMT, 14 of which show broad emission. The remaining 21 spectra with broad emission have been selected from the Data Release 5 of the Sloan Digital Sky Survey. The most plausible origin of broad line emission is the evolution of massive stars and their interaction with the circumstellar and interstellar medium. The broad emission with the lowest H$\alpha$ luminosities (10^36 - 10^39 erg/s) is likely produced in circumstellar envelopes around hot Ofp/WN9 and/or LBV stars. The broad emission with the highest Halpha luminosities (10^40 - 10^42 erg/s) probably arises from type IIp or type IIn supernovae (SNe). It can also come from active galactic nuclei (AGN) containing intermediate-mass black holes, although we find no strong evidence for hard non-thermal radiation in our sample galaxies. The oxygen abundance in the host galaxies with SN candidates is low and varies in the range 12 + log O/H = 7.36 - 8.31. However, type IIn SN / AGN candidates are found only in galaxies with 12 + log O/H < 7.99. Spectroscopic monitoring of these type IIn SN / AGN candidates over a time scale of several years is necessary to distinguish between the two possibilities.Comment: 50 pages, 6 figures. Accepted for publication in the Astrophysical Journa

Quantitative analysis of WC stars: Constraints on neon abundances from ISO/SWS spectroscopy

Neon abundances are derived in four Galactic WC stars -- gamma Vel (WR11, WC8+O7.5III), HD156385 (WR90, WC7), HD192103 (WR135, WC8), and WR146 (WC5+O8) - using mid-infrared fine structure lines obtained with ISO/SWS. Stellar parameters for each star are derived using a non-LTE model atmospheric code (Hillier & Miller 1998) together with ultraviolet (IUE), optical (INT, AAT) and infrared (UKIRT, ISO) spectroscopy. In the case of gamma Vel, we adopt results from De Marco et al. (2000), who followed an identical approach. ISO/SWS datasets reveal the [NeIII] 15.5um line in each of our targets, while [NeII] 12.8um, [SIV] 10.5um and [SIII] 18.7um are observed solely in gamma Vel. Using a method updated from Barlow et al. (1988) to account for clumped winds, we derive Ne/He=3-4x10^-3 by number, plus S/He=6x10^-5 for gamma Vel. Neon is highly enriched, such that Ne/S in gamma Vel is eight times higher than cosmic values. However, observed Ne/He ratios are a factor of two times lower than predictions of current evolutionary models of massive stars. An imprecise mass-loss and distance were responsible for the much greater discrepancy in neon content identified by Barlow et al. Our sample of WC5--8 stars span a narrow range in T* (=55--71kK), with no trend towards higher temperature at earlier spectral type, supporting earlier results for a larger sample by Koesterke & Hamann (1995). Stellar luminosities range from 100,000 to 500,000 Lo, while 10^-5.1 < Mdot/(Mo/yr) < 10^-4.5, adopting clumped winds, in which volume filling factors are 10%. In all cases, wind performance numbers are less than 10, significantly lower than recent estimates. Carbon abundances span 0.08 < C/He < 0.25 by number, while oxygen abundances remain poorly constrained.Comment: 16 pages,7 figures accepted for MNRA

Characteristic QSO Accretion Disk Temperatures from Spectroscopic Continuum Variability

Using Sloan Digital Sky Survey (SDSS) quasar spectra taken at multiple epochs, we find that the composite flux density differences in the rest frame wavelength range 1300-6000 AA can be fit by a standard thermal accretion disk model where the accretion rate has changed from one epoch to the next (without considering additional continuum emission components). The fit to the composite residual has two free parameters: a normalizing constant and the average characteristic temperature $\bar{T}^*$. In turn the characteristic temperature is dependent on the ratio of the mass accretion rate to the square of the black hole mass. We therefore conclude that most of the UV/optical variability may be due to processes involving the disk, and thus that a significant fraction of the UV/optical spectrum may come directly from the disk.Comment: 31 pages, 8 figure

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.Comment: 18 pages, 15 figures, accepted for publication in A&A; high-resolution PDF version available also at http://www.mpifr.de/staff/jgroh/etacar.htm

Toward a simulation approach for alkene ring-closing metathesis : scope and limitations of a model for RCM

A published model for revealing solvent effects on the ring-closing metathesis (RCM) reaction of di-Et diallylmalonate 7 has been evaluated over a wider range of conditions, to assess its suitability for new applications. Unfortunately, the model is too flexible and the published rate consts. do not agree with exptl. studies in the literature. However, by fixing the values of important rate consts. and restricting the concn. ranges studied, useful conclusions can be drawn about the relative rates of RCM of different substrates, precatalyst concn. can be simulated accurately and the effect of precatalyst loading can be anticipated. Progress has also been made toward applying the model to precatalyst evaluation, but further modifications to the model are necessary to achieve much broader aims