Radiation Hydrodynamic Simulations Research Articles

Multiphase Gas Nature in the Sub-parsec Region of the Active Galactic Nuclei. III. Eddington Ratio Dependence on the Structures of Dusty and Dust-free Outflows

We investigated the influence of the Eddington ratio on sub-parsec-scale outflows in active galactic nuclei (AGNs) with supermassive black holes (SMBHs) masses of 107 M ⊙ using two-dimensional radiation hydrodynamics simulations. When the range of the Eddington ratio is γ Edd > 10−3, the radiation force exceeds the gas pressure, leading to stronger outflows and larger dust sublimation radius. Although the sub-parsec-scale outflows are a time-dependent phenomenon, our simulations demonstrated that the radial distributions can be well explained by the steady solutions of the spherically symmetric stellar winds. The dynamic structure of sub-parsec-scale outflows is influenced by the dust sublimation radius and the critical radii determined by the dynamical equilibrium condition. Although significantly affecting the outflow velocity, the Eddington ratio exerts minimal effects on temperature and number density distribution. Furthermore, our analytical solutions highlight the importance of the dust sublimation scale as a crucial determinant of terminal velocity and column density in dusty outflows. Through comparisons of our numerical model with the obscuring fraction observed in nearby AGNs, we reveal insights into the Eddington ratio dependence and the tendency toward the large obscuring fraction of the dusty and dust-free gases. The analytical solutions are expected to facilitate an understanding of the dynamical structure and radiation structures along the line of sight and their viewing angles from observations of ionized outflows.

Dimming events of evolved stars due to clouds of molecular gas. Scenarios based on 3D radiation-hydrodynamics simulations with CO5BOLD

The dramatic dimming episode of the red supergiant Betelgeuse in 2019 and 2020, caused by a partial darkening of the stellar disk, has highlighted gaps in the understanding of the evolution of massive stars. We analyzed numerical models to investigate the processes behind the formation of dark surface patches and the associated reduction in the disk-integrated stellar light. With the CO5BOLD code, we performed global 3D radiation-hydrodynamical simulations of evolved stars, including convection in the stellar interior, self-excited pulsations, and the resulting atmospheric dynamics with strong radiative shocks. We attribute dimming phenomena to obscuring clouds of cool gas in the lower atmosphere, forming according to three different scenarios. One process transports material outward in a strong shock, similar to what occurs in 1D simulations of radially pulsating asymptotic giant branch (AGB) stars. However, in 3D models, deviations from spherical symmetry of the shock front can lead to further local density enhancements. Another mechanism is triggered by a large convective upflow structure, in combination with exceptionally strong radial pulsations. This induces Rayleigh-Taylor instabilities, causing plumes of material to be sent outward into the atmosphere. The third and rarest scenario involves large-amplitude convective fluctuations, leading to enhanced flows in deep downdrafts, which rebound and send material outward. In all cases, the dense gas above the stellar surface cools and darkens rapidly in visible light. AGB stars show localized dark patches regularly during intermediate phases of their large-amplitude pulsations, while more massive stars will only intermittently form such patches during luminosity minima. The episodic levitation of dense gas clumps above the stellar surface, followed by the formation of complex molecules in the cooling gas and possibly dust grains at a later stage, can account for the dark patches and strong dimming events of supergiant stars such as Betelgeuse.

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

X-ray phase-contrast imaging of strong shocks on OMEGA EP.

The ongoing improvement in laser technology and target fabrication is opening new possibilities for diagnostic development. An example is x-ray phase-contrast imaging (XPCI), which serves as an advanced x-ray imaging diagnostic in laser-driven experiments. In this work, we present the results of the XPCI platform that was developed at the OMEGA EP Laser-Facility to study multi-Mbar single and double shocks produced using a kilojoule laser driver. Two-dimensional radiation-hydrodynamic simulations agree well with the shock progression and the spherical curvature of the shock fronts. It is demonstrated that XPCI is an excellent method to determine with high accuracy the front position of a trailing shock wave propagating through an expanding CH plasma that was heated by a precursor Mbar shock wave. The interaction between the rarefaction wave and the shock wave results in a clear signature in the radiograph that is well reproduced by radiation-hydrodynamic simulations.

Spectral reconstruction for radiation hydrodynamic simulations of galaxy evolution

Radiation from stars and active galactic nuclei (AGN) plays an important role in galaxy formation and evolution, and profoundly transforms the intergalactic, circumgalactic, and interstellar medium (IGM, CGM, and ISM). On-the-fly radiative transfer (RT) has started being incorporated in cosmological simulations, but the complex evolving radiation spectra are often crudely approximated with a small number of broad bands with piece-wise constant intensity and a fixed photo-ionisation cross-section. Such a treatment is unable to capture the changes to the spectrum as light is absorbed while it propagates through a medium with non-zero opacity. This can lead to large errors in photo-ionisation and heating rates. In this work we present a novel approach of discretising the radiation field at discrete photon energies, at the edges of the typically used photo-ionising bands, in order to capture the power-law slope of the radiation field. In combination with power-law approximations for the photo-ionisation cross-sections, this model allows us to self-consistently combine radiation from sources with different spectra and accurately follow the ionisation states of primordial and metal species through time. The method is implemented in GASOLINE2 in connection with TREVR2, a fast reverse ray tracing algorithm with 𝒪(Nactive log2 N) scaling. We compare our new piece-wise power-law reconstruction to the piece-wise constant method in calculating the primordial chemistry photo-ionisation and heating rates under an evolving UV background (UVB) and stellar spectrum, and find that our method reduces errors significantly, by up to two orders of magnitude in the case of HeII ionisation. We apply our new spectral reconstruction method in RT post-processing of a cosmological zoom-in simulation, MUGS2 g1536, including radiation from stars and a live UVB, and find a significant increase in total neutral hydrogen (HI) mass in the ISM and the CGM due to shielding of the UVB and a low escape fraction of the stellar radiation. This demonstrates the importance of RT and an accurate spectral approximation in simulating the CGM-galaxy ecosystem.

Multidimensional Radiation Hydrodynamics Simulations of SN 1987A Shock Breakout

Shock breakout is the first electromagnetic signal from supernovae (SNe), which contains important information on the explosion energy and the size and chemical composition of the progenitor star. This paper presents the first two-dimensional (2D) multiwavelength radiation hydrodynamics simulations of SN 1987A shock breakout by using the CASTRO code with the OPAL opacity table considering eight photon groups from infrared to X-ray. To investigate the impact of the pre-SN environment of SN 1987A, we consider three possible circumstellar medium environments: a steady wind, an eruptive mass loss, and the existence of a companion star. In sum, the resulting breakout light curve has an hour-long duration and a peak luminosity of ∼4 × 1046 erg s−1, with a decay rate of ∼3.5 mag hr−1 in X-ray. The dominant band transits to UV around 3 hr after the initial breakout, and its luminosity has a decay rate of ∼1.5 mag hr−1 that agrees well with the observed shock breakout tail. The detailed features of breakout emission are sensitive to the pre-explosion environment. Furthermore, our 2D simulations demonstrate the importance of multidimensional mixing and its impacts on shock dynamics and radiation emission. The mixing emerging from the shock breakout may lead to a global asymmetry of SN ejecta and affect its later SN remnant formation.

Beryllium–tungsten graded density inner shells in double shell capsules for improved hydrodynamic stability

The outer surface of the high-Z inner shell in the double shell configuration of inertial confinement fusion experiments experiences Rayleigh–Taylor instability growth during the implosion process due to inverted density and pressure gradients between a highly compressed foam interstitial layer and the accelerating dense inner shell. Graded density layers have long been known to reduce instability growth rates. In this study, we employ high-fidelity radiation hydrodynamic simulations to demonstrate this improved stability when grading beryllium into tungsten. We first characterize the response to L-band preheat of these layers using a newly calibrated radiation drive. While graded layer capsules suffer reduced performance (here, measured as DD neutron yield from a CD foam fuel) in 1D simulations due to reduced kinetic energy coupling and reduced fuel compression, they suffer less of a performance drop when 2D instabilities are accounted for. With the improved stability of graded layers, we explore the performance of capsules with larger fuel radii and thinner shells as a preliminary study to find new designs in which graded layers produce the highest yields.

Sequential formation of supermassive stars and heavy seed BHs through the interplay of cosmological cold accretion and stellar radiative feedback

ABSTRACT Supermassive stars (SMSs) and heavy seed black holes, as their remnants, are promising candidates for supermassive black hole (SMBH) progenitors, especially for ones observed in the early universe $z\simeq 8.5-10$ by recent JWST observations. Expected cradles of SMSs are the atomic cooling haloes ($M_{\rm halo}\simeq 10^7\ \mathrm{M}_{\odot }$), where ‘cold accretion’ emerges and possibly forms SMSs. We perform a suit of cosmological radiation hydrodynamics simulations and investigate star formation after the emergence of cold accretion, solving radiative feedback from stars inside the halo. We follow the mass growth of the protostars for $\sim 3\ \mathrm{Myr}$, resolving the gas inflow down to $\sim 0.1\ \mathrm{pc}$ scales. We discover that, after cold accretion emerges, multiple SMSs of $m_{\star }\gtrsim 10^5\ \mathrm{M}_{\odot }$ form at the halo centre with the accretion rates maintained at $\dot{m}_{\star }\simeq 0.04\ \mathrm{M}_\odot \mathrm{yr}^{-1}$ for $\lesssim 3\ \mathrm{Myr}$. Cold accretion supplies gas at a rate of $\dot{M}_{\rm gas}\gtrsim 0.01-0.1\ \mathrm{M}_\odot \mathrm{yr}^{-1}$ from outside the halo virial radius to the central gas disc. Gravitational torques from spiral arms transport gas further inwards, which feeds the SMSs. Radiative feedback from stars suppresses $\mathrm{H}_2$ cooling and disc fragmentation, while photoevaporation is prevented by a dense envelope, which attenuates ionizing radiation. Our results suggest that cold accretion can bring efficient BH mass growth after seed formation in the later universe. Moreover, cold accretion and gas migration inside the central disc increase the mass concentration and provide a promising formation site for the extremely compact stellar clusters observed by JWST.

Role of nonlocal heat transport on the laser ablative Rayleigh-Taylor instability

Abstract Ablative Rayleigh–Taylor instability (ARTI) and nonlocal heat transport are the critical problems in laser-driven inertial confinement fusion, while their coupling with each other is not completely understood yet. Here the ARTI in the presence of nonlocal heat transport is studied self-consistently for the first time theoretically and by using radiation hydrodynamic simulations. It is found that the nonlocal heat flux generated by the hot electron transport tends to attenuate the growth of instability, especially for short wavelength perturbations. A linear theory of the ARTI coupled with the nonlocal heat flux is developed, and a prominent stabilization of the ablation front via the nonlocal heat flux is found, in good agreement with numerical simulations. This effect becomes more significant as the laser intensity increases. Our results should have important references for the target designing for inertial confinement fusion.

Spectral behavior and expansion dynamics of Gd plasma generated by dual-pulse laser irradiation.

Laser-produced gadolinium plasma (Gd-LPP) emerges as a promising candidate for next-generation nanolithography light sources. In this study, a dual laser pulse scheme was implemented to achieve a narrow spectral peak. By varying the pre-main pulse delay and pre-pulse laser energy, optimal conditions of 40 ns delay and 50 mJ energy were identified to improve spectral purity. Radiation hydrodynamics simulations revealed that the improved spectral purity stems from a flatter density gradient at the ablation front and a lower average electron density in the EUV emission region. Additionally, reheating the pre-formed plasma with a short main pulse mitigated plasma squeezing, resulting in an even lower electron density and thus improved spectral purity. Our findings suggest that spectral narrowing in the dual-pulse scheme, essential for better matching with multilayer reflection bandwidths, can be optimized through precise control of pre-pulse energy, pre-main delay, and main-pulse duration.

Modeling the hundreds-of-nanoseconds-long irradiation of tin droplets with a 2 µm-wavelength laser for future EUV lithography

Abstract The properties of an extreme ultraviolet (EUV) source driven by hundreds-of-nanoseconds-long laser pulses of λ laser = 2 µm wavelength are investigated through radiation-hydrodynamic simulations. We show that single-pulse irradiation of 30 µm-diameter tin droplets can generate in-band energies ∼ 30 mJ directed in the 2π sr solid angle subtended by the collector mirror, yielding energies ∼100 mJ produced from 45 µm-diameter droplets. Time-integrated in-band EUV emission profiles, generated by taking Abel transforms of the local net in-band emissivity, reveal EUV source sizes in the axial (laser) direction × radial direction of ∼600 (1000) × 100 (150) µm2 for 30 (45) µm-diameter droplets. We find that approximately 74% of the total in-band emission is produced during the first half of the total propulsion distance irrespective of droplet size or laser intensity. Furthermore, we propose a one-dimensional analytical propulsion model to qualitatively explain the simulated droplet trajectory and to predict the time taken for the droplet to vaporize. These findings offer motivation for the development of high-power EUV sources based on single-pulse, hundreds-of-nanoseconds-long irradiation of tin droplets.

Radiation Hydrodynamic Simulations of Massive Stars in Gas-rich Environments: Accretion of AGN Stars Suppressed by Thermal Feedback

Massive stars may form in or be captured into active galactic nuclei (AGN) disks. Recent 1D studies employing stellar-evolution codes have demonstrated the potential for rapid growth of such stars through accretion up to a few hundred solar masses. We perform 3D radiation hydrodynamic simulations of moderately massive stars’ envelopes in order to determine the rate and critical radius R crit of their accretion process in an isotropic gas-rich environment in the absence of luminosity-driven mass loss. We find that in the “fast-diffusion” regime where characteristic radiative diffusion speed c/τ is faster than the gas sound speed c s , the accretion rate is suppressed by feedback from gravitational and radiative advection energy flux, in addition to the stellar luminosity. Alternatively, in the “slow-diffusion” regime where c/τ < c s , due to adiabatic accretion, the stellar envelope expands quickly to become hydrostatic and further net accretion occurs on thermal timescales in the absence of self-gravity. When the radiation entropy of the medium is less than that of the star, however, this hydrostatic envelope can become more massive than the star itself. Within this subregime, the self-gravity of the envelope excites runaway growth. Applying our results to realistic environments, moderately massive stars (≲100M ⊙) embedded in AGN disks typically accrete in the fast-diffusion regime, leading to a reduction of steady-state accretion rate 1–2 orders of magnitudes lower than expected by previous 1D calculations and R crit smaller than the disk scale height, except in the opacity window at temperature T ∼ 2000 K. Accretion in slow diffusion regime occurs in regions with very high density ρ ≳ 10−9 g cm−3, and needs to be treated with caution in 1D long-term calculations.

Pre-peak Emission in Tidal Disruption Events

The rising part of a tidal disruption event light curve provides unique insight into early emission and the onset of accretion. Various mechanisms are proposed to explain the pre-peak emission, including shocks from debris interaction and reprocessing of disk emission. We study the pre-peak emission and its influence on the gas circularization by a series of gray radiation hydrodynamic simulations with varying black hole mass. We find that, given a super-Eddington fallback rate of 10ṀEdd , the stream–stream collision can occur multiple times and drive strong outflows of up to 9ṀEdd . By dispersing gas to ≳100r s , the outflow can delay gas circularization and leads to sub-Eddington accretion rates during the first few stream–stream collisions. The stream–stream collision shock and circularization shock can sustain a luminosity of ∼1044 erg s−1 for days. The luminosity is generally sub-Eddington and shows a weak correlation with accretion rate at early times. The outflow is optically thick, yielding a reprocessing layer with a size of ∼1014 cm and photospheric temperature of ∼4 × 104 K.

Three-dimensional reconstruction of laser-direct-drive inertial confinement fusion hot-spot plasma from x-ray diagnostics on the OMEGA laser facility (invited).

A deep-learning convolutional neural network (CNN) is used to infer, from x-ray images along multiple lines of sight, the low-mode shape of the hot-spot emission of deuterium-tritium (DT) laser-direct-drive cryogenic implosions on OMEGA. The motivation of this approach is to develop a physics-informed 3-D reconstruction technique that can be performed within minutes to facilitate the use of the results to inform changes to the initial target and laser conditions for the subsequent implosion. The CNN is trained on a 3D radiation-hydrodynamic simulation database to relate 2D x-ray images to 3D emissivity at stagnation. The CNN accounts for the lack of an absolute spatial reference and the different bands of photon energies in the x-ray images. While previous work [O. M. Mannion etal., Phys. Plasmas 28, 042701 (2021) and A. Lees etal., Phys. Rev. Lett. 127, 105001 (2021)] studied the effect of mode-1 asymmetries on implosion performance using nuclear diagnostics, this work focuses on the effect of mode 2 inferred from x-ray diagnostics on implosion performance. A current analysis of 19 DT cryogenic implosions indicates there is an upper limit of ∼20% reduction in the neutron yield caused by an ℓ = 2 amplitude for ℓ2/ℓ0 ≤ 0.32. These conclusions are supported by 2D simulations.

The Primary Flare Following a Stellar Collision in a Galactic Nucleus

High-velocity stellar collisions near supermassive black holes may result in a complete disruption of the stars. The initial disruption can have energies on par with supernovae and power a very fast transient. In this work, we examine the primary flare that follows the initial transient, which arises when streams of gas from the disrupted stars travel around the central black hole and collide with each other on the antipodal side with respect to the original collision. We present a simple analytic estimate for the properties of the flare, which depends on the distance of the collision from the central black hole and on the center of mass velocity of the colliding stars. We also present the first-of-their-kind radiation-hydrodynamics simulations of a few examples of stellar collisions and postcollision flow of the ejected gas and calculate the expected bolometric light curves. We find that such postcollision flares are expected to be similar to flares that arise in tidal disruptions events of single stars.

3D Radiation-hydrodynamical Simulations of Shadows on Transition Disks

Shadows are often observed in transition disks, which can result from obscuring by materials closer to the star, such as a misaligned inner disk. While shadows leave apparent darkened emission as observational signatures, they have significant dynamical impact on the disk. We carry out 3D radiation-hydrodynamical simulations to study shadows in transition disks and find that the temperature drop due to the shadow acts as an asymmetric driving force, leading to spirals in the cavity. These spirals have zero pattern speed following the fixed shadow. The pitch angle is given by tan−1(c s /v ϕ ) (6° if h/r = 0.1). These spirals transport mass through the cavity efficiently, with α ∼ 10−2 in our simulation. Besides spirals, the cavity edge can also form vortices and flocculent streamers. When present, these features could disturb the shadow-induced spirals. By carrying out Monte Carlo radiative transfer simulations, we show that these features resemble those observed in near-infrared scattered light images. In the vertical direction, the vertical gravity is no longer balanced by the pressure gradient alone. Instead, an azimuthal convective acceleration term balances the gravity–pressure difference, leading to azimuthally periodic upward and downward gas motion reaching 10% of the sound speed, which can be probed by Atacama Large Millimeter/submillimeter Array line observations.

Element Formation in Radiation-hydrodynamics Simulations of Kilonovae

Understanding the details of r-process nucleosynthesis in binary neutron star merger (BNSM) ejecta is key to interpreting kilonova observations and identifying the role of BNSMs in the origin of heavy elements. We present a self-consistent, two-dimensional, ray-by-ray radiation-hydrodynamic evolution of BNSM ejecta with an online nuclear network (NN) up to a timescale of days. For the first time, an initial numerical relativity ejecta profile composed of the dynamical component and spiral-wave and disk winds is evolved including detailed r-process reactions and nuclear heating effects. A simple model for the jet energy deposition is also included. Our simulation highlights that the common approach of relating in postprocessing the final nucleosynthesis yields to the initial thermodynamic profile of the ejecta can lead to inaccurate predictions. Moreover, we find that neglecting the details of the radiation-hydrodynamic evolution of the ejecta in nuclear calculations can introduce deviations of up to 1 order of magnitude in the final abundances of several elements, including very light and second r-process peak elements. The presence of a jet affects element production only in the innermost part of the polar ejecta, and it does not alter the global nucleosynthesis results. Overall, our analysis shows that employing an online NN improves the reliability of nucleosynthesis and kilonova light-curve predictions.

A GPU based 3D raytracing algorithm for DUED laser fusion code

Abstract These days, graphical processing units (GPUs) deliver performance&amp;#xD;comparable to that of hundreds of CPU cores. This level of performance allows&amp;#xD;certain classes of simulations to be run in-house on a standard consumer workstation,&amp;#xD;eliminating the need for a cluster. In this paper, it is shown that medium-resolution,&amp;#xD;2D radiation hydrodynamics simulations for laser-driven inertial confinement fusion&amp;#xD;with realistic 3D laser raytracing can now be conducted on a single consumer device.&amp;#xD;A novel raytracing module has indeed been developed for the 2D Lagrangian radiation-&amp;#xD;hydro-nuclear code DUED to leverage the computational power of GPUs. By&amp;#xD;employing 3D raytracing, more realistic investigations of laser-driven plasmas become&amp;#xD;feasible, with a particular focus on perturbations resulting from non-uniform laser&amp;#xD;irradiation.

A Thermodynamic Criterion for the Formation of Circumplanetary Disks

The formation of circumplanetary disks is central to our understanding of giant planet formation, influencing their growth rate during the post-runaway phase and observability while embedded in protoplanetary disks. We use three-dimensional global multifluid radiation hydrodynamics simulations with the FARGO3D code to define the thermodynamic conditions that enable circumplanetary disk formation around Jovian planets on wide orbits. Our simulations include stellar irradiation, viscous heating, static mesh refinement, and active calculation of opacity based on multifluid dust dynamics. We find a necessary condition for the formation of circumplanetary disks in terms of a mean cooling time: When the cooling time is at least 1 order of magnitude shorter than the orbital timescale, the specific angular momentum of the gas is nearly Keplerian at scales of one-third of the Hill radius. We show that the inclusion of multifluid dust dynamics favors rotational support because dust settling produces an anisotropic opacity distribution that favors rapid cooling. In all our models with radiation hydrodynamics, specific angular momentum decreases as time evolves, in agreement with the formation of an inner isentropic envelope due to compressional heating. The isentropic envelope can extend up to one-third of the Hill radius and shows negligible rotational support. Thus, our results imply that young gas giant planets may host spherical isentropic envelopes, rather than circumplanetary disks.

Early-Time Harmonic Generation from a Single-Mode Perturbation Driven by X-Ray Ablation.

X-ray ablation dynamics of the planar foil with preimposed sinusoidal ripples is investigated at the SG 100kJ Laser Facility. A significant fraction of the second harmonics is observed and identified at the beginning of the ablative drive when the amplitude of the perturbation is within the linear regime. With radiation-hydrodynamic simulations and a developed simple model, we can reveal that such a novel phenomenon is due to the fact that a sustained deformation of the ablation front is initiated since the ablation pressure is directed to the normal direction of the perturbation surface. We find that this deformation dominates the early-time perturbation evolution and results in a specific stage in addition to the traditional ablative Richtmyer-Meshkov instability phase. Our results can be applicable to various regions such as implosions in inertial confinement fusion and the dynamics of molecular clouds in astrophysics.