Application of the local self-consistent-field method to core-ionized and core-excited molecules, polymers, and proteins: True orthogonality between ground and excited states

Initial Maximum Overlap Method Embedded with Extremely Localized Molecular Orbitals for Core-Ionized States of Large Systems

Despite great advances in X-ray absorption spectroscopy for the investigation of small molecule electronic structure, the application to biosystems of experimental techniques developed within this research field remains a challenge. To partially circumvent the problem, users resort to theoretical methods to interpret or predict the X-ray absorption spectra of large molecules. To accomplish this task, only low-cost computational strategies can be exploited. For this reason, some of them are single Slater determinant wavefunction approaches coupled with multiscale embedding techniques designed to treat large systems of biological interest. Therefore, in this work, we propose to apply the recently developed IMOM/ELMO embedding method to the determination of core-ionized states. The IMOM/ELMO technique resulted from the combination of the single Slater determinant Δself-consistent-field-initial maximum overlap approach (ΔSCF-IMOM) with the QM/ELMO (quantum mechanics/extremely localized molecular orbital) embedding strategy, a method where only the chemically relevant region of the examined system is treated at fully quantum chemical level, while the rest is described through transferred and frozen extremely localized molecular orbitals (ELMOs). The IMOM/ELMO technique was initially validated by computing core-ionization energies for small molecules, and it was afterwards exploited to study larger biosystems. The obtained results are in line with those reported in previous studies that applied alternative ΔSCF approaches. This makes us envisage a possible future application of the proposed method to the interpretation of X-ray absorption spectra of large molecules.

Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: a Comparison between Energy- and Variance-Based Approaches

Self-consistent field methods for excited states offer an attractive low-cost route to study not only excitation energies but also properties of excited states. Here, we present the generalization of two self-consistent field methods, the maximum overlap method (MOM) and the σ-SCF method, to calculate excited states in strong magnetic fields and investigate their stability and accuracy in this context. These methods use different strategies to overcome the well-known variational collapse of energy-based optimizations to the lowest solution of a given symmetry. The MOM tackles this problem in the definition of the orbital occupations to constrain the self-consistent field procedure to converge on excited states, while the σ-SCF method is based on the minimization of the variance instead of the energy. To overcome the high computational cost of the variance minimization, we present a new implementation of the σ-SCF method with the resolution of identity approximation, allowing the use of large basis sets, which is an important requirement for calculations in strong magnetic fields. The accuracy of these methods is assessed by comparison with the benchmark literature data for He, H₂, and CH⁺. The results reveal severe limitations of the variance-based scheme, which become more acute in large basis sets. In particular, many states are not accessible using variance optimization. Detailed analysis shows that this is a general feature of variance optimization approaches due to the masking of local minima in the optimization. In contrast, the MOM shows promising performance for computing excited states under these conditions, yielding results consistent with available benchmark data for a diverse range of electronic states.

Relativistic EOM-CCSD for Core-Excited and Core-Ionized State Energies Based on the Four-Component Dirac-Coulomb(-Gaunt) Hamiltonian

We report an implementation of the core-valence separation approach to the four-component relativistic Hamiltonian-based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD) for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects of the different CVS-EOM-CCSD variants and the use of different Hamiltonians based on the exact two-component (X2C) framework on the energies of different core-ionized and -excited states in halogen- (CH₃I, HX, and X^-, X = Cl-At) and xenon-containing (Xe, XeF₂) species. Our results show that the X2C molecular mean-field approach [Sikkema, J.; J. Chem. Phys. 2009, 131, 124116], based on four-component Dirac-Coulomb mean-field calculations (²DC^M), is capable of providing core excitations and ionization energies that are nearly indistinguishable from the reference four-component energies for up to and including fifth-row elements. We observe that two-electron integrals over the small-component basis sets lead to non-negligible contributions to core binding energies for the K and L edges for atoms such as iodine or astatine and that the approach based on Dirac-Coulomb-Gaunt mean-field calculations (²DCG^M) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals.

The nuts and bolts of core-hole constrainedab initiosimulation forK-shell x-ray photoemission and absorption spectra

X-ray photoemission (XPS) and near edge x-ray absorption fine structure (NEXAFS) spectroscopy play an important role in investigating the structure and electronic structure of materials and surfaces.Ab initiosimulations provide crucial support for the interpretation of complex spectra containing overlapping signatures. Approximate core-hole simulation methods based on density functional theory (DFT) such as the delta-self-consistent-field (ΔSCF) method or the transition potential (TP) method are widely used to predictK-shell XPS and NEXAFS signatures of organic molecules, inorganic materials and metal-organic interfaces at reliable accuracy and affordable computational cost. We present the numerical and technical details of our variants of the ΔSCF and TP method (coined ΔIP-TP) to simulate XPS and NEXAFS transitions. Using exemplary molecules in gas-phase, in bulk crystals, and at metal-organic interfaces, we systematically assess how practical simulation choices affect the stability and accuracy of simulations. These include the choice of exchange-correlation functional, basis set, the method of core-hole localization, and the use of periodic boundary conditions (PBC). We particularly focus on the choice of aperiodic or periodic description of systems and how spurious charge effects in periodic calculations affect the simulation outcomes. For the benefit of practitioners in the field, we discuss sensible default choices, limitations of the methods, and future prospects.

Block-Localized Excitation for Excimer Complex and Diabatic Coupling

We describe a block-localized excitation (BLE) method to carry out constrained optimization of block-localized orbitals for constructing valence bond-like, diabatic excited configurations using multistate density functional theory (MSDFT). The method is an extension of the previous block-localized wave function method through a fragment-based ΔSCF approach to optimize excited determinants within a molecular complex. In BLE, both the number of electrons and the electronic spin of different fragments in a whole system can be constrained, whereas electrostatic, exchange, and polarization interactions among different blocks can be fully taken into account of. To avoid optimization collapse to unwanted states, a ΔSCF projection scheme and a maximum overlap of wave function approach have been presented. The method is illustrated by the excimer complex of two naphthalene molecules. With a minimum of eight spin-adapted configurational state functions, it was found that the inversion of L_a^- and L_b^- states near the optimal structure of the excimer complex is correctly produced, which is in quantitative agreement with DMRG-CASPT2 calculations and experiments. Trends in the computed transfer integrals associated with excited-state energy transfer both in the singlet and triplet states are discussed. The results suggest that MSDFT may be used as an efficient approach to treat intermolecular interactions in excited states with a minimal active space (MAS) for interpretation of the results and for dynamic simulations, although the selection of a small active space is often system dependent.

Targeting G-quadruplexes with Organic Dyes: Chelerythrine-DNA Binding Elucidated by Combining Molecular Modeling and Optical Spectroscopy

The DNA-binding of the natural benzophenanthridine alkaloid chelerythrine (CHE) has been assessed by combining molecular modeling and optical absorption spectroscopy. Specifically, both double-helical (B-DNA) and G-quadruplex sequences—representative of different topologies and possessing biological relevance, such as telomeric or regulatory sequences—have been considered. An original multiscale protocol, making use of molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations, allowed us to compare the theoretical and experimental circular dichroism spectra of the different DNA topologies, readily providing atomic-level details of the CHE–DNA binding modes. The binding selectivity towards G-quadruplexes is confirmed by both experimental and theoretical determination of the binding free energies. Overall, our mixed computational and experimental approach is able to shed light on the interaction of small molecules with different DNA conformations. In particular, CHE may be seen as the building block of promising drug candidates specifically targeting G-quadruplexes for both antitumoral and antiviral purposes.

On the orthogonality of states with approximate wavefunctions

A variational solution to the eigenvalue problem for the Hamiltonian H, with orthogonality restrictions on eigenvectors of H to the vector H ∣ Φ₀〉, where ∣Φ₀〉 is an approximate ground-state vector, is proposed as a means to calculate excited states. The asymptotic projection (AP) method proposed previously is further developed and applied to solve this problem in a simple way. We demonstrate that the AP methodology does not require an evaluation of the matrix elements of operator H², whereas conventional approaches-such as the elimination of off-diagonal Lagrange multipliers method, projection operator techniques, and other methods-do. It is shown, based on the results obtained for the single-electron molecular ions H₂⁺, HeH²⁺, and H₃²⁺, that applying the new method to determine excited-state wavefunctions yields the upper bounds for excited-state energies. We demonstrate that regardless of whether the orthogonality constraint for states (〈Φ| Φ₀〉 = 0) is applied, the zero-coupling constraint (〈Φ| H| Φ₀〉 = 0) is imposed, or both of these restrictions are enforced simultaneously, practically the same excited-state energy is obtained if the basis set is almost complete. For the systems considered here, all schemes are capable of giving a sub-μhartree level of accuracy for the ground and excited states computed with different basis sets.

Photophysics of chlorin e6: from one- and two-photon absorption to fluorescence and phosphorescence

Linear and non-linear optical properties of a known photosensitizer producing singlet oxygen, chlorin e6, have been studied, including dynamics effects.

From non-covalent binding to irreversible DNA lesions: nile blue and nile red as photosensitizing agents

We report a molecular modeling study, coupled with spectroscopy experiments, on the behavior of two well known organic dyes, nile blue and nile red, when interacting with B-DNA. In particular, we evidence the presence of two competitive binding modes, for both drugs. However their subsequent photophysical behavior is different and only nile blue is able to induce DNA photosensitization via an electron transfer mechanism. Most notably, even in the case of nile blue, its sensitization capabilities strongly depend on the environment resulting in a single active binding mode: the minor groove. Fluorescence spectroscopy confirms the presence of competitive interaction modes for both sensitizers, while the sensitization via electron transfer, is possible only in the case of nile blue.

Two-photon-absorption DNA sensitization via solvated electron production: unraveling photochemical pathways by molecular modeling and simulation

Infrared optical properties and photochemistry of a potential therapeutic dye interacting with DNA: a multiscale computational approach.

Hydrogen abstraction by photoexcited benzophenone: consequences for DNA photosensitization

We report a computational investigation of the hydrogen abstraction (H-abstraction) induced by triplet benzophenone (³BP) on thymine nucleobase and backbone sugar.

Gentlest ascent dynamics for calculating first excited state and exploring energy landscape of Kohn-Sham density functionals

We develop the gentlest ascent dynamics for Kohn-Sham density functional theory to search for the index-1 saddle points on the energy landscape of the Kohn-Sham density functionals. These stationary solutions correspond to excited states in the ground state functionals. As shown by various examples, the first excited states of many chemical systems are given by these index-1 saddle points. Our novel approach provides an alternative, more robust way to obtain these excited states, compared with the widely used ΔSCF approach. The method can be easily generalized to target higher index saddle points. Our results also reveal the physical interest and relevance of studying the Kohn-Sham energy landscape.

Theoretical Investigation of the Geometries and UV-vis Spectra of Poly(l-glutamic acid) Featuring a Photochromic Azobenzene Side Chain

The geometries and UV-vis spectra of azobenzene dyes grafted as a side chain on poly(l-glutamic acid) have been investigated using a combination of quantum mechanics/molecular mechanics (QM/MM) and time-dependent density functional theory (TD-DFT) methods at the TD-PBE0/6-311+G(d,p)//B3LYP/6-311G(d,p):Amber ff99 level of theory. The influence of the secondary structure of the polypeptide on the electronic properties of both the trans and cis conformations of azobenzene dyes has been studied. It turns out that the grafted dyes exhibit a red-shift of the π → π* absorption energies mainly due to the auxochromic shift induced by the peptidic group used to link the chromophoric unit to the polypeptide and that specific interactions between the glutamic side chain and the azobenzene moiety lead to a large blue-shift of the n → π* transition.

Photophysics of acetophenone interacting with DNA: why the road to photosensitization is open

Deoxyribonucleic acid photosensitization, i.e. the photoinduced electron- or energy-transfer of chromophores interacting with DNA, is a crucial phenomenon that triggers important DNA lesions such as pyrimidine dimerization, even upon absorption of relatively low-energy radiation. Oxidative lesions may also be produced via the photoinduced production of reactive oxygen species. Aromatic ketones, and acetophenone in particular, are well known for their sensitization effects. In this contribution we model the structural and dynamical properties of the acetophenone/DNA aggregates as well as their spectroscopic and photophysical properties using high-level hybrid quantum mechanics/molecular mechanics methods. We show that the key steps of the photochemistry of acetophenone in gas phase are conserved in the macromolecular environment and thus an ultrafast singlet-triplet conversion of acetophenone is expected prior to the transfer to DNA.

Projected hybrid orbitals: a general QM/MM method

A projected hybrid orbital (PHO) method was described to model the covalent boundary in a hybrid quantum mechanical and molecular mechanical (QM/MM) system. The PHO approach can be used in ab initio wave function theory and in density functional theory with any basis set without introducing system-dependent parameters. In this method, a secondary basis set on the boundary atom is introduced to formulate a set of hybrid atomic orbtials. The primary basis set on the boundary atom used for the QM subsystem is projected onto the secondary basis to yield a representation that provides a good approximation to the electron-withdrawing power of the primary basis set to balance electronic interactions between QM and MM subsystems. The PHO method has been tested on a range of molecules and properties. Comparison with results obtained from QM calculations on the entire system shows that the present PHO method is a robust and balanced QM/MM scheme that preserves the structural and electronic properties of the QM region.

Multireference space without first solving the configuration interaction problem

We further develop an idea to generate a compact multireference space without first solving the configuration interaction problem previously proposed for the ground state (GS) (Glushkov, Chem. Phys. Lett. 1995, 244, 1). In the present contribution, our attention is focused on low-lying excited states (ESs) with the same symmetry as the GS which can be adequately described in terms of an high-spin open-shell formalism. Two references Møller-Plesset (MP) like perturbation theory for ESs is developed. It is based on: (1) a main reference configuration constructed from the parent molecular orbitals adjusted to a given ES and (2) secondary double excitation configuration built on the GS like orbitals determined by the Hartree-Fock equations subject to some orthogonality constraints. It is shown how to modify the MP zeroth-order Hamiltonian so that the reference configurations and corresponding excitations are eigenfunctions of it and are compatible with orthogonality conditions for the GS and ES. Intruder states appearance is also discussed. The proposed scheme is applied to the GS, ES, and excitation energies of small molecules to illustrate and calibrate our calculations.

Orthogonality constrained density functional theory for electronic excited states

We report a novel scheme for computing electronic excitation energies within the framework of density functional theory (DFT) based on a time-independent variational formulation of DFT. The excited state density functional is recast as a Kohn-Sham functional, which is further simplified by an adiabatic approximation of the exchange-correlation functional. Under the adiabatic approximation, the minimization of the excited state Kohn-Sham functional is shown to be equivalent to a ground state DFT computation augmented with orthogonality constraints with respect to the ground state Kohn-Sham determinant. An algorithm for the optimization of the energy subject to orthogonality constraints, which does not suffer from variational collapse, is described and implemented. A benchmark test set containing 28 organic molecules (Schreiber, M. J. Chem. Phys. 2008, 128, 134110) was used to assess the quality of the excitation energies obtained. Two novel approaches to spin-adapt the resulting excitation energies are discussed and found to provide results with error metrics similar to those of time-dependent DFT. Similarities and differences with respect to other time-independent DFT approaches are highlighted and some of the advantages of our scheme-including the ability to correctly describe charge-transfer excitations-are critically assessed.

Theoretical modeling of large molecular systems. Advances in the local self consistent field method for mixed quantum mechanics/molecular mechanics calculations

Molecular mechanics methods can efficiently compute the macroscopic properties of a large molecular system but cannot represent the electronic changes that occur during a chemical reaction or an electronic transition. Quantum mechanical methods can accurately simulate these processes, but they require considerably greater computational resources. Because electronic changes typically occur in a limited part of the system, such as the solute in a molecular solution or the substrate within the active site of enzymatic reactions, researchers can limit the quantum computation to this part of the system. Researchers take into account the influence of the surroundings by embedding this quantum computation into a calculation of the whole system described at the molecular mechanical level, a strategy known as the mixed quantum mechanics/molecular mechanics (QM/MM) approach. The accuracy of this embedding varies according to the types of interactions included, whether they are purely mechanical or classically electrostatic. This embedding can also introduce the induced polarization of the surroundings. The difficulty in QM/MM calculations comes from the splitting of the system into two parts, which requires severing the chemical bonds that link the quantum mechanical subsystem to the classical subsystem. Typically, researchers replace the quantoclassical atoms, those at the boundary between the subsystems, with a monovalent link atom. For example, researchers might add a hydrogen atom when a C-C bond is cut. This Account describes another approach, the Local Self Consistent Field (LSCF), which was developed in our laboratory. LSCF links the quantum mechanical portion of the molecule to the classical portion using a strictly localized bond orbital extracted from a small model molecule for each bond. In this scenario, the quantoclassical atom has an apparent nuclear charge of +1. To achieve correct bond lengths and force constants, we must take into account the inner shell of the atom: for an sp(3) carbon atom, we consider the two core 1s electrons and treat that carbon as an atom with three electrons. This results in an LSCF+3 model. Similarly, a nitrogen atom with a lone pair of electrons available for conjugation is treated as an atom with five electrons (LSCF+5). This approach is particularly well suited to splitting peptide bonds and other bonds that include carbon or nitrogen atoms. To embed the induced polarization within the calculation, researchers must use a polarizable force field. However, because the parameters of the usual force fields include an average of the induction effects, researchers typically can obtain satisfactory results without explicitly introducing the polarization. When considering electronic transitions, researchers must take into account the changes in the electronic polarization. One approach is to simulate the electronic cloud of the surroundings by a continuum whose dielectric constant is equal to the square of the refractive index. This Electronic Response of the Surroundings (ERS) methodology allows researchers to model the changes in induced polarization easily. We illustrate this approach by modeling the electronic absorption of tryptophan in human serum albumin (HSA).

On orthogonality constrained multiple core-hole states and optimized effective potential method

An attempt to construct a multiple core-hole state within the optimized effective potential (OEP) methodology is presented. In contrast to the conventional Δ-self-consistent field method for hole states, the effects of removing an electron is achieved using some orthogonality constraints imposed on the orbitals so that a Slater determinant describing a hole state is constrained to be orthogonal to that of a neutral system. It is shown that single, double, and multiple core-hole states can be treated within a unified framework and can be easily implemented for atoms and molecules. For this purpose, a constrained OEP method proposed earlier for excited states (Glushkov and Levy, J. Chem. Phys. 2007, 126, 174106) is further developed to calculate single and double core ionization energies using a local effective potential expressed as a direct mapping of the external potential. The corresponding equations, determining core-hole orbitals from a one-particle Schrödinger equation with a local potential as well as correlation corrections derived from the second-order many-body perturbation theory are given. One of the advantages of the present direct mapping formulation is that the effective potential, which plays the role of the Kohn-Sham potential, has the symmetry of the external potential. Single and double core ionization potentials computed with the presented scheme were found to be in agreement with data available from experiment and other calculations. We also discuss core-hole state local potentials for the systems studied.

Doubly, triply, and multiply excited states from a constrained optimized effective potential method

This article further develops and applies a constrained optimized effective potential (COEP) approach for the practical calculations of doubly and multiply excited states of atoms and molecules. The COEP method uses the time-independent theory of pure excited states and implements a simple asymptotic projection method to take orthogonality constraints into account. We show that, in contrast with the common time-dependent density functional method, the COEP methodology is capable of treating doubly, triply, and multiply excited states and can be easily applied to both atoms and molecules. In particular, doubly excited energies of each state are calculated through a constrained minimization procedure including constraints that make its Slater determinantal functions orthogonal to those of the ground and all lower-lying doubly excited states. The performance of the proposed method is examined by calculations of doubly excited state energies for the He atom and H(2) molecule at exchange-only and exchange-correlation level of approximation.

On the suitability of strictly localized orbitals for hybrid QM/MM calculations

In the QM/MM method we have developed (LSCF/MM), the QM and the MM parts are held together by means of strictly localized bonding orbitals (SLBOs). Generally these SLBOs are derived from localized bond orbitals (LBOs) that undergo tails deletion, resulting in a nonpredictable change of their properties. An alternative set of SLBOs is provided by the extremely localized molecular orbitals (ELMOs) approach, where the orbitals are rigorously localized on some prefixed atoms without tails on the other atoms of the molecule. A comparative study of SLBOs arising from various localization schemes and ELMOs is presented to test the reliability and the transferability of these functions within the Local Self-Consistent Field (LSCF) framework. Two types of chemical bonds were considered: C--C and C--O single bonds. The localized functions are obtained on the ethane and the methanol molecules, and are tested on beta-alanine and diethyl ether molecules. Moreover, the various protonation forms of beta-alanine have been investigated to illustrate how well the polarity variation of the chemical bond can be handled throughout a chemical process. At last, rotation energy profiles around C--C and C--O bonds are reproduced for butane and fluoromethanol. Energetic, geometric, as well as electronic factors all indicate that ELMO functions are much more transferable from one molecule to another, leading to results closer to the usual SCF reference than any other calculations involving any other localized orbitals. When the shape of the orbital is the most important factor then ELMO functions will perform as well as any other localized orbital.