Accurate Potential Energy Curves Research Articles

Electronic structure, oscillator strength, and rovibrationally resolved photodissociation of 24MgH

Aims. A series of high-precision calculations for the electronic structure of MgH have been reported in the past two decades; however, most of them d not include the core-valence correlation and still exhibit distinct differences. Furthermore, the latest high-precision results have not been applied to the studies of photodissociation dynamics. The primary motivations of this paper are to calculate a more precise electronic structure of MgH consering a core-valence correlation and to prove the photodissociation cross-sections. Methods. The electronic structure of MgH is investigated by multi-reference configuration interaction calculations with Davson correction (MRCI+Q). We performed two different sets of calculations to investigate the core-valence correlation and, as a result, obtained accurate potential energy curves (PECs) and transition dipole moments (TDMs). An extrapolation procedure was also employed to eliminate the error of basis set. Then, the photodissociation cross-sections were calculated using high-precision PECs and TDMs. Results. The PECs and TDMs of the five lowest doublet electronic states, X2Σ+, B′2Σ+, E2Σ+, A2Π, and C2Π, are obtained from calculations including core-valence correlation, termed as CV-MRCI, while PECs of the ten lowest doublet states and three quartet states are also obtained from NCV-MRCI calculations without core-valence correlation. The spectroscopic constants and band oscillator strengths are also proved with high precision levels. The equilibrium Re and vertical excitation energy Te are only 0.1% different from the measurements. Based on the CV-MRCI results, the rovibrationally resolved photodissociation cross-sections for transitions from X2Σ+ to the other four states, as well as the total local thermodynamic equilibrium cross-sections for temperatures up to 10000 K, are calculated.

Quantum calculation of the pressure broadening of H and K lines of Ca+ calcium ion in a gaz of ground-state He helium atoms

In this theoretical study, we have conducted calculations to investigate the pressure broadening of the H ( 4 s 2 S 1 / 2 − 4 p 2 P 1 / 2 ) and the K ( 4 s 2 S 1 / 2 − 4 p 2 P 3 / 2 ) lines of the Ca + calcium ion in a gas of ground-state He helium atoms. The focus of our work is to determine the profiles of the pressure-broadened resonance lines and analyse their behaviour under the influence of binary collisions between Ca + ions and He atoms. To achieve this, we constructed accurate potential energy curves for the ground X 2 Σ 1 / 2 + and the excited D 2 Π 1 / 2 , D 2 Π 3 / 2 and E 2 Σ 1 / 2 + states. Additionally, we constructed the corresponding transition dipole moments based on the data obtained through high-level ab initio methods, including SA-CASSCF, MRCI, and SO coupling, with Davidson and BSSE corrections. The study spans a range of temperatures from 4000 K to 12,000 K and wavelengths from 200 nm to 500 nm . The theoretical profile is characterised by dominant free-free transitions. In the vicinity of the 436 nm wavelength position, there is a satellite peak in the red wing, attributed to the contribution of D 2 Π 1 / 2 ⟵ X 2 Σ 1 / 2 + and D 2 Π 3 / 2 ⟵ X 2 Σ 1 / 2 + transitions. Our findings are in good agreement with previous theoretical results.

Calculating Potential Energy Surfaces with Quantum Computers by Measuring Only the Density Along Adiabatic Transitions.

We show that chemically accurate potential energy surfaces (PESs) can be generated from quantum computers by measuring only the density along an adiabatic transition between different molecular geometries. In lieu of using phase estimation, the energy is evaluated by performing line-integration using the inverted real-space time-dependent density functional theory Kohn-Sham (KS) potential obtained from the geometry-varying densities of the full wave function. The accuracy of this method depends on the validity of the adiabatic evolution itself and the potential inversion process (which is theoretically exact but can be numerically unstable), whereas the total evolution time is the defining factor for the precision of phase estimation. We examine the method with a one-dimensional system of two electrons for both the ground and first triplet states in first quantization, as well as the ground state of three- and four-electron systems in second quantization. It is shown that few accurate measurements can be utilized to obtain chemical accuracy across the full potential energy curve, with a shorter propagation time than may be required using phase estimation for a similar accuracy. We also show that an accurate potential energy curve can be calculated by making many imprecise density measurements (using a few shots) along the time evolution and smoothing the resulting density evolution. Finally, it is important to note that the method is able to classically provide a check of its own accuracy by comparing the density resulting from a time-independent KS calculation using the inverted potential with the measured density. This can be used to determine whether longer adiabatic evolution times are required to satisfy the adiabatic theorem.

Energy-dependent photoion angular distributions in two-body Coulomb explosions of molecules.

We experimentally study two-body Coulomb explosions of CO2, O2, and CH3Cl molecules in intense femtosecond laser pulses. We observe an obvious variation in the ionic angular distribution of the fragments with respect to the kinetic energy releases (KERs). Using a classical model based on abinitio potential energy curves, we find that the dependence of the ionic angular distribution on the KER is relevant to the fact that the accurate potential energy deviates significantly from the value determined by applying the Coulomb interaction approximation at a relatively small internuclear distance of the molecule. We show that the KER-dependent ionic angular distribution provides an effective way to determine the critical internuclear distance at which the Coulomb interaction approximation holds or breaks down without relying on the knowledge of the accurate potential energy curves.

Transport cross sections and collision integrals for N(4S)–O(3P) and N(4S)–O(1D) interactions in high-temperature air plasmas

Collisions between nitrogen (N) and oxygen (O) play a crucial role in determining transport coefficients in high-temperature atmospheres of Earth and planetary. In this study, the momentum transfer, viscosity, third-moment, and fourth-moment transport cross sections for the N(4S)–O(3P) and N(4S)–O(1D) interactions are reported in the collision energy range of 10−6–10 Hartree based on the classical and semiclassical methods. The new and accurate potential energy curves for N–O interactions, which are used to provide the input for calculations of the cross sections, are calculated based on the state-of-the-art ab initio method. The classical and semiclassical collision integrals are provided at 300–50 000 K, and the results support the calculation of transport coefficients in a third-order approximation. In particular, the collision data for the N(4S)–O(1D) interaction based on ab initio points are reported for the first time. The calculated transport cross sections and collision integrals are helpful for studies of modeling the high-temperature air plasmas.

Benchmark Study of the Electronic States of the LiRb Molecule: Ab Initio Calculations with the Fock Space Coupled Cluster Approach.

Accurate potential energy curves (PECs) are determined for the twenty-two electronic states of LiRb. In contrast to previous studies, the applied approach relies on the first principle calculations involving correlation among all electrons. The current methodology is founded on the multireference coupled cluster (CC) scheme constructed within the Fock space (FS) formalism, specifically for the (2,0) sector. The FS methodology is established within the framework of the intermediate Hamiltonian formalism and offers an intruder-free, efficient computational scheme. This method has a distinctive feature that, when applied to the doubly ionized system, provides the characteristics of the neutral case. This proves especially beneficial when investigating PECs in situations where a closed-shell molecule dissociates into open-shell fragments, yet its double positive ion forms closed-shell species. In every instance, we successfully computed continuous PECs spanning the entire range of interatomic distances, from the equilibrium to the dissociation limit. Moreover, the spectroscopic characteristic of various electronic states is presented, including relativistic effects. Relativistic corrections included at the third-order Douglas-Kroll level have a non-negligible effect on the accuracy of the determined spectroscopic constants.

Prediction of vibrational spectrum and thermodynamic properties for phosphorus mononitride

In this work, an accurate potential energy curve (PEC) for the ground electronic state of phosphorus mononitride (PN) molecule has been determined from a variationally improved Hulburt-Hirschfelder (VIHH) oscillator model in conjunction with the experimental spectral constants (De,ωe,ωexe,Be,αe,re). We have numerically solved the Schrödinger equation for the VIHH potential using the LEVEL program, obtaining the pure vibrational spectrum that converges to the dissociation limit. In addition, the partition functions of PN molecule are calculated using the full rovibrational energies. Ultimately, thermodynamic properties like molar heat capacity, entropy, enthalpy, and Gibbs free energy were calculated for the PN molecule and show good agreement with those data from the NIST (National Institute of Standards and Technology) database.

Ab initio calculation of the potential energy curves and spectroscopic properties of low-lying electronic states of CH+ molecular cation

The combined use of the MRCI-SD(T) method with the cc-pV6Z basis set has enabled us to determine the potential energy curves of the molecular cation whose fundamental state lies above the corresponding state of its neutral parent at the dissociation limit. With the accurate potential energy curves thus obtained, we have calculated spectroscopic constants which show good agreement with the experimental and other theoretical values available in the literature. To account for perturbations where necessary, spin-orbit coupling effects have been introduced. An understanding of the polarity of studied in the various electronic states is made possible by calculating the dipole moments. In addition, transition dipole moments were computed prior to the obtention of molecular parameters such as oscillator strength, Einstein coefficients, radiative probabilities and Franck–Condon factors. Overall good accuracy is observed.

Surface characterization using Friction Force Microscopy and the Jarzynski equality

The Jarzynski equality is a fundamental result from non-equilibrium statistical mechanics that provides the difference in free energy between two states from the work done on the system by processes that can be carried out far from equilibrium. We evaluate the applicability of Jarzynski equality to map the potential energy of a model graphene surface using data from simulated Friction Force Microscopy (FFM). We model the scanning process of the FFM using the Prandtl-Tomlinson model and Langevin dynamics. By varying the simulation parameters, we verify the “stick–slip” and thermolubricity friction regimes, as well as the crossover between them. We then calculate the surface potential energy using the Jarzynski equality for these regimes. A new method for properly evaluating the free energy of the cantilever is introduced. We observe that the applicability of Jarzynski's equality is linked to the friction regimes: For thermolubricity, a very accurate potential energy curve is obtained for relatively few repetitions, but for the “stick–slip” movement, it is only possible to use Jarzynski’s equation in a small fraction of the scanning distance. For the crossover regime, it is possible to obtain a relatively accurate potential energy curve for a sufficiently large number of sampling repetitions.

High-accuracy Rb2+ interaction potentials based on coupled-cluster calculations

This work discusses a protocol for constructing highly accurate potential energy curves (PECs) for the lowest two states of Rb$_{2}^+$, i.e. $X\,{}^2{\Sigma}{_g^+}$ and $(1) {}^2\Sigma{_u^+}$, using an additivity scheme based on coupled-cluster theory. The approach exploits the findings of our previous work [J. Schnabel, L. Cheng and A. K\"ohn, J. Chem. Phys. 155, 124101 (2021)] to avoid the unphysical repulsive long-range barrier occurring for symmetric molecular ions when perturbative estimates of higher-order cluster operators are employed. Furthermore, care was taken to reproduce the physically correct exchange splitting of the $X {}^2{\Sigma}{_g^+}$ and $(1) {}^2{\Sigma}{_u^+}$ PECs. The accuracy of our computational approach is benchmarked for ionization energies of Rb and for spectroscopic constants as well as vibrational levels of the $a {}^3{\Sigma}{_u^+}$ triplet state of Rb\textsubscript{2}. We study high-level correlation contributions, high-level relativistic effects and inner-shell correlation contributions and find very good agreement with experimental reference values for the atomic ionization potential and the binding energy of Rb$_{2}$ in the $a\,{}^3{\Sigma}{_u^+}$ triplet state. Our final best estimate for the binding energy of the Rb$_{2}^+$ $X {}^2{\Sigma}{_g^+}$ state including zero-point vibrational contributions is $D_0 = 6179\,\mathrm{cm}^{-1}$ with an estimated error bound of $\mathcal{O}(\pm 30\,\mathrm{cm}^{-1})$. This value is smaller than the experimentally inferred lower bond of $D_0\ge 6307.5\,\mathrm{cm}^{-1}$ [Bellos et al., Phys. Rev. A 87, 012508 (2013)] and will require further investigation. For the $(1) {}^2{\Sigma}{_u^+}$ state a shallow potential with $D_0 = 78.4\,\mathrm{cm}^{-1}$ and an error bound of $\pm 9\,\mathrm{cm}^{-1}$ is computed.

Intensities of KCs [formula omitted] band system up to dissociation threshold: an interplay between spin-orbit, hyperfine and rovibronic coupling effects

The relative intensity distribution in the rotationally resolved laser-induced fluorescence spectra belonging to the E(4)1Σ+→(a3Σ1+,X1Σ+) band systems of the KCs molecule was analyzed. The experimental intensities in doublet P,R progressions assigned to spin-allowed E→X and spin-forbidden E→a transitions up to their common ground dissociation limit were described in the framework of a coupled-channels (CC) deperturbation model applied for the interacting X1Σ+ and a3Σ+ states. The CC intensity simulation was based solely on fixed electronic structure parameters as functions of the internuclear distance R, namely: accurate empirical potential energy curves for all three states, ab initio estimates for matrix elements A(R) of the hyperfine structure (HFS), and transition dipole moments dEX(R) and dEa(R). A comparison between the measured intensities and their theoretical counterparts demonstrates a strong competition between different intramolecular interactions. A weak spin-orbit coupling of the upper E(4)1Σ+ state with the remote 3Π states is responsible for appearance of the E→a vibrational bands for the intermediate va-values. In turn, the HFS coupling between X1Σ+ and a3Σ+ states leads to peculiarities in E→(X,a) intensities, which are pronounced for high vX/a-values in the vicinity of K(42S)+Cs(62S) dissociation threshold. Both adiabatic ro-vibrational and non-adiabatic electronic-rotational interactions explain the abrupt deviation of some observed P/R intensity ratios from the expected Hönl-London factors.

Potential energy curves for electronic states of the sodium dimer with multireference coupled cluster calculations

Accurate potential energy curves (PECs) were obtained for the 26 lowest lying electronic states of the sodium dimer. The current approach is based on the first principle size-extensive multireference coupled cluster method formulated in the (2,0) sector of the Fock space (FS). The (2,0) sector of FS provides the description of the states obtained by the attachment of two electrons to the reference system. This makes it possible to adopt the doubly ionised sodium dimer as a reference with its very concrete advantage in the calculations of the PECs, i.e. it dissociates into closed shell fragments (i.e. sodium cations) hence the restricted Hartree–Fock reference can be used in the whole range of interatomic distances. This results in the correct description of the interatomic potentials for all bond lengths. The theoretical PECs stay very close to the experimental curves wherever the latter are available and the computed spectroscopic parameters reproduce the experiment with very good accuracy. Relativistic corrections included at the second-order Douglass–Kroll level have a non-negligible effect on the accuracy of computed excitation and dissociation energies.

The vibrational and hyperfine spectroscopy toward laser cooling 87Sr35Cl

We theoretically investigate the possibility of laser cooling 87Sr35Cl molecule in accordance with vibrational and hyperfine spectroscopy. The potential energy curves and dipole moment of the X2Σ+1/2, A2Π1/2,3/2 and B2Σ+ states are calculated using ab initio method and the spectroscopic parameters are in good agreement with the experimental data. On account of the accurate potential energy curves and the transition dipole moment, the Franck – Condon factors and radiative lifetimes are predicted. Comparing the conditions of laser cooling candidate molecules, the A2Π 1/2 ↔ X2Σ+1/2 transition is selected as the laser cooling cycle system. In order to obtain an approximately closed cooling cycle system, we employed matrix element algorithm to calculated the hyperfine spectroscopy and branching ratios of the 87Sr35Cl molecule. Furthermore, an electro-optical modulator (EOM) is designed including six hyperfine levels of the ground state X2Σ+1/2 (v = 0, N = 1).

Erratum: "Accurate semiempirical potential energy curves for the a3Σ+-state of NaCs, KCs, and RbCs" [J. Chem. Phys. 154, 154304 (2021)

Erratum: "Accurate semiempirical potential energy curves for the a3Σ+-state of NaCs, KCs, and RbCs" [J. Chem. Phys. 154, 154304 (2021)

Theoretical study on spectroscopic properties of 8 Λ-S and 23 Ω states for BH molecule

In this work, the potential energy curves of eight low electronic states (X<sup>1</sup>Σ<sup>+</sup>, a<sup>3</sup>Π, A<sup>1</sup>Π, b<sup>3</sup>Σ<sup>-</sup>, 2<sup>3</sup>Π, 1<sup>3</sup>Σ<sup>+</sup>, 1<sup>5</sup>Σ<sup>-</sup>, and 1<sup>5</sup>Π) and twenty-three Ω states of BH molecule, and the transition dipole moments among the <inline-formula><tex-math id="M10">\begin{document}$ {\text{X}}{}^{\text{1}}{\Sigma}_{{{\text{0}}^ + }}^ + $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M10.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M10.png"/></alternatives></inline-formula>, <inline-formula><tex-math id="M11">\begin{document}$ {{\text{a}}^{\text{3}}}{\Pi_{{{\text{0}}^ + }}} $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M11.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M11.png"/></alternatives></inline-formula>, a<sup>3</sup>Π<sub>1</sub>, a<sup>3</sup>Π<sub>2</sub>, and A<sup>1</sup>Π<sub>1</sub> states are calculated by using the internally contracted multireference configuration interaction (icMRCI) method. In order to obtain the accurate potential energy curve, the errors caused by single and double electron excitation, core-valence correlation effects, relativistic effects and basis set truncation are corrected. The spectral and transition data of BH molecule are in good agreement with the available theoretical and experimental data. The calculation results show that the A<sup>1</sup>Π<sub>1</sub>(<i>υ′</i> = 0-2, <i>J′</i> = 1, +) →<inline-formula><tex-math id="M12">\begin{document}$ {\text{X}}{}^{\text{1}}{\Sigma}_{{{\text{0}}^ + }}^ + $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M12.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M12.png"/></alternatives></inline-formula>(<i>υ′′</i> = 0-2, <i>J′′</i> = 1, –) transition has large Einstein <i>A</i>-coefficient, weighted absorption oscillator strength, and highly diagonal vibrational branching ratio<i> R</i><sub><i>υ′υ′′</i></sub>, and the excited state A<sup>1</sup>Π<sub>1</sub>(<i>υ′</i> = 0, 1) have short spontaneous radiation lifetimes. Moreover, the effects of <inline-formula><tex-math id="M13">\begin{document}$ {{\text{a}}^{\text{3}}}{\Pi_{{{\text{0}}^ + }}} $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M13.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M13.png"/></alternatives></inline-formula>and a<sup>3</sup>Π<sub>1</sub> states on A<sup>1</sup>Π<sub>1</sub>(<i>υ′</i> = 0) ↔ <inline-formula><tex-math id="M14">\begin{document}$ {\text{X}}{}^{\text{1}}{\Sigma}_{{{\text{0}}^ + }}^ + $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M14.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M14.png"/></alternatives></inline-formula>(<i>υ′′</i> = 0) cycle transition can be ignored. Therefore, according to the A<sup>1</sup>Π<sub>1</sub>(<i>υ′</i><sub> </sub>= 0-1, <i>J′</i> = 1, +) ↔ <inline-formula><tex-math id="M15">\begin{document}$ {\text{X}}{}^{\text{1}}{\Sigma}_{{{\text{0}}^ + }}^ + $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M15.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="10-20220038_M15.png"/></alternatives></inline-formula>(<i>υ′′</i><sub> </sub>= 0-3, <i>J</i>′′ = 1, –) cycle transition, we propose to apply one main cooling laser (<i>λ</i><sub>00</sub> = 432.45 nm) and two repumping lasers (<i>λ</i><sub>10</sub> = 479.67 nm and <i>λ</i><sub>21</sub> = 481.40 nm) to laser cooling BH molecules, and evaluation of the cooling effect.

Ab initio calculation of the ground and first excited states of the lithium dimer

Based on a high level ab initio calculation which is carried out with the multireference configuration interaction method under the aug-cc-pVXZ (AVXZ) basis sets, X = T, Q, 5, the accurate potential energy curves (PECs) of the ground state and the first excited state of Li2 are constructed. By fitting the ab initio potential energy points with the Murrell–Sorbie potential function, the analytic potential energy functions (APEFs) are obtained. The molecular bond length at the equilibrium (R e ), the potential well depth (D e ), and the spectroscopic constants (B e , ω e , α e , and ω e χ e ) for the state and the state are deduced from the APEFs. The vibrational energy levels of the two electronic states are obtained by solving the time-independent Schrödinger equation with the Fourier grid Hamiltonian method. All the spectroscopic constants and the vibrational levels agree well with the experimental results. The Franck–Condon factors (FCFs) corresponding to the transitions from the vibrational level (v′ = 0) of the ground state to the vibrational levels (v″ = 0–74) of the first excited state have been calculated. The FCF for the vibronic transition of (v″ = 0) ←(v′ = 0) is the strongest. These PECs and corresponding spectroscopic constants provide reliable theoretical references to both the spectroscopic and the molecular dynamic studies of the Li2 dimer.

Limitations of perturbative coupled-cluster approximations for highly accurate investigations of Rb2+

We reveal limitations of several standard coupled-cluster (CC) methods with perturbation-theory based noniterative or approximate iterative treatments of triple excitations when applied to the determination of highly accurate potential energy curves (PECs) of ionic dimers, such as the XΣg+2 electronic ground state of Rb2 +. Such computations are of current interest for the understanding of ion-atom interactions in the ultracold regime. We demonstrate that these CC methods lead to an unphysical long-range barrier for the Rb2 + system. The barrier is small but spoils the long-range behavior of the PEC. The effect is also found for other X2 + systems, such as X = Li, Na, and K. Calculations using a flexible framework for obtaining leading perturbative triples corrections derived using an analytic CC singles and doubles energy derivative formulation demonstrate that the origin of this problem lies in the use of T̂3 amplitudes obtained from approximate CC singles, doubles, and triples amplitude equations. It is shown that the unphysical barrier is related to a symmetry instability of the underlying Hartree-Fock mean-field solution, leading to orbitals representing two +0.5-fold charged ions in the limit of separated fragments. This, in turn, leads to a wrong 1/R asymptote of the interaction potential computed by perturbation-based CC approximations. Physically meaningful perturbative corrections in the long-range tail of the PEC may instead be obtained using symmetry-broken reference determinants.

Highly Accurate Potential Energy Curves for the Hydrogen Molecular Ion

AbstractPotential energy curves of the hydrogen molecular ion in the Born‐Oppenheimer approximation are computed by means of the Riccati‐Padé method (RPM). The convergence properties of the method are analysed for different states. The equilibrium internuclear distance, as well as the corresponding electronic plus nuclear energy, and the associated separation constants, are computed to 40 digits of accuracy for several bound states. For the ground state the same parameters are computed with more than 100 digits of accuracy. Additional benchmark values of the electronic energy at different internuclear distances are given for several additional states. The software implementation of the RPM is given under a free software license. The results obtained in the present work are the most accurate available so far, and further additional benchmarks are made available through the software provided.

Radial molecular property functions of CH in its ground electronic state

Effective potential energy curves of the ground electronic states of isotopomers of the methylidyne radical CH are constructed by morphing accurate ab initio potential energy curves within the framework of Jenč’s reduced potential energy curve (RPC) approach. The morphing is performed by fitting the RPC parameters to available experimental ro-vibrational data. In turn, the resulting potential energy curves are used to determine the radial property functions describing the spin-orbit and fine and hyperfine energy patterns of the spin-ro-vibrational energy levels by morphing their ab initio approximants or by a direct fitting to experiment. The resulting functions enable reliable predicting outside prospected energy regions, and a quantitative probing of the impacts of external perturbations on the observed spectra.

The stability of e+H- 2.

The recently discovered positronic molecule e+H- 2 [J. Charry et al., Angew. Chem., Int. Ed. 57, 8859-8864 (2018)] has a new type of bond, the single-positron bond. We studied its stability using quantum Monte Carlo techniques. We computed an accurate potential energy curve of the reaction H- + PsH → e+H- 2 → H2 + Ps- to establish its global stability with respect to all possible dissociation channels and to define the range of its local stability. We showed that the e+H- 2 system is stable with respect to the dissociation into H- + PsH, with a binding energy of 23.5(1) mhartree. For R < 3.2 bohrs, the system is unstable, and it decays into H2 + Ps-. There are no other bound structures for R < 3.2 bohrs. We discuss possible routes to its experimental production.