Modeling Optical Spectra of Large Organic Systems Using Real-Time

Aug 16, 2017 - We present an implementation of a time-dependent semiempirical method (INDO/S) in NWChem using real-time (RT) propagation to address, i...
0 downloads 10 Views 1MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Modeling Optical Spectra of Large Organic Systems using RealTime Propagation of Semiempirical Effective Hamiltonians Soumen Ghosh, Amity Andersen, Laura Gagliardi, Christopher J. Cramer, and Niranjan Govind J. Chem. Theory Comput., Just Accepted Manuscript • DOI: 10.1021/acs.jctc.7b00618 • Publication Date (Web): 16 Aug 2017 Downloaded from http://pubs.acs.org on August 22, 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 23

Journal of Chemical Theory and Computation

Draft of August 12, 2017 Prepared for Journal of Chemical Theory and Computation 1 2 3

Modeling Optical Spectra of Large Organic Systems using Real-Time Propagation of Semiempirical Effective Hamiltonians 7

6

5

4

Soumen Ghosh†, Amity Andersen§, Laura Gagliardi*†, Christopher J. Cramer*†, and Niranjan Govind*§ 10

9

8



1

Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA. 14

13

12

§

Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99338, USA 16

15 17 19

18

Abstract: 21

20

We present an implementation of a time-dependent semiempirical method (INDO/S) in NWChem using real-time (RT) propagation to address, in principle, the entire spectrum of valence electronic excitations. Adopting this model, we study the UV/Vis spectra of medium-sized systems like P3B2, f-coronene, and in addition much larger systems like ubiquitin in the gas phase and the betanin chromophore in the presence of two explicit solvents (water and methanol). RT-INDO/S provides qualitatively and often quantitatively accurate results when compared with RT- TDDFT or experimental spectra. Even though we only consider the INDO/S Hamiltonian in this paper, our implementation provides a framework for performing electron dynamics in large systems using semiempirical Hartree-Fock (HF) Hamiltonians in general. 3

32

31

30

29

28

27

26

25

24

23

2

Introduction 35

34

In the last decade, performing large-scale static and dynamic calculations has become essential in different areas of chemical research,1 including, e.g., modeling the explicit solvation of chromophores,2-5 the active sites of enzymes,6,7 homogeneous clusters,8-11 proteins,12-16 photochemistry,17,18 and the properties of crystals,19,20 polymers,21-23 and nanomaterials.24,25. Performing ab initio electronic structure calculation on large systems has become significantly more practical, thanks to advances in computational hardware and new efficient algorithms. Density functional theory (DFT) calculations are now regularly performed for systems of up to O(103) atoms (using local functionals), although more stringent limitations are associated with post-Hartree-Fock (HF) methods that take increasing account of electron correlation. Even in favorable instances, however, such large-scale calculations tend to be restricted to molecular (or supermolecular) geometries that are static. In many instances, it is important not only to model interactions of the system with its environment but also its dynamical properties. When explicit solvation is included, many relevant chemical systems require from 103 to 104 atoms for a physical treatment, and such systems are not yet trivially accessible with ab initio DFT and post-HF methods. 5

54

53

52

51

50

49

48

47

46

45

4

43

42

41

40

39

38

37

36

Combined quantum mechanical/molecular mechanical (QM/MM) methods have often proven effective for the study of large systems,26 but the accuracy of QM/MM methods is strongly 57

56 58 59 60

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 2 of 23

2 1 3

2

dependent on the chosen size of the QM region,5,27-29 the accuracy of the electronic structure methods used, and the minimization of boundary effects.30 Treating an entire system at a single level of theory avoids such complications, and approximate electronic structure methods employing effective Hamiltonians have long been the subject of extensive study owing to their computational efficiency. Such methods were first developed by neglecting selected electronelectron repulsion integrals, and computing the remaining terms in the Hamiltonian by correlating them with different physical properties of relevant chemical systems. The first semiempirical methods to appear in the literature were restricted to π orbitals in organic systems, e.g., the Hückel31 and Pariser-Parr-Pople (PPP)32-34 methods. Subsequently, Hoffmann and co-workers extended the Hückel method to describe both π and σ electronic bonding.35 Pople and coworkers then developed all-valence-electron semiempircal methods based on HF theory, including three semiempirical approximations: complete neglect of differential overlap (CNDO),36,37 intermediate neglect of differential overlap (INDO)38 and neglect of diatomic differential overlap (NDDO)39. Among these approximate methods, INDO received special attention from Zerner and co-workers, who parametrized the INDO (INDO/S) method specifically to compute accurately the UV/Vis spectra of organic molecules at the configuration interaction including single excitations (CIS) as well as within the random phase approximation.40-44 INDO/S parameters are now available for transition metals,45,46 actinides,47 and lanthanides48. Modifications to the INDO/S method, INDO/X, have also been reported.49 An alternative effective Hamiltonian model formally derived from density functional theory, namely the density functional tight binding (DFTB) model,50 has also been extensively developed by a number of groups51-57. 29

28

27

26

25

24

23

2

21

20

19

18

17

16

15

14

13

12

1

10

9

8

7

6

5

4

30

Considering electron dynamics, several quantum mechanical approaches are available. Within the Born-Oppenheimer (BO) approximation, these include, e.g., direct integration of the Schrodinger equation for very small systems (e.g., H2)58, real-time configuration interaction singles (CIS),59,60 real-time time-dependent density functional theory (RT-TDDFT),61-65 and real-time timedependent Hartree-Fock (RT-TDHF)66,67. In the framework of semiempirical methods, wave packet dynamics with an extended Hückel Hamiltonian68 has been applied to gas phase systems,69 molecules in solution,70 and dye-sensitized semiconductor interfaces71-76 (recently, this latter model has been combined with nuclear dynamics as well77,78). Real-time dynamics approaches formulated for DFTB have also been used to study electron dynamics79 as have analogs formulated for semiempirical HF methods,80-83 with both applied to the study of both ground and excited states in small- to medium-sized systems. In the specific case of DFTB, interesting recent applications include the study of non-adiabatic electron-hole dynamics in nanoscale solar cell materials.84-86 45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

46

Regarding the prediction of UV/Vis spectra, the CIS method has been extensively explored with semiempirical methods,87,88 with the INDO/S Hamiltonian having been used in particular to study excited states in organic electronics,89-91 conjugated polymers,92-95 biological molecules,96-99 and nanoparticles.100,101 Other approaches applied to study very large systems include combining CIS with PPP to study the UV/Vis spectra of organic dyes102,103 and combining linear-response timedependent (LR-TD) theory with DFTB104,105. Since traditional implementations of CIS and LRTD approaches are typically “bottom-up” approaches, they are mostly used to study the lowest excited states and can become prohibitively expensive for large systems and systems with high densities of states. Real-time approaches mitigate this issue as they can be used to capture broad 58

57

56

5

54

53

52

51

50

49

48

47

59 60 ACS Paragon Plus Environment

Page 3 of 23

Journal of Chemical Theory and Computation

3 1 3

2

spectral regions depending on the size of the time step and the length of the simulation. 65,106-108 More recently real-time approaches have also been used to study excited-state absorption.109-111 In this context, we should also mention other viable approaches that can address wide absorption spectra, the complex polarization,112-114 damped response approaches115,116 and multishift linear solvers117. Other approaches include the self-consistent constricted variational DFT approach by Ziegler and co-workers118, the simplified Tamm-Dancoff119 and simplified TDDFT reported by Grimme and co-workers120. and the recent symmetric Lanczos algorithm and kernel polynomial method (KPM) due to Brabec and co-workers121. 14

13

12

1

10

9

8

7

6

5

4

In this paper, we present an implementation of RT-INDO/S, in the NWChem,122 computational chemistry program for the modeling UV/Vis spectra of large molecular complexes. We first outline the general methodology for real-time density matrix propagation with a semiempirical Hamiltonian. Next, we apply the approach to medium size molecular complexes like P3B2 (C72H38N20), and f-coronene (C108H42N12) in the gas phase. Finally, to highlight our approach, we calculate the UV/Vis absorption spectra of tyrosine isolated in the gas phase and present in the protein ubiquitin, and also of the betanin red pigment molecular complex in explicit aqueous and alcohol (methanol) solutions. We show that, without any modification of the existing parameters, RT-INDO/S provides qualitatively and often quantitatively accurate results when compared with experiment or RT- TDDFT results and it is thus applicable to a large class of systems. 26

25

24

23

2

21

20

19

18

17

16

15

27 29

28

Theory: 31

30

Before discussing the INDO/S approximation, we briefly discuss ab-initio HF theory to highlight changes in the equations due to the approximations. 3

32 34

Ab-initio Hartree-Fock (HF) Theory: 36

35

In Hartree-Fock theory123, for a closed-shell system, the molecular orbitals (MO), 𝜓𝑖 are considered as linear combination of atomic orbitals (LCAO) 𝜙𝜇 , 39

38

37

40

𝜓𝑖 = ∑ 𝐶𝑖𝜇 𝜙𝜇 42

𝜇

43

41

(1)

where 𝐶𝑖𝜇 is the linear expansion coefficient. Molecular orbitals are optimized by solving the matrix equation, 46

45

4

47

𝐹𝐶 = 𝑆𝐶𝐸 49

48

(2)

where F and S are the Fock and overlap matrices, C is the matrix of linear expansion coefficients and E is the diagonal matrix containing the MO eigenvalues, respectively. Matrix elements of S are given by, 52

51

50

53 54

𝑆𝜇𝜐 = ∫ 𝜙𝜇∗ (𝑟1 )𝜙𝜐 (𝑟1 )𝑑𝑟1 5

(3)

56

and those of F are given by, 57 58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 4 of 23

4 1 3

2

1 𝐹𝜇𝜐 = 𝐻𝜇𝜐 + ∑ 𝑃𝜆𝜎 [(𝜇𝜐|𝜆𝜎) − (𝜇𝜆|𝜐𝜎)] 2 5

4

(4)

𝜆𝜎

7

6

where the one-electron contribution is given by 8

atoms

9

𝐻𝜇𝜐 = 1

10

∫ 𝜙𝜇∗ (𝑟1 ) [−

1 2 𝑍𝐴 ∇ + ∑ ] 𝜙 (𝑟 )𝑑𝑟 (5) 2 𝑟1 − 𝑅𝐴 𝜐 1 1 𝐴

13

12

where ZA is the nuclear charge of atom A, 𝑟1 is the electronic position and RA is the atomic position and the two-electron electron-repulsion integrals are computed as 15

14 16

(𝜇𝜐|𝜆𝜎) = ∬ 𝜙𝜇∗ (𝑟1 )𝜙𝜐 (𝑟1 ) 19

18

17

1 ∗ 𝜙 (𝑟 )𝜙 (𝑟 )𝑑𝑟 𝑑𝑟 𝑟12 𝜆 2 𝜎 2 1 2

(6)

where 𝑟12 is the interelectronic distance. The one electron reduced density matrix is given by 20 21

𝑀𝑂

2

∗ 𝑃𝜆𝜎 = ∑ 𝐶𝜇𝑖 𝐶𝑖𝜎

23 24

(7)

𝑖

26

25

where the index i runs over only the occupied molecular orbitals. 27 28 30

29

INDO/S Approximation: 31

The INDO/S method is based on the INDO/1 method38 and can be summarized as follows: 3

32

(i) Overlap integrals (eq. 3) are neglected unless 𝜇 = 𝜐. Given S = 1, this reduces Eq. 2 to 35

34

𝐹𝐶 = 𝐶𝐸 37

36

(8)

(ii) Integrals of type (𝜇𝜐|𝜆𝜎) are set equal to zero unless 𝜇 = 𝜐 and 𝜆 = 𝜎, in which case, if the two pairs of basis functions are located on different atoms 40

39

38

(𝜇𝜐|𝜆𝜎) = 𝛿𝜇𝜐 𝛿𝜆𝜎 𝛾𝐴𝐵 (μ on atom A and  on atom B) (9) 43

42

41

where the Mataga-Nishimoto-Weiss formula40 takes 45

4

𝛾𝐴𝐵 = 48

47

46

𝑓𝛾 𝑅𝐴𝐵 +

2𝑓𝛾 (𝛾𝐴𝐴 +𝛾𝐵𝐵)

(10)

and 𝛾𝐴𝐴 = 𝐼𝐴 − 𝐴𝐴 , where IA and AA are the ionization potential and electron affinity of atom A,40,41 respectively, and Weiss suggested a value of 1.2 bohrs∙hartree for 𝑓𝛾 parameter. 𝑅𝐴𝐵 is the distance between atom A and B. Single-center electron-repulsion integrals in the INDO model are expressed through Slater-Condon parameters40 54

53

52

51

50

49

(𝑠𝑠|𝑠𝑠) = (𝑠𝑠|𝑝𝑝) = 𝐹 0 5 56 57 58 59 60 ACS Paragon Plus Environment

Page 5 of 23

Journal of Chemical Theory and Computation

5 1 2 3

1 (𝑠𝑝|𝑠𝑝) = 𝐺 1 (𝑠𝑝) 3 4 5 6

(𝑝𝑥 𝑝𝑥 |𝑝𝑥 𝑝𝑥 ) = 𝐹 0 (𝑝𝑝) + 7 8 9

4 25

𝐹 2 (𝑝𝑝) (11)

(𝑝𝑥 𝑝𝑥 |𝑝𝑦 𝑝𝑦 ) = 𝐹 0 (𝑝𝑝) − 10 1

2 2 𝐹 (𝑝𝑝) 25

3

(𝑝𝑥 𝑝𝑦 |𝑝𝑥 𝑝𝑦 ) = 25 𝐹 2 (𝑝𝑝) 12 13 14 16

15

(iii) Diagonal one-electron matrix elements for atomic orbitals ϕμ centered on atom A are computed as 18

17 19

𝑐𝑜𝑟𝑒 𝐻𝜇𝜇 = 𝑈𝜇𝜇 − ∑ 𝑍𝐵 𝛾𝐴𝐵

20 21

(12)

𝐵≠𝐴

23

2

where 𝑈𝜇𝜇 is the one-center core integrals. 25

24

(iv) Two-center off-diagonal one-electron matrix elements are calculated using the weighted overlap approach 27

26 28

𝑐𝑜𝑟𝑒 𝐻𝜇𝜐 = 𝑓𝜇𝜐 (𝛽𝜇0 + 𝛽𝜐0 )𝑆𝜇𝜐

30

29

(13)

where the parameter fμυ is equal to 1.267 and 0.585 for σ and π overlap41, respectively, 𝑆𝜇𝜐 is computed using Slater orbitals (i.e., atomic orbitals are not assumed to be orthonormal for this step), and a single βA is used for each element and semiempirically adjusted to yield best agreement with reference spectroscopic data. 35

34

3

32

31

36

Using the above approximations, F can be defined for closed-shell system as, 38

1 𝐹𝜇𝜇 = 𝑈𝜇𝜇 − ∑ 𝑍𝐵 𝛾𝐴𝐵 + ∑ 𝑃𝜐𝜐 [𝛾𝐴𝐴 − (𝜇𝜐|𝜇𝜐)] + ∑ ∑ 𝑃𝜆𝜆 𝛾𝐴𝐵 2 40

39

37

𝐵≠𝐴

42

41

𝜐

𝐵≠𝐴 𝜆∈𝐵

1

𝐹𝜇𝜐 = 2 𝑃𝜇𝜐 (3(𝜇𝜐|𝜇𝜐) − 𝛾𝐴𝐴 ) (𝜇, 𝜐 ∈ 𝐴) 4

43

(𝜇, 𝜐 ∈ 𝐴, 𝜆 ∈ 𝐵)

(14)

𝐹𝜇𝜆 = 𝑓𝜇𝜐 (𝛽𝜇0 + 𝛽𝜐0 )𝑆𝜇𝜐 − 𝑃𝜇𝜆 𝛾𝐴𝐵 (𝜇 ∈ 𝐴, 𝜆 ∈ 𝐵 ) 47

46

45

The INDO/S approximation uses a Slater-type minimal basis set. All parameters used in this paper are reported in the supporting information (Table S1). Equation 14 can be easily extended for openshell systems and our implementation can perform such calculations. However, all the systems reported in this paper are closed-shell. 51

50

49

48

53

52

Time-Dependent Fock Matrix: Unlike ground state DFT calculations where the Fock and density matrices are purely real (with real basis functions), both become complex quantities when propagated in time as we shall see shortly. Using Eq. 4, the time-dependent Fock matrix in the presence of an applied time-dependent field is given by, 57

56

5

54

58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 6 of 23

6 1 2 3

𝑐𝑜𝑟𝑒 𝐹𝜇𝜐 (𝑡) = 𝐻𝜇𝜐 + 𝐺𝜇𝜐 [𝑃(𝑡)] − 𝐷𝜇𝜐 𝐸(𝑡)

5

4

(15)

where 𝐺𝜇𝜐 is the two-electron component of the Fock matrix element, 𝐷𝜇𝜐 is an element of the dipole matrix, 7

6 8

𝑥 𝐷𝜗𝜇 = ∫ 𝜑𝜇∗ (𝑟)𝑥𝜑𝜗 (𝑟) 𝑑𝑟

9 1

10

(16)

and E is the external electric field. 13

12

For a delta function electric field pulse, 14

𝐸(𝑡) = 𝛿(𝑡)𝐸0 16

15

(17)

where 𝐸0 is the strength of the delta electric field. Under the influence of the weak perturbation, the difference between the instantaneous (𝑃(𝑡)) and initial density (𝑃0 ) matrices can be written as, 19

18

17

𝑊(𝑡) = 𝑃(𝑡) − 𝑃0 = −𝑖[𝐸 ∙ 𝐷, 𝑃0 ] (18) 20 2

21

Time Propagation: Time propagation of the density matrix is governed by the von Neumann equation, 24

23 25

𝑖 26 28

27

𝜕𝑃 = [𝐹(𝑡), 𝑃(𝑡)] 𝜕𝑡

(19)

During time propagation, the density matrix is in general complex. In the case of the Fock matrix, it can be shown that the contribution to the imaginary part comes exclusively from the exact exchange part of the Fock matrix.65 32

31

30

29

In semiempirical HF methods, the commutator [𝐹0 , 𝑃0 ] (where 𝐹0 and 𝑃0 are the Fock and density matrices before perturbation, respectively) is small but nonzero owing to the zero differential overlap approximation. In such cases, the evolution equation can be re-written as80,82, 36

35

34

3

37 38

𝑖 39 41

40

𝑖 43

42

𝜕(𝑃0 + 𝑊(𝑡)) = [𝐹(𝑃0 + 𝑊(𝑡)), (𝑃0 + 𝑊(𝑡))] 𝜕𝑡

𝜕𝑊(𝑡) = [𝐹(𝑃0 + 𝑊(𝑡)), (𝑃0 + 𝑊(𝑡))] − [𝐹0 , 𝑃0 ] 𝜕𝑡

(20) (21)

Using a small scaling factor 𝑔 (equal to 10-5 times the norm of W), the von Neumann operator, within linear-response, is constructed as, 46

45

4

48

47 49

𝐿𝑊(𝑡) =

[𝐹(𝑃0 + 𝑔𝑊(𝑡)), (𝑃0 + 𝑔𝑊(𝑡))] − [𝐹0 , 𝑃0 ] 𝜕𝑊(𝑡) = −𝑖 𝜕𝑡 𝑔

(22)

50

where L is called the Liouville superoperator. This way the response of the system can be studied under the influence of a weak delta function pulse.59 54

53

52

51

Chebyshev Propagator: A general overview of different time propagation schemes can be found in reviews by Kosloff124 and by Castro and coworkers.125 Here we propagate Eq. 22 using the iterative Chebyshev algorithm126 57

56

5

58 59 60 ACS Paragon Plus Environment

Page 7 of 23

Journal of Chemical Theory and Computation

7 1 2

𝐿 𝑊(𝑡) = 𝑒 𝐿∆𝑡 𝑊 𝑠𝑡𝑎𝑟𝑡 = ∑(2 − 𝛿𝑛0 )𝐽𝑛 (𝛼∆𝑡)𝑇𝑛 ( ) 𝑊 𝑠𝑡𝑎𝑟𝑡 = ∑(2 − 𝛿𝑛0 )𝐽𝑛 (𝛼∆𝑡)𝜁𝑛 𝛼 5

4

3

𝑛

6

(23)

𝑛

where 𝐽𝑛 are Bessel functions of the first kind of order n, ∆𝑡 is the time step, and α is a parameter required to scale eigenvalues of the Fock matrix into the range [−1,1]. In our implementation 10

9

8

7

3 𝛼 = (𝜀𝑚𝑎𝑥 − 𝜀𝑚𝑖𝑛 ) (24) 2 13

12

1

where 𝜀𝑚𝑎𝑥 and 𝜀𝑚𝑖𝑛 are the maximum and minimum eigenvalues of the Fock matrix, respectively. Usually, α needs to be slightly bigger than the spectral range of the Fock matrix for efficient propagation. 𝜁𝑛 is the modified Chebyshev series80, 17

16

15

14

𝐿 𝜁𝑛 = 𝑇𝑛 ( ) 𝑊 𝑠𝑡𝑎𝑟𝑡 𝛼 21

20

19

18

(25)

The modified Chebyshev polynomials can be evaluated recursively using following equations, 2

𝜁0 = 𝑊 𝑠𝑡𝑎𝑟𝑡 23 25

24

𝜁1 = 26

𝐿𝑊 𝑠𝑡𝑎𝑟𝑡

27 28

𝜁𝑛 = 2 29 31

30

𝛼

(26)

𝐿𝜁𝑛−1 + 𝜁𝑛−2 𝛼

From the above equations, it is clear that the Chebyshev propagator only involves matrix multiplications. However, in principle, determination of the parameter α requires on-the-fly diagonalization of the Fock matrix. For weak or no external perturbations, where fluctuations in the eigenspectrum are small, only a single diagonalization step is required to obtain the scaling factor. In our implementation the Chebyshev expansion is truncated to an accuracy of 10-10. 38

37

36

35

34

3

32

Dipole Moment: The time-dependent dipole moment is computed in the atomic orbital (AO) basis according to 41

40

39

𝜇(𝑡) = 𝑇𝑟[𝐷𝑃(𝑡)] 43

42

(27)

Real-time calculations are performed for each polarization of the external field pulse and the complex on-diagonal polarizability tensor is computed as 45

4

47

46

𝛼𝑖𝑖 = 49

48

𝜇𝑖 (𝜔) 𝑖 = 𝑥, 𝑦, 𝑧 𝐸𝑖 (𝜔)

(28)

where 𝜇𝑖 (𝜔) and 𝐸𝑖 (𝜔) are the Fourier transforms of the dipole moment and electric fields in the i direction. The absorption cross section formula is obtained from127,62, 65 52

51

50

53

𝜎(𝜔) = 56

5

54

4𝜋𝜔 𝐼𝑚[𝛼(𝜔)] (29) 𝑐

and the resulting dipole strength function is 57 58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 8 of 23

8 1 2 3

1 𝑆(𝜔) = 𝑇𝑟[𝜎(𝜔)] 3 4 6

5

(30)

Computational Details: 8

7

All calculations were performed with a development version of NWChem. Geometries for P3B2 and f-coronene were obtained from reference 108 and the geometry for ubiquitin was obtained from reference 104, respectively. RT-TDDFT and LR-TDDFT results with B3LYP/6-31G(d) method for P3B2 and f-coronene were obtained from reference 108. We have also performed INDO/S-CIS calculations within a (20e,20o) active space using the ZINDO-MN128 code to compare with RT-INDO/S results for P3B2 and f-coronene. INDO/S-CIS results are also reported in the supporting information for comparison. We note that the trends of our RT-INDO/S results compared to INDO/S-CIS are consistent with INDO/S-RPA results of Baker and Zerner and coworkers42,43 when compared with INDO/S-CIS, suggesting that RT and RPA calculations at the INDO/S level would show agreement similar to that we have previously observed with DFT (where RT-TDDFT is formally equivalent to RPA, but we do not have access to an INDO/S-RPA implementation). The gas-phase geometry of isolated tyrosine was optimized at the B3LYP129/ccpVTZ130 level. QM/MM calculations were performed to prepare the solvated betanin dye system in two solvents, water and methanol. Solvent partitions of 3080-water and 1263-methanol molecules were assigned as the MM regions, and the betanin molecule was treated as the QM region, respectively. The water molecules were assigned SPC/E force field parameters131, while the methanol molecules were assigned AMBER ff99 parameters132. For the betanin molecule, the Def2-SVP basis set133 and hybrid PBE0134 functional was used. Initially, both the water-solvated and methanol-solvated betanin system geometries were optimized. After a 1 ps of equilibration with the Berendsen thermostat (NVT)135 at 298 K and a time step of 0.25 fs, QM/MM molecular dynamics simulations were performed for 10 ps. 10 snapshots of the trajectory were extracted to generate cluster models having a 7 Å thick shell of solvent molecules. Cluster generation and visualization were done with VMD136. 38

37

36

35

34

3

32

31

30

29

28

27

26

25

24

23

2

21

20

19

18

17

16

15

14

13

12

1

10

9

For all systems considered, starting from the ground state, RT-INDO/S calculations were carried out under the influence of a weak delta function electric field pulse with a field strength of 0.0001 a.u. A time step of 0.1 a.u. (0.0024fs) was used and the system was propagated for 10000 time steps (24 fs). Since infinite time propagation is not practical, we multiply the time-dependent dipole moment by a damping factor 𝑒 −𝑡/𝜏 with τ = 200 a.u. before taking the Fourier transform in order to produce the absorption spectra. This has the effect of introducing artificial and identical lifetimes for all of the excited states. 46

45

4

43

42

41

40

39

47 48

Results and Discussion: 50

49

First, we have performed calculations on two medium size organic molecules, P3B2 and fcoronene, composed of 130 and 162 atoms, respectively. We compare our results with previous RT-TDDFT and LR-TDDFT calculations.108 As mentioned before, performing LR-TDDFT or INDO/S-CIS (see, e.g., Figures S1 and S2) calculations for the same energy range for which realtime calculations are performed would require the calculation of several thousands of roots. Realtime methods, on the other hand, are able to obtain spectra for a broad energy range because the 58

57

56

5

54

53

52

51

59 60 ACS Paragon Plus Environment

Page 9 of 23

Journal of Chemical Theory and Computation

9 1 3

2

delta function electric field pulse samples all frequencies and the resolution of a spectral window depends only on the time step and length of propagation. For instance, to resolve an excitation of energy ω, one requires a time step < π / ω, and typically a tenth of that for accuracy. In both cases, RT-INDO/S spectra show qualitative agreement with both RT-TDDFT and LR-TDDFT (Figures 2 and 3). However, the intensities of different peaks predicted by RT-INDO/S are different from RT-TDDFT, which is not unexpected for approximate Hamiltonians and the limitations of the Slater-type minimal basis sets. 1

10

9

8

7

6

5

4

12 13 14 15 16 17 18 19 20

Figure 1. Molecular structures of (a) P3B2, (b) f-coronene. Blue, light blue and white colors represent nitrogen, carbon and hydrogen, respectively. 23

2

21

24 25 26 27 28 29 30 31 32 3 34 35 36 37 39

38

Figure 2. RT-INDO/S, LR-TDDFT and RT-TDDFT spectra of P3B2. LR-TDDFT and RT-TDDFT calculations were performed using the B3LYP/6-31G(d)137 method. Only 200 roots were calculated with LR-TDDFT. LR-TDDFT and RT-TDDFT spectra were adapted with permission from reference 108. Copyright 2015 American Chemical Society. 4

43

42

41

40

45 46 47 48 49 50 51 52 53 54 5 56 57 58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 10 of 23

10 1 2 3 4 5 6 7 8 9 10 1 12 13 14 15 16 17 19

18

Figure 3. RT-INDO/S, LR-TDDFT and RT-TDDFT spectra of f-coronene. LR-TDDFT and RTTDDFT calculations were performed using B3LYP/6-31G(d) method. Only 40 roots were calculated with LR-TDDFT. LR-TDDFT and RT-TDDFT spectra were adapted with permission from reference 108. Copyright 2015 American Chemical Society. 25

24

23

2

21

20

While P3B2 and f-coronene are fairly large systems for RT-TDDFT calculations, performed with the B3LYP hybrid exchange-correlation functional, for INDO/S these systems are relatively small in size. To evaluate the performance of the RT-INDO/S method for considerably larger organic systems, we have also calculated the UV/Vis absorption spectrum of a small protein ubiquitin. Ubiquitin is a regulatory protein consisting of 76 amino acids. UV/Vis spectra for the entire protein in the gas phase and in aqueous solution are available and the low energy part of the spectrum includes a distinct feature associated with the single tyrosine chromophore in ubiquitin. QM/MM138 and TD-DFTB104,105 calculations have been performed for the entire protein and the tyrosine molecule in the gas phase. QM/MM calculation predicted a redshift of ~4 nm (~0.06 eV) upon embedding the chromophore in the protein environment. However, TD-DFTB calculations predict that low energy part of the spectra is not dominated only by tyrosine spectra. Other excitations observed in the low energy part of the spectrum appear due to the underestimation of charge transfer excitations,105 a well-known issue in DFTB method139,140. Ubiquitin has 1231 atoms, so that even in a small energy range, there can be thousands of single orbital transitions, making linear response and CI approaches prohibitively expensive to capture all relevant roots. However, RT-INDO/S calculations for tyrosine in the gas phase and in the ubiquitin environment are practical and predict a blueshift of ~0.04 eV in the spectra (Figure 4) in near quantitative agreement with experiment. That suggests that this method is potentially suitable for studying environmental effects in large protein molecules. In addition, it is important to note that unlike TD-DFTB calculation, low energy part of the RT-INDO/S spectrum is only dominated by tyrosine. This is expected, since, semiempirical HF methods do not severely underestimate charge transfer excitations like DFTB methods.88,141,142 52

51

50

49

48

47

46

45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

30

29

28

27

26

53 54 5 56 57 58 59 60 ACS Paragon Plus Environment

Page 11 of 23

Journal of Chemical Theory and Computation

11 1 2 3 4 5 6 7 8 9 10 1 12 13 14 15 16 17 18 19 20 21 2 24

23

Figure 4. RT-INDO/S spectra of tyrosine in gas phase and in ubiquitin. Experimentally observed red shift of 0.06 eV in tyrosine spectra at 4.49 eV when embedded in ubiquitin protein environment is qualitatively (~0.04 eV) reproduced with the RT-INDO/S method. Red and black vertical lines indicate the peaks in tyrosine spectra in gas phase and ubiquitin spectra respectively. 30

29

28

27

26

25

Another area where the RT-INDO/S method may be particularly useful is for the computation of the UV/Vis spectra of organic dyes in explicit solvents. Betanin is an organic dye, well known for its application in the food, cosmetic, and pharmaceutical industries.143,144 The photophysical properties of betanin in aqueous and alcoholic solutions have been studied.145 Gas phase UV/Vis spectra of betanin computed with RT-INDO/S shows a maximum at 2.24 eV which agrees fairly well with an LR-TDDFT (PBE0/6-31+G(d,p)146 ) excitation energy of 2.35 eV (Figure S3). In both water and methanol solutions, the experimental UV/Vis absorption spectra of betanin exhibit a maximum at 2.32 eV.145 It has been suggested recently that large numbers of solvent molecules can be required to get convergence of excitation energy in a truncated model.27,29 To model the UV/Vis spectra of betanin in aqueous and methanol solutions, we have treated the complex and the surrounding solvent explicitly. We have extracted large solvated clusters, 837-861 atoms for water and 816-882 atoms for methanol [Figures 5b, 5d], from initial QM/MM simulations described earlier, respectively. 46

45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

47

RT-INDO/S calculations were performed on clusters extracted from 10 snapshots for each of the two solvents and averaged to generate a composite spectrum. Results obtained from individual spectra are reported in the supporting information. Our calculations predict maxima at 2.31 eV and 2.34 eV in methanol and water (Figure 6) in near quantitative agreement with experiment.145 These results indicate that the RT-INDO/S method, with existing parameters, can be used to accurately calculate solvatochromic shifts in the spectra of organic dyes in explicit solvents. 5

54

53

52

51

50

49

48

56 57 58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 12 of 23

12 1 2 3 4 5 6 7 8 9 10 1 12 13 14 15 16 17 18 19 20 21 2

Figure 5. Betanin in explicit (a) water during QM/MM simulation, (b) water for the truncated model, (c) methanol during QM/MM simulation, (d) methanol for the truncated model. Blue, light blue, red and white colors represent nitrogen, carbon, oxygen and hydrogen respectively. 26

25

24

23

27 28 29 30 31 32 3 34 35 36 37 38 39 40 41 42 43 4 45 46 47 48 49 50 51 52 53 5

54

Figure 6. RT-INDO/S spectra of betanin in methanol and water solution. 10 different snapshots are considered for each solvent. The averaged intensity is plotted for the final spectra. 57

56 58 59 60

ACS Paragon Plus Environment

Page 13 of 23

Journal of Chemical Theory and Computation

13 1 3

2

While we have demonstrated the ability of the RT-INDO/S method to predict UV/Vis spectrum for medium and large organic molecular complexes and also to compute solvatochromatic shifts, we note that the INDO/S model was developed essentially to compute vertical excitation energies and ionization potentials and is not suitable for optimizing geometries or dynamics, whether on the ground or excited state. Thus, prior to application of the RT-INDO/S model, appropriate geometries/trajectories must be obtained using other standard quantum or molecular mechanical techniques (or taken from experiment). In addition, our current implementation of RT-INDO/S is only appropriate for single excitations. Doubly and other higher excited states are not currently accessible. In addition, since minimal basis sets are used, Rydberg excited states are not accessible with this model. This problem may be overcome by including diffuse basis functions and further parametrization.147 Even with the limitations mentioned above, we anticipate that RT-INDO/S will be a useful technique for the study of the electronic excitations in large-scale systems. In addition, the real-time approach described here combined with semiempirical Hamiltonians may prove useful for the study of charge and exciton dynamics in large systems. 68-78,84-86 The INDO/S Hamiltonian has been previously used to calculate parameters in different model Hamiltonians to perform exciton dynamics successfully in different chemical systems.97,148,149 INDO/S has also been found to be accurate for calculating parameters to model charge transfer energies and electronic coupling, related to charge transfer rates in chemical systems. 83,150,151 Despite these successes, the INDO/S Hamiltonian has not been used in dynamical studies in combination with real-time approaches. We are currently investigating this as well as the possibility of using other effective Hamiltonians for dynamical applications. 29

28

27

26

25

24

23

2

21

20

19

18

17

16

15

14

13

12

1

10

9

8

7

6

5

4

30

Conclusion: 31 32

We have reported an implementation of a real-time INDO/S method for the computation of UV/Vis spectra for medium-sized systems like P3B2 and f-coronene, and large systems like gasphase ubiquitin and the betanin chromophore in explicit water and methanol solvents. It is encouraging that, with no adjustment to the underlying INDO/S parameters, UV/Vis spectra predicted at the RT-INDO/S level are in qualitative, and in the largest systems quantitative, agreement with either more accurate RT-TDDFT or experimental spectra. The low computational cost of this method makes this method quite general and opens up the possibility to study electronic excitations in large-scale systems in a fully quantum mechanical framework, including evaluating UV/Vis solvatochromic shifts for chromophores in explicit solvents and assessing environmental effects on the absorption spectra of chromophores embedded in macromolecules. Given the robust parameterization of INDO/S for transition-metal containing compounds, there is reason to expect similar utility of our RT-INDO/S approach when applied to such cases as well, although this remains to be assessed in further detail. 50

49

48

47

46

45

4

43

42

41

40

39

38

37

36

35

34

3

Acknowledgments: 52

51

The authors thank Dr. David Bowman for helpful discussions. This work was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences and the Office of Advanced Scientific Computing Research through the Scientific Discovery through Advanced Computing (SciDAC) program under Award Numbers DE-SC0008666 (S.G., C.J.C., L.G.) and KC-030106062653 58

57

56

5

54

53

59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 14 of 23

14 1 3

2

(A.A., N.G.). The research was performed using The Minnesota Supercomputing Institute (MSI) at the University of Minnesota and EMSL, a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research and located at the Pacific Northwest National Laboratory (PNNL). PNNL is operated by Battelle Memorial Institute for the United States Department of Energy under DOE contract number DE-AC05-76RL1830. 10

9

8

7

6

5

4

Supporting Information 12

1

INDO/S parameters, INDOS-CIS spectrum of P3B2 and f-coronene, RT-INDO/S and LR-TDDFT spectrum of betanin, details of RT-INDO/S spectra for different snapshots of solvated betanin molecule and molecular geometries for all of the systems discussed herein are available in the supporting information. 18

17

16

15

14

13

References 19 21

20

Akimov, A. V.; Prezdo, O. V. Large-Scale Computations in Chemistry: A Bird’s Eye View of a Vibrant Field. Chem. Rev. 2015, 115, 5797-5890. 2 Pomogaev, V.; Pomogaeva, A.; Aoki, Y. Absorption Spectra of Estradiol and Tryptophan Constructed by the Statistical and Elongation Methods. J. Phys. Chem. A 2009, 113, 1429-1433. 3 Kistler, K. A.; Matsika, S. Solvatochromic Shifts of Uracil and Cytosine Using a Combined Multireference Configuration Interaction/Molecular Dynamics Approach and the Fragment Molecular Orbital Method. J. Phys. Chem. A 2009, 113, 12396-12403. 4 Bose, S.; Chakrabarty, S.; Ghosh, D. Effect of Solvation on Electron Detachment and Excitation Energies of a Green Fluorescent Protein Chromophore Variant. J. Phys. Chem. B 2016, 120, 4410−4420. 5 Isborn, C. M.; Gotz, A. W.; Clarke, M. A.; Walker, M. A. Walker, R. C.; Martinez, T. J. Electronic Absorption Spectra from MM and ab Initio QM/MM Molecular Dynamics: Environmental Effects on the Absorption Spectrum of Photoactive Yellow Protein. J. Chem. Theory Comput. 2012, 8, 5092-5106. 6 Åqvist, J.; Arieh, W. Simulation of enzyme reactions using valence bond force fields and other hybrid quantum/classical approaches. Chem. Rev. 1993, 93, 2523−2544. 7 Shaik, S.; Cohen, S.; Wang, Y.; Chen, H.; Kumar, D.; Thiel, W. P450 Enzymes: Their Structure, Reactivity, and Selectivity—Modeled by QM/MM Calculations. Chem. Rev. 2010, 110, 949− 1017. 8 Nair, N. N.; Bredow, T.; Jug, K. Molecular dynamics implementation in MSINDO: Study of silicon clusters. J. Comput. Chem. 2004, 25, 1255−1263. 9 Rincón, L.; Gonzalez, C. A. Extended Hückel Tight-Binding Calculations of Electronic Resonances in Linear Chains of Gold Atoms and Clusters. J. Phys. Chem. C 2010, 114, 20734−20740. 10 Cabrera-Trujillo, J. M.; Robles, J. Theoretical study of the structural and electronic properties of two-dimensionally polymerized fullerene clusters with 2, 3, 4, and 7 C60 molecules. Molecules. Phys. Rev. B 2001, 64, 165408. 11 Saha, A.; Raghavachari, K. Dimers of Dimers (DOD): A New Fragment-Based Method Applied to Large Water Clusters. J. Chem. Theory Comput. 2014, 10, 58−67. 1

5

54

53

52

51

50

49

48

47

46

45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

30

29

28

27

26

25

24

23

2

56 57 58 59 60 ACS Paragon Plus Environment

Page 15 of 23

Journal of Chemical Theory and Computation

15 1 2 4

3 12

Lee, T. S.; York, D. M.; Yang, W. Linear‐scaling semiempirical quantum calculations for macromolecules. J. Chem. Phys. 1996, 105, 2744−2750. 13 Mei, Y.; Ji, C.; Zhang, J. Z. H. A new quantum method for electrostatic solvation energy of protein. J. Chem. Phys. 2006, 125, 094906. 14 Amari, S.; Aizawa, M.; Zhang, J.; Fukuzawa, K.; Mochizuki, Y.; Iwasawa, Y.; Nakata, K.;Chuman, H.; Nakano, T. VISCANA:  Visualized Cluster Analysis of Protein−Ligand Interaction Based on the ab Initio Fragment Molecular Orbital Method for Virtual Ligand Screening. J. Chem. Inf. Comput. Sci. 2006, 46, 221-330. 14

13

12

1

10

9

8

7

6

5

15

Du, S.; Sakurai, M. Multivariate analysis of properties of amino acid residues in proteins from a viewpoint of functional site prediction. Chem. Phys. Lett. 2010, 488, 81-85. 16 Fedorov, D. G.; Alexeev, Y.; Kitaura, K. Geometry Optimization of the Active Site of a Large System with the Fragment Molecular Orbital Method. J. Phys. Chem. Lett. 2011, 2, 282-288. 20

19

18

17

16

15

17

23

Arora, P.; Slipchenko, L. V.; Webb, S. P.; DeFusco, A.; Gordon, M. S. Solvent-Induced Frequency Shifts: Configuration Interaction Singles Combined with the Effective Fragment Potential Method. J. Phys. Chem. A 2010, 114, 6742-6750. 24

18

2

21

Yoo, S.; Zahariev, F.; Sok, S.; Gordon, M. S. Solvent effects on optical properties of molecules: a combined time-dependent density functional theory/effective fragment potential approach. J. Chem. Phys. 2008, 129, 144112. 19 Netzloff, H. M.; Collins, M. A. Ab Initio Energies of Nonconducting Crystals by Systematic Fragmentation. J. Chem. Phys. 2007, 127, 134113. 20 Babu, K.; Ganesh, V.; Gadre, S. R.; Ghermani, N. E. Tailoring approach for exploring electron densities and electrostatic potentials of molecular crystals. Theor. Chem. Acc. 2004, 111, 255−263. 21 Korchowiec, J.; Gu, F. L.; Aoki, Y. Elongation Method for Linear Scaling SCF Calculations of Polymers. J. Comput. Methods Sci. Eng. 2006, 6, 189−200. 22 Korchowiec, J.; de Silva, P.; Makowski, M.; Gu, F. L.; Aoki, Y. Elongation Cutoff Technique at Kohn-Sham Level of Theory. Int. J. Quantum Chem. 2010, 110, 2130−2139. 23 Xie, P.; Orimoto, Y.; Aoki, Y. An Efficient Local Molecular Dynamics Polymerization Simulation Combined with an Ab Initio MO Method. Materials 2013, 6, 870−885. 24 Fedorov, D. G.; Avramov, P. V.; Jensen, J. H.; Kitaura, K. Analytic gradient for the adaptive frozen orbital bond detachment in the fragment molecular orbital method. Chem. Phys. Lett. 2009, 477, 169-175. 25 Otaki, H.; Ando, K. The role of intermolecular hydrogen bond on dielectric properties in hydrogen-bonded material 5-bromo-9-hydroxyphenalenone: theoretical investigation. Phys. Chem. Chem. Phys. 2011, 13, 10719-10728. 26 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. 27 Provorse, M. R.; Peev, T.; Xiong, C.; Isborn, C. M. Convergence of Excitation Energies in Mixed Quantum and Classical Solvent: Comparison of Continuum and Point Charge Models. J. Phys. Chem. B 2016, 120, 12148-12159. 28 Zuehlsdorff, T. J.; Haynes, P. D.; Hanke, F.; Payne, M. C.; Hine, N. D. M. Solvent Effects on Electronic Excitations of an Organic Chromophore. J. Chem. Theory Comput. 2016, 12, 1853– 1861. 57

56

5

54

53

52

51

50

49

48

47

46

45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

30

29

28

27

26

25

58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 16 of 23

16 1 2 4

3 29

50

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, 13, 2159-2171. 30 Brunk, E.; Rothlisberger, U. Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. Chem. Rev. 2015, 115, 6217−6263. 31 Hückel, E. Quantentheoretische Beiträge zum Benzolproblem. Z Phys. 1931, 70, 204-286. 32 Praiser, R.; Parr, R. A Semi‐Empirical Theory of the Electronic Spectra and Electronic Structure of Complex Unsaturated Molecules. I. J. Chem. Phys. 1953, 21, 466-471. 33 Praiser, R.; Parr, R. A Semi‐Empirical Theory of the Electronic Spectra and Electronic Structure of Complex Unsaturated Molecules. II. J. Chem. Phys. 1953, 21, 767-776. 34 Pople, J. A. Electron interaction in unsaturated hydrocarbons. Trans. Faraday Soc. 1953, 49, 1375-1385. 35 Hoffmann, R. An Extended Hueckel Theory. I. Hydrocarbons. J. Chem. Phys. 1963, 39, 13971412. 36 Pople, J. A.; Segal, G. A. Approximate Self-consistent Molecular Orbital Theory. II. Calculations with Complete Neglect of Differential Overlap. J. Chem. Phys. 1965, 43, S136-S149. 37 Pople, J. A.; Segal, G. A. Approximate Self-consistent Molecular Orbital Theory. III. CNDO results for AB2 and AB3 System. J. Chem. Phys. 1966, 44, 3289-3296. 38 Pople, J. A.; Beveridge, D. L.; Dobosh, P. A. Approximate Self‐Consistent Molecular‐Orbital Theory. V. Intermediate Neglect of Differential Overlap. J. Chem. Phys. 1967, 47, 2026-2033. 39 Pople, J. A.; Santry, D. P.; Segal, G. A. Approximate Self-consistent Molecular Orbital Theory. I. Invariant Procedures. J. Chem. Phys. 1965, 43, S129-S135. 40 Ridley J.; Zerner M. C. An intermediate neglect of differential overlap technique for spectroscopy: Pyrrole and the azines. Theor. Chim. Acta. 1973, 32, 111–134. 41 Ridley J.; Zerner M. C. Triplet states via intermediate neglect of differential overlap: Benzene, pyridine and the diazines. Theor. Chim. Acta. 1976, 42, 223-236. 42 Baker, J. D.; Zerner, M. C. Characterization of the random phase approximation with the intermediate neglect of differential overlap Hamiltonian for electronic spectroscopy. J. Phys. Chem. 1991, 95, 8614-8619. 43 Baker, J. D.; Zerner, M. C. Applications of the random phase approximation with the INDO/S Hamiltonian: UV-VIS spectra of free base porphin. Chem. Phys. Lett. 1990, 175, 192-196. 44 Parkinson, W. A.; Zerner, M. C. Calculation of molecular polarizabilities for large systems within the random phase approximation. Chem. Phys. Lett. 1987, 139, 563-570. 45 Bacon A. D.; Zerner M. C. An intermediate neglect of differential overlap theory for transition metal complexes: Fe, Co and Cu chlorides. Theor. Chim. Acta. 1979, 53, 21-54. 46 Zerner M. C.; Loew, G. H.; Kirchner, R. F.; Mueller-Westerhof, U. T. An intermediate neglect of differential overlap technique for spectroscopy of transition-metal complexes. Ferrocene. J. Am. Chem. Soc. 1980, 102, 589–599. 51

47

49

48

47

46

45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

30

29

28

27

26

25

24

23

2

21

20

19

18

17

16

15

14

13

12

1

10

9

8

7

6

5

Cory, M. G.; Köstlmeier, S.; Kotzian, M.; Rösch, N.; Zerner M. C. An intermediate neglect of differential overlap technique for actinide compounds. J. Chem. Phys. 1994, 100, 1353-1365. 48 Kotzian, M.; Rösch, N.; Zerner, M. C. Intermediate neglect of differential overlap spectroscopic studies on lanthanide complexes. Theor. Chim. Acta. 1992, 81, 201-222. 49 Voityuk, A.A. INDO/X: A New Semiempirical Method for Excited States of Organic and Biological Molecules. J. Chem. Theory Comput. 2014,10, 4950-4958. 58

57

56

5

54

53

52

59 60 ACS Paragon Plus Environment

Page 17 of 23

Journal of Chemical Theory and Computation

17 1 2 4

3 50

Porezag, D.; Frauenheim, T.; Köhler, T.; Seifert, G.; Kaschner, R. Construction of tight-bindinglike potentials on the basis of density-functional theory: Application to carbon. Phys. Rev. B 1995, 51, 12947−12957. 8

7

6

5

51

Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, T.; Suhai, S.; Seifert, G. Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties. Phys. Rev. B: Condens. Matter Mater. Phys. 1998, 58, 7260. 52 Niehaus, T. A.; Suhai, S.; Della Sala, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, T. Tightbinding approach to time-dependent density-functional response theory. Phys. Rev. B: Condens. Matter Mater. Phys. 2001, 63, 085108. 53 Krüger, T.; Elstner, M.; Schiffels, P.; Frauenheim, T. Validation of the density-functional based tight-binding approximation method for the calculation of reaction energies and other data. J. Chem. Phys. 2005, 122, 114110. 54 Aradi, B.; Hourahine, B.; Frauenheim, T. DFTB+, a Sparse Matrix-Based Implementation of the DFTB Method. J. Phys. Chem. A 2007, 111, 5678−5684. 55 Elstner, M. SCC-DFTB: What Is the Proper Degree of Self-Consistency? J. Phys. Chem. A 2007, 111, 5614−5621. 56 Seifert, G. Tight-Binding Density Functional Theory:  An Approximate Kohn−Sham DFT Scheme. J. Phys. Chem. A 2007, 111, 5609−5613. 57 Niehaus, T.; della Sala, F. Range separated functionals in the density functional based tightbinding method: Formalism. Phys. Status Solidi B 2012, 249, 237-244. 58 Harumiya, K.; Kawata, I.; Kono, H.; Fujimura, Y. Exact two-electron wave packet dynamics of H2H2 in an intense laser field: Formation of localized ionic states H+H−. J. Chem. Phys. 2000, 113, 8953–8960. 32

31

30

29

28

27

26

25

24

23

2

21

20

19

18

17

16

15

14

13

12

1

10

9

59

3

Krause, P.; Klamroth, T.; Saalfrank, P. Time-dependent configuration-interaction calculations of laser-pulse-driven many-electron dynamics: Controlled dipole switching in lithium cyanide. J. Chem. Phys. 2005, 123, 074105. 60 Greenman, L.; Ho, P. J.; Pabst, S.; Kamarchik, E.; Mazziotti, D. A.; Santra, R. Implementation of the time-dependent configuration-interaction singles method for atomic strong-field processes. Phys. Rev. A 2010, 82, 023406. 40

39

38

37

36

35

34

61

Yabana, K.; Bertsch, G. F. Time-dependent local-density approximation in real time. Phys. Rev. B 1996, 54, 4484–4487. 62 Castro, A.; Appel, H.; Oliveira, M.; Rozzi, C. A.; Andrade, X.; Lorenzen, F.; Marques, M. A. L.; Gross, E. K. U.; Rubio, A. Octopus: a tool for the application of time-dependent density functional theory. Phys. Status Solidi B 2006, 243, 2465–2488. 63 Baer, R.; Neuhauser, D. Real-time linear response for time-dependent density-functional theory. J. Chem. Phys. 2004, 121, 9803-9807. 64 Theilhaber, J. Ab initio simulations of sodium using time-dependent density-functional theory. Phys. Rev. B 1992, 46, 12990–13003. 65 Lopata, K.; Govind, N. Modeling Fast Electron Dynamics with Real-Time Time-Dependent Density Functional Theory: Application to Small Molecules and Chromophores. J. Chem. Theory Comput. 2011, 7, 1344–1355. 66 Kulander, K. C. Time-dependent Hartree-Fock theory of multiphoton ionization: Helium. Phys. Rev. A 1987, 36, 2726–2738. 57

56

5

54

53

52

51

50

49

48

47

46

45

4

43

42

41

58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 18 of 23

18 1 2 4

3 67

Li, X.; Smith, S. M.; Markevitch, A. N.; Romanov, D. A.; Levis, R. J.; Schlegel, H. B. A timedependent Hartree–Fock approach for studying the electronic optical response of molecules in intense fields. Phys. Chem. Chem. Phys. 2005, 7, 233–239. 8

7

6

5

68

Hoff, D. A.; Rego, L. G. C. Chapter 4: Modelling electron quantum dynamics in large molecular systems. In Chemical Modelling; Joswig, J.-O., Springborg, M., Eds.; 2013; Vol. 10, pp 102−126. 69 Rego, L. G. C.; Hames, B. C.; Mazon, K. T.; Joswig, J. Intramolecular Polarization Induces Electron–Hole Charge Separation in Light-Harvesting Molecular Triads. J. Phys. Chem. C 2014, 118, 126−134. 70 Hoff, D. A.; Silva, R.; Rego, L. G. C. Subpicosecond Dynamics of Metal-to-Ligand ChargeTransfer Excited States in Solvated [Ru(bpy)3]2+ Complexes. J. Phys. Chem. C 2011, 115, 15617−15626. 71 Rego, L. G. C., Batista, V. S. Quantum Dynamics Simulations of Interfacial Electron Transfer in Sensitized TiO2 Semiconductors. J. Am. Chem. Soc., 2003, 125, 7989-7997. 72 Rego, L. G. C.; da Silva, R.; Freire, J. A.; Snoeberger, R. C.; Batista, V. S. Visible Light Sensitization of TiO2 Surfaces with Alq3 Complexes. J. Phys. Chem. C 2010, 114, 1317−1325. 73 Bowman, D. N.; Jakubikova, E. Low-Spin versus High-Spin Ground State in Pseudo-Octahedral Iron Complexes. Inorg. Chem. 2012, 51, 6011−6019. 74 Mukherjee, S.; Bowman, D. N.; Jakubikova, E. Cyclometalated Fe(II) Complexes as Sensitizers in Dye-Sensitized Solar Cells. Inorg. Chem. 2015, 54, 560−569. 75 Bowman, D. N.; Blew, J. H.; Tsuchiya, T.; Jakubikova, E. Elucidating Band-Selective Sensitization in Iron(II) Polypyridine-TiO2 Assemblies. Inorg. Chem. 2013, 52, 8621−8628. 76 Bowman, D. N.; Mukherjee, S.; Barnes, L. J.; Jakubikova, E. Linker dependence of interfacial electron transfer rates in Fe (II)-polypyridine sensitized solar cells. J. Phys.: Condens. Matter 2015, 27, 134205. 77 Torres, A.; Oliboni, R. S.; Rego, L. G. C. Vibronic and Coherent Effects on Interfacial Electron Transfer Dynamics. J. Phys. Chem. Lett. 2015, 6, 4927−4935. 78 da Silva, R.; Bortolini, G.; Torres, A.; Rego. L. G. C. A Nonadiabatic Excited State Molecular Mechanics/Extended Hückel Ehrenfest Method. J. Phys. Chem. C 2016, 120, 27688−27698. 79 Jakowski, J.; Morokuma, K. Liouville–von Neumann molecular dynamics J. Chem. Phys. 2009, 130, 224106. 80 Bartell, L. A.; Wall, M. R.; Neuhauser, D. A time-dependent semiempirical approach to determining excited states. J. Chem. Phys. 2010, 132, 234106. 81 Arnsten, C.; Lopata, K.; Wall, M. R.; Bartell, L.; Neuhauser, D. Modeling molecular effects on plasmon transport: Silver nanoparticles with tartrazine. J. Chem. Phys. 2011, 134, 084101. 82 Lopata, K.; Reslan, R.; Kowalska, M.; Neuhauser, D.; Govind, N.; Kowalski, K. Excited-State Studies of Polyacenes: A Comparative Picture Using EOMCCSD, CR-EOMCCSD(T), RangeSeparated (LR/RT)-TDDFT, TD-PM3, and TD-ZINDO. J. Chem. Theory Comput. 2011, 7, 3686– 3693. 83 Reslan, R.; Lopata, K.; Arntsen, C.; Govind, N.; Neuhauser, D. Electron transfer beyond the static picture: A TDDFT/TD-ZINDO study of a pentacene dimer. J. Chem. Phys. 2012, 137, 22A502. 84 Pal, S.; Trivedi, D. J.; Akimov, A. V.; Aradi, B.; Frauenheim, T.; Prezhdo, O. V. Nonadiabatic Molecular Dynamics for Thousand Atom Systems: A Tight-Binding Approach toward PYXAID. J. Chem. Theory Comput. 2016, 12, 1436−1448. 57

56

5

54

53

52

51

50

49

48

47

46

45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

30

29

28

27

26

25

24

23

2

21

20

19

18

17

16

15

14

13

12

1

10

9

58 59 60 ACS Paragon Plus Environment

Page 19 of 23

Journal of Chemical Theory and Computation

19 1 2 4

3 85

Long, R.; Fang, W.; Prezdo, O. V. Strong Interaction at the Perovskite/TiO2 Interface Facilitates Ultrafast Photoinduced Charge Separation: A Nonadiabatic Molecular Dynamics Study. J. Phys. Chem. C 2017, 121, 3797–3806. 86 Akimov, A. V. Nonadiabatic Molecular Dynamics with Tight-Binding Fragment Molecular Orbitals. J. Chem. Theory Comput. 2016, 12, 5719−5736. 87 Tretiak, S.; Mukamel, S. Density Matrix Analysis and Simulation of Electronic Excitations in Conjugated and Aggregated Molecules. Chem. Rev. 2002, 102, 3171-3212. 88 Voityuk, A. A. Intermediate neglect of differential overlap for spectroscopy. WIREs Comput. Mol. Sci. 2013, 3, 515−527. 89 Quarti C., Fazzi D., Del Zoppo M. A computational investigation on singlet and triplet exciton couplings in acene molecular crystals. Phys. Chem. Chem. Phys. 2011, 13, 18615–18625. 90 Nguyen, K. A.; Day, P. N.; Pachter, R.; Tretiak, S.; Chernyak, V.; Mukamel, S. Analysis of Absorption Spectra of Zinc Porphyrin, Zinc meso-Tetraphenylporphyrin, and Halogenated Derivatives. J Phys. Chem. A 2002, 106, 10285–10293. 21

20

19

18

17

16

15

14

13

12

1

10

9

8

7

6

5

91

Rand, B. P.; Cheyns, D., Vasseur, K.; Giebink, N. C.; Mothy, S.; Yi, Y. P.; Coropceanu, V.; Beljonne, D.; Cornil, J.; Bredas, J. L.; Genoe, J. The Impact of Molecular Orientation on the Photovoltaic Properties of a Phthalocyanine/Fullerene Heterojunction. Adv. Funct. Mater. 2012, 22, 2987–2995. 92 Hutchison, G. R.; Ratner, M. A.; Marks, T. J. Accurate Prediction of Band Gaps in Neutral Heterocyclic Conjugated Polymers. J. Phys. Chem. A 2002, 106, 10596–10605. 93 de Melo, J. S.; Pina, J.; Rodrigues, L.; M, Becker R. S. A comprehensive study of the spectral and photophysical properties of arylthiophenes. J. Photochem. Photobiol. A Chem. 2008, 194, 67– 75. 94 Beenken, W. J. D. Photo-induced charge transfer in fullerene–oligothiophene dyads – A quantum-chemical study. Chem Phys 2009, 357, 144–150. 95 Pourtois, G.; Beljonne, D.; Cornil, J.; Ratner, M. A.; Bredas, J. L. Photoinduced ElectronTransfer Processes along Molecular Wires Based on Phenylenevinylene Oligomers:  A QuantumChemical Insight. J Am. Chem. Soc. 2002, 124, 4436–4447. 38

37

36

35

34

3

32

31

30

29

28

27

26

25

24

23

2

96

39

Voityuk, A. A.; Michel-Beyerle, M. E.; Rösch, N. Quantum chemical modeling of structure and absorption spectra of the chromophore in green fluorescent proteins. Chem. Phys. 1998, 231, 13– 25. 97 Cory, M. G.; Zerner, M. C.; Xu, X. C.; Shulten, K. Electronic Excitations in Aggregates of Bacteriochlorophylls. J. Phys. Chem. B 1998, 102, 7640–7650. 98 Broo, A.; Holmén, A. Calculations and Characterization of the Electronic Spectra of DNA Bases Based on ab Initio MP2 Geometries of Different Tautomeric Forms. J. Phys. Chem. A 1997, 101, 3589–3600. 49

48

47

46

45

4

43

42

41

40

99

Starikov, E. B.; Cuniberti, G.; Tanaka, S. Conformation Dependence of DNA Exciton Parentage. J. Phys. Chem. B 2009, 113, 10428–10435. 100 Gieseking, R. L.; Ratner, M. A.; Schatz, G.C. Semiempirical Modeling of Ag Nanoclusters: New Parameters for Optical Property Studies Enable Determination of Double Excitation Contributions to Plasmonic Excitation. J. Phys. Chem. A 2016, 120, 4542-4549. 101 Gieseking, R. L.; Ratner, M. A.; Schatz, G. C. Quantum Mechanical Identification of Quadrupolar Plasmonic Excited States in Silver Nanorods. J. Phys. Chem. A 2016, 120, 9324– 9329. 58

57

56

5

54

53

52

51

50

59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 20 of 23

20 1 2 4

3 102

Lobaugh, J., Rossky, P. J. Computer Simulation of the Excited State Dynamics of Betaine-30 in Acetonitrile. J. Phys. Chem. A. 1999, 103, 9432-9447. 103 Lobaugh, J., Rossky, P. J. Solvent and Intramolecular Effects on the Absorption Spectrum of Betaine-30. J. Phys. Chem. A. 2000, 104, 899-907. 104 Rüger, R.; Lenthe, E. V.; Lu, Y.; Frenzel, J.; Heine, T.; Visscher, L. Efficient Calculation of Electronic Absorption Spectra by Means of Intensity-Selected Time-Dependent Density Functional Tight Binding. J. Chem. Theory Comput. 2015, 11, 157−167. 105 Rüger, R.; Lenthe, E. V.; Lu, Y.; Frenzel, J.; Heine, T.; Visscher, L. Correction to Efficient Calculation of Electronic Absorption Spectra by Means of Intensity-Selected Time-Dependent Density Functional Tight Binding. J. Chem. Theory Comput. 2017, DOI: 10.1021/acs.jctc.7b00564. 106 Lopata, K.; Govind, N. Near and Above Ionization Electronic Excitations with Non-Hermitian Real-Time Time-Dependent Density Functional Theory. J. Chem. Theory Comput. 2013, 9, 49394946. 107 Wang, Y.; Lopata, K.; Chambers, S. A.; Govind, N.; Sushko, P.V. Optical Absorption and Band Gap Reduction in (Fe1–xCrx)2O3 Solid Solutions: A First-Principles Study. J. Phys. Chem. C 2013, 117, 25504-25512. 108 Tussupbayev, S.; Govind, N.; Lopata, K.; Cramer, C. J. Comparison of Real-Time and LinearResponse Time-Dependent Density Functional Theories for Molecular Chromophores Ranging from Sparse to High Densities of States. J. Chem. Theory Comput. 2015, 11, 1102−1109. 109 Fischer, S. A.; Cramer, C. J.; Govind, N. Excited State Absorption from Real-Time TimeDependent Density Functional Theory. J. Chem. Theory Comput. 2015, 11, 4294-4303. 110 Fischer, S. A.; Cramer, C. J.; Govind, N. Excited-State Absorption from Real-Time TimeDependent Density Functional Theory: Optical Limiting in Zinc Phthalocyanine. J. Chem. Phys. Lett. 2016, 7, 1387-1391. 111 Nguyen, T. S.; Koh, J. H.; Lefelhocz, S.; Parkhill, J. Black-Box, Real-Time Simulations of Transient Absorption Spectroscopy. J. Chem. Phys. Lett. 2016, 7, 1590-1595. 112 Norman, P.; Bishop, D. M.; Jensen, H. J. A.; Oddershede, J. Near-resonant absorption in the time-dependent self-consistent field and multiconfigurational self-consistent field approximations. J. Chem. Phys. 2001, 115, 10323−10334. 113 Norman, P.; Bishop, D. M.; Jensen, H. J. A.; Oddershede, J. J. Chem. Phys. 2005, 123, 194103. 114 Kauczor, J.; Norman, P. Nonlinear response theory with relaxation: The first-order hyperpolarizability. J. Chem. Theory Comput. 2014, 10, 2449−2455. 115 Jensen, L.; Autschbach, J.; Schatz, G. C. Finite lifetime effects on the polarizability within time-dependent density-functional theory. J. Chem. Phys. 2005, 122, 224115. 116 Devarajan, A.; Gaenko, A.; Autschbach, J. Two-component relativistic density functional method for computing nonsingular complex linear response of molecules based on the zeroth order regular approximation. J. Chem. Phys. 2009, 130, 194102. 117 Hubener, H.; Giustino, F. Linear optical response of finite systems using multishift linear system solvers.J. Chem. Phys. 2014, 141, 044117. 118 Ziegler, T.; Krykunov, M.; Cullen, J. The implementation of a self-consistent constricted variational density functional theory for the description of excited states. J. Chem. Phys. 2012, 136, 124107. 119 Grimme, S. Simplified Tamm−Dancoff density functional approach for the electronic excitation spectra of very large molecules. J. Chem. Phys. 2013, 138, 244104. 57

56

5

54

53

52

51

50

49

48

47

46

45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

30

29

28

27

26

25

24

23

2

21

20

19

18

17

16

15

14

13

12

1

10

9

8

7

6

5

58 59 60 ACS Paragon Plus Environment

Page 21 of 23

Journal of Chemical Theory and Computation

21 1 2 4

3 120

Bannwarth, C.; Grimme, S. A simplified time-dependent densityfunctional theory approach for electronic ultraviolet and circular dichroism spectra of very large molecules. Comput. Theor. Chem. 2014, 1040, 45−53. 121 Brabec, J.; Lin, L.; Shao, M.; Govind, N.; Yang, C.; Saad, Y.; Ng, E. G. Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT. J. Chem. Theory Comput. 2015, 11, 5197-5208. 122 Valiev, M.; Bylaska, E. J.; Govind, N.; Kowalski, K.; Straatsma,T. P.; Van Dam, H. J. J.; Wang, D.; Nieplocha, J.; Apra, E.; Windus, T. L.; de Jong, W. NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations. Comput. Phys. Commun. 2010, 181, 1477−1489. 123 Roothaan, C. C. J. New Developments in Molecular Orbital Theory. Rev. Mod. Phys. 1951, 23, 69-89. 124 Kosloff, R. Time-dependent quantum-mechanical methods for molecular dynamics. J. Phys. Chem. 1988, 92, 2087–2100. 125 Castro, A.; Marques, M. A. L.; Rubio, A. Propagators for the time-dependent Kohn-Sham equations. J. Chem. Phys. 2004, 121, 3425–3433. 126 Tal-Ezer, H.; Kosloff, R. An accurate and efficient scheme for propagating the time dependent Schrödinger equation. J. Chem. Phys. 1984, 81, 3967-3971. 127 Castro, A.; Marques, M. A. L.; Alonso, J. A.; Bertsch, G. F.; Yabana, K.; Rubio, A. Can optical spectroscopy directly elucidate the ground state of C20? J. Chem. Phys. 2002, 116, 1930-1933. 128 Zerner, M. C.; Ridley, J. E.; Bacon, A. D.; Edwards, W. D.; Head, J. D.; McKelvey, J.; Culberson, J. C.; Knappe, P.; Cory, M. G.; Weiner, B.; Baker, J. D.; Parkinson, W. A.; Kannis, D.; Yu, J.; Roesch, N.; Kotzian, M.; Tamm, T.; Karelson, M. M.; Zheng, X.; Pearl, G.; Broo, A.; Albert, K.; Cullen, J. M.; Li, J.; Hawkins, G. D.; Thompson, J. D.; Kelly, C. P.; Liotard, D. A.; Cramer, C. J.; Truhlar, D. G. ZINDO-MN1.2, Quantum Theory Project; University of Florida, Gainesville, and Department of Chemistry, University of Minnesota, Minneapolis, 2005. 129 Becke, A. D. Density‐functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. 130 Dunning Jr., T. H. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007-1023. 131 Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P. The missing term in effective pair potentials. J. Phys. Chem. 1987, 91, 6269-6271. 132 Case, D. A.; Babin, V.; Berryman, J. T.; Betz, R. M.; Cai, Q.; Cerutti, D. S.; Cheatham III, T. E.; Darden, T. A.; Duke, R. E.; Gohlke, H.; Goetz, A. W.; Gusarov, S.; Homeyer, N.; Janowski, P.; Kaus, J.; Kolossváry, I.; Kovalenko, A.; Lee, T. S.; LeGrand, S.; Luchko, T.; Luo, R.; Madej, B.; Merz, K. M.; Paesani, F.; Roe, D. R.; Roitberg, A.; Sagui, C.; Salomon-Ferrer, R.; Seabra, G.; Simmerling, C. L.; Smith, W.; Swails, J.; Walker, R. C.; Wang, J.; Wolf, R. M.; Wu, X.; Kollman, P. A. AMBER 14, University of California: San Francisco, 2014. 133 Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297-3305. 134 Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158–6170. 135 Berendsen, H. J. C.; Postma, J. P. M.; Gunsteren, W. F. v.; DiNola, A.; Haak, J. R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684-3690. 57

56

5

54

53

52

51

50

49

48

47

46

45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

30

29

28

27

26

25

24

23

2

21

20

19

18

17

16

15

14

13

12

1

10

9

8

7

6

5

58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 22 of 23

22 1 2 4

3 136

6

Humphrey, W.; Dalke, A.; Schulten, K. J. VMD: Visual molecular dynamics. Mol. Graphics 1996, 14, 33-38. 7

137

5

Hariharan, P. C.; Pople, J. A. The influence of polarization functions on molecular orbital hydrogenation energies. Theor. Chim. Acta. 1973, 28, 213-222. 138 Bellina, B.; Compagnon, I.; Joly, L.; Albrieux, F.; Allouche, A.; Bertorelle, F.; Lemoine, J.; Antoine, R.; Dugourd, P. UV spectroscopy of entire proteins in the gas phase. Int. J. Mass Spectrom. 2010, 297, 36−40 Special Issue: Ion Spectroscopy. 139 Kranz, J. J.; Elstner, M.; Aradi, B.; Frauenheim, T.; Lutsker, V.; Garcia, A. D.; Niehaus, T. A. Time-Dependent Extension of the Long-Range Corrected Density Functional Based TightBinding Method. J. Chem. Theory Comput. 2017, 13, 1737-1747. 140 Humeniuk, A.; Mitric, R. Long-range correction for tight-binding TD-DFT. J. Chem. Phys. 2015, 143, 134120. 141 Metcalfe, R. A.; Lever, A. B. P. Tetraammineruthenium(II) and -ruthenium(III) Complexes of o-Benzoquinone Diimine and Their Redox Series. Inorg. Chem. 1997, 36, 4762-4771. 142 Kogh-Jespersen, K.; Westbrook, J. D.; Potenza, J. A.; Schugar, H. J. Molecular and Electronic Structure of (L-Histidinato)pentammineruthenium(III) Chloride Monohydrate, ( NH3)5RuIII( his)Cl3.H2O. X-ray Structure, Single-Crystal Polarized Charge-Transfer Spectra, and ab Initio and Semiempirical Molecular Orbital Calculations. J. Am. Chem. Soc. 1987, 109, 7025-7031. 143 Esatbeyoglu, T.; Wagner, A. E.; Schini-Kerth, V. B.; Rimbach, G. Betanin—A food colorant with biological activity. Mol. Nutr. Food Res. 2015, 59, 36-47. 144 Stintzing, F. C.; Carle, R. Functional properties of anthocyanins and betalains in plants, food, and in human nutrition. Trends Food Sci. Technol. 2004, 15, 19-38. 145 Wendel, M.; Nizinski, S.; Dorota, T.; Starzak, K.; Szot, D.; Prukala, D.; Sikorski, M.; Wybraniec, S.; Burdzinski, G. Time-resolved spectroscopy of the singlet excited state of betanin in aqueous and alcoholic solutions. Phys. Chem. Chem. Phys., 2015, 17, 18152-18158. 146 Hehre, W. J.; Ditchfield, R.; Pople, J. A. Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules. J. Chem. Phys. 1972, 56, 2257-2261. 147 Head, J. Inclusion of Rydberg orbitals in semiempirical calculations on excited states? Int. J. Quant. Chem. 2003, 95, 580-592. 148 Hu, X.; Ritz, T.; Damjanovic, A.; Schulten, K. Pigment Organization and Transfer of Electronic Excitation in the Photosynthetic Unit of Purple Bacteria. J. Phys. Chem. B 1997, 101, 3854-3871. 149 Emanuele, E.; Zakrzewska, K.; Markovitsi, D.; Lavery, R.; Milie, P. Exciton States of Dynamic DNA Double Helices: Alternating dCdG Sequences. J. Phys. Chem. B 2005, 109,16109-16118. 150 Voityuk, A. A. Assessment of semiempirical methods for the computation of charge transfer in DNA p-stacks. Chem. Phys. Lett. 2006, 427, 177-180. 49

48

47

46

45

4

43

42

41

40

39

38

37

36

35

34

3

32

31

30

29

28

27

26

25

24

23

2

21

20

19

18

17

16

15

14

13

12

1

10

9

8

50 151

51

Voityuk, A. A. Electronic couplings and on-site energies for hole transfer in DNA: Systematic quantum mechanical/molecular dynamic study. J. Chem. Phys. 2008, 128, 115101. 53

52 54 5 56 57 58 59 60

ACS Paragon Plus Environment

Page 23 of 23

Journal of Chemical Theory and Computation

23 1 2 3 5

4

For Table of Contents only 6 7 8 9 10 1 12 13 14 15 16 17 18 19 20 21 2 23 24 25 26 27 28 29 30 31 32 3 34 35 36 37 38 39 40 41 42 43 4 45 46 47 48 49 50 51 52 53 54 5 56 57 58 59 60 ACS Paragon Plus Environment