A Hybrid Density Functional Theory/Molecular Mechanics Approach for Linear Response Properties in Heterogeneous Environments.

We introduce a density functional theory/molecular mechanical approach for computation of linear response properties of molecules in heterogeneous environments, such as metal surfaces or nanoparticles embedded in solvents. The heterogeneous embedding environment, consisting from metallic and nonmetallic parts, is described by combined force fields, where conventional force fields are used for the nonmetallic part and capacitance-polarization-based force fields are used for the metallic part. The presented approach enables studies of properties and spectra of systems embedded in or placed at arbitrary shaped metallic surfaces, clusters, or nanoparticles. The capability and performance of the proposed approach is illustrated by sample calculations of optical absorption spectra of thymidine absorbed on gold surfaces in an aqueous environment, where we study how different organizations of the gold surface and how the combined, nonadditive effect of the two environments is reflected in the optical absorption spectrum.

[1]  Luca Frediani,et al.  A combined quantum mechanics/molecular mechanics study of the one- and two-photon absorption in the green fluorescent protein. , 2012, Physical chemistry chemical physics : PCCP.

[2]  K. Mikkelsen,et al.  Charge transfer excitation energies in pyridine–silver complexes studied by a QM/MM method , 2009 .

[3]  Luca Frediani,et al.  The Dalton quantum chemistry program system , 2013, Wiley interdisciplinary reviews. Computational molecular science.

[4]  Nathalie Reuter,et al.  Frontier Bonds in QM/MM Methods: A Comparison of Different Approaches , 2000 .

[5]  V. Barone,et al.  Singlet excited-state behavior of uracil and thymine in aqueous solution: a combined experimental and computational study of 11 uracil derivatives. , 2006, Journal of the American Chemical Society.

[6]  G. Schatz,et al.  Discrete dipole approximation for calculating extinction and Raman intensities for small particles with arbitrary shapes , 1995 .

[7]  M. Field,et al.  A Generalized Hybrid Orbital (GHO) Method for the Treatment of Boundary Atoms in Combined QM/MM Calculations , 1998 .

[8]  Branislav Jansik,et al.  Calculations of static and dynamic polarizabilities of excited states by means of density functional theory. , 2004, The Journal of chemical physics.

[9]  Anders Wallqvist,et al.  A molecular dynamics study of polarizable water , 1989 .

[10]  Paweł Sałek,et al.  Density-functional theory of linear and nonlinear time-dependent molecular properties , 2002 .

[11]  V. Barone,et al.  Toward reliable density functional methods without adjustable parameters: The PBE0 model , 1999 .

[12]  C. Morari,et al.  Electronic Properties of DNA Nucleosides Adsorbed on a Au(100) Surface , 2012 .

[13]  F. Porta,et al.  Au–thymine, thymidine and thymidine 5′-monophosphate nanoparticles: chemical characterisation and cellular uptake studies into U87 cancer cells , 2012 .

[14]  Roberto Improta,et al.  Quantum mechanical computations and spectroscopy: from small rigid molecules in the gas phase to large flexible molecules in solution. , 2008, Accounts of chemical research.

[15]  Lasse Jensen,et al.  A discrete interaction model/quantum mechanical method to describe the interaction of metal nanoparticles and molecular absorption. , 2011, The Journal of chemical physics.

[16]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[17]  D. Neuhauser,et al.  Molecular nanopolaritonics: cross manipulation of near-field plasmons and molecules. I. Theory and application to junction control. , 2007, The Journal of chemical physics.

[18]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .

[19]  Paweł Sałek,et al.  Restricted density functional theory of linear time-dependent properties in open-shell molecules , 2003 .

[20]  Erik Lindahl,et al.  Association dynamics and linear and nonlinear optical properties of an N-acetylaladanamide probe in a POPC membrane. , 2013, Journal of the American Chemical Society.

[21]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[22]  Vincenzo Barone,et al.  Linear Response Theory and Electronic Transition Energies for a Fully Polarizable QM/Classical Hamiltonian. , 2012, Journal of chemical theory and computation.

[23]  K. Mikkelsen,et al.  An atomic capacitance–polarizability model for the calculation of molecular dipole moments and polarizabilities , 2001 .

[24]  K. Mikkelsen,et al.  Heterogeneous solvation: An ab initio approach , 2001 .

[25]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[26]  Jacopo Tomasi,et al.  Enhanced response properties of a chromophore physisorbed on a metal particle , 2001 .

[27]  G. Mie Beiträge zur Optik trüber Medien, speziell kolloidaler Metallösungen , 1908 .

[28]  L. Jensen,et al.  A discrete interaction model/quantum mechanical method for describing response properties of molecules adsorbed on metal nanoparticles. , 2010, The Journal of chemical physics.

[29]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[30]  L. Jensen,et al.  A discrete interaction model/quantum mechanical method for simulating surface-enhanced Raman spectroscopy. , 2012, The Journal of chemical physics.

[31]  E. Davidson,et al.  One- and two-electron integrals over cartesian gaussian functions , 1978 .

[32]  A. Mayer Formulation in terms of normalized propagators of a charge-dipole model enabling the calculation of the polarization properties of fullerenes and carbon nanotubes , 2007 .

[33]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[34]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[35]  Luca Frediani,et al.  Excitation energies in solution: the fully polarizable QM/MM/PCM method. , 2011, The journal of physical chemistry. B.

[36]  L. Jensen,et al.  Electrostatic Interaction Model for the Calculation of the Polarizability of Large Noble Metal Nanoclusters , 2008 .

[37]  Roland Lindh,et al.  Local properties of quantum chemical systems: the LoProp approach. , 2004, The Journal of chemical physics.

[38]  Lasse Jensen,et al.  Theoretical studies of plasmonics using electronic structure methods. , 2011, Chemical reviews.

[39]  D. Truhlar,et al.  QM/MM: what have we learned, where are we, and where do we go from here? , 2007 .

[40]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[41]  Jan H. Jensen,et al.  QM/MM Boundaries Across Covalent Bonds: A Frozen Localized Molecular Orbital-Based Approach for the Effective Fragment Potential Method , 2000 .

[42]  H. Ågren,et al.  Encapsulation Influence on EPR Parameters of Spin-Labels: 2,2,6,6-Tetramethyl-4-methoxypiperidine-1-oxyl in Cucurbit[8]uril. , 2012, Journal of chemical theory and computation.

[43]  Paweł Sałek,et al.  Density functional theory of nonlinear triplet response properties with applications to phosphorescence , 2003 .

[44]  Xie,et al.  Single molecule emission characteristics in near-field microscopy. , 1995, Physical review letters.

[45]  Tiffany R Walsh,et al.  GolP-CHARMM: First-Principles Based Force Fields for the Interaction of Proteins with Au(111) and Au(100). , 2013, Journal of chemical theory and computation.

[46]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[47]  Jacob Kongsted,et al.  Density functional self-consistent quantum mechanics/molecular mechanics theory for linear and nonlinear molecular properties: Applications to solvated water and formaldehyde. , 2007, The Journal of chemical physics.

[48]  Kurt V. Mikkelsen,et al.  The QM/MM approach for wavefunctions, energies and response functions within self-consistent field and coupled cluster theories , 2002 .

[49]  K. Prince,et al.  Photoemission Study of Thymidine Adsorbed on Au(111) and Cu(110) , 2010 .

[50]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[51]  Marcella Iannuzzi,et al.  Simulation of Adsorption Processes at Metallic Interfaces: An Image Charge Augmented QM/MM Approach. , 2013, Journal of chemical theory and computation.

[52]  Alessandro Laio,et al.  A Hamiltonian electrostatic coupling scheme for hybrid Car-Parrinello molecular dynamics simulations , 2002 .

[53]  G. Schatz,et al.  Many-body theory of surface-enhanced Raman scattering , 2008, 0809.3207.

[54]  D. Sánchez-Portal,et al.  The SIESTA method for ab initio order-N materials simulation , 2001, cond-mat/0111138.