A computed line list for the H2D+ molecular ion

A comprehensive, calculated line list of frequencies and transition probabilities for the singly deuterated isotopologue of H3+, H2D+, is presented. The line list, called ST1, contains over 22 million rotational-vibrational transitions occurring between more than 33 thousand energy levels; it covers frequencies up to 18500 cm-1. All energy levels with rotational quantum number, J, up to 20 are considered, making the line list useful for temperatures up to at least 3000 K. About 15% of these levels are fully assigned with approximate rotational and vibrational quantum numbers. The list is calculated using a previously proposed, high accuracy, ab initio model and consistency checks are carried out to test and validate the results. These checks confirm the accuracy of the list. A temperature-dependent partition function, valid over a more extended temperature range than those previously published, and cooling function are presented. Temperature-dependent synthetic spectra in the frequency range 0 - 10000 cm-1 are also given.

💡 Research Summary

This paper presents ST1, a comprehensive computed line list for the singly deuterated isotopologue of H₃⁺, H₂D⁺. The authors generated over 22 million rotational‑vibrational transitions linking more than 33 000 energy levels, covering frequencies up to 18 500 cm⁻¹ and rotational quantum numbers J ≤ 20. Approximately 15 % of the levels are fully assigned with approximate vibrational (v₁, v₂, v₃) and rotational (J, Kₐ, K_c) quantum numbers; the remainder are labeled only by J and symmetry (ortho/para).

The calculations were performed with the DVR3D suite, which solves the exact nuclear‑motion Hamiltonian within the Born‑Oppenheimer approximation. The underlying potential energy surface (PES) is the global H₃⁺ surface of Polyansky et al. (2000), itself based on ultra‑high‑accuracy ab‑initio points (Cencek et al., 1998) and supplemented by high‑energy data. An adiabatic correction specific to H₂D⁺ (Polyansky & Tennyson, 1999) and vibrational‑mass scaling were applied to mitigate non‑adiabatic effects (μ_H = 1.0075372 u, μ_D = 2.0138140 u for vibrations; μ_H = 1.00727647 u, μ_D = 2.01355320 u for rotations). The nuclear‑motion problem was expressed in Jacobi coordinates (r₁, r₂, θ) with Morse‑type radial functions for the H–H bond and spherical oscillators for the D‑center‑of‑mass distance. Convergence tests (e.g., J = 3 and J = 15) ensured that level energies are accurate to better than 0.01 cm⁻¹ for essentially all states considered.

Transition intensities were obtained using the ab‑initio dipole moment surface of Röthse et al. (1994), evaluated with a 50‑point Gauss‑Legendre quadrature. The final ST1 data set consists of two files (levels and transitions) formatted analogously to the BT2 water line list, facilitating immediate integration into existing astrophysical radiative‑transfer codes. The total computational effort amounted to roughly 8 000 CPU hours on 32‑ and 64‑bit Linux machines.

Validation was carried out against a broad set of laboratory measurements (Shy et al., 1981; Amano & Watson, 1984; Főster et al., 1986; Fárník et al., 2002; Asvany et al., 2007). Statistical analysis of 9–73 lines per source shows mean frequency deviations ranging from –0.014 to +0.242 cm⁻¹ with standard deviations of 0.02–0.09 cm⁻¹, indicating a substantial improvement over earlier theoretical lists (e.g., Miller et al., 1989). Moreover, Einstein A‑coefficients derived from ST1 reproduce the relative Einstein B‑coefficients measured by Asvany et al. within experimental uncertainties, confirming the reliability of the computed line strengths.

Using the ST1 energy levels, the authors computed the temperature‑dependent partition function Q(T) for 5 K ≤ T ≤ 4000 K. Comparison with the classic Sidhu et al. (1992) values shows excellent agreement below ~1200 K, but significant divergence at higher temperatures where the inclusion of many high‑lying levels in ST1 raises Q(T) appreciably. A fifth‑order polynomial fit to Q(T) is provided for convenient use in modeling. A cooling function, derived from the summed spontaneous emission rates, is also presented, highlighting H₂D⁺’s potential role in the thermal evolution of primordial gas clouds.

Synthetic spectra were generated for the 0–10 000 cm⁻¹ range at several temperatures, demonstrating the appearance of strong rotational lines (e.g., the 372 GHz 1₁₀–1₁₁ transition) and vibrational bands (ν₁, ν₂, ν₃). These spectra are directly comparable with observations from facilities such as ALMA, Herschel, and future far‑infrared missions, enabling the use of H₂D⁺ as a diagnostic of cold, dense interstellar environments, deuterium fractionation, and even as a probe of dark‑matter clumps.

In summary, the ST1 line list provides an unprecedentedly complete and accurate spectroscopic dataset for H₂D⁺, extending the applicability of H₂D⁺ modeling to temperatures up to at least 3000 K. Its high‑quality frequencies, Einstein A‑coefficients, partition function, and cooling function make it a valuable resource for astrochemical networks, radiative‑transfer simulations, and the interpretation of both current and forthcoming astronomical observations. Future work may expand this approach to other H₃⁺ isotopologues (e.g., D₂H⁺) and incorporate further experimental refinements, thereby establishing a comprehensive spectroscopic foundation for the entire H₃⁺ family.