Showing papers on "Path integral molecular dynamics published in 2004"

Quantum statistics and classical mechanics: Real time correlation functions from ring polymer molecular dynamics

Abstract: We propose an approximate method for calculating Kubo-transformed real-time correlation functions involving position-dependent operators, based on path integral (Parrinello-Rahman) molecular dynamics. The method gives the exact quantum mechanical correlation function at time zero, exactly satisfies the quantum mechanical detailed balance condition, and for correlation functions of the form C(Ax)(t) and C(xB)(t) it gives the exact result for a harmonic potential. It also works reasonably well at short times for more general potentials and correlation functions, as we illustrate with some example calculations. The method provides a consistent improvement over purely classical molecular dynamics that is most apparent in the low-temperature regime.

Theoretical study on photoexcitation dynamics of the K atom attached to helium clusters and the solvation structures of K*Hen exciplexes

Abstract: The dynamics of the K atom attached to the 300-atom helium cluster upon photoexcitation has been studied theoretically based on the quantum molecular dynamics method. With the use of a hybrid method between time-dependent quantum dynamics and semiclassical path integral centroid molecular dynamics, the calculation includes the electronic state mixing for the K(2P) states as well as the quantum motions of helium atoms. It was found from the calculated trajectories that the K atom mostly desorbs by itself from the surface of the helium cluster without depending on the types of initial electronic excited states. It was also found that the K*He exciplex can be formed only as a minor product, following the initial excitation to the Π-type electronic state alone. Helium solvation structures of the K*Hen exciplexes were also investigated by path integral molecular dynamics simulations. The result implies that stable exciplexes can be formed up to n = 6 on the lowest excited state until the first solvation shell around the 4p orbital of K* is completed. Meanwhile, it was found that only two helium atoms can be bound to the K(2P) atom on the second lowest excited state.

Path integral hybrid Monte Carlo algorithm for correlated Bose fluids.

Abstract: Path integral hybrid Monte Carlo (PIHMC) algorithm for strongly correlated Bose fluids has been developed. This is an extended version of our previous method [S. Miura and S. Okazaki, Chem. Phys. Lett. 308, 115 (1999)] applied to a model system consisting of noninteracting bosons. Our PIHMC method for the correlated Bose fluids is constituted of two trial moves to sample path-variables describing system coordinates along imaginary time and a permutation of particle labels giving a boundary condition with respect to imaginary time. The path-variables for a given permutation are generated by a hybrid Monte Carlo method based on path integral molecular dynamics techniques. Equations of motion for the path-variables are formulated on the basis of a collective coordinate representation of the path, staging variables, to enhance the sampling efficiency. The permutation sampling to satisfy Bose–Einstein statistics is performed using the multilevel Metropolis method developed by Ceperley and Pollock [Phys. Rev. Lett. 56, 351 (1986)]. Our PIHMC method has successfully been applied to liquid helium-4 at a state point where the system is in a superfluid phase. Parameters determining the sampling efficiency are optimized in such a way that correlation among successive PIHMC steps is minimized.

Theoretical study on isotope and temperature effect in hydronium ion using ab initio path integral simulation

Abstract: We have applied the ab initio path integral molecular dynamics simulation to study hydronium ion and its isotopes, which are the simplest systems for hydrated proton and deuteron. In this simulation, all the rotational and vibrational degrees of freedom are treated fully quantum mechanically, while the potential energies of the respective atomic configurations are calculated “on the fly” using ab initio quantum chemical approach. With the careful treatment of the ab initio electronic structure calculation by relevant choices in electron correlation level and basis set, this scheme is theoretically quite rigorous except for Born–Oppenheimer approximation. This accurate calculation allows a close insight into the structural shifts for the isotopes of hydronium ion by taking account of both quantum mechanical and thermal effects. In fact, the calculation is shown to be successful to quantitatively extract the geometrical isotope effect with respect to the Walden inversion. It is also shown that this leads to the isotope effect on the electronic structure as well as the thermochemical properties.

Quantum effect on the internal proton transfer and structural fluctuation in the H+ 5 cluster.

Abstract: The thermal equilibrium state of H+(5) is investigated by means of an ab initio path integral molecular dynamics (PIMD) method, in which degrees of freedom of both nuclei and electrons at finite temperature are quantized within the adiabatic approximation. The second-order Moller-Plesset force field has been employed for the present ab initio PIMD. At 5-200 K, H+(5) is shown to have the structure that the proton is surrounded by the two H(2) units without any exchange of an atom between the central proton and the H(2) unit. At 5 K, the quantum tunneling of the central proton occurs more easily when the distance between the two H(2) units is shortened. At the high temperature of 200 K, the central proton is more delocalized in space between the two H(2) units, with less correlation with the stretching of the distance between the two H(2) units. As for the rotation of the H(2) units around the C(2) axis of H+(5) , the dihedral angle distribution is homogeneous at all temperatures, suggesting that the two H(2) units freely rotate around the C(2) axis, while this quantum effect on the rotation of the H(2) units becomes more weakened with increasing temperature. The influence of the structural fluctuation of H+(5) on molecular orbital energies has been examined to conclude that the highest occupied molecular orbital-lowest unoccupied molecular orbital energy gap is largely reduced with the increase of temperature because of the spatial expansion of the whole cluster.

Ab initio centroid path integral molecular dynamics: Application to vibrational dynamics of diatomic molecular systems

Abstract: An ab initio centroid molecular dynamics (CMD) method is developed by combining the CMD method with the ab initio molecular orbital method. The ab initio CMD method is applied to vibrational dynamics of diatomic molecules, H2 and HF. For the H2 molecule, the temperature dependence of the peak frequency of the vibrational spectral density is investigated. The results are compared with those obtained by the ab initio classical molecular dynamics method and exact quantum mechanical treatment. It is shown that the vibrational frequency obtained from the ab initio CMD approaches the exact first excitation frequency as the temperature lowers. For the HF molecule, the position autocorrelation function is also analyzed in detail. The present CMD method is shown to well reproduce the exact quantum result for the information on the vibrational properties of the system.

Effect of temperature on the formation of electronic bound states in an expanded bcc hydrogenoid crystal: A restricted path-integral molecular dynamics simulation

Abstract: We have used the restricted path-integral molecular dynamics method to study the correlated electronic structure of a half-filled expanded three-dimensional hydrogenoid body-centered cubic lattice at finite temperatures. Starting from a paramagnetic metallic state with electron gas character, we find that bound electrons form upon expansion of the lattice. The bound electrons are spatially localized with their center for the motion of gyration located on ionic positions. The region of coexistence of bound and unbound states in the temperature-density plane is reminiscent of that associated with a first-order transition. At constant temperature, the number of bound electrons increases monotonously with decreasing density. The width of the region of coexistence narrows with increasing temperature.

Formation of bound states in expanded metal studied via path integral molecular dynamics

Abstract: The usefulness of the restricted path integral molecular dynamics method for the study of strongly correlated electrons is demonstrated by studying the formation of bound electronic states in a half-filled expanded three-dimensional hydrogenoid body-centred cubic lattice at finite temperature. Starting from a metallic state with one-component plasma character, we find that bound electrons form upon expansion of the lattice. The bound electrons are spatially localized with their centre for the motion of gyration located at ionic positions. The number of bound electrons increases monotonically with decreasing density.