Pierre Lesaffre

Bio: Pierre Lesaffre is an academic researcher from École Normale Supérieure. The author has contributed to research in topics: Physics & Turbulence. The author has an hindex of 12, co-authored 23 publications receiving 3360 citations. Previous affiliations of Pierre Lesaffre include University of Paris & Centre national de la recherche scientifique.

Modules for Experiments in Stellar Astrophysics (MESA)

Abstract: 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/).

Modules for Experiments in Stellar Astrophysics (MESA)

Abstract: 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 libraries for a wide range of applications in computational stellar astrophysics. A newly designed 1-D stellar evolution module, MESA star, 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. MESA star 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. Independently usable modules provide equation of state, opacity, nuclear reaction rates, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own public interface. Examples include comparisons to other codes and show evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets; the complete evolution of a 1 Msun star from the pre-main sequence to a cooling white dwarf; the Solar sound speed profile; the evolution of intermediate mass stars through the thermal pulses on the He-shell burning AGB phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; evolutionary tracks of massive stars from the pre-main sequence to the onset of core collapse; stars undergoing Roche lobe overflow; and accretion onto a neutron star. Instructions for downloading and installing MESA can be found on the project web site (this http URL).

Sulphur-bearing molecules in diffuse molecular clouds: new results from SOFIA/GREAT and the IRAM 30 m telescope

Abstract: We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the SH 1383 GHz 2 Π3/2 J = 5/2 ← 3/2 lambda doublet toward the star-forming regions W31C, G29.96–0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the EMIR receivers on the IRAM 30 m telescope at Pico Veleta to detect the H2 S 110 −101 (169 GHz), CS J = 2−1 (98 GHz) and SO 32 −21 (99 GHz) transitions. Upper limits on the H3 S+ 10 −00 (293 GHz) transition were also obtained at the IRAM 30 m. In nine foreground absorption components detected towards these sources, the inferred column densities of the four detected molecules showed relatively constant ratios, with N (SH) /N (H2 S) in the range 1.1−3.0, N (CS) /N (H2 S) in the range 0.32−0.61, and N (SO) /N (H2 S) in the range 0.08−0.30. The column densities of the sulphur-bearing molecules are very well correlated amongst themselves, moderately well correlated with CH (a surrogate tracer for H2 ), and poorly correlated with atomic hydrogen. The observed SH/H2 ratios – in the range 5 to 26 × 10-9 – indicate that SH (and other sulphur-bearing molecules) account for ≪ 1% of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the column densities of H2 S and SO by a factor of about 10; more elaborate shock models – in which account is taken of the velocity drift, relative to H2 , of SH molecules produced by the dissociative recombination of H3 S+ – reduce this discrepancy to a factor ~3.

Ortho-H2 and the age of prestellar cores

Abstract: 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

Abstract: 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.

A massive pulsar in a compact relativistic binary

Abstract: Many physically motivated extensions to general relativity (GR) predict significant deviations at energies present in massive neutron stars. We report the measurement of a 2.01 \(\pm \) 0.04 solar mass (M\(_\odot \)) pulsar in a 2.46-h orbit around a 0.172 \(\pm \) 0.003 M\(_\odot \) white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detection experiments. Additionally, the system strengthens recent constraints on the properties of dense matter and provides novel insight to binary stellar astrophysics and pulsar recycling.

Modules for Experiments in Stellar Astrophysics (MESA): Planets, Oscillations, Rotation, and Massive Stars

Abstract: We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA star. Improvements in MESA star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 M ? stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA star solves the fully coupled stellar structure and composition equations, and we show how this has improved the scaling of MESA's calculational speed on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit that packages all the required components needed to form a unified, maintained, and well-validated build environment for MESA. We also highlight a few tools developed by the community for rapid visualization of MESA star results.

Modules for experiments in stellar astrophysics (mesa): binaries, pulsations, and explosions

Abstract: We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

The CO-to-H2 Conversion Factor

Abstract: CO line emission represents the most accessible and widely used tracer of the molecular interstellar medium. This renders the translation of observed CO intensity into total H2 gas mass critical to understand star formation and the interstellar medium in our Galaxy and beyond. We review the theoretical underpinning, techniques, and results of efforts to estimate this CO-to-H2 “conversion factor,” XCO, in different environments. In the Milky Way disk, we recommend a conversion factor XCO = 2×10 20 cm −2 (K km s −1 ) −1 with ±30% uncertainty. Studies of other “normal galaxies” return similar values in Milky Way-like disks, but with greater scatter and systematic uncertainty. Departures from this Galactic conversion factor are both observed and expected. Dust-based determinations, theoretical arguments, and scaling relations all suggest that XCO increases with decreasing metallicity, turning up sharply below metallicity ≈ 1/3–1/2 solar in a manner consistent with model predictions that identify shielding as a key parameter. Based on spectral line modeling and dust observations, XCO appears to drop in the central, bright regions of some but not all galaxies, often coincident with regions of bright CO emission and high stellar surface density. This lower XCO is also present in the overwhelmingly molecular interstellar medium of starburst galaxies, where several lines of evidence point to a lower CO-to-H2 conversion factor. At high redshift, direct evidence regarding the conversion factor remains scarce; we review what is known based on dynamical modeling and other arguments. Subject headings: ISM: general — ISM: molecules — galaxies: ISM — radio lines: ISM

Solving Ordinary Differential Equations

Abstract: Initial-value ordinary differential equation solution via variable order Adams method (SIVA/DIVA) package is collection of subroutines for solution of nonstiff ordinary differential equations. There are versions for single-precision and double-precision arithmetic. Requires fewer evaluations of derivatives than other variable-order Adams predictor/ corrector methods. Option for direct integration of second-order equations makes integration of trajectory problems significantly more efficient. Written in FORTRAN 77.