Theoretical calculations of zero‐temperature absorption spectra of Li in solid H2

Spectroscopy of Alkali Atoms in Solid Matrices of Rare Gases: Experimental Results and Theoretical Analysis

We present an experimental and theoretical investigation of the spectroscopy of dilute alkali atoms in a solid matrix of inert gases at cryogenic temperatures, specifically Rubidium atoms in a solid Argon or Neon matrix, and related aspects of the interaction energies between the alkali atoms and the atoms of the solid matrix. The system considered is relevant for matrix isolation spectroscopy, and it is at the basis of a recently proposed detector of cosmological axions, exploiting magnetic-type transitions between Zeeman sublevels of alkali atoms in a magnetic field, tuned to the axion mass, assumed in the meV range. Axions are one of the supposed constituents of the dark matter (DM) of the Universe. This kind of spectroscopy could be also relevant for the experimental search of new physics beyond the Standard Model, in particular the search of violations of time-reversal or parity-charge-conjugation (CP) symmetry. In order to efficiently resolve the axion-induced transition in alkali-doped solid matrices, it is necessary to reduce as much as possible the spectral linewidth of the electronic transitions involved. The theoretical investigation presented in this paper aims to estimate the order of magnitude of the inhomogeneous contribution to the linewidth due to the alkali–matrix interactions (Coulomb/exchange and dispersion), and to compare the theoretical results with our experimental measurements of spectra of dilute Rubidium atoms in Argon and Neon solid matrix. The comparison of the expected or measured spectral linewidths will be important for selecting the most appropriate combination of alkali atoms and matrix inert elements to be used in the proposed axion detection scheme. It is finally suggested that dilute Lithium atoms diffused in a cold parahydrogen solid matrix could be, overall, a good system upon which the proposed detector could be based.

Equation of state and first principles prediction of the vibrational matrix shift of solid parahydrogen

We generate the equation of state (EOS) of solid parahydrogen (para-H₂) using a path-integral Monte Carlo (PIMC) simulation based on a highly accurate first-principles adiabatic hindered rotor potential energy curve for the para-H₂ dimer. The EOS curves for the fcc and hcp structures of solid para-H₂ near the equilibrium density show that the hcp structure is the more stable of the two, in agreement with experiment. To accurately reproduce the structural and energy properties of solid para-H₂, we eliminated by extrapolation the systematic errors associated with the choice of simulation parameters used in the PIMC calculation. We also investigate the temperature dependence of the EOS curves, and the invariance of the equilibrium density with temperature is satisfyingly reproduced. The pressure as a function of density and the compressibility as a function of pressure are both calculated using the obtained EOS and are compared with previous simulation results and experiments. We also report the first ever a priori prediction of a vibrational matrix shift from first-principles two-body potential functions, and its result for the equilibrium state agrees well with experiment.

Coinage metal exciplexes with helium atoms: a theoretical study of M*(2L)He(n) (M = Cu, Ag, Au; L = P,D)

The structure and energetics of exciplexes M*((2)L)He(n) (M = Cu, Ag and Au; L = P and D) in their vibrational ground state are studied by employing diffusion Monte Carlo (DMC). Interaction potentials between the excited coinage metals and He atoms are built using the Diatomics-in-Molecule (DIM) approach and ab initio potential curves for the M((2)L)-He dimers. Extending our previous work [Cargnoni et al., J. Phys. Chem. A, 2011, 115, 7141], we computed the dimer potential for Au in the (2)P and (2)D states, as well for Cu and Ag in the (2)D state, employing basis set superposition error-corrected Configuration Interaction calculations. We found that the (2)Π potential correlating with the (2)P state of Au is substantially less binding than for Ag and Cu, a trend well supported by the M(+) ionic radiuses. Conversely, the interaction potentials between a (n - 1)d(9)ns(2 2)D metal and He present a very weak dependency on M itself or the projection of the angular momentum along the dimer axis. This is due to the screening exerted by the ns(2) electrons on the hole in the (n - 1)d shell. Including the spin-orbit coupling perturbatively in the DIM energy matrix has a major effect on the lowest potential energy surface of the (2)P manifold, the one for Cu allowing the formation of a "belt" of five He atoms while the one for Au being completely repulsive. Conversely, spin-orbit coupling has only a weak effect on the (2)D manifold due to the nearly degenerate nature of the diatomic potentials. Structural and energetic results from DMC have been used to support experimental indications for the formation of metastable exciplexes or the opening of non-radiative depopulation channels in bulk and cold gaseous He.

Infrared-active spin-orbit transitions of halogen atom dopants in solid parahydrogen: the role of trapping site geometry

We present theoretical calculations of the (2)P(1/2) ← (2)P(3/2) spin-orbit transition of Cl dopants embedded as substitutional impurities in solid parahydrogen (pH2) matrices. In the lower-energy (2)P(3/2) spin-orbit level, the Cl atom's electron density distribution is anisotropic, and slightly distorts the geometry of the atom's trapping site. This distortion leads to a blue shift in the spin-orbit transition energy; the blue shift is enhanced when we account for the large-amplitude zero point motions of the pH2 molecules surrounding the Cl dopant. We also show that the intensity of the transition depends on the geometry of the trapping site. In the gas phase, the (2)P(1/2) ← (2)P(3/2) atomic transition is electric dipole forbidden. However, when the Cl atom resides in trapping sites that mimic the hexagonal close packed morphology of pure solid pH2, the transition becomes electric dipole allowed through interaction-induced transition dipole moments. These transition dipole moments originate in the anisotropic electron density distribution of the lower-energy (2)P(3/2) spin-orbit level.

Ca impurity in small mixed (4)He-(3)He clusters

The structure of small mixed helium clusters doped with one calcium atom has been determined within the diffusion Monte Carlo framework. The results show that the calcium atom sits at the (4)He-(3)He interface. This is in agreement with previous studies, both experimental and theoretical, performed for large clusters. A comparison between the results obtained for the largest cluster we have considered for each isotope shows a clear tendency of the Ca atom to reside in a deep dimple at the surface of the cluster for (4)He clusters, and to become fully solvated for (3)He clusters. We have calculated the absorption spectrum of Ca around the 4s4p<--4s(2) transition and have found that it is blueshifted from that of the free-atom transition by an amount that depends on the size and composition of the cluster.

The ultrafast structural response of solid parahydrogen: A complementary experimental/simulation investigation

We present a complete characterization, based on femtosecond pump-probe spectroscopy and molecular dynamics simulations, of the ultrafast dynamics of electronic bubble formation in solid parahydrogen upon impulsive excitation of impurity-doped sites, which correlate with the lowest Rydberg state of the NO impurity. The high temporal resolution of the experiment allows us to identify three time scales in the structural dynamics. A first ultrafast expansion (<150 fs), associated with the release of approximately 80% of the excess energy available to the system after excitation, is accompanied by a transient narrowing of the spatial distribution of the first shell of H2 molecules around the impurity. In a subsequent stage (up to approximately 800 fs), the cavity expansion slows down, and energy starts to flow irreversibly into the crystal. Finally, the lattice undergoes a slow structural reorganization at the impurity site (5-10 ps). A weak low-frequency recurrence, probably associated with an elastic response of the crystal, is observed at approximately 10 ps. The absence of polarization dependence indicates that the dynamics is largely dominated by translational (radial) motions of the molecules surrounding NO and not by the rotational motion of the impurity. Molecular dynamics simulations with temperature corrections, to mimic zero-point fluctuations, fully support the experimental results and show that the bubble model is suited to describe the dynamics of the system. It appears that the response of the medium around the impurity at short times is typical of a liquid solvent rather than that of a solid.

The OH radical-H2O molecular interaction potential

The OH radical is one of the most important oxidants in the atmosphere due to its high reactivity. The study of hydrogen-bonded complexes of OH with the water molecules is a topic of significant current interest. In this work, we present the development of a new analytical functional form for the interaction potential between the rigid OH radical and H(2)O molecules. To do this we fit a selected functional form to a set of high level ab initio data. Since there is a low-lying excited state for the H(2)O.OH complex, the impact of the excited state on the chemical behavior of the OH radical can be very important. We perform a potential energy surface scan using the CCSD(T)/aug-cc-pVTZ level of electronic structure theory for both excited and ground states. To model the physics of the unpaired electron in the OH radical, we develop a tensor polarizability generalization of the Thole-type all-atom polarizable rigid potential for the OH radical, which effectively describes the interaction of OH with H(2)O for both ground and excited states. The stationary points of (H(2)O)(n)OH clusters were identified as a benchmark of the potential.

Path integral Monte Carlo simulation of the absorption spectra of an Al atom embedded in helium

We use a multilevel path integral Monte-Carlo (PIMC) method to simulate the arrangement of He atoms around a single Al atom doped in a He cluster. High-level ab initio Al-He pair potentials and a Balling and Wright pairwise Hamiltonian model are used to describe the full potential and the electronic asymmetry arising from the open-shell character of the Al atom in its ground and excited electronic states. Our calculations show that the doping of the Al 3p electron strongly influences the He packing. The results of the PIMC simulation are used to predict the electronic excitation spectrum of an Al atom embedded in He clusters. With inclusion of tail corrections for the ground and excited states potentials, the calculated 3d<--3p spectrum agrees reasonably well with the experimental spectrum. The blueshift of the calculated spectrum associated with the 4s<--3p transition of solvated Al is about 25 nm (2000 cm-1) larger than seen in experiments on Al embedded in bulk liquid He. We predict that the spectrum associated with the 4p<--3p transition will be blueshifted by approximately 7000 cm-1 (nearly 1 eV).

Predicting atomic dopant solvation in helium clusters: The MgHen case

We present a quantum Monte Carlo study of the solvation and spectroscopic properties of the Mg-doped helium clusters MgHe(n) with n=2-50. Three high-level [MP4, CCSD(T), and CCSDT] MgHe interaction potentials have been used to study the sensitivity of the dopant location on the shape of the pair interaction. Despite the similar MgHe well depth, the pair distribution functions obtained in the diffusion Monte Carlo simulations markedly differ for the three pair potentials, therefore indicating different solubility properties for Mg in He(n). Moreover, we found interesting size effects for the behavior of the Mg impurity. As a sensitive probe of the solvation properties, the Mg excitation spectra have been simulated for various cluster sizes and compared with the available experimental results. The interaction between the excited 1P Mg atom and the He moiety has been approximated using the diatomics-in-molecules method and the two excited 1pi and 1sigma MgHe potentials. The shape of the simulated MgHe50 spectra shows a substantial dependency on the location of the Mg impurity, and hence on the MgHe pair interaction employed. To unravel the dependency of the solvation behavior on the shape of the computed potentials, exact density-functional theory has been adapted to the case of doped He(n) and various energy distributions have been computed. The results indicate the shape of the repulsive part of the MgHe potential as an important cause of the different behaviors.

Lattice response of quantum solids to an impulsive local perturbation

The lattice response of solid para-H2 to an impulsive electronic excitation was studied using femtosecond pump-probe spectroscopy. The evolution of an electronic bubble in the crystal, created upon excitation of the A(3ssigma) Rydberg state of an NO impurity, was followed in real time, with a resolution of 100 fs. The experimental results, interpreted in connection with molecular dynamics simulations with quantum corrections, indicate the presence of three stages in the dynamics: a sub-100 fs "adiabatic" phase, a 0.5-1 ps phase, corresponding to the interaction of the first with the next shells driven by the bubble expansion, and a 5 ps phase, corresponding to a slow rearrangement of the environment surrounding the impurity. These findings indicate that the lattice response in solid para-H2 resembles that of a liquid.

Zero-temperature equation of state of quasi-one-dimensional H2

We have studied molecular hydrogen in a pure 1D geometry and inside a narrow carbon nanotube by means of the diffusion Monte Carlo method. The one dimensionality of H2 in the nanotube is well maintained in a large density range, this system being closer to an ideal 1D fluid than liquid 4He in the same setup. H2 shares with 4He the existence of a stable liquid phase and a quasicontinuous liquid-solid transition at very high linear densities.

Evidence for Superfluidity in Para-Hydrogen Clusters Inside Helium-4 Droplets at 0.15 Kelvin

A linear carbonyl sulfide (OCS) molecule surrounded by 14 to 16 para-hydrogen (pH(2)) molecules, or similar numbers of ortho-deuterium (oD(2)) molecules, within large helium-4 ((4)He) droplets and inside mixed (4)He/(3)He droplets was investigated by infrared spectroscopy. In the pure (4)He droplets (0.38 kelvin), both systems exhibited spectral features that indicate the excitation of angular momentum around the OCS axis. In the colder (4)He/(3)He droplets (0.15 kelvin), these features remained in the oD(2) cluster spectra but disappeared in the pH(2) spectra, indicating that the angular momentum is no longer excited. These results are consistent with the onset of superfluidity, thereby providing the first evidence for superfluidity in a liquid other than helium.