Benchmarking the completely renormalised equation-of-motion coupled-cluster approaches for vertical excitation energies

Unraveling the effect of aromaticity for the dynamics of excited states of single benzene fluorophores

The photophysical properties of a series of recently synthesized single benzene fluorophores were investigated using ensemble density functional theory calculations. The energetic stability of the ground and excited state species were counterposed against the aromaticity index derived from local vibrational modes. It was found that the large Stokes shift of the fluorophores (up to ca. 5800 cm  - 1 ) originates from the effect of electron donating and electron withdrawing substituents rather than π -delocalization and related (anti-)aromaticity. On the basis of nonadiabatic molecular dynamics simulations, the absence of fluorescence from one of the regioisomers was explained by the occurrence of easily accessible S  1 /S  0 conical intersections below the vertical excitation energy level. It is demonstrated in the manuscript that the analysis of local mode force constants and the related aromaticity index represent a useful tool for the characterization of π -delocalization effects in π -conjugated compounds.

Aufbau Suppressed Coupled Cluster Theory for Electronically Excited States

We introduce an approach to improve single-reference coupled cluster theory in settings where the Aufbau determinant is absent from or plays only a small role in the true wave function. Using a de-excitation operator that can be efficiently hidden within a similarity transform, we create a coupled cluster wave function in which de-excitations work to suppress the Aufbau determinant and produce wave functions dominated by other determinants. Thanks to an invertible and fully exponential form, the approach is systematically improvable, size consistent, size extensive, and, interestingly, size intensive in a granular way that should make the adoption of some ground state techniques, such as local correlation, relatively straightforward. In this initial study, we apply the general formalism to create a state-specific method for orbital-relaxed, singly excited states. We find that this approach matches the accuracy of similar-cost equation-of-motion methods in valence excitations while offering improved accuracy for charge transfer states. We also find the approach to be more accurate than excited-state-specific perturbation theory in both types of states.

Leveling the Mountain Range of Excited-State Benchmarking through Multistate Density Functional Theory

The performance of multistate density functional theory (MSDFT) with nonorthogonal state interaction (NOSI) is assessed for 100 vertical excitation energies against the theoretical best estimates extracted to the full configuration interaction accuracy on the database developed by Loos et al. in 2018 (Loos2018). Two optimization techniques, namely, block-localized excitation and target state optimization, are examined along with two ways of estimating the transition density functional (TDF) for the correlation energy of the Hamiltonian matrix density functional. The results from the two optimization methods are similar. It was found that MSDFT-NOSI using the spin-multiplet degeneracy constraint for the TDF of spin-coupling interaction, along with the M06-2X functional, yields a root-mean-square error (RMSE) of 0.22 eV, which performs noticeably better than time-dependent density functional theory (DFT) at an RMSE of 0.43 eV using the same functional and basis set on the Loos2018 database. In comparison with wave function theory, NOSI has smaller errors than CIS(D∞), LR-CC2, and ADC(3) all of which have an RMSE of 0.28 eV, but somewhat greater than STEOM-CCSD (RMSE of 0.14 eV) and LR-CCSD (RMSE of 0.11 eV) wave function methods. In comparison with Kohn-Sham (KS) DFT calculations, the multistate DFT approach has little double counting of correlation. Importantly, there is no noticeable difference in the performance of MSDFT-NOSI on the valence, Rydberg, singlet, triplet, and double-excitation states. Although the use of another hybrid functional PBE0 leads to a greater RMSE of 0.36 eV, the deviation is systematic with a linear regression slope of 0.994 against the results with M06-2X. The present benchmark reveals that density functional approximations developed for KS-DFT for the ground state with a noninteracting reference may be adopted in MSDFT calculations in which the state interaction is key.

Excited-State-Specific Pseudoprojected Coupled-Cluster Theory

We present an excited-state-specific coupled-cluster approach in which both the molecular orbitals and cluster amplitudes are optimized for an individual excited state. The theory is formulated via a pseudoprojection of the traditional coupled-cluster wavefunction that allows correlation effects to be introduced atop an excited-state mean field starting point. The approach shares much in common with ground-state CCSD, including size extensivity and an N6 cost scaling. Preliminary numerical tests show that, when augmented with N5 cost perturbative corrections for key terms, the method can improve over excited-state-specific second-order perturbation theory in valence, charge transfer, and Rydberg states.

Unitary Coupled Cluster: Seizing the Quantum Moment

Shallow, CNOT-efficient quantum circuits are crucial for performing accurate computational chemistry simulations on current noisy quantum hardware. Here, we explore the usefulness of noniterative energy corrections, based on the method of moments of coupled-cluster theory, for accelerating convergence toward full configuration interaction. Our preliminary numerical results relying on iteratively constructed ansätze suggest that chemically accurate energies can be obtained with substantially more compact circuits, implying enhanced resilience to gate and decoherence noise.

The vertical excitation energies and a lifetime of the two lowest singlet excited states of the conjugated polyenes from C2 to C22: Ab initio, DFT, and semiclassical MNDO-MD simulations

Electronic excited states in the series of polyene molecules were explored. Optimal ground-state geometry was used for the evaluation of vertical excitation energies. Results of a chosen set of functionals were compared to post-HF methods (EOM-CCSD, NEVPT2, CASPT2, and MRCI). In addition, the semiempirical OM2/MNDO method using MRCISD computational level was confronted with the above-mentioned techniques. Despite the fact that the first excited state has a significant double-excitation character some functionals were able to qualitatively determine the correct state order (where the lowest excited state has a A g - character). The most successful functionals in transition energies predictions were PBE, TPSS and BLYP in Tamm-Dancoff approach (TDA), which had the smallest root-mean-square deviation (RMSD) scoring towards the experimental values. Regarding RMSD scoring, the OM2/MNDO method performed fairly well, too. Besides absorption spectra, lifetimes of the first two excited states were estimated based on a stochastic approach exploring a swarm of OM2/MNDO hopping dynamics using the Tully fewest switch algorithm for each molecule. The longest lifetime of the first excited state (S₁ ) was found for decapentaene (about 5 ps). Further elongation of the conjugated chain caused a mild decrease of this value to ca 1.5 ps for docosaundecaene.

Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.

A hybrid coupled cluster-machine learning algorithm: Development of various regression models and benchmark applications

The iterative solution of the coupled cluster equations exhibits a synergistic relationship among the various cluster amplitudes. The iteration scheme is analyzed as a multivariate discrete time propagation of nonlinearly coupled equations, which is dictated by only a few principal cluster amplitudes. These principal amplitudes usually correspond to only a few valence excitations, whereas all other cluster amplitudes are enslaved and behave as auxiliary variables [Agarawal et al., J. Chem. Phys. 154, 044110 (2021)]. We develop a coupled cluster-machine learning hybrid scheme where various supervised machine learning strategies are introduced to establish the interdependence between the principal and auxiliary amplitudes on-the-fly. While the coupled cluster equations are solved only to determine the principal amplitudes, the auxiliary amplitudes, on the other hand, are determined via regression as unique functionals of the principal amplitudes. This leads to significant reduction in the number of independent degrees of freedom during the iterative optimization, which saves significant computation time. A few different regression techniques have been developed, which have their own advantages and disadvantages. The scheme has been applied to several molecules in their equilibrium and stretched geometries, and our scheme, with all the regression models, shows a significant reduction in computation time over the canonical coupled cluster calculations without unduly sacrificing the accuracy.

Internal Conversion between Bright (11Bu+) and Dark (21Ag-) States in s-trans-Butadiene and s-trans-Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1¹B_u⁺ and 2¹A_g^-, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N_o}, and DEA-EOMCC(4p2h){N_u} or DEA-EOMCC(3p1h,4p2h){N_u}, consistently predict that the 1¹B_u⁺/2¹A_g^- crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1¹B_u⁺ and 2¹A_g^- potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT2, do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR(4,4), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

An approximate coupled cluster theory via nonlinear dynamics and synergetics: The adiabatic decoupling conditions

The coupled cluster iteration scheme is analyzed as a multivariate discrete time map using nonlinear dynamics and synergetics. The nonlinearly coupled set of equations to determine the cluster amplitudes are driven by a fraction of the entire set of cluster amplitudes. These driver amplitudes enslave all other amplitudes through a synergistic inter-relationship, where the latter class of amplitudes behave as the auxiliary variables. The driver and the auxiliary variables exhibit vastly different time scales of relaxation during the iteration process to reach the fixed points. The fast varying auxiliary amplitudes are small in magnitude, while the driver amplitudes are large, and they have a much longer time scale of relaxation. Exploiting their difference in relaxation time scale, we employ an adiabatic decoupling approximation, where each of the fast relaxing auxiliary modes is expressed as a unique function of the principal amplitudes. This results in a tremendous reduction in the independent degrees of freedom. On the other hand, only the driver amplitudes are determined accurately via exact coupled cluster equations. We will demonstrate that the iteration scheme has an order of magnitude reduction in computational scaling than the conventional scheme. With a few pilot numerical examples, we would demonstrate that this scheme can achieve very high accuracy with significant savings in computational time.

Excited states in the adiabatic connection fluctuation-dissipation theory: Recovering missing correlation energy from the negative part of the density response spectrum

The adiabatic connection (AC) theory offers an alternative to the perturbation theory methods for computing correlation energy in the multireference wavefunction framework. We show that the AC correlation energy formula can be expressed in terms of the density linear response function as a sum of components related to positive and negative parts of the transition energy spectrum. Consequently, generalization of the adiabatic connection fluctuation-dissipation theory to electronically excited states is obtained. The component of the linear response function related to the negative-transition energy enters the correlation energy expression with an opposite sign to that of the positive-transition part and is non-negligible in the description of excited states. To illustrate this, we analyze the approximate AC model in which the linear response function is obtained in the extended random phase approximation (ERPA). We demonstrate that AC can be successfully combined with the ERPA for excited states, provided that the negative-excitation component of the response function is rigorously accounted for. The resulting AC0D model, an extension of the AC0 scheme introduced in our earlier works, is applied to a benchmark set of singlet excitation energies of organic molecules. AC0D constitutes a significant improvement over AC0 by bringing the excitation energies of the lowest excited states to a satisfactory agreement with theoretical best estimates, which parallels or even exceeds the accuracy of the n-electron valence state perturbation theory method. For higher excitations, AC0D is less reliable due to the gradual deterioration of the underlying ERPA linear response.

Improved Description of Charge-Transfer Potential Energy Surfaces via Spin-Component-Scaled CC2 and ADC(2) Methods

The molecular level understanding of electronic transport properties depends on the reliable theoretical description of charge-transfer (CT)-type electronic states. In this paper, the performance of spin-component-scaled variants of the popular CC2 and ADC(2) methods is evaluated for CT states, following benchmark strategies of earlier studies that revealed a compromised accuracy of the unmodified models. In addition to statistics on the accuracy of vertical excitation energies at equilibrium and infinite separation of bimolecular complexes, potential energy surfaces of the ammonia-fluorine complex are also reported. The results show the capability of spin-component-scaled approaches to reduce the large errors of their regular counterparts to a significant extent, outperforming even the coupled-cluster single and double method in many cases. The cost-effective scaled-opposite-spin variants are found to provide a remarkably good agreement with the CCSDT-3 reference data, thereby being recommended methods of choice in the study of charge-transfer states.

A hybrid approach to excited-state-specific variational Monte Carlo and doubly excited states

We extend our hybrid linear-method/accelerated-descent variational Monte Carlo optimization approach to excited states and investigate its efficacy in double excitations. In addition to showing a superior statistical efficiency when compared to the linear method, our tests on small molecules show good energetic agreement with benchmark methods. We also demonstrate the ability to treat double excitations in systems that are too large for a full treatment by using selected configuration interaction methods via an application to 4-aminobenzonitrile. Finally, we investigate the stability of state-specific variance optimization against collapse to other states' variance minima and find that symmetry, Ansatz quality, and sample size all have roles to play in achieving stability.

Isoenergetic two-photon excitation enhances solvent-to-solute excited-state proton transfer

Two-photon excitation (TPE) is an attractive means for controlling chemistry in both space and time. Since isoenergetic one- and two-photon excitations (OPE and TPE) in non-centrosymmetric molecules are allowed to reach the same excited state, it is usually assumed that they produce similar excited-state reactivity. We compare the solvent-to-solute excited-state proton transfer of the super photobase FR0-SB following isoenergetic OPE and TPE. We find up to 62% increased reactivity following TPE compared to OPE. From steady-state spectroscopy, we rule out the involvement of different excited states and find that OPE and TPE spectra are identical in non-polar solvents but not in polar ones. We propose that differences in the matrix elements that contribute to the two-photon absorption cross sections lead to the observed enhanced isoenergetic reactivity, consistent with the predictions of our high-level coupled-cluster-based computational protocol. We find that polar solvent configurations favor greater dipole moment change between ground and excited states, which enters the probability for TPE as the absolute value squared. This, in turn, causes a difference in the Franck-Condon region reached via TPE compared to OPE. We conclude that a new method has been found for controlling chemical reactivity via the matrix elements that affect two-photon cross sections, which may be of great utility for spatial and temporal precision chemistry.

A New Benchmark Set for Excitation Energy of Charge Transfer States: Systematic Investigation of Coupled Cluster Type Methods

The numerous existing publications on benchmarking quantum chemistry methods for excited states rarely include Charge Transfer (CT) states, although many interesting phenomena in, e.g., biochemistry and material physics involve the transfer of electrons between fragments of the system. Therefore, it is timely to test the accuracy of quantum chemical methods for CT states, as well. In this study we first propose a new benchmark set consisting of dimers having low-energy CT states. On this set, the vertical excitation energy has been calculated with Coupled Cluster methods including triple excitations (CC3, CCSDT-3, CCSD(T)(a)*), as well as with methods including full or approximate doubles (CCSD, STEOM-CCSD, CC2, ADC(2), EOM-CCSD(2)). The results show that the popular CC2 and ADC(2) methods are much less accurate for CT states than for valence states. On the other hand, EOM-CCSD seems to have similar systematic overestimation of the excitation energies for both types of states. Among the triples methods the novel EOM-CCSD(T)(a)* method including noniterative triple excitations is found to stand out with its consistently good performance for all types of states, delivering essentially EOM-CCSDT quality results.

Recent developments in the general atomic and molecular electronic structure system

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized.

Quantum computation solves a half-century-old enigma: Elusive vibrational states of magnesium dimer found

The high-lying vibrational states of the magnesium dimer (Mg₂), which has been recognized as an important system in studies of ultracold and collisional phenomena, have eluded experimental characterization for half a century. Until now, only the first 14 vibrational states of Mg₂ have been experimentally resolved, although it has been suggested that the ground-state potential may support five additional levels. Here, we present highly accurate ab initio potential energy curves based on state-of-the-art coupled-cluster and full configuration interaction computations for the ground and excited electronic states involved in the experimental investigations of Mg₂. Our ground-state potential unambiguously confirms the existence of 19 vibrational levels, with ~1 cm^-1 root mean square deviation between the calculated rovibrational term values and the available experimental and experimentally derived data. Our computations reproduce the latest laser-induced fluorescence spectrum and provide guidance for the experimental detection of the previously unresolved vibrational levels.

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules

Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.

Proton Abstraction Mediates Interactions between the Super Photobase FR0-SB and Surrounding Alcohol Solvent

We report on the motional and proton transfer dynamics of the super photobase FR0-SB in the series of normal alcohols C1 (methanol) through C8 (n-octanol) and ethylene glycol. Steady-state and time-resolved fluorescence data reveal that the proton abstraction dynamics of excited FR0-SB depend on the identity of the solvent and that the transfer of the proton from solvent to FR0-SB*, forming FR0-HSB⁺*, fundamentally alters the nature of interactions between the excited molecule and its surroundings. In its unprotonated state, solvent interactions with FR0-SB* are consistent with slip limit behavior, and in its protonated form, intermolecular interactions are consistent with a much stronger interaction of FR0-HSB⁺* with the deprotonated solvent RO^-. We understand the excited-state population dynamics in the context of a kinetic model involving a transition state wherein FR0-HSB⁺* is still bound to the negatively charged alkoxide, prior to solvation of the two charged species. Data acquired in ethylene glycol confirm the hypothesis that the rotational diffusion dynamics of FR0-SB* are largely mediated by solvent viscosity while proton transfer dynamics are mediated by the lifetime of the transition state. Taken collectively, our results demonstrate that FR0-SB* extracts solvent protons efficiently and in a predictable manner, consistent with a ca. 3-fold increase in dipole moment upon photoexcitation as determined by ab initio calculations based on the equation-of-motion coupled-cluster theory.

Benchmarks for Electronically Excited States with CASSCF Methods

The accuracy of three different complete active space (CAS) self-consistent field (CASSCF) methods is investigated for the electronically excited-state benchmark set of Schreiber , M. ; et al. J. Chem. Phys. 2008 , 128 , 134110 . Comparison of the CASSCF linear response (LR) methods MC-RPA and MC-TDA and the state-averaged (SA) CASSCF method is made for 122 singlet excitation energies and 69 oscillator strengths. Of all CASSCF methods, when considering the complete test set, MC-RPA performs best for both excitation energies and oscillator strengths with a mean absolute error (MAE) of 0.74 eV and 51%, respectively. MC-TDA and SA-CASSCF show a similar accuracy for the excitation energies with a MAE of ∼1 eV with respect to more accurate coupled cluster (CC3) excitation energies. The opposite trend is observed for the subset of n → π* excitation energies for which SA-CASSCF exhibits the least deviations (MAE 0.65 eV). By looking at s-tetrazine in more detail, we conclude that better performance for the n → π* SA-CASSCF excitation energies can be attributed to a fortunate error compensation. For oscillator strengths, SA-CASSCF performs worst for the complete test set (MAE 100%) as well as for the subsets of n → π* (MAE 192%) and π → π* excitations (MAE 84.9%). In general, CASSCF gives the worst performance for excitation energies of all excited-state ab initio methods considered so far due to lacking the major part of dynamic electron correlation, though MC-RPA and TD-DFT (BP86) show similar performance. Among all LR-type methods, LR-CASSCF oscillator strengths are the ones with the least accuracy for the same reason. As state-specific orbital relaxation effects are accounted for in LR-CASSCF, oscillator strengths are significantly more accurate than those of MS-CASPT2. Our findings should encourage further developments of response theory-based multireference methods with higher accuracy and feasibility.

The ring-opening channel and the influence of Rydberg states on the excited state dynamics of furan and its derivatives

One important relaxation pathway for photo-excited five-membered heterocyclic organic molecules is ring-opening via a dissociative πσ^* state. In this study, we investigate the influence of this pathway in furan and several hydrogenated and methylated derivatives by combining time-resolved photoelectron spectroscopy with time-dependent density functional theory and coupled cluster calculations. We find strong experimental evidence that the ring-opening channel is the major relaxation channel in furan, 2,3-dihydrofuran, and 2-methylfuran (2-MF). In 2,5-dimethylfuran (25-DMF), however, we observe that the molecules relax either via a π3s Rydberg state or through a direct return to the ground state by undergoing ring-puckering motions. From the supporting calculations, for 2-MF and 25-DMF, we predict that there is strong mixing between the πσ^* state and the π3s Rydberg state along the ring opening pathway. However, in 25-DMF, no crossing between the πσ^*/π3s state and the initially excited ππ^* state can be found along the ring opening coordinate, effectively blocking this channel.

Combined quantum-mechanical molecular mechanics calculations with NWChem and AMBER: Excited state properties of green fluorescent protein chromophore analogue in aqueous solution

Combined quantum mechanical molecular mechanics (QM/MM) calculations have become a popular methodology for efficient and accurate description of large molecular systems. In this work we introduce our development of a QM/MM framework based on two well-known codes-NWChem and AMBER. As an initial application area we are focused on excited state properties of small molecules in an aqueous phase using an analogue of the green fluorescent protein (GFP) chromophore as a particular test case. Our approach incorporates high level coupled cluster theory for the analysis of excited states providing a reliable theoretical analysis of effects of an aqueous solvation environment on the photochemical properties of the GFP chromophore. Using a systematic approach, which involves comparison of gas phase and aqueous phase results for different protonation states and conformations, we resolve existing uncertainties regarding the theoretical interpretation of experimental data. We observe that the impact of aqueous environment on charged states generally results in blue shifts of the absorption spectra, but the magnitude of the effect is sensitive to both protonation state and conformation and can be rationalized based on charge movement into the area of higher/lower external electrostatic potentials. At neutral pH levels the experimentally observed absorption signal is most likely coming from the phenol protonated form. Our results also show that the high level electron correlated method is essential for a proper description of excited states of GFP. © 2017 Wiley Periodicals, Inc.