Accelerating Realtime TDDFT with Block-Orthogonalized Manby

Jul 19, 2017 - We find that by mixing levels of detail the embedded mean-field theory scheme is a simple, accurate, and effective way to accelerate RT...
0 downloads 0 Views 23MB Size
Subscriber access provided by UNIV OF YORK

Article

Accelerating Realtime TDDFT with BlockOrthogonalized Manby-Miller Embedding Theory Kevin J. Koh, Triet S. Nguyen-Beck, and John Parkhill J. Chem. Theory Comput., Just Accepted Manuscript • DOI: 10.1021/acs.jctc.7b00494 • Publication Date (Web): 19 Jul 2017 Downloaded from http://pubs.acs.org on July 20, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Chemical Theory and Computation is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Accelerating Realtime TDDFT with Block-Orthogonalized Manby-Miller Embedding Theory Kevin J. Koh, Triet S. Nguyen-Beck, and John Parkhill∗ 251 Nieuwland Science Hall, Notre Dame, IN 46556 E-mail: [email protected]

Abstract Realtime time-dependent density-functional theory (RT-TDDFT) is one of the most practical techniques available to simulate electronic dynamics of molecules and materials. Promising applications of RT-TDDFT to study non-linear spectra and energy transport, demand simulations of large solvated systems over long timescales, which are computationally quite costly. In this paper, we apply an embedding technique developed for ground-state SCF methods by Manby and Miller to accelerate realtime TDDFT. We assess the accuracy and speed of these approximations by studying the absorption spectra of solvated and covalently split chromophores. Our embedding approach is also compared with less accurate, less costly QM/MM charge embeddings. We find that by mixing levels of detail the embedded mean-field theory scheme is a simple, accurate, and effective way to accelerate RT-TDDFT simulations.

1 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1

Introduction

Realtime time-dependent density functional theory (RT-TDDFT) has become a popular choice for propagating electronic densities of medium-sized systems in response to one-body time-dependent perturbations. It can be used to simulate materials which transform light into electronic energy. 1–3 With multiple propagations, the technique can be used to simulate non-linear spectra with promising advantages over response-theory based approaches. 4–14 Although RT-TDDFT is a relatively affordable approach, it is still prohibitively expensive to routinely reach picosecond timescales, especially with functionals that feature a highaccuracy treatment of exchange. 15 In this work, we incorporate a ground-state DFT embedding approach to RT-TDDFT, showing that it is a simple and accurate way to obtain high quality results at nearly the cost of a low-quality propagation. Subsystem embedding has gained momentum as a computationally efficient strategy for modeling electronic structures of complex systems at mixed levels of detail, methods including QM/MM, 16,17 density-matrix embedding theory, 18 dynamical mean-field theory, 19 and other density-based embedding methods. 20–34 Manby and Miller have introduced several variations of embedding theory, 29,35 all of which produce effective Fock matrices for a subsystem-of-interest including effects of electronic environments at lower levels of detail. In this study, we apply the DFT embedding techniques developed by the groups of Manby and Miller to improve the efficiency of realtime TDDFT. We employ the block orthogonalized partitioning variant 29 that was introduced to produce more physical embedded density matrices. We show that when the method, hereafter referred to as the block-orthogonalized Manby-Miller embedding (BOMME), is used in place of the ordinary Fock build for RTTDDFT, explicit electronic solvation dynamics can be treated at little additional cost. The accuracy of absorption spectra produced using this multiresolution dynamics approximation is in good agreement with high-level calculations performed on the whole system and experimental spectra. The advantage of this approach is the possibility of obtaining highlevel quality simulations with nearly the speed of a low-level approximation. We present an 2 ACS Paragon Plus Environment

Page 2 of 20

Page 3 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

open-source implementation where 90% of the ideal speedup is achieved for just eight water molecules solvating a benzene. Investigating the accuracy of the technique further, we also compare with QM/MM point-charge embedding to demonstrate that the BOMME spectra are significantly closer to high-level calculations and experimental results than point-charge embedding. We examine the possibility of using the technique to mix levels of detail within a single molecule as well, and also discuss how to avoid artifacts from partitioning of the system.

2

Methods

The basic formalism of RT-TDDFT 36–39 and the BOMME technique 29 are described in previous work. We will discuss a density matrix-based realtime implementation, intended for studying non-equilibrium electronic relaxation. Realtime TDDFT and Hartree-Fock lead to a Liouville equation for the one-particle density matrix: 40,41 i γ˙ = − [Fˆ (γ), γ] h ¯

(1)

where γ is the time-dependent one-electron density matrix, and F is the time-dependent Fock matrix of HF (or KS) theory expressed in an orthonormal basis. The Fock matrix can be written as: F [γ] = h0 + G[γ] in which h0 is the core Hamiltonian, and G contains the two-electron contributions:

G[γ] = J[γ] + cx K[γ] + Vxc [γ]

(2)

J and K are the density-dependent Coulomb and exchange operators, Vxc is the KS exchangecorrelation contribution to the Fock matrix, 42 and cx is the fraction of the exact exchange for hybrid KS theories. The bottleneck of integrating Eq.1 comes from the formation of J, Vxc and especially K which is done several times a femtosecond. Techniques which decrease

3 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 20

the cost of the Fock build can easily save days of propagation time. 43–52 The BOMME method 29,35 describes a subsystem with a high-level Fock matrix and the remaining degrees of freedom with a less expensive Fock matrix by reducing the quality of basis and exchange. The block orthogonalized (BO) partitioning scheme uses a projected basis in place of the conventional atomic-orbital (AO) partitioning to define the high- and low-level components of the system. Ding and coworkers have shown that this reduces artifacts in properties related to the embedding approximation:

˜ Low [˜ H˜0 = U T H0 U , G γ ] = U T GLow [γ]U , and γ˜ = U T γU

(3)

In Eq.3, U is the transformation matrix from the nonorthogonal AO basis set to the BO basis set:

 AA

 I U=

−P

AB

IBB

0

  

In this blocked matrix, the identity matrices of subsystems A and B, IAA and IBB , have the dimensions of na and nb which are the number of basis functions in each subblock, and P AB is the projection matrix, PAB = (SAA )−1 SAB , in which SAB is the AO overlap between the subsystems. The effective Fock matrix in the BO scheme is then constructed as:

˜ Low [˜ ˜ High [˜ ˜ Low [˜ F˜ [˜ γ ] = h˜0 + G γ ] + (G γ AA ] − G γ AA ])

(4)

˜ High and G ˜ Low represent the where quantities in tildes are expressed in the BO basis. G two-electron contributions to the Fock matrix from the high- and low-level of theories, and γ˜ AA denotes the density matrix blocks that belong to the subsystem with high-level theory (hereafter called the AA block). The composition of the Fock matrix formed by BOMME is illustrated in Figure 1. We prefer the simplest scheme to calculate EEX [γ AA ], which only

4 ACS Paragon Plus Environment

Page 5 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

considers the exact exchange interaction within the AA block:

EEX0 = EEX [γ AA ] = −

1 X AA AA (µκ|υλ)γµυ γκλ 4 µυλκ∈A

(5)

Ding et al. also discussed two more complex schemes, but determined that EEX0 has the best balance between accuracy and speed, and we find it sufficient for producing accurate spectra. We have implemented RT-TDDFT propagation with BOMME method in the PySCF

F

high = H,J + Vxc +

F

= H,J + Vxc

K

High-Level Theory

Low-Level Theory

low

Approximate Fock Matrix

~ F

= H,J + Vxc +

K

BOMME

Figure 1: The Fock matrix composition in the BOMME scheme compared with its counterparts in purely high- and low-level theory. chemistry package, 53 integrating the equation of motion with a midpoint step. 34,39,54,55 Pseudocodes and a link to the GitHub repository for the code used in this work are included in the Supporting Information (SI). As demonstrated in Figure 2, the method achieves significant speedup compared with simulations of purely high-level theory, and becomes relatively affordable as more explicit solvent molecules are included in the simulation. It is important to note that the combined BOMME/RT-TDDFT scheme as described above only evolves electrons in real time, and does not propagate explicit nuclear trajectories. Consequently, the method is inappropriate for electronic dynamics that involve significant rearrangement of the nuclei. Some more costly propagations and QM/MM calculations used our implementations in the Q-Chem quantum chemistry package. 56 All propagations 5 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Calculation time ratio

1.0■ 1.0■■ Calculation time ratio

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

0.8



0.8

● ■ ◆

0.6● 0.6● 0.4

0.4 ●

0.2

0.2



◆◆

0.0 0.0 1 12





● ■ ◆





37413 hr/ps

Page 6 of 20





BOMME Manby-Miller High-level High Theory Theory Low-level Low Theory Theory





23









34

● ◆

● ◆

2134 hr/ps

45 56 67 78 # water # water molecule molecule

89

● ◆

● ◆

10 9

10

Figure 2: Speed performance (wall-hours per picosecond of simulation) with and without embedding is demonstrated on a benzene surrounded with an increasing number of water molecules. Ratios of simulation time are calculated with respect to the high-level theory scheme. The functional/basis set combinations are PBE0/6-31G* on high-level theory (red line and red shaded region), PBE/STO-3G on low-level theory (green line and green shaded region), and mixed bases on BOMME method (blue line and red/green shaded region). were performed on ground-state DFT optimized geometries. Each system experiences an impulsive electric field with a strength of 0.01 a.u. and duration of 0.07 atomic time units. In the figures provided with each spectrum, the atoms belong to high- and low-level theory regions are specified by the red and green shaded regions, respectively.

3 3.1

Results Mixed Functional Propagation

Figure 3 shows the absorption spectra obtained by performing RT-TDDFT on a system of one methanol molecule surrounded by three water molecules treated at the PBE0 level and PBE level, respectively. This simple system demonstrates the overall features of the BOMME method. Stability of dipole moments, density trace, and total entropy are shown in Figures S3, S4, and S5, respectively (SI). Note that the entropy of the high- and lowlevel subsystems are not constant, since electrons and coherence are exchanged between these regions throughout the propagation. As expected for embedding, the absorption spectrum of the whole system is a mixture of the high-quality methanol 6 ACS Paragon Plus Environment

BOMME BOMME High-level Theory High-level Theory Low-level Theory Low-level Theory Isolated Methanol Isolated Methanol

Absorption (arb. units)

2

1

0

8

2

1

0

89

9 10

10 11

11 12 13 14 12 13 14 15 Energy Energy (eV) (eV)

BOMME BOMME High-level Theory High-level Theory Low-level Theory Low-level Theory Isolated Methanol Isolated Methanol

(b)

Absorption (arb. units)

(a)

Absorption (arb. units)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

15 16

16 17

2

1

Absorption (arb. units)

Page 7 of 20

6

4

2

0 0 17 18 8 18 918

10

20 15 11 1912 13 14 EnergyEnergy (eV) (eV)

16 2117

18

22

Figure 3: Absorption spectra of one methanol with 3 water molecules: (a) Low-energy region, and (b) high-energy region. Note that the spectra are plotted at different scales for clarity. High theory employs PBE0, and low-level theory employs PBE. Both methods used a mixed basis of 6-31G*/STO-3G on methanol and water, respectively. spectrum and low-level water spectra. The BOMME spectrum visibly matches peaks of the high-level spectrum in the low-energy region (up to 18 eV), and it matches peaks of the low-level spectrum in the higher-energy region. All peaks of the BOMME method under 18 eV, other than peaks at 13.44 eV and 17.32 eV, correspond to the PBE0/6-31G* methanol peaks, whereas absorption above 18 eV and peaks at 13.44 eV and 17.32 eV correspond to PBE/STO-3G water peaks. For this particular geometry the solvation shift is small. A more strongly coupled geometry was also examined (Section 3 of SI). At this geometry, a bright singlet of the isolated methanol is at 11.85 eV, whereas the high-level theory places this state at 11.66 eV, and BOMME at 11.61 eV. No sort of unphysical collapse, instability, or artifactual transitions are observed and the dipole induced by the impulse agrees well with the high-level calculation. 29 Stepping up the system size, we examined the errors of the BOMME method that increase with the size of environment or polarizability of the solvated chromophore. We obtained the absorption spectra of a benzene with one water molecule (Figure 4.a) and ten water molecules (Figure 4.b). In Figure 4.a, the lowest π → π ∗ transition energy of BOMME, high-, and low-level theory are 7.37 eV, 7.40 eV and 8.80 eV, respectively. The good agreement is maintained in Figure 4.b where the peaks are 7.24 eV, 7.35 eV and 8.74 eV, respectively. 7 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

40

40

0

5

Absorption

30

10

0.5

0.0

20

10 8 8 10 Energy (eV) Energy (eV)

6 6

12 10

15 14

Absorption (arb. units)

10

1.0

BOMME Manby-Miller Manby-Miller High Theory High-level Low Theory Low-level

1.0

0.5

0

5 0.0 12

10 8 Energy (eV) Energy (eV)

6

10

15

12

Figure 4: Absorption spectra of benzene with (a) one and (b) ten water molecules. Highlevel theory used PBE0/6-31G*, low-level theory used PBE/STO-3G, and isolated methanol used PBE0/6-31G*. In the presence of larger number of environment, the BOMME method still performs close to the high-level theory as shown in Figure 4.b while the computing cost remains close to low-level theory (Figure 2). We also tested an intramolecular embedding, which we expected to be significantly

Absorption (arb. units)

1.5

1.0

0.5

0.0

1.0

0.5

0.0

Absorption (arb. units)

20

1.5 (b)

BOMME Manby-Miller Manby-Miller High Theory High-level Low Theory Low-level

Absorption (arb. units)

30

Absorption (arb. units)

1.5 (a) Absorption

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 20

65

1.0

Manby-Miller Manby-Miller BOMME High-level Theory High-level Theory Low-level Theory Low-level Theory

0.5

0.0 5

68 7 10 7 8 6 EnergyEnergy (eV) (eV) Energy (eV)

128

9

9

Figure 5: Absorption spectra of 3,4-Dimethylpyrazole. High-level theory used PBE0/6-31G*, and low-level theory used PBE/6-31g*. more challenging and possibly reveal false absorption peaks due to non-physical motion between high- and low-level regions. We chose 3,4-dimethylpyrazole as the first test system. Indeed, naively applying the same propagation described so far lead to spurious peaks in the BOMME method. Analyzing the time-dependent electron density, we found that, as expected, characteristic DFT over-delocalization of the density from high- to low-level regions 8 ACS Paragon Plus Environment

Page 9 of 20

caused these poles. The exciting pulse delocalized the electron and hole density into the two different regions. 57 We found a simple solution to this problem was to only perturb the high-quality block, by restricting all fields to the BO high-quality block. After a modified perturbation was applied on the system, the resulting spectra has no visible artifacts. The first significant excited-state peak is at 7.43 eV, 7.88 eV and 7.98 eV for low-level theory, the BOMME method, and high-level theory, respectively (Figure 5). Simply applying the perturbation to only the high-level region appears to effectively

0.5

0.0

1.0

1.0

0.5

0.5

Absorption (arb. units)

1.0

Absorption (arb. units)

1.5 Absorption (arb. units)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

0.0 0.0 0 60

Manby-Miller Manby-Miller BOMME High-level High-level Theory Theory Theory Low-level Low-level Theory Theory Theory

2 210 8 Energy Energy (eV) (eV) Energy (eV)

12 4

4

Figure 6: Absorption spectra of all-trans-retinal. High-level theory used PBE0 and low-level theory used PBE. All methods used mixed basis of 6-31G*/STO-3G. solve the problem of spurious embedding absorbances. Figure 6 shows the absorption spectra of all-trans-retinal. For this larger intramolecular example, BOMME also functions well when we choose a reasonable definition of high- and low-level regions that do not cleave the optically active parts of the chromophore. There is a low-level treatment of exchange energy between high- and low-level regions; therefore, when mixing exchange functionals, users should plot the electron and hole densities and verify that both densities are contained in the high-quality region. If either quasi-particle leaks into the low-quality region, the coupling between regions is physically too strong to treat accurately with BOMME, and the user should choose a different partitioning scheme. The code we have posted to GitHub is able to generate densities for this purpose.

9 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

3.2

BOMME vs. point-charge embedding

In this section, we compared the BOMME method with point-charge QM/MM embedding for a few systems, to demonstrate the value of explicit electron solvation. Figure S6 shows the low-lying absorption spectrum of benzene and ten water molecules generated by purely highlevel theory, BOMME, and QM/MM. The solvation shift of the lowest lying state in highlevel theory and BOMME were roughly -0.17 eV, whereas the QM/MM spectrum is actually blue-shifted +0.02 eV. Figure 7 shows the absorption spectra of a neutral GFP chromophore obtained using the BOMME method, QM/MM method and experiment. Point-charge embedding produces a first bright state approximately 40 nm away from the experimental peak, which is at 370 nm, whereas BOMME is in fortuitously good agreement with the experiment. In both cases, the solvation shift is captured well by BOMME, whereas point-charge embedding is totally inadequate. BOMME Experiment QM/MM Isolated GFP

1.0

Absorption (arb. units) Absorption (arb. units)

0.5 1.0

Absorption (arb. units)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

BOMME BOMME Experiment Experiment QM/MM QM/MM Isolated GFP GFP Isolated

1.0

0.0 0.5 300 0.5

350

400 Wavelength (nm)

450

500

Figure 7: Absorption spectra of a neutral GFP chromophore with ten water molecules calculated with BOMME 0.0 (6-31G/PBE0, STO-3G/PBE) and QM/MM, and of an isolated 0.0 300 300 350 350 400 400 450 450 500 500 GFP neutral chromophore (PBE0/6-31G),Wavelength compared with experimental spectrum from Ref. (nm) (nm) Wavelength 58.

4

Conclusions

In this work, we investigate the usefulness of BOMME applied to RT-TDDFT. The BO variation of the BOMME scheme is found to produce spectra that are in good agreement 10 ACS Paragon Plus Environment

Page 10 of 20

Page 11 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

with high-level calculations and experimental data at a cost much closer to the low-level calculations. We have shown that the method works as both intermolecular and intramolecular embedding schemes, and usefully captures the shift of a solvated chromophore. However, the case of intramolecular embedding requires careful attention to the way density moves in the propagation, and to verify that the electron and hole do not significantly leak into the low-quality region. It is also advantageous to only perturb the high-quality region of the spectrum, to avoid driving electrons non-physically across the boundary between highand low-quality regions. Because of its simplicity, the method can easily be combined with a multiple timestepping scheme to obtain further speedups; 59,60 we may explore these directions in future work. It can also be combined with our methods for electronic dissipation to study electronic relaxation phenomena and non-linear spectra.

5

Associated Content

Pseudocode and schematic for BOMME, additional plots for a system of one methanol and three water molecules, and spectra for a system of one benzene and ten water molecules. This material is available free of charge via the Internet at http://pubs.acs.org.

6

Author Information

Corresponding author: John Parkhill, E-mail: [email protected]

7

Notes

The authors declare no competing financial interest.

11 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Acknowledgement We thank the University of Notre Dame’s College of Science and Department of Chemistry and Biochemistry for generous start-up funding. T. S. N.-B. acknowledges the support of the NSF Graduate Research Fellowship under grant DGE-1313583 and the Patrick and Jana Eilers Graduate Student Fellowship for Energy Related Research.

References (1) Lopata, K.; Govind, N. Modeling Fast Electron Dynamics with Real-Time TimeDependent Density Functional Theory: Application to Small Molecules and Chromophores. J. Chem. Theory Comput 2011, 7, 1344–1355. (2) Tussupbayev, S.; Govind, N.; Lopata, K.; Cramer, C. J. Comparison of Real-Time and Linear-Response Time-Dependent Density Functional Theories for Molecular Chromophores Ranging from Sparse to High Densities of States. J. Chem. Theory Comput 2015, 11, 1102–1109. (3) Peng, B.; Lingerfelt, D. B.; Ding, F.; Aikens, C. M.; Li, X. Real-Time TDDFT Studies of Exciton Decay and Transfer in Silver Nanowire Arrays. J. Phys. Chem. C 2015, 119, 6421–6427. (4) Fischer, S. A.; Cramer, C. J.; Govind, N. Excited State Absorption from Real-Time Time-Dependent Density Functional Theory. J. Chem. Theory Comput. 2015, 11, 4294–4303. (5) Nguyen, T. S.; Koh, J. H.; Lefelhocz, S.; Parkhill, J. Black-Box, Real-Time Simulations of Transient Absorption Spectroscopy. J. Phys. Chem. Lett. 2016, 7, 1590–1595. (6) Nguyen, T. S.; Parkhill, J. Nonradiative Relaxation in Real-Time Electronic Dynamics OSCF2: Organolead Triiodide Perovskite. J. Phys. Chem. A 2016, 120, 6880–6887. 12 ACS Paragon Plus Environment

Page 12 of 20

Page 13 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

(7) Isborn, C. M.; Li, X. Modeling the Doubly Excited State with Time-Dependent Hartree–Fock and Density Functional Theories. J. Chem. Phys. 2008, 129, 204107. (8) Takimoto, Y.; Isborn, C.; Eichinger, B.; Rehr, J.; Robinson, B. Frequency and Solvent Dependence of Nonlinear Optical Properties of Molecules. J. Phys. Chem. C 2008, 112, 8016–8021. (9) Isborn, C. M.; Li, X. Singlet- Triplet Transitions in Real-Time Time-Dependent Hartree- Fock/Density Functional Theory. J. Chem. Theory Comput 2009, 5, 2415– 2419. (10) Habenicht, B. F.; Tani, N. P.; Provorse, M. R.; Isborn, C. M. Two-electron Rabi Oscillations in Real-Time Time-Dependent Density-Functional Theory. J. Chem. Phys. 2014, 141, 184112. (11) Provorse, M. R.; Isborn, C. M. Electron Dynamics with Real-Time Time-Dependent Density Functional Theory. Int. J. Quantum Chem. 2016, 116, 739–749. (12) Lopata, K.; Govind, N. Near and Above Ionization Electronic Excitations with NonHermitian Real-Time Time-Dependent Density Functional Theory. J. Chem. Theory Comput 2013, 9, 4939–4946. (13) DePrince III, A. E.; Pelton, M.; Guest, J. R.; Gray, S. K. Emergence of Excited-State Plasmon Modes in Linear Hydrogen Chains from Time-Dependent Quantum Mechanical Methods. Phys. Rev. Lett. 2011, 107, 196806. (14) Nascimento, D. R.; DePrince, A. E. Linear Absorption Spectra from Explicitly TimeDependent Equation-of-Motion Coupled-Cluster Theory. J. Chem. Theory Comput 2016, 12, 5834–5840. (15) Cui, G.; Fang, W.; Yang, W. Reformulating Time-Dependent Density Functional The-

13 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ory with Non-Orthogonal Localized Molecular Orbitals. Phys. Chem. Chem. Phys. 2010, 12, 416–421. (16) Warshel, A.; Levitt, M. Theoretical Studies of Enzymic Reactions: Dielectric, Electrostatic and Steric Stabilization of the Carbonium Ion in the Reaction of Lysozyme. J. Mol. Biol. 1976, 103, 227–249. (17) Senn, H. M.; Thiel, W. QM/MM Methods for Biomolecular Systems. Angew. Chem., Int. Ed. 2009, 48, 1198–1229. (18) Knizia, G.; Chan, G. K.-L. Density Matrix Embedding: A Simple Alternative to Dynamical Mean-Field Theory. Phys. Rev. Lett. 2012, 109, 186404. (19) Georges, A.; Kotliar, G.; Krauth, W.; Rozenberg, M. J. Dynamical Mean-Field Theory of Strongly Correlated Fermion Systems and the Limit of Infinite Dimensions. Rev. Mod. Phys. 1996, 68, 13. (20) Govind, N.; Wang, Y. A.; Carter, E. A. Electronic-Structure Calculations by FirstPrinciples Density-Based Embedding of Explicitly Correlated Systems. J. Chem. Phys. 1999, 110, 7677–7688. (21) Neugebauer, J.; Jacob, C. R.; Wesolowski, T. A.; Baerends, E. J. An Explicit Quantum Chemical Method for Modeling Large Solvation Shells Applied to Aminocoumarin C151. J. Phys. Chem. A 2005, 109, 7805–7814. (22) Goodpaster, J. D.; Ananth, N.; Manby, F. R.; Miller III, T. F. Exact Nonadditive Kinetic Potentials for Embedded Density Functional Theory. J. Chem. Phys. 2010, 133, 084103. (23) Goodpaster, J. D.; Barnes, T. A.; Miller III, T. F. Embedded Density Functional Theory for Covalently Bonded and Strongly Interacting Subsystems. J. Chem. Phys. 2011, 134, 164108. 14 ACS Paragon Plus Environment

Page 14 of 20

Page 15 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

(24) Manby, F. R.; Stella, M.; Goodpaster, J. D.; Miller III, T. F. A Simple, Exact DensityFunctional-Theory Embedding Scheme. J. Chem. Theory Comput. 2012, 8, 2564. (25) Barnes, T. A.; Goodpaster, J. D.; Manby, F. R.; Miller III, T. F. Accurate Basis Set Truncation for Wavefunction Embedding. J. Chem. Phys. 2013, 139, 024103. (26) Goodpaster, J. D.; Barnes, T. A.; Manby, F. R.; Miller III, T. F. Accurate and Systematically Improvable Density Functional Theory Embedding for Correlated Wavefunctions. J. Chem. Phys. 2014, 140, 18A507. (27) Bennie, S. J.; Stella, M.; Miller III, T. F.; Manby, F. R. Accelerating Wavefunction in Density-Functional-Theory Embedding by Truncating the Active Basis Set. J. Chem. Phys. 2015, 143, 024105. (28) Wesolowski, T. A.; Shedge, S.; Zhou, X. Frozen-Density Embedding Strategy for Multilevel Simulations of Electronic Structure. Chem. Rev. 2015, 115, 5891–5928. (29) Ding, F.; Manby, F. R.; Miller III, T. F. Embedded Mean-Field Theory with BlockOrthogonalized Partitioning. J. Chem. Theory Comput. 2017, 13, 1605–1615. (30) Ramos, P.; Pavanello, M. Performance of Frozen Density Embedding for Modeling Hole Transfer Reactions. J. Phys. Chem. B 2015, 119, 7541–7557. (31) Pavanello, M.; Neugebauer, J. Modelling Charge Transfer Reactions with the Frozen Density Embedding Formalism. J. Chem. Phys. 2011, 135, 234103. (32) Kretchmer, J. S.; Chan, G. K. A Real-Time Extension of Density Matrix Embedding Theory for Non-Equilibrium Electron Dynamics. arXiv Preprint 2016, arXiv:1609.07678. (33) Zheng, B.-X.; Kretchmer, J. S.; Shi, H.; Zhang, S.; Chan, G. K.-L. Cluster Size Convergence of the Density Matrix Embedding Theory and Its Dynamical Cluster Formulation:

15 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

A Study with an Auxiliary-Field Quantum Monte Carlo Solver. Phys. Rev. B 2017, 95, 045103. (34) Li, X.; Tully, J. C.; Schlegel, H. B.; Frisch, M. J. Ab Initio Ehrenfest Dynamics. J. Chem. Phys. 2005, 123, 084106. (35) Fornace, M. E.; Lee, J.; Miyamoto, K.; Manby, F. R.; Miller III, T. F. Embedded Mean-Field Theory. J. Chem. Theory Comput. 2015, 11, 568–580. (36) McLachlan, A. D.; Ball, M. A. Time-Dependent Hartree-Fock Theory for Molecules. Rev. Mod. Phys. 1964, 36, 844–855. (37) Casida, M. E. Time-Dependent Density-Functional Theory for Molecules and Molecular Solids. J. Mol. Struct.: THEOCHEM 2009, 914, 3 – 18. (38) Furche, F.; Ahlrichs, R. Adiabatic Time-Dependent Density Functional Methods for Excited State Properties. J. Chem. Phys. 2002, 117, 7433–7447. (39) Nguyen, T. S.; Parkhill, J. Nonadiabatic Dynamics for Electrons at Second-Order: Real-Time TDDFT and OSCF2. J. Chem. Theory Comput. 2015, 11, 2918–2924. (40) Kohn, W.; Sham, L. J. Self-Consistent Equations including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133. (41) Parr, R. G. Horizons of Quantum Chemistry; Springer, 1980; pp 5–15. (42) Pople, J. A.; Gill, P. M.; Johnson, B. G. Kohn-Sham Density-Functional Theory within a Finite Basis Set. Chem. Phys. Lett. 1992, 199, 557–560. (43) Alml¨of, J.; Faegri, K.; Korsell, K. Principles for a Direct SCF Approach to LICAO-MO ab-Initio Calculations. J. Comput. Chem. 1982, 3, 385–399. (44) H¨aser, M.; Ahlrichs, R. Improvements on the Direct SCF Method. J. Comput. Chem. 1989, 10, 104–111. 16 ACS Paragon Plus Environment

Page 16 of 20

Page 17 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

(45) White, C. A.; Head-Gordon, M. Derivation and Efficient Implementation of the Fast Multipole Method. J. Chem. Phys. 1994, 101, 6593–6605. (46) Strain, M. C.; Scuseria, G. E.; Frisch, M. J. Achieving Linear Scaling for the Electronic Quantum Coulomb Problem. Science 1996, 271, 51. (47) Challacombe, M.; Schwegler, E.; Alml¨of, J. Fast Assembly of the Coulomb Matrix: A Quantum Chemical Tree Code. J. Chem. Phys. 1996, 104, 4685–4698. (48) Schwegler, E.; Challacombe, M. Linear Scaling Computation of the Hartree–Fock Exchange Matrix. J. Chem. Phys. 1996, 105, 2726–2734. (49) Burant, J. C.; Scuseria, G. E.; Frisch, M. J. A Linear Scaling Method for Hartree-Fock Exchange Calculations of Large Molecules. J. Chem. Phys. 1996, 105, 8969–8972. (50) Izmaylov, A. F.; Scuseria, G. E.; Frisch, M. J. Efficient Evaluation of Short-Range Hartree-Fock Exchange in Large Molecules and Periodic Systems. J. Chem. Phys. 2006, 125, 104103. (51) Rudberg, E.; Salek, P. Efficient Implementation of the Fast Multipole Method. J. Chem. Phys. 2006, 125, 084106. (52) Liang, W.; Chapman, C. T.; Li, X. Efficient First-Principles Electronic Dynamics. J. Chem. Phys. 2011, 134, 184102. (53) Sun, Q.; Berkelbach, T. C.; Blunt, N. S.; Booth, G. H.; Guo, S.; Li, Z.; Liu, J.; McClain, J.; Sharma, S.; Wouters, S.; Chan, G. K.-L. The Python-Based Simulations of Chemistry Framework (PySCF). arXiv Preprint 2017, arXiv:1701.08223. (54) Micha, D. A. Time-Dependent Many-Electron Treatment of Electronic Energy and Charge Transfer in Atomic Collisions. J. Phys. Chem. A 1999, 103, 7562–7574. (55) Castro, A.; Marques, M. A. L.; Rubio, A. Propagators for the Time-Dependent KohnSham Equations. J. Chem. Phys. 2004, 121, 3425–3433. 17 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(56) Shao, Y.; Gan, Z.; Epifanovsky, E.; Gilbert, A. T.; Wormit, M.; Kussmann, J.; Lange, A. W.; Behn, A.; Deng, J.; Feng, X.; Ghosh, D.; Goldey, M.; Horn, P. R.; Jacobson, L. D.; Kaliman, I.; Khaliullin, R. Z.; Ku´s, T.; Landau, A.; Liu, J.; Proynov, E. I.; Rhee, Y. M.; Richard, R. M.; Rohrdanz, M. A.; Steele, R. P.; Sundstrom, E. J.; III, H. L. W.; Zimmerman, P. M.; Zuev, D.; Albrecht, B.; Alguire, E.; Austin, B.; Beran, G. J. O.; Bernard, Y. A.; Berquist, E.; Brandhorst, K.; Bravaya, K. B.; Brown, S. T.; Casanova, D.; Chang, C.-M.; Chen, Y.; Chien, S. H.; Closser, K. D.; Crittenden, D. L.; Diedenhofen, M.; Jr., R. A. D.; Do, H.; Dutoi, A. D.; Edgar, R. G.; Fatehi, S.; Fusti-Molnar, L.; Ghysels, A.; Golubeva-Zadorozhnaya, A.; Gomes, J.; Hanson-Heine, M. W.; Harbach, P. H.; Hauser, A. W.; Hohenstein, E. G.; Holden, Z. C.; Jagau, T.-C.; Ji, H.; Kaduk, B.; Khistyaev, K.; Kim, J.; Kim, J.; King, R. A.; Klunzinger, P.; Kosenkov, D.; Kowalczyk, T.; Krauter, C. M.; Lao, K. U.; Laurent, A. D.; Lawler, K. V.; Levchenko, S. V.; Lin, C. Y.; Liu, F.; Livshits, E.; Lochan, R. C.; Luenser, A.; Manohar, P.; Manzer, S. F.; Mao, S.-P.; Mardirossian, N.; Marenich, A. V.; Maurer, S. A.; Mayhall, N. J.; Neuscamman, E.; Oana, C. M.; Olivares-Amaya, R.; O’Neill, D. P.; Parkhill, J. A.; Perrine, T. M.; Peverati, R.; Prociuk, A.; Rehn, D. R.; Rosta, E.; Russ, N. J.; Sharada, S. M.; Sharma, S.; Small, D. W.; Sodt, A.; Stein, T.; St¨ uck, D.; Su, Y.-C.; Thom, A. J.; Tsuchimochi, T.; Vanovschi, V.; Vogt, L.; Vydrov, O.; Wang, T.; Watson, M. A.; Wenzel, J.; White, A.; Williams, C. F.; Yang, J.; Yeganeh, S.; Yost, S. R.; You, Z.-Q.; Zhang, I. Y.; Zhang, X.; Zhao, Y.; Brooks, B. R.; Chan, G. K.; Chipman, D. M.; Cramer, C. J.; III, W. A. G.; Gordon, M. S.; Hehre, W. J.; Klamt, A.; III, H. F. S.; Schmidt, M. W.; Sherrill, C. D.; Truhlar, D. G.; Warshel, A.; Xu, X.; Aspuru-Guzik, A.; Baer, R.; Bell, A. T.; Besley, N. A.; Chai, J.-D.; Dreuw, A.; Dunietz, B. D.; Furlani, T. R.; Gwaltney, S. R.; Hsu, C.-P.; Jung, Y.; Kong, J.; Lambrecht, D. S.; Liang, W.; Ochsenfeld, C.; Rassolov, V. A.; Slipchenko, L. V.; Subotnik, J. E.; Voorhis, T. V.; Herbert, J. M.; Krylov, A. I.; Gill, P. M.; Head-Gordon, M. Advances in Molecular Quantum Chemistry Contained

18 ACS Paragon Plus Environment

Page 18 of 20

Page 19 of 20

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

in the Q-Chem 4 Program Package. Mol. Phys. 2015, 113, 184–215. (57) Milanese, J. M.; Provorse, M. R.; Alameda, E.; Isborn, C. M. Convergence of Computed Aqueous Absorption Spectra with Explicit Quantum Mechanical Solvent. J. Chem. Theory Comput. 2017, (58) Nielsen, S. B.; Lapierre, A.; Andersen, J. U.; Pedersen, U. V.; Tomita, S.; Andersen, L. H. Absorption Spectrum of the Green Fluorescent Protein Chromophore Anion in Vacuo. Phys. Rev. Lett. 2001, 87, 228102. (59) Tuckerman, M.; Berne, B. J.; Martyna, G. J. Reversible Multiple Time Scale Molecular Dynamics. J. Chem. Phys. 1992, 97, 1990–2001. (60) McClean, J. R.; Parkhill, J. A.; Aspuru-Guzik, A. Feynmans Clock, a New Variational Principle, and Parallel-in-Time Quantum Dynamics. Proc. Natl. Acad. Sci. U.S.A. 2013, 110, E3901–E3909.

19 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

8

For Table of Contents use only RT-TDDFT with Manby-Miller Embedding 1.0

0.5 0.5

0.5

Absorption (arb. units)

1.0 1.0

Absorption (arb. units) Absorption (arb. units)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

0.0 300 0.0 0.0 300 300

BOMME Experiment BOMME BOMME QM/MM Experiment Experiment Isolated GFP QM/MM QM/MM Isolated GFP GFP Isolated

350 400 450 350 350 400 (nm) 450 450 400 Wavelength Wavelength (nm) (nm) Wavelength

~ F

H

J,Vxc

K

500 500 500

Figure 8: The embedded mean-field theory scheme and the absorption spectra of a GFP neutral chromophore solvated in water computed with BOMME/RT-TDDFT and QM/MM.

20 ACS Paragon Plus Environment

Page 20 of 20