Persistent molecular superfluid response in doped para-hydrogen clusters

Theory cracks old data: Rovibrational energy levels of orthoH2-CO derived from experiment

Measurements of rovibrational spectra of clusters provide physical insight only if spectral lines can be assigned to pairs of quantum states, and further insight is obtained if one can deduce the quantitative energy-level pattern. Both steps can be so difficult that some measured spectra remain unassigned, one example is orthoH2-CO. To extend the scope of spectroscopic insights, we propose to use theoretical information in interpretation of spectra. We first performed high accuracy, full-dimensional calculations of the orthoH2-CO spectrum, at the highest practically achievable levels of electronic structure theory and quantum nuclear dynamics. Then, an iterative, theory-guided method developed here allowed us to fully interpret the spectrum of orthoH2-CO, extending the range of van der Waals clusters for which spectroscopy can provide physical insights.

MLRNet: Combining the Physics-Motivated Potential Models with Neural Networks for Intermolecular Potential Energy Surface Construction

A physics-based machine learning model called MLRNet has been developed to construct the high-accuracy two-body intermolecular potential energy surface (IPES). The outputs of the neural network are integrated into the physically realistic Morse/long-range (MLR) function, which ensures that the MLRNet has meaningful extrapolation at both short and long ranges and solves the asymptotic problem in common neural network potential (NNP) models. The neural network representation of the MLR parameters is more flexible and more efficient than the polynomial expansion in the conventional mdMLR model, especially for systems containing nonrigid monomer(s). The present work illustrates the basic framework of the current MLRNet model, including (i) how to combine the physically meaningful MLR function with different possible NN structures, (ii) the preservation of permutation symmetry, and (iii) the predetermination of the long-range function u_LR. We choose two realistic systems to demonstrate the performance of MLRNet: the three-dimensional IPES of CO₂-He including the CO₂ antisymmetric vibration Q₃ and the six-dimensional IPES of the H₂O-Ar system. In both cases, the fitting errors of the MLRNet are several times smaller than those of the conventional mdMLR model. Both short-range and long-range extrapolation tests were performed to illustrate the extrapolation ability of the MLRNet and its damping function version. Moreover, for the 6-D H₂O-Ar system, the MLRNet only needs 1596 trainable parameters, which is almost equal to the number needed for the 5-D mdMLR model (1509) and half that needed for the PIP-NN model (3501) within similar accuracy, which illustrates the model efficiency in high-dimensional IPES fitting.

Neural network interaction potentials for para-hydrogen with flexible molecules

The study of molecular impurities in para-hydrogen ( pH2) clusters is key to push forward our understanding of intra- and intermolecular interactions, including their impact on the superfluid response of this bosonic quantum solvent. This includes tagging with only one or very few pH2, the microsolvation regime for intermediate particle numbers, and matrix isolation with many solvent molecules. However, the fundamental coupling between the bosonic pH2 environment and the (ro-)vibrational motion of molecular impurities remains poorly understood. Quantum simulations can, in principle, provide the necessary atomistic insight, but they require very accurate descriptions of the involved interactions. Here, we present a data-driven approach for the generation of impurity⋯ pH2 interaction potentials based on machine learning techniques, which retain the full flexibility of the dopant species. We employ the well-established adiabatic hindered rotor (AHR) averaging technique to include the impact of the nuclear spin statistics on the symmetry-allowed rotational quantum numbers of pH2. Embedding this averaging procedure within the high-dimensional neural network potential (NNP) framework enables the generation of highly accurate AHR-averaged NNPs at coupled cluster accuracy, namely, explicitly correlated coupled cluster single, double, and scaled perturbative triples, CCSD(T*)-F12a/aVTZcp, in an automated manner. We apply this methodology to the water and protonated water molecules as representative cases for quasi-rigid and highly flexible molecules, respectively, and obtain AHR-averaged NNPs that reliably describe the corresponding H2O⋯ pH2 and H3O+⋯ pH2 interactions. Using path integral simulations, we show for the hydronium cation, H3O+, that umbrella-like tunneling inversion has a strong impact on the first and second pH2 microsolvation shells. The automated and data-driven nature of our protocol opens the door to the study of bosonic pH2 quantum solvation for a wide range of embedded impurities.

Electronic spectroscopy of Mg-phthalocyanine embedded in cold hydrogen clusters produced by a pulsed nozzle

Cold clusters of molecular hydrogen were created using a pulsed nozzle. The thermodynamical states of the clusters were characterized by measuring the cluster beam velocity and the laser-induced fluorescence (LIF) spectra of embedded molecules. Two distinct velocity components were identified in the beam that originates from different clustering mechanisms. The fast velocity component corresponds to the expansion of H₂ from the gas phase, while the slow velocity component corresponds to the expansion from the liquid phase. The velocity distribution of these two components showed no significant difference between the expansions of para and normal hydrogen. In this study, LIF spectroscopy of single Mg-phthalocyanine molecules embedded in the H₂ clusters consisting of 10⁵ H₂ molecules was used to investigate the properties of the fast component. The observed peak frequencies of the LIF signals, compared to those observed in helium droplets, were used to infer the possible presence of the liquid phase in the fast component of the H₂ clusters below 5 K. The shift, linewidth, and splitting in the spectra, which strongly depend on the ortho/para ratio, are attributed to the local configurations of hydrogen in the vicinity of the probe molecules.

Quantum Microsolvation of Protonated Methane with ^{4}He: Large-Amplitude Motion Heavily Influences Bosonic Exchange

Quantum simulations of small CH_{5}^{+}·^{4}He_{n} complexes disclose significant and antagonistic impact of small-amplitude local vibrational motion vs large-amplitude global fluxional motion within the CH_{5}^{+} impurity on helium in real and permutation space. While the former significantly enhances bosonic exchange in the surrounding ^{4}He microsolvation shell compared to the rigid-body reference, the latter greatly suppresses long permutation cycles, which is traced back to the different nature of these quantum fluctuations. Therefore, it is expected that the resulting impact on local superfluidity is generic for fluctuating impurities in bosonic environments.

Suppression of Parahydrogen Superfluidity in a Doped Nanoscale Bose Fluid Mixture

Helium (^{4}He) nanodroplets provide a unique environment to observe the microscopic origins of superfluidity. The search for another superfluid substance has been an ongoing quest in the field of quantum fluids. Nearly two decades ago, experiments on doped parahydrogen (p-H_{2}) clusters embedded in ^{4}He droplets displayed anomalous spectroscopic signatures that were interpreted as a sign of the superfluidity of p-H_{2} [S. Grebenev et al., Science 289, 1532 (2000)SCIEAS0036-807510.1126/science.289.5484.1532]. Here, we observe, using first-principles quantum Monte Carlo simulations, a phase separation between a symmetric and localized p-H_{2} core and ^{4}He shells. The p-H_{2} core has minimal superfluid response. These findings are consistent with the recorded spectra but not with their original interpretation, and lead us to conclude that doped p-H_{2} clusters form a nonsuperfluid core in ^{4}He droplets.

Lieb-Liniger-like model of quantum solvation in CO-(4)HeN clusters

Small (4)He clusters doped with various molecules allow for the study of "quantum solvation" as a function of cluster size. A peculiarity of quantum solvation is that, as the number of (4)He atoms is increased from N = 1, the solvent appears to decouple from the molecule which, in turn, appears to undergo free rotation. This is generally taken to signify the onset of "microscopic superfluidity." Currently, little is known about the quantum mechanics of the decoupling mechanism, mainly because the system is a quantum (N + 1)-body problem in three dimensions which makes computations difficult. Here, a one-dimensional model is studied in which the (4)He atoms are confined to revolve on a ring and encircle a rotating CO molecule. The Lanczos algorithm is used to investigate the eigenvalue spectrum as the number of (4)He atoms is varied. Substantial solvent decoupling is observed for as few as N = 5 (4)He atoms. Examination of the Hamiltonian matrix, which has an almost block diagonal structure, reveals increasingly weak inter-block (solvent-molecule) coupling as the number of (4)He atoms is increased. In the absence of a dopant molecule the system is similar to a Lieb-Liniger (LL) gas and we find a relatively rapid transition to the LL limit as N is increased. In essence, the molecule initially-for very small N-provides a central, if relatively weak, attraction to organize the cluster; as more (4)He atoms are added, the repulsive interactions between the identical bosons start to dominate as the solvation ring (shell) becomes more crowded which causes the molecule to start to decouple. For low N, the molecule pins the atoms in place relative to itself; as N increases the atom-atom repulsion starts to dominate the Hamiltonian and the molecule decouples. We conclude that, while the notion of superfluidity is a useful and correct description of the decoupling process, a molecular viewpoint provides complementary insights into the quantum mechanism of the transition from a molecular cluster to a quantum solvated molecule.

Microwave spectroscopy of the seeded binary and ternary clusters CO-(pH2)2, CO-pH2-He, CO-HD, and CO-(oD2)(N=1,2)

We report the Fourier transform microwave spectra of the a-type J = 1-0 transitions of the binary and ternary CO-(pH2)2, CO-pH2-He, CO-HD, and CO-(oD2)N=1,2 clusters. In addition to the normal isotopologue of CO for all clusters, we observed the transitions of the minor isotopologues, (13)C(16)O, (12)C(18)O, and (13)C(18)O, for CO-(pH2)2 and CO-pH2-He. All transitions lie within 335 MHz of the experimentally or theoretically predicted values. In comparison to previously reported infrared spectra [Moroni et al., J. Chem. Phys. 122, 094314 (2005)], we are able to tentatively determine the vibrational shift for CO-pH2-He, in addition to its b-type J = 1-0 transition frequency. The a-type frequency of CO-pH2-He is similar to that of CO-He2 [Surin et al., Phys. Rev. Lett. 101, 233401 (2008)], suggesting that the pH2 molecule has a strong localizing effect on the He density. Perturbation theory analysis of CO-oD2 reveals that it is approximately T-shaped, with an anisotropy of the intermolecular potential amounting to ∼9 cm(-1).

Size dependent transition to solid hydrogen and argon clusters probed via spectroscopy of PTCDA embedded in helium nanodroplets

Complexes made of either Ar(N) or (H2)N clusters (N = 1-170) and a single PTCDA molecule (3,4,9,10-perylene-tetracarboxylic-dianhydride) are assembled inside helium droplets and spectroscopically studied via laser-induced fluorescence spectroscopy. The frequency shift and line-broadening are analyzed as a function of N and of the pick-up order of the PTCDA and cluster material in order to track liquid or solid properties of the clusters. For argon, the solid phase is observed for N > 10 above which the pick-up order dramatically influences the localization of the chromophore with respect to the Ar cluster. If the droplets are doped first with Ar, the chromophore remains on the surface of a solid cluster whereas for the reversed pick-up order the molecule is surrounded by an argon shell. At N < 10 wetting and the formation of the first solvation shell are observed. For para-hydrogen, a transition to the solid is observed at N ~ 20-25, confirming previous theoretical predictions on the existence of a liquid-like phase at such small sizes, even below the bulk hydrogen freezing temperature.

Rotationally adiabatic pair interactions of para- and ortho-hydrogen with the halogen molecules F2, Cl2, and Br2

The present work is concerned with the weak interactions between hydrogen and halogen molecules, i.e., the interactions of pairs H2-X2 with X = F, Cl, Br, which are dominated by dispersion and quadrupole-quadrupole forces. The global minimum of the four-dimensional (4D) coupled cluster with singles and doubles and perturbative triples (CCSD(T)) pair potentials is always a T shaped structure where H2 acts as the hat of the T, with well depths (De) of 1.3, 2.4, and 3.1 kJ/mol for F2, Cl2, and Br2, respectively. MP2/AVQZ results, in reasonable agreement with CCSD(T) results extrapolated to the basis set limit, are used for detailed scans of the potentials. Due to the large difference in the rotational constants of the monomers, in the adiabatic approximation, one can solve the rotational Schrödinger equation for H2 in the potential of the X2 molecule. This yields effective two-dimensional rotationally adiabatic potential energy surfaces where pH2 and oH2 are point-like particles. These potentials for the H2-X2 complexes have global and local minima for effective linear and T-shaped complexes, respectively, which are separated by 0.4-1.0 kJ/mol, where oH2 binds stronger than pH2 to X2, due to higher alignment to minima structures of the 4D-pair potential. Further, we provide fits of an analytical function to the rotationally adiabatic potentials.

Path-integral Monte Carlo method for Rényi entanglement entropies

We introduce a quantum Monte Carlo algorithm to measure the Rényi entanglement entropies in systems of interacting bosons in the continuum. This approach is based on a path-integral ground state method that can be applied to interacting itinerant bosons in any spatial dimension with direct relevance to experimental systems of quantum fluids. We demonstrate how it may be used to compute spatial mode entanglement, particle partitioned entanglement, and the entanglement of particles, providing insights into quantum correlations generated by fluctuations, indistinguishability, and interactions. We present proof-of-principle calculations and benchmark against an exactly soluble model of interacting bosons in one spatial dimension. As this algorithm retains the fundamental polynomial scaling of quantum Monte Carlo when applied to sign-problem-free models, future applications should allow for the study of entanglement entropy in large-scale many-body systems of interacting bosons.

Abnormal superfluid fraction of harmonically trapped few-fermion systems

Superfluidity is a fascinating phenomenon that, at the macroscopic scale, leads to dissipationless flow and the emergence of vortices. While these macroscopic manifestations of superfluidity are well described by theories that have their origin in Landau's two-fluid model, our microscopic understanding of superfluidity is far from complete. Using analytical and numerical ab initio approaches, this Letter determines the superfluid fraction and local superfluid density of small harmonically trapped two-component Fermi gases as a function of the interaction strength and temperature. At low temperature, we find that the superfluid fraction is, in certain regions of the parameter space, negative. This counterintuitive finding is traced back to the symmetry of the system's ground state wave function, which gives rise to a diverging quantum moment of inertia I(q). Analogous abnormal behavior of I(q) has been observed in even-odd nuclei at low temperature. Our predictions can be tested in modern cold atom experiments.

Microscopic superfluidity in 4He clusters stirred by a rotating impurity molecule

The effective moment of inertia of a CO impurity molecule in 4HeN and p-(H2)N solvent clusters initially increases with N but then commences a nonclassical decrease at N=4 (4He) or N=6 (p-H2). This suggests molecule-solvent decoupling and a transition to microscopic superfluidity. However, the quantum decoupling mechanism has not been elucidated. To understand the decoupling mechanism, a one-dimensional model is introduced in which the 4He atoms are confined to a ring. This model captures the physics and shows that decoupling happens primarily because of bosonic solvent-solvent repulsion. Quantum Monte Carlo and basis set calculations suggest that the system can be modeled as a stirred Tonks-Girardeau gas. This allows the N-particle time-dependent Schrödinger equation to be solved directly. Computations of the integrated particle current reveal a threshold for stirring and current generation, indicative of superfluidity.

Microscopic molecular superfluid response: theory and simulations

Since its discovery in 1938, superfluidity has been the subject of much investigation because it provides a unique example of a macroscopic manifestation of quantum mechanics. About 60 years later, scientists successfully observed this phenomenon in the microscopic world though the spectroscopic Andronikashvili experiment in helium nano-droplets. This reduction of scale suggests that not only helium but also para-H2 (pH2) can be a candidate for superfluidity. This expectation is based on the fact that the smaller number of neighbours and surface effects of a finite-size cluster may hinder solidification and promote a liquid-like phase. The first prediction of superfluidity in pH2 clusters was reported in 1991 based on quantum Monte Carlo simulations. The possible superfluidity of pH2 was later indirectly observed in a spectroscopic Andronikashvili experiment in 2000. Since then, a growing number of studies have appeared, and theoretical simulations have been playing a special role because they help guide and interpret experiments. In this review, we go over the theoretical studies of pH2 superfluid clusters since the experiment of 2000. We provide a historical perspective and introduce the basic theoretical formalism along with key experimental advances. We then present illustrative results of the theoretical studies and comment on the possible future developments in the field. We include sufficient theoretical details such that the review can serve as a guide for newcomers to the field.

Probing the Superfluid Response of para-Hydrogen with a Sulfur Dioxide Dopant

We recently presented the first attempt at using an asymmetric top molecule (para-water) to probe the superfluidity of nanoclusters (of para-hydrogen) [ Zeng , T. ; Li , H. ; Roy , P.-N. J. Phys. Chem. Lett. 2013 , 4 , 18 - 22 ]. Unfortunately, para-water could not be used to probe the para-hydrogen superfluid response. We now report a theoretical simulation of sulfur dioxide rotating in para-hydrogen clusters and show that this asymmetric top can serve as a genuine probe of superfluidity. With this probe, we predict that as few as four para-hydrogen molecules are enough to form a superfluid cluster, the smallest superfluid system to date. We also propose the concept of "exchange superfluid fraction" as a more precise measurement. New superfluid scenarios brought about by an asymmetric top dopant and potential experimental measurements are discussed.

Simulating Asymmetric Top Impurities in Superfluid Clusters: A para-Water Dopant in para-Hydrogen

We present the first simulation study of bosonic clusters doped with an asymmetric top molecule. The path-integral Monte Carlo method with the latest methodological advance in treating rigid-body rotation [Noya, E. G.; Vega, C.; McBride, C. J. Chem. Phys.2011, 134, 054117] is employed to study a para-water impurity in para-hydrogen clusters with up to 20 para-hydrogen molecules. The growth pattern of the doped clusters is similar in nature to that of pure clusters. The para-water molecule appears to rotate freely in the cluster. The presence of para-water substantially quenches the superfluid response of para-hydrogen with respect to the space-fixed frame.