Modeling of Protostellar Clouds and their Observational Properties

A physical model and two-dimensional numerical method for computing the evolution and spectra of protostellar clouds are described. The physical model is based on a system of magneto-gasdynamical equations, including ohmic and ambipolar diffusion, and a scheme for calculating the thermal and ionization structure of a cloud. The dust and gas temperatures are determined during the calculations of the thermal structure of the cloud. The results of computing the dynamical and thermal structure of the cloud are used to model the radiative transfer in continuum and in molecular lines. We presented the results for clouds in hydrostatic and thermal equilibrium. The evolution of a rotating magnetic protostellar cloud starting from a quasi-static state is also considered. Spectral maps for optically thick lines of linear molecules are analyzed. We have shown that the influence of the magnetic field and rotation can lead to a redistribution of angular momentum in the cloud and the formation of a characteristic rotational velocity structure. As a result, the distribution of the velocity centroid of the molecular lines can acquire an hourglass shape. We plan to use the developed program package together with a model for the chemical evolution to interpret and model observed starless and protostellar cores.

💡 Research Summary

The paper presents a comprehensive two‑dimensional numerical framework for modeling the dynamical, thermal, and ionization structure of protostellar clouds, and for synthesizing their observable signatures in both continuum and molecular line emission. The core of the physical model is the set of non‑ideal magnetohydrodynamic (MHD) equations that include explicit terms for Ohmic resistivity and ambipolar diffusion. Ohmic diffusion dominates in the dense, weakly ionized core where electron–ion collisions impede current flow, while ambipolar diffusion governs the coupling between neutral gas and magnetic fields in the more diffuse envelope where ion–neutral drag is the primary transport mechanism. By retaining both terms, the authors capture the realistic transition from a well‑coupled magnetic field in the outer layers to a strongly decoupled field in the central region, a feature that is often omitted in ideal MHD treatments.

Thermal balance is treated self‑consistently for gas and dust. Gas heating and cooling processes include line cooling (CO, H₂O, CS, etc.), collisional ionization, and exothermic chemical reactions; cooling is balanced against cosmic‑ray heating, photoelectric heating on grains, and compressional heating. Dust temperature is obtained from the equilibrium between absorption of external radiation (interstellar radiation field) and re‑emission in the infrared, while gas–dust collisional coupling is explicitly calculated, allowing the two temperatures to converge in high‑density regions. The ionization structure is derived from a reduced chemical network that accounts for cosmic‑ray ionization, UV photoionization, and thermal ionization, providing electron and ion abundances needed for the non‑ideal MHD terms.

The dynamical and thermal solutions are then fed into radiative transfer modules. Continuum radiative transfer is solved using a 2‑D ray‑tracing scheme that accounts for anisotropic dust opacity, while molecular line transfer employs a non‑local large‑velocity‑gradient (LVG) approximation adapted to the 2‑D velocity field. This dual approach enables the generation of synthetic spectral maps for both optically thick and thin transitions of linear molecules (e.g., CO J = 2‑1, HCO⁺ J = 1‑0). The authors demonstrate that, in rotating, magnetized clouds, the velocity centroid maps acquire a characteristic hourglass shape: the central region shows a reversal of the line‑of‑sight velocity due to magnetic braking, while the outer parts retain the original rotational signature. This morphology is absent in purely hydrostatic or purely rotating models, highlighting the diagnostic power of combined kinematic and magnetic effects.

Two main classes of models are explored. The first is a static, thermally balanced sphere that serves as a benchmark; it reproduces the expected spherical symmetry in density, temperature, and line profiles. The second class follows the evolution of a quasi‑static, rotating, magnetized cloud. Initially, the cloud is near hydrostatic equilibrium, but as ambipolar diffusion allows magnetic flux to slip outward, angular momentum is redistributed outward by magnetic braking. The core contracts, flattens into a disk‑like structure, and the magnetic field lines become pinched, forming an hourglass configuration. The resulting temperature distribution shows dust temperatures of ~12 K and gas temperatures of ~10 K in the dense core, with a gradual rise toward the envelope. Synthetic line maps reveal broadened, asymmetric profiles in the inner region and narrower, symmetric profiles outward, matching observations of starless cores and early Class 0 objects.

The authors conclude that their integrated modeling platform, which couples non‑ideal MHD, detailed thermal balance, and radiative transfer, can reproduce key observational signatures of starless and protostellar cores. They outline a future program to incorporate a full time‑dependent chemical evolution module, enabling direct comparison with high‑resolution interferometric data from facilities such as ALMA and NOEMA. This will allow the community to constrain the relative importance of magnetic fields, rotation, and chemistry in the earliest phases of star formation, bridging the gap between theory and observation.