arxiv: 1907.03367 · v1 · pith:SN5ZQP6Bnew · submitted 2019-07-07 · 🌌 astro-ph.EP · astro-ph.SR

ExoMol line list XXXIV: A Rovibrational Line List for Phosphinidene (PH) in its X\,{}³Sigma^- and a\,{}¹Delta Electronic States

Jonathan Langleben , Jonathan Tennyson , Sergei N. Yurchenko , Peter Bernath This is my paper

Pith reviewed 2026-05-25 00:58 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.SR

keywords phosphinidenePHrovibrational line listpotential energy curvesdipole moment curvesexoplanet atmospheresmolecular spectroscopy

0

0 comments

The pith

A rovibronic line list for phosphinidene is built from fitted potential curves and ab initio dipoles for use in hot atmospheres.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper produces a line list for the X 3Σ− ground state and a 1Δ excited state of 31PH by fitting potential energy curves to measured transition frequencies. These curves, together with Born-Oppenheimer and spin-splitting corrections, are inserted into the nuclear Schrödinger equation to obtain 2528 rovibronic states. Line intensities are then obtained from separate ab initio dipole-moment curves. The resulting LaTY list supplies 65055 transitions up to 24500 cm−1 and J = 80 and is shown to generate emission and absorption spectra at temperatures up to 4000 K.

Core claim

Empirical potential energy curves fitted to laboratory transition frequencies reproduce those frequencies with a root-mean-square error of 0.01 cm−1; when combined with ab initio dipole-moment curves the same curves yield a complete set of Einstein A coefficients for the X–X and a–a systems of 31PH.

What carries the argument

The LaTY line list, generated by solving the nuclear Schrödinger equation on empirically adjusted PECs plus correction terms and scaling the resulting wavefunctions by ab initio dipole-moment curves.

If this is right

Emission and absorption spectra of PH can be simulated directly at any temperature up to 4000 K.
The list supplies complete data for both the triplet ground state and the singlet excited state through J = 80.
Transitions are provided up to 24500 cm−1, covering the wavelength range needed for cool stellar and exoplanet atmosphere models.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same fitting procedure could be applied to other isotopologues once laboratory frequencies become available.
Detection of PH features in exoplanet transmission spectra would become feasible once the line list is incorporated into atmospheric retrieval codes.

Load-bearing premise

Ab initio dipole-moment curves, once paired with the fitted potentials, give intensities that are accurate enough for astrophysical modeling without further empirical correction.

What would settle it

Laboratory or astronomical spectra whose measured absolute line intensities differ by more than the expected uncertainty from the intensities computed in the LaTY list.

read the original abstract

A rovibronic line list for the ground ($X$ $^3\Sigma^-$) and first excited ($a$ $^1\Delta$) states of phosphinidene, $^{31}$PH, is computed. The line list is designed for studies of exoplanetary and cool stellar atmospheres with temperatures up to 4000 K. A combination of empirical and ab initio data are used to produce the line list: potential energy curves (PECs) are fitted using experimental transition frequencies; these transitions are reproduced with a root mean square error of 0.01 cm$^{-1}$. The nuclear Schr\"{o}dinger is solved using these PECs plus Born-Oppenheimer and spin splitting correction terms. Line intensities and Einstein $A$ coefficients are computed using ab initio Dipole Moment Curves (DMC) $X$-$X$ and $a$-$a$. The resulting LaTY line list, which contains 65055 transitions for 2528 rovibronic states up to 24500 cm$^{-1}$ and $J=80$ is used to simulate spectra in emission and absorption at a range of temperatures. The line list is made available in electronic form at the CDS and ExoMol databases.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Standard ExoMol line list for PH with tight frequency fit but intensities left on untested ab initio DMCs.

read the letter

The paper delivers a new line list for 31PH in the X 3Σ− and a 1Δ states, with 65055 transitions up to J=80 and 24500 cm−1. Frequencies come from fitted PECs plus corrections and match experiment to 0.01 cm−1 RMS. That part is straightforward and useful for anyone who needs positions for PH in cool atmospheres up to 4000 K. The list is deposited in the usual ExoMol and CDS places, so it slots directly into existing tools. What they do well is stick to the series method without claiming more than the data support. The nuclear Schrödinger solution and the coverage look complete for the stated range. The soft spot is the intensities. Einstein A values rest on ab initio X-X and a-a dipole curves with no reported comparison to measured strengths, no scaling, and no cross-checks against other electronic-structure methods. The stress-test note is right on this: the energy levels are anchored, but the absolute scale for opacities is an assumption. For applications that only need line positions this is minor; for retrievals that weight line strengths it is the main uncertainty. This is for molecular spectroscopists and exoplanet atmosphere modelers who already use ExoMol data and now have one more species covered. It is not a broad methodological advance, but the specific data product fills a documented gap. The work shows clear, reproducible steps and honest engagement with the literature on this narrow task. I would send it to peer review so referees can check the intensity section and any supplementary validation tables.

Referee Report

1 major / 0 minor

Summary. The manuscript computes a rovibronic line list (LaTY) for 31PH in the X 3Σ− and a 1Δ states for exoplanet and cool stellar atmosphere modeling up to 4000 K. Potential energy curves are fitted to experimental transition frequencies (RMS error 0.01 cm−1); the nuclear Schrödinger equation is solved with these PECs plus Born-Oppenheimer and spin-splitting corrections. Einstein A coefficients and intensities are obtained from separate ab initio dipole moment curves for X-X and a-a transitions. The list contains 65055 transitions for 2528 states up to 24500 cm−1 and J=80; simulated emission and absorption spectra are presented. Data are deposited at CDS and ExoMol.

Significance. If the intensities hold, the line list would supply a practical tool for phosphorus chemistry in astrophysical environments, extending the ExoMol series. The empirical PEC fitting that reproduces observed frequencies to 0.01 cm−1 RMS is a clear technical strength and supports reliable line positions.

major comments (1)

[Abstract and intensity-calculation section] Abstract and intensity-calculation section: intensities and Einstein A coefficients are derived solely from ab initio DMCs with no reported comparison to measured line strengths, no empirical scaling, no cross-method benchmarks, and no uncertainty estimates, even though the list is explicitly intended for quantitative modeling at 4000 K. Frequency accuracy alone does not establish the intensity scale required by the central claim.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful review and for highlighting the strength of the empirical PEC fitting. We address the single major comment below.

read point-by-point responses

Referee: [Abstract and intensity-calculation section] Abstract and intensity-calculation section: intensities and Einstein A coefficients are derived solely from ab initio DMCs with no reported comparison to measured line strengths, no empirical scaling, no cross-method benchmarks, and no uncertainty estimates, even though the list is explicitly intended for quantitative modeling at 4000 K. Frequency accuracy alone does not establish the intensity scale required by the central claim.

Authors: We agree that the manuscript reports no experimental line-strength comparisons, empirical scaling, or explicit uncertainty estimates for the intensities. No measured intensities for the X–X or a–a transitions of PH appear in the literature, precluding direct validation. The DMCs were obtained at the MRCI/aug-cc-pV5Z level. In the revised manuscript we will add a dedicated paragraph in the intensity section that (i) states the absence of experimental benchmarks, (ii) references accuracy assessments of the same level of theory for chemically similar diatomics (e.g., NH, OH), and (iii) supplies conservative uncertainty estimates (approximately 10–20 % for strong lines) based on those benchmarks. We will also revise the abstract to indicate that intensities are the best currently available from theory rather than claiming quantitative accuracy without qualification. revision: yes

Circularity Check

0 steps flagged

No circularity: frequencies fitted to external data, intensities from independent ab initio DMCs

full rationale

The derivation fits PECs to experimental transition frequencies (RMSE 0.01 cm^{-1}), solves the nuclear Schrödinger equation with those PECs plus corrections, and computes Einstein A coefficients from separate ab initio X-X and a-a DMCs. No step reduces by construction to its own inputs, no fitted parameter is relabeled as a prediction, and no load-bearing claim rests on self-citation chains or imported uniqueness theorems. The line list is a direct forward computation from these distinct empirical and ab initio sources; the central output does not tautologically reproduce the fit inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The claim rests on fitting a small number of parameters in the PECs to experimental frequencies, standard quantum-mechanical solution of the nuclear Schrödinger equation, and the domain assumption that ab initio DMCs give usable intensities. No new particles or forces are introduced.

free parameters (1)

PEC fitting parameters
Potential energy curves are fitted to experimental transition frequencies; exact number and functional form not stated in abstract

axioms (2)

domain assumption Born-Oppenheimer approximation with added correction terms suffices for the nuclear motion
Invoked when solving the nuclear Schrödinger equation after fitting PECs
domain assumption Ab initio dipole moment curves yield intensities accurate for the target temperature range
Used to compute Einstein A coefficients and line intensities

pith-pipeline@v0.9.0 · 5790 in / 1411 out tokens · 33635 ms · 2026-05-25T00:58:52.818431+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.