Interference, fluctuation, and alternation of electron tunneling in protein media. 2. Non-condon theory for the energy gap dependence of electron transfer rate

Developing the quantum transition rate theory of Prezhdo and Rossky (J. Chem. Phys. 1997, 107, 5863), we produced a new non-Condon theory of the rate of electron transfer (ET) which happens through a protein medium with conformational fluctuation. The new theory is expressed by a convolution form of the power spectrum for the autocorrelation function of the electronic tunneling matrix element T(DA)(t) with quantum correction and the ordinary Franck-Condon factor. The new theory satisfies the detailed balance condition for the forward and backward ET rates. The ET rate formula is divided into two terms of elastic and inelastic tunneling mechanisms on the mathematical basis. The present theory is applied to the ET from Bph(-) to Q(A) in the reaction center of Rhodobacter sphaeroides. Numerical calculations of T(DA)(t) were made by a combined method of molecular dynamics simulations and quantum chemistry calculations. We showed that the normalized autocorrelation function of T(DA)(t) is almost expressed by exponential forms. The calculated energy gap law of the ET rate is nearly Marcus' parabola in most of the normal region and around the maximum region, but it does not decay substantially in the inverted region, which is called the anomalous inverted region. We also showed that the energy gap law at the high uphill energy gap in the normal region is elevated considerably from the Marcus' parabola, which is called the anomalous normal region. Those anomalous energy gap laws are due to the inelastic tunneling mechanism which works actively at the energy gap far from zero. We presented an empirical formula for easily calculating the non-Condon ET rate, which is usable by many researchers. We provided experimental evidence for the anomalous inverted region which was basically reproduced by the present theory. The present theory was extensively compared with the previous non-Condon theories.     

Mechanism of intramolecular electron transfer in the photoexcited Zn-substituted cytochrome c: theoretical and experimental perspective

Photoinduced electron transfer (ET) in zinc-substituted cytochrome c (Zn-cyt c) has been utilized in many studies on the long-range ET in protein. Attempting to understand its ET mechanism in terms of electronic structure of the molecule, we have calculated an all-electron wave function for the ground-state of Zn-cyt c on the basis of density functional theory (DFT). The four molecular orbitals (MOs) responsible for excitation by UV-vis light (Gouterman's 4-orbitals) are assigned on the basis of the excited states of chromophore model for Zn-porphine complex calculated with the time-dependent DFT method. ET rates between each Gouterman's 4-orbitals and other MOs were estimated using Fermi's golden rule. It appeared that the two occupied MOs of the 4-orbitals show exclusively higher ET rate from/to particular MOs that localize on outermost amino acid residues (Lys 7 or Asn 54), respectively, whereas ET rates involving the two unoccupied MOs of the 4-orbitals are much slower. These results imply that the intramolecular ET in photoexcited Zn-cyt c is governed by the hole transfer through occupied MOs. The couplings of MOs between zinc porphyrin core and specific amino acid residues on the protein surface have been demonstrated in Zn-cyt c immobilized on an Au electrode via carboxylic acid group-terminated self-assembled monolayer. The Zn-cyt c-modified electrode showed photocurrents responsible for photoillumination. The action spectrum of the photocurrent was identical with the absorption spectrum of Zn-cyt c, indicating photoinduced electron conduction via occupied MOs. The voltage dependence of the photocurrent appeared to be linear and bidirectional like a photoconductor, which strongly supports the intramolecular ET mechanism in Zn-cyt c proposed on the basis of the theoretical calculations.     

An efficient, fragment-based electronic structure method for molecular systems: self-consistent polarization with perturbative two-body exchange and dispersion

We report a fragment-based electronic structure method, intended for the study of clusters and molecular liquids, that incorporates electronic polarization (induction) in a self-consistent fashion but treats intermolecular exchange and dispersion interactions perturbatively, as post-self-consistent field corrections, using a form of pairwise symmetry-adapted perturbation theory. The computational cost of the method scales quadratically as a function of the number of fragments (monomers), but could be made to scale linearly by exploiting distance-dependent thresholds. Extensive benchmark calculations are reported using the S22 database of high-level ab initio binding energies for dimers, and we find that average errors can be reduced to <1 kcal/mol with a suitable choice of basis set. Comparison to ab initio benchmarks for water clusters as large as (H(2)O)(20) demonstrates that the method recovers ?90% of the binding energy in these systems, at a tiny fraction of the computational cost. As such, this approach represents a promising path toward accurate, systematically improvable, and parameter-free simulation of molecular liquids.     

Interference, fluctuation, and alternation of electron tunneling in protein media. 1. Two tunneling routes in photosynthetic reaction center alternate due to thermal fluctuation of protein conformation

Electron tunneling routes for the electron transfer from the bacteriopheophytin anion to the primary quinone in the bacterial photosynthetic reaction center of Rhodobactor sphaeroides are investigated by a combined method of molecular dynamics simulations for the protein conformation fluctuation and quantum chemical calculations for the electronic states of the donor, acceptor, and protein medium. The analysis of the tunneling route is made by mapping interatomic electron tunneling currents for each protein conformation. We found that there are two dominant routes mainly passing through Trp(M252) (Trp route) or mainly passing through Met(M218) (Met route). Actual electron tunneling pathways alternate between the two routes, depending on the protein conformation which varies with time. When either the Trp route or the Met route dominates, the electron tunneling matrix element /T(DA)/ becomes large. When both the Trp route and the Met route dominate, /T(DA)/ becomes very small due to the destructive interference of the electron tunneling currents between the two routes. We found that a linear relationship exists between the value of /T(DA)/ and the inverse of the degree of destructive interference Q for a wide range of values (ca. 3-10(3) for Q). A similar relationship was also found previously for electron transfer in ruthenium-modified azurins, suggesting that this relationship holds true in general. From these results, we are led to the conclusion that /T(DA)/ cannot exceed a maximum value at Q = 1, even if much variation of /T(DA)/ happens due to the fluctuation of protein conformation. We also conclude that the property of the electron transfer alternates between constructive and destructive interference, due to the fluctuation of protein conformation. It is impossible to keep a system in either constructive or destructive interference because thermal fluctuation of protein conformation takes place.     

Fragment-Based Ab Initio Molecular Dynamics Simulation for Combustion

We develop a fragment-based ab initio molecular dynamics (FB-AIMD) method for efficient dynamics simulation of the combustion process. In this method, the intermolecular interactions are treated by a fragment-based many-body expansion in which three- or higher body interactions are neglected, while two-body interactions are computed if the distance between the two fragments is smaller than a cutoff value. The accuracy of the method was verified by comparing FB-AIMD calculated energies and atomic forces of several different systems with those obtained by standard full system quantum calculations. The computational cost of the FB-AIMD method scales linearly with the size of the system, and the calculation is easily parallelizable. The method is applied to methane combustion as a benchmark. Detailed reaction network of methane reaction is analyzed, and important reaction species are tracked in real time. The current result of methane simulation is in excellent agreement with known experimental findings and with prior theoretical studies.     

Spectral-product methods for electronic structure calculations

Progress is reported in development, implementation, and application of a spectral method for ab initio studies of the electronic structure of matter. In this approach, antisymmetry restrictions are enforced subsequent to construction of the many-electron Hamiltonian matrix in a complete orthonormal spectral-product basis. Transformation to a permutation-symmetry representation obtained from the eigenstates of the aggregate electron antisymmetrizer is seen to enforce the requirements of the Pauli principle ex post facto, and to eliminate the unphysical (non-Pauli) states spanned by the product representation. Results identical with conventional use of prior antisymmetrization of configurational state functions are obtained in applications to many-electron atoms. The development provides certain advantages over conventional methods for polyatomic molecules, and, in particular, facilitates incorporation of fragment information in the form of Hermitian matrix representatives of atomic and diatomic operators which include the non-local effects of overall electron antisymmetry. An exact atomic-pair expression is obtained in this way for polyatomic Hamiltonian matrices which avoids the ambiguities of previously described semi-empirical fragment-based methods for electronic structure calculations. Illustrative applications to the well-known low-lying doublet states of the H₃ molecule in a minimal-basis-set demonstrate that the eigensurfaces of the antisymmetrizer can anticipate the structures of the more familiar energy surfaces, including the seams of intersection common in high-symmetry molecular geometries. The calculated H₃ energy surfaces are found to be in good agreement with corresponding valence-bond results which include all three-center terms, and are in general accord with accurate values obtained employing conventional high-level computational-chemistry procedures. By avoiding the repeated evaluations of the many-centered one- and two-electron integrals required in construction of polyatomic Hamiltonian matrices in the antisymmetric basis states commonly employed in conventional calculations, and by performing the required atomic and atomic-pair calculations once and for all, the spectral-product approach may provide an alternative potentially efficient ab initio formalism suitable for computational studies of adiabatic potential energy surfaces more generally. Contribution to the Mark S. Gordon 65th Birthday Festschrift Issue.     

Fragment-based quantum mechanical methods for periodic systems with Ewald summation and mean image charge convention for long-range electrostatic interactions

We describe an Ewald-summation method to incorporate long-range electrostatic interactions into fragment-based electronic structure methods for periodic systems. The present method is an extension of the particle-mesh Ewald technique for combined quantum mechanical and molecular mechanical (QM/MM) calculations, and it has been implemented into the explicit polarization (X-Pol) potential to illustrate the computational details. As in the QM/MM-Ewald method, the X-Pol-Ewald approach is a linear-scaling electrostatic method, in which the short-range electrostatic interactions are determined explicitly in real space and the long-range Ewald pair potential is incorporated into the Fock matrix as a correction. To avoid the time-consuming Fock matrix update during the self-consistent field procedure, a mean image charge (MIC) approximation is introduced, in which the running average with a user-chosen correlation time is used to represent the long-range electrostatic correction as an average effect. Test simulations on liquid water show that the present X-Pol-Ewald method takes about 25% more CPU time than the usual X-Pol method using spherical cutoff, whereas the use of the MIC approximation reduces the extra costs for long-range electrostatic interactions by 15%. The present X-Pol-Ewald method provides a general procedure for incorporating long-range electrostatic effects into fragment-based electronic structure methods for treating biomolecular and condensed-phase systems under periodic boundary conditions.     

Spin–orbit and correlation effects in platinum hydride (PtH)

The low-lying electronic states of PtH were studied by all-electron one- and two-component variational calculations on the multireference CI levels. The orbital optimization is performed within a one-component formalism, whereas the further refinement of the wave functions follows two different schemes: The most demanding approach introduces spin–orbit coupling in the CI optimization step, giving a simultaneous treatment of electron correlation and spin–orbit coupling. The second, considerably less demanding approach, corresponds almost to a perturbational treatment, introducing spin–orbit coupling as a final step after the CI optimization by diagonalizing the resulting Hamiltonian matrix over CI states. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 68: 53–64, 1998     

Generalized relativistic effective core potentials for superheavy elements

The generalized relativistic effective core potentials (GRECPs) for calculations of electronic structure and physical-chemical properties of compounds containing superheavy elements (Z ≥ 104) are presented. Features of accounting for the finite nuclear size effects which are unusually large for superheavy elements are discussed in details. Accuracy of the GRECPs is analyzed in atomic calculations compared to all-electron studies with the Dirac-Coulomb-Breit Hamiltonian. Applications of the GRECP method in molecular and cluster calculations are surveyed.     

A localized molecular-orbital assembler approach for Hartree-Fock calculations of large molecules

We describe an alternative fragment-based method, the localized molecular-orbital assembler method, for Hartree-Fock (HF) calculations of macromolecules. In this approach, a large molecule is divided into many small-size fragments, each of which is capped by its local surroundings. Then the conventional HF calculations are preformed on these capped fragments (or subsystems) and the canonical molecular orbitals of these systems are transferred into localized molecular orbitals (LMOs). By assembling the LMOs of these subsystems into a set of LMOs of the target molecule, the total density matrix of the target molecule is constructed and correspondingly the HF energy or other molecular properties can be approximately computed. This approach computationally achieves linear scaling even for medium-sized systems. Our test calculations with double-zeta and polarized double-zeta basis sets demonstrate that the present approach is able to reproduce the conventional HF energies within a few millihartrees for a broad range of molecules.     

Ab initio derivation of multi-orbital extended Hubbard model for molecular crystals

From configuration interaction (CI) ab initio calculations, we derive an effective two-orbital extended Hubbard model based on the gerade (g) and ungerade (u) molecular orbitals (MOs) of the charge-transfer molecular conductor (TTM-TTP)I(3) and the single-component molecular conductor [Au(tmdt)(2)]. First, by focusing on the isolated molecule, we determine the parameters for the model Hamiltonian so as to reproduce the CI Hamiltonian matrix. Next, we extend the analysis to two neighboring molecule pairs in the crystal and we perform similar calculations to evaluate the inter-molecular interactions. From the resulting tight-binding parameters, we analyze the band structure to confirm that two bands overlap and mix in together, supporting the multi-band feature. Furthermore, using a fragment decomposition, we derive the effective model based on the fragment MOs and show that the staking TTM-TTP molecules can be described by the zig-zag two-leg ladder with the inter-molecular transfer integral being larger than the intra-fragment transfer integral within the molecule. The inter-site interactions between the fragments follow a Coulomb law, supporting the fragment decomposition strategy.     

Linear-scaling fixed-node diffusion quantum Monte Carlo: accounting for the nodal information in a density matrix-based scheme

A reformulation of the fixed-node diffusion quantum Monte Carlo method (FN-DQMC) in terms of the N-particle density matrix is presented, which allows us to reduce the computational effort to linear for the evaluation of the local energy. The reformulation is based on our recently introduced density matrix-based approach for a linear-scaling variational QMC method [J. Kussmann et al., Phys. Rev. B. 75, 165107 (2007)]. However, within the latter approach of using the positive semi-definite N-particle trial density (rhoN T(R)=mid R:Psi(T)(R)mid R:(2)), the nodal information of the trial function is lost. Therefore, a straightforward application to the FN-DQMC method is not possible, in which the sign of the trial function is usually traced in order to confine the random walkers to their nodal pockets. As a solution, we reformulate the FN-DQMC approach in terms of off-diagonal elements of the N-particle density matrix rhoN T(R;R'), so that the nodal information of the trial density matrix is obtained. Besides all-electron moves, a scheme to perform single-electron moves within N-PDM QMC is described in detail. The efficiency of our method is illustrated for exemplary calculations.     

Independent trajectory implementation of the semiclassical Liouville method: application to multidimensional reaction dynamics

We describe an independent trajectory implementation of semiclassical Liouville method for simulating quantum processes using classical trajectories. In this approach, a single ensemble of trajectories describes all semiclassical density matrix elements of a coupled electronic state problem, with the ensemble evolving classically under a single reference Hamiltonian chosen on the basis of physical grounds. In this paper, we introduce an additional uncoupled trajectory approximation, allowing the members of the ensemble to evolve independently of one another and eliminating the major computational costs of our previous coupled trajectory implementation. The accuracy of the method is demonstrated for model one-dimensional problems. In addition, the approach is applied to the chemical reaction dynamics of a collinear triatomic system, yielding excellent agreement with exact calculations. This method allows molecular dynamics involving coupled electronic surfaces to be modeled with essentially the same effort as classical molecular dynamics and ensemble averaging.     

Finite group theory for large systems. 3. Symmetry-generation of reduced matrix elements for icosahedral C(20) and C(60) molecules

This paper uses symmetry-generation to simplify the determination of Hamiltonian reduced matrix elements. It is part of a series on using computers to apply finite group theory to quantum mechanical calculations on large systems. Symmetry-generation is an expression of the whole molecule as a sum of symmetry transformations on a smaller reference structure. Then on a suitably-conditioned symmetry-adapted basis, the reduced matrix elements of the Hamiltonian are averages of certain elements of the simpler reference structure matrix. The smaller the reference structure, the greater is the computational savings. Single atom reference structures are used here for the Hückel treatment of icosahedral C(20) and C(60) fullerenes. The analytical power of this approach is illustrated by determining the two bond lengths of C(60) from spectral data.     

Two-body Reduced Density Matrix Reconstruction for Van der Waals Systems

A method for the calculation of the electronic energy of a correlated system is presented. This approach is based on the reconstruction of the total two-body reduced density matrix by doing separate configurations interaction calculations on fragments. The method has been tested on Van der Waals systems and has been implemented by considering restrictive N-representability conditions. It is shown that the computational strategy presented in this work can describe with good accuracy weak dispersion interactions, and considerably lowers the size-consistency error of a classical configuration interaction calculation.     

On relativistic effects in ground state potential curves of Zn2, Cd2, and Hg2 dimers. A CCSD(T) study

The ground state potential curves of the Zn2, Cd2, and Hg2 dimers calculated at different levels of theory are presented and compared with each other as well as with experimental and other theoretical studies. The calculations at the level of Dirac-Coulomb Hamiltonian (DCH), 4-component spin-free Hamiltonian, nonrelativistic Lévy-Leblond Hamiltonian and at the level of simple Coulombic correction to DCH are presented. The potential curves are calculated in an all-electron supermolecular approach including the correction to basis set superposition error (BSSE). Electron correlation is treated at the coupled cluster level including single and double excitations and noniterative triple corrections, CCSD(T). In addition, simulations of the temperature dependence of dynamic viscosities in the low-density limit using the obtained ground state potential curves are presented.     

Low-temperature electronic transport through macromolecules and characteristics of intramolecular electron transfer

Long-distance electron transfer (ET) plays an important part in many biological processes. Also, fundamental understanding of ET processes could give grounds for designing miniaturized electronic devices. So far, experimental data on the ET mostly concern ET rates which characterize ET processes as a whole. Here, we develop a different approach which could provide more information about intrinsic characteristics of the long-range intramolecular ET. A starting point of the studies is an obvious resemblance between ET processes and electric transport through molecular wires placed between metallic contacts. Accordingly, the theory of electronic transport through molecular wires is applied to analyze characteristics of a long-range electron transfer through molecular bridges. Assuming a coherent electron tunneling to be a predominant mechanism of ET at low temperatures, it is shown that low-temperature current-voltage characteristics could exhibit a special structure, and the latter contains information concerning intrinsic features of the intramolecular ET. Using the Buttiker dephasing model within the scattering matrix formalism, we analyze the effect of dephasing on the electron transmission function and current-voltage curves.     

A new fragment-based approach for calculating electronic excitation energies of large systems

We present a new fragment-based scheme to calculate the excited states of large systems without necessity of a Hartree-Fock (HF) solution of the whole system. This method is based on the implementation of the renormalized excitonic method [M. A. Hajj et al., Phys. Rev. B 72, 224412 (2005)] at ab initio level, which assumes that the excitation of the whole system can be expressed by a linear combination of various local excitations. We decomposed the whole system into several blocks and then constructed the effective Hamiltonians for the intra- and inter-block interactions with block canonical molecular orbitals instead of widely used localized molecular orbitals. Accordingly, we avoided the prerequisite HF solution and the localization procedure of the molecular orbitals in the popular local correlation methods. Test calculations were implemented for hydrogen molecule chains at the full configuration interaction, symmetry adapted cluster/symmetry adapted cluster configuration interaction, HF/configuration interaction singles (CIS) levels and more realistic polyene systems at the HF/CIS level. The calculated vertical excitation energies for lowest excited states are in reasonable accordance with those determined by the calculations of the whole systems with traditional methods, showing that our new fragment-based method can give good estimates for low-lying energy spectra of both weak and moderate interaction systems with economic computational costs.