Potential energy and dipole moment surfaces of H(3) (-) molecule

Superconductivity in Dilute Hydrides of Ammonia under Pressure

The past decade has witnessed great progress in predicting and synthesizing polyhydrides that exhibit superconductivity under pressure. Dopants allow these compounds to become metals at pressures lower than those required to metallize elemental hydrogen. Here, we show that by combining the fundamental planetary building blocks of molecular hydrogen and ammonia, conventional superconducting compounds can be formed at high pressure. Through extensive theoretical calculations, we predict metallic metastable structures with NHn (n = 10, 11, 24) stoichiometries that are based on NH4+ superalkali cations and complex hydrogenic lattices. The hydrogen atoms in the molecular cation contribute to the superconducting mechanism, and the estimated superconducting critical temperatures, Tc's, are comparable to the highest values computed for the alkali metal polyhydrides. The largest calculated (isotropic Eliashberg) Tc is ∼180 K for Pnma-NH10 at 300 GPa. Our results suggest that other molecular cations can be mixed with hydrogen under pressure, yielding superconducting compounds.

Polyatomic radiative association by quasiclassical trajectory calculations: Formation of HCN and HNC molecules in H + CN collisions

We have developed the polyatomic extension of the established [M. Gustafsson, J. Chem. Phys. 138, 074308 (2013)] classical theory of radiative association in the absence of electronic transitions. The cross section and the emission spectrum of the process is calculated by a quasiclassical trajectory method combined with the classical Larmor formula which can provide the radiated power in collisions. We have also proposed a Monte Carlo scheme for efficient computation of ro-vibrationally quantum state resolved cross sections for radiative association. Besides the method development, the global potential energy and dipole surfaces for H + CN collisions have been calculated and fitted to test our polyatomic semiclassical method.

Analysis of QED and non-adiabaticity effects on the rovibrational spectrum of H3 + using geometry-dependent effective nuclear masses

The influence of QED effects (including one- and two-electron Lamb-shift, Araki-Sucher term, one-loop self-energy, and finite nuclear size correction) together with non-adiabatic effects on the rovibrational bound states of H₃ ⁺ has been investigated. Non-adiabaticity is modeled by using geometry-dependent effective nuclear masses together with only one single potential energy surface. In conclusion, for rovibrational states below 20 000 cm^-1, QED and relativistic effects do nearly compensate, and a potential energy surface based on Born-Oppenheimer energies and diagonal adiabatic corrections has nearly the same quality as the one including relativity with QED; the deviations between the two approaches for individual rovibrational states are mostly below 0.02 cm^-1. The inclusion of non-adiabatic effects is important, and it reduces deviations from experiments mostly below 0.1 cm^-1.

Magic numbers, quantum delocalization, and orientational disordering in anionic hydrogen and deuterium clusters

The Diffusion Monte Carlo (DMC) method was applied to anionic hydrogen clusters H^-(H₂)_n (n = 1-16, 32) and their deuterated analogs using a polarizable all-atom potential energy surface (PES) developed by Calvo and Yurtsever. For the hydrogen clusters, the binding energy ΔE_n appears to be a smooth function of the cluster size n, thus contradicting the previous claim that n = 12 is a "magic number" cluster. The structures of the low energy minima of the PES for these clusters belong to the icosahedral motif with the H₂ molecules aligned toward the central H^- ion. However, their ground state wavefunctions are highly delocalized and resemble neither the structures of the global nor local minima. Moreover, the strong nuclear quantum effects result in a nearly complete orientational disordering of the H₂ molecules. For the deuterium clusters, the ground state wavefunctions are localized and the D₂ molecules are aligned toward the central D^- ion. However, their structures are still characterized as disordered and, as such, do not display size sensitivity. In addition, DMC simulations were performed on the mixed H^-(H₂)_n(D₂)_p clusters with (n, p) = (6, 6) and (16, 16). Again, in contradiction to the previous claim, we found that the "more quantum" H₂ molecules prefer to reside farther from the central H^- ion than the D₂ molecules.

Quantum state-to-state study for (H−(D−),HD) collisions on two potential energy surfaces

Revealing the reaction mechanisms of the H⁻/D⁻ + HD reaction – an exact quantum dynamics study on two potential energy surfaces.

Low frequency vibrational anharmonicity and nuclear spin effects of Cl-(H2) and Cl-(D2)

Low frequency combination bands of ³⁵Cl^-(H₂) and ³⁵Cl^-(D₂) have been measured in the region between 600 and 1100 cm^-1 by infrared predissociation spectroscopy in a cryogenic 22-pole ion trap using a free electron laser at the FELIX Laboratory as a tunable light source. The ³⁵Cl^-(H₂) (³⁵Cl^-(D₂)) spectrum contains three bands at 773 cm^-1 (620 cm^-1), 889 cm^-1 (692 cm^-1), and 978 cm^-1 (750 cm^-1) with decreasing intensity toward higher photon energies. Comparison of the experimentally determined transition frequencies with anharmonic vibrational self-consistent field and vibrational configuration interaction calculations suggests the assignment of the combination bands v₁ + v₂, 2v₁ + v₂, and 3v₁ + v₂ for ³⁵Cl^-(H₂) and 2v₁ + v₂, 3v₁ + v₂, and 4v₁ + v₂ for ³⁵Cl^-(D₂), where v₁ is the ³⁵Cl^-⋯H₂ stretching fundamental and v₂ is the Cl^-(H₂) bend. The observed asymmetric temperature dependent line shape of the v₁ + v₂ transition can be modeled by a series of ∑⁺-∏ ro-vibrational transitions, when substantially decreasing the rotational constant in the vibrationally excited state by 35%. The spectrum of ³⁵Cl^-(D₂) shows a splitting of 7 cm^-1 for the strongest band which can be attributed to the tunneling of the ortho/para states of D₂.

Investigation of Nonadiabatic Effects for the Vibrational Spectrum of a Triatomic Molecule: The Use of a Single Potential Energy Surface with Distance-Dependent Masses for H3. J Phys Chem A 2017;121:7016-7030. [PMID: 28820589 DOI: 10.1021/acs.jpca.7b04703]

On the basis of first-principles, the influence of nonadiabatic effects on the vibrational bound states of H₃⁺ has been investigated using distance-dependent reduced masses and only one single potential energy surface. For these new vibrational calculations, potentials based on explicitly correlated wave functions are used where, in addition, adiabatic corrections and relativistic contributions are taken into account. For the first time, several different fully distance-dependent reduced mass surfaces in three dimensions have been incorporated in the vibrational calculations.

Anionic Hydrogen Cluster Ions as a New Form of Condensed Hydrogen

We report the first experimental observation of negatively charged hydrogen and deuterium cluster ions, H_{n}^{-} and D_{n}^{-}, where n≥5. These anions are formed by an electron addition to liquid helium nanodroplets doped with molecular hydrogen or deuterium. The ions are stable for at least the lifetime of the experiment, which is several tens of microseconds. Only anions with odd values of n are detected, and some specific ions show anomalously high abundances. The sizes of these "magic number" ions suggest an icosahedral framework of H_{2} (D_{2}) molecules in solvent shells around a central H^{-} (D^{-}) ion. The first three shells, which contain a total of 44 H_{2} or D_{2} molecules, appear to be solidlike, but thereafter a more liquidlike arrangement of the H_{2} (D_{2}) molecules is adopted.

Synthesis of sodium polyhydrides at high pressures

The only known compound of sodium and hydrogen is archetypal ionic NaH. Application of high pressure is known to promote states with higher atomic coordination, but extensive searches for polyhydrides with unusual stoichiometry have had only limited success in spite of several theoretical predictions. Here we report the first observation of the formation of polyhydrides of Na (NaH₃ and NaH₇) above 40 GPa and 2,000 K. We combine synchrotron X-ray diffraction and Raman spectroscopy in a laser-heated diamond anvil cell and theoretical random structure searching, which both agree on the stable structures and compositions. Our results support the formation of multicenter bonding in a material with unusual stoichiometry. These results are applicable to the design of new energetic solids and high-temperature superconductors based on hydrogen-rich materials.

The only known compound of sodium and hydrogen is ionic NaH, but theory predicts the existence of polyhydrides at high pressure. Here, the authors report observations of the formation of polyhydrides above 40 GPa and 2000 K, supporting the idea of multicentre bonding in a material with unusual stoichiometry.

Compressed cesium polyhydrides: Cs+ sublattices and H3(-) three-connected nets

The cesium polyhydrides (CsH(n), n > 1) are predicted to become stable, with respect to decomposition into CsH and H2, at pressures as low as 2 GPa. The CsH3 stoichiometry is found to have the lowest enthalpy of formation from CsH and H2 between 30 and 200 GPa. Evolutionary algorithms predict five distinct, mechanically stable, nearly isoenthalpic CsH3 phases consisting of H3(–) molecules and Cs+ atoms. The H3(–) sublattices in two of these adopt a hexagonal three-connected net; in the other three the net is twisted, like the silicon sublattice in the α-ThSi2 structure. The former emerge as being metallic below 100 GPa in our screened hybrid density functional theory calculations, whereas the latter remain insulating up to pressures greater than 250 GPa. The Cs+ cations in the most-stable I4(1)/amd CsH3 phase adopt the positions of the Cs atoms in Cs-IV, and the H3(–) molecules are found in the (interstitial) regions which display a maximum in the electron density.

Potential energy and dipole moment surfaces of HCO- for the search of H- in the interstellar medium

Potential energy and permanent dipole moment surfaces of the electronic ground state of formyl negative ion HCO(-) are determined for a large number of geometries using the coupled-cluster theory with single and double and perturbative treatment of triple excitations ab initio method with a large basis set. The obtained data are used to construct interpolated surfaces, which are extended analytically to the region of large separations between CO and H(-) with the multipole expansion approach. We have calculated the energy of the lowest rovibrational levels of HCO(-) that should guide the spectroscopic characterization of HCO(-) in laboratory experiments. The study can also help to detect HCO(-) in the cold and dense regions of the interstellar medium where the anion could be formed through the association of abundant CO with still unobserved H(-).

A systematic investigation of the ground state potential energy surface of H3+

Based on different ab initio electronic structure calculations (CI-R12 and Gaussian Geminals) of the Born-Oppenheimer electronic energy E(BO) of H(3)(+) from high to highest quality, we build up a potential energy surface which represents a highly reliable form of the topology of the whole potential region, locally and globally. We use the CI-R12 method in order to get within reasonable CPU-time a relatively dense grid of energy points. We demonstrate that CI-R12 is good enough to give an accurate surface, i.e., Gaussian Geminals are not absolutely necessary. For different types of potential energy surface fits, we performed variational calculations of all bound vibrational states, including resonances above the dissociation limit, for total angular momentum J = 0. We clarify the differences between different fits of the energy to various functional forms of the potential surface. Small rms-values (<1 cm(-1)) of the fit do not provide precise information about the interpolatory behaviour of the fit functions.

Rubidium polyhydrides under pressure: emergence of the linear H3(-) species

The structures of compressed rubidium polyhydrides, RbH(n) with n>1, and their evolution under pressure are studied using density functional theory calculations. These phases, which start to stabilize at only P = 2 GPa, consist of Rb(+) cations and one or more of the following species: H(-) anions, H(2) molecules, and H(3)(-) molecules. The latter motif, the simplest example of a three-center four-electron bond, is found in the most stable structures, RbH(5) and RbH(3) , which metallize above 200 GPa. At the highest pressures studied, our evolutionary searches find an RbH(6) phase which contains polymeric (H(3)(-))(∞) chains that show signs of one-dimensional liquid-like behavior.

Structure and stability of the negative hydrogen molecular ion

We present the results of a Coulomb explosion experiment that allows for the imaging of the rovibrational wave function of the metastable H2- ion. Our measurements confirm the predicted large internuclear separation of 6 a.u., and they show that the ion decays by autodetachment rather than by spontaneous dissociation. Imaging of the resulting H2 products reveals a large angular momentum of J = 25 ± 2, quantifying the rotation that leads to the metastability of this most fundamental molecular anion.