1
|
Zane AR, Curotto E. Electrolyte clusters as hydrogen sponges: diffusion Monte Carlo simulations. Phys Chem Chem Phys 2022; 24:26094-26101. [PMID: 36268906 DOI: 10.1039/d2cp03658d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
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.
Collapse
Affiliation(s)
- A R Zane
- Department of Chemistry & Physics Arcadia University, 450 S. Easton Rd., Glenside, PA 19038, USA.
| | - E Curotto
- Department of Chemistry & Physics Arcadia University, 450 S. Easton Rd., Glenside, PA 19038, USA.
| |
Collapse
|
2
|
Curotto E. Ion-Stockmayer clusters: Minima, classical thermodynamics, and variational ground state estimates of Li+(CH3NO2)n(n= 1–20). J Chem Phys 2015; 143:214301. [DOI: 10.1063/1.4936587] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- E. Curotto
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| |
Collapse
|
3
|
Wolf S, Curotto E. Ring polymer dynamics for rigid tops with an improved integrator. J Chem Phys 2014; 141:024116. [DOI: 10.1063/1.4887460] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
4
|
Abstract
We formulate an extension of the ring polymer dynamics approach to curved spaces using stereographic projection coordinates. We test the theory by simulating the particle in a ring, T(1), mapped by a stereographic projection using three potentials. Two of these are quadratic, and one is a nonconfining sinusoidal model. We propose a new class of algorithms for the integration of the ring polymer Hamilton equations in curved spaces. These are designed to improve the energy conservation of symplectic integrators based on the split operator approach. For manifolds, the position-position autocorrelation function can be formulated in numerous ways. We find that the position-position autocorrelation function computed from configurations in the Euclidean space R(2) that contains T(1) as a submanifold has the best statistical properties. The agreement with exact results obtained with vector space methods is excellent for all three potentials, for all values of time in the interval simulated, and for a relatively broad range of temperatures.
Collapse
Affiliation(s)
- S Wolf
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | | |
Collapse
|
5
|
Curotto E, Mella M. Quantum Monte Carlo simulations of selected ammonia clusters (n = 2–5): Isotope effects on the ground state of typical hydrogen bonded systems. J Chem Phys 2010; 133:214301. [DOI: 10.1063/1.3506027] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
|
6
|
Asare E, Musah AR, Curotto E, Freeman DL, Doll JD. The thermodynamic and ground state properties of the TIP4P water octamer. J Chem Phys 2009; 131:184508. [DOI: 10.1063/1.3259047] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
7
|
Lubombo C, Curotto E, Janeiro Barral PE, Mella M. Thermodynamic properties of ammonia clusters (NH3)n n=2–11: Comparing classical and quantum simulation results for hydrogen bonded species. J Chem Phys 2009; 131:034312. [DOI: 10.1063/1.3159398] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
8
|
Curotto E, Freeman DL, Doll JD. A stereographic projection path integral study of the coupling between the orientation and the bending degrees of freedom of water. J Chem Phys 2008; 128:204107. [DOI: 10.1063/1.2925681] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Luan T, Curotto E, Mella M. Importance sampling for quantum Monte Carlo in manifolds: Addressing the time scale problem in simulations of molecular aggregates. J Chem Phys 2008; 128:164102. [DOI: 10.1063/1.2898539] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Avilés MW, McCandless ML, Curotto E. Stereographic projection path integral simulations of (HCl)n clusters (n=2–5): Evidence of quantum induced melting in small hydrogen bonded networks. J Chem Phys 2008; 128:124517. [DOI: 10.1063/1.2837802] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
Xie D, Ran H, Zhou Y. Potential energy surfaces and predicted infrared spectra for van der Waals complexes: dependence on one intramolecular vibrational coordinate. INT REV PHYS CHEM 2007. [DOI: 10.1080/01442350701437926] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
12
|
Avilés MW, Curotto E. Stereographic projection diffusion monte carlo (SPDMC) algorithms for molecular condensed matter. J Phys Chem A 2007; 111:2610-8. [PMID: 17388350 DOI: 10.1021/jp066827i] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We develop and test three algorithms for diffusion Monte Carlo simulations in non-Euclidean manifolds. The methods are based on the construction of the "velocity" distribution by rejection techniques and are capable of functioning in a broad class of non-Euclidean spaces generated by holonomic constraints. The formulation of the propagator for non-Euclidean manifolds avoids the use of Lagrange multipliers; it is derived instead from the Feynman quantization in manifolds proposed by DeWitt. The manifolds are mapped onto Rd by using stereographic projection coordinates. Numerical tests are conducted for the particle in a ring of unit radius subjected to a sinusoidal potential, for the electron in the field of an infinitely massive proton, and for a water molecule modeled as an asymmetric top subjected to an external field.
Collapse
Affiliation(s)
- Michael W Avilés
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | | |
Collapse
|
13
|
Langley SF, Curotto E, Freeman DL, Doll JD. Rigid quantum Monte Carlo simulations of condensed molecular matter: Water clusters in the n=2→8 range. J Chem Phys 2007; 126:084506. [PMID: 17343457 DOI: 10.1063/1.2484229] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The numerical advantage of quantum Monte Carlo simulations of rigid bodies relative to the flexible simulations is investigated for some simple systems. The results show that if high frequency modes in molecular condensed matter are predominantly in the ground state, the convergence of path integral simulations becomes nonuniform. Rigid body quantum parallel tempering simulations are necessary to accurately capture thermodynamic phenomena in the temperature range where the dynamics are influenced by intermolecular degrees of freedom; the stereographic projection path integral adapted for quantum simulations of asymmetric tops is a significantly more efficient strategy compared with Cartesian coordinate simulations for molecular condensed matter under these conditions. The reweighted random series approach for stereographic path integral Monte Carlo is refined and implemented for the quantum simulation of water clusters treated as an assembly of rigid asymmetric tops.
Collapse
Affiliation(s)
- Stephen F Langley
- Department of Chemistry, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | | | | | | |
Collapse
|
14
|
Avilés MW, Gray PT, Curotto E. Stereographic projection path-integral simulations of (HF)n clusters. J Chem Phys 2006; 124:174305. [PMID: 16689569 DOI: 10.1063/1.2192773] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We perform several quantum canonical ensemble simulations of (HF)(n) clusters. The HF stretches are rigid, and the stereographic projection path-integral method is employed for the simulation in the resulting curved configuration space. We make use of the reweighted random series techniques to accelerate the convergence of the path-integral simulation with respect to the number of path coefficients. We develop and test estimators for the total energy and heat capacity based on a finite difference approach for non-Euclidean spaces. The quantum effects at temperatures below 400 K are substantial for all sizes. We observe interesting thermodynamic behaviors in the quantum simulations of the octamer and the heptamer.
Collapse
Affiliation(s)
- Michael W Avilés
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038, USA
| | | | | |
Collapse
|
15
|
Pav DM, Curotto E. Parallel tempering simulations of the 13-center Lennard-Jones dipole-dipole cluster (muD=0-->0.5 a.u.). J Chem Phys 2005; 123:144301. [PMID: 16238385 DOI: 10.1063/1.2049279] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate the thermodynamic behavior of the thirteen center uniform Lennard-Jones dipole-dipole cluster [(LJDD)(13)] for a wide range of dipole moment strengths. We find a relatively wide range of potential parameters where solid-solid coexistence manifests itself. Using structural characterization methods we determine the shape of the few isomers that contribute to the solid-solid coexistence region. The thermal distributions of the size of the net dipole moment are broad even at the coldest temperatures of the simulation where the (LJDD)(13) cluster is solid.
Collapse
Affiliation(s)
- Dawn M Pav
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038, USA
| | | |
Collapse
|
16
|
Curotto E. A reweighted random series method for stereographic projection path integrals. J Chem Phys 2005; 123:134102. [PMID: 16223270 DOI: 10.1063/1.2036970] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
A set of general reweighted random series methods for metric affine spaces is developed. The extension of the theorems to metric affine spaces demands the introduction of a configuration-independent reference metric tensor; this geometric object is used to treat the path expansion coefficients beyond the core path, in both the partial averaging and the reweighted random series approach. Numerical tests are conducted by simulating a particle in a ring. The reweighted random series results show better convergence properties and better statistical quality at a fraction of the cost compared with the related partial averaging simulation.
Collapse
Affiliation(s)
- E Curotto
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038, USA.
| |
Collapse
|
17
|
Avilés MW, Curotto E. Partial averaging and the centroid virial estimator for stereographic projection path-integral simulations in curved spaces. J Chem Phys 2005; 122:164109. [PMID: 15945674 DOI: 10.1063/1.1884109] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We develop and test three different partial averaging theories for the stereographic projection path integral in curved spaces. Additionally, we adapt and test the centroid virial estimator for the kinetic energy in curved spaces. We tested both a confining as well as a nonconfining potential for the particle in a ring. All three partial averaging theories are demonstrated to converge linearly in the asymptotic region with k(-2)max, the number of Fourier coefficients. We use an error estimator to determine the optimal parameters for an extrapolation to infinite kmax. We verify that the centroid virial estimator (derived from the primitive DeWitt path-integral formula) converges to the kinetic energy for both confining and nonconfining potentials.
Collapse
Affiliation(s)
- Michael W Avilés
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038, USA
| | | |
Collapse
|
18
|
Oppenheimer CA, Curotto E. Parameter space minimization methods: Applications to Lennard-Jones–dipole-dipole clusters. J Chem Phys 2004; 121:6226-39. [PMID: 15446915 DOI: 10.1063/1.1786916] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The morphology of the uniform Lennard-Jones-dipole-dipole cluster with 13 centers (LJDD)13 is investigated over a relatively wide range of values of the dipole moment. We introduce and compare several necessary modifications of the basin-hopping algorithm for global optimization to improve its efficiency. We develop a general algorithm for T=0 Brownian dynamics in curved spaces, and a graph theoretical approach necessary for the elimination of dissociated states. We find that the (LJDD)13 cluster has icosahedral symmetry for small to moderate values of the dipole moment. As the dipole moment increases, however, its morphology shifts to an hexagonal antiprism, and eventually to a ring.
Collapse
Affiliation(s)
- Craig A Oppenheimer
- Department of Chemistry and Physics, Arcadia University, Glenside, Pennsylvania 19038-3295, USA
| | | |
Collapse
|