We discuss the cosmological simulation code GADGET-2, a new massively parallel TreeSPH code, capable of following a collisionless fluid with the N-body method, and an ideal gas by means of smoothed particle hydrodynamics (SPH). Our implementation of SPH manifestly conserves energy and entropy in regions free of dissipation, while allowing for fully adaptive smoothing lengths. Gravitational forces are computed with a hierarchical multipole expansion, which can optionally be applied in the form of a TreePM algorithm, where only short-range forces are computed with the ‘tree’ method while long-range forces are determined with Fourier techniques. Time integration is based on a quasi-symplectic scheme where long-range and short-range forces can be integrated with different time-steps. Individual and adaptive short-range time-steps may also be employed. The domain decomposition used in the parallelization algorithm is based on a space-filling curve, resulting in high flexibility and tree force errors that do not depend on the way the domains are cut. The code is efficient in terms of memory consumption and required communication bandwidth. It has been used to compute the first cosmological N-body simulation with more than 10 10 dark matter particles, reaching a homogeneous spatial dynamic range of 10 5 per dimension in a three-dimensional box. It has also been used to carry out very large cosmological SPH simulations that account for radiative cooling and star formation, reaching total particle numbers of more than 250 million. We present the algorithms used by the code and discuss their accuracy and performance using a number of test problems. GADGET-2 is publicly released to the research community. Ke yw ords: methods: numerical ‐ galaxies: interactions ‐ dark matter.

/pdf/the-cosmological-simulation-code-gadget-2-48ba4eyb8k.pdf

The Cosmological simulation code GADGET-2

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.

https://iopscience.iop.org/article/10.1088/0067-0049/220/1/15/pdf

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

amount of energy, a tiny fraction of which is sufficient to explode the star as a supernova. The authors examine our current understanding of the lives and deaths of massive stars, with special attention to the relevant nuclear and stellar physics. Emphasis is placed upon their post-helium-burning evolution. Current views regarding the supernova explosion mechanism are reviewed, and the hydrodynamics of supernova shock propagation and ‘‘fallback’’ is discussed. The calculated neutron star masses, supernova light curves, and spectra from these model stars are shown to be consistent with observations. During all phases, particular attention is paid to the nucleosynthesis of heavy elements. Such stars are capable of producing, with few exceptions, the isotopes between mass 16 and 88 as well as a large fraction of still heavier elements made by the r and p processes.

/pdf/the-evolution-and-explosion-of-massive-stars-57jc6k6s1i.pdf

The evolution and explosion of massive stars

A new N-body and hydrodynamical code, called RAMSES, is presented. It has been designed to study structure formation in the universe with high spatial resolution. The code is based on Adaptive Mesh Renement (AMR) technique, with a tree-based data structure allowing recursive grid renements on a cell-by-cell basis. The N-body solver is very similar to the one developed for the ART code (Kravtsov et al. 1997), with minor dierences in the exact implementation. The hydrodynamical solver is based on a second-order Godunov method, a modern shock-capturing scheme known to compute accurately the thermal history of the fluid component. The accuracy of the code is carefully estimated using various test cases, from pure gas dynamical tests to cosmological ones. The specic renement strategy used in cosmological simulations is described, and potential spurious eects associated with shock waves propagation in the resulting AMR grid are discussed and found to be negligible. Results obtained in a large N-body and hydrodynamical simulation of structure formation in a low density CDM universe are reported, with 256 3 particles and 4:1 10 7 cells in the AMR grid, reaching a formal resolution of 8192 3 .A convergence analysis of dierent quantities, such as dark matter density power spectrum, gas pressure power spectrum and individual haloe temperature proles, shows that numerical results are converging down to the actual resolution limit of the code, and are well reproduced by recent analytical predictions in the framework of the halo model.

/pdf/cosmological-hydrodynamics-with-adaptive-mesh-refinement-a-2696ypxntf.pdf

Cosmological hydrodynamics with adaptive mesh refinement - A new high resolution code called RAMSES

Aims. Many topical astrophysical research areas, such as the properties of planet host stars, the nature of the progenitors of different types of supernovae and gamma ray bursts, and the evolution of galaxies, require complete and homogeneous sets of stellar models at different metallicities in order to be studied during the whole of cosmic history. We present here a first set of models for solar metallicity, where the effects of rotation are accounted for in a homogeneous way.Methods. We computed a grid of 48 different stellar evolutionary tracks, both rotating and non-rotating, at Z  = 0.014, spanning a wide mass range from 0.8 to 120 M ⊙ . For each of the stellar masses considered, electronic tables provide data for 400 stages along the evolutionary track and at each stage, a set of 43 physical data are given. These grids thus provide an extensive and detailed data basis for comparisons with the observations. The rotating models start on the zero-age main sequence (ZAMS) with a rotation rate υ ini /υ crit  = 0.4. The evolution is computed until the end of the central carbon-burning phase, the early asymptotic giant branch (AGB) phase, or the core helium-flash for, respectively, the massive, intermediate, and both low and very low mass stars. The initial abundances are those deduced by Asplund and collaborators, which best fit the observed abundances of massive stars in the solar neighbourhood. We update both the opacities and nuclear reaction rates, and introduce new prescriptions for the mass-loss rates as stars approach the Eddington and/or the critical velocity. We account for both atomic diffusion and magnetic braking in our low-mass star models.Results. The present rotating models provide a good description of the average evolution of non-interacting stars. In particular, they reproduce the observed main-sequence width, the positions of the red giant and supergiant stars in the Hertzsprung-Russell (HR) diagram, the observed surface compositions and rotational velocities. Very interestingly, the enhancement of the mass loss during the red-supergiant stage, when the luminosity becomes supra-Eddington in some outer layers, help models above 15−20 M ⊙ to lose a significant part of their hydrogen envelope and evolve back into the blue part of the HR diagram. This result has interesting consequences for the blue to red supergiant ratio, the minimum mass for stars to become Wolf-Rayet stars, and the maximum initial mass of stars that explode as type II−P supernovae.

/pdf/grids-of-stellar-models-with-rotation-i-models-from-0-8-to-3bjuk0gnjq.pdf

Grids of stellar models with rotation - I. Models from 0.8 to 120 M⊙ at solar metallicity (Z = 0.014)

Predictive r-process calculations

PARAMESH: A Parallel, Adaptive Mesh Refinement Toolkit and Performance of the ASCI/FLASH code

FLASH: A Multidimensional Hydrodynamics Code for Modeling Astrophysical Thermonuclear Flashes

Computational skills are required across all astronomy disciplines Many students enter degree programs without sufficient skills to solve computational problems in their core classes or contribute immediately to research We recommend advocacy for computational literacy, familiarity with fundamental software carpentry skills, and mastery of basic numerical methods by the completion of an undergraduate degree in Astronomy 
We recommend the AAS Education Task Force advocate for a significant increase in computational literacy 
We encourage the AAS to modestly fund efforts aimed at providing Open Education Resources (OER) that will significantly impact computational literacy in astronomy education

Frank Timmes

Papers

Predictive r-process calculations

PARAMESH: A Parallel, Adaptive Mesh Refinement Toolkit and Performance of the ASCI/FLASH code

FLASH: A Multidimensional Hydrodynamics Code for Modeling Astrophysical Thermonuclear Flashes

The Importance of Computation in Astronomy Education

Tracing the Cosmic Star Formation History to its Beginnings: GRBs as Tools