Comparison of basis set effects and the performance of ab initio and DFT methods for probing equilibrium fluctuations

The electronic absorption and emission spectra of large molecules reflect the extent and timescale of electron‐vibration coupling and therefore the extent and timescale of relaxation/reorganization in response to a perturbation. In this paper, we present a comparison of the calculated absorption and emission spectra of NADH in liver alcohol dehydrogenase (LADH), using quantum mechanical/molecular mechanical methods, in which we vary the QM component. Specifically, we have looked at the influence of basis set (STO‐3G, 3‐21G*, 6‐31G*, CC‐pVDZ, and 6‐311G**), as well as the influence of applying the DFT TD‐B3LYP and ab initio TD‐HF and CIS methods to the calculation of absorption/emission spectra and the reorganization energy (Stokes shift). The ab initio TD‐HF and CIS methods reproduce the experimentally determined Stokes shift and spectral profiles to a high level of agreement, while the TD‐B3LYP method significantly underestimates the Stokes shift, by 45%. We comment on the origin of this problem and suggest that ab initio methods may be naturally more suited to predicting molecular behavior away from equilibrium geometries. © 2006 Wiley Periodicals, Inc. J Comput Chem 28: 478–490, 2007

[1]  I. Gould,et al.  Optical properties of solvated molecules calculated by a QMMM method Chlorophyll a and bacteriochlorophyll a , 1997 .

[2]  Dennis R. Salahub,et al.  Dynamic polarizabilities and excitation spectra from a molecular implementation of time‐dependent density‐functional response theory: N2 as a case study , 1996 .

[3]  Nicholas C. Handy,et al.  On the determination of excitation energies using density functional theory , 2000 .

[4]  M. Head‐Gordon,et al.  Time-Dependent Density Functional Study of the Electronic Excited States of Polycyclic Aromatic Hydrocarbon Radical Ions , 2003 .

[5]  D. Chandler,et al.  Introduction To Modern Statistical Mechanics , 1987 .

[6]  H. Eklund,et al.  3 Alcohol Dehydrogenases , 1975 .

[7]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[8]  M. E. Casida.  Time-Dependent Density Functional Response Theory for Molecules , 1995 .

[9]  N. Handy,et al.  Study of excited states of furan and pyrrole by time-dependent density functional theory , 2002 .

[10]  S. Mukamel.  Fluorescence and absorption of large anharmonic molecules - spectroscopy without eigenstates , 1985 .

[11]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[12]  R. Kubo.  GENERALIZED CUMULANT EXPANSION METHOD , 1962 .

[13]  Michael J. Frisch,et al.  Toward a systematic molecular orbital theory for excited states , 1992 .

[14]  B. Plapp,et al.  Alternative pathways and reactions of benzyl alcohol and benzaldehyde with horse liver alcohol dehydrogenase. , 1993, Biochemistry.

[15]  Mark Earl Casida,et al.  In Recent Advances in Density-Functional Methods , 1995 .

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

[17]  Carole Van Caillie,et al.  Geometric derivatives of density functional theory excitation energies using gradient-corrected functionals , 2000 .

[18]  J. Pople,et al.  Self-consistent molecular orbital methods. 21. Small split-valence basis sets for first-row elements , 2002 .

[19]  D. Klug,et al.  A Quantum Mechanical/Molecular Mechanical Approach to Relaxation Dynamics: Calculation of the Optical Properties of Solvated Bacteriochlorophyll-a , 1999 .

[20]  Luis Serrano-Andrés,et al.  Does density functional theory contribute to the understanding of excited states of unsaturated organic compounds , 1999 .

[21]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[22]  EDWIN C. Webb.  The Enzymes , 1961, Nature.

[23]  J. B. Collins,et al.  Self‐consistent molecular orbital methods. XVII. Geometries and binding energies of second‐row molecules. A comparison of three basis sets , 1976 .

[24]  G. Scuseria,et al.  An efficient implementation of time-dependent density-functional theory for the calculation of excitation energies of large molecules , 1998 .

[25]  Alexander D. MacKerell,et al.  A molecular mechanics force field for NAD+ NADH, and the pyrophosphate groups of nucleotides , 1997, J. Comput. Chem..

[26]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[27]  Roger D. Amos,et al.  Geometric derivatives of excitation energies using SCF and DFT , 1999 .

[28]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[29]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[30]  S. Mukamel.  Principles of Nonlinear Optical Spectroscopy , 1995 .

[31]  D. Chong.  Recent Advances in Density Functional Methods Part III , 2002 .

[32]  R. Ahlrichs,et al.  Treatment of electronic excitations within the adiabatic approximation of time dependent density functional theory , 1996 .

[33]  David R. Klug,et al.  Large and Fast Relaxations inside a Protein: Calculation and Measurement of Reorganization Energies in Alcohol Dehydrogenase , 2002 .

[34]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. I. Use of Gaussian Expansions of Slater‐Type Atomic Orbitals , 1969 .

[35]  S. Ramaswamy,et al.  Binding of formamides to liver alcohol dehydrogenase. , 1997, Biochemistry.

[36]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[37]  M. Petersilka,et al.  Excitation energies from time-dependent density-functional theory. , 1996 .

[38]  Dennis R. Salahub,et al.  Molecular excitation energies to high-lying bound states from time-dependent density-functional response theory: Characterization and correction of the time-dependent local density approximation ionization threshold , 1998 .

[39]  R. Kubo.  The fluctuation-dissipation theorem , 1966 .