Anharmonic Infrared Spectroscopy through the Fourier Transform of Time Correlation Function Formalism in ONETEP.

Density functional theory molecular dynamics (DFT-MD) provides an efficient framework for accurately computing several types of spectra. The major benefit of DFT-MD approaches lies in the ability to naturally take into account the effects of temperature and anharmonicity, without having to introduce any ad hoc or a posteriori corrections. Consequently, computational spectroscopy based on DFT-MD approaches plays a pivotal role in the understanding and assignment of experimental peaks and bands at finite temperature, particularly in the case of floppy molecules. Linear-scaling DFT methods can be used to study large and complex systems, such as peptides, DNA strands, amorphous solids, and molecules in solution. Here, we present the implementation of DFT-MD IR spectroscopy in the ONETEP linear-scaling code. In addition, two methods for partitioning the dipole moment within the ONETEP framework are presented. Dipole moment partitioning allows us to compute spectra of molecules in solution, which fully include the effects of the solvent, while at the same time removing the solvent contribution from the spectra.

[1]  D. Vanderbilt,et al.  Exponential decay properties of Wannier functions and related quantities. , 2001, Physical review letters.

[2]  M. Alderton,et al.  Distributed multipole analysis , 2006 .

[3]  Chris-Kriton Skylaris,et al.  Introducing ONETEP: linear-scaling density functional simulations on parallel computers. , 2005, The Journal of chemical physics.

[4]  J. Dziedzic,et al.  Linear-scaling calculation of Hartree-Fock exchange energy with non-orthogonal generalised Wannier functions. , 2013, The Journal of chemical physics.

[5]  M. Gaigeot.  Infrared spectroscopy of the alanine dipeptide analog in liquid water with DFT-MD. Direct evidence for P(II)/beta conformations. , 2010, Physical chemistry chemical physics : PCCP.

[6]  M. Scheffler,et al.  Coupled cluster benchmarks of water monomers and dimers extracted from density-functional theory liquid water: the importance of monomer deformations. , 2009, The Journal of chemical physics.

[7]  Vanderbilt,et al.  Generalization of the density-matrix method to a nonorthogonal basis. , 1994, Physical review. B, Condensed matter.

[8]  Arash A. Mostofi,et al.  Nonorthogonal generalized Wannier function pseudopotential plane-wave method , 2002 .

[9]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[10]  Anthony J. Stone,et al.  Distributed multipole analysis, or how to describe a molecular charge distribution , 1981 .

[11]  D. Marx,et al.  Quantum corrections to classical time-correlation functions: hydrogen bonding and anharmonic floppy modes. , 2004, The Journal of chemical physics.

[12]  N. Marzari,et al.  Maximally-localized Wannier Functions: Theory and Applications , 2011, 1112.5411.

[13]  M. Robinson,et al.  Accurate ionic forces and geometry optimization in linear-scaling density-functional theory with local orbitals , 2011, 1103.5869.

[14]  Adam Wasserman,et al.  On the foundations of chemical reactivity theory. , 2007, The journal of physical chemistry. A.

[15]  Stefan Grimme,et al.  Semiempirical GGA‐type density functional constructed with a long‐range dispersion correction , 2006, J. Comput. Chem..

[16]  Nevena Todorova,et al.  Dimensionality of Carbon Nanomaterials Determines the Binding and Dynamics of Amyloidogenic Peptides: Multiscale Theoretical Simulations , 2013, PLoS Comput. Biol..

[17]  P. Tavan,et al.  Vibrational spectra from atomic fluctuations in dynamics simulations. I. Theory, limitations, and a sample application. , 2004, The Journal of chemical physics.

[18]  Alanine polypeptide structural fingerprints at room temperature: what can be gained from non-harmonic Car-Parrinello molecular dynamics simulations. , 2008, The journal of physical chemistry. A.

[19]  S. Corcelli,et al.  Molecular dynamics investigation of the vibrational spectroscopy of isolated water in an ionic liquid. , 2014, The journal of physical chemistry. B.

[20]  R. Zwanzig,et al.  Time-Correlation Functions and Transport Coefficients in Statistical Mechanics , 1965 .

[21]  Peter D. Haynes,et al.  Localised spherical-wave basis set for O(N) total-energy pseudopotential calculations , 1997 .

[22]  Chris-Kriton Skylaris,et al.  Including dispersion interactions in the ONETEP program for linear-scaling density functional theory calculations , 2009, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[23]  L. H. D. L. Peña.  On the zero temperature limit of the Kubo-transformed quantum time correlation function , 2014 .

[24]  Walter Kohn.  NEARSIGHTEDNESS OF ELECTRONIC MATTER , 2008 .

[25]  T. Keyes,et al.  Time correlation function and finite field approaches to the calculation of the fifth order Raman response in liquid xenon. , 2006, The Journal of chemical physics.

[26]  Gustavo E. Scuseria,et al.  Linear scaling conjugate gradient density matrix search as an alternative to diagonalization for first principles electronic structure calculations , 1997 .

[27]  Robert Rein,et al.  On Physical Properties and Interactions of Polyatomic Molecules: With Application to Molecular Recognition in Biology , 1973 .

[28]  Yuanqin Yu,et al.  New C-H Stretching Vibrational Spectral Features in the Raman Spectra of Gaseous and Liquid Ethanol † , 2007 .

[29]  M. Tuckerman.  Statistical Mechanics: Theory and Molecular Simulation , 2010 .

[30]  M. Gaigeot.  Theoretical spectroscopy of floppy peptides at room temperature. A DFTMD perspective: gas and aqueous phase. , 2010, Physical chemistry chemical physics : PCCP.

[31]  Arash A. Mostofi,et al.  Density kernel optimization in the ONETEP code , 2008 .

[32]  S R Elliott,et al.  Spatial decay of the single-particle density matrix in insulators: analytic results in two and three dimensions. , 2002, Physical review letters.

[33]  Barbara Kirchner,et al.  Computing vibrational spectra from ab initio molecular dynamics. , 2013, Physical chemistry chemical physics : PCCP.

[34]  Kohn,et al.  Density functional and density matrix method scaling linearly with the number of atoms. , 1996, Physical review letters.

[35]  Arash A. Mostofi,et al.  Minimal parameter implicit solvent model for ab initio electronic-structure calculations , 2011, 1112.1268.

[36]  S. Saito,et al.  Ultrafast energy relaxation and anisotropy decay of the librational motion in liquid water: A molecular dynamics study. , 2009, The Journal of chemical physics.

[37]  Sohrab Ismail-Beigi,et al.  LOCALITY OF THE DENSITY MATRIX IN METALS, SEMICONDUCTORS, AND INSULATORS , 1999 .

[38]  S. Goedecker.  DECAY PROPERTIES OF THE FINITE-TEMPERATURE DENSITY MATRIX IN METALS , 1998 .

[39]  Generalized Wannier functions: A comparison of molecular electric dipole polarizabilities , 2012, 1203.4371.

[40]  T. Fox,et al.  Large‐scale DFT calculations in implicit solvent—A case study on the T4 lysozyme L99A/M102Q protein , 2013 .

[41]  Jonathan W Essex,et al.  Free energies of binding from large-scale first-principles quantum mechanical calculations: application to ligand hydration energies. , 2013, The journal of physical chemistry. B.

[42]  Chris-Kriton Skylaris,et al.  Hybrid MPI-OpenMP Parallelism in the ONETEP Linear-Scaling Electronic Structure Code: Application to the Delamination of Cellulose Nanofibrils. , 2014, Journal of chemical theory and computation.

[43]  Renato Lemus,et al.  Vibrational excitations in H2O in the framework of a local model , 2004 .

[44]  T. Thonhauser,et al.  van der Waals density functional study of energetic, structural, and vibrational properties of small water clusters and ice I h , 2011, 1106.4830.

[45]  Minhaeng Cho,et al.  Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: comparisons with IR and vibrational circular dichroism spectra. , 2008, The Journal of chemical physics.

[46]  C. Anfinsen.  Principles that govern the folding of protein chains. , 1973, Science.

[47]  B. Brooks,et al.  A 500 ps molecular dynamics simulation study of interleukin-1 beta in water. Correlation with nuclear magnetic resonance spectroscopy and crystallography. , 1992, Journal of molecular biology.

[48]  Anthony J Stone,et al.  Distributed Multipole Analysis:  Stability for Large Basis Sets. , 2005, Journal of chemical theory and computation.

[49]  R. Mcweeny.  Some Recent Advances in Density Matrix Theory , 1960 .

[50]  J. D. Cloizeaux.  Analytical Properties of n-Dimensional Energy Bands and Wannier Functions , 1964 .