---
title: 'SEDONA-GesaRaT: 3D NLTE Supernova Modeling'
url: https://www.emergentmind.com/topics/sedona-gesarat
type: topic
---

# SEDONA-GesaRaT: 3D NLTE Supernova Modeling

Searching arXiv for the core SEDONA-GesaRaT paper and closely related SEDONA methodology papers.
Searching arXiv for related work on the integral-based technique and SEDONA-based multidimensional radiative transfer.
SEDONA-GesaRaT is a SEDONA-based radiative-transfer framework for supernova modeling in which the dominant computational bottleneck of non-local thermodynamic equilibrium atomic physics is replaced by a learned surrogate, the atomic physics neural network (APNN). In the reported implementation, the framework retains SEDONA’s Monte Carlo transport machinery, couples it to the integral-based technique (IBT) for spectropolarimetric extraction, and is demonstrated on the three-dimensional Type Ia supernova model N100, including spatially resolved linear-polarization data cubes. Its central claim is practical rather than terminological: direct \(3\)-D, time-dependent, NLTE spectropolarimetric simulation becomes feasible at a total cost of about \(\sim 3000\) core-hours, where earlier workflows were effectively restricted to either \(1\)-D NLTE or \(3\)-D LTE calculations at similar computational budgets [2507.11767].

## 1. Definition and scientific setting

SEDONA-GesaRaT is presented as a modified SEDONA framework rather than a wholly separate transport code. The scientific problem it addresses is the long-standing tradeoff in supernova radiative transfer between multidimensional geometry and realistic atomic kinetics. In \(3\)-D Monte Carlo transport, the cost of packet propagation is already large; in NLTE, the additional requirement to solve the atomic level-population problem zone by zone and timestep by timestep raises the cost sharply. The paper states that, in SEDONA, enabling the traditional NLTE solver for selected elements increases computation cost by roughly a factor of \(\sim 100\) relative to LTE. This matters because supernova remnants, spectra, and spectropolarimetry all indicate intrinsically asymmetric ejecta, while \(1\)-D NLTE and \(3\)-D LTE each discard one of the two essential ingredients: realistic atomic physics or realistic geometry [2507.11767].

Within that context, SEDONA-GesaRaT is defined by an internal substitution. The transport solver, packet propagation, and interaction treatment remain those of SEDONA, but the costly atomic-physics stage is replaced by APNN, a neural-network surrogate trained on \(119\) one-dimensional SN Ia radiative-transfer simulations. The framework is therefore aimed at accelerating the NLTE material-property evaluation that supplies the transport coefficients, not at replacing Monte Carlo transport itself. This distinction is central to its identity.

## 2. Framework composition and workflow

The framework consists of three functional layers: the original SEDONA transport engine, APNN as the atomic-physics surrogate, and IBT as the spectropolarimetric extraction method. The resulting workflow is local-condition \(\rightarrow\) transport-coefficient prediction \(\rightarrow\) Monte Carlo propagation \(\rightarrow\) low-noise observer extraction.

| Component | Function | Basis |
|---|---|---|
| SEDONA | Time-dependent multidimensional Monte Carlo radiative transfer | [2507.11767] |
| APNN | Predicts \(\bar{k}\), \(\bar{j}\), and \(\bar{k}_e\) from local state variables | [2507.11767] |
| IBT | High-S/N spectropolarimetric extraction using integral packets | [2409.07729] |

Operationally, ejecta structure is supplied on a grid, radioactive decay deposition and packet transport proceed as in SEDONA, and each zone and timestep provide the local radiation estimator \(\bar{J}_\nu\), density, and composition. Instead of explicitly solving the rate-equation problem for atomic level populations and then recomputing emissivities and opacities, APNN maps local conditions directly to the absorption or extinction coefficient spectrum \(\bar{k}\), the emissivity spectrum \(\bar{j}\), and the electron-scattering coefficient \(\bar{k}_e\) or \(\bar{\sigma}\). Monte Carlo transport then uses those coefficients exactly where the original SEDONA solver would have used conventionally computed NLTE quantities [2507.11767].

A common misconception is that SEDONA-GesaRaT is an AI replacement for radiative transfer in its entirety. The paper is explicit that this is not the case. The machine-learning component is inserted inside the source-function and atomic-property calculation stage, whereas transport itself remains the original SEDONA Monte Carlo scheme. This preserves the transport formalism while relocating the computational savings to the atomic microphysics.

## 3. APNN and the learned NLTE surrogate

The expensive physical task that APNN bypasses is the determination of atomic level populations under statistical equilibrium. Rather than solving directly for \(n_i\) from the coupled transition-rate system, the network learns an operational mapping
\[
(\rho,\ \text{composition},\ \bar{J}_\nu)\longmapsto(\bar{k}_\nu,\ \bar{j}_\nu,\ \bar{k}_e).
\]
This reframes NLTE atomic physics as a supervised regression problem whose outputs are the transport coefficients actually required by SEDONA [2507.11767].

The network inputs are a rescaled density, a normalized mean-intensity spectrum, a rescaled integrated radiation energy, and a \(28\)-element abundance vector from H to Ni. The appendix states these explicitly as \(\rho_{\rm new}\), \(\bar{J}_{\rm nu}\), \(\bar{E}\), and \(N_{\rm elem}\). The outputs are not trained in raw form. Because \(\bar{k}\) and \(\bar{j}\) span many orders of magnitude, each spectrum is decomposed into a minimum, a maximum, and a normalized spectral shape; electron scattering is predicted through a separate rescaled scalar \(\sigma_e\). This output factorization is intended to compress dynamic range and separate amplitude from morphology.

Architecturally, APNN is a convolutional neural network with \(33\) convolution layers, \(5\) fully connected layers, and \(13{,}480{,}395\) trainable parameters. The \(2305\)-point input spectrum \(\bar{J}_{\rm nu}\) passes through four down-sampling blocks, each block containing four sub-blocks of convolution, layer normalization, and ReLU followed by max pooling. The encoded spectral features are concatenated with \(N_{\rm elem}\), \(\bar{E}\), and \(\rho_{\rm new}\), passed through the fully connected stack to predict the scale parameters \(j_{\rm hi}, j_{\rm lo}, k_{\rm hi}, k_{\rm lo}\), and then decoded through four up-sampling blocks plus a convolution layer to reconstruct the spectral shapes \(j_{\rm nu}\) and \(k_{\rm nu}\). The design is therefore encoder-decoder-like and tailored to structured frequency-dependent outputs.

Training data are generated from \(119\) one-dimensional SN Ia SEDONA simulations using HESMA ejecta models, evolved from day \(1\) to day \(60\). At each radial zone and timestep, the recorded quantities are elemental abundances, plasma density, \(\bar{J}_\nu\), extinction coefficients, and emissivities. The frequency grid spans \(5\times10^{13}\) to \(5\times10^{15}\,\mathrm{Hz}\) with logarithmic spacing \(\Delta\nu/\nu=0.002\), giving \(2305\) frequency points. Four atomic-physics targets are used: CMF\_100 LTE, CMF\_All Ca NLTE, CMF\_All Si, Ca NLTE, and CMF\_All Si, S, Ca NLTE. The \(119\) models are split into \(80\) training, \(20\) validation, and \(19\) testing cases. Training uses Adam, mean squared error loss, batch size \(64\), learning rate \(10^{-4}\), and \(200\) epochs; training time is about \(2\) days on an NVIDIA A100 GPU, with implementation in PyTorch and deployment through LibTorch [2507.11767].

The paper emphasizes that APNN is not universal. It is trained on SN Ia ejecta only, tied to the specific \(2305\)-point frequency grid, and learns only the atomic recipes represented in its training set. This restricts immediate portability, but also defines the domain within which the reported acceleration claims are meaningful.

## 4. Transport formalism, polarization, and IBT

SEDONA-GesaRaT retains SEDONA’s time-dependent Monte Carlo transport for supernova ejecta. The simulations include gamma-ray emission from radioactive isotopes, gamma-ray interactions through Compton scattering and photoionization or photoelectric processes, and optical, UV, and IR transport with Thomson scattering, bound-bound transitions, bound-free transitions, and free-free transitions. The ejecta are assumed to undergo homologous expansion, while line transport uses the Sobolev approximation and expansion-opacity formalism [2507.11767].

Several transport coefficients and interaction probabilities are stated explicitly. The comoving-frame Thomson extinction coefficient is
\[
\bar{k}_e=\sigma_T n_e,
\]
and the Sobolev optical depth for an individual line is
\[
\tau_{\rm Sob}=\frac{h}{4\pi}\left(N_lB_{lu}-N_uB_{ul}\right)t_{\rm exp}c.
\]
The effective interaction probability is written as
\[
\tau_{\rm eff}=1-e^{-\tau_{\rm Sob}},
\]
and the total bound-bound extinction in a frequency bin is
\[
\bar{k}_{bb}(\bar{\nu},\bar{\nu}+\Delta\bar{\nu})=
\frac{(\bar{\nu}+0.5\Delta\bar{\nu})\sum_i\tau_{{\rm eff},i}}
{\Delta\bar{\nu}\,t_{\rm exp}\,c}.
\]
These expressions clarify that the framework accelerates coefficient generation but does not alter the radiative-transfer formalism itself.

For spectropolarimetry, the reported observables are the standard linear-polarization fractions \(Q/I\) and \(U/I\). The companion methodological paper on IBT defines it as an observer-extraction algorithm based on integral packets rather than event-spawned virtual packets. Implemented in \(3\)-D SEDONA, IBT was reported to be faster than the event-based technique by a factor of \(10\)–\(30\) and to improve signal-to-noise by a factor of \(5\)–\(10\) at fixed Monte Carlo quanta, while also passing spherical-symmetry, mirror-symmetry, and cross-comparison tests against direct counting and event-based extraction [2409.07729]. In SEDONA-GesaRaT, IBT is the mechanism that makes the low-noise spatially resolved polarization products practical once the atomic-physics bottleneck is removed.

## 5. Validation, N100 application, and computational performance

Validation is reported at two levels. First, after APNN is integrated into SEDONA-GesaRaT, the authors run one-dimensional time-dependent tests on the held-out \(19\)-model test set and report that the spectral time series reproduce the original SEDONA results from early to late phases with high accuracy, with predicted and conventional absorption or emission spectra overlapping visually. Second, they compare SEDONA-GesaRaT against original SEDONA in \(3\)-D using the simpler CMF\_100 LTE setup on the N100 model, because a direct \(3\)-D NLTE-to-\(3\)-D NLTE comparison with the traditional solver would be prohibitively expensive [2507.11767].

For the \(3\)-D LTE benchmark, the setup is a \(50\times50\times50\) Cartesian grid, \(4\times10^6\) packets per timestep, evolution from day \(1\) to day \(60\), and execution on a \(48\)-core node. The reported wall-clock times are \(\sim 73.3\) hr for original SEDONA and \(\sim 56.3\) hr for SEDONA-GesaRaT. Flux spectra and linear polarization are generally consistent from \(2000\)–\(50000\) Å and from \(5\)–\(50\) days, with the main discrepancies occurring in Si II \(5640\) Å around day \(17\) and the Ca II IR triplet around day \(41\).

The paper’s principal scientific demonstration is the \(3\)-D NLTE application to the delayed-detonation SN Ia model N100. The ejecta are mapped to a \(50^3\) Cartesian grid with maximum ejecta velocity \(28600\ {\rm km\ s^{-1}}\). The showcased viewing direction is
\[
\mu=0.5,\qquad \varphi=3.4455.
\]
The framework computes time-dependent \(3\)-D spectra, linear spectropolarimetry, and spatially resolved spectropolarimetry data cubes. The resolved image calculation uses \(4\times10^6\) packets per timestep and a \(20\times20\) observer image grid spanning
\[
-32000\,t_{\rm exp}\ {\rm km\,s^{-1}} \quad \text{to} \quad +32000\,t_{\rm exp}\ {\rm km\,s^{-1}}
\]
in both projected axes. The total runtime for this \(3\)-D NLTE spectropolarimetry simulation is \(63\) hours on a \(48\)-core node, or \(3024\) core-hours [2507.11767].

The physical outputs are not limited to global spectra. At \(17\)–\(18\) days after explosion, putting S in NLTE enhances the emission feature of S II \(6715\) Å and the absorption features of S II \(5454\) Å and \(5640\) Å. Putting Si in NLTE enhances the P-Cygni profile of Si II \(6347/6371\) Å, enhances the absorption feature of Si II \(4128\) Å, and suppresses the linear polarization signal near \(6100\) Å. The spatially resolved polarization maps at \(3697\) Å and \(3800\) Å are interpreted as tracing Ca II H\&K line formation in velocity slices of roughly \(17000\)–\(20000\ {\rm km\ s^{-1}}\) and \(10000\)–\(13000\ {\rm km\ s^{-1}}\), respectively, and their morphology broadly overlaps the Ca density distribution in those slices.

The per-zone timing benchmark summarizes the acceleration more starkly:

| Atomic treatment | Core-seconds per zone |
|---|---:|
| CMF\_100 LTE | 1.19 |
| CMF\_All Ca NLTE | 57.92 |
| CMF\_All Si, Ca NLTE | 65.86 |
| CMF\_All Si, S, Ca NLTE | 72.09 |
| APNN | 0.17 |

These figures imply that APNN is about \(7\) times faster than the reduced LTE calculation and roughly \(341\)–\(424\) times faster than the explicit NLTE variants. The paper summarizes this as a reduction from about \(\sim 50\) core-seconds to \(0.17\) core-seconds per zone, that is, roughly two to three orders of magnitude [2507.11767].

## 6. Limits, interpretation, and research role

The principal limitations are sharply defined. APNN is trained only on SN Ia models; application to core-collapse supernovae would require new training data. The learned surrogate is tied to a specific \(2305\)-point frequency grid and to selected atomic recipes rather than a general all-species NLTE treatment. Residual feature-level inaccuracies remain, particularly in S II \(5640\) Å and parts of the Ca II IR triplet. The transport still relies on the expansion-opacity formalism, whose reliability degrades when a bin contains multiple optically thick lines. More broadly, improvements may require better atomic data libraries, more stable or more complete NLTE solvers, more diverse training sets, and refined network architectures [2507.11767].

These limits also define the correct interpretation of the framework. SEDONA-GesaRaT is not a proof that NLTE atomic physics has been solved generically by machine learning; it is a domain-specific acceleration strategy for a constrained but astrophysically important class of SN Ia problems. Its significance lies in changing the feasible scale of multidimensional radiative-transfer studies. The paper explicitly connects the reported efficiency to systematic surveys of \(3\)-D explosion models, parameter studies over viewing angle, composition, and asymmetry, model libraries for interpreting observed spectropolarimetry, and integration with earlier AI-assisted inversion approaches. A plausible implication is that the framework repositions detailed \(3\)-D NLTE spectropolarimetry from a boutique calculation to a surveyable regime, but only within the bounds of its training domain and transport approximations.

In that sense, SEDONA-GesaRaT occupies a specific place in the SEDONA lineage. It leaves the Monte Carlo transport core intact, incorporates IBT for practical observer extraction, and relocates the dominant computational burden from direct atomic solving to supervised surrogate evaluation. The result is a SEDONA derivative whose novelty is architectural and computational: direct \(3\)-D NLTE spectropolarimetric simulation of the N100 model, including spatially resolved linear polarization data cubes, at a reported cost of \(\sim 3000\) core-hours [2507.11767].

Source: https://www.emergentmind.com/topics/sedona-gesarat