On the convergence of diffusion Monte Carlo in non-Euclidean spaces. II. Diffusion with sources and sinks

Electrolyte clusters as hydrogen sponges: diffusion Monte Carlo simulations

We carry out Diffusion Monte Carlo simulations of up to five hydrogen molecules aggregated with two Stockmayer clusters that solvate a single lithium ion. The first one contains six point dipole solvent particles with parameters tuned to emulate nitromethane. The second cluster is a relative large system investigated recently [G. DiEmma, S. Kalette, and E. Curotto, Chem. Phys. Lett. 2019, 725, 80-86]. In both cases we find that the aggregated hydrogen molecules perturb significantly the ground state of the host cluster and form a distorted tetrahedral cage around the Li⁺ ion. The fifth hydrogen molecule is absorbed by the larger Stockmayer cluster while remaining in the proximity of the solvated charge.

Diffusion Monte Carlo simulations of gas phase and adsorbed D2-(H2)n clusters

We have computed ground state energies and analyzed radial distributions for several gas phase and adsorbed D₂(H₂)_n and HD(H₂)_n clusters. An external model potential designed to mimic ionic adsorption sites inside porous materials is used [M. Mella and E. Curotto, J. Phys. Chem. A 121, 5005 (2017)]. The isotopic substitution lowers the ground state energies by the expected amount based on the mass differences when these are compared with the energies of the pure clusters in the gas phase. A similar impact is found for adsorbed aggregates. The dissociation energy of D₂ from the adsorbed clusters is always much higher than that of H₂ from both pure and doped aggregates. Radial distributions of D₂ and H₂ are compared for both the gas phase and adsorbed species. For the gas phase clusters, two types of hydrogen-hydrogen interactions are considered: one based on the assumption that rotations and translations are adiabatically decoupled and the other based on nonisotropic four-dimensional potential. In the gas phase clusters of sufficiently large size, we find the heavier isotopomer more likely to be near the center of mass. However, there is a considerable overlap among the radial distributions of the two species. For the adsorbed clusters, we invariably find the heavy isotope located closer to the attractive interaction source than H₂, and at the periphery of the aggregate, H₂ molecules being substantially excluded from the interaction with the source. This finding rationalizes the dissociation energy results. For D₂-(H₂)_n clusters with n≥12, such preference leads to the desorption of D₂ from the aggregate, a phenomenon driven by the minimization of the total energy that can be obtained by reducing the confinement of (H₂)₁₂. The same happens for (H₂)₁₃, indicating that such an effect may be quite general and impact on the absorption of quantum species inside porous materials.

Structural and dynamic inhomogeneities induced by curvature gradients in elliptic colloidal halos of paramagnetic particles

Paramagnetic colloidal particles distributed along an ellipse are used as a model system to study the effects of curvature gradients on the structure and dynamics of colloids in curved manifolds. Unlike what happens for circular and spherical systems, in the present case, the equilibrium one-particle distribution function displays inhomogeneities due to the changing curvature along the ellipse. The ensuing effects on the two-body correlations are also analyzed, leading to the observation of anisotropic and long-ranged effects. Another noticeable consequence is the slowing down of the self-diffusion of these particles, which for large eccentricities may induce metastable states; this is evaluated by means of the time-dependent self-distribution.

Quest for Inexpensive Hydrogen Isotopic Fractionation: Do We Need 2D Quantum Confining in Porous Materials or Are Rough Surfaces Enough? The Case of Ammonia Nanoclusters

We study the adsorption energetics and quantum properties of the molecular hydrogen isotopes H₂, D₂, and T₂ onto the surface of rigid ammonia nanoclusters with quantum simulations and accurate model potential energy surfaces (PES). A highly efficient diffusion Monte Carlo (DMC) algorithm for rigid rotors allowed us to accurately define zero-point adsorption energies for the three isotopes, as well as the degree of translational and rotational delocalization that each affords on the surface. From the data emerges that the quantum adsorption energy (E_ads) of T₂ can be up to twice the one of H₂ at 0 K, suggesting the possibility of exploiting some form of solid ammonia to selectivity separate hydrogen isotopes at low temperatures (≃20 K). This is discussed by focusing on the structural motif that may be more effective for the task. The analysis of the contributions to E_ads, however, surprisingly indicates that the average kinetic energy (E^kin) and rotation energy (E_rot^kin) of T₂ can also be, respectively, 2 times and 20 times higher than those of H₂; this finding markedly deviates from what is predicted for hydrogen molecules inside carbon nanotubes (CNT) or metallic-organic frameworks (MOF), where E^kin and E_rot^kin is higher for H₂ due to the unavoidable effects of confinement and hindrance to its rotational motion. The rationale for these differences is provided by the geometrical distributions for the rigid rotors, which reveal an increasingly stronger coupling between rotational and translational degrees of freedom upon increasing the isotopic mass. This effect has never been observed before on adsorbing surfaces (e.g., graphite) and is induced by a strongly anisotropic and anharmonic bowl-like potential experienced by the rotors.

Smart darting diffusion Monte Carlo: Applications to lithium ion-Stockmayer clusters

In a recent investigation [K. Roberts et al., J. Chem. Phys. 136, 074104 (2012)], we have shown that, for a sufficiently complex potential, the Diffusion Monte Carlo (DMC) random walk can become quasiergodic, and we have introduced smart darting-like moves to improve the sampling. In this article, we systematically characterize the bias that smart darting moves introduce in the estimate of the ground state energy of a bosonic system. We then test a simple approach to eliminate completely such bias from the results. The approach is applied for the determination of the ground state of lithium ion-n-dipoles clusters in the n = 8-20 range. For these, the smart darting diffusion Monte Carlo simulations find the same ground state energy and mixed-distribution as the traditional approach for n < 14. In larger systems we find that while the ground state energies agree quantitatively with or without smart darting moves, the mixed-distributions can be significantly different. Some evidence is offered to conclude that introducing smart darting-like moves in traditional DMC simulations may produce a more reliable ground state mixed-distribution.

On the convergence of diffusion Monte Carlo in non-Euclidean spaces. I. Free diffusion

We develop a set of diffusion Monte Carlo algorithms for general compactly supported Riemannian manifolds that converge weakly to second order with respect to the time step. The approaches are designed to work for cases that include non-orthogonal coordinate systems, nonuniform metric tensors, manifold boundaries, and multiply connected spaces. The methods do not require specially designed coordinate charts and can in principle work with atlases of charts. Several numerical tests for free diffusion in compactly supported Riemannian manifolds are carried out for spaces relevant to the chemical physics community. These include the circle, the 2-sphere, and the ellipsoid of inertia mapped with traditional angles. In all cases, we observe second order convergence, and in the case of the sphere, we gain insight into the function of the advection term that is generated by the curved nature of the space.