Path integral Monte Carlo applications to quantum fluids in confined geometries

Carbon Nanotubes Immersed in Superfluid Helium: The Impact of Quantum Confinement on Wetting and Capillary Action

A recent experimental study [ Ohba, Sci. Rep. 2016, 6, 28992 ] of gas adsorption on single-walled carbon nanotubes at temperatures between 2 and 5 K reported a quenched propagation of helium through carbon nanotubes with diameters below 7 Å despite the small kinetic diameter of helium atoms. After assessing the performance of a potential model for the He-nanotube interaction via ab initio calculations with density functional theory-based symmetry adapted perturbation theory, we apply orbital-free helium density functional theory to show that the counterintuitive experimental result is a consequence of the exceptionally high zero-point energy of helium and its tendency to form spatially separated layers of helium upon adsorption at the lowest temperatures. Helium filling factors are derived for a series of carbon nanotubes and compared to the available experimental data.

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.

Spectral shifts and helium configurations in 4He(N)-tetracene clusters

Spectral shifts of electronic transitions of tetracene in helium droplets are investigated in a theoretical study of (4)He(N)-tetracene clusters with 1 < or = N < or = 150. Utilizing a pairwise interaction for the S(0) state of tetracene with helium that is extended by semiempirical terms to construct a potential for the S(1) state of tetracene with helium, the spectral shift is calculated from path integral Monte Carlo calculations of the helium equilibrium properties with tetracene in the S(0) and S(1) states at T = 0 and at T = 0.625 K. The calculated spectral shifts are in quantitative agreement with available experimental measurements for small values of N (< or = 8) at T approximately 0.4 K and show qualitative agreement for larger N (10-20). The extrapolated value of the spectral shift in large droplets (N approximately 10(4)) is approximately 90% of the experimentally measured value. We find no evidence of multiple configurations of helium for any cluster size for either the S(0) or S(1) state of tetracene. These results suggest that the observed spectral splitting of electronic transitions of tetracene in large helium droplets is not due to the coexistence of static metastable helium densities, unlike the situation previously analyzed for the phthalocyanine molecule.

Multiple solvation configurations around phthalocyanine in helium droplets

Recent measurements of the emission spectrum of phthalocyanine solvated in superfluid helium nanodroplets exhibit a constant 10.3 cm(-1) splitting of each emission line relative to the absorption spectrum. This splitting has been attributed to two distinct helium environments near the surface of the phthalocyanine molecule. Rigid-body path-integral Monte Carlo provides a means of investigating the origin of the splitting on a detailed microscopic level. Path-integral Monte Carlo simulations of 4He(N)-phthalocyanine at 0.625 K with N ranging from 24 to 150 show two distinct helium configurations. One configuration is commensurate with the molecular substrate and the other is a triangular lattice. We investigate the energetics of these two configurations and use a method for calculating electronic spectral shifts for aromatic molecule-rare-gas clusters due to dispersive interactions to estimate the spectral splitting that would arise from the two helium configurations seen for N=150. The results are in reasonable agreement with the experimentally measured splitting, supporting the existence of two distinct local helium environments near the surface of the molecule in the nanodroplets.

Low variance energy estimators for systems of quantum Drude oscillators: Treating harmonic path integrals with large separations of time scales

In the effort to develop atomistic models capable of accurately describing nanoscale systems with complex interfaces, it has become clear that simple treatments with rigid charge distributions and dispersion coefficients selected to generate bulk properties are insufficient to predict important physical properties. The quantum Drude oscillator model, a system of one-electron pseudoatoms whose "pseudoelectrons" are harmonically bound to their respective "pseudonuclei," is capable of treating many-body polarization and dispersion interactions in molecular systems on an equal footing due to the ability of the pseudoatoms to mimic the long-range interactions that characterize real materials. Using imaginary time path integration, the Drude oscillator model can, in principle, be solved in computer operation counts that scale linearly with the number of atoms in the system. In practice, however, standard expressions for the energy and pressure, including the commonly used virial estimator, have extremely large variances that require untenably long simulation times to generate converged averages. In this paper, low-variance estimators for the internal energy are derived, in which the large zero-point energy of the oscillators does not contribute to the variance. The new estimators are applicable to any system of harmonic oscillators coupled to one another (or to the environment) via an arbitrary set of anharmonic interactions. The variance of the new estimators is found to be much smaller than standard estimators in three example problems, a one-dimensional anharmonic oscillator and quantum Drude models of the xenon dimer and solid (fcc) xenon, respectively, yielding 2-3 orders of magnitude improvement in computational efficiency.

Spatial Delocalization in para-H2 Clusters

We applied the quantum path integral Monte Carlo method for the study of (para-H)N (N = 5-33) clusters at T = 2 K, exploring static and dynamic order, which originates from the effects of zero-point energy, kinetic energy, and thermal fluctuations in quantum clusters. Information on dynamic structure was inferred from the asymptotic tails of the cage correlation function calculated from the centroid Monte Carlo trajectory. The centroid cage correlation function decays to zero for large clusters (N = 15-33), manifesting the interchange of molecules between different solvation shells, with statistically diminishing back interchange. Further evidence for the floppiness of para-hydrogen clusters emerges from the Monte Carlo evolution of the centroid of a tagged molecule, which exhibits significant changes in the list of its first and second solvation shells due to the interchange of molecules between these shells.