Journal
JOURNAL OF CHEMICAL PHYSICS
Volume 131, Issue 9, Pages -Publisher
AMER INST PHYSICS
DOI: 10.1063/1.3216520
Keywords
electric potential; polymers; potential energy functions; quantum theory
Funding
- European Union [MTKD-CT-029583]
- U. S. Office of Naval Research [N000140510460]
- U. K. Engineering and Physical Sciences Research Council [E01741X]
- E. U. [205066]
- Engineering and Physical Sciences Research Council [EP/E01741X/1] Funding Source: researchfish
- EPSRC [EP/E01741X/1] Funding Source: UKRI
Ask authors/readers for more resources
A quantum simulation of an imaginary time path integral typically requires around n times more computational effort than the corresponding classical simulation, where n is the number of ring polymer beads (or imaginary time slices) used in the calculation. It is however possible to improve on this estimate by decomposing the potential into a sum of slowly and rapidly varying contributions. If the slowly varying contribution changes only slightly over the length scale of the ring polymer, it can be evaluated on a contracted ring polymer with fewer than the full n beads (or equivalently on a lower order Fourier decomposition of the imaginary time path). Here we develop and test this idea for systems with polarizable force fields. The development consists of iterating the induction on the contracted ring polymer and applying an appropriate transformation to obtain the forces on the original n beads. In combination with a splitting of the Coulomb potential into its short- and long-range parts, this results in a method with little more than classical computational effort in the limit of large system size. The method is illustrated with simulations of liquid water at 300 K and hexagonal ice at 100 K using a recently developed flexible and polarizable Thole-type potential energy model.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available