The Fast & Fourier: Using neural operators to speed up radiative modeling of protoplanetary disks
Abstract
Edge-on protoplanetary disks provide a spatially resolved view of the vertical distribution of icy planetary building blocks, but inferring their relative abundances requires radiative-transfer (RT) modeling with multiple anisotropic scattering. Each model costs 101-103 CPU hours, which precludes principled Bayesian inference and homogeneous, population-scale analyses. We develop and validate a fast and reusable surrogate for Bayesian inference with ice-bearing edge-on disk models. We train a Fourier Neural Operator on a grid of RADMC-3D spectra. The network maps disk parameters and analytic physical input fields to the spectrum, using mixed-resolution spectra. We use the surrogate in an explicit-likelihood nested sampler and validate it on held-out simulations and synthetic recovery tests with TARP coverage. Finally, we use it to retrieve the disk structure and ice abundances of ESO-Hα 569 from its spectral energy distribution. The surrogate reproduces the RADMC-3D spectra with a median error of 4.84%. We recover synthetic parameters without detectable bias and with nominal coverage. For ESO-Hα 569, the posterior favors a massive (Mgas=0.028+0.015−0.011 M⊙), nearly edge-on (i=79.8+5.2−2.4 deg) disk consistent with the image-plus-SED model of Wolff et al. (2017). Its ice abundances and cavity size remain weakly constrained by broadband data alone. The marginal spectrum-evaluation speed-up is ∼109-1010. Neural-operator emulation makes full posterior inference feasible within the adopted model and training prior. Key posterior points should still be checked with high-fidelity RT models. By drastically cutting evaluation time, this approach enables for the first time homogeneous population studies of ices in disks, which new JWST surveys can capitalize on to map the range of compositions within planetary nurseries.
Related
- The Fast and Fourier: Using neural operators to speed up radiative modeling of protoplanetary disks
- 3D radiative transfer modeling of scattering polarization with partial frequency redistribution. I. Verification and disk-center results for the solar Ca i line at 4227,Å
- Molecules in the shadows of the protoplanetary disk HD 143006: Evidence for C/O ≳ 1 from simultaneously modeling fourteen ALMA lines
- Three-dimensional Radiative-magnetohydrodynamic Modeling of Particle-beam Heating of the Solar Atmosphere
- Determining elemental composition in laboratory meteorite ablation spectra through radiative transfer modeling
BibTeX
@misc{welzel2026fast,
title = {{The Fast & Fourier: Using neural operators to speed up radiative modeling of protoplanetary disks}},
author = {Lukas Welzel and Melissa K. McClure},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.12632},
eprint = {2610.12632},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.12632},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - The Fast & Fourier: Using neural operators to speed up radiative modeling of protoplanetary disks
AU - Lukas Welzel
AU - Melissa K. McClure
PY - 2026
DA - 2026/10/08
PB - arXiv
DO - 10.48550/arxiv.2610.12632
UR - https://arxiv.org/abs/2610.12632
AB - Edge-on protoplanetary disks provide a spatially resolved view of the vertical distribution of icy planetary building blocks, but inferring their relative abundances requires radiative-transfer (RT) modeling with multiple anisotropic scattering. Each model costs 10^1-10^3 CPU hours, which precludes principled Bayesian inference and homogeneous, population-scale analyses. We develop and validate a fast and reusable surrogate for Bayesian inference with ice-bearing edge-on disk models. We train a Fourier Neural Operator on a grid of RADMC-3D spectra. The network maps disk parameters and analytic physical input fields to the spectrum, using mixed-resolution spectra. We use the surrogate in an explicit-likelihood nested sampler and validate it on held-out simulations and synthetic recovery tests with TARP coverage. Finally, we use it to retrieve the disk structure and ice abundances of ESO-Hα 569 from its spectral energy distribution. The surrogate reproduces the RADMC-3D spectra with a median error of 4.84%. We recover synthetic parameters without detectable bias and with nominal coverage. For ESO-Hα 569, the posterior favors a massive (M_gas=0.028^+0.015_−0.011 M_⊙), nearly edge-on (i=79.8^+5.2_−2.4 deg) disk consistent with the image-plus-SED model of Wolff et al. (2017). Its ice abundances and cavity size remain weakly constrained by broadband data alone. The marginal spectrum-evaluation speed-up is ∼10^9-10^10. Neural-operator emulation makes full posterior inference feasible within the adopted model and training prior. Key posterior points should still be checked with high-fidelity RT models. By drastically cutting evaluation time, this approach enables for the first time homogeneous population studies of ices in disks, which new JWST surveys can capitalize on to map the range of compositions within planetary nurseries.
ER -