scholarly journals Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: Using density functional theory to accelerate wavefunction methods

2016 ◽  
Vol 144 (5) ◽  
pp. 054111 ◽  
Author(s):  
V. Kapil ◽  
J. VandeVondele ◽  
M. Ceriotti
2021 ◽  
Author(s):  
Xi Xu ◽  
Zehua Chen ◽  
Yang Yang

Nuclear quantum effects play a crucial role in many chemical and biological systems involving hydrogen atoms yet are difficult to include in practical molecular simulations. In this Letter, we combine our recently developed methods of constrained nuclear-electronic orbital density functional theory (cNEO-DFT) and constrained minimized energy surface molecular dynamics (CMES-MD) to create a new method for accurately and efficiently describing nuclear quantum effects in molecular simulations. Using this new method, dubbed cNEO-MD, the vibrational spectra of a set of small molecules are calculated and compared with those from conventional ab initio molecular dynamics (AIMD) as well as from experiments. With the same formal scaling, cNEO-MD greatly outperforms AIMD in describing the vibrational modes with significant hydrogen motion characters, demonstrating the promise of cNEO-MD for simulating chemical and biological systems with significant nuclear quantum effects.


2021 ◽  
Vol 22 (10) ◽  
pp. 5220
Author(s):  
Jarosław J. Panek ◽  
Joanna Zasada ◽  
Bartłomiej M. Szyja ◽  
Beata Kizior ◽  
Aneta Jezierska

The O-H...N and O-H...O hydrogen bonds were investigated in 10-hydroxybenzo[h]quinoline (HBQ) and benzo[h]quinoline-2-methylresorcinol complex in vacuo, solvent and crystalline phases. The chosen systems contain analogous donor and acceptor moieties but differently coupled (intra- versus intermolecularly). Car–Parrinello molecular dynamics (CPMD) was employed to shed light onto principle components of interactions responsible for the self-assembly. It was applied to study the dynamics of the hydrogen bonds and vibrational features as well as to provide initial geometries for incorporation of quantum effects and electronic structure studies. The vibrational features were revealed using Fourier transformation of the autocorrelation function of atomic velocity and by inclusion of nuclear quantum effects on the O-H stretching solving vibrational Schrödinger equation a posteriori. The potential of mean force (Pmf) was computed for the whole trajectory to derive the probability density distribution and for the O-H stretching mode from the proton vibrational eigenfunctions and eigenvalues incorporating statistical sampling and nuclear quantum effects. The electronic structure changes of the benzo[h]quinoline-2-methylresorcinol dimer and trimers were studied based on Constrained Density Functional Theory (CDFT) whereas the Electron Localization Function (ELF) method was applied for all systems. It was found that the bridged proton is localized on the donor side in both investigated systems in vacuo. The crystalline phase simulations indicated bridged proton-sharing and transfer events in HBQ. These effects are even more pronounced when nuclear quantization is taken into account, and the quantized Pmf allows the proton to sample the acceptor area more efficiently. The CDFT indicated the charge depletion at the bridged proton for the analyzed dimer and trimers in solvent. The ELF analysis showed the presence of the isolated proton (a signature of the strongest hydrogen bonds) only in some parts of the HBQ crystal simulation. The collected data underline the importance of the intramolecular coupling between the donor and acceptor moieties.


Author(s):  
Lijuan Meng ◽  
Jinlian Lu ◽  
Yujie Bai ◽  
Lili Liu ◽  
Tang Jingyi ◽  
...  

Understanding the fundamentals of chemical vapor deposition bilayer graphene growth is crucial for its synthesis. By employing density functional theory calculations and classical molecular dynamics simulations, we have investigated the...


Author(s):  
Alberto Rodríguez-Fernández ◽  
Laurent Bonnet ◽  
Pascal Larrégaray ◽  
Ricardo Díez Muiño

The dissociation process of hydrogen molecules on W(110) was studied using density functional theory and classical molecular dynamics.


Sign in / Sign up

Export Citation Format

Share Document