Pierre Lesaffre

MODULES FOR EXPERIMENTS IN STELLAR ASTROPHYSICS (MESA)

Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source, robust, efficient, thread-safe libraries for a wide range of applications in computational stellar astrophysics. A one-dimensional stellar evolution module, MESAstar, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very low mass to massive stars, including advanced evolutionary phases. MESAstar solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. State-of-the-art modules provide equation of state, opacity, nuclear reaction rates, element diffusion data, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own explicitly defined public interface to facilitate independent development. Several detailed examples indicate the extensive verification and testing that is continuously performed and demonstrate the wide range of capabilities that MESA possesses. These examples include evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets to very old ages; the complete evolutionary track of a 1 M ☉ star from the pre-main sequence (PMS) to a cooling white dwarf; the solar sound speed profile; the evolution of intermediate-mass stars through the He-core burning phase and thermal pulses on the He-shell burning asymptotic giant branch phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; the complete evolutionary tracks of massive stars from the PMS to the onset of core collapse; mass transfer from stars undergoing Roche lobe overflow; and the evolution of helium accretion onto a neutron star. MESA can be downloaded from the project Web site (http://mesa.sourceforge.net/).

ORTHO-H2 AND THE AGE OF INTERSTELLAR DARK CLOUDS

Interstellar dark clouds are the sites of star formation. Their main component, dihydrogen, exists under two states, ortho and para. H2 is supposed to form in the ortho:para ratio (OPR) of 3:1 and to subsequently decay to almost pure para-H2 (OPR ≤ 0.001). Only if the H2 OPR is low enough, will deuteration enrichment, as observed in the cores of these clouds, be efficient. The second condition for strong deuteration enrichment is the local disappearance of CO, which freezes out onto grains in the core formation process. We show that this latter condition does not apply to DCO+, which, therefore, should be present all over the cloud. We find that an OPR ≥ 0.1 is necessary to prevent DCO+ large-scale apparition. We conclude that the inevitable decay of ortho-H2 sets an upper limit of ~6 million years to the age of starless molecular clouds under usual conditions.

MRI channel flows and their parasites

Local simulations of the magnetorotational instability (MRI) in accretion discs can exhibit recurrent coherent structures called channel flows. The formation and destruction of these structures may play a role in the development and saturation of MRI-induced turbulence, and consequently help us understand the time-dependent accretion behaviour of certain astrophysical objects. Previous investigations have revealed that channel solutions are attacked by various parasitic modes, foremost of which is an analogue of the Kelvin-Helmholtz instability. We revisit these instabilities and show how they relate to the classical instabilities of plasma physics, the kink and pinch modes. However, we argue that in most cases channels emerge from developed turbulence and are eventually destroyed by turbulent mixing, not by the parasites. The exceptions are clean isolated channels, which appear in systems near criticality or which emerge from low amplitude initial conditions. These structures can achieve large amplitudes and are only then destroyed, giving rise to eruptive behaviour.

Strong CH+ J = 1-0 emission and absorption in DR21

We report the first detection of the ground-state rotational transition of the methylidyne cation CH^+ towards the massive star-forming region DR 21 with the HIFI instrument onboard the Herschel satellite. The line profile exhibits a broad emission line, in addition to two deep and broad absorption features associated with the DR 21 molecular ridge and foreground gas. These observations allow us to determine a ^(12)CH^(+)J = 1–0 line frequency of ν = 835 137 ± 3 MHz, in good agreement with a recent experimental determination. We estimate the CH^+ column density to be a few 10^(13) cm^(-2) in the gas seen in emission, and >10^(14) cm^(-2) in the components responsible for the absorption, which is indicative of a high line of sight average abundance [CH^+] /[H] > 1.2 × 10^(-8). We show that the CH^+ column densities agree well with the predictions of state-of-the-art C-shock models in dense UV-illuminated gas for the emission line, and with those of turbulent dissipation models in diffuse gas for the absorption lines.

Ortho-H2 and the age of prestellar cores

Prestellar cores form from the contraction of cold gas and dust material in dark clouds before they collapse to form protostars. Several concurrent theories exist to describe this contraction but they are currently difficult to distinguish. One major difference is the timescale involved in forming the prestellar cores: some theories advocate nearly free-fall speed via, e.g., rapid turbulence decay, while others can accommodate much longer periods to let the gas accumulate via, e.g., ambipolar diffusion. To tell the difference between these theories, measuring the age of prestellar cores could greatly help. However, no reliable clock currently exists. We present a simple chemical clock based on the regulation of the deuteration by the abundance of ortho‐H2 that slowly decays away from the ortho-para statistical ratio of 3 down to or less than 0.001. We use a chemical network fully coupled to a hydrodynamical model that follows the contraction of a cloud, starting from uniform density, and reaches a density profile typical of a prestellar core. We compute the N2D + /N2H + ratioalong the density profile. The disappearance of ortho-H2 is tied tothe duration of the contraction and the N2D + /N2H + ratio increases in the wake of the ortho-H2 abundance decrease. By adjusting the time of contraction, we obtain different deuteration profiles that we can compare to the observations. Our model can test fast contractions (from 10 4 to 10 6 cm −3 in ∼0.5 My) and slow contractions (from 10 4 to 10 6 cm −3 in ∼5 My). We have tested the sensitivity of the models to various initial conditions. The slowcontraction deuteration profile is approximately insensitive to these variations, while the fast-contraction deuteration profile shows significant variations. We found that, in all cases, the deuteration profile remains clearly distinguishable whether it comes from the fast collapse or the slow collapse. We also study the para-D2H + /ortho-H2D + ratio and find that its variation is not monotonic, so it does not discriminate between models. Applying this model to L183 (=L134N), we find that the N2D + /N2H + ratio would be higher than unity for evolutionary timescales of a few megayears independently of other parameters, such as cosmic ray ionization rate or grain size (within reasonable ranges). A good fit to the observations is only obtained for fast contraction (≤0.7 My from the beginning of the contraction and ≤4 My from the birth of the molecular cloud based on the need to keep a high ortho-H2 abundance when the contraction starts ‐ ortho-H2/para-H2 ≥ 0.2 ‐ to match the observations). This chemical clock therefore rules out slow contraction in L183 and steady-state chemical models, since steady state is clearly not reached here. This clock should be applied to other cores to help distinguish slow and fast contraction theories over a large sample of cases.

UV-driven chemistry in simulations of the interstellar medium I. Post-processed chemistry with the Meudon PDR code

Context. Observations have long demonstrated the molecular diversity of the di use interstellar medium (ISM). Only now, with the advent of high-performance computing, does it become possible for numerical simulations of astrophysical fluids to include a treatment of chemistry, in order to faithfully reproduce the abundances of the many observed species, and especially that of CO, which is used as a proxy for molecular hydrogen. When applying photon-dominated region (PDR) codes to describe the UV-driven chemistry of uniform density cloud models, it is found that the observed abundances of CO are not well reproduced. Aims. Our main purpose is to estimate the e ect of assuming uniform density on the line-of-sight in PDR chemistry models, compared to a more realistic distribution for which total gas densities may well vary by several orders of magnitude. A secondary goal of this paper is to estimate the amount of molecular hydrogen which is not properly traced by the CO (J = 1! 0) line, the so-called ”dark molecular gas”. Methods. We use results from a magnetohydrodynamical (MHD) simulation as a model for the density structures found in a turbulent di use ISM with no star-formation activity. The Meudon PDR code is then applied to a number of lines of sight through this model, to derive their chemical structures. Results. It is found that, compared to the uniform density assumption, maximal chemical abundances for H2, CO, CH and CN are increased by a factor 2 4 when taking into account density fluctuations on the line of sight. The correlations between column densities of CO, CH and CN with respect to those of H2 are also found to be in better overall agreement with observations. For instance, at N(H2) & 2 10 20 cm 2 , while observations suggest that d[logN(CO)]=d[logN(H2)] ’ 3:07 0:73, we find d[logN(CO)]=d[logN(H2)] ’ 14 when assuming uniform density, and d[logN(CO)]=d[logN(H2)] ’ 5:2 when including density fluctuations.

H2, H3+ and the age of molecular clouds and prestellar cores

Measuring the age of molecular clouds and prestellar cores is a difficult task that has not yet been successfully accomplished although the information is of paramount importance to help in understanding and discriminating between different formation scenarios. Most chemical clocks suffer from unknown initial conditions and are therefore difficult to use. We propose a new approach based on a subset of deuterium chemistry that takes place in the gas phase and for which initial conditions are relatively well known. It relies primarily on the conversion of H3+ into H2D+ to initiate deuterium enrichment of the molecular gas. This conversion is controlled by the ortho/para ratio of H2 that is thought to be produced with the statistical ratio of 3 and subsequently slowly decays to an almost pure para-H2 phase. This slow decay takes approximately 1 Myr and allows us to set an upper limit on the age of molecular clouds. The deuterium enrichment of the core takes longer to reach equilibrium and allows us to estimate the time necessary to form a dense prestellar core, i.e. the last step before the collapse of the core into a protostar. We find that the observed abundance and distribution of DCO+ and N2D+ argue against quasi-static core formation and favour dynamical formation on time scales of less than 1 Myr. Another consequence is that ortho-H2 remains comparable to para-H2 in abundance outside the dense cores.

Tracing extended low-velocity shocks through SiO emission case study of the W43-MM1 ridge

Aims. Previous literature suggests that the densest structures in the interstellar medium form through colliding flows, but patent evidence of this process is still missing. Recent literature proposes using SiO line emission to trace low-velocity shocks associated with cloud formation through collision. In this paper we investigate the bright and extended SiO(2−1) emission observed along the ∼5 pc-long W43-MM1 ridge to determine its origin. Methods. We used high angular resolution images of the SiO(2−1) and HCN(1−0) emission lines obtained with the IRAM plateau de Bure (PdBI) interferometer and combined with data from the IRAM 30 m radiotelescope. These data were complemented by a Herschel column density map of the region. We performed spectral analysis of SiO and HCN emission line profiles to identify protostellar outflows and spatially disentangle two velocity components associated with low- and high-velocity shocks. Then, we compared the low-velocity shock component to a dedicated grid of one-dimensional (1D) radiative shock models. Results. We find that the SiO emission originates from a mixture of high-velocity shocks caused by bipolar outflows and lowvelocity shocks. Using SiO and HCN emission lines, we extract seven bipolar outflows associated with massive dense cores previously identified within the W43-MM1 mini-starburst cluster. Comparing observations with dedicated Paris-Durham shock models constrains the velocity of the low-velocity shock component from 7 to 12 km s−1 . Conclusions. The SiO arising from low-velocity shocks spreads along the complete length of the ridge. Its contribution represents at least 45% and up to 100% of the total SiO emission depending on the area considered. The low-velocity component of SiO is most likely associated with the ridge formation through colliding flows or cloud-cloud collision.

Turbulent energy dissipation and intermittency in ambipolar diffusion magnetohydrodynamics

The dissipation of kinetic and magnetic energy in the interstellar medium (ISM) can proceed through viscous, Ohmic or ambipolar diffusion (AD). It occurs at very small scales compared to the scales at which energy is presumed to be injected. This localized heating may impact the ISM evolution but also its chemistry, thus providing observable features. Here, we perform 3D spectral simulations of decaying magnetohydrodynamic turbulence including the effects of AD. We find that the AD heating power spectrum peaks at scales in the inertial range, due to a strong alignment of the magnetic and current vectors in the dissipative range. AD affects much greater scales than the AD scale predicted by dimensional analysis. We find that energy dissipation is highly concentrated on thin sheets. Its probability density function follows a lognormal law with a power-law tail which hints at intermittency, a property which we quantify by use of structure function exponents. Finally, we extract structures of high dissipation, defined as connected sets of points where the total dissipation is most intense and we measure the scaling exponents of their geometric and dynamical characteristics: the inclusion of AD favours small sizes in the dissipative range.

Uncertainties and robustness of the ignition process in type Ia supernovae

Context. It is widely accepted that the onset of the explosive carbon burning in the core of a carbon-oxygen white dwarf (CO WD) triggers the ignition of a type Ia supernova (SN Ia). The features of the ignition are among the few free parameters of the SN Ia explosion theory. Aims. We explore the role for the ignition process of two different issues: firstly, the ignition is studied in WD models coming from different accretion histories. Secondly, we estimate how a different reaction rate for C-burning can affect the ignition. Methods. Two-dimensional hydrodynamical simulations of temperature perturbations in the WD core (“bubbles”) are performed with the FLASH code. In order to evaluate the impact of the C-burning reaction rate on the WD model, the evolution code FLASH_THE_TORTOISE from Lesaffre et al. (2006, MNRAS, 368, 187) is used. Results. In different WD models a key role is played by the different gravitational acceleration in the progenitor’s core. As a consequence, the ignition is disfavored at a large distance from the WD center in models with a larger central density, resulting from the evolution of initially more massive progenitors. Changes in the C reaction rate at T < 5 × 10 8 K slightly influence the ignition density in the WD core, while the ignition temperature is almost unaffected. Recent measurements of new resonances in the C-burning reaction rate (Spillane et al. 2007, Phys. Rev. Lett., 98, 122501) do not affect the core conditions of the WD significantly. Conclusions. This simple analysis, performed on the features of the temperature perturbations in the WD core, should be extended in the framework of the state-of-the-art numerical tools for studying the turbulent convection and ignition in the WD core. Future measurements of the C-burning reactions cross section at low energy, though certainly useful, are not expected to affect our current understanding of the ignition process dramatically.