Falk Herwig

Bio: Falk Herwig is an academic researcher from University of Victoria. The author has contributed to research in topics: Nucleosynthesis & Stellar evolution. The author has an hindex of 59, co-authored 287 publications receiving 14062 citations. Previous affiliations of Falk Herwig include Los Alamos National Laboratory & Keele University.

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

Evolution of Asymptotic Giant Branch Stars

Abstract: ▪ Abstract The current status of modeling the evolution and nucleosynthesis of asymptotic giant branch (AGB) stars is reviewed. The principles of AGB evolution have been investigated in recent years leading to improved and refined models, for example with regard to hot-bottom burning or the third dredge-up. The postprocessing s-process model yields quantitative results that reproduce many observations. However, these and most other processes in AGB stars are intimately related to the physics of stellar mixing. Mixing in AGB stars is currently not well-enough understood for accurate yield predictions. Several constraints and methods are available to improve the models. Some regimes of AGB evolution have not yet been studied in sufficient detail. These include the super-AGB stars and AGB stars at extremely low or ultra low metallicity.

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).

The Weak s-Process in Massive Stars and its Dependence on the Neutron Capture Cross Sections

Abstract: The slow neutron capture process in massive stars (weak s process) produces most of the s-process isotopes between iron and strontium. Neutrons are provided by the 22Ne(?,n)25Mg reaction, which is activated at the end of the convective He-burning core and in the subsequent convective C-burning shell. The s-process-rich material in the supernova ejecta carries the signature of these two phases. In the past years, new measurements of neutron capture cross sections of isotopes beyond iron significantly changed the predicted weak s-process distribution. The reason is that the variation of the Maxwellian-averaged cross sections (MACS) is propagated to heavier isotopes along the s path. In the light of these results, we present updated nucleosynthesis calculations for a 25 M ? star of Population I (solar metallicity) in convective He-burning core and convective C-burning shell conditions. In comparison with previous simulations based on the Bao et?al. compilation, the new measurement of neutron capture cross sections leads to an increase of s-process yields from nickel up to selenium. The variation of the cross section of one isotope along the s-process path is propagated to heavier isotopes, where the propagation efficiency is higher for low cross sections. New 74Ge, 75As, and 78Se MACS result in a higher production of germanium, arsenic, and selenium, thereby reducing the s-process yields of heavier elements by propagation. Results are reported for the He core and for the C shell. In shell C-burning, the s-process nucleosynthesis is more uncertain than in the He core, due to higher MACS uncertainties at higher temperatures. We also analyze the impact of using the new lower solar abundances for CNO isotopes on the s-process predictions, where CNO is the source of 22Ne, and we show that beyond Zn this is affecting the s-process yields more than nuclear or stellar model uncertainties considered in this paper. In particular, using the new updated initial composition, we obtain a high s-process production (overproduction higher than 16O, ~100) for Cu, Ga, Ge, and As. Using the older abundances by Anders & Grevesse, also Se, Br, Kr, and Rb are efficiently produced. Our results have important implications in explaining the origin of copper in the solar abundance distribution, pointing to a prevailing contribution from the weak s-process in agreement with spectroscopic observations and Galactic chemical evolution calculations. Because of the improvement due to the new MACS for nickel and copper isotopes, the nucleosynthesis of copper is less affected by nuclear uncertainties compared to heavier s-process elements. An experimental determination of the 63Ni MACS is required for a further improvement of the abundance prediction of copper. The available spectroscopic observations of germanium and gallium in stars are also discussed, where most of the cosmic abundances of these elements derives from the s-process in massive stars.

The supernova channel of super-AGB stars

Abstract: We study the late evolution of solar metallicity stars in the transition region between white dwarf formation and core collapse. This includes the super-asymptotic giant branch (super-AGB, SAGB) stars, which ignite carbon burning and form an oxygen-neon (ONe) core. SAGB star cores may grow to the Chandrasekhar mass because of continued H- and He-shell burning, ending as core-collapse supernovae. From stellar evolution models we find that the initial mass range for SAGB evolution is -->7.5–9.25 M☉. We perform calculations with three different stellar evolution codes to judge the robustness of our results. The mass range significantly depends on the treatment of semiconvective mixing and convective overshooting. To consider the effect of a large number of thermal pulses, as expected in SAGB stars, we construct synthetic SAGB models that are calibrated through stellar evolution simulations. The synthetic model enables us to compute the evolution of the main properties of SAGB stars from the onset of thermal pulses until the core reaches the Chandrasekhar mass or is uncovered by the stellar wind. Thereby, we differentiate the stellar initial mass ranges that produce ONe WDs from that leading to electron-capture SNe. The latter is found to be -->9.0–9.25 M☉ for our fiducial model, implying that electron-capture SNe would constitute about 4% of all SNe in the local universe. The error in this determination due to uncertainties in the third dredge-up efficiency and AGB mass-loss rate could lead to about a doubling of the number of electron-capture SNe, which provides a firm upper limit to their contribution to all supernovae of ~20%.

The Observation of Gravitational Waves from a Binary Black Hole Merger

Abstract: On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0×10(-21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ. The source lies at a luminosity distance of 410(-180)(+160) Mpc corresponding to a redshift z=0.09(-0.04)(+0.03). In the source frame, the initial black hole masses are 36(-4)(+5)M⊙ and 29(-4)(+4)M⊙, and the final black hole mass is 62(-4)(+4)M⊙, with 3.0(-0.5)(+0.5)M⊙c(2) radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.

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

Standard Solar Composition

Abstract: We review the current status of our knowledge of the chemical composition of the Sun, essentially derived from the analysis of the solar photospheric spectrum. The comparison of solar and meteoritic abundances confirms that there is a very good agreement between the two sets of abundances. They are used to construct a Standard Abundance Distribution.

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.