Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDiscretized Wiener-Khinchin theorem for Fourier-Laplace transformation: application to molecular simulations
129
0
0.0
(
0
)
Authors
Akira Koyama
Ask ChatGPT about the research
No Arabic abstract
The Wiener-Khinchin theorem for the Fourier-Laplace transformation (WKT-FLT) provides a robust method to calculate numerically single-side Fourier transforms of arbitrary autocorrelation functions from molecular simulations. However, the existing WKT-FLT equation produces two artifacts in the output of the frequency-domain relaxation function. In addition, these artifacts are more apparent in the frequency-domain response function converted from the relaxation function. We find the sources of these artifacts that are associated with the discretization of the WKT-FLT equation. Taking these sources into account, we derive the new discretized WKT-FLT equations designated for both the frequency-domain relaxation and response functions with the artifacts removed. The use of the discretized WKT-FLT equations is illustrated by a flow chart of an on-the-fly algorithm. We also give application examples of the discretized WKT-FLT equations for computing dynamic structure factor and wave-vector-dependent dynamic susceptibility from molecular simulations.
rate research
Read More
The classical Wiener-Khinchin theorem (WKT), which can extract spectral information by classical interferometers through Fourier transform, is a fundamental theorem used in many disciplines. However, there is still need for a quantum version of WKT, which could connect correlated biphoton spectral information by quantum interferometers. Here, we extend the classical WKT to its quantum counterpart, i.e., extended WKT (e-WKT), which is based on two-photon quantum interferometry. According to the e-WKT, the difference-frequency distribution of the biphoton wavefunctions can be extracted by applying a Fourier transform on the time-domain Hong-Ou-Mandel interference (HOMI) patterns, while the sum-frequency distribution can be extracted by applying a Fourier transform on the time-domain NOON state interference (NOONI) patterns. We also experimentally verified the WKT and e-WKT in a Mach-Zehnder interference (MZI), a HOMI and a NOONI. This theorem can be directly applied to quantum spectroscopy, where the spectral correlation information of biphotons can be obtained from time-domain quantum interferences by Fourier transform. This may open a new pathway for the study of light-matter interaction at the single photon level.
Recent developments in path integral methodology have significantly reduced the computational expense of including quantum mechanical effects in the nuclear motion in ab initio molecular dynamics simulations. However, the implementation of these developments requires a considerable programming effort, which has hindered their adoption. Here we describe i-PI, an interface written in Python that has been designed to minimise the effort required to bring state-of-the-art path integral techniques to an electronic structure program. While it is best suited to first principles calculations and path integral molecular dynamics, i-PI can also be used to perform classical molecular dynamics simulations, and can just as easily be interfaced with an empirical forcefield code. To give just one example of the many potential applications of the interface, we use it in conjunction with the CP2K electronic structure package to showcase the importance of nuclear quantum effects in high pressure water.
Computer simulations of model systems are widely used to explore striking phenomena in promising applications spanning from physics, chemistry, biology, to materials science and engineering. The long range electrostatic interactions between charged particles constitute a prominent factor in determining structures and states of model systems. How to efficiently calculate electrostatic interactions in model systems subjected to partial or full periodic boundary conditions has been a grand challenging task. In the past decades, a large variety of computational schemes have been proposed, among which the Ewald summation method is the most reliable route to accurately deal with electrostatic interactions in model systems. In addition, extensive effort has been done to improve computational efficiency of the Ewald summation based methods. Representative examples are approaches based on cutoffs, reaction fields, multi-poles, multi-grids, and particle-mesh schemes. We sketched an ENUF method, an abbreviation for the Ewald summation method based on Non-Uniform fast Fourier transform technique, and have implemented this method in particle-based simulation packages to calculate electrostatic energies and forces at micro- and mesoscopic levels. Extensive computational studies of conformational properties of polyelectrolytes, dendrimer-membrane complexes, and ionic fluids demonstrated that the ENUF method and its derivatives conserve both energy and momentum to floating point accuracy, and exhibit a computational complexity of $mathcal{O}(Nlog N)$ with optimal physical parameters. These ENUF based methods are attractive alternatives in molecular simulations where high accuracy and efficiency of simulation methods are needed to accelerate calculations of electrostatic interactions at extended spatiotemporal scales.
A discretized version of the adiabatic theorem is described with the help of a rule relating a Hermitian operator to its expectation value and variance. The simple initial operator X with known ground state is transformed in a series of N small steps into a more complicated final operator Z with unknown ground state. Each operator along the discretised path in the space of Hermitian matrices is used to measure the state, initially the ground state of X. Measurements similar to the Zeno effect or Renningers negative measurements modify the state incrementally. This process eventually leads to an eigenstate combination of Z. In the limit of vanishing step size the state stays with overwhelming probability in the ground state of each of the N observables.
A fitting scheme is proposed to obtain effective potentials from Car-Parrinello molecular dynamics (CPMD) simulations. It is used to parameterize a new pair potential for silica. MD simulations with this new potential are done to determine structural and dynamic properties and to compare these properties to those obtained from CPMD and a MD simulation using the so-called BKS potential. The new potential reproduces accurately the liquid structure generated by the CPMD trajectories, the experimental activation energies for the self-diffusion constants and the experimental density of amorphous silica. Also lattice parameters and elastic constants of alpha-quartz are well-reproduced, showing the transferability of the new potential.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
