An official website of the United States government
DL_POLY Quantum 2.1 is introduced here as a highly modular, sustainable, and scalable general-purpose molecular dynamics (MD) simulation software for large-scale long-time MD simulations of condensed phase and interfacial systems with the essential nuclear quantum effects (NQEs) included. The new release improves upon version 2.0 through the introduction of several emerging real-time path integral (PI) methods, including fast centroid molecular dynamics (f-CMD) and fast quasi-CMD (f-QCMD) methods, as well as our recently introduced hybrid CMD (h-CMD) method for the accurate and efficient simulation of vibrational infrared spectra. Several test cases, including liquid bulk water at 300 K and ice Ih at 150 K, are used to showcase the performance of different implemented PI methods in simulating the infrared spectra at both ambient conditions and low temperatures where NQEs become more apparent. Additionally, using different salt-in-water (i.e., dilute) and water-in-salt (i.e., concentrated) lithium bis(trifluoromethanesulfonyl)imide (Li-TFSI) aqueous electrolyte solutions, we demonstrate the applicability of our recently introduced h-CMD method implemented in DL_POLY Quantum 2.1 for the large scale simulation of infrared (IR) spectra of complex heterogeneous systems. We show that h-CMD can overcome the curvature problem of CMD and the artificial broadening of T-RPMD for the accurate simulation of the vibrational spectra of complex, heterogeneous systems with NQEs included.
more »
« less
h-CMD: An efficient hybrid fast centroid and quasi-centroid molecular dynamics method for the simulation of vibrational spectra
Developing efficient path integral (PI) methods for atomistic simulations of vibrational spectra in heterogeneous condensed phases and interfaces has long been a challenging task. Here, we present the h-CMD method, short for hybrid centroid molecular dynamics, which combines the recently introduced fast quasi-CMD (f-QCMD) method with fast CMD (f-CMD). In this scheme, molecules that are believed to suffer more seriously from the curvature problem of CMD, e.g., water, are treated with f-QCMD, while the rest, e.g., solid surfaces, are treated with f-CMD. To test the accuracy of the newly introduced scheme, the infrared spectra of the interfacial D2O confined in the archetypal ZIF-90 framework are simulated using h-CMD compared to a variety of other PI methods, including thermostatted ring-polymer molecular dynamics (T-RPMD) and partially adiabatic CMD as well as f-CMD and experiment as reference. Comparisons are also made with classical MD, where nuclear quantum effects are neglected entirely. Our detailed comparisons at different temperatures of 250–600 K show that h-CMD produces O–D stretches that are in close agreement with the experiment, correcting the known curvature problem and redshifting of the stretch peaks of CMD. h-CMD also corrects the known issues associated with too artificially dampened and broadened spectra of T-RPMD, which leads to missing the characteristic doublet feature of the interfacial confined water, rendering it unsuitable for these systems. The new h-CMD method broadens the applicability of f-QCMD to heterogeneous condensed phases and interfaces, where defining curvilinear coordinates for the entire system is not feasible.
more »
« less
- Award ID(s):
- 2302618
- PAR ID:
- 10564558
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 162
- Issue:
- 1
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The vibrational spectra of condensed and gas-phase systems are influenced by thequantum-mechanical behavior of light nuclei. Full-dimensional simulations of approximate quantum dynamics are possible thanks to the imaginary time path-integral (PI) formulation of quantum statistical mechanics, albeit at a high computational cost which increases sharply with decreasing temperature. By leveraging advances in machine-learned coarse-graining, we develop a PI method with the reduced computational cost of a classical simulation. We also propose a simple temperature elevation scheme to significantly attenuate the artifacts of standard PI approaches as well as eliminate the unfavorable temperature scaling of the computational cost. We illustrate the approach, by calculating vibrational spectra using standard models of water molecules and bulk water, demonstrating significant computational savings and dramatically improved accuracy compared to more expensive reference approaches. Our simple, efficient, and accurate method has prospects for routine calculations of vibrational spectra for a wide range of molecular systems - with an explicit treatment of the quantum nature of nuclei.more » « less
-
Abstract For a free group F r F_{r} of finite rank r ≥ 2 r\geq 2 and a non-trivial element w ∈ F r w\in F_{r} , the primitivity rank π ( w ) \pi(w) is the smallest rank of a subgroup H ≤ F r H\leq F_{r} such that w ∈ H w\in H and 𝑤 is not primitive in 𝐻 (if no such 𝐻 exists, one puts π ( w ) = ∞ \pi(w)=\infty ).The set of all subgroups of F r F_{r} of rank π ( w ) \pi(w) containing 𝑤 as a non-primitive element is denoted by Crit ( w ) \operatorname{Crit}(w) .These notions were introduced by Puder (2014).We prove that there exists an exponentially generic subset V ⊆ F r V\subseteq F_{r} such that, for every w ∈ V w\in V , we have π ( w ) = r \pi(w)=r and Crit ( w ) = { F r } \operatorname{Crit}(w)=\{F_{r}\} .more » « less
-
Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium–nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium–nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.more » « less
-
F or c e d at a f or a fl a p pi n g f oil e n er g y h ar v e st er wit h a cti v e l e a di n g e d g e m oti o n o p er ati n g i n t h e l o w r e d u c e d fr e q u e n c y r a n g e i s c oll e ct e d t o d et er mi n e h o w l e a di n g e d g e m oti o n aff e ct s e n er g y h ar v e sti n g p erf or m a n c e. T h e f oil pi v ot s a b o ut t h e mi dc h or d a n d o p er at e s i n t h e l o w r e d u c e d fr e q u e n c y r a n g e of 𝑓𝑓 𝑓𝑓 / 𝑈𝑈 ∞ = 0. 0 6 , 0. 0 8, a n d 0. 1 0 wit h 𝑅𝑅 𝑅𝑅 = 2 0 ,0 0 0 − 3 0 ,0 0 0 , wit h a pit c hi n g a m plit u d e of 𝜃𝜃 0 = 7 0 ∘ , a n d a h e a vi n g a m plit u d e of ℎ 0 = 0. 5 𝑓𝑓 . It i s f o u n d t h at l e a di n g e d g e m oti o n s t h at r e d u c e t h e eff e cti v e a n gl e of att a c k e arl y t h e str o k e w or k t o b ot h i n cr e a s e t h e lift f or c e s a s w ell a s s hift t h e p e a k lift f or c e l at er i n t h e fl a p pi n g str o k e. L e a di n g e d g e m oti o n s i n w hi c h t h e eff e cti v e a n gl e of att a c k i s i n cr e a s e d e arl y i n t h e str o k e s h o w d e cr e a s e d p erf or m a n c e. I n a d diti o n a di s cr et e v ort e x m o d el wit h v ort e x s h e d di n g at t h e l e a di n g e d g e i s i m pl e m e nt f or t h e m oti o n s st u di e d; it i s f o u n d t h at t h e m e c h a ni s m f or s h e d di n g at t h e l e a di n g e d g e i s n ot a d e q u at e f or t hi s p ar a m et er r a n g e a n d t h e m o d el c o n si st e ntl y o v er pr e di ct s t h e a er o d y n a mi c f or c e s.more » « less
