Astrophysical Simulations Research Articles

Analysis the Principle and Realization for Cosmic Dust Simulation

Astrophysics simulations are becoming crucial to the understanding of how the universe came to be the way it is today. Cosmic dusts are quite important in these astrophysics’ simulations, although they have only recently been introduced to the simulations. This study offers analysis on the different parts of current cosmic dust simulations. Astrophysics simulations and cosmic dusts’ role in this simulation was first discussed. Simulation requirements on cosmic dusts were then evaluated, both in the software and the hardware. The realization of cosmic dusts simulations was then discussed, focusing on the creation, growth, and destruction of cosmic dusts and briefly discussing the algorithmic process of incorporating calculations and presenting data. At last, limitations on current cosmic dusts simulations and some future prospects were proposed. These results provide an organization of the different prospects of current progresses in cosmic dust simulations, which could be beneficial as an overview of the topic and a base for which directions for future researches could be suggested.

Comparison and Applications for the State-of-art Astrophysics Simulations: Evidence from THESAN and illustrisTNG

The understandings of the formation of galaxy have been remarkably improved into a new level due to the simulations in astrophysics. Early models using simple N-body simulations were limited due to computational power. Later on, the development of algorithms had enhanced accuracy and scalability. In the modern era (2000s-present), due to the high-resolution simulation supported by projects including the MillenniumTNG and THESAN significantly enhanced the detailed studies of the formation and evolution of galaxy. This study compares the differences of application and comparison between the simulations of the two projects, the MillenniumTNG and THESAN, specific in the concept of high-redshift. While the two projects dive into deep using different methods and taking different angles, they share the same overarching goals of having further understanding of the early universe and the high-redshift galaxies. In particular, the two projects complement each other through their simulations and observations, pursuing better insights into the cosmos at large.

Cosmological Radiation Simulation Analysis: Taking Arepo-RT as an Example

With the fast-growing progress in computational physics, astrophysical simulations are made possible and continually upgraded with novel physical models. Contemporarily, many modern simulations feature the capability of cosmological radiation. In order to collate the recent improvements in this field, the numerical solver Arepo-RT is discussed as an important example. This study reorganizes and interprets the main physical models and numerical techniques used in Arepo-RT and similar solvers, with some details simplified to suit common readers. This paper describes simply and briefly the finite volume method, the radiative transfer model, the radiation source model, the Godunov-type schemes and the time integration, all of which are frequently used in modern radiation simulations. Apart from the realizations of the solvers, the results and limitations are also outlined in this paper, including the justification of reduced speed of light approximation, the results compared with JWST, and the issue of missing photons. These results shed light on guiding further improvements for Arepo-RT.

Comparing Performance in the Codes AREPO, GIZMO and GAMER-2

Modern astrophysical simulation is the main aspect of human getting insight from the cosmos. Numerical simulations became irreplaceable tools for generating various phases of the universe, the formation of galaxies, and other astrophysical phenomena. As the advances of computational power, validity and diversity of astrophysical modeling have been growing. In addition to the development of the processor performance, the diversity of the software community is also rising therefore, research between various simulation codes is worthwhile. It allowed researchers to know the nature of different simulation techniques linked to their basic algorithms. This study summarizes the basics of common simulation methods used by researchers including N-body Simulations, Smoothed Particle Hydrodynamics, Monte Carlo Methods, Finite Difference and Finite Element Methods. Furthermore, this research analyzes the basics of AREPO, GIZMO, and GAMER-2 codes, and compared the performance of these codes. Last, according to the analysis, this study discusses and demonstrates the advantages and limitations of these codes. These results shed light on guiding further exploration of cosmology simulations.

Numerical studies of hydrodynamic relativistic astrophysical jet: The modification of vorticity transport equation

In this Letter, we have carried out the two-dimensional numerical simulation of axisymmetric relativistic jet in cylindrical coordinates by employing higher order finite volume method in PLUTO [Mignone et al., “PLUTO: A numerical code for computational astrophysics,” Astrophys. J. Suppl. Ser. 170, 228 (2007)] solver. The modified vorticity transport equation has been proposed for relativistic flow by taking the curl of the momentum equation, which shows significant change in the baroclinic vorticity production term due to relativistic effect. Both mathematical analysis and numerical results show that the vorticity production term due to baroclinic torque is heavily influenced due to the presence of specific enthalpy gradient and square of Lorentz factor gradient in a relativistic fluid flow.

AHKASH: a new Hybrid particle-in-cell code for simulations of astrophysical collisionless plasma.

We introduce Astrophysical Hybrid-Kinetic simulations with the flash code ([Formula: see text]) - a new Hybrid particle-in-cell (PIC) code developed within the framework of the multiphysics code flash. The new code uses a second-order accurate Boris integrator and a predictor-predictor-corrector algorithm for advancing the Hybrid-kinetic equations, using the constraint transport method to ensure that magnetic fields are divergence-free. The code supports various interpolation schemes between the particles and grid cells, with post-interpolation smoothing to reduce finite particle noise. We further implement a [Formula: see text] method to study instabilities in weakly collisional plasmas. The new code is tested on standard physical problems such as the motion of charged particles in uniform and spatially varying magnetic fields, the propagation of Alfvén and whistler waves, and Landau damping of ion acoustic waves. We test different interpolation kernels and demonstrate the necessity of performing post-interpolation smoothing. We couple the turbgen turbulence driving module to the new Hybrid PIC code, allowing us to test the code on the highly complex physical problem of the turbulent dynamo. To investigate steady-state turbulence with a fixed sonic Mach number, it is important to maintain isothermal plasma conditions. Therefore, we introduce a novel cooling method for Hybrid PIC codes and provide tests and calibrations of this method to keep the plasma isothermal. We describe and test the 'hybrid precision' method, which significantly reduces (by a factor [Formula: see text]) the computational cost, without compromising the accuracy of the numerical solutions. Finally, we test the parallel scalability of the new code, showing excellent scaling up to 10,000 cores.

Tetrahedral grids in Monte Carlo radiative transfer

Context. To understand the structures of complex astrophysical objects, 3D numerical simulations of radiative transfer processes are invaluable. For Monte Carlo radiative transfer, the most common radiative transfer method in 3D, the design of a spatial grid is important and non-trivial. Common choices include hierarchical octree and unstructured Voronoi grids, each of which has advantages and limitations. Tetrahedral grids, commonly used in ray-tracing computer graphics, can be an interesting alternative option. Aims. We aim to investigate the possibilities, advantages, and limitations of tetrahedral grids in the context of Monte Carlo radiative transfer. In particular, we want to compare the performance of tetrahedral grids to other commonly used grid structures. Methods. We implemented a tetrahedral grid structure, based on the open-source library TetGen, in the generic Monte Carlo radiative transfer code SKIRT. Tetrahedral grids can be imported from external applications or they can be constructed and adaptively refined within SKIRT. We implemented an efficient grid traversal method based on Plücker coordinates and Plücker products. Results. The correct implementation of the tetrahedral grid construction and the grid traversal algorithm in SKIRT were validated using 2D radiative transfer benchmark problems. Using a simple 3D model, we compared the performance of tetrahedral, octree, and Voronoi grids. With a constant cell count, the octree grid outperforms the tetrahedral and Voronoi grids in terms of traversal speed, whereas the tetrahedral grid is poorer than the other grids in terms of grid quality. All told, we find that the performance of tetrahedral grids is relatively poor compared to octree and Voronoi grids. Conclusions. Although the adaptively constructed tetrahedral grids might not be favourable in most media representative of astrophysical simulation models, they still form an interesting unstructured alternative to Voronoi grids for specific applications. In particular, they might prove useful for radiative transfer post-processing of hydrodynamical simulations run on tetrahedral or unstructured grids.

Dissipated correction map method with trapezoidal rule for the simulations of gravitational waves from spinning compact binary

ABSTRACT The correction map method means extended phase-space algorithm with correction map. In our research, we have developed a correction map method, specifically the dissipated correction map method with trapezoidal rule, for numerical simulations of gravitational waves from spinning compact binary systems. This new correction map method, denoted as $CM3$, has shown remarkable performance in various simulation results, such as phase-space distance, dissipated energy error, and gravitational waveform, closely resembling the high-order precision implicit Gaussian algorithm. When compared with the previously used mid-point map which is denoted as $C_2$, the $CM3$ consistently exhibits a closer alignment with the highly accurate Gaussian algorithm in waveform evolution and orbital trajectory analysis. Through detailed comparisons and analyses, it is evident that $CM3$ outperforms other algorithms, including $CM2$ and $C_2$ mentioned in this paper, in terms of accuracy and precision in simulating spinning compact binary systems. The incorporation of the trapezoidal rule and the optimization with a scale factor $\gamma$ have significantly enhanced the performance of $CM3$, making it a promising method for future numerical simulations in astrophysics. With the groundbreaking detection of gravitational waves by the LIGO/VIRGO collaboration, interest in this research domain has soared. Our work contributes valuable insights for the application of matched filtering techniques in the analysis of gravitational wave signals, enhancing the precision and reliability of these detection.

Consistent Eccentricities for Gravitational-wave Astronomy: Resolving Discrepancies between Astrophysical Simulations and Waveform Models

Detecting imprints of orbital eccentricity in gravitational-wave (GW) signals promises to shed light on the formation mechanisms of binary black holes. To constrain the formation mechanisms, distributions of eccentricity derived from numerical simulations of astrophysical formation channels are compared to the estimates of eccentricity inferred from GW signals. We report that the definition of eccentricity typically used in astrophysical simulations is inconsistent with the one used while modeling GW signals, with the differences mainly arising due to the choice of reference frequency used in both cases. We also posit a prescription for calculating eccentricity from astrophysical simulations, by evolving ordinary differential equations obtained from post-Newtonian theory and using the dominant (ℓ = m = 2) mode’s frequency as the reference frequency; this ensures consistency in the definitions. On comparing the existing eccentricities of the binaries present in the Cluster Monte Carlo catalog of globular cluster simulations with the eccentricities calculated using the prescription presented here, we find a significant discrepancy at e ≳ 0.2; this discrepancy becomes worse with increasing eccentricity. We note the implications this discrepancy has for existing studies and recommend that care be taken when comparing data-driven constraints on eccentricity to expectations from astrophysical formation channels.

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10−3 to 10−1, which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

Emulating the interstellar medium chemistry with neural operators

Context. The study of galaxy formation and evolution critically depends on our understanding of the complex photo-chemical processes that govern the evolution and thermodynamics of the interstellar medium (ISM). In a computational sense, resolving the chemistry is among the weightiest tasks in cosmological and astrophysical simulations. Aims. Astrophysical simulations can include photo-chemical models that allow for a wide range of densities (n), abundances of different species (ni/n) and temperature (T), and plausible evolution scenarios of the ISM under the action of a radiation field (F) with different spectral shapes and intensities. The evolution of such a non-equilibrium photo-chemical network relies on implicit, precise, computationally costly, ordinary differential equations (ODE) solvers. Here, we aim to substitute such procedural solvers with fast, pre-trained emulators based on neural operators. Methods. We emulated a non-equilibrium chemical network up to H2 formation (9 species, 52 reactions) by adopting the DeepONet formalism, namely: by splitting the ODE solver operator that maps the initial conditions and time evolution into a tensor product of two neural networks (named branch and trunk). We used KROME to generate a training set, spanning −2 < log(n/cm−3) ≤ 3.5, log(20) ≤ log(T/K) ≤ 5.5, −6 ≤ log(ni/n) < 0, and adopting an incident radiation field, F, sampled in 10 energy bins with a continuity prior. We separately trained the solver for T and each ni for ≃4.34 GPUhrs. Results. Compared with the reference solutions obtained by KROME for single-zone models, the typical precision obtained is of the order of 10−2, that is, it is 10 times better when using a training that is 40 times less costly, with respect to previous emulators that only considered a fixed F. DeepONet also performs well for T and ni outside the range of the training sample. Furthermore, the emulator aptly reproduces the ion and temperature profiles of photo dissociation regions as well; namely, by giving errors that are comparable to the typical difference between various photo-ionization codes. The present model achieves a speed-up of a factor of 128× with respect to stiff ODE solvers. Conclusions. Our neural emulator represents a significant leap forward in the modelling of ISM chemistry, offering a good balance of precision, versatility, and computational efficiency. Nevertheless, further work is required to address the challenges represented by the extrapolation beyond the training time domain and the removal of potential outliers.

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth-order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

Challenges Modeling the Low-Luminosity Type Iax Supernovae

Numerical models allow the investigation of phenomena that cannot exist in a laboratory. Computational simulations are therefore essential for advancing our knowledge of astrophysics, however, the very nature of simulation requires making assumptions that can substantially affect their outcome. Here, we present the challenges faced when simulating dim thermonuclear explosions, Type Iax supernovae. This class of dim events produce a slow moving, sparse ejecta that presents challenges for simulation. We investigate the limitations of the equation of state and its applicability to the expanding, cooling ejecta. We also discuss how the“fluff”, i.e. the low-density gas on the grid in lieu of vacuum, inhibits the ejecta as it expands. We explore how the final state of the simulation changes as we vary the character of the burning, which influences the outcome of the explosion. These challenges are applicable to a wide range of astrophysical simulations, and are important to discuss and overcome as a community.

A Framework for Exploring Nuclear Physics Sensitivity in Numerical Simulations

We describe the AMReX-Astrophysics framework for exploring the sensitivity of astrophysical simulations to the details of a nuclear reaction network, including the number of nuclei, choice of reaction rates, and approximations used. This is explored by modeling a simple detonation with the Castro simulation code. The entire simulation methodology is open-source and GPU-enabled.

Ookami: An A64FX Computing Resource

We present a look at Ookami, a project providing community access to a testbed supercomputer with the ARM-based A64FX processors developed by a collaboration between RIKEN and Fujitsu and deployed in the Japanese supercomputer Fugaku. We provide an overview of the project and details of the hardware, and describe the user base and education/training program. We present highlights from previous performance studies of two astrophysical simulation codes and present a strong scaling study of a full 3D supernova simulation as an example of the the machine’s capability.

A new capacitive inductive system design for LASER-induced kilotesla magnetic field generation

This research focuses on exploring the nanosecond laser-driven coil systems capable of generating kT magnetic fields and the diverse applications of this system. Through investigating the effects of laser parameters and coil structures, the aim of this study is to unveil the physics of these generated intense magnetic fields. The outcomes gained from this research give an important and fundamental understanding on high magnetic field production, informing the development in laser-driven systems. The implications of this study extend to plasma physics, astrophysics simulations and fusion research. Furthermore, the study explains the advantages and applications of these intense magnetic fields and includes measurements of laser pulse powers according to coil materials.

Constraining the Onset Density for the QCD Phase Transition with the Neutrino Signal from Core-collapse Supernovae

The occurrence of a first-order hadron–quark matter phase transition at high baryon densities is investigated in astrophysical simulations of core-collapse supernovae, to decipher yet incompletely understood properties of the dense matter equation of state (EOS) using neutrinos from such cosmic events. It is found that the emission of a nonstandard second neutrino burst, dominated by electron antineutrinos, is not only a measurable signal for the appearance of deconfined quark matter but also reveals information about the state of matter at extreme conditions encountered at the supernova (SN) interior. To this end, a large set of spherically symmetric SN models is investigated, studying the dependence on the EOS and the stellar progenitor. General relativistic neutrino-radiation hydrodynamics is employed featuring three-flavor Boltzmann neutrino transport and a microscopic hadron-quark hybrid matter EOS class. Therefore, the DD2 relativistic mean-field hadronic model is employed, and several variations of it, and the string-flip model for the description of deconfined quark matter. The resulting hybrid model covers a representative range of onset densities for the phase transition and latent heats. This facilitates the direct connection between intrinsic signatures of the neutrino signal and properties of the EOS. In particular, a set of linear relations has been found empirically. These potentially provide a constraint for the onset density of a possible QCD phase transition from the future neutrino observation of the next galactic core-collapse SN, if a millisecond electron anti-neutrino burst is present around or less than 1 s.

The AGORA High-resolution Galaxy Simulations Comparison Project. V. Satellite Galaxy Populations in a Cosmological Zoom-in Simulation of a Milky Way–Mass Halo

We analyze and compare the satellite halo populations at z ∼ 2 in the high-resolution cosmological zoom-in simulations of a 1012 M ⊙ target halo (z = 0 mass) carried out on eight widely used astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Arepo-t, and Gizmo) for the AGORA High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near z = 2 for each code (hereafter “z ∼ 2”) at which the eight simulations are in the same stage in the target halo’s merger history. After identifying the matched pairs of halos between the CosmoRun simulations and the DMO simulations, we discover that each CosmoRun halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at z ∼ 2, the number of satellite galaxies is significantly fewer than that of dark matter halos in all participating AGORA simulations and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called “missing satellite problem” is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in AGORA and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution (≲100 proper pc at z = 2). We also compare other properties such as the stellar mass–halo mass relation and the mass–metallicity relation. Our work highlights the value of comparison studies such as AGORA, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms.

Mancha3D Code: Multipurpose Advanced Nonideal MHD Code for High-Resolution Simulations in Astrophysics

The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic (MHD) processes in solar/stellar atmospheres. The code includes nonideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high-quality realistic simulations of magnetoconvection, as well as idealized simulations of particular processes, such as wave propagation, instabilities or energetic events. The paper summarizes the equations and methods used in the Mancha3D (Multifluid (-purpose -physics -dimensional) Advanced Non-ideal MHD Code for High resolution simulations in Astrophysics 3D) code. It also describes its numerical stability and parallel performance and efficiency. The code is based on a finite difference discretization and a memory-saving Runge–Kutta (RK) scheme. It handles nonideal effects through super-time-stepping and Hall diffusion schemes, and takes into account thermal conduction by solving an additional hyperbolic equation for the heat flux. The code is easily configurable to perform different kinds of simulations. Several examples of the code usage are given. It is demonstrated that splitting variables into equilibrium and perturbation parts is essential for simulations of wave propagation in a static background. A perfectly matched layer (PML) boundary condition built into the code greatly facilitates a nonreflective open boundary implementation. Spatial filtering is an important numerical remedy to eliminate grid-size perturbations enhancing the code stability. Parallel performance analysis reveals that the code is strongly memory bound, which is a natural consequence of the numerical techniques used, such as split variables and PML boundary conditions. Both strong and weak scalings show adequate performance up to several thousands of processors (CPUs).

TREVR2: illuminating fast N log2N radiative transfer

ABSTRACT We present trevr2 (Tree-based REVerse Ray Tracing 2), a fast, general algorithm for computing the radiation field, suitable for both particle and mesh codes. It is designed to self-consistently evolve chemistry for zoomed-in astrophysical simulations, such as cosmological galaxies with both internal sources and prescribed background radiation, rather than large periodic volumes. Light is propagated until absorbed, with no imposed speed limit other than those due to opacity changes (e.g. ionization fronts). trevr2 searches outward from receiving gas in discrete directions set by the healpix algorithm (unlike its slower predecessor trevr), accumulating optical depth and adding the flux due to sources combined into progressively larger tree cells with distance. We demonstrate Nactivelog2N execution time with absorption and many sources. This allows multiband radiative transfer costs comparable to tree-based gravity and hydrodynamics, and the usual speed-up when active particles evolve on individual time-steps. Sources embedded in non-homogeneous absorbing material introduce systematic errors. We introduce transmission averaging instead of absorption averaging which dramatically reduces these systematic effects. We outline other ways to address systematics including an explicit complex source model. We demonstrate the overall performance of the method via a set of astrophysical test problems.