Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Effects of molecular rotation on densities in doped He4 clusters

Mehul V. Patel, Alexandra Viel, Francesco Paesani, Patrick Huang, K. Birgitta Whaley

Published Feb 28, 2003• 37 citations
Publisher

Abstract

The effects of including rotational degrees of freedom on helium solvation densities in molecule-doped helium clusters are investigated for a variety of molecules. Helium densities and cluster energetics are calculated with diffusion Monte Carlo methods. The rotationally induced changes in the helium density distributions are examined and quantified with a theoretical estimator applicable to molecules of arbitrary symmetry. This analysis leads to a discussion of adiabatic following of molecular rotation in a solvating helium environment. We make a detailed comparative study of the effect of molecular rotation as a function of four impurity molecules with varying mass and symmetry: SF6, OCS, HCN, and benzene (C6H6). We find that even for the heaviest rotors, only a fraction of the solvating helium density adiabatically follows the molecular motion in the quantum ground state. For the lightest molecule, HCN, a negligible degree of adiabatic following is found. A discussion of the various definitions is presented to clarify the meaning of adiabatic following, and its applicability to dynamical models of quantum rotation in helium droplets is evaluated in light of the quantitative findings of incomplete adiabatic following established here.

DOI
10.1063/1.1545106

Related papers

Rotational excitations of N2O in small helium clusters and the role of Bose permutation symmetry

Francesco Paesani, K. Birgitta Whaley

We present a detailed study of the energetics, structures, and Bose properties of small clusters of (4)He containing a single nitrous oxide (N(2)O) molecule, from N=1 (4)He up to sizes corresponding to completion of the first solvation shell around N(2)O (N=16 (4)He). Ground state properties are calculated using the importance-sampled rigid-body diffusion Monte Carlo method, rotational excited state calculations are made with the projection operator imaginary time spectral evolution method, and Bose permutation exchange and associated superfluid properties are calculated with the finite temperature path integral method. For N< or =5 the helium atoms are seen to form an equatorial ring around the molecular axis, at N=6 helium density starts to occupy the second (local) minimum of the N(2)O-He interaction at the oxygen side of the molecule, and N=9 is the critical size at which there is onset of helium solvation all along the molecular axis. For N> or =8 six (4)He atoms are distributed in a symmetric, quasirigid ring around N(2)O. Path integral calculations show essentially complete superfluid response to rotation about the molecular axis for N> or =5, and a rise of the perpendicular superfluid response from zero to appreciable values for N> or =8. Rotational excited states are computed for three values of the total angular momentum, J=1-3, and the energy levels fitted to obtain effective spectroscopic constants that show excellent agreement with the experimentally observed N dependence of the effective rotational constant B(eff). The non-monotonic behavior of the rotational constant is seen to be due to the onset of long (4)He permutation exchanges and associated perpendicular superfluid response of the clusters for N> or =8. We provide a detailed analysis of the role of the helium solvation structure and superfluid properties in determining the effective rotational constants.

The Journal of Chemical Physics•Sep 8, 2004•59 citations

Rotational fluctuation of molecules in quantum clusters. II. Molecular rotation and superfluidity in OCS-doped helium-4 clusters

Shinichi Miura

In this paper, quantum fluctuations of a carbonyl sulfide molecule in helium-4 clusters are studied as a function of cluster size N in a small-to-large size regime (2<or=N<or=64). The molecular rotation of the dopant shows nonmonotonic size dependence in the range of 10<or=N<or=20, reflecting the density distribution of the helium atoms around the molecule. The size dependence on the rotational constant shows a plateau for N>or=20, which is larger than the experimental nanodroplet value. Superfluid response of the doped cluster is found to show remarkable anisotropy especially for N<or=20. The superfluid fraction regarding the axis perpendicular to the molecular axis shows a steep increase at N=10, giving the significant enhancement of the rotational fluctuation of the molecule. On the other hand, the superfluid fraction regarding the axis parallel to the molecular axis reaches 0.9 at N=5, arising from the bosonic exchange cycles of the helium atoms around the molecular axis. The anisotropy in the superfluid response is found to be the direct consequence of the configurations of the bosonic exchange cycles.

The Journal of Chemical Physics•Mar 20, 2007•28 citations

Rotational Superfluidity in Small Helium Droplets

David Mateo, Frisly Gonzalez, Jussi Eloranta

The first minimum appearing in molecular rotational constants as a function of helium droplet size has been previously associated with the onset of superfluidity in these finite systems. We investigate this relationship by bosonic density functional theory calculations of classical molecular rotors (OCS, N2O, CO, and HCN) interacting with the surrounding helium. The calculated rotational constants are in fair agreement with the existing experimental data, demonstrating the applicability of the theoretical model. Inspection of the spatial evolution of the global phase and density shows the increase in the rotational constant after the first minimum correlates with continuous coverage of the molecule by helium and the appearance of angular phase coherence rather than completion of the first solvent shell. We assign the observed phenomenon to quantum phase transition between a localized state and one-dimensional superfluid, which represents the onset of rotational superfluidity in small helium droplets.

The Journal of Physical Chemistry A•Aug 6, 2014•24 citations

Rotational fluctuation of molecules in quantum clusters. II. Molecular rotation and superfluidity in OCS-doped helium-4 clusters

伸一 三浦, Shinichi Miura

In this paper, quantum fluctuations of a carbonyl sulfide molecule in helium-4 clusters are studied as a function of cluster size N in a small-to-large size regime (2 ≤ N ≤ 64). The molecular rotation of the dopant shows nonmonotonic size dependence in the range of 10≤N≤20, reflecting the density distribution of the helium atoms around the molecule. The size dependence on the rotational constant shows a plateau for N≥20, which is larger than the experimental nanodroplet value. Superfluid response of the doped cluster is found to show remarkable anisotropy especially for N ≤20. The superfluid fraction regarding the axis perpendicular to the molecular axis shows a steep increase at N = 10, giving the significant enhancement of the rotational fluctuation of the molecule. On the other hand, the superfluid fraction regarding the axis parallel to the molecular axis reaches 0.9 at N=5, arising from the bosonic exchange cycles of the helium atoms around the molecular axis. The anisotropy in the superfluid response is found to be the direct consequence of the configurations of the bosonic exchange cycles. © 2007 American Institute of Physics.

•Jan 1, 2007•0 citations

Computational Studies of Microscopic Superfluidity in the 4He Clusters

Angeline R. Wairegi

The physics that result in the decoupling of a molecule from a bosonic solvent at 0 K are studied. Fixed-node diffusion Monte Carlo (FNDMC) coupled with a Genetic Algorithm is used to perform simulations of the bosonic droplets doped with various molecules. The efficacy and accuracy of this approach is tested on a strongly coupled 2-dimensional quartic oscillator with excellent results. This algorithm is then applied to 4He-CO and 4He-HCN clusters respectively in an effort to determine the factors that result in the onset of microscopic superfluidity. The decoupling of the doped molecule from the bosonic solvent is found to be, primarily, a result of the combined effect of the repulsive interaction between the helium atoms and bose symmetry. The effects of rotor size versus molecular anisotropy in a NH3 molecule seeded into a 4He droplet is studied as well. Simulations are done using the accurate rotational constants (B0=9.945 cm-1, C0=6.229 cm-1) and using "fudged" versions of the rotational constants (Bfudged=0.9945 cm-1, Cfudged=0.6229 cm-1) for the |0011〉state. The simulations done with the fudged rotational constants experience a slightly smaller reduction than those done using the accurate rotational constants. This is attributed to the importance of molecular anisotropy versus the size of larger rotational constants in molecules whose rotational constants fall in an intermediate regime.

•Oct 8, 2021•0 citations

Quantum Simulations of Helium Clusters with Ionic and Open Shell Dopants

Mirjana Mladenović, Marius Lewerenz, Petr Slavíček

Many experiments on pure and doped helium clusters result in the production of a broad distribution of charged fragment ions after electron impact or laser ionisation. These ion yield distributions occasionally exhibit distinct stability patterns. We use the diffusion quantum Monte Carlo technique to study these stability patterns and structural features of helium clusters with open shell atomic and molecular dopants. These simulations require reliable many-body potentials which we construct from high level ab initio CCSD(T) calculations for several electronic states mixed by spin-orbit coupling and including non-additive interactions arising from induction. For small clusters numerically exact calculations of rovibrational properties are used to establish the quality of our potential surfaces. The approach is illustrated with applications to I$^{q}$@He$_n$, q=-1,0,+1,+2 clusters which appear in ionisation and photodissociation experiments. Spin-orbit coupling due to the iodine atom is found to have a profound effect on the structure of the mixed clusters. Radial and angular helium density profiles are computed and show that even in the presence of pronounced radial structuring exchange between 'shells' remains possible and angular disorder remains significant. The identification of a well defined ionic snowball core appears incorrect. We present first results for CO$^+$ ions in helium clusters. He-CO$^+$ is an interesting astrophysical collision system but its interest for helium cluster studies is the similarity of the CO$^+$ rotational constant with the one of neutral CO. Our present understanding of rotation inside helium clusters relies on studies of molecules where changing the molecule implies changing dynamical parameters and the interaction potential. The CO/CO$^+$ case allows to study the specific effect of changing only the interaction energy. CO$^+$@He$_n$ should be relatively easy to make by cluster ionisation or aggregation in ion drift tubes and its strong dipole and vibrational transition moment make it an easy target for high resolution spectroscopy. We have computed accurate ab initio surfaces for the two lowest electronic states of He-CO$^+$ to predict rovibrational spectroscopic and collisional properties. A many-body model using these surfaces is used to study larger CO$^+$He$_n$ clusters.

•Nov 28, 2011•0 citations
The effects of including rotational degrees of freedom on helium solvation densities in molecule-doped helium clusters are investigated for a variety of molecules. Helium densities and cluster energetics are calculated with diffusion Monte Carlo methods. The rotationally induced changes in the helium density distributions are examined and quantified with a theoretical estimator applicable to molecules of arbitrary symmetry. This analysis leads to a discussion of adiabatic following of molecular rotation in a solvating helium environment. We make a detailed comparative study of the effect of molecular rotation as a function of four impurity molecules with varying mass and symmetry: SF6, OCS, HCN, and benzene (C6H6). We find that even for the heaviest rotors, only a fraction of the solvating helium density adiabatically follows the molecular motion in the quantum ground state. For the lightest molecule, HCN, a negligible degree of adiabatic following is found. A discussion of the various definitions is presented to clarify the meaning of adiabatic following, and its applicability to dynamical models of quantum rotation in helium droplets is evaluated in light of the quantitative findings of incomplete adiabatic following established here.
We present a detailed study of the energetics, structures, and Bose properties of small clusters of (4)He containing a single nitrous oxide (N(2)O) molecule, from N=1 (4)He up to sizes corresponding to completion of the first solvation shell around N(2)O (N=16 (4)He). Ground state properties are calculated using the importance-sampled rigid-body diffusion Monte Carlo method, rotational excited state calculations are made with the projection operator imaginary time spectral evolution method, and Bose permutation exchange and associated superfluid properties are calculated with the finite temperature path integral method. For N< or =5 the helium atoms are seen to form an equatorial ring around the molecular axis, at N=6 helium density starts to occupy the second (local) minimum of the N(2)O-He interaction at the oxygen side of the molecule, and N=9 is the critical size at which there is onset of helium solvation all along the molecular axis. For N> or =8 six (4)He atoms are distributed in a symmetric, quasirigid ring around N(2)O. Path integral calculations show essentially complete superfluid response to rotation about the molecular axis for N> or =5, and a rise of the perpendicular superfluid response from zero to appreciable values for N> or =8. Rotational excited states are computed for three values of the total angular momentum, J=1-3, and the energy levels fitted to obtain effective spectroscopic constants that show excellent agreement with the experimentally observed N dependence of the effective rotational constant B(eff). The non-monotonic behavior of the rotational constant is seen to be due to the onset of long (4)He permutation exchanges and associated perpendicular superfluid response of the clusters for N> or =8. We provide a detailed analysis of the role of the helium solvation structure and superfluid properties in determining the effective rotational constants.
In this paper, quantum fluctuations of a carbonyl sulfide molecule in helium-4 clusters are studied as a function of cluster size N in a small-to-large size regime (2<or=N<or=64). The molecular rotation of the dopant shows nonmonotonic size dependence in the range of 10<or=N<or=20, reflecting the density distribution of the helium atoms around the molecule. The size dependence on the rotational constant shows a plateau for N>or=20, which is larger than the experimental nanodroplet value. Superfluid response of the doped cluster is found to show remarkable anisotropy especially for N<or=20. The superfluid fraction regarding the axis perpendicular to the molecular axis shows a steep increase at N=10, giving the significant enhancement of the rotational fluctuation of the molecule. On the other hand, the superfluid fraction regarding the axis parallel to the molecular axis reaches 0.9 at N=5, arising from the bosonic exchange cycles of the helium atoms around the molecular axis. The anisotropy in the superfluid response is found to be the direct consequence of the configurations of the bosonic exchange cycles.
The first minimum appearing in molecular rotational constants as a function of helium droplet size has been previously associated with the onset of superfluidity in these finite systems. We investigate this relationship by bosonic density functional theory calculations of classical molecular rotors (OCS, N2O, CO, and HCN) interacting with the surrounding helium. The calculated rotational constants are in fair agreement with the existing experimental data, demonstrating the applicability of the theoretical model. Inspection of the spatial evolution of the global phase and density shows the increase in the rotational constant after the first minimum correlates with continuous coverage of the molecule by helium and the appearance of angular phase coherence rather than completion of the first solvent shell. We assign the observed phenomenon to quantum phase transition between a localized state and one-dimensional superfluid, which represents the onset of rotational superfluidity in small helium droplets.
In this paper, quantum fluctuations of a carbonyl sulfide molecule in helium-4 clusters are studied as a function of cluster size N in a small-to-large size regime (2 ≤ N ≤ 64). The molecular rotation of the dopant shows nonmonotonic size dependence in the range of 10≤N≤20, reflecting the density distribution of the helium atoms around the molecule. The size dependence on the rotational constant shows a plateau for N≥20, which is larger than the experimental nanodroplet value. Superfluid response of the doped cluster is found to show remarkable anisotropy especially for N ≤20. The superfluid fraction regarding the axis perpendicular to the molecular axis shows a steep increase at N = 10, giving the significant enhancement of the rotational fluctuation of the molecule. On the other hand, the superfluid fraction regarding the axis parallel to the molecular axis reaches 0.9 at N=5, arising from the bosonic exchange cycles of the helium atoms around the molecular axis. The anisotropy in the superfluid response is found to be the direct consequence of the configurations of the bosonic exchange cycles. © 2007 American Institute of Physics.
The physics that result in the decoupling of a molecule from a bosonic solvent at 0 K are studied. Fixed-node diffusion Monte Carlo (FNDMC) coupled with a Genetic Algorithm is used to perform simulations of the bosonic droplets doped with various molecules. The efficacy and accuracy of this approach is tested on a strongly coupled 2-dimensional quartic oscillator with excellent results. This algorithm is then applied to 4He-CO and 4He-HCN clusters respectively in an effort to determine the factors that result in the onset of microscopic superfluidity. The decoupling of the doped molecule from the bosonic solvent is found to be, primarily, a result of the combined effect of the repulsive interaction between the helium atoms and bose symmetry. The effects of rotor size versus molecular anisotropy in a NH3 molecule seeded into a 4He droplet is studied as well. Simulations are done using the accurate rotational constants (B0=9.945 cm-1, C0=6.229 cm-1) and using "fudged" versions of the rotational constants (Bfudged=0.9945 cm-1, Cfudged=0.6229 cm-1) for the |0011〉state. The simulations done with the fudged rotational constants experience a slightly smaller reduction than those done using the accurate rotational constants. This is attributed to the importance of molecular anisotropy versus the size of larger rotational constants in molecules whose rotational constants fall in an intermediate regime.
Many experiments on pure and doped helium clusters result in the production of a broad distribution of charged fragment ions after electron impact or laser ionisation. These ion yield distributions occasionally exhibit distinct stability patterns. We use the diffusion quantum Monte Carlo technique to study these stability patterns and structural features of helium clusters with open shell atomic and molecular dopants. These simulations require reliable many-body potentials which we construct from high level ab initio CCSD(T) calculations for several electronic states mixed by spin-orbit coupling and including non-additive interactions arising from induction. For small clusters numerically exact calculations of rovibrational properties are used to establish the quality of our potential surfaces. The approach is illustrated with applications to Iq^{q}q@Hen_nn​, q=-1,0,+1,+2 clusters which appear in ionisation and photodissociation experiments. Spin-orbit coupling due to the iodine atom is found to have a profound effect on the structure of the mixed clusters. Radial and angular helium density profiles are computed and show that even in the presence of pronounced radial structuring exchange between 'shells' remains possible and angular disorder remains significant. The identification of a well defined ionic snowball core appears incorrect. We present first results for CO+^++ ions in helium clusters. He-CO+^++ is an interesting astrophysical collision system but its interest for helium cluster studies is the similarity of the CO+^++ rotational constant with the one of neutral CO. Our present understanding of rotation inside helium clusters relies on studies of molecules where changing the molecule implies changing dynamical parameters and the interaction potential. The CO/CO+^++ case allows to study the specific effect of changing only the interaction energy. CO+^++@Hen_nn​ should be relatively easy to make by cluster ionisation or aggregation in ion drift tubes and its strong dipole and vibrational transition moment make it an easy target for high resolution spectroscopy. We have computed accurate ab initio surfaces for the two lowest electronic states of He-CO+^++ to predict rovibrational spectroscopic and collisional properties. A many-body model using these surfaces is used to study larger CO+^++Hen_nn​ clusters.