Quantum-induced solid-solid transitions and melting in the Lennard-Jones LJ38 cluster

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.

Exact bound rovibrational spectra of the neon tetramer

Exact quantum dynamics calculations are performed for the bound rovibrational states of the neon tetramer (Ne₄) in its ground electronic state, using pair-wise Lennard-Jones potentials and the ScalIT suite of parallel codes. The vibrational states separate into a low-lying group mostly localized to a single potential well and a higher-energy delocalized group lying above the isomerization threshold-with such a structure serving as a testament to the "intermediate" quantum nature of the Ne₄ system. To accurately and efficiently represent both groups of states, the phase-space optimized discrete variable representation (PSO-DVR) approach was used, as implemented in the ScalIT code. The resultant 1D PSO effective potentials also shed significant light on the quantum dynamics of the system. All vibrational states were computed well up into the isomerization band and labeled up to the classical isomerization threshold-defined as the addition of the classical energy of a single bond, ϵ = 24.7 cm^-1, to the quantum ground state energy. Rovibrational energy levels for all total angular momentum values in the range J = 1-5 were also computed, treating all Coriolis coupling exactly.

Magic numbers, quantum delocalization, and orientational disordering in anionic hydrogen and deuterium clusters

The Diffusion Monte Carlo (DMC) method was applied to anionic hydrogen clusters H^-(H₂)_n (n = 1-16, 32) and their deuterated analogs using a polarizable all-atom potential energy surface (PES) developed by Calvo and Yurtsever. For the hydrogen clusters, the binding energy ΔE_n appears to be a smooth function of the cluster size n, thus contradicting the previous claim that n = 12 is a "magic number" cluster. The structures of the low energy minima of the PES for these clusters belong to the icosahedral motif with the H₂ molecules aligned toward the central H^- ion. However, their ground state wavefunctions are highly delocalized and resemble neither the structures of the global nor local minima. Moreover, the strong nuclear quantum effects result in a nearly complete orientational disordering of the H₂ molecules. For the deuterium clusters, the ground state wavefunctions are localized and the D₂ molecules are aligned toward the central D^- ion. However, their structures are still characterized as disordered and, as such, do not display size sensitivity. In addition, DMC simulations were performed on the mixed H^-(H₂)_n(D₂)_p clusters with (n, p) = (6, 6) and (16, 16). Again, in contradiction to the previous claim, we found that the "more quantum" H₂ molecules prefer to reside farther from the central H^- ion than the D₂ molecules.