Radiation Hydrodynamics Model Research Articles

Global classical solutions to a multidimensional radiation hydrodynamics model with symmetry and large initial data

AbstractAs a first stage to study the global large solutions of the radiation hydrodynamics model with viscosity and thermal conductivity in the high‐dimensional space, we study the problems in high dimensions with some symmetry, such as the spherically or cylindrically symmetric solutions. Specifically, we will study the global classical large solutions to the radiation hydrodynamics model with spherically or cylindrically symmetric initial data. The key point is to obtain the strict positive lower and upper bounds of the density and the lower bound of the temperature . Compared with the Navier–Stokes equations, these estimates in the present paper are more complicated due to the influence of the radiation. To overcome the difficulties caused by the radiation, we construct a pointwise estimate between the radiative heat flux and the temperature by studying the boundary value problem of the corresponding ordinary differential equation. And we consider a general heat conductivity: if ; if . This can be viewed as the first result about the global classical large solutions of the radiation hydrodynamics model with some symmetry in the high‐dimensional space.

A Multiwavelength Survey of Nearby M Dwarfs: Optical and Near-ultraviolet Flares and Activity with Contemporaneous TESS, Kepler/K2, Swift, and HST Observations

We present a comprehensive multiwavelength investigation into flares and activity in nearby M dwarf stars. We leverage the most extensive contemporaneous data set obtained through the Transiting Exoplanet Sky Survey, Kepler/K2, the Neil Gehrels Swift Observatory, and the Hubble Space Telescope, spanning the optical and near-ultraviolet (NUV) regimes. In total, we observed 213 NUV flares on 24 nearby M dwarfs, with ∼27% of them having detected optical counterparts, and found that all optical flares had NUV counterparts. We explore NUV/optical energy fractionation in M dwarf flares. Our findings reveal a slight decrease in the ratio of optical to NUV energies with increasing NUV energies, a trend in agreement with prior investigations on G–K stars’ flares at higher energies. Our analysis yields an average NUV fraction of flaring time for M0–M3 dwarfs of 2.1%, while for M4–M6 dwarfs it is 5%. We present an empirical relationship between NUV and optical flare energies and compare to predictions from radiative hydrodynamic and blackbody models. We conducted a comparison of the flare frequency distribution (FFDs) of NUV and optical flares, revealing that the FFDs of both NUV and optical flares exhibit comparable slopes across all spectral subtypes. NUV flares on stars affect the atmospheric chemistry, the radiation environment, and the overall potential to sustain life on any exoplanets they host. We find that early and mid-M dwarfs (M0–M5) have the potential to generate NUV flares capable of initiating abiogenesis.

Global Classical Large Solutions for the Radiation Hydrodynamics Model in Unbounded Domains

Global Classical Large Solutions for the Radiation Hydrodynamics Model in Unbounded Domains

M3DIS - A grid of 3D radiation-hydrodynamics stellar atmosphere models for stellar surveys. I. Procedure, validation, and the Sun

Large-scale stellar surveys, such as SDSS-V, 4MOST, WEAVE, and PLATO, require accurate atmospheric models and synthetic spectra of stars for accurate analyses of fundamental stellar parameters and chemical abundances. The primary goal of our work is to develop a new approach to solve radiation-hydrodynamics (RHD) and generate model stellar spectra in a self-consistent and highly efficient framework. We build upon the Copenhagen legacy RHD code, the MULTI3D non-local thermodynamic equilibrium (NLTE) code, and the DISPATCH high-performance framework. The new approach allows us to calculate 3D RHD models of stellar atmospheres on timescales of a few thousand CPU hours and to perform subsequent spectrum synthesis in local thermodynamic equilibrium (LTE) or NLTE for the desired physical conditions within the parameter space of FGK-type stars. We compare the 3D RHD solar model with other available models and validate its performance against solar observations, including the centre-to-limb variation of intensities and key solar diagnostic lines of H and Fe. We show that the performance of the new code allows to overcome the main bottleneck in 3D NLTE spectroscopy and enables calculations of multi-dimensional grids of synthetic stellar observables for comparison with modern astronomical observations.

On the role of target mass in extreme ultraviolet light generation from CO2-driven tin plasmas for nanolithography

Target conditioning is a crucial ingredient of high-power extreme ultraviolet (EUV) source operation in state-of-the-art nanolithography. It involves deforming tin microdroplets into tens of nanometer-thin sheets, sheets which are subsequently irradiated by intense CO2 laser radiation to form a hot, EUV-emitting plasma. Recent experiments have found that a substantial fraction of the initial droplet mass is lost in the deformation phase through fragmentation. The goal of the present study is to investigate, using radiation-hydrodynamic modeling, how variations in the sheet mass affect EUV source power and the laser-to-in-band conversion efficiency (CE). It is found that high-mass sheets can “feed” the plasma with sufficient mass to sustain the production of in-band-emitting charge states over the course of laser irradiation. Low-mass sheets, on the contrary, cannot supply enough mass to sustain this production over the pulse, thus leading to a reduction in in-band power and CE. The dependence of CE on laser energy and target thickness is quantified, and a rather weak reduction of CE with increasing laser energy for high-mass sheets is identified.

Circumnuclear Multiphase Gas in the Circinus Galaxy. VI. Detectability of Molecular Inflow and Atomic Outflow

Recent submillimeter observations have revealed signs of parsec-scale molecular inflow and atomic outflow in the nearest Seyfert 2 galaxy, the Circinus galaxy. To verify the gas kinematics suggested by these observations, we performed molecular and atomic line transfer calculations based on a physics-based 3D radiation-hydrodynamic model, which has been compared with multiwavelength observations in this paper series. The major-axis position–velocity diagram (PVD) of CO(3–2) reproduces the observed faint emission at the systemic velocity, and our calculations confirm that this component originates from failed winds falling back to the disk plane. The minor-axis PVD of [C i](3 P 1–3 P 0), when created using only the gas with positive radial velocities, presents a sign of blueshifted and redshifted offset peaks similar to those in the observation, suggesting that the observed peaks indeed originate from the outflow, but that the model may lack outflows as strong as those in the Circinus galaxy. Similar to the observed HCN(3–2), the similar dense-gas tracer HCO+(3–2) can exhibit nuclear spectra with inverse P-Cygni profiles with ∼0.5 pc beams, but the line shape is azimuthally dependent. The corresponding continuum absorbers are inflowing clumps at 5–10 pc from the center. To detect significant absorption with a high probability, the inclination must be fairly edge-on (≳85°), and the beam size must be small (≲1 pc). These results suggest that HCN or HCO+ and [C i] lines are effective for observing parsec-scale inflows and outflows, respectively.

An improved methodology for modeling short pulse buried layer x-ray emission spectra

Radiation-hydrodynamic and spectroscopic modeling are important aspects of high energy density experimental design. In this paper, we improve the performance and capabilities over those obtainable with a previous methodology used for simulating x-ray emission spectra from buried layer targets heated by short pulse lasers. The improvement incorporates post-processing HYDRA radiation-hydrodynamic output with a non-local thermodynamic equilibrium atomic-kinetics radiation transport code, Cretin. Each code uses an independent radiation field which allows decoupling HYDRA's radiation group structure from Cretin's spectral output to improve the speed and flexibility of the design methodology. The execution time decreases from 2–3 days to a few hours while the flexibility of the improved methodology allows for performing sensitivity studies including a comparison of steady-state and time-dependent atomic kinetics and differences in the radiation group structure.

Cosmic variance suppression in radiation-hydrodynamic modelling of the reionization-era 21-cm signal

ABSTRACT The 21-cm line emitted by neutral hydrogen is the most promising probe of the Epoch of Reionization (EoR). Multiple radio interferometric instruments are on the cusp of detecting its power spectrum. It is therefore essential to deliver robust theoretical predictions, enabling sound inference of the coeval Universe properties. The nature of this signal traditionally required the modelling of $\mathcal {O}(10^{7-8} \, {\rm Mpc}^3)$ volumes to suppress the impact of cosmic variance. However, the recently proposed Fixed & Paired (F&P) approach uses carefully crafted simulation pairs to achieve equal results in smaller volumes. In this work, we thoroughly test the applicability of and improvement granted by this technique to different observables of the 21-cm signal from the EoR. We employ radiation-magneto-hydrodynamics simulations to ensure the most realistic physical description of this epoch, greatly improving over previous studies using a seminumerical approach without accurate galaxy formation physics and radiative transfer. We estimate the statistical improvement granted by the F&P technique on predictions of the skewness, power spectrum, bispectrum, and ionized regions size distribution of the 21-cm signal at redshift 7 ≤ z ≤ 10 (corresponding to ${\ge}80{{\ \rm per\ cent}}$ of the gas being neutral). We find that the effective volume of F&P simulations is at least 3.5 times larger than traditional simulations. This directly translates into an equal improvement in the computational cost (in terms of time and memory). Finally, we confirm that a combination of different observables like skewness, power spectrum, and bispectrum across different redshifts can be utilized to maximize the improvement.

3D non-LTE modeling of the stellar center-to-limb variation for transmission spectroscopy studies

Context. Transmission spectroscopy is one of the most powerful techniques used to characterize transiting exoplanets, since it allows for the abundance of the atomic and molecular species in the planetary atmosphere to be measured. However, stellar lines may bias the determination of such abundances if their center-to-limb variations (CLVs) are not properly accounted for. Aims. This paper aims to show that three-dimensional (3D) radiation hydrodynamic models and the assumption of non-local ther-modynamic equilibrium (non-LTE) line formation are required for an accurate modeling of the stellar CLV of the Na I D1 and K I resonance lines on transmission spectra. Methods. We modeled the CLV of the Na I D1 and K I resonance lines in the Sun with 3D non-LTE radiative transfer. The synthetic spectra were compared to solar observations with high spatial and spectral resolution, including new data collected with the CRISP instrument at the Swedish 1-m Solar Telescope between µ = 0.1 and µ = 1.0. Results. Our 3D non-LTE modeling of the Na I D1 resonance line at 5896 Å and the K I 7699 Å resonance line in the Sun is in good agreement with the observed CLV in the solar spectrum. Moreover, the simulated CLV curve for a Jupiter-Sun system inferred with a 3D non-LTE analysis shows significant differences from the one obtained from a 1D atmosphere. The latter does indeed tend to overestimate the amplitude of the transmission curve by a factor that is on the same order of magnitude as a planetary absorption depth (i.e., up to 0.2%). Conclusions. This work highlights the fact that to correctly characterize exoplanetary atmospheres, 3D non-LTE synthetic spectra ought to be used to estimate the stellar CLV effect in transmission spectra of solar-like planet hosts. Moreover, since different spectral lines show different CLV curves for the same geometry of the planet-star system, it is fundamental to model the CLV individually for each line of interest. The work will be extended to other lines and FGK-type stars, allowing for synthetic high-resolution spectra to mitigate the stellar contamination of low-resolution planetary spectra, for example, those drawn from JWST.

Bridging the Gap between Luminous Red Novae and Common Envelope Evolution: The Role of Recombination Energy and Radiation Force

Luminous red novae and their connection to common envelope evolution (CEE) remain elusive in astrophysics. Here, we present a radiation hydrodynamic model capable of simulating the light curves of material ejected during a CEE. For the first time, the radiation hydrodynamic model incorporates complete recombination physics for hydrogen and helium. The radiation hydrodynamic equations are solved with Guangqi. With time-independent ejecta simulations, we show that the peaks in the light curves are attributed to radiation-dominated ejecta, while the extended plateaus are produced by matter-dominated ejecta. To showcase our model’s capability, we fit the light curve of AT 2019zhd. The central mass object of 6 M ⊙ is assumed based on observations and scaling relations. Our model demonstrates that the ejecta mass of AT 2019zhd falls within the range of 0.04–0.1 M ⊙. Additionally, we demonstrate that recombination energy and radiation force acceleration significantly impact the light curves, whereas dust formation has a limited effect during the peak and plateau phases.

3D NLTE modelling of Y and Eu

Context. Abundances of s- and r-process elements in Sun-like stars constrain nucleosynthesis in extreme astrophysical events, such as compact binary mergers and explosions of highly magnetised rapidly rotating massive stars. Aims. We measure solar abundances of yttrium (Y) and europium (Eu) using 3D non-local thermal equilibrium (NLTE) models. We use the model to determine the abundance of Y, and also explore the model’s ability to reproduce the solar centre-to-limb variation of its lines. In addition, we determine the Eu abundance using solar disc-centre and integrated flux spectra. Methods. We developed an NLTE model of Eu and updated our model of Y with collisional data from detailed quantum-mechanical calculations. We used the IAG spatially resolved high-resolution solar spectra to derive the solar abundances of Y across the solar disc and of Eu for integrated flux and at disc centre using a set of carefully selected lines and a 3D radiation-hydrodynamics model of the solar atmosphere. Results. We find 3D NLTE solar abundances of A(Y)3D NLTE = 2.30 ± 0.03stat ± 0.07syst dex based on observations at all angles and A(Eu) = 0.57 ± 0.01stat ± 0.06syst dex based on the integrated flux and disc-centre intensity. 3D NLTE modelling offers the most consistent abundances across the solar disc, and resolves the problem of severe systematic bias in Y and Eu abundances inherent to 1D LTE, 1D NLTE, and 3D LTE modelling.

Reconstruction and analysis of transient evolution images of laser-produced plasma plumes.

We introduce a method for the analysis and simulation of transient images of laser-produced plasma (LPP) plumes. This method comprises three steps: (i) calculating the two-dimensional distribution of plasma parameters using a radiation hydrodynamics model, (ii) constructing radiation paths through ray tracing, and (iii) solving the radiation transport equation along these paths. In our simulations, we have meticulously considered factors that could influence the imaging results, including the quantum efficiency to different radiation wavelengths, the imaging lens' transmittance, the target surface's reflectivity, and the absorption, emission, and scattering quantum effect of the detector processes occurring in the plasma. We applied this method to analyze and simulate the transient images of aluminum plasma plumes in a background air environment at a pressure of 2000 Pa. The results demonstrate that our method not only produces simulated images that align with experimental results but also provides a reliable distribution of plasma state parameters and clearly identifies the ion species radiating in different bands. Given its capability in transient image reconstruction and its adaptability as a tool for spectral simulation and analysis of LPPs, we believe this method holds significant potential for spectral diagnostics in fields such as laser-induced breakdown spectroscopy, extreme ultraviolet lithography sources, and high-energy-density physics, among others.

Can Proton Beam Heating Flare Models Explain Sunquakes?

Solar Dynamics Observatory (SDO)/Helioseismic and Magnetic Imager (HMI) observations reveal a class of solar flares with substantial energy and momentum impacts in the photosphere, concurrent with white-light emission and helioseismic responses, known as sunquakes. Previous radiative hydrodynamic modeling has demonstrated the challenges of explaining sunquakes in the framework of the standard flare model of “electron beam” heating. One of the possibilities to explain the sunquakes and other signatures of the photospheric impact is to consider additional heating mechanisms involved in solar flares, for example via flare-accelerated protons. In this work, we analyze a set of single-loop Fokker–Planck and radiative hydrodynamics RADYN+FP simulations where the atmosphere is heated by nonthermal power-law-distributed proton beams which can penetrate deeper than the electron beams into the low atmospheric layers. Using the output of the RADYN models, we calculate synthetic Fe i 6173 Å line Stokes profiles and from those the line-of-sight observables of the SDO/HMI instrument, as well as the 3D helioseismic response, and compare them with the corresponding observational characteristics. These initial results show that the models with proton beam heating can produce the enhancement of the HMI continuum observable and explain qualitatively the generation of sunquakes. The continuum observable enhancement is evident in all models but is more prominent in ones with E c ≥ 500 keV. In contrast, the models with E c ≤ 100 keV provide a stronger sunquake-like helioseismic impact according to the 3D acoustic modeling, suggesting that low-energy (deka- and hecto-keV) protons have an important role in the generation of sunquakes.

Spectro-polarimetric Properties of Sunquake Sources in X1.5 Flare and Evidence for Electron and Proton Beam Impacts

The first significant sunquake event of Solar Cycle 25 was observed during the X1.5 flare of 2022 May 10, by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory. We perform a detailed spectro-polarimetric analysis of the sunquake photospheric sources, using the Stokes profiles of the Fe i 6173 Å line, reconstructed from the HMI linear and circular polarized filtergrams. The results show fast variations of the continuum emission with rapid growth and slower decay lasting 3–4 minutes, coinciding in time with the hard X-ray impulses observed by the Konus instrument on board the Wind spacecraft. The variations in the line core appeared slightly ahead of the variations in the line wings, showing that the heating started in the higher atmospheric layers and propagated downward. The most significant feature of the line profile variations is the transient emission in the line core in three of the four sources, indicating intense, impulsive heating in the lower chromosphere and photosphere. In addition, the observed variations of the Stokes profiles reflect transient and permanent changes in the magnetic field strength and geometry in the sunquake sources. Comparison with the radiative hydrodynamics models shows that the physical processes in the impulsive flare phase are substantially more complex than those predicted by proton and electron beam flare models currently presented in the literature.

Global Regularity for a Radiation Hydrodynamics Model with Viscosity and Thermal Conductivity

Global Regularity for a Radiation Hydrodynamics Model with Viscosity and Thermal Conductivity

Space-time behavior for radiative hydrodynamics model with or without heat conduction

We consider space-time behaviors of smooth solutions for the radiative hydrodynamics system with or without heat conduction in the whole space \(R^3\) by using Green's function method. This result exhibits the generalized Huygens' principle as the classical compressible Navier-Stokes equations [3,26], which is different from the Hamer model for radiating gases in [36].
 For more information see https://ejde.math.txstate.edu/Volumes/2023/60/abstr.html

Three-component modelling of O-rich AGB star winds

Stellar winds of cool and pulsating asymptotic giant branch (AGB) stars enrich the interstellar medium with large amounts of processed elements and various types of dust. We present the first study on the influence of gas-to-dust drift on ab initio simulations of stellar winds of M-type stars driven by radiation pressure on forsterite particles. Our study is based on our radiation hydrodynamic model code T-800 that includes frequency-dependent radiative transfer, dust extinction based on Mie scattering, grain growth and ablation, gas-to-dust drift using one mean grain size, a piston that simulates stellar pulsations, and an accurate high spatial resolution numerical scheme. To enable this study, we calculated new gas opacities based on the EXOMOL database, and we extended the model code to handle the formation of minerals that may form in M-type stars. We determine the effects of drift by comparing drift models to our new and extant non-drift models. Three out of four new drift models show high drift velocities, 87–310 km s−1. Our new drift model mass-loss rates are 1.7–13 per cent of the corresponding values of our non-drift models, but compared to the results of two extant non-drift models that use the same stellar parameters, these same values are 0.33–1.5 per cent. Meanwhile, a comparison of other properties such as the expansion velocity and grain size show similar values. Our results, which are based on single-component forsterite particles, show that the inclusion of gas-to-drift is of fundamental importance in stellar wind models driven by such transparent grains. Assuming that the drift velocity is insignificant, properties such as the mass-loss rate may be off from more realistic values by a factor of 50 or more.

Asymptotic stability of viscous contact wave to a radiation hydrodynamic limit model

This paper is concerned with the large time behavior of the solutions for 1D radiation hydrodynamic limit model without viscosity and its asymptotic stability of the viscous contact discontinuity wave under the smallness assumption of the strength of the contact wave and initial perturbations. The present pressure includes a fourth-order term about the absolute temperature from radiation effect which brings the main difficulty. Furthermore, the dissipative of the system is weaker for the lack of viscosity. All these make the problem more challenging. The prove is mainly based on the energy method, including normal and radial directions energy estimates.

Turbulence-supported Massive Star Envelopes

The outer envelopes of massive (M ≳ 10 M ⊙) stars exhibit large increases in opacities from forests of lines and ionization transitions (particularly from iron and helium) that trigger near-surface convection zones. One-dimensional (1D) models predict density inversions and supersonic motions that must be resolved with computationally intensive three-dimensional (3D) radiation hydrodynamic (RHD) modeling. Only in the last decade have computational tools advanced to the point where ab initio 3D models of these turbulent envelopes can be calculated, enabling us to present five 3D RHD Athena++ models (four previously published and one new 13 M ⊙ model). When convective motions are subsonic, we find excellent agreement between 3D and 1D velocity magnitudes, stellar structure, and photospheric quantities. However, when convective velocities approach the sound speed, hydrostatic balance fails as the turbulent pressure can account for 80% of the force balance. As predicted by Henyey, we show that this additional pressure support leads to a modified temperature gradient, which reduces the superadiabaticity where convection is occurring. In addition, all five models display significant overshooting from the convection in the Fe convection zone. As a result, the turbulent velocities at the surface are indicative of those in the Fe zone. There are no confined convection zones as seen in 1D models. In particular, helium convection zones seen in 1D models are significantly modified. Stochastic low-frequency brightness variability is also present in the 13 M ⊙ model with comparable amplitude and characteristic frequency to observed stars.

Large time behavior of solutions to a diffusion approximation radiation hydrodynamics model

This paper concerns the large time behavior of solutions to a diffusion approximation radiation hydrodynamics model when the initial data is a small perturbation around an equilibrium state. The global-in-time well-posedness of solutions is achieved in Sobolev spaces depending on the Littlewood-Paley decomposition technique together with certain elaborate energy estimates in frequency space. Moreover, the optimal decay rates of solutions are also derived provided the initial data satisfy an additional L1 condition. Meanwhile, the same decay rates of the solutions to the approximation system without thermal conductivity could also be established.