Quantum structural fluctuation in para-hydrogen clusters revealed by the variational path integral method

Ground state chemical potential of parahydrogen clusters of size N = 21-40

We report the ground state chemical potential of parahydrogen clusters between N = 21-40 calculated using the Langevin equation Path Integral Ground State method. There has been much debate in the past whether the chemical potential size evolution in this region is jagged (indicating magic number cluster sizes) or if it is smooth (indicating some quantum melting below 1 K). We compare to previous diffusion Monte Carlo and Path Integral Ground State (PIGS) results, including very recent Variational Path Integral Molecular Dynamics (VPIMD) calculations [S. Miura, J. Chem. Phys. 148, 102333 (2018)]. We find that the ground state chemical potential is not a smooth curve and that magic number clusters are present, consistent with VPIMD and PIGS Monte Carlo results.

A path integral ground state approach for asymmetric top rotors with nuclear spin symmetry: Application to water chains

A path integral ground state (PIGS) approach for the simulation of asymmetric top rotors is presented. The method is based on Monte Carlo sampling of angular degrees of freedom. A symmetry-adapted rotational density matrix is used to account for nuclear spin statistics. To illustrate the method, ground-state properties of collections of para-water molecules confined to a one-dimensional lattice are computed. Those include energetic and structural observables. An advantage of the PIGS method is that expectation values can be obtained directly since the square of the wavefunction is sampled during a simulation. To benchmark the method, ground state energies and orientational distributions are computed using exact diagonalization for a single para-water molecule in an external field using a finite basis of symmetric top eigenfunctions. Benchmark results are also provided for N = 2 para-water molecules pinned to lattice sites at various distances to sample the crossover from hydrogen bonding to the dipole-dipole interaction regime. Excellent agreement between the PIGS results and the finite basis set calculations is observed. A thorough analysis of the convergence in terms of the imaginary time propagation length and systematic Trotter error is performed. The PIGS approach is then applied to a chain of N = 11 water molecules, and an equation of state is constructed in terms of the intermolecular separation. Ordering effects are also studied, and a transition between hydrogen bonding to dipole-dipole alignment is observed. The method is scalable and can also be applied in higher dimensions.

The Loss of Size Sensitivity in para-Hydrogen Clusters Due to the Strong Quantum Delocalization

para-Hydrogen (pH₂)_N clusters have been the focus of numerous computational studies. Originally motivated by the possibility of observing superfluidity, these studies also revealed rich and complex structural properties of (pH₂)_N. However, their structural analysis was typically limited to attempts to identify "magic number clusters" by computing their ground state energies E_N and the chemical potential μ_N = E_N-E_N-1 as a function of N. This was followed by structural analysis based on an ill-defined radial density profile. Surprisingly, however, there were remarkable discrepancies between the results reported in the literature for cluster sizes beyond approximately N = 25, and this ambiguity remained unsettled until now. In the present paper, we apply the diffusion Monte Carlo method to resolve inconsistencies in cluster sizes within the range (N = 24-28). Here, we try to avoid speculations based on the highly demanding energy calculations whose numerical accuracy harbors ambiguity. Instead, we focus on the direct and unambiguous structural analysis of the ground state wavefunctions, which supports the conclusion that the clusters are structurally the same in the size range considered. That is, there are no magic number clusters at least in the range N = 24-28, contrary to what some of the previous publications have suggested. This lack of size sensitivity of para-hydrogen clusters is a direct consequence of the strong quantum delocalization in these systems.