Adaptive Configuration Interaction for Computing Challenging

Sep 11, 2017 - We introduce and analyze various approaches for computing excited electronic states using our recently developed adaptive configuration...
1 downloads 18 Views 1MB Size
Subscriber access provided by UNIVERSITY OF THE SUNSHINE COAST

Article

Adaptive configuration interaction for computing challenging electronic excited states with tunable accuracy Jeffrey B. Schriber, and Francesco A. Evangelista J. Chem. Theory Comput., Just Accepted Manuscript • DOI: 10.1021/acs.jctc.7b00725 • Publication Date (Web): 11 Sep 2017 Downloaded from http://pubs.acs.org on September 14, 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 39

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

Adaptive configuration interaction for computing challenging electronic excited states with tunable accuracy Jeffrey B. Schriber and Francesco A. Evangelista∗ Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University Atlanta, Georgia, 30322, USA E-mail: [email protected]

Abstract

2

We introduce and analyze various approaches to compute excited electronic states using our 4

recently developed adaptive configuration interaction (ACI) method [J.B. Schriber and F.A. Evangelista, J. Chem. Phys. 144, 161106 (2016)]. These ACI methods aim at describing

6

multiple electronic states with equal accuracy, including challenging cases like multi-electron, charge transfer, and near-degenerate states. We develop both state-averaged and state-specific

8

approaches to compute excited states whose absolute energy error can be tuned by a userspecified energy error threshold, σ. State averaged schemes are found to be more efficient, in

10

that they obtain all states simultaneously in one computation, but lose some degree of statewise tunability. State-specific algorithms allow for direct control of the error of each state,

12

though states must be computed sequentially. We compare each method using methylene, LiF, and all-trans polyene benchmark data.

14

1. INTRODUCTION Perhaps the biggest difficulty in modeling electronic excited states is in quantifying varying de-

16

grees of dynamic and static electron correlation equally well among states of different character.

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

This challenge is exacerbated, for example, in multielectron excited states, which frequently occur 18

in chemically relevant systems containing transition metals or highly conjugated hydrocarbons. 1–3 Conventional excited-state electronic structure methods, including time-dependent density func-

20

tional theory (TD-DFT) 4,5 and the equation-of-motion coupled cluster hierarchy, 6,7 often cannot efficiently compute multielectron excited states due to their considerable multireference character

22

in which double or higher excitations dominate. Another example of challenging excited states are avoided crossings and conical intersections, which play a fundamental role in non-radiative pro-

24

cesses. 8–10 In these cases, two or more electronic states with potentially different character become degenerate or near-degenerate, a situation that is difficult to describe with single-reference methods

26

and state-specific multireference approaches, especially when the states involved are the groundand first-excited state. In this work, we generalize our previously developed adaptive configura-

28

tion interaction (ACI), 11 which is very well-suited for treating large numbers of strongly correlated electrons, to the general description of excited states, including multielectron and near-degenerate

30

states. In order to describe electronic states in which static correlation plays a dominant role, a mul-

32

tireference method is typically required. The most common choice is to build a multideterminantal reference wave function by permuting a chosen number of electrons in a set of selected active orbitals.

34

While these complete active space (CAS) wave functions, with (CASSCF) or without (CASCI) orbital optimization, are good references for excited state computations, the number of variational

36

parameters in CASCI wave functions grows combinatorially with the number of electrons and orbitals and quickly becomes prohibitive. Therefore, efficient alternatives to CASCI are required to

38

describe the excitations of large, chemically-relevant molecular systems. Wave function factorization techniques, 12–18 stochastic approaches, 19–28 symmetry projection

40

schemes, 29–35 and various truncated configuration interaction (CI) methods 11,36–56 have risen as the most economical, near-exact representations of compressed full CI (FCI) or CASCI wave functions.

42

Moreover, these techniques have allowed accurate descriptions of unprecedentedly large systems with significant multireference character by finding an efficient, sub-combinatorial scaling of the number

44

of parameters with respect to the size of the active space. The density matrix renormalization group (DMRG) 12–18 and full CI quantum Monte Carlo (FCIQMC) 20,22 have been particularly successful in

46

describing the ground states of large systems by correlating to near exactness active spaces beyond 2 ACS Paragon Plus Environment

Page 2 of 39

Page 3 of 39

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 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

66

for example, the renormalized basis states are in general not optimal for multiple roots, and can be sequentially optimized for each root in a state-specific scheme. Recently introduced alternatives

68

include derivations of Tamm-Dancoff (TDA) or random phase approximations (RPA) from a DMRG matrix product state reference, 65 in addition to time-dependent formalisms. 70 For FCIQMC, state-

70

specific algorithms can be achieved by using a Gaussian projector shifted near the energy of the desired root, 66 or by projecting out all previously computed lower roots from the solution of the

72

current root. 68,71,72 However, near-degenerate electronic states can pose a significant challenge for state-specific methods, and typically an excited state method that computes states simultaneously

74

is required. Alternatively, methods that compute states simultaneously usually incorporate some form of

76

state-averaging, in which a common basis is used to describe multiple states. In DMRG, this averaging is realized by renormalizing basis states with respect to multiple roots, though this procedure

78

generally requires more basis states to achieve accuracy comparable with ground state DMRG calculations. 63 Coe and Patterson introduced a modified Monte Carlo CI (MCCI), which stochastically

80

samples Hilbert space by selecting determinants based on a criterion averaged over roots of interests, and then forms and diagonalizes a matrix within the selected space. 73 Blunt et al. recently proposed

82

a modified FCIQMC method in which all states are computed simultaneously, and orthogonalized with a Graham-Schmidt procedure at each iteration. 68 Finally, Rodríguez-Guzmán and co-workers

84

presented a method in which ground and excited states are expanded in terms of nonorthogonal variationally-optimized symmetry-projected configurations. 74

86

Another method well-suited to accurately describe very large excited states is the variational 2RDM (v-2RDM) method, 75–80 which forgoes wave functions entirely and directly optimizes the one

88

and two-particle reduced density matrices subjected to N -representability conditions. The v-2RDM method can treat the lowest energy state of a given symmetry and total spin with high accuracy, 80,81

90

and its ability to compute excited states of the same symmetry based on the equations of motion of the 1-RDM has also been recently reported. 82

92

Two of the most commonly used selected CI methods, CIPSI 37 and MRD-CI, 38,83 both have been modified to compute excited states using a single space of selected determinants. 39,60,84,85

94

Though the definition of the selection importance criterion is different for each method, both add determinants to the model space based on the maximum value of the selection criterion with respect 4 ACS Paragon Plus Environment

Page 4 of 39

Page 5 of 39

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

96

Journal of Chemical Theory and Computation

to all electronic states. This approach is equivalent to adding determinants to the model space that are important to at least one state being studied. While the success of both methods has prompted

98

numerous excited state studies, 86–92 [See Ref. 93 and references therein] there has been a lack of comparison of numerical results to determine the most effective way to compute excited states with

100

selected CI methods. In this work, we provide a detailed analysis of both simultaneous and sequential ACI methods for

102

computing excited state wave functions. Although our analysis is focused on ACI, our comparisons of various excited state methods is applicable to selected CI methods in general. One of the hallmark

104

features of ACI is its ability to compute ground state energies to a user-defined accuracy. 11 Similar in spirit to the recently proposed variance-matching variational Monte Carlo method of Robinson and

106

Neuscamman, 94 which computes all states to a similar quality as monitored by the wave function’s energy variance, we seek an adaptation of ACI such that any excited state energy can be efficiently

108

computed with a specifiable error. As shown in this work, ACI can be generalized to compute excitation energies with state-specific error control and near perfect error cancellation.

110

2. THEORY 2.1. Review of ACI.

112

We begin with a brief overview of ACI. 11 ACI is an iterative selected CI

method which generates and optimal space of Slater determinants such that the total error in the ground state energy approximately matches a user defined parameter, σ: |E0FCI − E0ACI (σ)| ≈ σ,

114

(1)

where E0FCI and E0ACI are the FCI and ACI ground state energies, respectively. The procedure is as follows:

116

1. At iteration k, we define a reference space, P (k) = {Φµ }, and associated reference wave function, (k)

ΨP =

X

c(k) µ Φµ ,

(2)

µ∈P (k)

118

(k)

where expansion coefficients cµ are obtained by diagonalizing the Hamiltonian in the P (k) space. The procedure is started with a reference space, P (0) , comprised of all determinants 5 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

120

Page 6 of 39

in a small CASCI (< 1000 determinants). 2. All singly and doubly excited determinants are generated from P (k) to form the first order

122

interacting space, F (k) . An estimate of the energy contribution to the reference space, ǫ(ΦI ), is computed for each determinant ΦI ∈ F (k) . This energy importance measure is defined as

124

(k)

the variational correlation energy of a two-state system involving ΨP and ΦI , and takes the form, 95 1 ǫ(ΦI ) = ∆I − 2

q  2 2 ∆I + 4VI ,

(3)

where the Epstein–Nesbet denominator (∆I ) and coupling (VI ) are defined as (k) ˆ (k) ˆ |ΦI i , ∆I = hΨP | H |ΨP i − hΦI | H (k) ˆ VI = hΨP | H |ΦI i .

126

(4) (5)

The generation of all singles and doubles has a computational cost that scales as |F (k) | ≈ |P (k) |NO2 NV2 , where |P (k) | is the number of determinants in the set P (k) , while NO and NV

128

are the number of occupied and virtual orbitals, respectively. The quantity VI is computed as VI =

ˆ c(k) µ hΦµ | H |ΦI i ,

X

(6)

µ∈P (k) 130

by an algorithm that generates all single and double excitations out of the determinants in P (k) . To reduce the memory necessary to store the vector VI , we apply screening and store

132

(k) ˆ |ΦI i | > τV , where τV is a only those elements for which one of the contributions |cµ hΦµ | H

user-specified screening threshold. Unless otherwise noted, calculations reported in this work 134

use a screening threshold τV = 10−12 Eh . 3. We define the secondary space, Q(k) , as the smallest possible subset of F (k) such that the

136

accumulation of error estimates for sampled determinants not in Q(k) is approximately equal to σ, the energy selection parameter. In practice, this truncation of F (k) is performed by

138

accumulating the ǫ(ΦI ) estimates in increasing order until σ is reached, X

|ǫ(ΦI )| ≤ σ.

ΦI ∈F (k) \Q(k)

6 ACS Paragon Plus Environment

(7)

Page 7 of 39

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

The total model space M (k) is then defined as the union of the secondary space and the 140

reference, M (k) = Q(k) ∪ P (k) .

(8)

In our current implementation, this step has a computational cost that scales as |F (k) | ln |F (k) | 142

due to sorting of the error estimates |ǫ(ΦI )|. In principle, sorting can be avoided by using an optimization algorithm that identifies an error threshold η such that the difference between

144

σ and the sum of the error estimates |ǫ(ΦI )| for the determinants with |ǫ(ΦI )| < η is zero. This can be achieved by introducing a function f (η),

f (η) =

X

(9)

|ǫ(ΦI )|θ(η − |ǫ(ΦI )|) − σ,

ΦI ∈F (k)

146

148

and seeking the largest η such that f (η) ≤ 0. The function f (η) is monotonous (albeit not P continuous) and has one zero if σ ≤ ΦI ∈F (k) |ǫ(ΦI )| (otherwise, all determinants should be included). Therefore, η may be found by using the bisection method with cost proportional to |F (k) |.

150

4. The Hamiltonian is formed and diagonalized in the model space, M (k) , yielding a total energy (k)

EM and wave function, (k)

X

|ΨM i =

(k)

(10)

CI |ΦI i .

ΦI ∈M (k) 152

(k)

The procedure ends when EM converges, which occurs when the determinantal makeup of M (k) is identical between iterations [M (k) = M (k+1) ]. This diagonalization is usually the

154

most expensive step in the ACI procedure and has a cost between |M (k) | and |M (k) |NO2 NV2 . 5. If the energy is not converged, the updated reference P (k+1) space is generated from the

156

most important subset of the M (k) space. Specifically, P (k+1) is formed by storing determi(k)

(k)

nants in M (k) with the largest |CI |2 values, until the accumulation of the |CI |2 values of 158

determinants kept reaches a second user-defined value, X

|Cµ(k) |2 ≤ 1 − γσ,

Φµ ∈P (k+1)

7 ACS Paragon Plus Environment

(11)

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 8 of 39

where γ is a user-specified constant with units (energy)−1 . Following our previous work, 11 160

we set γ = 1 Eh−1 . This step requires sorting of the determinant coefficients and scales as |M (k) | ln |M (k) |. However, as noted above, sorting may also be avoided in this step. Analo-

162

gous to the initiator approximation in FCIQMC, 22 this step greatly increases the efficiency of ACI by only allowing determinants important in the model space to generate candidate

164

determinants in F (k) . The aimed selection strategy gives the user nearly perfect control over the error for ground state

166

calculations. This control is achieved in step 3 by accurately identifying determinants that need to be explicitly treated in the diagonalization procedure and ignoring those determinants whose cumu-

168

lative energy contribution is smaller than the desired error. Another recently proposed selected CI method, heat-bath CI (HCI), 52,62 uses a very similar algorithm to ACI with the important excep-

170

tion that its generation of the model space is done without sampling any unselected determinants. Specifically, it selects determinants individually using a threshold ǫ1 , only adding a determinant ΦI

172

to the model space for which the inequality, ∀ΦJ ∈ P (k) ,

|HIJ cJ | > ǫ1

(12)

holds. Sorting of the integrals allows the selection of determinants to be done with perfect efficiency, 174

meaning that no determinants are sampled that are not added to the model space. Additionally, HCI continually grows the model space, so that at iteration k + 1 all the model space determinants

176

are retained, P (k+1) = M (k) ,

(13)

rather than truncating it as is done with ACI. The resulting method is extremely efficient, but does 178

not provide an inherent means of estimating the total energy error. The predictive control over the error in ACI relies on a fairly complete sampling of unselected determinants, and a potentially

180

fruitful development in ACI would be to somehow adopt the heat-bath sampling method while still preserving the predictive power of the current method. In the following section, we generalize the

182

ground state ACI method to determinant spaces capable of describing excited state wave functions.

8 ACS Paragon Plus Environment

Page 9 of 39

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

2.2. Excited State Methods in ACI. 184

In principle, excited state wave functions can be obtained

nearly automatically by solving for the lowest eigenvalues of the Hamiltonian in the ground state determinant space (M0 ) optimized by ACI. However, the determinant set M0 is biased towards

186

the ground state and provides no way to control the accuracy of higher eigenvalues. One solution is a state-averaged (SA) approach in which the definition of the selection importance criterion,

188

ǫ(ΦI ), is modified to reflect importance with respect to multiple roots, resulting in a single space of determinants capable of describing all roots of interest. Alternatively, an algorithm may be preferred

190

wherein determinant selection is done separately for each root, allowing each state to be described by potentially unique determinant spaces. We expand upon both approaches below.

192

For all methods considered, we generalize our definition of the energy importance criterion, ǫ(ΦI ), to apply to any arbitrary state n, q  1 2 , ∆I,n − ∆2I,n + 4VI,n 2

(14)

(k) (k) ˆ ˆ |ΦI i , |ΨP,n i − hΦI | H ∆I,n = hΨP,n | H

(15)

ǫ(ΦI , n) = where,

(k) ˆ VI,n = hΨP,n | H |ΦI i ,

194

(16)

(k)

and |ΨP,n i refers to the n-th eigenstate of the P (k) space. For converged wave functions, the (k) superscript will be dropped.

196

By combining different definitions of selection criteria and model spaces, it is possible to formulate several excited state methods. In this work we consider five excited state ACI methods, which,

198

for convenience are summarized in Table 1. We separate these approaches into two categories: methods that use a single model space for all states, and methods that use or combine different

200

model spaces. What follows is a detailed description of each excited state method.

2.2.1. State-averaged methods based on a single model space. 202

The goal of a state-

averaged ACI (SA-ACI) procedure is to obtain a single compact space of determinants capable of describing multiple electronic states with controlled accuracy. Algorithmically, the only differ-

9 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

204

Page 10 of 39

ence in SA-ACI compared to the ground state algorithm is the definition of the energy importance criterion, ǫ(ΦI , n), where we propose two ways to define a determinant’s importance in lowering the

206

energy of multiple states. One option is to simply take the average importance of a determinant (|ǫ(ΦI , n)|) among N roots of interest,

ǫ¯(ΦI ) =

N −1 1 X |ǫ(ΦI , n)|. N

(17)

n=0

208

In addition, the reference space is updated by selecting determinants based on an average of the expansion coefficient, −1  1 NX  |Cµ,n |2 ≤ 1 − γσ, N (k+1)

X

Φµ ∈P 210

(18)

n=0

where Cµ,n is the CI coefficient from root n corresponding to the reference space determinant Φµ . This averaged procedure, denoted SA-ACIavg , only requires one optimization to obtain multiple

212

roots, which are rigorously orthogonal as they are eigenvectors of a common Hamiltonian. However, the near-exact energy control in ground state ACI is exchanged for a control over the average error

214

among roots. The loss of state-specific error control may become potentially hazardous if one state has a

216

significantly different determinantal makeup with respect to the other computed roots, where the optimized space will become biased towards those roots described by similar determinants. This

218

biasing can be remedied either by adding user-specified weights to each state, or by defining the selection importance criterion as a maximum function rather than an average. Specifically, we can

220

redefine ǫ¯(ΦI ) as, ǫ¯(ΦI ) =

max

0≤n≤N −1

(19)

|ǫ(ΦI , n)|,

and similarly the maximum of the expansion coefficient among roots is then used to truncate the 222

model space M to update the reference. This procedure, denoted as SA-ACImax , will build a final model space of determinants which are important to at least one root of interest rather than all to

224

roots of interest in averaged way.

2.2.2. Methods based on state-specific model spaces. 226

Perhaps the most attractive feature

of ground state ACI is that the error is predicted by σ to remarkable accuracy. SA-ACI methods 10 ACS Paragon Plus Environment

Page 11 of 39

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

achieve this in an averaged way because all states are computed with the same set of determinants. 228

To better preserve error prediction in higher roots, we propose a second class of ACI methods which form separate model spaces Mn with corresponding wave functions ΨMn for each root n of

230

interest, a procedure we define as unconstrained ACI. Unfortunately, the model space wave functions, ΨMn , are in general not mutually orthogonal when optimized separately. In fact, the degree of

232

nonorthogonality between states can be viewed as a measure of the inaccuracy of individually optimized wave functions, since all states become orthogonal in the FCI limit (σ → 0). We propose

234

three methods which attempt to preserve state-wise energy control while maintaining orthogonality among roots.

236

Two of the following methods begin with an unconstrained ACI calculation, yielding separate determinant spaces for each root. In order to correctly converge each state to the right solution,

238

a reasonable initial guess is required. To achieve correct convergence to each root, we start each computation with three or four warmup iterations of SA-ACI, select the desired root, and then

240

optimize it individually. To avoid the possibility of root flipping, we do all optimizations within a root-following scheme

242

State Combined ACI. The N non-orthogonal states are first optimized separately with unconstrained ACI to produce a set of converged model spaces {Mn }, with which we can define a final

244

model space MSC as the union of all individually selected spaces,

MSC = M0 ∪ M1 ∪ · · · ∪ MN −1 .

(20)

Finally, the matrix representation of the Hamiltonian in the MSC space can be diagonalized to 246

yield orthogonal eigenstates for the N roots. This state-combined procedure (SC-ACI) can also be viewed as a more rigorous way to remove any determinant space biasing in SA-ACIavg . While the

248

energy errors derived from the individually converged model subspaces {Mn } are well controlled by σ, the final energies computed from the final model space MSC may not be. Since this final model

250

space will have more determinants than any individually optimized space (unless all spaces Mn are identical), then the energies of the final model space will be lower than their individually optimized

252

counterparts. This energy lowering effect should not be considered entirely a good consequence, as the lowering may be inconsistent among roots, resulting in poorer error cancellation. Nonetheless, 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

254

Page 12 of 39

σ remains an approximate upper bound to the expected error for each state. Multistate ACI. Another approach to overcome the nonorthogonality problem in uncon-

256

strained ACI is to simply recouple the nonorthogonal solutions to yield a multistate ACI (MS-ACI). Specifically, we form the multi-state Hamiltonian in the basis of individually optimized states, and

258

then solve the resulting generalized eigenvalue problem for N roots,

HC = SCE,

(21)

where the Hamiltonian (H) and overlap (S) matrices are defined as ˆ |ΨMn i (H)mn = hΨMm | H

260

m, n = 0, . . . , N − 1,

(22)

and (S)mn = hΨMm |ΨMn i

m, n = 0, . . . , N − 1,

(23)

while E = diag(E0 , . . . , EN −1 ) is a diagonal matrix. Though simple, this approach is vulnerable to 262

numerical instability when the eigenvalues of the overlap matrix, S, are nearly zero, corresponding to near linear-dependencies in the basis of individually optimized states. Near-degenerate electronic

264

states may also pose a significant challenge for the initial computation of nonorthogonal wave functions, however, and issues in correctly identifying the root to be optimized and root-flipping could

266

preclude the use of an ACI method. Additionally, nonorthogonal wave functions create difficulties for certain properties and computations that require multiple wave functions, such as transition

268

dipole moments. Orthogonality-constrained ACI. The final excited state method, orthogonality constrained

270

ACI (OC-ACI), sequentially optimizes determinant spaces Mn in the orthogonal complement of the space spanned by the previous roots. This strategy has been successfully applied in the context

272

of density functional theory to obtain both valence 96 and core-excited states. 97 In this procedure, the selection and screening for P (k) and M (k) spaces is identical to the unconstrained ACI, but the

274

computation of the energies and wave functions for these spaces differs. Specifically, in the optimization of the n-th root, we define the projector onto the space com-

12 ACS Paragon Plus Environment

Page 13 of 39

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

Table 1 Comparison of excited state ACI methods.

Method

Energy importance criterion

SA-ACIavg

1 X ǫ(ΦI , n) N +1

Model space

N

Single average space (Mavg )

n=0

SA-ACImax SC-ACI MS-ACI OC-ACI

276

Single average space (Mmax )

max ǫ(ΦI , n) 0≤n≤N        ǫ(ΦI , n)      

n Mn )

Union of individual model spaces (MSC =

S

Union of individual model spaces (MMS =

S

n Mn )

Different spaces per state (M0 ⊆ M1 ⊆ M2 . . .)

ˆ M,n ) plementary to the first n − 1 roots (Q

ˆ M,n = 1 − Q

n−1 X

|ΨM,m i hΨM,m | ,

(24)

m=0

where |ΨM,m i is the optimized wave function for root m computed from the Mm model space. Thus, 278

this projector makes the n-th root orthogonal to all n−1 roots in the intersection of their subspaces. At each iteration, the P (k) and M (k) space wave functions are computed by projecting out the converged wave functions of the previous roots, ˆQ ˆ M,n |Ψ (k) i = EP,n |Ψ (k) i ˆ M,n H Q P ,n P ,n

(25)

ˆ M,n H ˆQ ˆ M,n |Ψ (k) i = EM,n |Ψ (k) i , Q M ,n M ,n

(26)

where the lowest eigenstate of this projected Hamiltonian yields the desired excited state wave 280

function. Note that this procedure is equivalent to the Gram-Schmidt orthogonalization of an excited state wave function with respect to lower roots, with the requirement that the Gram-

282

Schmidt orthogonalized remains variationally optimal. In practice, Eqs. (25)-(26) are implemented by projecting out eigenvectors for the first n − 1 roots from the trial vectors in the Davidson–Liu

284

algorithm. 98,99 The excited state ACI methods are all developed in the freely available open-source software,

286

Forte, 100 written as a plugin to Psi4. 101,102

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

Page 14 of 39

3. RESULTS AND DISCUSSION

288

3.1. Methylene.

We compare all ACI excited state methods using benchmark data on the

three lowest singlet A1 states of methylene. These states provide a rigorous test for the various ACI 290

excited state methods due to their different determinantal character. 103–108 The 1 1 A1 state is closedshell with a single dominant configuration, while the 2 1 A1 state has considerable double excitation

292

character and the 3 1 A1 is a multi-configurational open-shell singlet biradical. The geometries used were taken from Ref. 105 and are of C2v symmetry. The cc-pVDZ basis was used, 109 augmented

294

with additional s functions for both the carbon and hydrogen atoms, as outlined in Ref. 105, and no orbitals were frozen, resulting in a FCI space equivalent to a CAS(8e,27o). Table 2 Methylene vertical excitation energies (eV) obtained with truncated CI, EOM-CCSD, CC3, and ACI methods. All computations use restricted Hartree–Fock orbitals and the modified cc-pVDZ basis described in Ref. 105 with no frozen orbitals, resulting in a CAS(8e,27o).

All values shown are errors with respect to FCI, for which the excitation energy is instead shown. Number of variational Excitation Energy Error (eV) 1 Method parameters (Np ) 2 A1 3 1 A1 CISD 3,607 2.639 1.990 CISDT 83,689 1.113 0.035 CISDTQ 1,020,759 0.046 0.032 EOM-CCSD CC3

3,607 83,689

1.456 0.471

−0.005 −0.004

0.056 0.021 0.028 −0.016 −0.007

0.220 0.095 0.063 0.026 0.008

54,575 112,563 70,253 23,563/24,401/47,049 23,563/24,094/47,875

0.003 0.002 0.003 < 0.001 < 0.001

0.027 0.015 0.013 0.008 < 0.001

77,145,700

4.655

6.514

σ = 0.25 eV SA-ACIavg SA-ACImax SC-ACI MS-ACI OC-ACI

6279 12684 8871 1513/2525/6250 1513/2059/6282 σ = 0.054 eV

SA-ACIavg SA-ACImax SC-ACI MS-ACI OC-ACI FCI

14 ACS Paragon Plus Environment

Page 15 of 39

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

296

Journal of Chemical Theory and Computation

In the limit σ → 0, all ACI methods are equivalent to FCI. However, the rates at which each method converges to FCI as σ is reduced will likely be different, as each method uses a

298

different selection algorithm. Table 2 shows the error in excitation energy for the two excited states studied computed with conventional CI truncated to double (SD), triple (SDT), and quadruple

300

(SDTQ) excitations, equation of motion coupled cluster with singles and doubles (EOM-CCSD), 110 approximate coupled cluster with triples (CC3), 111 and excited-state ACI methods. For CI and CC

302

methods, the computed 2 1 A1 excitation energy converges to the correct result much faster than for the 2 1 A1 vertical excitation energy. This result is well understood considering that the 3 1 A1

304

state, though a multireference open-shell singlet, is most significantly composed of determinants with excitation rank up to triples. For the 2 1 A1 state, however, higher excitations are required,

306

and the inclusion of triples in CC methods is required even to achieve sub-eV accuracy. For the ACI methods, this trend does not hold—the 2 1 A1 state is more easily described than is

308

the 3 1 A1 state, and only slightly. Regardless of excitation rank, the 2 1 A1 state is largely dominated by a single determinant, so the ACI methods converge very quickly to the FCI result for this typically

310

challenging case. With a threshold of σ = 0.054 eV, all ACI methods produce excitation energy errors on the order of 1 meV or better for the 2 1 A1 state. The 3 1 A1 excitation energy is more

312

challenging since it is multideterminantal, and for the OC-ACI and MS-ACI methods it requires additional determinants to achieve the same level of accuracy. Among the ACI methods, the OC-

314

ACI shows the best performance for methylene as it gives ∼1-8 meV accuracy while requiring much fewer determinants than the other methods. Even for the modest σ values shown, we see that ACI,

316

particularly OC-ACI, can compute excited state energies more accurately and efficiently than the CI and CC methods in our example.

318

To further understand the varying efficiencies of each ACI method, we can analyze the absolute error of each state with respect to the required number of determinants, shown in Figure 2. Note

320

that the efficiencies of SA-ACIavg , SA-ACImax , and SC-ACI are very similar for all three states, meaning that they produce comparable absolute errors with a similar number of determinants.

322

In contrast, OC-ACI and MS-ACI show an improved efficiency for the 1 1 A1 and 2 1 A1 states. These lowest two states can be treated with fewer determinants than the 3 1 A1 state to achieve

324

a similar absolute error, but this property cannot be exploited by SA-ACI and SC-ACI as they employ a single set of determinants for all states. In general, the efficiencies of ACI methods which 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 ACS Paragon Plus Environment

Page 16 of 39

Page 17 of 39

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

Table 3 Comparison of LiF avoided crossing curves. We report the average error for the lowest two Σ+ states, the non paralellism error (NPE)a , the standard deviation (σ∆E ) with respect to FCI, and the average number of variational parameters (Np ) required for each computation. For OC-ACI computations, the two numbers represent the average number of determinants required for ground and excited states, respectively. All values are reported in mEh , and computed from seven bond distances from 8.5 a.u. to 14.5 a.u. incremented by 1 a.u. in order to enable comparison with PMC-SD and SA-MCCI. We do not show data for MS-ACI and SC-ACI as they diverge near the avoided crossing. All ACI computations used restricted Hartree–Fock orbitals and the custom basis set defined in Ref. 112, with the two 1s-like molecular orbitals frozen on each atom resulting in a CAS(8,27).

Method CISDT CISDTQ

Average Error (mEh ) X 1 1 Σ+ 2 1 Σ+ 14.075 15.584 1.295 1.273

NPE (mEh ) X 2 1 Σ+ 9.825 10.231 1.275 1.378

σ∆E 4.400 0.593

Np 76,365 978,423

1 Σ+

SA-ACIavg (50) SA-ACIavg (10) SA-ACIavg (1) SA-ACIavg (0.5)

42.853 10.243 0.891 0.448

48.873 10.052 0.830 0.410

13.896 1.262 0.089 0.055

21.110 1.314 0.109 0.071

7.368 0.437 0.047 0.029

488 3,341 57,798 95,711

OC-ACI(50) OC-ACI(10) OC-ACI(1) OC-ACI(0.5)

46.185 10.721 0.963 0.471

46.150 9.685 0.871 0.442

2.845 1.043 0.102 0.038

11.205 1.593 0.077 0.041

2.234 0.721 0.055 0.020

246/261 1576/1929 37,986/38,675 65,144/67,980

– 1.571

– 1.143

– 1.000

– 2.000

0.586 0.497

2734 5 × 106

b SA-MCCI c PMC-SD a

NPE is computed as the minimum error along the curve subtracted from the maximum error along the curve. Values taken from Ref. 73, Np refers to configuration state functions. c Values computed from data in Ref. 71, Np refers to the approximate number of walkers in projector Monte Carlo simulation. b

near-degenerate electronic states. 112–114 Near equilibrium, the ground state wave function is qualita332

tively single-determinantal with ionic character, and the first excited state is an open-shell biradical singlet dominated by low-weighted pairs of covalent configurations, requiring a multireference treat-

334

ment. At the avoided crossing point, the ionic and covalent character swaps between these two states, thus requiring a balanced treatment of strong and weak correlation to achieve good accuracy

336

for both states. 115,116 All LiF computations used the Li (9s 5p)/[4s 2p]; F (9s 6p 1d)/[4s 3p 1d] basis defined in Ref.

338

112 and restricted Hartree–Fock orbitals with the 1s-like molecular orbitals frozen on each atom, resulting in the explicit treatment of 8 electrons in 27 orbitals. To ensure proper convergence, we use

340

a γ value of 0.1 mEh−1 . The FCI solution requires 7.7 × 107 determinants and predicts the avoided 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 ACS Paragon Plus Environment

Page 18 of 39

Page 19 of 39

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

360

Journal of Chemical Theory and Computation

SA-ACIavg methods. As shown in the previous example, the SA-ACImax method has an efficiency that is comparable or smaller than that of the SA-ACIavg scheme, therefore, we do not report results

362

for it. SA-ACImax results are not reported since this In Figure 3, we report absolute energy errors in

364

mEh for CISDT, CISDTQ, OC-ACI, and SA-ACIavg for σ = 10 and 1 mEh . CISDT and CISDTQ fail to describe the multireference character of the covalent state with the same accuracy as the ionic

366

state, while both states are described to similar accuracy in the ACI methods. For σ = 10 mEh , the OC-ACI(10) and SA-ACIavg (10) produce errors near 10 mEh throughout the curve, with deviations

368

up to about 2 mEh near the avoided crossing point. Decreasing σ to 1 mEh gives nearly identical results for the two ACI methods, which both consistently yield an error near the preselected value

370

without noticeable deviations at this scale. In Table 3 we report the number of variational parameters in each computation in addition to

372

the average error of each state, the non-parallelism error (NPE) and the standard deviation (σ∆E ) of each curve with respect to FCI. The OC-ACI(10) and SA-MCCI produce an avoided crossing

374

curve of similarly good quality using a modest number of parameters. To achieve a comparably accurate PES, PMC-SD requires each state to be computed to within 1 mEh , while all ACI methods

376

can produce a given σ∆E value by computing each state with an order of magnitude less accuracy. The systematic improvability allows the ACI methods to compute the avoided crossing with NPEs

378

on the order of 40 µEh . The OC-ACI and SA-ACIavg give reasonably similar energy errors for all chosen σ values, but

380

with OC-ACI requiring fewer determinants without exception. Furthermore, the OC-ACI shows consistently better NPEs for all σ values shown. This result is somewhat surprising in that the SA-

382

ACIavg is very well suited for this particular problem. While the determinantal character for each state is different and swaps at the avoided crossing point, the total determinantal makeup of both

384

states in combination remains fairly unchanged throughout dissociation. Despite the determinantal basis for the SA-ACIavg Hamiltonian remaining fairly constant, it still shows greater NPEs with

386

respect to OC-ACI because it lacks the direct state-specific error control of OC-ACI.

3.3. 388

Extended polyenes.

The trans-polyene series (Cn Hn+2 ) is a set of nearly one-dimensional

conjugated carbon chains, which contain numerous low-lying π-π ∗ excited states often with double 19 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

excitation character. As the chain length increases, the energies of π bonding and anti bonding or390

bital pairs converge, which further enhances the strong correlation effects present in these low-lying excited states. 117–119

392

In the following sections, we first compare all ACI methods with available FCI benchmarks on octatetraene to establish the capacity of each method to describe these low-lying transitions. Next,

394

we apply the OC-ACI to a series of extended polyenes.

3.3.1. Comparison of Methods. 396

To investigate the performance of each ACI method on

large numbers of dense states, we studied the eight lowest singlets of octatetraene, the all-trans polyene with eight carbon atoms and four double bonds. The geometry was obtained from a DFT

398

optimization at the B3LYP/cc-pVDZ level of theory using C2h symmetry. 120 The ACI and CASCI computations employed a CAS(8,16) π-double active space in which all eight π electrons were

400

correlated in sixteen π bonding and anti bonding orbitals, all within the cc-pVDZ basis. To achieve fast convergence in ACI computations, we use Pipek-Mezey split-localized orbitals 121,122 in which

402

the active π orbitals occupied in the initial RHF reference wave function are localized separately from the active π orbitals unoccupied in the reference. Additionally, the ACI computations were

404

performed without symmetry for better localization, though we report the results using C2h labels. The CAS(8,16) is within the scope of CASCI calculations, allowing direct comparison of absolute

406

and excitation energies computed with ACI. Table 4 shows the absolute and excitation energy errors with respect to CASCI of the lowest 8 states of octatetraene for all ACI methods computed with

408

σ = 10 and 1 mEh . Similarly with the LiF example, the unconstrained ACI computation required for MS-ACI and SC-ACI suffered from incurable variational collapse for σ = 10 mEh . In this case,

410

the 2 1 Bu state collapses into the 1 1 Bu state, so MS-ACI is not feasible, but we do report SC-ACI values determined from the seven roots the unconstrained ACI did converge. For σ = 1 mEh , all

412

correct roots were found in the unconstrained ACI. In Table 4, we see that OC-ACI and MS-ACI consistently give the lowest errors in excitation

414

energy, 0.24 and 0.20 mEh respectively, compared to SC-ACI (0.37 mEh ), SA-ACIavg (0.77 mEh ) , and SA-ACImax (0.56 mEh ) for σ = 1 mEh . However, SC-ACI consistently gives the lowest absolute

416

error for each state, and for many states the SA-ACI methods give lower absolute errors than do the state-specific ACI methods. Most notably, the ground state absolute energy errors for SC-ACI 20 ACS Paragon Plus Environment

Page 20 of 39

Page 21 of 39

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

Table 4 Comparison of excited state ACI methods for the lowest eight singlet electronic states of octatetraene (C8 H10 ) computed with σ = 1 and 10 mEh . All computations use split-localized orbitals, with an active space including all sixteen π orbitals defined in the cc-pVDZ basis. We report both the excitation energy errors (∆Eex ) and absolute energy errors (∆Eabs ) with respect to FCI in mEh .

X 1 Ag

3 1 Bu

5 1 Ag

Average of |∆E|

9.41 12.13 4.27 −1.85

∆Eex (σ = 10 mEh ) 6.75 9.40 12.06 6.12 12.69 12.52 3.45 4.49 4.32 −0.51 −0.62 −2.15

8.99 11.76 3.76 −2.20

7.43 9.93 4.10 −0.22

9.43 11.16 4.19 1.44

0.85 0.60 0.40 −0.22 −0.30

0.71 0.57 0.37 −0.20 −0.24

∆Eex (σ = 1 mEh ) 0.68 0.73 0.91 0.50 0.54 0.56 0.35 0.37 0.38 −0.18 −0.18 −0.22 −0.20 −0.23 −0.26

0.75 0.56 0.36 −0.23 −0.27

0.77 0.59 0.38 −0.14 −0.19

0.77 0.56 0.37 0.20 0.24

15.30 16.63 8.61 10.13

12.75 15.81 7.92 10.81

∆Eabs (σ = 10 mEh ) 10.09 12.74 9.79 16.36 7.10 8.14 12.15 12.04

15.39 16.20 7.97 10.51

12.33 15.45 7.41 10.46

10.76 13.61 7.76 12.44

11.59 13.44 7.32 11.40

1.12 0.81 0.59 0.92 0.88

0.98 0.78 0.56 0.95 0.94

∆Eabs (σ = 1 mEh ) 0.95 0.99 0.70 0.75 0.54 0.56 0.96 0.96 0.99 0.96

1.17 0.77 0.57 0.93 0.93

1.02 0.77 0.55 0.91 0.91

1.04 0.79 0.57 1.00 1.00

0.94 0.70 0.52 0.97 0.97

1 1 Bu

SA-ACIavg SA-ACImax SC-ACI OC-ACI

11.97 12.96 4.95 −2.53

SA-ACIavg SA-ACImax SC-ACI MS-ACI OC-ACI

SA-ACIavg 3.34 SA-ACImax 3.67 SC-ACI 3.65 OC-ACI 12.66

SA-ACIavg SA-ACImax SC-ACI MS-ACI OC-ACI

418

4 1 Ag

2 1 Ag

0.26 0.21 0.19 1.14 1.18

State 2 1 Bu 3 1 Ag

(0.19 mEh ), SA-ACIavg (0.26 mEh ), and SA-ACImax (0.21 mEh ) are significantly lower than the corresponding error of OC-ACI (1.18 mEh ) and MS-ACI (1.14 mEh ), which are much closer to

420

the value predicted by σ. It is this over-stabilization of the ground state in SC-ACI and SA-ACI methods which causes the large errors in the excitation energies. In fact, OC-ACI and MS-ACI

422

both give an average absolute energy error of 0.97 mEh , close to σ = 1.0 mEh , with little deviation from each state, resulting in very accurate excitation energies.

424

As mentioned with the previous methylene example, the state-specific ACI methods perform better for relative energies as they allow each state to be described by a different number of de-

426

terminants. Figure 4 shows the absolute error with respect to CASCI of each of the eight states 21 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 ACS Paragon Plus Environment

Page 22 of 39

Page 23 of 39

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

Table 5 Vertical excitation energies (eV) for the polyene series (C2n H2n+2 ) computed with various methods. All ACI, CASCI, and DRMG methods use a CAS(2n,2n), and the NOCI calculation was performed with 10 nonorthogonal determinants as described in Ref 129.

State 2 1 A− g 1 1 Bu− 3 1 A− g

CASCI 4.97 6.13 6.82

OC-ACI(10) 4.96 6.12 6.86

DMRG PAO-CASCI 64 6.33 7.49 7.95

DMRG CASSCF 64 4.69 5.88 6.60

3.97 4.94 5.77

3.97 4.95 5.82

5.40 6.30 7.01

3.76 4.74 5.59

4.10

6 (12,12)

2 1 A− g 1 1 Bu− 3 1 A− g

3.19 3.98 5.12

3.42 4.18 4.81

3.43 4.18 4.85

4.90 5.60 6.28

3.25 4.03 4.78

3.46

8 (16,16)

2 1 A− g 1 1 Bu− 3 1 A− g

2.50 3.10 3.99

3.12 3.71 4.31

4.60 5.15 5.71

2.93 3.57 4.20

2.94

10(20,20)

2 1 A− g 1 1 Bu− 3 1 A− g

2.04 2.51 3.11

12(24,24)

2 1 A− g 1 1 Bu− 3 1 A− g

3.04 3.45 3.89

4.42 4.85 5.31

2.73 3.25 3.78

n CAS 4 (8,8)

446

NOCI 129 5.09

CASCIMRMP 118 4.26 5.30 7.20

1.70 2.05 2.45

n = 12, we use a prescreening threshold value of τV = 10−8 in building the external space F (k) which introduces an error negligible with respect to the chosen value of σ = 10 mEh . For these polyenes,

448

we report the excited states with C2h symmetry labels in addition to the particle/hole symmetry label in order to distinguish states characterized by ionic (+) or covalent (−) valence structures,

450

as is done conventionally (see Ref. 64 for more details). Moreover, this labeling distinguishes the 1 1 Bu− state, which we compute, from the 1 1 Bu+ state, which we do not compute, despite being

452

the second excited state in experimental results. We do not compute the 1 1 Bu+ state because it is largely described by π − σ excitations outside of our selected active space.

454

1 − 1 − Table 1 shows vertical excitation energies for the 2 1 A− g , 1 Bu , and 3 Ag states as computed

with CASCI, OC-ACI(10), DMRG with a projected atomic orbital basis (DMRG-PAO-CASCI), 64 456

DMRG with orbital opitmization (DMRG-CASSCF), 64 non-orthogonal CI (NOCI), 129 and the multireferece Møller–Plesset method with a CASCI reference (CASCI-MRMP). 118 The DMRG meth-

458

ods used the same active space as the OC-ACI, and CASCI-MRMP used a subset of the π valence

23 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 ACS Paragon Plus Environment

Page 24 of 39

Page 25 of 39

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

4. CONCLUSIONS

478

In this study, we have detailed five methods for computing challenging excited states using the ACI formalism. These excited state methods iteratively optimize either a single space of determinants or

480

separate determinant spaces in describing multiple electronic states. We rigorously compared each method using challenging excited-state benchmarks characterized by strong correlation and double

482

excitation character, including methylene excited states, the LiF avoided crossing, and many excited states of large polyenes. In these comparisons, we find that the OC-ACI method is generally the most

484

effective and widely applicable because it shows the best error prediction and tunability while still maintaining orthogonality among solutions. We attribute its success to its flexibility; by allowing

486

separate determinant spaces to describe each electronic state, the method can use the appropriate determinants for each state such that their preselected errors are the same. We demonstrate the

488

ability of OC-ACI to compute large polyenes, with up to a CAS(24,24), with remarkable error cancellation, in addition to computing excited states with accuracy competitive with DMRG.

490

In its current form, OC-ACI is well-suited to compute excitation energies of small molecules to high accuracy. Despite its remarkable error cancellation, computations in larger basis sets including

492

extrapolations to the complete basis limit are required for experimental predictability. For small systems, OC-ACI could be a very valuable tool in the benchmarking of excited states, where single-

494

reference methods like CC3 currently represent the highest applied level of theory. In application to larger systems (>100 orbitals), OC-ACI and SA-ACI are most useful in generating reference wave

496

functions which then need to be connected to theories of dynamical correlation to enable comparison to experiment. We are actively working on these extensions.

498

While this article was in revision, we noticed a similar study posted to the arXiv 130 in which the aforementioned HCI is adapted for excited states. This work employs an algorithm similar to

500

SA-ACI in addition to a semi-stochastic multireference perturbation theory, all which allow the authors to accurately compute excitation energies using large basis sets.

502

ACKNOWLEDGEMENTS This work was supported by Emory University startup funds and by a Research Fellowship of the

504

Alfred P. Sloan Foundation. 25 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

ASSOCIATED CONTENT

506

Supporting information.

We provide Cartesian coordinates of methylene and all polyenes (text

file) in xyz format. For all reported excitation energies in the article, we also provide the absolute 508

energies for the corresponding states. Finally, for the extended polyene computations with SO-ACI, we report RHF energies, absolute SO-ACI energies, and the required number of determinants for

510

each state. This information is available free of charge via the Internet at http://pubs.acs.org.

References 512

(1) González, L.; Escudero, D.; Andrés, L. S. Progress and challenges in the calculation of electronic excited states. ChemPhysChem 2012, 13, 28–51.

514

(2) Daniel, C. Photochemistry and photophysics of transition metal complexes: quantum chemistry. Coord. Chem. Rev. 2015, 282-283, 19–32.

516

(3) Dreuw, A.; Head-Gordon, M. Single-reference ab initio methods for the calculation of excited states of large molecules. Chem. Rev. 2005, 105, 4009–4037.

518

(4) Runge, E.; Gross, E. K. U. Density-functional theory for time-dependent systems. Phys. Rev. Lett. 1984, 52, 997–1000.

520

(5) Burke, K.; Werschnik, J.; Gross, E. K. U. Time-dependent density functional theory: Past, present, and future. J. Chem. Phys. 2005, 123, 062206.

522

(6) Stanton, J. F.; Bartlett, R. J. The equation of motion coupled-cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited

524

state properties. J. Chem. Phys. 1993, 98, 7029–7039. (7) Christiansen, O.; Koch, H.; Jørgensen, P. The second-order approximate coupled cluster

526

singles and doubles model CC2. Chem. Phys. Lett. 1995, 243, 409–418. (8) Matsika, S.; Krause, P. Nonadiabatic events and conical intersections. Ann. Rev. Phys. Chem.

528

2011, 62, 621–643.

26 ACS Paragon Plus Environment

Page 26 of 39

Page 27 of 39

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

(9) Blancafort, L. Photochemistry and photophysics at extended seams of conical intersection. 530

ChemPhysChem 2014, 15, 3166–3181. (10) Zobel, J. P.; Nogueira, J. J.; González, L. The IPEA dilemma in CASPT2. Chemical Science

532

2017, 8, 1482–1499. (11) Schriber, J. B.; Evangelista, F. A. Communication: An adaptive configuration interaction

534

approach for strongly correlated electrons with tunable accuracy. J. Chem. Phys. 2016, 144, 161106.

536

(12) White, S. R.; Martin, R. L. Ab initio quantum chemistry using the density matrix renormalization group. J. Chem. Phys. 1999, 110, 4127–4130.

538

(13) Chan, G. K.-L.; Head-Gordon, M. Highly correlated calculations with a polynomial cost algorithm: A study of the density matrix renormalization group. J. Chem. Phys. 2002, 116,

540

4462–4476. (14) Chan, G. K.-L.; Dorando, J. J.; Ghosh, D.; Hachmann, J.; Neuscamman, E.; Wang, H.;

542

Yanai, T. In Frontiers in Quantum Systems in Chemistry and Physics; Wilson, S., Grout, P., Maruani, J., Delgado-Barrio, G., Piecuch, P., Eds.; Springer Netherlands, 2008; pp 49–65.

544

(15) Chan, G. K.-L.; Sharma, S. The density matrix renormalization group in quantum chemistry. Annu. Rev. Phys. Chem. 2011, 62, 465–481.

546

(16) Sharma, S.; Chan, G. K.-L. Spin-adapted density matrix renormalization group algorithms for quantum chemistry. J. Chem. Phys. 2012, 136, 124121.

548

(17) Olivares-Amaya, R.; Hu, W.; Nakatani, N.; Sharma, S.; Yang, J.; Chan, G. K.-L. The ab-initio density matrix renormalization group in practice. J. Chem. Phys. 2015, 142, 034102.

550

(18) Sharma, S. A general non-Abelian density matrix renormalization group algorithm with application to the C2 dimer. J. Chem. Phys. 2015, 142, 024107.

552

(19) Greer, J. C. Estimating full configuration interaction limits from a Monte Carlo selection of the expansion space. J. Chem. Phys. 1995, 103, 1821–1828.

27 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

554

(20) Booth, G. H.; Thom, A. J.; Alavi, A. Fermion Monte Carlo without fixed nodes: A game of life, death, and annihilation in Slater determinant space. J. Chem. Phys. 2009, 131, 054106.

556

(21) Purwanto, W.; Zhang, S.; Krakauer, H. Excited state calculations using phaseless auxiliaryfield quantum Monte Carlo: Potential energy curves of low-lying C2 singlet states. J. Chem.

558

Phys. 2009, 130, 094107. (22) Cleland, D.; Booth, G. H.; Alavi, A. Communications: Survival of the fittest: Accelerating

560

convergence in full configuration-interaction quantum Monte Carlo. J. Chem. Phys. 2010, 132, 041103.

562

(23) Coe, J. P.; Paterson, M. J. Development of Monte Carlo configuration interaction: natural orbitals and second-order perturbation theory. J. Chem. Phys. 2012, 137, 204108.

564

(24) Giner, E.; Scemama, A.; Caffarel, M. Using perturbatively selected configuration interaction in quantum Monte Carlo calculations. Can. J. Chem. 2013, 91, 879–885.

566

(25) Petruzielo, F. R.; Holmes, A. A.; Changlani, H. J.; Nightingale, M. P.; Umrigar, C. J. Semistochastic projector Monte Carlo method. Phys. Rev. Lett. 2012, 109, 230201.

568

(26) Ten-no, S. Stochastic determination of effective Hamiltonian for the full configuration interaction solution of quasi-degenerate electronic states. J. Chem. Phys. 2013, 138, 164126.

570

(27) Li Manni, G.; Smart, S. D.; Alavi, A. Combining the complete active space self-consistent field method and the full configuration interaction quantum monte carlo within a super-

572

CI framework, with application to challenging metal-porphyrins. J. Chem. Theor. Comput. 2016, 12, 1245–1258.

574

(28) Giner, E.; Assaraf, R.; Toulouse, J. Quantum Monte Carlo with reoptimised perturbatively selected configuration-interaction wave functions. Mol. Phys. 2016, 114, 910–920.

576

(29) Jiménez-Hoyos, C. A.; Henderson, T. M.; Tsuchimochi, T.; Scuseria, G. E. Projected Hartree–Fock theory. J. Chem. Phys. 2012, 136, 164109.

28 ACS Paragon Plus Environment

Page 28 of 39

Page 29 of 39

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

578

Journal of Chemical Theory and Computation

(30) Jiménez-Hoyos, C. A.; Rodríguez-Guzmán, R.; Scuseria, G. E. Multi-component symmetryprojected approach for molecular ground state correlations. J. Chem. Phys. 2013, 139,

580

204102. (31) Samanta, K.; Jiménez-Hoyos, C. A.; Scuseria, G. E. Exploring copper oxide cores using the

582

projected Hartree–Fock method. J. Chem. Theory Comput. 2012, 8, 4944–4949. (32) Rivero, P.; Jiménez-Hoyos, C. A.; Scuseria, G. E. Entanglement and polyradical character

584

of polycyclic aromatic hydrocarbons predicted by projected Hartree–Fock theory. J. Phys. Chem. B 2013, 117, 12750–12758.

586

(33) Degroote, M.; Henderson, T. M.; Zhao, J.; Dukelsky, J.; Scuseria, G. E. Polynomial similarity transformation theory: A smooth interpolation between coupled cluster doubles and projected

588

BCS applied to the reduced BCS Hamiltonian. Phys. Rev. B 2016, 93, 125124. (34) Qiu, Y.; Henderson, T. M.; Scuseria, G. E. Communication: Projected Hartree–Fock theory

590

as a polynomial similarity transformation theory of single excitations. J. Chem. Phys. 2016, 145, 111102.

592

(35) Wahlen-Strothman, J. M.; Henderson, T. M.; Hermes, M. R.; Degroote, M.; Qiu, Y.; Zhao, J.; Dukelsky, J.; Scuseria, G. E. Merging symmetry projection methods with coupled cluster

594

theory: Lessons from the Lipkin model Hamiltonian. J. Chem. Phys. 2017, 146, 054110. (36) Bender, C. F.; Davidson, E. R. Studies in configuration interaction: The first-row diatomic

596

hydrides. Phys. Rev. 1969, 183, 23. (37) Huron, B.; Malrieu, J. P.; Rancurel, P. Iterative perturbation calculations of ground and

598

excited state energies from multiconfigurational zeroth-order wavefunctions. J. Chem. Phys. 1973, 58, 5745–5759.

600

(38) Buenker, R. J.; Peyerimhoff, S. D. Individualized configuration selection in CI calculations with subsequent energy extrapolation. Theor. Chim. Acta 1974, 35, 33–58.

602

(39) Evangelisti, S.; Daudey, J.-P.; Malrieu, J. P. Convergence of an improved CIPSI algorithm. Chem. Phys. 1983, 75, 91–102. 29 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

604

(40) Harrison, R. J. Approximating full configuration interaction with selected configuration interaction and perturbation theory. J. Chem. Phys. 1991, 94, 5021–5031.

606

(41) García, V. M.; Castell, O.; Caballol, R.; Malrieu, J. P. An iterative difference-dedicated configuration interaction. Proposal and test studies. Chem. Phys. Lett. 1995, 238, 222–229.

608

(42) Neese, F. A spectroscopy oriented configuration interaction procedure. J. Chem. Phys. 2003, 119, 9428–9443.

610

(43) Nakatsuji, H.; Ehara, M. Iterative CI general singles and doubles (ICIGSD) method for calculating the exact wave functions of the ground and excited states of molecules. J. Chem.

612

Phys. 2005, 122, 194108. (44) Abrams, M. L.; Sherrill, C. D. Important configurations in configuration interaction and

614

coupled-cluster wave functions. Chem. Phys. Lett. 2005, 412, 121–124. (45) Hu, W.; Chan, G. Excited-state geometry optimization with the density matrix renormaliza-

616

tion group, as applied to polyenes. J. Chem. Theor. Comput. 2015, 11, 3000–3009. (46) Bytautas, L.; Ruedenberg, K. A priori identification of configurational deadwood. Chem.

618

Phys. 2009, 356, 64–75. (47) Roth, R. Importance truncation for large-scale configuration interaction approaches. Phys.

620

Rev. C 2009, 79, 064324. (48) Evangelista, F. A. Adaptive multiconfigurational wave functions. J. Chem. Phys. 2014, 140,

622

054109. (49) Knowles, P. J. Compressive sampling in configuration interaction wavefunctions. Mol. Phys.

624

2015, 113, 1655–1660. (50) Liu, W.; Hoffmann, M. R. iCI: Iterative CI toward full CI. J. Chem. Theor. Comput. 2016,

626

12, 1169–1178. (51) Tubman, N. M.; Lee, J.; Takeshita, T. Y.; Head-Gordon, M.; Whaley, K. B. A deterministic

628

alternative to the full configuration interaction quantum Monte Carlo method. J. Chem. Phys. 2016, 145, 044112. 30 ACS Paragon Plus Environment

Page 30 of 39

Page 31 of 39

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

630

Journal of Chemical Theory and Computation

(52) Holmes, A.; Tubman, N.; Umrigar, C. J. Heat-bath configuration interaction: An efficient selected CI algorithm inspired by heat-bath sampling. J. Chem. Theor. Comput. 2016, 12,

632

3674–3680. (53) Zhang, T.; Evangelista, F. A. A deterministic projector configuration interaction approach

634

for the ground state of quantum many-body systems. J. Chem. Theor. Comput. 2016, 12, 4326–4337.

636

(54) Boschen, J. S.; Theis, D.; Ruedenberg, K.; Windus, T. L. Correlation energy extrapolation by many-body expansion. J. Phys. Chem. A 2017, 121, 836–844.

638

(55) Chien, A. D.; Zimmerman, P. M. Recovering dynamic correlation in spin flip configuration interaction through a difference dedicated approach. J. Chem. Phys. 2017, 146, 014103.

640

(56) Coe, J. P.; Paterson, M. J. Open-shell systems investigated with Monte Carlo configuration interaction. Int. J. Quant. Chem. 2016, 116, 1772–1782.

642

(57) Cimiraglia, R. Second order perturbation correction to CI energies by use of diagrammatic techniques: An improvement to the CIPSI algorithm. J. Chem. Phys. 1985, 83, 1746–1749.

644

(58) Caballol, R.; Malrieu, J. P. Direct selected configuration interaction using a hole-particle formalism. Chem. Phys. Lett. 1992, 188, 543–549.

646

(59) Hanrath, M.; Engels, B. New algorithms for an individually selecting MR-CI program. Chem. Phys. 1997, 225, 197–202.

648

(60) Angeli, C.; Persico, M. Multireference perturbation CI II. Selection of the zero-order space. Theor. Chem. Acc. 1997, 98, 117–128.

650

(61) McClean, J. R.; Aspuru-Guzik, A. Compact wavefunctions from compressed imaginary time evolution. RSC Advances 2015, 5, 102277.

652

(62) Sharma, S.; Holmes, A. A.; Jeanmairet, G.; Alavi, A.; Umrigar, C. J. Semistochastic heatbath configuration interaction method: Selected configuration interaction with semistochastic

654

perturbation theory. J. Chem. Theory Comput. 2017, 13, 1595–1604.

31 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

(63) Dorando, J. J.; Hachmann, J.; Chan, G. K.-L. Targeted excited state algorithms. J. Chem. 656

Phys. 2007, 127, 084109. (64) Ghosh, D.; Hachmann, J.; Yanai, T.; Chan, G. K.-L. Orbital optimization in the density

658

matrix renormalization group, with applications to polyenes and β-carotene. J. Chem. Phys. 2008, 128, 144117.

660

(65) Nakatani, N.; Wouters, S.; Van Neck, D.; Chan, G. K.-L. Linear response theory for the density matrix renormalization group: Efficient algorithms for strongly correlated excited

662

states. J. Chem. Phys. 2014, 140, 024108. (66) Booth, G. H.; Chan, G. K.-L. Communication: Excited states, dynamic correlation functions

664

and spectral properties from full configuration interaction quantum Monte Carlo. J. Chem. Phys. 2012, 137, 191102.

666

(67) Humeniuk, A.; Mitrić, R. Excited states from quantum Monte Carlo in the basis of Slater determinants. J. Chem. Phys. 2014, 141, 194104.

668

(68) Blunt, N. S.; Smart, S. D.; Booth, G. H.; Alavi, A. An excited-state approach within full configuration interaction quantum Monte Carlo. J. Chem. Phys. 2015, 143, 134117.

670

(69) Blunt, N. S.; Alavi, A.; Booth, G. H. Krylov-projected quantum Monte Carlo method. Phys. Rev. Lett. 2015, 115, 050603.

672

(70) Ronca, E.; Li, Z.; Jiménez-Hoyos, C. A.; Chan, G. K.-L. Time-step targeting time-dependent and dynamical density matrix renormalization group algorithms with ab initio Hamiltonians.

674

2017, arXiv:1706.09537. (71) Ohtsuka, Y.; Nagase, S. Projector Monte Carlo method based on Slater determinants: Test

676

application to singlet excited states of H2 O and LiF. Chem. Phys. Lett. 2010, 485, 367–370. (72) Blunt, N. S.; Booth, G. H.; Alavi, A. Density matrices in full configuration interaction quan-

678

tum Monte Carlo: Excited states, transition dipole moments, and parallel distribution. J. Chem. Phys. 2017, 146, 244105.

32 ACS Paragon Plus Environment

Page 32 of 39

Page 33 of 39

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

680

Journal of Chemical Theory and Computation

(73) Coe, J. P.; Paterson, M. J. State-averaged Monte Carlo configuration interaction applied to electronically excited states. J. Chem. Phys. 2013, 139, 154103.

682

(74) Rodríguez-Guzmán, R.; Jiménez-Hoyos, C. A.; Schutski, R.; Scuseria, G. E. Multireference symmetry-projected variational approaches for ground and excited states of the one-

684

dimensional Hubbard model. Phys. Rev. B 2013, 87, 235129. (75) Mazziotti, D. A. Pursuit of n-representability for the contracted Schrödinger equation through

686

density-matrix reconstruction. Phys. Rev. A 1999, 60, 3618–3626. (76) Mazziotti, D. A. First-order semidefinite programming for the direct determination of two-

688

electron reduced density matrices with application to many-electron atoms and molecules. J. Chem. Phys. 2004, 121, 10957.

690

(77) Mazziotti, D. A. Variational reduced-density-matrix method using three-particle nrepresentability conditions with application to many-electron molecules. Phys. Rev. A 2006,

692

74, 032501. (78) Mazziotti, D. A. Parametrization of the two-electron reduced density matrix for its direct

694

calculation without the many-electron wave function. Phys. Rev. Lett. 2008, 101, 253002. (79) DePrince III, A. E.; Mazziotti, D. A. Open-shell molecular electronic states from the para-

696

metric two-electron reduced-density-matrix method. J. Chem. Phys. 2009, 130, 164109. (80) Fosso-Tande, J.; Nascimento, D. R.; DePrince III, A. E. Accuracy of two-particle n-

698

representability conditions for describing different spin states and the singlet–triplet gap in the linear acene series. Mol. Phys. 2016, 114, 423–430.

700

(81) Fosso-Tande, J.; Nguyen, T.-S.; Gidofalvi, G.; DePrince III, A. E. Large-Scale Variational Two-Electron Reduced-Density-Matrix-Driven Complete Active Space Self-Consistent Field

702

Methods. J. Chem. Theor. Comput. 2016, 12, 2260–2271. (82) Jeffcoat, D. B.; DePrince III, A. E. N-representability-driven reconstruction of the two-

704

electron reduced-density matrix for a real-time time-dependent electronic structure method. J. Chem. Phys. 2014, 141, 214104. 33 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

706

(83) Buenker, R. J.; Phillips, R. A.; Krebs, S.; Liebermann, H.-P.; Alekseyev, A. B.; Funke, P. The Wuppertal multireference configuration interaction (MRD-CI) program system. Theor.

708

Chem. Acc. 2014, 133, 1–12. (84) Angeli, C.; Cimiraglia, R.; Persico, M.; Toniolo, A. Multireference perturbation CI I. Ex-

710

trapolation procedures with CAS or selected zero-order spaces. Theor. Chem. Acc. 1997, 98, 57–63.

712

(85) Krebs, S.; Buenker, R. J. Multireference configuration interaction study of the mixed Valence–Rydberg character of the C2 H4 1(π,π*)âĂĽV state. J. Chem. Phys. 1998, 106, 7208–

714

7214. (86) Banichevich, A.; Peyerimhoff, S. D. Theoretical study of the ground and excited states of

716

ozone in its symmetric nuclear arrangement. Chem. Phys. 1993, 174, 93–109. (87) Hogreve, H. The ground state geometry of the Ne+3 trimer. Chem. Phys. Lett. 1993, 215,

718

72–80. (88) Mawhinney, R. C.; Bruna, P. J.; Grein, F. Multireference configuration interaction studies on

720

metastable states of the dication BN2+ . J. Chem. Phys. 1995, 103, 8944–8954. (89) Hogreve, H. Ab initio study of the dication carbon trimer C2+ 3 . J. Chem. Phys. 1998, 102,

722

3281–3291. (90) Mennucci, B.; Toniolo, A.; Tomasi, J. Ab initio study of the electronic excited states in

724

4-(n,n-dimethylamino) benzonitrile with inclusion of solvent effects:âĂĽ the internal charge transfer process. J. Am. Chem. Soc. 2000, 122, 10621–10630.

726

(91) Żuchowski, P. S.; Guerout, R.; Dulieu, O. Ground- and excited-state properties of the polar and paramagnetic RbSr molecule: A comparative study. Phys. Rev. A 2014, 90, 012507.

728

(92) Chaieb, M.; Habli, H.; Mejrissi, L.; Oujia, B.; Gadéa, F. X. Ab initio spectroscopic study for the NaRb molecule in ground and excited states. Int. J. Quant. Chem. 2014, 114, 731–747.

34 ACS Paragon Plus Environment

Page 34 of 39

Page 35 of 39

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

730

Journal of Chemical Theory and Computation

(93) Szalay, P. G.; Müller, T.; Gidofalvi, G.; Lischka, H.; Shepard, R. Multiconfiguration selfconsistent field and multireference configuration interaction methods and applications. Chem.

732

Rev. 2012, 112, 108–181. (94) Robinson, P. J.; Neuscamman, E. Excitation variance matching with limited configuration

734

interaction expansions in variational monte carlo. 2017, arXiv:1705.04856. (95) Assfeld, X.; Almlöf, J. E.; Truhlar, D. G. Degeneracy-corrected perturbation theory for elec-

736

tronic structure calculations. Chem. Phys. Lett. 1995, 241, 438–444. (96) Evangelista, F. A.; Shushkov, P.; Tully, J. C. Orthogonality Constrained Density Functional

738

Theory for Electronic Excited States. J. Phys. Chem. A 2013, 117, 7378–7392. (97) Derricotte, W. D.; Evangelista, F. A. Simulation of X-ray absorption spectra with orthogonal-

740

ity constrained density functional theory. Phys. Chem. Chem. Phys. 2015, 17, 14360–14374. (98) Davidson, E. R. The iterative calculation of a few of the lowest eigenvalues and corresponding

742

eigenvectors of large real-symmetric matrices. J. Comp. Phys. 1975, 17, 87–94. (99) Liu, B. The simultaneous expansion for the solution of several of the lowest eigenvalues and

744

corresponding eigenvectors of large real-symmetric matrices; Numerical Algorithms in Chemistry: Algebraic Method, 1978.

746

(100) Forte, a suite of quantum chemistry methods for strongly correlated electrons. For the current version, see https://github.com/evangelistalab/forte, 2017.

748

(101) Turney, J. M.; Simmonett, A. C.; Parrish, R. M.; Hohenstein, E. G.; Evangelista, F. A.; Fermann, J. T.; Mintz, B. J.; Burns, L. A.; Wilke, J. J.; Abrams, M. L.; others, Psi4: an

750

open-source ab initio electronic structure program. WIREs: Comp. Mol. Sci. 2012, 2, 556– 565.

752

(102) Parrish, R. M.; Burns, L. A.; Smith, D. G. A.; Simmonett, A. C.; DePrince III, A. E.; Hohenstein, E. G.; Bozkaya, U.; Sokolov, A. Y.; Di Remigio, R.; Richard, R. M.; Gonthier, J. F.;

754

James, A. M.; McAlexander, H. R.; Kumar, A.; Saitow, M.; Wang, X.; Pritchard, B. P.; Verma, P.; Schaefer, H. F.; Patkowski, K.; King, R. A.; Valeev, E. F.; Evangelista, F. A.; 35 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

756

Turney, J. M.; Crawford, T. D.; Sherrill, C. D. Psi4 1.1: An Open-Source Electronic Structure Program Emphasizing Automation, Advanced Libraries, and Interoperability. J. Chem.

758

Theory Comput. 2017, Article ASAP. (103) Schaefer, H. F. Methylene: A paradigm for computational quantum chemistry. Science 1986,

760

231, 1100–1107. (104) Sherrill, C. D.; Leininger, M. L.; Van Huis, T. J.; Schaefer, H. F. Structures and vibrational

762

frequencies in the full configuration interaction limit: Predictions for four electronic states of methylene using a triple-zeta plus double polarization (TZ2P) basis. J. Chem. Phys. 1998,

764

108, 1040–1049. (105) Koch, H.; Christiansen, O.; Jørgensen, P.; Olsen, J. Excitation energies of BH, CH 2 and

766

Ne in full configuration interaction and the hierarchy CCS, CC2, CCSD and CC3 of coupled cluster models. Chem. Phys. Lett. 1995, 244, 75–82.

768

(106) Zimmerman, P. M.; Toulouse, J.; Zhang, Z.; Musgrave, C. B.; Umrigar, C. J. Excited states of methylene from quantum Monte Carlo. J. Chem. Phys. 2009, 131, 124103.

770

(107) Reynolds, P. J.; Dupuis, M.; Lester Jr, W. A. Quantum Monte Carlo calculation of the singlet–triplet splitting in methylene. J. Chem. Phys. 1985, 82, 1983–1990.

772

(108) Sherrill, C. D.; Van Huis, T. J.; Yamaguchi, Y.; Schaefer, H. F. Full configuration interaction benchmarks for the X 3B1, a 1 A1 , b 1 B1 and c 1 A1 states of methylene. J. Mol. Struct.

774

(Theochem) 1997, 400, 139–156. (109) Dunning Jr, T. H. Gaussian basis sets for use in correlated molecular calculations. I. The

776

atoms boron through neon and hydrogen. J. Chem. Phys. 1998, 90, 1007–1023. (110) Stanton, J. F.; Bartlett, R. J. The equation of motion coupled-cluster method. A systematic

778

biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties. J. Chem. Phys. 1998, 98, 7029–7039.

780

(111) Christiansen, O.; Koch, H.; Jorgensen, P. Response functions in the CC3 iterative triple excitation model. J. Chem. Phys. 1995, 103, 7429–7441. 36 ACS Paragon Plus Environment

Page 36 of 39

Page 37 of 39

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

782

Journal of Chemical Theory and Computation

(112) Bauschlicher Jr, C. W.; Langhoff, S. R. Full configuration-interaction study of the ionic–neutral curve crossing in LiF. J. Chem. Phys. 1988, 89, 4246–4254.

784

(113) Hanrath, M. Multi-reference coupled-cluster study of the ionic-neutral curve crossing LiF. Mol. Phys. 2008, 106, 1949–1957.

786

(114) Varandas, A. J. C. Accurate ab initio potential energy curves for the classic Li-F ioniccovalent interaction by extrapolation to the complete basis set limit and modeling of the

788

radial nonadiabatic coupling. J. Chem. Phys. 2009, 131, 124128. (115) Malrieu, J. P.; Heully, J.-L.; Zaitsevskii, A. Multiconfigurational second-order perturbative

790

methods: Overview and comparison of basic properties. Theoret. Chim. Acta 1995, 90, 167– 187.

792

(116) Finley, J.; Malmqvist, P.-Å.; Roos, B. O.; Serrano-Andrés, L. The multi-state CASPT2 method. Chem. Phys. Lett. 1998, 288, 299–306.

794

(117) Barford, W.; Bursill, R. J.; Lavrentiev, M. Y. Density-matrix renormalization-group calculations of excited states of linear polyenes. Phys. Rev. B 2001, 63, 195108.

796

(118) Kurashige, Y.; Nakano, H.; Nakao, Y.; Hirao, K. The π-π* excited states of long linear polyenes studied by the CASCI-MRMP method. Chem. Phys. Lett. 2004, 400, 425–429.

798

(119) Starcke, J. H.; Wormit, M.; Schirmer, J.; Dreuw, A. How much double excitation character do the lowest excited states of linear polyenes have? Chem. Phys. 2006, 329, 39–49.

800

(120) Becke, A. D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100.

802

(121) Pipek, J.; Mezey, P. G. A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions. J. Chem. Phys. 1989, 90,

804

4916–4926. (122) Bytautas, L.; Ivanic, J.; Ruedenberg, K. Split-localized orbitals can yield stronger configura-

806

tion interaction convergence than natural orbitals. J. Chem. Phys. 2003, 119, 8217–8224.

37 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

(123) Mikhailov, I. A.; Tafur, S.; Masunov, A. E. Double excitations and state-to-state transition 808

dipoles in π- π* excited singlet states of linear polyenes: Time-dependent density-functional theory versus multiconfigurational methods. Phys. Rev. A 2008, 77, 012510.

810

(124) Tavan, P.; Schulten, K. Electronic excitations in finite and infinite polyenes. Phys. Rev. B 1987, 36, 4337.

812

(125) Schmidt, M.; Tavan, P. Electronic excitations in long polyenes revisited. J. Chem. Phys. 2012, 136, 124309.

814

(126) Hachmann, J.; Cardoen, W.; Chan, G. K.-L. Multireference correlation in long molecules with the quadratic scaling density matrix renormalization group. J. Chem. Phys. 2006, 125,

816

144101. (127) Boguslawski, K. Targeting excited states in all-trans polyenes with electron-pair states. J.

818

Chem. Phys. 2016, 145, 234105. (128) Sokolov, A. Y.; Guo, S.; Ronca, E.; Chan, G. K.-L. Time-dependent N-electron valence

820

perturbation theory with matrix product state reference wavefunctions for large active spaces and basis sets: Applications to the chromium dimer and all-transpolyenes. J. Chem. Phys.

822

2017, 146, 244102. (129) Sundstrom, E. J.; Head-Gordon, M. Non-orthogonal configuration interaction for the calcu-

824

lation of multielectron excited states. J. Chem. Phys. 2014, 140, 114103. (130) Holmes, A. A.; Umrigar, C. J.; Sharma, S. Excited states using semistochastic heat-bath

826

configuration interaction. arXiv preprint arXiv:1708.03456v1 2017,

38 ACS Paragon Plus Environment

Page 38 of 39

Page 39 of 39

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 ACS Paragon Plus Environment