Abstract
In path integral molecular dynamics (PIMD) simulations, additional vibrations appear beyond the physical vibrations. In harmonic approximation, the frequencies of these internal modes can be determined from the physical frequencies. We show that this formal effect of the path integral simulations on the vibrations can be considered as a convolution if we use the square of the frequency as an independent variable. This convolution can be represented as a matrix multiplication. The potential of the formalism is demonstrated in two applications. We present an alternative method to determine the power spectrum of thermostats used in PIMD simulations. We also show that in simple anharmonic model systems, the physical frequencies can be obtained from ring polymer molecular dynamics simulations by deconvolution.
Supplementary materials
Title
Supporting Information
Description
Analysis of the effect of the lambda parameter on the deconvolution.
Actions
Title
Input files
Description
CP2K input files for the PIMD/RPMD simulation of OH and water molecule.
Actions
Supplementary weblinks
Title
Wolfram code to solve the path integral weight function problem
Description
This code demonstrates how to calculate the weight functions for path integral simulations. It works for small number of beads.
Actions
View Title
Fortran code to solve the path integral weight function problem
Description
This code is more effective than the Wolfram code.
Actions
View Title
Wolfram code for the deconvolution
Description
Wolfram code for the deconvolution of the vibrational density of states determined from path integral simulations
Actions
View