PubMed Health⌕ Search

Biomedical subjects

Michael J Frisch

Publications and source records attributed to Michael J Frisch.

16 recordsLinked to original sources

Time-dependent density functional theory investigation of the absorption, fluorescence, and phosphorescence spectra of solvated coumarins.

Using time-dependent density functional theory (TD-DFT) and the polarizable continuum model, we have computed the electronic transitions of a large panel of coumarin dyes in their enol, keto, cationic, and anionic forms. Several processes have been studied: absorption, fluorescence, 0-0 phosphorescence, and triplet-triplet excitations. For each process, detailed comparison with experimental data has been carried out. Using the PBE06-31+G(d) scheme, it turns out that for a given electronic transition the experimental shifts resulting from the substitution of the coumarin core are nicely reproduced. Indeed, once a simple statistical correction is applied, the mean absolute errors on the absorption and fluorescence wavelengths are limited to 8 nm (0.09 eV) and 9 nm (0.07 eV), respectively. A valuable correlation between the experimental and theoretical phosphorescence auxochromic displacements has also been unravelled. The differences between the wavelengths of the various electronic processes of a given dye tend to be fairly predicted, especially for the fluorescence-phosphoresence shifts that are strongly overestimated by TD-DFT.

Journal Article↗

Efficient evaluation of short-range Hartree-Fock exchange in large molecules and periodic systems.

We present an efficient algorithm for the evaluation of short-range Hartree-Fock exchange energies and geometry gradients in Gaussian basis sets. Our method uses a hierarchy of screening levels to eliminate negligible two-electron integrals whose evaluation is the fundamental computational bottleneck of the procedure. By applying our screening technique to the Heyd-Scuseria-Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] short-range Coulomb hybrid density functional, we achieve a computational efficiency comparable with that of standard nonhybrid density functional calculations.

Journal Article↗

A restricted-open-shell complete-basis-set model chemistry.

A restricted-open-shell model chemistry based on the complete basis set-quadratic Becke3 (CBS-QB3) model is formulated and denoted ROCBS-QB3. As the name implies, this method uses spin-restricted wave functions, both for the direct calculations of the various components of the electronic energy and for extrapolating the correlation energy to the complete-basis-set limit. These modifications eliminate the need for empirical corrections that are incorporated in standard CBS-QB3 to compensate for spin contamination when spin-unrestricted wave functions are used. We employ an initial test set of 19 severely spin-contaminated species including doublet radicals and both singlet and triplet biradicals. The mean absolute deviation (MAD) from experiment for the new ROCBS-QB3 model (3.6+/-1.5 kJ mol(-1)) is slightly smaller than that of the standard unrestricted CBS-QB3 version (4.8+/-1.5 kJ mol(-1)) and substantially smaller than the MAD for the unrestricted CBS-QB3 before inclusion of the spin correction (16.1+/-1.5 kJ mol(-1)). However, when applied to calculate the heats of formation at 298 K for the moderately spin-contaminated radicals in the G2/97 test set, ROCBS-QB3 does not perform quite as well as the standard unrestricted CBS-QB3, with a MAD from experiment of 3.8+/-1.6 kJ mol(-1) (compared with 2.9+/-1.6 kJ mol(-1) for standard CBS-QB3). ROCBS-QB3 performs marginally better than standard CBS-QB3 for the G2/97 set of ionization energies with a MAD of 4.1+/-0.1 kJ mol(-1) (compared with 4.4+/-0.1 kJ mol(-1)) and electron affinities with a MAD of 3.9+/-0.2 kJ mol(-1) (compared with 4.3+/-0.2 kJ mol(-1)), but the differences in MAD values are comparable to the experimental uncertainties. Our overall conclusion is that ROCBS-QB3 eliminates the spin correction in standard CBS-QB3 with no loss in accuracy.

Journal Article↗

A state-specific polarizable continuum model time dependent density functional theory method for excited state calculations in solution.

An effective state specific (SS) model for the inclusion of solvent effects in time dependent density functional theory (TD-DFT) computations of excited electronic states has been developed and coded in the framework of the so-called polarizable continuum model (PCM). Different relaxation time regimes can be treated thus giving access to a number of different spectroscopic properties together with solvent relaxation energies of paramount relevance in electron transfer processes. SS and conventional linear response (LR) models have been compared for two benchmark systems (coumarin 153 and formaldehyde in different solvents) and in the limiting simple case of a dipolar solute embedded in a spherical cavity. The results point out the complementarity of LR and SS approaches and the advantages of the latter model especially for polar solvents. The favorable scaling properties of PCM-TD-DFT models in both SS and LR variants and their availability in effective quantum mechanical codes pave the route for the computation of reliable spectroscopic properties of large molecules of technological and/or biological interest in their natural environments.

Biophysics↗

The convergence of complete active space self-consistent-field configuration interaction including all single and double excitation energies to the complete basis set limit.

Examination of the convergence of full valence complete active space self-consistent-field configuration interaction including all single and double excitation (CASSCF-CISD) energies with expansion of the one-electron basis set reveals a pattern very similar to the convergence of single determinant energies. Calculations on the lowest four singlet states and the lowest four triplet states of N(2) with the sequence of n-tuple-zeta augmented polarized (nZaP) basis sets (n=2, 3, 4, 5, and 6) are used to establish the complete basis set limits. Full configuration-interaction (CI) and core electron contributions must be included for very accurate potential energy surfaces. However, a simple extrapolation scheme that has no adjustable parameters and requires nothing more demanding than CAS(10e(-),8orb)-CISD/3ZaP calculations gives the R(e), omega(e), omega(e)X(e), T(e), and D(e) for these eight states with rms errors of 0.0006 Angstrom, 4.43 cm(-1), 0.35 cm(-1), 0.063 eV, and 0.018 eV, respectively.

Journal Article↗

Is the protein surrounding the active site critical for hydrogen peroxide reduction by selenoprotein glutathione peroxidase? An ONIOM study.

In this ONIOM(QM:MM) study, we evaluate the role of the protein surroundings in the mechanism of H2O2 reduction catalyzed by the glutathione peroxidase enzyme, using the whole monomer (3113 atoms in 196 amino acid residues) as a model. A new optimization scheme that allows the full optimization of transition states for large systems has been utilized. It was found that in the presence of the surrounding protein the optimized active site structure bears a closer resemblance to the one in the X-ray structure than that without the surrounding protein. H2O2 reduction occurs through a two-step mechanism. In the first step, the selenolate anion (E-Se(-)) formation occurs with a barrier of 16.4 kcal/mol and is endothermic by 12.0 kcal/mol. The Gln83 residue plays the key role of the proton abstractor, which is in line with the experimental suggestion. In the second step, the O-O bond is cleaved, and selenenic acid (R-Se-OH) and a water molecule are formed. The calculated barrier for this process is 6.0 kcal/mol, and it is exothermic by 80.9 kcal/mol. The overall barrier of 18.0 kcal/mol for H2O2 reduction is in reasonable agreement with the experimentally measured barrier of 14.9 kcal/mol. The protein surroundings has been calculated to exert a net effect of only 0.70 kcal/mol (in comparison to the "active site only" model including solvent effects) on the overall barrier, which is most likely due to the active site being located at the enzyme surface.

Binding Sites↗

Geometries and properties of excited states in the gas phase and in solution: theory and application of a time-dependent density functional theory polarizable continuum model.

In this paper we present the theory and implementation of analytic derivatives of time-dependent density functional theory (TDDFT) excited states energies, both in vacuo and including solvent effects by means of the polarizable continuum model. The method is applied to two case studies: p-nitroaniline and 4-(dimethyl)aminobenzonitrile. For both molecules PCM-TDDFT is shown to be successful in supporting the analysis of experimental data with useful insights for a better understanding of photophysical and photochemical pathways in solution.

Journal Article↗

Singlet excited-state behavior of uracil and thymine in aqueous solution: a combined experimental and computational study of 11 uracil derivatives.

The excited-state properties of uracil, thymine, and nine other derivatives of uracil have been studied by steady-state and time-resolved spectroscopy. The excited-state lifetimes were measured using femtosecond fluorescence upconversion in the UV. The absorption and emission spectra of five representative compounds have been computed at the TD-DFT level, using the PBE0 exchange-correlation functional for ground- and excited-state geometry optimization and the Polarizable Continuum Model (PCM) to simulate the aqueous solution. The calculated spectra are in good agreement with the experimental ones. Experiments show that the excited-state lifetimes of all the compounds examined are dominated by an ultrafast (<100 fs) component. Only 5-substituted compounds show more complex behavior than uracil, exhibiting longer excited-state lifetimes and biexponential fluorescence decays. The S(0)/S(1) conical intersection, located at CASSCF (8/8) level, is indeed characterized by pyramidalization and out of plane motion of the substituents on the C5 atom. A thorough analysis of the excited-state Potential Energy Surfaces, performed at the PCM/TD-DFT(PBE0) level in aqueous solution, shows that the energy barrier separating the local S(1) minimum from the conical intersection increases going from uracil through thymine to 5-fluorouracil, in agreement with the ordering of the experimental excited-state lifetime.

Models, Molecular↗

Nonresonant optical activity of isolated organic molecules.

Cavity ring-down polarimetry (CRDP) has been exploited to interrogate the nonresonant optical activity (or circular birefringence) of prototypical organic compounds in the vapor phase, thereby revealing the intrinsic chiro-optical response evoked from isolated (solvent-free) molecules. Specific polarization rotation parameters have been measured at two distinct excitation wavelengths (355 nm and 633 nm) for a variety of gas-phase species drawn from the terpene, epoxide, and alkane/alkene families, with complementary solution-phase polarimetric studies serving to highlight the pronounced influence of solute-solvent interactions. Time-dependent linear response calculations performed at high levels of density functional theory have been enlisted to unravel the structural and electronic origins for observed behavior. Aside from elucidating the complex solvation processes that mediate chiro-optical phenomena taking place in condensed media, this study affords a critical assessment for emerging ab initio predictions of nonresonant optical activity and for their promising ability to assist in the determination of absolute molecular stereochemistry.

Epichlorohydrin↗

Ab initio Ehrenfest dynamics.

We present an ab initio direct Ehrenfest dynamics scheme using a three time-step integrator. The three different time steps are implemented with nuclear velocity Verlet, nuclear-position-coupled midpoint Fock integrator, and time-dependent Hartree-Fock with a modified midpoint and unitary transformation algorithm. The computational cost of the ab initio direct Ehrenfest dynamics presented here is found to be only a factor of 2-4 larger than that of Born-Oppenheimer (BO) dynamics. As an example, we compute the vibration of the NaCl molecule and the intramolecular torsional motion of H2C=NH2+ by Ehrenfest dynamics compared with BO dynamics. For the vibration of NaCl with an initial kinetic energy of 1.16 eV, Ehrenfest dynamics converges to BO dynamics with the same vibrational frequency. The intramolecular rotation of H2C=NH2+ produces significant electronic excitation in the Ehrenfest trajectory. The amount of nonadiabaticity, suggested by the amplitude of the coherent progression of the excited and ground electronic states, is observed to be directly related to the strength of the electron-nuclear coupling. Such nonadiabaticity is seen to have a significant effect on the dynamics compared with the adiabatic approximation.

Journal Article↗

The convergence of complete active space self-consistent-field energies to the complete basis set limit.

Examination of the convergence of full valence complete active space self-consistent-field energies with expansion of the one-electron basis set reveals a pattern very similar to the convergence of single determinant Hartree-Fock energies. Calculations on 26 molecular examples with the sequence of ntuple-zeta augmented polarized (nZaP) basis sets (n=2, 3, 4, 5, and 6) are used to evaluate complete basis set extrapolation schemes. The most effective extrapolation reduces the rms one-electron basis set truncation errors from 3.03, 0.58, and 0.12 mhartree to 0.23, 0.05, and 0.014 mhartree for the 3ZaP, 4ZaP, and 5ZaP basis sets, respectively.

Journal Article↗

Excited states of porphyrin isomers and porphycene derivatives: a SAC-CI study.

Excited states of free-base porphyrin isomers, porphycene (Pc), corrphycene (Cor), and hemiporphycene (hPc), were studied by the Symmetry-Adapted Cluster (SAC)/SAC-Configuration Interaction (CI) method. The absorption peaks of the porphyrin isomers were assigned on the basis of the SAC-CI spectra. The X, Y, X', and Y' bands of the porphyrin isomers, which have weak intensities, are identified. The differences in the Q-band absorptions among the isomers were clearly explained by the four-orbital model. In Cor and hPc, the wave function of the B-band corresponds to the mixture of the four-orbital excitations and the optically forbidden excitation of free-base porphin (P), due to the molecular symmetry lowering in the isomers. The B-band character is described by the five-orbital model in Pc and the six-orbital model in Cor and hPc. Two tetrazaporphycenes and two ring-extended (dibenzo) porphycenes were designed, and the Q-band transition moment was successfully controlled. These examples show that the control of the four-orbital energy levels is the guiding principle for pigment design in porphyrin compounds.

Electrochemistry↗

Effect of time-dependent basis functions and their superposition error on atom-centered density matrix propagation (ADMP): connections to wavelet theory of multiresolution analysis.

We present a rigorous analysis of the primitive Gaussian basis sets used in the electronic structure theory. This leads to fundamental connections between Gaussian basis functions and the wavelet theory of multiresolution analysis. We also obtain a general description of basis set superposition error which holds for all localized, orthogonal or nonorthogonal, basis functions. The standard counterpoise correction of quantum chemistry is seen to arise as a special case of this treatment. Computational study of the weakly bound water dimer illustrates that basis set superposition error is much less for basis functions beyond the 6-31+G(*) level of Gaussians when structure, energetics, frequencies, and radial distribution functions are to be calculated. This result will be invaluable in the use of atom-centered Gaussian functions for ab initio molecular dynamics studies using Born-Oppenheimer and atom-centered density matrix propagation.

Journal Article↗

Analytical energy gradient of the symmetry-adapted-cluster configuration-interaction general-R method for singlet to septet ground and excited states.

A method of calculating analytical energy gradients of the singlet and triplet excited states, ionized states, electron-attached states, and high-spin states from quartet to septet states by the symmetry-adapted-cluster configuration-interaction general-R method is developed and implemented. This method is a powerful tool in the studies of geometries, dynamics, and properties of the states of molecules in which not only one-electron processes but also two- and multielectron processes are involved. The performance of the present method was confirmed by calculating the geometries and the spectroscopic constants of the diatomic and polyatomic molecules in various electronic states involving the ground state and the one- to three-electron excited states. The accurate descriptions were obtained for the equilibrium geometries, vibrational frequencies, and adiabatic excitation energies, which show the potential usefulness of the present method. The particularly interesting applications were to the C' 1Ag state of acetylene, the A 2Deltau and B 2Sigmau+ states of CNC and the 4B1 and a 4Piu states of N3 radical.

Journal Article↗

Geometry optimization with QM/MM, ONIOM, and other combined methods. I. Microiterations and constraints.

Hybrid energy methods such as QM/MM and ONIOM, that combine different levels of theory into one calculation, have been very successful in describing large systems. Geometry optimization methods can take advantage of the partitioning of these calculations into a region treated at a quantum mechanical (QM) level of theory and the larger, remaining region treated by an inexpensive method such as molecular mechanics (MM). A series of microiterations can be employed to fully optimize the MM region for each optimization step in the QM region. Cartesian coordinates are used for the MM region and are chosen so that the internal coordinates of the QM region remain constant during the microiterations. The coordinates of the MM region are augmented to permit rigid body translation and rotation of the QM region. This is essential if any atoms in the MM region are constrained, but it also improves the efficiency of unconstrained optimizations. Because of the microiterations, special care is needed for the optimization step in the QM region so that the system remains in the same local valley during the course of the optimization. The optimization methodology with microiterations, constraints, and step-size control are illustrated by calculations on bacteriorhodopsin and other systems.

Journal Article↗