Chemistry in AGB stars: successes and challenges

Emission and absorption line observations of molecules in late-type stars are a vital component in our understanding of stellar evolution, dust formation and mass loss in these objects. The molecular composition of the gas in the circumstellar envelopes of AGB stars reflects chemical processes in gas whose properties are strong functions of radius with density and temperature varying by more than ten and two orders of magnitude, respectively. In addition, the interstellar UV field plays a critical role in determining not only molecular abundances but also their radial distributions. In this article, I shall briefly review some recent successful approaches to describing chemistry in both the inner and outer envelopes and outline areas of challenge for the future.

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 Last Minutes of Oxygen Shell Burning in a Massive Star

We present the first 3D simulation of the last minutes of oxygen shell burning in an 18 solar mass supernova progenitor up to the onset of core collapse. A moving inner boundary is used to accurately model the contraction of the silicon and iron core according to a 1D stellar evolution model with a self-consistent treatment of core deleptonization and nuclear quasi-equilibrium. The simulation covers the full solid angle to allow the emergence of large-scale convective modes. Due to core contraction and the concomitant acceleration of nuclear burning, the convective Mach number increases to ~0.1 at collapse, and an l=2 mode emerges shortly before the end of the simulation. Aside from a growth of the oxygen shell from 0.51 to 0.56 solar masses due to entrainment from the carbon shell, the convective flow is reasonably well described by mixing length theory, and the dominant scales are compatible with estimates from linear stability analysis. We deduce that artificial changes in the physics, such as accelerated core contraction, can have precarious consequences for the state of convection at collapse. We argue that scaling laws for the convective velocities and eddy sizes furnish good estimates for the state of shell convection at collapse and develop a simple analytic theory for the impact of convective seed perturbations on shock revival in the ensuing supernova. We predict a reduction of the critical luminosity for explosion by 12--24% due to seed asphericities for our 3D progenitor model relative to the case without large seed perturbations.

From stellar to planetary composition: Galactic chemical evolution of Mg/Si mineralogical ratio

Aims. The main goal of this work is to study element ratios that are important for the formation of planets of different masses. Methods. We study potential correlations between the existence of planetary companions and the relative elemental abundances of their host stars. We use a large sample of FGK-type dwarf stars for which precise Mg, Si, and Fe abundances have been derived using HARPS high-resolution and high-quality data. Results. A first analysis of the data suggests that low-mass planet host stars show higher [Mg/Si] ratios, while giant planet hosts present [Mg/Si] that is lower than field stars. However, we found that the [Mg/Si] ratio significantly depends on metallicity through Galactic chemical evolution. After removing the Galactic evolution trend only the difference in the [Mg/Si] elemental ratio between low-mass planet hosts and non-hosts was present in a significant way. These results suggest that low-mass planets are more prevalent around stars with high [Mg/Si]. Conclusions. Our results demonstrate the importance of Galactic chemical evolution and indicate that it may play an important role in the planetary internal structure and composition. The results also show that abundance ratios may be a very relevant issue for our understanding of planet formation and evolution.

Researches on the evolution of the stellar systems ...

v. 1. On the universality of the law of gravitation and on the orbits and general characteristics of binary stars.--v. 2. The capture theory of cosmical evolution, founded on dynamical principles and illustrated by phenomena observed in the spiral nebulae, the planetary system, the double and multiple stars and clusters and the star-clouds of the Milky way.

The Evolution of Massive Young Stellar Objects in the Large Magellanic Cloud

This thesis presents an analysis of the largest catalog to date of infrared spectra of massive young stellar objects in the Large Magellanic Cloud. Evidenced by their very different spectral features, the luminous objects span a range of evolutionary states from those most embedded in their natal molecular material to those that have dissipated and ionized their surroundings to form compact HII regions and photodissociation regions. We quantify the contributions of the various spectral features using the statistical method of principal component analysis. Using this analysis, we classify the YSO spectra into several distinct groups based upon their dominant spectral features: silicate absorption (S Group), silicate absorption and fine-structure line emission (SE), polycyclic aromatic hydrocarbon (PAH) emission (P Group), PAH and fine-structure line emission (PE), and only fine-structure line emission (E). Based upon the relative numbers of sources in each category, we are able to estimate the amount of time massive YSOs spend in each evolutionary stage. We find that approximately 50% of the sources have ionic fine-structure lines, indicating that a compact HII region forms about half-way through the YSO lifetime probed in our study. Of the 277 YSOs we collected spectra for, 41 have ice absorption features, indicating they are surrounded by cold ice-bearing dust particles. We have decomposed the shape of the ice features to probe the composition and thermal history of the ice. We find that most the CO2 ice is embedded a polar ice matrix that has been thermally processed by the embedded YSO. The amount of thermal processing may be correlated with the luminosity of the YSO. Using the Australia Telescope Compact Array, we imaged the dense gas around a subsample of our sources in the HII complexes N44, N105, N113, and N159 using HCO+ and HCN as dense gas tracers. We find that the molecular material in star forming environments is highly clumpy, with clumps that range from subparsec to ~2 parsecs in size and with masses between 10^2 to 10^4 solar masses. We find that there are varying levels of star formation in the clumps, with the lower-mass clumps tending to be without massive YSOs. These YSO-less clumps could either represent an earlier stage of clump to the more massive YSO-bearing ones or clumps that will never form a massive star. Clumps with massive YSOs at their centers have masses larger than those with massive YSOs at their edges, and we suggest that the difference is evolutionary: edge YSO clumps are more advanced than those with YSOs at their centers. Clumps with YSOs at their edges may have had a significant fraction of their mass disrupted or destroyed by the forming massive star. We find that the strength of the silicate absorption seen in YSO IR spectra feature is well-correlated with the on-source HCO+ and HCN flux densities, such that the strength of the feature is indicative of the embeddedness of the YSO. We estimate that ~40% of the entire spectral sample has strong silicate absorption features, implying that the YSOs are embedded in circumstellar material for about 40% of the time probed in our study.

Orbital evolution and orbital phase resolved spectroscopy of the HMXB pulsar 4U 1538-52 with RXTE-PCA and BeppoSAX

We report here results from detailed timing and spectral studies of the high mass X-ray binary pulsar 4U 1538-52 over several binary periods using observations made with the Rossi X-ray Timing Explorer (RXTE) and BeppoSAX satellites. Pulse timing analysis with the 2003 RXTE data over two binary orbits confirms an eccentric orbit of the system. Combining the orbitial parameters determined from this observation with the earlier measurements we did not find any evidence of orbital decay in this X-ray binary. We have carried out orbital phase resolved spectroscopy to measure changes in the spectral parameters with orbital phase, particularly the absorption column density and the iron line flux. The RXTE-PCA spectra in the 3-20 keV energy range were fitted with a power law and a high energy cut-off along with a Gaussian line at similar to 6.4 keV, whereas the BeppoSAX spectra needed only a power law and Gaussian emission line at similar to 6.4keV in the restricted energy range of 0.3-10.0 keV. An absorption along the line of sight was included for both the RXTE and BeppoSAX data. The variation of the free spectral parameters over the binary orbit was investigated and we found that the variation of the column density of absorbing material in the line of sight with orbital phase is in reasonable agreement with a simple model of a spherically symmetric stellar wind from the companion star.

Angular momentum transport and evolution of lopsided galaxies

The surface brightness distribution in the majority of stellar galactic discs falls off exponentially. Often what lies beyond such a stellar disc is the neutral hydrogen gas whose distribution also follows a nearly exponential profile at least for a number of nearby disc galaxies. Both the stars and gas are commonly known to host lopsided asymmetry especially in the outer parts of a galaxy. The role of such asymmetry in the dynamical evolution of a galaxy has not been explored so far. Following Lindblad's original idea of kinematic density waves, we show that the outer part of an exponential disc is ideally suitable for hosting lopsided asymmetry. Further, we compute the transport of angular momentum in the combined stars and gas disc embedded in a dark matter halo. We show that in a pure star and gas disc, there is a transition point where the free precession frequency of a lopsided mode, Omega - kappa, changes from retrograde to prograde and this in turn reverses the direction of angular momentum flow in the disc leading to an unphysical behaviour. We show that this problem is overcome in the presence of a dark matter halo, which sets the angular momentum flow outwards as required for disc evolution, provided the lopsidedness is leading in nature. This, plus the well-known angular momentum transport in the inner parts due to spiral arms, can facilitate an inflow of gas from outside perhaps through the cosmic filaments.

Relativistic velocity - potential hydrodynamics and stellar stability

The equations of relativistic, perfect-fluid hydrodynamics are cast in Eulerian form using six scalar "velocity-potential" fields, each of which has an equation of evolution. These equations determine the motion of the fluid through the equation

U_ʋ=µ^-1 (ø,_ʋ + αβ,_ʋ + ƟS,_ʋ).

Einstein's equations and the velocity-potential hydrodynamical equations follow from a variational principle whose action is

I = (R + 16π p) (-g)^1/2 d⁴x,

where R is the scalar curvature of spacetime and p is the pressure of the fluid. These equations are also cast into Hamiltonian form, with Hamiltonian density –T₀⁰ (-g^oo)^-1/2.

The second variation of the action is used as the Lagrangian governing the evolution of small perturbations of differentially rotating stellar models. In Newtonian gravity this leads to linear dynamical stability criteria already known. In general relativity it leads to a new sufficient condition for the stability of such models against arbitrary perturbations.

By introducing three scalar fields defined by

ρ ᵴ = ∇λ + ∇x(xi + ∇xɣi)

(where ᵴ is the vector displacement of the perturbed fluid element, ρ is the mass-density, and i, is an arbitrary vector), the Newtonian stability criteria are greatly simplified for the purpose of practical applications. The relativistic stability criterion is not yet in a form that permits practical calculations, but ways to place it in such a form are discussed.

The VLT-FLAMES survey of massive stars: wind properties and evolution of hot massive stars in the Large Magellanic Cloud

We have studied the optical spectra of a sample of 28 O- and early B-type stars in the Large Magellanic Cloud, 22 of which are associated with the young star forming region N11. Our observations sample the central associations of LH9 and LH10, and the surrounding regions. Stellar parameters are determined using an automated fitting method ( Mokiem et al. 2005), which combines the stellar atmosphere code fastwind ( Puls et al. 2005) with the genetic algorithm based optimisation routine PIKAIA ( Charbonneau 1995). We derive an age of 7.0 +/- 1.0 and 3.0 +/- 1.0 Myr for LH9 and LH10, respectively. The age difference and relative distance of the associations are consistent with a sequential star formation scenario in which stellar activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type O2. From helium and hydrogen line fitting we find the hottest three of these stars to be similar to 49- 54 kK ( compared to similar to 45- 46 kK for O3 stars). Detailed determination of the helium mass fraction reveals that the masses of helium enriched dwarfs and giants derived in our spectroscopic analysis are systematically lower than those implied by non-rotating evolutionary tracks. We interpret this as evidence for efficient rotationally enhanced mixing leading to the surfacing of primary helium and to an increase of the stellar luminosity. This result is consistent with findings for SMC stars by Mokiem et al. ( 2006). For bright giants and supergiants no such mass discrepancy is found; these stars therefore appear to follow tracks of modestly or non-rotating objects. The set of programme stars was sufficiently large to establish the mass loss rates of OB stars in this Z similar to 1/2 Z(circle dot) environment sufficiently accurate to allow for a quantitative comparison with similar objects in the Galaxy and the SMC. The mass loss properties are found to be intermediate to massive stars in the Galaxy and SMC. Comparing the derived modified wind momenta D-mom as a function of luminosity with predictions for LMC metallicities by Vink et al. ( 2001) yields good agreement in the entire luminosity range that was investigated, i.e. 5.0

First stellar abundances in NGC 6822 from VLT-UVES and Keck-Hires spectroscopy

We have obtained the first high-resolution spectra of individual stars in the dwarf irregular galaxy NGC 6822. The spectra of the two A-type supergiants were obtained at the Very Large Telescope and Keck Observatories, using the Ultraviolet-Visual Echelle Spectrograph and the High Resolution Echelle Spectrometer, respectively. A detailed model atmospheres analysis has been used to determine their atmospheric parameters and elemental abundances. The mean iron abundance from these two stars is [[Fe/H]] = -0.49 +/- 0.22 (+/- 0.21),(6) with Cr yielding a similar underabundance, [[Cr/H]] = -0.50 +/- 0.20 (+/- 0.16). This confirms that NGC 6822 has a metallicity that is slightly higher than that of the SMC and is the first determination of the present-day iron group abundances in NGC 6822. The mean stellar oxygen abundance, 12 + log (O/H) = 8.36 +/- 0.19 (+/- 0.21), is in good agreement with the nebular oxygen results. Oxygen has the same underabundance as iron, [[O/ Fe]] = + 0.02 +/- 0.20 (+/- 0.21). This O/Fe ratio is very similar to that seen in the Magellanic Clouds, which supports the picture that chemical evolution occurs more slowly in these lower mass galaxies, although the O/Fe ratio is also consistent with that observed in comparatively metal-poor stars in the Galactic disk. Combining all of the available abundance observations for NGC 6822 shows that there is no trend in abundance with galactocentric distance. However, a subset of the highest quality data is consistent with a radial abundance gradient. More high-quality stellar and nebular observations are needed to confirm this intriguing possibility.

The Rossiter-McLaughlin effect reloaded: Probing the 3D spin-orbit geometry, differential stellar rotation, and the spatially-resolved stellar spectrum of star-planet systems

When a planet transits its host star, it blocks regions of the stellar surface from view; this causes a distortion of the spectral lines and a change in the line-of-sight (LOS) velocities, known as the Rossiter-McLaughlin (RM) effect. Since the LOS velocities depend, in part, on the stellar rotation, the RM waveform is sensitive to the star-planet alignment (which provides information on the system’s dynamical history). We present a new RM modelling technique that directly measures the spatially-resolved stellar spectrum behind the planet. This is done by scaling the continuum flux of the (HARPS) spectra by the transit light curve, and then subtracting the infrom the out-of-transit spectra to isolate the starlight behind the planet. This technique does not assume any shape for the intrinsic local profiles. In it, we also allow for differential stellar rotation and centre-to-limb variations in the convective blueshift. We apply this technique to HD 189733 and compare to 3D magnetohydrodynamic (MHD) simulations. We reject rigid body rotation with high confidence (>99% probability), which allows us to determine the occulted stellar latitudes and measure the stellar inclination. In turn, we determine both the sky-projected (λ ≈ −0.4 ± 0.2◦) and true 3D obliquity (ψ ≈ 7+12 −4 ◦ ). We also find good agreement with the MHD simulations, with no significant centre-to-limb variations detectable in the local profiles. Hence, this technique provides a new powerful tool that can probe stellar photospheres, differential rotation, determine 3D obliquities, and remove sky-projection biases in planet migration theories. This technique can be implemented with existing instrumentation, but will become even more powerful with the next generation of high-precision radial velocity spectrographs.

Modeling the Rossiter-McLaughlin Effect: Impact of the Convective Center-to-limb Variations in the Stellar Photosphere

bservations of the Rossiter–McLaughlin (RM) effect provide information on star–planet alignments, which can inform planetary migration and evolution theories. Here, we go beyond the classical RM modeling and explore the impact of a convective blueshift that varies across the stellar disk and non-Gaussian stellar photospheric profiles. We simulated an aligned hot Jupiter with a four-day orbit about a Sun-like star and injected center-to-limb velocity (and profile shape) variations based on radiative 3D magnetohydrodynamic simulations of solar surface convection. The residuals between our modeling and classical RM modeling were dependent on the intrinsic profile width and v sin i; the amplitude of the residuals increased with increasing v sin i and with decreasing intrinsic profile width. For slowly rotating stars the center-to-limb convective variation dominated the residuals (with amplitudes of 10 s of cm s−1 to ~1 m s−1); however, for faster rotating stars the dominant residual signature was due a non-Gaussian intrinsic profile (with amplitudes from 0.5 to 9 m s−1). When the impact factor was 0, neglecting to account for the convective center-to-limb variation led to an uncertainty in the obliquity of ~10°–20°, even though the true v sin i was known. Additionally, neglecting to properly model an asymmetric intrinsic profile had a greater impact for more rapidly rotating stars (e.g., v sin i = 6 km s−1) and caused systematic errors on the order of ~20° in the measured obliquities. Hence, neglecting the impact of stellar surface convection may bias star–planet alignment measurements and consequently theories on planetary migration and evolution.

GHASP: an H alpha kinematic survey of spiral and irregular galaxies - IX. The near-infrared, stellar and baryonic Tully-Fisher relations

We studied, for the first time, the near-infrared, stellar and baryonic Tully-Fisher relations for a sample of field galaxies taken from a homogeneous Fabry-Perot sample of galaxies [the Gassendi HAlpha survey of SPirals (GHASP) survey]. The main advantage of GHASP over other samples is that the maximum rotational velocities were estimated from 2D velocity fields, avoiding assumptions about the inclination and position angle of the galaxies. By combining these data with 2MASS photometry, optical colours, HI masses and different mass-to-light ratio estimators, we found a slope of 4.48 +/- 0.38 and 3.64 +/- 0.28 for the stellar and baryonic Tully-Fisher relation, respectively. We found that these values do not change significantly when different mass-to-light ratio recipes were used. We also point out, for the first time, that the rising rotation curves as well as asymmetric rotation curves show a larger dispersion in the Tully-Fisher relation than the flat ones or the symmetric ones. Using the baryonic mass and the optical radius of galaxies, we found that the surface baryonic mass density is almost constant for all the galaxies of this sample. In this study we also emphasize the presence of a break in the NIR Tully-Fisher relation at M(H,K) similar to -20 and we confirm that late-type galaxies present higher total-to-baryonic mass ratios than early-type spirals, suggesting that supernova feedback is actually an important issue in late-type spirals. Due to the well-defined sample selection criteria and the homogeneity of the data analysis, the Tully-Fisher relation for GHASP galaxies can be used as a reference for the study of this relation in other environments and at higher redshifts.

Diagnostic tool to analyse colour-magnitude diagrams of poorly populated stellar concentrations

The dynamical processes that lead to open cluster disruption cause its mass to decrease. To investigate such processes from the observational point of view, it is important to identify open cluster remnants (OCRs), which are intrinsically poorly populated. Due to their nature, distinguishing them from field star fluctuations is still an unresolved issue. In this work, we developed a statistical diagnostic tool to distinguish poorly populated star concentrations from background field fluctuations. We use 2MASS photometry to explore one of the conditions required for a stellar group to be a physical group: to produce distinct sequences in a colour-magnitude diagram (CMD). We use automated tools to (i) derive the limiting radius; (ii) decontaminate the field and assign membership probabilities; (iii) fit isochrones; and (iv) compare object and field CMDs, considering the isochrone solution, in order to verify the similarity. If the object cannot be statistically considered as a field fluctuation, we derive its probable age, distance modulus, reddening and uncertainties in a self-consistent way. As a test, we apply the tool to open clusters and comparison fields. Finally, we study the OCR candidates DoDz 6, NGC 272, ESO 435 SC48 and ESO 325 SC15. The tool is optimized to treat these low-statistic objects and to separate the best OCR candidates for studies on kinematics and chemical composition. The study of the possible OCRs will certainly provide a deep understanding of OCR properties and constraints for theoretical models, including insights into the evolution of open clusters and dissolution rates.