Extraction of spin-averaged rovibrational transition frequencies in HD + for the determination of fundamental constants

We present a comprehensive analysis of all currently available high-accuracy frequency measurements of rotational and rovibrational transitions in the hydrogen molecular ion HD + . Our analysis utilises the theoretically calculated hyperfine structure to extract the values of three spin-averaged transition frequencies through a global linear least-squares adjustment that takes into account theory-induced correlations between the different transitions. We subsequently use the three spin-averaged transition frequencies as input data in a second adjustment which employs precise theoretical expressions for the transition frequencies, written as a function of the proton, deuteron and electron relative atomic masses, the Rydberg constant, and the proton and deuteron charge radii. Our analysis shows that the HD + data may significantly improve the value of the electron relative atomic mass and the proton-electron mass ratio, in particular if combined with recent high-precision measurements of particle atomic masses and mass ratios obtained from Penning traps.


Introduction
Fundamental physical constants are determined from precision measurements, often performed on single particles or simple atomic quantum systems, by adjusting the values of the constants such that theoretical predictions match the experimental observations as closely as possible.Such adjustments, implemented via least-squares optimization of the constants of interest with respect to large bodies of input data, are carried out by standards organizations such as the Committee on Data for Science and Technology (CODATA) and the Atomic Mass Data Center (AMDC).from laser spectroscopy of both electronic and muonic hydrogen and deuterium, complemented by high-precision quantum electrodynamics (QED) theoretical predictions (see Ref. [1] and references therein).By contrast, the relative atomic masses of particles like the electron and light nuclei such as the proton and deuteron have traditionally been determined using Penning-trap experiments and (for the electron) theoretical calculations of bound-state g-factors; for recent examples see Refs.[2][3][4][5][6][7].In 1976 Wing et al. already recognised that vibrational spectroscopy of hydrogen molecular ions, such as H + 2 and HD + , could also be used for determinations of particle mass ratios provided that both the measurements and the theory of these three-body, two-center quantum systems would be sufficiently accurate [8].This concept was extended to laser spectroscopy of antiprotonic helium, pHe, by Hori et al. who combined their results with precise ab initio theory for a meaningful determination of the antiproton-electron mass ratio which under the assumption of charge, parity and time reversal (CPT) invariance must be the same as the proton-electron mass ratio, m p /m e [9][10][11].As a result, pHe data were included in the 2006 and 2010 CODATA adjustments where they contributed primarily to the value of the electron relative atomic mass, A r (e) [12,13].However, the 2014 and 2018 CODATA adjustments no longer considered pHe after A r (e) was determined using Penning-trap measurements and QED theory with over an order of magnitude smaller uncertainty [2].
Over the past two decades, advances in laser and terahertz spectroscopy of HD + stored in radiofrequency (rf) ion traps and sympathetically cooled by laser-cooled ions, considerably improved the accuracy, culminating in the measurement of various rovibrational transition frequencies with uncertainties in the low parts-per-trillion (ppt) range [14][15][16].Meanwhile, the precision of theoretical QED predictions of relevant transition frequencies in HD + improved to below 10 ppt, such that the uncertainty of theoretically predicted frequencies is now limited by (primarily) the uncertainty of the 2018 CODATA value of m p /m e [17][18][19].As anticipated by Wing et al. [8] and more recently also by Karr et al. [20], the availability of theoretical and experimental transition frequencies with uncertainties in the low-ppt range enabled a determination of m p /m e with a precision of about 20 ppt [14][15][16].These values represented not only the most precise determinations of m p /m e to date, but they were also found to be consistent with a recent precise measurement of the proton relative atomic mass, A r (p) [3,4], that deviated significantly from determinations of A r (p) included in the 2018 CODATA adjustment [1].Furthermore, several researchers have pointed out that results from HD + and pHe form a link between 'atomic' fundamental constants such as R ∞ , r p and r d , and particle mass ratios [21][22][23].
The recent advances in HD + experiments and theory and their demonstrated potential for the determination of fundamental constants warrant a closer investigation to prepare available HD + data for inclusion in future CODATA adjustments.One crucial issue is that, in general, atomic and molecular spectroscopic data include line shifts due to hyperfine interactions, whereas CODATA adjustments have historically considered only energy differences between hyperfine centroids.This implies that experimental HD + spectroscopic data must undergo a preparatory analysis in order to remove hyperfine shifts from measured rovibrational transition frequencies.Here, we will use the term 'hyperfine components' to indicate rovibrational transitions between individual hyperfine levels.In principle, the spin-averaged transition frequency can be determined from the measured transition frequencies of a sufficient number of hyperfine components.However, in the experiments conducted so far, the number of measured hyperfine components has been too small, and no other independent experimental data on the relevant HD + hyperfine structure exist.To derive spin-averaged transition frequencies from limited sets of measured frequencies of hyperfine components, several methods have been developed which overcome the lack of experimental data by using the theoretically predicted hyperfine level structure [14,24,25].In this work we adopt the approach of [24], which is based on a least-squares optimization and takes into account correlations between theoretical hyperfine coefficients.
This article is organised as follows.The state-of-the-art HD + hyperfine theory used for the extraction of spin-averaged transition frequencies is reviewed in Section 2, along with a discussion of correlated hyperfine theory uncertainties.Section 3 summarises the available experimental data (that include hyperfine shifts) related to three different rotational-vibrational transitions in HD + , followed by a global adjustment of the three spin-averaged rotational-vibrational transition frequencies in Section 4. In Section 5, the adjusted spin-averaged transition frequencies and their covariances are subsequently compared to theoretical expressions for these spin-averaged transition frequencies that are functional expressions of several fundamental constants, including A r (e), A r (p), and the deuteron relative atomic mass, A r (d).We use the corresponding observational equations in a second adjustment, aimed at determining the values of the fundamental constants.We consider various sets of input data and observational equations (including recent Penning-trap measurements of relative atomic masses), and compare values of mass ratios such as m p /m e derived from each adjustment, in order to illustrate the potential impact of HD + on the precision of fundamental particle mass ratios.Conclusions are presented in Section 6.

Hyperfine Hamiltonian
The hyperfine structure of a rovibrational state of HD + can be described by the following effective spin Hamiltonian [26]: where I d , I p , and s e are the spins of the deuteron, proton, and electron, respectively, and L is the total angular momentum excluding electron and nuclear spins (denoted by N in standard spectroscopic notation).Here and throughout the paper, spin and angular momentum operators are written in units of and are thus dimensionless.
The coefficients E k , k = 1, 2, . . ., 9 are then energies and depend on the rovibrational state (v, L), where v stands for the vibrational quantum number.This Hamiltonian includes all spin-dependent interactions that appear at the leading order m e c 2 α 4 , i.e. appear in the Breit-Pauli Hamiltonian, except for the proton-deuteron spin-spin contact interaction, which is negligibly small because of the strong Coulomb repulsion between the proton and deuteron.Additional spin-dependent interactions appear at higher orders (m e c 2 α 6 and above).
The largest missing term in Equation ( 1) is a coupling between proton and deuteron spins (proportional to I p •I d ), mediated by the electron.The coupling coefficient of this interaction was shown to be on the order of 100 Hz [25], which is much smaller than the current theoretical uncertainty of E 4 (see below).The effective spin Hamiltonian may thus be considered complete for our purposes.The largest terms in Equation ( 1) are the electron-proton (E 4 /h ∼ 900 MHz in the vibrational ground state) and electron-deuteron (E 5 /h ∼ 140 MHz) spin-spin contact interactions.They were calculated with relative uncertainties below 1 parts per million (ppm) [27].A key element to reach this accuracy is that the sum of nuclear structure and recoil corrections to E 4 (E 5 ), which we denote as 'nuclear correction' in the following, was determined from the difference between the experimental groundstate hyperfine splitting of the hydrogen (deuterium) atom and the corresponding nonrecoil, point-particle QED prediction, under the assumption that the nuclear correction is entirely described by a contact (delta-function) term.Uncertainties associated with nuclear corrections are thus suppressed, and the theoretical precision is limited by unevaluated higher-order QED terms: (i) nonrecoil corrections, dominated by the contribution of order α(Zα 2 ) ln(Zα)E F [28,29], and (ii) recoil corrections that deviate from the contact-term approximation, i.e. mainly the contribution of order (Zα) 2 (m/M )E F [30].Here, E F stands for the leading Fermi contribution of order m e c 2 α 4 .
The next largest terms are the spin-orbit (E 1 /h ∼ 30 MHz in the vibrational ground state) and spin-spin tensor interactions (E 6 /h ∼ 8 MHz, E 7 /h ∼ 1 MHz for L = 1).They were computed with few-ppm uncertainties through calculation of higherorder corrections at orders m e c 2 α 6 [31] and m e c 2 α 7 ln(α) [32,33].The uncertainty is dominated by unevaluated nonlogarithmic corrections at the m e c 2 α 7 order.
The remaining hyperfine coefficients are much smaller, ranging from E 8 /h ∼ 3 kHz to E 2 /h ∼ 30 kHz.They were calculated in the framework of the Breit-Pauli Hamiltonian accounting for the electron's anomalous magnetic moment [26].In general, the expected relative uncertainty of this approximation is on the order of α 2 , but in this case second-order contributions induced by the leading hyperfine interaction terms can give larger contributions.Due to this, the relative uncertainty of E 2 , E 3 , and E 8 is estimated to be 5α 2 [33].In the case of the E 9 term that describes the effect of the deuteron quadrupole moment, second-order terms are too small to significantly influence the uncertainty.The relative uncertainty of E 9 is estimated to be 98 ppm [33], which accounts for both the theoretical uncertainty and the uncertainty of the deuteron quadrupole moment from [34].
State-of-the art values of all the hyperfine coefficients, for rovibrational levels involved in published high-precision measurements, are given in Table 1.
In our least-squares adjustment of HD + experimental data (see Section 4), theoretical uncertainties in the hyperfine coefficients E k (v, L) are taken into account by introducing additive corrections δE th k (v, L) which are treated as adjusted constants, following the approach of CODATA [1].For each coefficient, an input datum δE k (v, L) is included with a zero value and an uncertainty equal to that given in Table 1, and an observational equation δE k (v, L) .= δE th k (v, L) is added to the adjustment.Beyond the uncertainties themselves, it is also necessary to estimate the correlation coefficients among the corrections δE k (v, L) due to common sources of uncertainty.As shown in [24], the assumed level of correlation can significantly influence the spinaveraged frequencies of rovibrational transitions determined from an incomplete set of measurements of their hyperfine components.Correlation coefficients are defined as follows: where u(X, Y ) is the covariance between quantities X and Y , and u(X) = u(X, X) is the uncertainty of X.
Here, we make assumptions similar to those presented in Appendix C of [33]: • Perfect positive correlations are assumed between δE 4 (v, L) and δE 4 (v , L ), as well as between δE 5 (v, L) and δE 5 (v , L ). • For the other hyperfine coefficients (k = 4, 5), correlation coefficients between δE k (v, L) and δE k (v , L ) are expected to be positive, but are difficult to estimate.
In our analysis, we thus vary these coefficients between 0 and 1 to investigate their influence on the spin-averaged frequencies, and additional uncertainties associated with these correlations.• Perfect positive correlations are assumed between δE 2 (v, L) and δE 3 (v, L), as well as between δE 6 (v, L) and δE 7 (v, L).This stems from the fact that uncalculated QED corrections in both coefficients are of the same nature.• The correlation coefficient between δE 4 (v, L) and δE 5 (v, L) is estimated [33] to be r(δE 4 (v, L), δE 5 (v, L)) = 0.4016.• No other correlations are assumed between δE k (v, L) and δE l (v, L).
• No correlations are assumed between δE k (v, L) and δE l (v , L ) with k = l and (v, L) = (v , L ), except for those that are implied by the above hypotheses, e.g.: Except for the unknown correlation between δE k (v, L) and δE k (v , L ) (k = 4, 5), we do not consider the effect of uncertainties in other correlation coefficients.Indeed, although these uncertainties are not strictly zero, they are estimated to be too small to have any significant impact on the spin-averaged frequencies.[27] (Ref.[33]).Those of E 2 /h, E 3 /h, E 8 /h, and E 9 /h are taken from [26].The value of E 9 /h has been updated using the latest determination of the deuteron's quadrupole moment [34].

Hyperfine splitting
The hyperfine structure of a rovibrational state (v, L) is obtained by diagonalizing the effective spin Hamiltonian (1).Given the relative sizes of the E k , the natural coupling scheme of angular momenta is [26] The matrix elements of the effective spin Hamiltonian in the basis of coupled states |vLF SJ are calculated using standard angular algebra procedures.J is an exact quantum number, whereas F and S are only approximate quantum numbers.A spin eigenstate is denoted by |vL F SJ , which is the state having the largest weight in its decomposition over the basis states.In the following, we omit the ˜symbols for simplicity.
Since the matrix of H eff is block diagonal with blocks of size 4 × 4 at most, its eigenvalues E hfs (v, L, F, S, J) can be expressed analytically as a function of the hyperfine interaction coefficients E k .These expressions are then linearised around the theoretical values of the coefficients using the sensitivity coefficients The theoretical hyperfine shifts E hfs and sensitivity coefficients γ k of all the levels involved in the measurements are given in Table 2.  2. Theoretical hyperfine shifts (column 3) and sensitivity coefficients (columns 4-12) for all the levels involved in high-precision measurements.Missing numbers for the (v, L) = (0, 0) state imply a zero value.The uncertainties of the hyperfine shifts are calculated assuming no correlation between δE k (v, L) and δE k (v , L ) for k = 4, 5, which provides an upper bound of the uncertainty.

Available experimental data
At the time of writing, ppt-range frequency measurements of three rotationalvibrational transitions in HD + had been reported in peer-reviewed publications.These are the (v, L): (0, 0) → (0, 1) rotational transition at 1.31 THz [14], the (v, L): (0, 0) → (1, 1) rovibrational transition at 58.6 THz [16], and the (v, L): (0, 3) → (9, 3) vibrational overtone at 415 THz [15].All frequency measurements were made on HD + ions stored in an rf trap, and sympathetically cooled to millikelvin temperatures via lasercooled Be + ions stored in the same trap.For further details on the experimental setups and the justification of the achieved measurement uncertainties we refer to the original publications.
Table 3 provides an overview of the frequencies (and their uncertainties) of the 10 measured hyperfine components belonging to the three rotational-vibrational transitions of interest.The measured frequencies are denoted by where a = (F, S, J) and a = (F , S , J ) are the initial and final hyperfine states.
In all cases the experimental uncertainties receive appreciable contributions from statistical as well as systematic sources of uncertainty.In the uncertainty evaluation of the (v, L): (0, 0) → (0, 1) and (v, L): (0, 0) → (1, 1) transitions, no mention is made of ] (where a = (F, S, J) and a = (F , S , J )), of various hyperfine components used in the adjustment of the spin-averaged transition frequencies (Section 4).The alternate labels for the hyperfine components in column 6 are those used in the original references (column 7).
possible correlated systematic errors [14,16].In the case of the (v, L): (0, 3) → (9, 3) transition, a correlation between two negligibly small (sub-0.01-kHz)uncertainties in the estimate of the Zeeman shift was identified [15].In what follows we assume that the uncertainties of all experimental transition frequencies are uncorrelated.

Adjustment details
The spin-averaged transition frequencies are determined from a least-squares adjustment of the available theoretical and experimental data presented in Secs. 2 and 3, respectively.We follow the fitting procedure described in Appendix E of [35].
All experimental data listed in Table 3 are included as input data for the adjustment.The associated "observational equations" (see [35] for details) read where the dotted equality sign means that the left and right hand sides are not equal in general (since the set of equations is overdetermined) but should agree within estimated uncertainties.Here, f SA [(v, L) → (v , L )] are the spin-averaged frequencies of the rovibrational transitions.Moreover, there are observational equations associated with the additive corrections δE th k (v, L) to theoretical hyperfine coefficients, as explained in Section 2.1: The adjusted constants are f SA [(v, L) → (v , L )] and δE th k (v, L).However, some of the δE k (v, L) are redundant because certain pairs of coefficients are perfectly correlated with each other (see Section 2.1).This leads to a singular covariance matrix, which cannot be inverted as required for the least-squares adjustment.To solve this, redundant variables are eliminated using the relationships: where u k (v, L) is the uncertainty of δE k (v, L).There are a total of N = 32 input data: 10 experimental frequencies (Table 3), and 22 additive energy corrections (Table 4).The total number of adjusted constants is M = 25: 3 spin-averaged transition frequencies and the 22 energy corrections.
Table 4. Input data for the additive energy corrections to account for the theoretical uncertainties of hyperfine interaction coefficients.

Effect of unknown hyperfine correlation coefficients
As discussed in Section 2.1, the level of correlation between some of the theoretical hyperfine coefficients (namely, δE k (v, L) and δE k (v , L ) with k = 4, 5) is not well known.We thus varied their correlation coefficients between 0 and 1 to study their influence on spin-averaged frequencies.These frequencies were found to be mostly sensitive to the correlation between δE 1 (v, L) and δE 1 (v , L ), with a smaller effect from the correlation between δE 6 (v, L) and δE 6 (v , L ) (which themselves are fully correlated with δE 7 (v, L) and δE 7 (v , L ), respectively).Correlations involving E 2 , E 3 , E 8 , and E 9 have a negligible impact.In order to estimate the additional shift, f corr (v, L, v , L ), and its uncertainty, δf corr (v, L, v , L ), caused by the unknown correlation coefficients, we carried out a Monte-Carlo simulation similar to [24], as further detailed in Appendix A. We find that the mean shift of the simulated spin-averaged transition frequencies is well accounted for when carrying out a single adjustment using correlation coefficients of 0.5.Furthermore, we determine the full widths of the simulated frequency distributions to be 2.09 Hz, 144 Hz, and 237 Hz for the rotational, v = 0 → 1, and v = 0 → 9 transitions, respectively.We conservatively define δf corr (v, L, v , L ) to be equal to the full width divided by two.As will be shown below, these correlation-induced uncertainties are considerably smaller than the frequency uncertainties resulting from the adjustment itself.The question is whether the uncertainties δf corr (v, L, v , L ) are correlated (which would introduce correlations between different transition frequencies f SA [(v, L) → (v , L )] in addition to the correlations that may result from the adjustment itself).We empirically find either strong correlations (i.e.r(δf corr (v, L, v , L ), δf corr (v , L , v , L )) = +1) or strong anticorrelations (i.e.r(δf corr (v, L, v , L ), δf corr (v , L , v , L )) = −1); see Appendix A. We can thus construct a 'correlation-induced' covariance matrix from the correlation-induced uncertainties defined above, in combination with the empirically found correlation coefficients.

Results and discussion
In Table 5, we compare our derived spin-averaged frequencies with those in the literature.These literature values are based on the same experimental data as used in our analysis; literature values based on specific subsets of the experimental data (such as the value presented in [25]) are not considered here.The rotational transition is shifted by about four standard deviations with respect to the literature value of [14].This is mainly due to the fact that we make different assumptions regarding correlations between theoretical hyperfine coefficients, especially those between E 4 (0, 0) and E 4 (0, 1).As justified in Ref. [33], we assume perfect positive correlation, whereas zero correlation was assumed in the initial analysis [14].The sensitivity of the rotational transition on the assumed level of correlation was already observed in [24] (see Fig 4 (b) in that reference).For the other transitions, the frequencies determined in this work are well within the error bars of the previous determinations.
As can be seen in Table 5, the improvement of hyperfine structure theory [33] has allowed the reduction of the uncertainties of the v = 0 → 1 and v = 0 → 9 transition frequencies with respect to their initial determinations [15,16].For the rotational transition, the precision has only slightly improved, because in this case 6 spin components have been measured, making the extracted frequency much less sensitive to theoretical uncertainties.
The normalised residuals of the adjustment (assuming a value of 0.5 for the unknown correlation coefficients) are shown in Figure 1.They reveal substantial discrepancies, with normalised residuals larger than two, between the experimental data and the theoretical hyperfine structure: • the hyperfine components A1, A3, and A5 of the rotational transition (with labels Ai defined in Table 3) have discrepancies of 3.0σ (121 Hz), 2.1σ (35 Hz), and −7.1σ (−427 Hz), respectively; • both components of the v = 0 → 9 transition, A9 and A10, exhibit discrepancies of −6.9σ (−4.5 kHz) and 6.1σ (3.6 kHz), respectively; • several theoretical spin coefficients, in particular E 1 (v = 0, L = 1) and E 6 (v = 0, L = 1), with deviations of 4.2σ (h×484 Hz) and 4.7σ (h×85 Hz), respectively.Deviations slightly above 2σ are also found in Discrepancies in the hyperfine structure of the L = 0 → 1 and v = 0 → 9 transitions had already been discussed in [27,33].Their origin is presently unknown; they could be due to a problem in the theory, in the experiments, or both.It is furthermore worth mentioning that experimental data on the hyperfine structure of H + 2 are in good agreement with theoretical predictions; see [33] and references therein.
We treat these deviations following the procedure used in the CODATA adjustments (see, e.g., [1]).A multiplicative expansion factor is applied to the initially assigned uncertainties of all (experimental and theoretical) input data, such that the absolute values of all normalised residuals are smaller than two.The required expansion factor is η = 3.56.As a result, the uncertainties of the spin-averaged frequencies are multiplied by the same factor (third line in Table 5).The expanded uncertainties are one order of magnitude larger than the frequency uncertainties caused by the unknown correlation coefficients.
The spin-averaged frequencies obtained by the least-squares adjustment have small but nonzero correlations between them due to correlations between δE k in different rovibrational levels.The magnitude of the corresponding correlation coefficients is less than 0.01, and their values (and in some case also their signs) depend on the choice of hyperfine correlation coefficients (Section 2.1).Therefore, the correlations between the three adjusted spin-averaged transition frequencies have an uncertainty due to the unknown hyperfine theory correlations.In Section 5 we show that these correlations have a negligible impact on the fundamental constants that may be derived from the HD + data, and we do not investigate this particular source of uncertainty further here.However, for completeness we do combine the covariance matrix resulting from the adjustment with the correlation-induced covariance matrix that was constructed as described in Section 4.2.The resulting (combined) uncertainties of the spin-averaged transition frequencies are listed in the last line of Table 5, and their (combined) correlation coefficients are presented in Table 8.Table 5. Spin-averaged transition frequencies f SA (in kHz).Previous determinations from the original publications [14][15][16] are given in the first line.The results of the least-squares adjustment performed in the present work are given in line 2.These results were obtained with an expansion factor η = 1 (i.e.no expansion factor applied).The spin-averaged frequencies in line 3 are obtained by applying an expansion factor η = 3.56 to the uncertainties of all (experimental and theoretical) input data.Finally, line 4 shows the frequencies for η = 3.56 and including the hyperfine-correlation-induced uncertainty; see text for details.The frequencies on line 4 are our recommended values.3 and 4.

Implications for the determination of particle masses
Using the spin-averaged frequencies determined in this work, we can now assess the potential contribution of HD + spectroscopy to the determination of fundamental constants, in particular the relative atomic masses A r (e), A r (p), A r (d), and related mass ratios.To do this, we use a least-squares adjustment following the same procedure as in the previous section, and compare three different scenarios: • Adjustment 1 only includes Penning trap measurements related to the electron, proton, and deuteron relative atomic masses; • Adjustment 2 only includes HD + spectroscopic data and values of fundamental constants (other than atomic masses) involved in the theoretical energy levels of HD + ; • Adjustment 3 includes all input data of Adjustments 1 and 2.
Theoretical rovibrational transition frequencies of HD + depend on [20]: the Rydberg constant (R ∞ ), nucleus-to-electron mass ratios λ p ≡ m p /m e and λ d ≡ m d /m e , the charge radii of the proton r p and deuteron r d , and the fine-structure constant α.The dependence on α is very weak because it only enters in the relativistic and QED corrections, and can be safely neglected without loss of precision.As noted in previous works, e.g., [14,17], the mass dependence can to a good approximation be reduced to a dependence on a single parameter λ pd ≡ µ pd /m e , where µ pd is the proton-deuteron reduced mass.This is exactly true in the adiabatic approximation, but, due to nonadiabatic effects, there is actually a much weaker dependence on a second mass ratio, for example λ d .In our Adjustment 2, we neglect this weak dependence, allowing λ pd to be determined from HD + spectroscopic data alone, as done in [14,16].Although λ pd is not an adjustable parameter in Adjustment 1, it is possible to infer a value of λ pd from the adjusted values of A r (e), A r (p), and A r (d).The relative impact of HD + spectroscopy and Penning-trap measurements on future values of particle mass ratios can then be assessed by comparing the determinations of λ pd (and their uncertainties) obtained from Adjustments 1 and 2. To further illustrate the effect of including HD + data, we also compare the values of λ p (= m p /m e ) inferred from Adjustments 1 and 3.
In the perspective of a global adjustment of fundamental constants, it is preferable to take the dependence of the theoretical HD + transition frequencies on both parameters λ pd , λ d into account in order to avoid any unnecessary loss of precision.Furthermore, for the global Adjustment 3, we have found it more convenient to parametrise the mass dependence of HD + transition frequencies in terms of the quantities that are treated as adjusted parameters in the CODATA adjustments [1], i.e., the relative atomic masses A r (e), A r (p), and A r (d).

Adjustment details
Let us briefly describe Adjustment 1.It includes the 2018 CODATA value of electron's relative atomic mass, for which no new measurement has been reported since then.For the proton relative atomic mass, we use the value of A r ( 1 H) from the latest (2020) Atomic Mass Evaluation (AME) [36,37], corrected for the electron mass and for the theoretical binding energy.This value takes into account recent high-precision Penning trap data [5,6], except for a very precise measurement of the deuteron-to-proton mass ratio, m d /m p , performed in 2021 [7].We thus also include the latter, together with the value of the deuteron mass deduced from the 2020 AME value of A r ( 2 H) (again corrected for the electron mass and the binding energy), improving the proton mass determination.Adjustment 1 thus comprises N = 4 input data and M = 3 adjusted parameters, the relative atomic masses of the electron, proton and deuteron.
The other adjustments include HD + spectroscopic data.Similarly to the hyperfine structure, theoretical uncertainties of spin-averaged transition frequencies are accounted for by introducing additive corrections δf th (v, L, v , L ) to the theoretical frequencies.They are treated as adjusted constants, and input data with a zero value and an uncertainty equal to the theoretical uncertainty of f th SA [(v, L) → (v , L )] are included (B1-B3 in Table 7).Their correlation coefficients (listed in Table 8) are estimated using the results of [19].Correlation coefficients are found to be close to 1, because the uncertainties are dominated by unevaluated QED terms, whereas numerical uncertainties are negligibly small.
In Adjustment 2, each of the three spin-averaged frequencies determined in the previous section (inputs A1-A3 in Table 7) is associated with an observational equation of the form Here, is the reference value of the theoretical transition frequency [19], calculated with 2018 CODATA values of fundamental constants.The reference (2018 CODATA) values of A r (e), R ∞ , r p , and r d are given in Table 7.Those of the proton and deuteron atomic masses are: A CODATA2018 As explained above, we furthermore have observational equations of the form with input value δf (v, L, v , L ) = 0, to deal with theoretical uncertainties.Finally, we require observational equations for R ∞ , r p , r d as well as relative atomic masses and mass ratios.They are written as etc., where the input data R CODATA2018 ∞ , r CODATA2018 p , . . .are found in Table 7 under labels C1-C7 (note that only C5-C7 are used in Adjustment 2).The correlation coefficients among these input data are given in Table 8.
The coefficients β, η in Eq. ( 8) depend on (v, L, v , L ), but for brevity this dependence is not indicated.β R∞ is the linear sensitivity coefficient to the Rydberg constant, and is simply given by In the dependence on the nuclear charge radii, following CODATA practice we only take into account the leading-order finite-size correction (see Eq. ( 6) of [38]).Higherorder finite-size corrections for the deuteron, which are included in the theoretical predictions [19], are taken as fixed, and the corresponding uncertainties are taken into account in δf .This yields a quadratic dependence on r p and r d , with sensitivity coefficients with i = p,d.a 0 is Bohr's radius, and x p (x d ) is the position vector of the electron with respect to the proton (deuteron).Brackets denote expectation values over the nonrelativistic wavefunction of the (v, L) rovibrational state, which are calculated numerically with high precision using variational wavefunctions [18].Finally, in order to evaluate the linear sensitivity coefficients with respect to the mass ratios, β λpd and β λd , it suffices to consider the nonrelativistic energy levels, which yields a relative accuracy on the order of α 2 for these coefficients.Following the approach of [39], they can be expressed as where X is the position of the proton with respect to the deuteron.The expectation values are calculated numerically with high precision using variational wavefunctions [39].For Adjustment 3, the mass dependence is parametrised in terms of the relative atomic masses, A r (e), A r (p), and A r (d).The observational equation (8) becomes Here, the observational equations, values, and correlation coefficients for δf are the same as for Adjustment 2. The sensitivity coefficients β Ar(i) (i = e, p, d) with respect to the relative atomic masses are given by The theoretical frequencies and all the sensitivity coefficients appearing in Eqs. ( 8) and ( 14) are given in Table 6.The dependence of the sensitivity coefficients on fundamental constants is weak, so that the values of Table 6 can be safely used during all iterations of the least-squares adjustment procedure.Note that the theoretical frequencies differ very slightly from those given in [19]   CODATA values of A r (e), A r (p) and A r (d), whereas the values of [19] were obtained with 2018 CODATA values of m p /m e and m d /m e .When switching between these two parametrizations, round-off errors in the published recommended values cause tiny shifts in the nonrelativistic energy levels.
In total, Adjustment 2 contains N = 9 input data and M = 7 adjusted parameters: λ pd , R ∞ , r p , r d , and the additive corrections for the three rovibrational transition frequencies.
Adjustment 3 combines the input data of the first two adjustments, with N = 13 input data and M = 9 adjusted parameters.

Results and discussion
The determinations of λ pd and of the proton-electron mass ratio, λ p , from the adjustments described above are shown in Table 9.
Adjustment 1 illustrates the potential precision permitted by Penning trap measurements alone.The relative uncertainty of λ pd is found to be 3.0 × 10 −11 , mainly limited by the electron's relative atomic mass that has a relative uncertainty of 2.9 × 10 −11 (see Table 7).The contribution from the proton and deuteron relative atomic masses to the uncertainty of λ pd is 7. Table 7. Input data for the adjustments discussed in Section 5. Entries C1 and C5-C7 also serve as reference values in equations ( 8) and (14).
Correlation coefficients r(A1, A2) = 0.00036 r(A1, A3) = 0.00644 r(A2, A3) = 0.00741 r(B1, B2) = 0.99570 r(B1, B3) = 0.95733 r(B2, B3) = 0.97998 r(C1, C5) = 0.00704 r(C1, C6) = −0.00133r(C1, C7) = 0.00317 r(C2, C3) = 0.31986 r(C5, C6) = 0.88592 r(C5, C7) = 0.90366 r(C6, C7) = 0.99165  Adjustment 2 provides an alternative determination of λ pd from HD + spectroscopic data.Its uncertainty is slightly smaller (2.1 × 10 −11 ), and the obtained values are in reasonable agreement, with a mild (1.6σ) tension.Since the precision of the first adjustment is limited by the electron's atomic mass, the smaller uncertainty from HD + data indicates that the impact of HD + on the global CODATA adjustment would primarily be to improve the precision of A r (e).This is verified in the last adjustment (Adjustment 3), which illustrates the potential improvement that can be obtained by combining Penning trap with HD + spectroscopy data.Whereas the proton and deuteron atomic masses are very precisely determined by mass spectrometry, HD + spectroscopy, being sensitive to the nucleus-to-electron mass ratios, provides a link to the electron mass, allowing to improve the precision of A r (e) and m p /m e to 1.8 × 10 −11 , which for the latter would represent an improvement by more than a factor of three over its 2018 CODATA uncertainty.It is worth noting that the adjustment reveals no significant tension in the data, with all normalised residuals being smaller than 1.2 in absolute value.The found values of A r (e) and m p /m e are furthermore in good agreement with the 2018 CODATA values.
We carried out Adjustment 3 taking into account the weak correlations between the adjusted experimental spin-averaged frequencies, discussed in Section 4.3.To assess the influence of these correlations, we repeated Adjustment 3 assuming zero correlation between the spin-averaged frequencies.The value of m p /m e obtained from this adjustment is identical to the one shown in the last lines of Table 9, indicating a negligible impact of these correlations on the determination of fundamental constants from recent HD + theory and measurements.

Conclusion
In conclusion, we have carried out a comprehensive analysis of the available highprecision spectroscopic data of rotational and vibrational transitions in HD + , with the aim of providing a set of spin-averaged rotational-vibrational transition frequencies that can serve as input data for future CODATA adjustments of the fundamental physical constants.A crucial part of the analysis involves the use of theoretically computed hyperfine structure to remove hyperfine shifts from the measured transition frequencies, taking into account the correlations between theroretical spin coefficients.The experimental input data are complemented by the most recent theoretical pre-dictions of the spin-averaged transition frequencies, and by linearised expressions that parametrise these theoretical predictions in terms of relevant fundamental constants.Uncertainties and correlation coefficients pertaining to the input data were also provided, which completes the set of observational equations and covariances needed for future CODATA adjustments.To illustrate the potential of the HD + data, we have carried out adjustments of the electron, proton and deuteron relative atomic mass using state-of-the-art data from Penning-trap measurements as well as the HD + data presented here.We found that the HD + data have a particularly large impact on the precision of the electron relative atomic mass which, together with the precise values of the proton and deuteron relative atomic masses from Penning-trap measurements, may improve the precision of m p /m e by more than a factor of three compared to the current 2018 CODATA value.frequency obtained for the (0.5, 0.5, 0.5, 0.5, 0.5) case is 0.05 Hz, −3 Hz, and −2 Hz (i.e. less than 3% of the full width).Figure A1 furthermore shows the dependence of the correlation-induced shift if all correlation coefficients are varied as (ξ, ξ, ξ, ξ, ξ), with ξ between 0 and 1. Comparing the shifts of the three transitions shows that the shift of the rotational transition is anticorrelated with the shift of the v = 0 → 1 and v = 0 → 9 transitions, whereas the shifts for the latter two are correlated.

Figure 1 .
Figure 1.Normalised residuals of the 32 input data for the adjustment of spin-averaged transition frequencies.Labels Ai and Bi follow those defined in Tables3 and 4.
= 2.013 553 212 745.The mass ratios λ pd and λ d are calculated using λ pd = A r (p)A r (d)/(A r (e)[A r (p) + A r (d)]) and λ d = A r (d)/A r (e) (and similar for the reference values λ pd0 and λ d0 ).

Table 1 .
Theoretical hyperfine coefficients of rovibrational states of HD + involved in high-precision measurements, in kHz.Missing numbers in the first line imply a zero value.The values of E 4 /h, E 5 /h (E 1 /h, E 6 /h, E 7 /h) were calculated in Ref.

Table 3 .
Experimentally determined transition frequencies, f exp because they are calculated with 2018

Table 6 .
Reference theoretical values and sensitivity coefficients of HD + transition frequencies.The theoretical frequencies are calculated using 2018 CODATA values of Ar(e), Ar(p), Ar(d), cR∞, rp, r d , and α.The sensitivity coefficient for the Rydberg constant β cR∞ is dimensionless.u stands for the relative atomic mass unit.

Table 8 .
Values of all nonzero correlation coefficients in the input data of Table7.

Table 9 .
Adjustment results.ur stands for relative uncertainty.