Search NASA⌕ Search

SEARCH · Search NASA

Results for “astrophysical simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Results from Binary Black Hole Simulations in Astrophysics Applications

Present and planned gravitational wave observatories are opening a new astronomical window to the sky. A key source of gravitational waves is the merger of two black holes. The Laser Interferometer Space Antenna (LISA), in particular, is expected to observe these events with signal-to-noise ratio's in the thousands. To fully reap the scientific benefits of these observations requires a detailed understanding, based on numerical simulations, of the predictions of General Relativity for the waveform signals. New techniques for simulating binary black hole mergers, introduced two years ago, have led to dramatic advances in applied numerical simulation work. Over the last two years, numerical relativity researchers have made tremendous strides in understanding the late stages of binary black hole mergers. Simulations have been applied to test much of the basic physics of binary black hole interactions, showing robust results for merger waveform predictions, and illuminating such phenomena as spin-precession. Calculations have shown that merging systems can be kicked at up to 2500 km/s by the thrust from asymmetric emission. Recently, long lasting simulations of ten or more orbits allow tests of post-Newtonian (PN) approximation results for radiation from the last orbits of the binary's inspiral. Already, analytic waveform models based PN techniques with incorporated information from numerical simulations may be adequate for observations with current ground based observatories. As new advances in simulations continue to rapidly improve our theoretical understanding of the systems, it seems certain that high-precision predictions will be available in time for LISA and other advanced ground-based instruments. Future gravitational wave observatories are expected to make precision.

Baker, John G.↗

Exascale models of stellar explosions: Quintessential multi-physics simulation

The ExaStar project aims to deliver an efficient, versatile, and portable software ecosystem for multi-physics astrophysics simulations run on exascale machines. The code suite is a component-based multi-physics toolkit, built on the capabilities of current simulation codes (in particular Flash-X and Castro), and based on the massively parallel adaptive mesh refinement framework AMReX. It includes modules for hydrodynamics, advanced radiation transport, thermonuclear kinetics, and nuclear microphysics. The code will reach exascale efficiency by building upon current multi- and many-core packages integrated into an orchestration system that uses a combination of configuration tools, code translators, and a domain-specific asynchronous runtime to manage performance across a range of platform architectures. The target science includes multi-physics simulations of astrophysical explosions (such as supernovae and neutron star mergers) to understand the cosmic origin of the elements and the fundamental physics of matter and neutrinos under extreme conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

ZERNIPAX: A fast and accurate Zernike polynomial calculator in Python

Zernike polynomials serve as an orthogonal basis on the unit disc, and have proven to be effective in optics simulations, astrophysics, and more recently in plasma simulations. Unlike Bessel functions, Zernike polynomials are inherently finite and smooth at the disc center (r=0), ensuring continuous differentiability along the axis. This property makes them particularly suitable for simulations, requiring no additional handling at the origin. We developed ZERNIPAX, an open-source Python package capable of utilizing CPU/GPUs, leveraging Google's JAX package and available on GitHub as well as the Python software repository PyPI. Furthermore, our implementation of the recursion relation between Jacobi polynomials significantly improves computation time compared to alternative methods by use of parallel computing while still performing more accurately for high-mode numbers.

Astrophysics↗

Laboratory Studies of the Optical Properties and Condensation Processes of Cosmic Dust Grains

A laboratory facility for conducting a variety of experiments on single isolated dust particles of astrophysical interest levitated in an electrodynamics balance has been developed at NASA/Marshall Space Flight Center. The objective of the research is to employ this experimental technique for studies of the physical and optical properties of individual cosmic dust grains of 0.1-100 micron size in controlled pressure/temperatures environments simulating astrophysical conditions. The physical and optical properties of the analogs of interstellar and interplanetary dust grains of known composition and size distribution will be investigated by this facility. In particular, we will carry out three classes of experiments to study the micro-physics of cosmic dust grains. (1) Charge characteristics of micron size single dust grains to determine the photoelectric efficiencies, yields, and equilibrium potentials when exposed to UV radiation. (2) Infrared optical properties of dust particles (extinction coefficients and scattering phase functions) in the 1-30 micron region using infrared diode lasers and measuring the scattered radiation. (3) Condensation experiments to investigate the condensation of volatile gases on colder nucleated particles in dense interstellar clouds and lower planetary atmospheres. The condensation experiments will involve levitated nucleus dust grains of known composition and initial mass (or m/q ratio), cooled to a temperature and pressure (or scaled pressure) simulating the astrophysical conditions, and injection of a volatile gas at a higher temperature from a controlled port. The increase in the mass due to condensation on the particle will be monitored as a function of the dust particle temperature and the partial pressure of the injected volatile gas. The measured data will permit determination of the sticking coefficients of volatile gases and growth rates of dust particles of astrophysical interest. Some preliminary results based on measurements of photoelectric emission and radiation pressure on single isolated 0.2 to 6.6 micron size silica particles exposed to UV radiation at 120-200 nm and green laser light at 532 nm are presented.

Abbas, M. M.↗

The AGORA High-resolution Galaxy Simulations Comparison Project. III. Cosmological Zoom-in Simulation of a Milky Way–mass Halo

In this work, we present a suite of high-resolution cosmological zoom-in simulations to z = 4 of a 10 12 M⊙ halo at z = 0, obtained using seven contemporary astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, and Gizmo) widely used in the numerical galaxy formation community. The physics prescriptions for gas cooling and heating and star formation are the same as the ones used in our previous Assembling Galaxies of Resolved Anatomy (AGORA) disk comparison but now account for the effects of cosmological processes such as the expansion of the universe, intergalactic gas inflow, and the cosmic ultraviolet background radiation emitted by massive stars and quasars. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. In the first two steps, we methodically calibrate the gas physics, such as cooling and heating, in simulations without star formation. In the third step, we seek agreement on the total stellar mass produced with the common star formation prescription used in the AGORA disk comparison, in stellar-feedback-free simulations. In the last calibration step, we activate stellar feedback, where each code group is asked to set the feedback prescription to as close to the most widely used one in its code community as possible, while aiming for convergence in the stellar mass at z = 4 to the values predicted by semiempirical models. After all the participating code groups successfully complete the calibration steps, we achieve a suite of cosmological simulations with similar mass assembly histories down to z = 4. With numerical accuracy that resolves the internal structure of a target halo (≲100 physical pc at z = 4), we find that the codes overall agree well with one another, e.g., in gas and stellar properties, but also show differences, e.g., in circumgalactic medium (CGM) properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, high-resolution cosmological zoom-in simulations can have robust and reproducible results. New code groups are invited to join and enrich this comparison by generating equivalent models or to test the code's compatibility on their own, by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the zoom-in simulations presented here will be presented in forthcoming reports from the AGORA Collaboration, including studies of the CGM, simulations by additional codes, and results at lower redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing with 3D Adaptive Mesh Refinement and Local Adaptive Time Stepping

General relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multifaceted optimizations that, collectively, accelerate computation by 2–5 orders of magnitude for a wide range of applications. First, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the three dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Second, we demonstrate that local adaptive time stepping on a logarithmic spherical-polar grid accelerates computation by a factor of ≲10 compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of ~10 9 zone cycles per second per node on 5400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate H-AMR's computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of 13,440 × 4608 × 8092 cells and a total of ≲22 billion cells and ~0.65 × 10 8 time steps, it is among the largest astrophysical simulations ever performed. We find that frame dragging by the black hole tears up the disk into two independently precessing subdisks. The innermost subdisk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗

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 10 12 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.

79 ASTRONOMY AND ASTROPHYSICS↗

Heasim and Skyback Simulation Tools and Their Application to the Hitomi Mission

We present an introduction to the heasim multimission observation and skyback background, high-energy pseudo Monte Carlo astrophysical simulation tools. Heasim may be used to accurately and efficiently construct flexible image transport system (FITS) event files for simple or composite sources with a wide range of standard and user-defined spatial, spectral, and temporal characteristics. Skyback is designed to enable users to assess the impact of background discrete and diffuse emission on prospective observations, and skyback output may be directly input into heasim. We present a brief overview of heasim and skyback input, algorithms, usage, and output. We also introduce the sxsbranch tool that computes Hitomi soft X-ray spectrometer resolution grade branching ratios, emphasizing its application to simulations. We include several examples of particular relevance to the Hitomi mission.

Loewenstein, Michael↗

A Search for PAHs in Astrophysical Environments

We present the results of a dedicated search for the spectral signatures in the visible range of neutral polycyclic aromatic hydrocarbons (PAHs) in astronomical observations representing various astrophysical environments, probing a total column of line of sight material corresponding to Av=50. Laboratory measurements of PAHs in simulated astrophysical conditions are now available (see contribution of Salama et al.) which provide for the first time the exact wavelengths for the spectral features of these molecules, as well as detailed information on the intrinsic line profiles and oscillator strengths. These measurements therefore allow a direct comparison to astronomical observations and an estimate of, or upper limit to, the abundance of individual PAHs in space. As the column densities for individual PAHs in interstellar or circumstellar lines of sight are expected to be very low, such a comparison and analysis requires astronomical observations at very high signal to noise. We present such a data set here for lines of sight representing diffuse clouds and circumstellar environments of carbon stars, and their comparison with gas phase spectra of a representative set of free, cold PAHs.

Salama, F.↗

Parallel Multiscale Algorithms for Astrophysical Fluid Dynamics Simulations

Our goal is to develop software libraries and applications for astrophysical fluid dynamics simulations in multidimensions that will enable us to resolve the large spatial and temporal variations that inevitably arise due to gravity, fronts and microphysical phenomena. The software must run efficiently on parallel computers and be general enough to allow the incorporation of a wide variety of physics. Cosmological structure formation with realistic gas physics is the primary application driver in this work. Accurate simulations of e.g. galaxy formation require a spatial dynamic range (i.e., ratio of system scale to smallest resolved feature) of 104 or more in three dimensions in arbitrary topologies. We take this as our technical requirement. We have achieved, and in fact, surpassed these goals.

Norman, Michael L.↗

Sugar and Sugar Derivatives in Residues Produced from the UV Irradiation of Astrophysical Ice Analogs

A large variety and number of organic compounds of prebiotic interest are known to be present in carbonaceous chondrites. Among them, one sugar (dihydroxyacetone) as well as several sugar acids, sugar alcohols, and other sugar derivatives have been reported in the Murchison and Murray meteorites. Their presence, along with amino acids, amphiphiles, and nucleobases strongly suggests that molecules essential to life can form abiotically under astrophysical conditions. This hypothesis is supported by laboratory studies on the formation of complex organic molecules from the ultraviolet (UV) irradiation of simulated astrophysical ice mixtures consisting of H2O, CO, CO2, CH3OH, CH4, NH3, etc., at low temperature. In the past 15 years, these studies have shown that the organic residues recovered at room temperature contain amino acids, amphiphiles, nucleobases, as well as other complex organics. However, no systematic search for the presence of sugars and sugar derivatives in laboratory residues have been reported to date, despite the fact that those compounds are of primary prebiotic significance. Indeed, only small (up to 3 carbon atoms) sugar derivatives including glycerol and glyceric acid have been detected in residues so far.

Ice Irradiation↗

Measurement of Characteristics of Micron Size Individual Dust Particles of Astrophysical Interest

A laboratory facility for levitating single isolated dust particles in an electrodynamic balance has been developed at NASA's Marshall Space Flight Center for conducting studies of the physical and optical properties of the analogs of interstellar and interplanetary dust grains of 0.2-20 micron size under controlled pressures/temperatures simulating astrophysical environments. We plan three classes of experiments using this facility: (1) Charge characteristics of micron size single dust grains: The photoelectric efficiencies, yields, and equilibrium potentials when exposed to UV radiation found from these measurements will provide much-needed photoelectric emission data for individual dust particles; (2) Infrared optical properties of dust particles: Specifically, we will determines the complex refractive indices, the extinction coefficients, the scattering phase functions, and the polarization properties of single dust grains of interest in interstellar environments, in the 1- 25 micron spectral region; (3) Condensation experiments to investigate the deposition of volatile gases on colder nucleated particles in dense interstellar clouds and lower planetary atmospheres: The measured data will permit determination of the sticking efficiencies of volatile gases of astrophysical interest. Brief descriptions of the experimental setup for the last two classes of experiments will be given. We will present results of measurements of photoelectric emission using 0.2-6.6 micron size silica particles exposed to UV radiation at 120-200 nm and also results of radiation pressure measurements using the same size silica particles and laser light at 5320 Angstrom.

Craven, P. D.↗

New Covariant Density Functionals of Nuclear Matter for Compact Star Simulations

We generate three families of extended covariant density functionals of nuclear matter that have varying slope of symmetry energy and skewness at nuclear saturation density, but otherwise share the same basic parameters (symmetry energy, compressibility, saturation parameters, etc.) with the standard DDME2, DD2, and MPE functionals. Tables of the parameters of these new density functionals are given, which can be straightforwardly used in DDME2, DD2, and MPE parameterization-based codes. Furthermore, we provide tables of a large number of equations of state (81 for each family) that can be used in astrophysical simulations to assess the impact of variations of the not-well-known slope of symmetry energy and skewness of nuclear systems on the astrophysics of compact objects. We also provide tables of computed integral parameters (mass, radius, and tidal deformability) that can be used, e.g., for modeling gravitational waveforms. Finally, for the extended DDME2-based parameterization, we implement a first-order phase transition to quark matter to obtain a family of equations of state that accommodates a phase transition to quark matter. Analogous tables of the equations of state and integral parameters are provided for this case as well.

79 ASTRONOMY AND ASTROPHYSICS↗

Laboratory Studies of the Optical Properties and Condensation Processes of Cosmic Dust Particles

A laboratory facility for levitating single isolated dust particles in an electrodynamics balance has been developing at NASA/Marshall Space Flight Center for conducting a variety of experimental, of astrophysical interest. The objective of this research is to employ this innovative experimental technique for studies of the physical and optical properties of the analogs of cosmic grains of 0.2-10 micron size in a chamber with controlled pressure/temperatures simulating astrophysical environments. In particular, we will carry out three classes of experiments to investigate the microphysics of the analogs of interstellar and interplanetary dust grains. (1) Charge characteristics of micron size single dust grains to determine the photoelectric efficiencies, yields, and equilibrium potentials when exposed to UV radiation. These measurements will provide the much-needed photoelectric emission data relating to individual particles as opposed to that for the bulk materials available so far. (2) Infrared optical properties of dust particles obtained by irradiating the particles with radiation from tunable infrared diode lasers and measuring the scattered radiation. Specifically, the complex refractive indices, the extinction coefficients, the scattering phase functions, and the polarization properties of single dust grains of interest in interstellar environments, in the 1-25 micron spectral region will be determined. (3) Condensation experiments to investigate the deposition of volatile gases on colder nucleated particles in dense interstellar clouds and lower planetary atmospheres. The increase in the mass or m/q ratio due to condensation on the particle will be monitored as a function of the dust particle temperature and the partial pressure of the injected volatile gas. The measured data wild permit determination of the sticking efficiencies of volatile gases of astrophysical interest. Preliminary results based on photoelectric emission experiments on 0.2-6.6 micron size silica particles exposed to UV radiation in the 120-200 nm spectral region will be presented.

Abbas, Mian M.↗

Mechanisms for the Formations of the Thymine Under Astrophysical Conditions and Implications for the Origin of Life

Nucleobases are the carriers of the genetic information in ribonucleic acid and deoxyribonucleic acid (DNA) for all life on Earth. Their presence in meteorites clearly indicates that compounds of biological importance can form via non-biological processes in extraterrestrial environments. Recent experimental studies have shown that the pyrimidine-based nucleobases uracil and cytosine can be easily formed from the ultraviolet irradiation of pyrimidine in H2O-rich ice mixtures that simulate astrophysical processes. In contrast, thymine, which is found only in DNA, is more difficult to form under the same experimental conditions, as its formation usually requires a higher photon dose. Earlier quantum chemical studies confirmed that the reaction pathways were favorable provided that several H2O molecules surrounded the reactants. However, the present quantum chemical study shows that the formation of thymine is limited because of the inefficiency of the methylation of pyrimidine and its oxidized derivatives in an H2O ice, as supported by the laboratory studies. Our results constrain the formation of thymine in astrophysical environments and thus the inventory of organic molecules delivered to the early Earth and have implications for the role of thymine and DNA in the origin of life.

Origin↗