Biomedical subjects

Troy Van Voorhis

Publications and source records attributed to Troy Van Voorhis.

8 recordsLinked to original sources

PubMed2006

Extracting electron transfer coupling elements from constrained density functional theory.

Constrained density functional theory (DFT) is a useful tool for studying electron transfer (ET) reactions. It can straightforwardly construct the charge-localized diabatic states and give a direct measure of the inner-sphere reorganization energy. In this work, a method is presented for calculating the electronic coupling matrix element (Hab) based on constrained DFT. This method completely avoids the use of ground-state DFT energies because they are known to irrationally predict fractional electron transfer in many cases. Instead it makes use of the constrained DFT energies and the Kohn-Sham wave functions for the diabatic states in a careful way. Test calculations on the Zn2+ and the benzene-Cl atom systems show that the new prescription yields reasonable agreement with the standard generalized Mulliken-Hush method. We then proceed to produce the diabatic and adiabatic potential energy curves along the reaction pathway for intervalence ET in the tetrathiafulvalene-diquinone (Q-TTF-Q) anion. While the unconstrained DFT curve has no reaction barrier and gives Hab approximately 17 kcal/mol, which qualitatively disagrees with experimental results, the Hab calculated from constrained DFT is about 3 kcalmol and the generated ground state has a barrier height of 1.70 kcal/mol, successfully predicting (Q-TTF-Q)- to be a class II mixed-valence compound.

Electrons↗

PubMed2006

Direct calculation of electron transfer parameters through constrained density functional theory.

It is shown that constrained density functional theory (DFT) can be used to access diabatic potential energy surfaces in the Marcus theory of electron transfer, thus providing a means to directly calculate the driving force and the inner-sphere reorganization energy. We present in this report an analytic expression for the forces in constrained DFT and their implementation in geometry optimization, a prerequisite for the calculation of electron transfer parameters. The method is then applied to study the symmetric mixed-valence complex tetrathiafulvalene-diquinone radical anion, which is observed experimentally to be a Robin-Day class II compound but found by DFT to be in class III. Constrained DFT avoids this pitfall of over-delocalization and provides a way to find the charge-localized structure. In another application, driving forces and inner-sphere reorganization energies are calculated for the charge recombination (CR) reactions in formanilide-anthraquinone (FA-AQ) and ferrocene-formanilide-anthraquinone (Fc-FA-AQ). While the two compounds have similar reorganization energies, the driving force in FA-AQ is 1 eV larger than in Fc-FA-AQ, in agreement with experimental observations and supporting the experimental conclusion that the anomalously long-lived FA-AQ charge-separated state arises because the electron transfer is in the Marcus inverted region.

Anions↗

PubMed2006

Advances in methods and algorithms in a modern quantum chemistry program package.

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.

Algorithms↗

PubMed2006

A coherent state approach to semiclassical nonadiabatic dynamics.

A semiclassical (SC) approximation to the quantum mechanical propagator for nonadiabatic systems is derived. Our derivation starts with an exact path integral expression that uses canonical coherent states for the nuclear degrees of freedom and spin coherent states for the electronic degrees of freedom. A stationary path approximation (SPA) is then applied to the path integral to obtain the SC approximation. The SPA results in complex classical trajectories of both nuclear and electronic degrees of freedom and a double ended boundary condition. The root search problem is solved using the previously proposed "real trajectory local search" algorithm. The SC approximation is tested on three simple one dimensional two-state systems proposed by Tully [J. Chem. Phys. 93, 1061 (1990)], and the SC results are compared to Ehrenfest and surface hopping predictions. Excellent agreement with quantum results is reached when the SC trajectory is far away from caustics. We discuss the origin of caustics in this SC formalism and the strengths and weaknesses of this approach.

Journal Article↗

PubMed2006

Accurate magnetic exchange couplings in transition-metal complexes from constrained density-functional theory.

We demonstrate an accurate method for extracting Heisenberg exchange-coupling constants (J) from density-functional theory (DFT) calculations. We note that the true uncoupled low-spin state of a given molecule should be identified with the ground state of the system subject to a constraint on the spin density of the atoms. Using an efficient optimization strategy for constrained DFT we obtain these states directly, leading to a simple, physically motivated formula for J. Our method only depends on state energies and their associated electron densities and assigns no unphysical meaning to the Kohn-Sham determinant or individual orbitals. We study several bimetallic transition-metal complexes and find that the constrained DFT approach is competitive with, if not better than, the best broken symmetry DFT results. The success of constrained DFT in these cases appears to result from a balanced elimination of self-interaction error and static correlation from the simulation.

Journal Article↗

PubMed2005

Density-matrix renormalization-group algorithms with nonorthogonal orbitals and non-Hermitian operators, and applications to polyenes.

We describe the theory and implementation of two extensions to the density-matrix renormalization-group (DMRG) algorithm in quantum chemistry: (i) to work with an underlying nonorthogonal one-particle basis (using a biorthogonal formulation) and (ii) to use non-Hermitian and complex operators and complex wave functions, which occur naturally in biorthogonal formulations. Using these developments, we carry out ground-state calculations on ethene, butadiene, and hexatriene, in a polarized atomic-orbital basis. The description of correlation in these systems using a localized nonorthogonal basis is improved over molecular-orbital DMRG calculations, and comparable to or better than coupled-cluster calculations, although we encountered numerical problems associated with non-Hermiticity. We believe that the non-Hermitian DMRG algorithm may further become useful in conjunction with other non-Hermitian Hamiltonians, for example, similarity-transformed coupled-cluster Hamiltonians.

Algorithms↗

PubMed2005

Fluctuation-dissipation theorem density-functional theory.

Using the fluctuation-dissipation theorem (FDT) in the context of density-functional theory (DFT), one can derive an exact expression for the ground-state correlation energy in terms of the frequency-dependent density response function. When combined with time-dependent density-functional theory, a new class of density functionals results that use approximations to the exchange-correlation kernel fxc as input. This FDT-DFT scheme holds promise to solve two of the most distressing problems of conventional Kohn-Sham DFT: (i) It leads to correlation energy functionals compatible with exact exchange, and (ii) it naturally includes dispersion. The price is a moderately expensive O(N6) scaling of computational cost and a slower basis set convergence. These general features of FDT-DFT have all been recognized previously. In this paper, we present the first benchmark results for a set of molecules using FDT-DFT beyond the random-phase approximation (RPA)-that is, the first such results with fxc not equal to 0. We show that kernels derived from the adiabatic local-density approximation and other semilocal functionals suffer from an "ultraviolet catastrophe," producing a pair density that diverges at small interparticle distance. Nevertheless, dispersion interactions can be treated accurately if hybrid functionals are employed, as is demonstrated for He2 and HeNe. We outline constraints that future approximations to fxc should satisfy and discuss the prospects of FDT-DFT.

Journal Article↗

PubMed2004

Semiclassical representations of electronic structure and dynamics.

We use a new formulation of the semiclassical coherent state propagator to derive and evaluate several different approximate representations of electron dynamics. For each representation we examine: (1) its ability to treat quantum effects and electron correlation, (2) its expected scaling with system size, and (3) the types of systems for which it can be used. We also apply two of the methods to a pair of model problems, namely the minimal basis electron dynamics in H2 and the magnetization dynamics in a cluster model of the Kagome lattice, in order to verify the feasibility of these approaches for realistic systems. Based on all these criteria, we find that the representation that takes the electron spins as the classical variables is particularly promising for the quantitative and qualitative description of large systems.

Journal Article↗