Photoionization of Co+ and electron-impact excitation of Co2 + using the Dirac R-matrix method

Modelling of massive stars and supernovae (SNe) plays a crucial role in understanding galaxies. From this modelling we can derive fundamental constraints on stellar evolution, mass-loss processes, mixing, and the products of nucleosynthesis. Proper account must be taken of all important processes that populate and depopulate the levels (collisional excitation, de-excitation, ionization, recombination, photoionization, bound–bound processes). For the analysis of Type Ia SNe and core collapse SNe (Types Ib, Ic and II) Fe group elements are particularly important. Unfortunately little data is currently available and most noticeably absent are the photoionization cross-sections for the Fe-peaks which have high abundances in SNe. Important interactions for both photoionization and electron-impact excitation are calculated using the relativistic Dirac atomic R-matrix codes (DARC) for low-ionization stages of Cobalt. All results are calculated up to photon energies of 45 eV and electron energies up to 20 eV. The wavefunction representation of Co III has been generated using GRASP0 by including the dominant 3d⁷, 3d⁶[4s, 4p], 3p⁴3d⁹ and 3p⁶3d⁹ configurations, resulting in 292 fine structure levels. Electron-impact collision strengths and Maxwellian averaged effective collision strengths across a wide range of astrophysically relevant temperatures are computed for Co III. In addition, statistically weighted level-resolved ground and metastable photoionization cross-sections are presented for Co II and compared directly with existing work.

The Status of Multi-Dimensional Core-Collapse Supernova Models

Models of neutrino-driven core-collapse supernova explosions have matured considerably in recent years. Explosions of low-mass progenitors can routinely be simulated in 1D, 2D, and 3D. Nucleosynthesis calculations indicate that these supernovae could be contributors of some lighter neutron-rich elements beyond iron. The explosion mechanism of more massive stars remains under investigation, although first 3D models of neutrino-driven explosions employing multi-group neutrino transport have become available. Together with earlier 2D models and more simplified 3D simulations, these have elucidated the interplay between neutrino heating and hydrodynamic instabilities in the post-shock region that is essential for shock revival. However, some physical ingredients may still need to be added/improved before simulations can robustly explain supernova explosions over a wide range of progenitors. Solutions recently suggested in the literature include uncertainties in the neutrino rates, rotation, and seed perturbations from convective shell burning. We review the implications of 3D simulations of shell burning in supernova progenitors for the ‘perturbations-aided neutrino-driven mechanism,’ whose efficacy is illustrated by the first successful multi-group neutrino hydrodynamics simulation of an 18 solar mass progenitor with 3D initial conditions. We conclude with speculations about the impact of 3D effects on the structure of massive stars through convective boundary mixing.

Dynamical stability of nascent neutron stars

This thesis presents a study of the dynamical stability of nascent neutron stars resulting from the accretion induced collapse of rapidly rotating white dwarfs.

Chapter 2 and part of Chapter 3 study the equilibrium models for these neutron stars. They are constructed by assuming that the neutron stars have the same masses, angular momenta, and specific angular momentum distributions as the pre-collapse white dwarfs. If the pre-collapse white dwarf is rapidly rotating, the collapsed object will contain a high density central core of size about 20 km, surrounded by a massive accretion torus extending to hundreds of kilometers from the rotation axis. The ratio of the rotational kinetic energy to gravitational binding energy, β, of these neutron stars is all found to be less than 0.27.

Chapter 3 studies the dynamical stability of these neutron stars by numerically evolving the linearized hydrodynamical equations. A dynamical bar-mode instability is observed when the β of the star is greater than the critical value β_d ≈ 0.25. It is expected that the unstable mode will persist until a substantial amount of angular momentum is carried away by gravitational radiation. The detectability of these sources is studied and it is estimated that LIGO II is unlikely to detect them unless the event rate is greater than 10^-6/year/galaxy.

All the calculations on the structure and stability of the neutron stars in Chapters 2 and 3 are carried out using Newtonian hydrodynamics and gravity. Chapter 4 studies the relativistic effects on the structure of these neutron stars. New techniques are developed and used to construct neutron star models to the first post-Newtonian (1PN) order. The structures of the 1PN models are qualitatively similar to the corresponding Newtonian models, but the values of β are somewhat smaller. The maximum β for these 1PN neutron stars is found to be 0.24, which is 8% smaller than the Newtonian result (0.26). However, relativistic effects will also change the critical value β_d. A detailed post-Newtonian stability analysis has yet to be carried out to study the relativistic effects on the dynamical stability of these neutron stars.

The massive binary companion star to the progenitor of supernova 1993J

The massive star that underwent a collapse of its core to produce supernova (SN)1993J was subsequently identified as a non-variable red supergiant star in images of the galaxy M81 taken before explosion(1, 2). It showed an excess in ultraviolet and B-band colours, suggesting either the presence of a hot, massive companion star or that it was embedded in an unresolved young stellar association1. The spectra of SN1993J underwent a remarkable transformation from the signature of a hydrogen-rich type II supernova to one of a helium-rich (hydrogen-deficient) type Ib(3, 4). The spectral and photometric peculiarities were best explained by models in which the 13�20 solar mass supergiant had lost almost its entire hydrogen envelope to a close binary companion(5, 6, 7), producing a 'type IIb' supernova, but the hypothetical massive companion stars for this class of supernovae have so far eluded discovery. Here we report photometric and spectroscopic observations of SN1993J ten years after the explosion. At the position of the fading supernova we detect the unambiguous signature of a massive star: the binary companion to the progenitor.

Ruling out a massive asymptotic giant-branch star as the progenitor of supernova 2005cs

We calculate the predicted UBVRIJHK absolute magnitudes for models of supernova progenitors and apply the result to the case of supernova 2005cs. We agree with previous results that the initial mass of the star was low, around 6 to 8 M-circle dot. However, such stars are thought to go through a second dredge-up to become asymptotic giant branch (AGB) stars. We show that had this occurred to the progenitor of 2005cs it would have been observed in JHK pre-explosion images. The progenitor was not detected in these bands and therefore we conclude that it was not an AGB star. Furthermore, if some AGB stars do produce supernovae they will have a clear signature in pre-explosion near-infrared images. Electron-capture supernovae are thought to occur in AGB stars, hence the implication is that 2005cs was not an electron-capture supernova but was the collapse of an iron core.

Blue luminous stars in nearby galaxies: Quantitative spectral analysis of M33 B-type supergiant stars

We present the detailed spectral analysis of a sample of M33 B-type supergiant stars, aimed at the determination of their fundamental parameters and chemical composition. The analysis is based on a grid of non-LTE metal line-blanketed model atmospheres including the effects of stellar winds and spherical extension computed with the code FASTWIND. Surface abundance ratios of C, N, and O are used to discuss the chemical evolutionary status of each individual star. The comparison of observed stellar properties with theoretical predictions of massive star evolutionary models shows good agreement within the uncertainties of the analysis. The spatial distribution of the sample allows us to investigate the existence of radial abundance gradients in the disk of M33. The comparison of stellar and H II region O abundances ( based on direct determinations of the electron temperature of the nebulae) shows good agreement. Using a simple linear radial representation, the stellar oxygen abundances result in a gradient of -0.0145 +/- 0.005 dex arcmin(-1) (or -0.06 +/- 0.02 dex kpc(-1)) up to a distance equal to similar to 1.1 times the isophotal radius of the galaxy. A more complex representation cannot be completely discarded by our stellar sample. The stellar Mg and Si abundances follow the trend displayed by O abundances, although with shallower gradients. These differences in gradient slope cannot be explained at this point. The derived abundances of the three alpha-elements yield solar metallicity in the central regions of the disk of M33. A comparison with recent planetary nebula data from Magrini and coworkers indicates that the disk of M33 has not suffered from a significant O enrichment in the last 3 Gyr.

X-shooting Herbig Ae/Be stars: Accretion probed by near-infrared He I emission

The Herbig Ae/Be stars are intermediate mass pre-main sequence stars that bridge the gap between the low mass T Tauri stars and the Massive Young Stellar Objects. In this mass range, the acting star forming mechanism switches from magnetically controlled accretion to an as yet unknown mechanism, but which is likely to be direct disk accretion onto the star. We observed a large sample of Herbig Ae/Be stars with X-shooter to address this issue from a multi-wavelength perspective. It is the largest such study to date, not only because of the number of objects involved, but also because of the large wavelength coverage from the blue to the near-infrared. This allows many accretion diagnostics to be studied simultaneously. By correlating the various properties with mass, temperature and age, we aim to determine where and whether the magnetically controlled mass accretion mechanism halts and the proposed direct disk accretion takes over. Here, we will give an overview of the background, present some observations and discuss our initial results. We will introduce a new accretion diagnostic for the research of Herbig Ae/Be stars, the HeI 1.083 µm line. © 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Finding exoplanets orbiting young active stars - I. Technique

Stellar activity, such as starspots, can induce radial velocity (RV) variations that can mask or even mimic the RV signature of orbiting exoplanets. For this reason RV exoplanet surveys have been unsuccessful when searching for planets around young, active stars and are therefore failing to explore an important regime which can help to reveal how planets form and migrate. This paper describes a new technique to remove spot signatures from the stellar line-profiles of moderately rotating, active stars (v sin i ranging from 10 to 50 km s(-1)). By doing so it allows planetary RV signals to be uncovered. We used simulated models of a G5V type star with differing dark spots on its surface along with archive data of the known active star HD 49933 to validate our method. The results showed that starspots could be effectively cleaned from the line-profiles so that the stellar RV jitter was reduced by more than 80 per cent. Applying this procedure to the same models and HD 49933 data, but with fake planets injected, enabled the effective removal of starspots so that Jupiter mass planets on short orbital periods were successfully recovered. These results show that this approach can be useful in the search for hot-Jupiter planets that orbit around young, active stars with a v sin i of similar to 10-50 km/s.

Interacting supernovae from photoionization-confined shells around red supergiant stars

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

Do exoplanets spin up their host stars?

We have collected initial evidence that tidal interaction between a late-type star and its close-in, massive planet can lead to a spin-up of the host star. We propose to explore this further by studying a small sample of proper motion pairs in which one of the stars is orbited by a Hot Jupiter. We will determine if the activity-estimated age appears to be strongly different for the two stars, which would indicate a tidal spin up of the Hot Jupiter host star. We propose to observe 4 such systems with Chandra/ACIS-S, and to perform a similar observation of one additional system with large angular separation using XMM-Newton/EPIC. The total proposed exposure times are 141 ks (Chandra) and 38 ks (XMM).

Planets spinning up their host stars

We have collected initial evidence that tidal interaction between a late-type star and its close-in, massive planet can lead to a spin-up of the host star. We propose to explore this further by studying a small sample of proper motion pairs in which one of the stars is orbited a Hot Jupiter. We will determine if the gyrochronal age is different for the two stars, which would indicate a tidal spin up of the planet host star. We propose to observe 3 such systems with XMM, and to perform similar Chandra observations of 3 more systems with angular separations

ON THE NATURE OF THE PROTOTYPE LUMINOUS BLUE VARIABLE AG CARINAE. II. WITNESSING A MASSIVE STAR EVOLVING CLOSE TO THE EDDINGTON AND BISTABILITY LIMITS

We show that the significantly different effective temperatures (T(eff)) achieved by the luminous blue variable AG Carinae during the consecutive visual minima of 1985-1990 (T(eff) similar or equal to 22,800 K) and 2000-2001 (T(eff) similar or equal to 17,000 K) place the star on different sides of the bistability limit, which occurs in line-driven stellar winds around T(eff) similar to 21,000 K. Decisive evidence is provided by huge changes in the optical depth of the Lyman continuum in the inner wind as T(eff) 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(-1) during 1985-1990 which, combined with the high luminosity (L(star) = 1.5 x 10(6) L(circle dot)), puts AG Car extremely close to the Eddington limit modified by rotation (Omega Gamma limit): for an inclination angle of 90 degrees, Gamma(Omega) greater than or similar to 1.0 for M(circle dot) less than or similar to 60. Based on evolutionary models and mass budget, we obtain an initial mass of similar to 100 M(circle dot) and a current mass of similar to 60-70 M(circle dot) for AG Car. Therefore, AG Car is close to, if not at, the Omega Gamma limit during visual minimum. Assuming M = 70 M(circle dot), 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.

A possible stellar metallic enhancement in post-T Tauri stars by a planetesimal bombardment

The photospheres of stars hosting planets have larger metallicity than stars lacking planets. This could be the result of a metallic star contamination produced by the bombarding of hydrogen-deficient solid bodies. In the present work we study the possibility of an earlier metal enrichment of the photospheres by means of impacting planetesimals during the first 20-30 Myr. Here we explore this contamination process by simulating the interactions of an inward migrating planet with a disc of planetesimal interior to its orbit. The results show the percentage of planetesimals that fall on the star. We identified the dependence of the planet's eccentricity (e(p)) and time-scale of migration (tau) on the rate of infalling planetesimals. For very fast migrations (tau= 10(2) and 10(3) yr) there is no capture in mean motion resonances, independently of the value of e(p). Then, due to the planet's migration the planetesimals suffer close approaches with the planet and more than 80 per cent of them are ejected from the system. For slow migrations (tau= 10(5)and 10(6) yr) the percentage of collisions with the planet decreases with the increase of the planet's eccentricity. For e(p) = 0 and 0.1 most of the planetesimals were captured in the 2:1 resonance and more than 65 per cent of them collided with the star. Whereas migration of a Jupiter mass planet to very short pericentric distances requires unrealistic high disc masses, these requirements are much smaller for smaller migrating planets. Our simulations for a slowly migrating 0.1 M-Jupiter planet, even demanding a possible primitive disc three times more massive than a primitive solar nebula, produces maximum [Fe/H] enrichments of the order of 0.18 dex. These calculations open possibilities to explain hot Jupiter exoplanet metallicities.

Stars of WIMPs

A class of boson-fermion stars, whose spin-0 and spin-1/2 constituents interact through a U(1) current-current term in the Lagrangian density, is analyzed. It is shown that it describes the low-energy behavior of a system of weakly interacting massive particles (WIMPs) from the leptonic sector of the minimal supersymmetric standard model. In this case the effective coupling constant A is related to the Fermi constant GF.

Free massive particles with total energy E

We analyze free elementary particles with a rest mass m and total energy Estars, and use Unruh-DeWitt-like detectors to calculate the associated particle detection rate in each case. The (mean) particle position is identified with the spatial average of the excitation probability of the detectors, which are supposed to cover the whole space. Our results are shown to be in harmony with general relativity classical predictions. Eventually we reconcile our conclusions with Earth-based experiments which are in good agreement with E ≥mc2.