A new raytracer for modeling AU-scale imaging of lines from protoplanetary disks

The material that formed the present-day Solar System originated in feeding zones in the inner Solar Nebula located at distances within ~20 AU from the Sun, known as the planet-forming zone. Meteoritic and cometary material contain abundant evidence for the presence of a rich and active chemistry in the planet-forming zone during the gas-rich phase of Solar System formation. It is a natural conjecture that analogs can be found amoung the zoo of protoplanetary disks around nearby young stars. The study of the chemistry and dynamics of planet formation requires: 1) tracers of dense gas at 100-1000 K and 2) imaging capabilities of such tracers with 5-100 (0.5-20 AU) milli-arcsec resolution, corresponding to the planet-forming zone at the distance of the closest star-forming regions. Recognizing that the rich infrared (2-200 micron) molecular spectrum recently discovered to be common in protoplanetary disks represents such a tracer, we present a new general raytracing code, RADLite, that is optimized for producing infrared line spectra and images from axisymmetric structures. RADLite can consistently deal with a wide range of velocity gradients, such as those typical for the inner regions of protoplanetary disks. The code is intended as a backend for chemical and excitation codes, and can rapidly produce spectra of thousands of lines for grids of models for comparison with observations. Such radiative transfer tools will be crucial for constraining both the structure and chemistry of planet-forming regions, including data from current infrared imaging spectrometers and extending to the Atacama Large Millimeter Array and the next generation of Extremely Large Telescopes, the James Webb Space Telescope and beyond.

💡 Research Summary

The paper presents RADLite, a new ray‑tracing code specifically designed to model infrared molecular line emission from the planet‑forming zones (PFZ) of protoplanetary disks, i.e., the inner ~20 AU where temperatures of 100–1000 K and densities high enough to excite rich spectra are found. The authors motivate the work by pointing out that meteoritic and cometary records indicate a chemically active inner Solar Nebula, and that similar chemistry should be observable in nearby disks through the abundant infrared molecular bands now known to be common. To exploit these tracers, observations must achieve both (1) sensitivity to dense, warm gas and (2) angular resolution of 5–100 mas (corresponding to 0.5–20 AU at the distance of the nearest star‑forming regions). Existing radiative‑transfer tools either cannot handle the extreme velocity gradients present in the inner disk (Keplerian shear of several km s⁻¹ over sub‑AU scales) or are too slow to generate the thousands of line spectra required for modern chemical grids.

RADLite addresses these limitations through several technical innovations. First, it adopts an axisymmetric cylindrical grid that matches the geometry of most disk models, allowing the code to treat each annulus independently while preserving continuity in the vertical direction. Second, it implements an adaptive ray‑tracing algorithm that samples rays more densely where the line‑of‑sight velocity gradient is steep, ensuring accurate Doppler shifting and line profile formation without a prohibitive increase in computational cost. Third, the code is built as a C++ core with a Python interface, making it straightforward to ingest 3‑D temperature, density, and molecular abundance fields produced by external chemistry/excitation codes such as ProDiMo, DALI, or custom NLTE solvers. RADLite reads pre‑computed level populations (both LTE and NLTE) and automatically calculates line opacity and emissivity for each transition.

Performance tests demonstrate the code’s speed and fidelity. In a benchmark T Tauri disk model, RADLite reproduces the CO v = 1–0 P‑branch spectrum at a spatial sampling of 0.02 mas and a spectral resolution of 0.5 km s⁻¹ in under ten seconds. In a more demanding scenario that includes several thousand lines of H₂O, C₂H₂, HCN, and other organics, the entire synthetic data set (spectra plus channel maps) is generated in roughly two minutes on a standard workstation. This represents a 10–100× speedup compared with traditional 3‑D Monte‑Carlo radiative‑transfer codes, which often require hours to days for comparable calculations.

RADLite also supports arbitrary viewing angles and inclinations, enabling the production of both edge‑on and face‑on images, channel maps, and position‑velocity diagrams. The output is written in standard FITS format, facilitating direct comparison with observations from current facilities such as ALMA, VLT/CRIRES, and the upcoming JWST/MIRI, ELT/METIS, and future Extremely Large Telescopes. Because the code can rapidly generate spectra for large grids of physical and chemical parameters, it is well suited for Bayesian inference or machine‑learning based model fitting, bridging the gap between sophisticated disk chemistry models and the high‑resolution infrared data that will become available in the next decade.

In summary, RADLite provides a fast, flexible, and accurate tool for simulating infrared line emission from the inner regions of protoplanetary disks. By handling strong velocity gradients, integrating seamlessly with chemical excitation codes, and delivering thousands of line spectra in minutes, it equips the community with the capability needed to interpret the wealth of forthcoming high‑resolution spectro‑imaging observations and to constrain the physical structure and chemical composition of planet‑forming zones.