Benchmarking Post-Hartree–Fock Methods To Describe the Nonlinear

Oct 7, 2016 - Hales, Joel M.; Matichak, Jonathan; Barlow, Stephen; Ohira, Shino; Yesudas, Kada; Bredas, Jean-Luc; Perry, Joseph W.; ..... K.J. Thorley...
0 downloads 0 Views 1MB Size
Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Article

Benchmarking Post-Hartree Fock Methods to Describe the Nonlinear Optical Properties of Polymethines: An investigation of the accuracy of Algebraic Diagrammatic Construction (ADC) approaches Stefan Knippenberg, Rebecca L. Gieseking, Dirk R. Rehn, Sukrit Mukhopadhyay, Andreas Dreuw, and Jean-Luc Bredas J. Chem. Theory Comput., Just Accepted Manuscript • DOI: 10.1021/acs.jctc.6b00615 • Publication Date (Web): 07 Oct 2016 Downloaded from http://pubs.acs.org on October 8, 2016

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 38

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

Benchmarking Post-Hartree Fock Methods to Describe the Nonlinear Optical Properties of Polymethines: An investigation of the accuracy of Algebraic Diagrammatic Construction (ADC) approaches

Stefan Knippenberg,a,* Rebecca L. Gieseking,b,c Dirk R. Rehn,d Sukrit Mukhopadhyay,b,$ Andreas Dreuw,e and Jean-Luc Brédasb,f

a

Division of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, Roslagstullsbacken 15, S-106 91 Stockholm, Sweden b

School of Chemistry and Biochemistry and

Center for Organic Materials for All-Optical Switching (COMAS), Georgia Institute of Technology, Atlanta, Georgia 30332-0400 c

Department of Chemistry, Northwestern University, 2145 Sheridan Rd, Evanston, Illinois 60208

d

Department of Physics, Chemistry and Biology, Linköping University, S-581 83 Linköping, Sweden e

Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany

f

Physical Science and Engineering Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Kingdom of Saudi-Arabia

* [email protected] $

Present address: The Dow Chemical Company, Midland, Michigan 48674.

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

Abstract

Third-order nonlinear optical (NLO) properties of polymethine dyes have been widely studied for applications such as all-optical switching. However, the limited accuracy of the current computational methodologies has prevented a comprehensive understanding of the nature of the lowest excited states and their influence on the molecular optical and NLO properties. Here, attention is paid to the lowest excited-state energies and their energetic ratio, as these characteristics impact the figure-of-merit for all-optical switching. For a series of model polymethines, we compare several Algebraic Diagrammatic Construction (ADC) schemes for the polarization propagator with approximate second-order Coupled Cluster (CC2) theory, the widely used INDO/MRDCI approach and the SAC-CI algorithm incorporating singles and doubles linked excitation operators (SAC-CI SD-R). We focus in particular on the ground-toexcited state transition dipole moments and the corresponding state dipole moments, since these quantities are found to be of utmost importance for an effective description of the third-order polarizability  and two-photon absorption spectra. A sum-over-states expression has been used, which is found to quickly converge. While ADC(3/2) has been found to be the most appropriate method to calculate these properties, CC2 performs poorly.

2

ACS Paragon Plus Environment

Page 2 of 38

Page 3 of 38

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

1. Introduction

Polymethines with an odd number of methine groups have been widely applied in nonlinear optics (NLO) over the past several decades.1-5 Unlike most π-conjugated systems, these positively or negatively charged closed-shell compounds inherently have a molecular thirdorder polarizability γ with a very large negative real part and a small imaginary part (Re(γ) and Im(γ), respectively).1-3,5 This makes them particularly suitable for applications such as alloptical switching (AOS), which requires molecules with a figure-of-merit (FOM), defined as |Re(γ)/Im(γ)| > 4π.6 However, the optical and NLO properties of polymethines have proven particularly challenging to compute accurately.3,7,8

An important characteristic of polymethines is a large energy gap between the first excited state e, which has a large one-photon absorption (OPA) cross-section, and the second excited state e’, which has a significant two-photon absorption (TPA) cross section.1,2,9 If the energy ħω of the incoming light signal is tuned such that Ege < 2ħω < Ege’, the magnitude of Re(γ) is significantly increased due to pre-resonant enhancement because of the small energy difference between ħω and Ege. Also, since 2ħω is energetically located between states e and e’, Im(γ) is expected to be small.5 The required FOM for all-optical switching has recently been achieved in thiopyrylium and selenopyrylium polymethines by leveraging these unique molecular properties.5,10 In contrast, only a very small FOM is generally possible for polyenes and most other π-conjugated systems because the first two excited states e and e’ are much more similar in energy,11-13 and as a result, 2ħω must be < Ege to avoid TPA losses. In these systems, the magnitude of Re(γ) is limited because little pre-resonant enhancement can be achieved without incurring substantial TPA losses (ħω has to be ] 2𝑛 𝑐 𝜖0 ℏ

(2)

where L denotes a local field factor; n, the refractive index; c, the speed of light in vacuum; and 𝜖0 , the permittivity of free space. Since typical formulations of the local field factor for solvated chromophores can lead to systematic errors,55 here we focus on simulations in vacuum, and thus L and n are set to 1. For more background to these formulas, we refer to the discussion in the Supplementary Information.

3. Computational details

The geometries of the four model streptocyanines with three, five, seven, and nine carbon atoms along the backbone and two dimethylamino end groups (Figure 1), are optimized using a density functional theory (DFT) approach with the long-range corrected B97x functional56 and Pople's 6-31G** basis set.57,58 While most standard DFT functionals over-delocalize the ground-state wavefunction in π-conjugated systems, the long-range corrected functionals typically provide a more accurate description of the ground-state wavefunction and degree of bond-length alternation (BLA);45 the B97x functional in particular has been shown to provide reliable geometries.59 We have recently shown that LRC functionals are in excellent agreement with CCSD(T) geometries for small polymethines,45 and the geometries used here are consistent with those previously used.20 We recall that the streptocyanines under investigation are closedshell singlet cations, with a C2v molecular point group symmetry.

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

Page 10 of 38

Using these geometries, the excited-state properties were computed at several levels of theory. The ADC computations were performed with Pople's 6-31G and 6-31G* basis sets as well as Dunning's cc-pVDZ basis set.60 Since ADC(2)-s calculations with the aug-cc-pVDZ basis set could not be completed due to convergence problems, the slightly smaller cc-pVDZ++ basis set, which excludes d-type diffuse functions on heavy atoms, was used instead. The difference between double and triple zeta basis sets is expected to be rather limited for the smaller closedshell cation compounds, while the computational demands become prohibitive towards the larger compounds.61-63 In this study, we do not aim to describe excited states with a profound Rydberg character. Focusing upon a comparison between different methods, the chosen moderate basis sets bear insight into the relevant →* excitations of the polymethines.

The electronic structure and absorption spectra of the streptocyanines have also been investigated using the symmetry adapted cluster configuration interaction (SAC-CI) SD-R method64,65 with the 6-31G* basis set. In our study, the maximum excitation level of the linked excitation operators was two. Neither the 9, 11, 13 and 15 core orbitals nor the corresponding number of uppermost virtual orbitals for 3C, 5C, 7C and 9C, respectively, are included in the active space.

Calculations have in addition been performed using the Intermediate Neglect of Differential Overlap (INDO) Hamiltonian66,67 with multireference determinant single- and doubleexcitation configuration interaction (MRDCI). Five reference determinants have been included, as well as an active space of 20 orbitals.21,68 The standard Mataga and Nishimoto parameterization69 was applied to express the Coulomb repulsion term; for the sake of comparison, the Ohno parameterization70 was also used. The nature of the excited states

10

ACS Paragon Plus Environment

Page 11 of 38

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

computed with INDO/MRDCI has been visualized with the ZOA program.71 Finally, RICC2/TZVP and CCSD/6-31G* calculations have also been performed.

For the ADC(2)-s, ADC(2)-x and ADC(3/2) calculations, a development version of the Q-Chem package has been used,72 while Gaussian09 has been employed for the SAC-CI calculations and the geometry optimizations.73 The INDO/MRDCI calculations were performed using the ZINDO code.74 The RI-CC2 calculations were performed with the Turbomole program 6.3.1,75 while the CCSD calculations were carried out using Dalton 2013.76,77

For the third-order polarizability , 300 excited states have been taken into account in the Sumover-States expression at the INDO/MRDCI level. For SAC-CI, as well as for the ADC(2)-s and ADC(2)-x levels, 5 excited states have been used. For ADC(3/2), the number of states has been truncated to the one-photon allowed first excited state and the next two two-photon allowed ones. Only the lowest excited states are found to be prominent contributors in the sumover-states expression of equation (1) through their coupling to the first excited state e, which is in itself strongly coupled to the ground state. In the static limit or at optical frequencies of light, this expression converges quickly because the denominator contains a product of three energy differences. In addition, a response theory approach has been used to compute the TPA cross-sections at the ADC(2)-s and ADC(2)-x levels. To compute the response-theory values, the two-photon absorption cross-section is expressed in terms of an inversion of the ADC matrix (see equation S4 in SI); the resulting equations are solved using iterative solvers for systems of linear equations and the DIIS (Direct Inversion of the Iterative Subspace) formalism as detailed in the SI.78

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

Page 12 of 38

4. Results and discussion

4.1. Excited-state energies

As can be inferred by examining the essential-state model for γ, an accurate evaluation of the NLO properties of polymethine dyes requires a proper evaluation of, at least, the lowest two excited states. The nature of the lowest excited states in streptocyanines is very similar for all computational approaches considered: the first excited state e has B2 symmetry and primarily involves a H → L transition, while the second and third excited states, e’ and e”, have A1 symmetry and involve linear combinations of the H-1 → L, H → L+1, and H2 → L2 electronic transitions. All of the orbitals involved in these transitions, see Figure 2, are π orbitals that extend across the entire length of the molecule. The relative contributions of these transitions to states e’ and e” vary among the levels of theory under consideration (see Table 1 for 5C and Table S1 in the Supplementary Information, SI, for data concerning the other three compounds). At the ADC(2)-x, the ADC(3/2), and INDO/MRDCI levels, the e’ state has a significant doubleexcitation character. For ADC(2)-s and CC2, these contributions are largely absent, which is not unexpected since these methodologies only incorporate zeroth-order orbital differences on the diagonal of the 2p-2h block. The e” state still has a considerable amount of doubleexcitation character, although it is less significant than in the e’ state. Among the methods studied, only SAC-CI yields a comparable amount of double excitation character for states e’ and e”.

12

ACS Paragon Plus Environment

Page 13 of 38

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

Figure 2. Frontier molecular orbitals of the four streptocyanines under consideration (HF/631G* data).

Table 1: Analysis and decomposition of the 𝒆, 𝒆′ and 𝒆" excited states for 𝟓𝑪, obtained at different levels of theory and the 6-31G* basis set. 𝐚 𝑒: 0.65 |𝐻 → 𝐿〉 𝑒′: 0.64 |𝐻 − 1 → 𝐿〉 𝑒": 0.64 |𝐻 → 𝐿 + 1〉 ADC(2)-x 𝑒: 0.64 |𝐻 → 𝐿〉 𝑒′: 0.47 |𝐻 − 1 → 𝐿〉 + 0.26 |𝐻 2 → 𝐿2 〉 𝑒" b : 0.50 |𝐻 → 𝐿 + 1〉 ADC(3/2) 𝑒: 0.64 |𝐻 → 𝐿〉 𝑒′: 0.42 |𝐻 − 1 → 𝐿〉 + 0.31 |𝐻 2 → 𝐿2 〉 + 0.27 |𝐻 → 𝐿 + 1〉 𝑒": 0.48 |𝐻 → 𝐿 + 1〉 + 0.41 |𝐻 − 1 → 𝐿〉 CC2 𝑒: 0.98 |𝐻 → 𝐿〉 𝑒′: 0.91 |𝐻 − 1 → 𝐿〉 + 0.39 |𝐻 → 𝐿 + 1〉 𝑒": 0.86 |𝐻 → 𝐿 + 1〉 + 0.44 |𝐻 − 1 → 𝐿〉 SAC-CI 𝑒: 0.94 |𝐻 → 𝐿〉 𝑒′: 0.66 |𝐻 − 1 → 𝐿〉 + 0.22 |𝐻 2 → 𝐿2 〉 + 0.15 |𝐻 → 𝐿 + 1〉 𝑒": 0.91 |𝐻 → 𝐿 + 1〉 + 0.22 |𝐻 − 1 → 𝐿〉 + 0.17 |𝐻 2 → 𝐿2 〉 INDO/MRDCI 𝑒: 0.90 |𝐻 → 𝐿〉 𝑒′: 0.33 |𝐻 − 1 → 𝐿〉 + 0.66|𝐻 → 𝐿 + 1〉 + 0.54|𝐻 2 → 𝐿2 〉 𝑒": 0.75 |𝐻 − 1 → 𝐿〉 + 0.49|𝐻 → 𝐿 + 1〉 ADC(2)-s

a

𝐻 and 𝐿 denote the highest occupied and the lowest unoccupied orbitals, respectively. The INDO/MRDCI results are obtained with the standard Mataga-Nishimoto parameterization. The INDO/MRDCI orbitals are found to be identical to the ones obtained by HF/6-31G*. b At the ADC(2)-x level, the 𝑒" excited state consists of 0.48|𝐻 → 𝐿 + 1〉 + 0.11|𝐻2 → 𝐿2 〉 for 7𝐶, while it amounts to 0.46|𝐻 → 𝐿 + 1〉 + 0.12|𝐻2 → 𝐿2 〉 for 9𝐶.

We now turn to the excited-state energies (see Table S2) and first focus on the energy Ege of state e, which decreases as molecular length and conjugation increase (Figure 3). We compare 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 38

the computed excited-state energies to the experimental absorption maxima in dilute solution.16,79 Although in general comparing the computed Ege to the experimental absorption spectra requires accounting for vibronic effects,80-82 the streptocyanine absorption spectra have narrow peaks where the 0-0 vibronic peaks dominate due to the non-bonding character of the frontier molecular orbitals and the resulting small geometric changes between the ground and excited states.3,16,46 Thus, a direct comparison of Ege to the absorption maxima is expected to be reasonable for these systems. We note that experimental absorption spectra have been measured for dilute solutions of the dyes in dichloromethane, methanol, and trifluoroethanol, and the absorption maxima vary by no more than 0.1 eV among these solvents; test calculations on 3C and 5C taking some of these solvents into account at ADC(2)-s/ptSS-PCM(PTD*) level83 along with the 6-31G* basis set show negligible shifts with respect to the gas phase values.

Among the methods under consideration, substantial variation in the Ege values is observed. Whereas CC2 and CCSD overestimate the experimental Ege by 0.4-0.5 eV as has been previously observed,8 INDO/MRDCI dramatically underestimates Ege, particularly for the shortest streptocyanines. SAC-CI yields more accurate Ege values, although the decrease in Ege with increasing molecular length is steeper than is seen experimentally in contrast with the other methods considered. Among the ADC approaches, the ADC(3/2) and ADC(2)-s results are within about 0.2 eV of the experimental Ege (4.01 [3.97], 3.02 [3.01], 2.44 [2.44] and 2.03 [1.96] eV for the experimental absorption maxima in methanol [trifluoroethanol])8,16,79, whereas ADC(2)-x underestimates Ege by about 0.5-0.6 eV. The Ege values obtained using polarized basis functions are up to 0.1 eV higher than those obtained using unpolarized ones, especially for SAC-CI and ADC(3/2). From these results, the conclusion can be drawn that a more complete treatment of double excitations does not necessarily result in a more accurate Ege: since the first excited state e can be described within the singly excited manifold, the accuracy 14

ACS Paragon Plus Environment

Page 15 of 38

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

of the ADC approaches does not correlate with the treatment of double excitations. Furthermore, the accuracy is not significantly improved when going from CC2 to CCSD.

Figure 3. Comparison between the lowest excited-state energies (Ege) of the four cyanines, as obtained by various methods (in eV). Solid lines correspond to the 6-31G* basis set, dashed lines correspond to the 6-31G basis set, and dotted lines correspond to the TZVP basis set. The experimental data in trifluoroethanol are taken from Refs 16 and 17.

Because of experimental challenges, to the best of our knowledge, the second excited state energy Ege’ has only been measured for 3C and 7C.46,47 Since the TPA of state e’ is close to the onset of OPA into state e, it is hard to observe the full vibronic structure of the TPA peak. However, the relatively sharp onset of the TPA peaks46,47 and the non-bonding character of the orbitals involved in state e’ suggest that the absorption is sufficiently narrow that a direct 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

Page 16 of 38

comparison between Ege’ and the TPA peak is valid. For the 7C molecule, the ADC(3/2) excitation energy is within 0.3 eV of the experimental value; in the 3C molecule, however, the value is overestimated by more than 1 eV. An inspection of the energy differences Ege’ [Ege”] between the ground state and the second [third] excited state e’ [e”] reveals the relative overestimation of the CC2 energies (see Figure 4), which is consistent with the overestimation of Ege at this level.

Figure 4: The comparison between the (a) second (Ege’) and (b) third excited state (Ege”) energies of the four cyanines obtained by various methods (in eV). Solid lines correspond to the 6-31G* basis set and dashed lines correspond to the 6-31G basis set. The experimental data in trifluoroethanol are taken from Ref. 46. 16

ACS Paragon Plus Environment

Page 17 of 38

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

As was seen for Ege, the INDO/MRCI second and third excited-state energies exhibit not only the lowest values but also a different slope. The slopes of the ADC(3/2) and ADC(2)-x excitedstate energy curves are found to be largely parallel. With respect to ADC(2)-s and CC2, however, the slopes of the former higher-order methods are found to be steeper with the elongation of the cyanine backbone. Especially for Ege’, the differences between the SAC-CI values and the ADC(2)-x results increase with the length of the cyanine backbone. This is a consequence of the large double-excitation character of the A1 states, which is somewhat more pronounced for e’ than for e”. In the case of Ege”, this effect is therefore less pronounced when considering the energy differences among the different ADC(2)-s, SAC-CI, and ADC(3/2) results for the two longest compounds. This finding is in line with our work on linear polyacenes, where it was found that the evolution of the 21Ag state energy is much steeper at the ADC(2)-x level than at the ADC(2)-s level, owing to the pronounced double excitation character of this state in the longer polyacenes.37

As described in the Introduction, the properties of the second excited state, in particular the energy gap between the first and second excited states, are key in evaluating the third-order NLO properties of the specific compound. The ratio Ege’/Ege of the first two excited-state energies is experimentally observed to be 1.71 for the 7C streptocyanine; ratios in the range 1.67-1.73 have also been measured for polymethines with thiazolium (5-9 carbons) and indolium (7-11 carbons) end groups. The consistency in this ratio across a range of lengths and end groups suggests that this ratio should be applicable to the longer streptocyanines. For shorter polymethines, the ratio may be smaller, with values of 1.45 observed for the 3C streptocyanine and 1.55 for the 5C polymethine with indolium end groups.

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

Page 18 of 38

At the INDO/MRDCI level, Ege’/Ege is quite large, particularly for the longer cyanines; in fact, Ege’ is underestimated relative to the experimental values despite this large ratio that comes about because Ege is much smaller than the experimental value. Although SAC-CI predicts the closest ratio compared to experiment for the 3C streptocyanine, the ratios for the larger cyanines are substantially overestimated due to the pronounced underestimation of Ege. In contrast, CC2 provides Ege’/Ege values close to 1.7 for the range of lengths considered but only partially captures the experimental decrease in the ratio for the shortest cyanines.

Figure 5. The ratio of the first and second excited state energies of the four cyanines. Solid lines correspond to the 6-31G* basis set and dashed lines correspond to the 6-31G basis set. The shown experimental values in trifluoroethanol are taken from Ref. 46; for longer cyanine molecules, the experimentally extrapolated ratio is ~1.7 (see Ref. 9).

Among the ADC approaches, there is likewise significant variation in the Ege’/Ege ratios. ADC(3/2)/6-31G* provides ratios very close to 1.7 value for the shorter cyanines, whereas both ADC(2) approaches overestimate the ratio. The ADC(2)-s ratios start to plateau for longer cyanines, while the ADC(2)-x values continue to increase nearly linearly. Within each scheme,

18

ACS Paragon Plus Environment

Page 19 of 38

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 ratios are smaller (and closer to 1.7) when polarization functions are considered: Ege slightly increases and Ege’ slightly decreases upon inclusion of polarization functions.

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

Page 20 of 38

4.2. Excited-state properties The third-order NLO properties depend critically on the transition dipole moments and state dipole moments among the lowest several streptocyanine electronic states (see equation (1)). The difference Δμeg between the state dipole moments of states e and g appears in the D term; due to the C2v symmetry, the streptocyanine Δμeg can only be non-zero along the short axis of the π system. As Δμeg is negligibly small, we will focus on the transition dipole moments μge, which dominate the N term, and the transition dipole moments between excited states e and e’ (μee’) and e” (μee”), which are important for the T term.

The transition dipole moment μge from state g to state e is roughly the same at all levels of calculation for each molecule. The absolute value of the x-component is  7.6 Debye for 3C,  10.7 Debye for 5C,  13.4 Debye for 7C, and 16.0 Debye for 9C at the ADC(3/2) level. The INDO/MRCI and CC2 methods slightly overestimate the μge values relative to ADC(3/2), while SAC-CI underestimates them especially for the longest compounds. Since this transition predominantly involves a H → L excitation and the spatial distributions of the HOMO and LUMO wavefunctions are consistent among all levels of theory, the small variations in μge among the applied methods are unsurprising (see Figure 6a and Table S2).

20

ACS Paragon Plus Environment

Page 21 of 38

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

Figure 6. Transition dipole moments (in Debye) from: (a) the ground state g to the first excited state e; (b) between the first excited state e and the second excited state e’; and (c) between the first excited state e and the third excited state e”, for the four cyanines. Solid lines correspond to the 6-31G* basis set and dashed lines correspond to the 6-31G basis set. 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

Page 22 of 38

From Figure 6 (b), it can be seen that at all levels of theory, ee’ increases with increasing streptocyanine length. The ee’ values strongly depend on the level of theory and, in particular, on the treatment of double excitations. The methods that do not describe double excitations beyond zeroth order, like SAC-CI and CC2, provide generally large values of ee’, whereas INDO/MRDCI, which treats single and double excitations on equal footing within the chosen active space, yields much smaller ee’ values and a very weak length dependence. A similar pattern can be observed when comparing the ADC methodologies. Within the ADC series, the value of ee’ sharply decreases at all molecular lengths as the treatment of double excitations improves from ADC(2)-s to ADC(2)-x and ADC(3/2). These observations are consistent with our recent work demonstrating that the double excitation has a critical role in decreasing the magnitude of ee’.21 Since the magnitude of Im(γ) and thus the TPA cross-section of the second excited state strongly depend on ee’ as will be detailed in the following section, these results highlight the importance of adequately treating double excitations to accurately model the streptocyanine NLO properties.

The transition dipole moments obtained at the CC2 level are determined through first order, while the ones based upon the ADC(2)-s, ADC(2)-x and ADC(3/2) methodologies are calculated through second order in perturbation theory. This is illustrated in Figure 6 (c), where CC2 is seen to exhibit a different pattern compared to both ADC(2) methods and to ADC(3/2). This evolution contrasts with that of the energies, which are consistently obtained through second order for CC2 as well as both ADC(2) methods; for ADC(3/2), the energy is obtained through third order. For the discussion of the third-order molecular polarizability  and the twophoton absorption cross-section, which strongly depend on the quality of the description of the 22

ACS Paragon Plus Environment

Page 23 of 38

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

first excited state dipole moments, it is thus more appropriate to focus on the performance of the ADC methods.

The transition dipole moments ee’ and ee” at the ADC(2)-s and ADC(2)-x levels exhibit the same characteristics, in the sense that both gradually increase as the streptocyanine backbone is extended (Figure 6). In contrast, the ADC(3/2) level, which gives a rather flat curve for ee’ around ~2.5-3 Debye, gives rise to a considerably steep curve and higher values when ee” is considered for the longer cyanines. The values for ee’ at the ADC(3/2) level are considerably smaller than both the ADC(2)-s and ADC(2)-x ones; however, for ee”, the ADC(3/2) results are intermediate.

4.3 Third-order polarizabilities (;-,,)

We now turn to a discussion of the third-order polarizabilities of the streptocyanines. We first consider the static (ħω = 0) values. As Im(γstat) is negligible, we focus here on Re(γstat) (see Table 2). Since the SOS expressions contain a product of four (state or transition) dipole moments in the numerator and three energies in the denominator, the value for Re(γstat) is expectedly very sensitive to the accuracy of the excited-state properties described in the previous section. Re(γ) converges rapidly with the number of excited states included in the summation within the SOS formalism. For ADC, differences due to varying the number of excited states (5 or 3, see computational details) are only of marginal importance. A consideration of 5 excited states at the INDO/MRDCI level introduces no more than a 25% error relative to a summation over 300 excited states. With an increasing molecular length, this error decreases and does not change the trends in the magnitude of Re(γ). This convergence is 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

Page 24 of 38

due to the form of the SOS expression and the unique excited-state structure of the polymethines. Because the denominator consists of a product of three energy differences, the magnitudes of the contributions to Re(γstat) decrease rapidly with increasing excited-state energy. Since the first excited state is energetically well-separated from the higher excited states and is more strongly optically allowed than any of the other low-lying states, its contribution typically dominates the magnitude of Re(γstat) and only a few two-photon-absorbing states have to be considered.

To rationalize the magnitude of Re(γstat), we consider the magnitude of the terms appearing in the essential-state model (Equation (1)). As we mentioned earlier, due to the C2v molecular symmetry, the state dipole moments can only be non-zero along the short axis of the π system. Thus, there is little difference in state dipole moments among the different states and the D term is negligibly small. We note that in longer polymethines where ground-state symmetrybreaking occurs, the D term has a much more significant contribution to Re(γ).44,84 For polymethines, the N term typically dominates; as μge falls within a narrow range for all of the levels of theory considered, the magnitude of the N term is primarily dependent on Ege. Since μge increases and Ege decreases with increasing molecular length, the evolution of the N term causes Re(γstat) to increase by more than two orders of magnitude as the molecular backbone is extended from three to nine carbons. The T term is usually smaller in magnitude than the N term; the relative magnitudes of these two terms depend on the energies and transition dipole moments of the TPA-active excited states. We note that the experimental Re(γstat) value is small and positive for the 3C streptocyanine with a perchlorate anion in trifluoroethanol.46 However, among our calculations, all approaches give negative Re(γstat) values, and only the SAC-CI value is comparable in magnitude to the experimental value (see Table 2); the difference in sign may be related to the essential-state approximation, which considers only a handful of excited 24

ACS Paragon Plus Environment

Page 25 of 38

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

states; we also neglect factors such as vibronic effects and the distribution of geometries present experimentally. For streptocyanine 5C, all methods give the same negative sign, but only the SAC-CI value is on the same order of magnitude as the experimental one.

Table 2: Real parts of the static third-order polarizability [𝜸𝒔𝒕𝒂𝒕 ] for the 3C, 5C, 7C and 9C cyanines. 𝐚 ADC(2)-s ADC(2)-x CC2 SAC-CI ADC(3/2) INDO/MRDCI Experiment c

6-31G*

3C -31 -81 -43 -2 -45 -141 +5

5C -274 -827 -267 -34 -404 -1360 b -53

7C -1337 -4957 -795 -565

9C -4582 -21302 +2443 -3333

-6885 -1224

-23609 -5082

a

𝛾𝑠𝑡𝑎𝑡 is given in 10−36 esu. All experimental and computed values are reported using a Taylor series expansion for 𝛾, which yields values a factor of 6 larger than the power series expansion commonly used in experimental studies.3,49 b Standard (Mataga-Nishimoto) parametrisation; the values for the Ohno parametrization is: 𝛾𝑠𝑡𝑎𝑡 = −384 ⋅ 10−36 esu. c In trifluoroethanol; taken from Ref. 46.

The relative magnitudes of Re(γstat) at the various levels of theory (Table 2) can be rationalized in terms of the excited-state properties discussed in the previous Section. The magnitude of Ege plays a particularly significant role in determining the magnitude of Re(γstat): INDO/MRDCI and ADC(2)-x underestimate Ege and thus substantially overestimate Re(γstat) relative to experiment. This difference is primarily due to the enhancement of the N term when Ege is small. The accuracy of description of the higher-lying excited-state properties has a smaller but still noticeable effect via the T term. For example, when comparing Re(γstat) for the 3C cyanine at the ADC(2)-s/6-31G and ADC(3/2)/6-31G levels, the first excited-state properties and thus the magnitudes of the N term are similar. However, since μee’ and μee” are several times larger

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

Page 26 of 38

at the ADC(2)-s/6-31G level, the T term is larger and overall Re(γstat) is less negative in the latter case.

Especially for the larger cyanines 7C and 9C, the CC2-calculated Re(stat) values deviate severely from the experimental ones and the ones obtained with all other considered methods. The unphysical, large positive value for 9C is a direct consequence of the dramatic overestimation of ee”. It can be seen that the quality of the transition dipole moments among the excited states e, e’, and e” are decisive for the third-order polarizability: the CC2 energies Ege, Ege’ and Ege” largely agree with the ADC(2)-s ones; however, the transition dipole moments and the Re(stat) values are markedly different.

4.4

Two-photon absorption

As described in the Introduction, the TPA cross-section δ at a photon energy ħω is directly proportional to Im(γ) at that energy. Im(γ) is negligibly small except at energies such that 2ħω approaches resonance with a (two-photon allowed) excited state. At these energies, the denominators (𝐸ge − 2ℏ𝜔 − 𝑖Γge ) in the D term or (𝐸ge′ − 2ℏ𝜔 − 𝑖Γge′ ) in the T term become vanishingly small. Here, we calculate the TPA cross-section δ for the lowest streptocyanine excited states (for each state x, at an energy ħω = 𝐸gx /2).

Because δ is exquisitely sensitive to the detuning of ħω from 𝐸ge , even slight deviations in the ratio 𝐸ge′ /𝐸ge can lead to substantial errors in δ; we recall that 𝐸ge′ /𝐸ge is expected to be close to 1.7 across a range of polymethine lengths. As a result, when 𝐸ge′ /𝐸ge is very close to 2, the computed δ is dramatically enhanced, which is not physically meaningful; if this ratio was 26

ACS Paragon Plus Environment

Page 27 of 38

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

observed experimentally, one-photon absorption would dominate the signal and TPA would be unobservable. Therefore, the exceptionally large δ values obtained using several of the computational approaches for the polymethines of 5 carbons and longer must be interpreted with care (see Tables 3-4). The chosen line broadening value  of 0.1 eV may also affect the accuracy; as the experimental OPA peaks have widths close to 0.1 eV and the TPA peaks are likewise comparably narrow as detailed in the previous Section,9,47 this value of  is expected to be physically reasonable. Since the computed values of δ span several orders of magnitude, small inaccuracies in  will not significantly affect the comparisons made here.

Table 3: Results of the two-photon absorption calculations obtained using response theory and the three-state model for 3C and 5C. 𝐚 Method Exc state ADC(2)-s/6-31G* 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 ADC(2)-x/6-31G* 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 ADC(3/2)/6-31G 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 ADC(3/2)/6-31G* 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 CC2/6-31G* 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 SAC-CI/6-31G* 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 INDO/MRDCI b 𝑒′ 𝐴1 𝑒" 𝐴1 a

3C Energy Resp. theo Three-state TPA Exc. state 4.18 (0.946) 1 0 𝑒 𝐵1 7.05 (0.037) 1186 1188 𝑒′ 𝐴1 7.38 (0.007) 842 849 𝑒" 𝐴1 3.46 (0.758) 1 0 𝑒 𝐵1 5.81 (0.027) 430 436 𝑒′ 𝐴1 54 56 6.57 c (0.000) 𝑒" 𝐴1 4.14 (0.914) 0 𝑒 𝐵1 6.92 (0.036) 168 𝑒′ 𝐴1 7.51 (0.001) 64 𝑒" 𝐴1 4.24 (0.907) 0 𝑒 𝐵1 6.89 (0.032) 115 𝑒′ 𝐴1 7.45 (0.001) 23 𝑒" 𝐴1 4.48 (1.014) 0 𝑒 𝐵1 7.23 (0.042) 117 𝑒′ 𝐴1 7.56 (0.000) 5 𝑒" 𝐴1 4.18 (0.829) 0 𝑒 𝐵1 6.21 (0.025) 613 𝑒′ 𝐴1 7.44 (0.012) 536 𝑒" 𝐴1 2.84 (0.674) 0 𝑒 𝐵1 4.84 (0.016) 94 𝑒′ 𝐴1 711 5.90 c (0.015) 𝑒" 𝐴1

5C Energy Resp. theo Three-state TPA 3.24 (1.424) 2 0 5.81 (0.028) 11168 11168 125647 6.36 c (0.005) 125651 2.51 (1.066) 1 0 4.53 (0.018) 3606 3638 5.46 (0.000) 105 3.12 (1.361) 0 5.49 (0.023) 774 9506 6.32 c (0.000) 3.24 (1.353) 0 5.44 (0.020) 296 2598 6.29 c (0.000) 3.54 (1.593) 0 5.97 (0.030) 2945 6.54 (0.000) 5479 3.10 (1.161) 0 5.01 (0.020) 4959 6.35 c (0.006) 39416 2.10 (0.990) 0 15796 4.11 c (0.015) 4.84 (0.002) 257

The excitation energies are given in eV, along with their oscillator strengths (in parentheses); the results of the TPA are given in Göppert-Mayer (GM). For 3C, at the CASPT2 level of theory and an ANO-S basis set with a contraction scheme of 4s3p1d/2s, for 𝑒 [ 1 𝐵2 ] and 𝑒′ [ 1 𝐴1 ] energies of 4.17 eV and 6.12 eV are reported. For 5C, using the same method and basis set, for 𝑒 [ 1 𝐵2 ] and 𝑒′ [ 1 𝐴1 ] energies of 2.95 eV and 4.68 eV are reported.19 For the e’ state of 3C, an experimental TPA cross section of 10 GM has been reported in trifluoroethanol.46 b Standard (Mataga-Nishimoto) parametrisation. c One half of the energy of this state is relatively close to the energy of the first excited state (ratio 1.9-2.1).

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

Page 28 of 38

Table 4: Results of the two-photon absorption calculations obtained using response theory and the three-state model for 7C and 9C. 𝐚 Method

Exc state

ADC(2)-s/6-31G*

𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1

ADC(2)-x/6-31G*

ADC(3/2)/6-31G

CC2/6-31G*

SAC-CI/6-31G* INDO/MRDCI b a

7C Resp. theo Three-state TPA Exc state

Energy 2.60 4.80 5.41 c 1.87 3.58 c 4.43 2.40 4.39 5.24 2.89 4.98 5.59 c 2.24 4.20 5.31 1.66 3.49 c 4.11

(1.789) (0.012) (0.015) (1.242) (0.011) (0.002) (1.651) (0.016) (0.001) (2.095) (0.019) (0.007) (1.291) (0.013) (0.013) (1.221) (0.015) (0.001)

3 53391 93496 2 50734 399

0 53522 93405 0 50976 0 3785 4672 0 13662 305380 0 157058 1542 0 4768 1024

𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1 𝑒 𝐵1 𝑒′ 𝐴1 𝑒" 𝐴1

9C Resp. Three-state TPA theo (2.111) 4 0 (0.003) 202831 203480 (0.027) 47642 47531 (1.353) 3 0 (0.006) 886215 886126 (0.007) 966 (1.863) 0 (0.012) 28253 (0.005) 8166 (2.588) 0 (0.010) 13153 (0.019) 442346102 (1.362) 0 (0.007) 444088 (0.027) 1318 (1.435) 0 (0.013) 3342 (0.008) 2863

Energy 2.17 4.09 4.73 1.44 2.93 c 3.69 1.91 3.63 c 4.44 2.46 4.28 4.90 c 1.71 3.61 4.56 1.38 3.06 3.64

The excitation energies are given in eV, along with their oscillator strengths (in parentheses); the results of the TPA are given in Göppert-Mayer (GM). For 7C, at the CASPT2 level of theory and an ANO-S basis set with a contraction scheme of 4s3p1d/2s, for 𝑒 [ 1 𝐵2 ] and 𝑒′ [ 1 𝐴1 ] energies of 2.32 eV and 3.56 eV are reported. For 9C, using the same method and basis set, for 𝑒 [ 1 𝐵2 ] and 𝑒′ [ 1 𝐴1 ] energies of 1.93 eV and 2.93 eV are reported.19 For the e’ state of 7C, an experimental TPA cross section of 557 GM has been reported in trifluoroethanol.46 b Standard (Mataga-Nishimoto) parametrisation. c One half of the energy of this state is relatively close to the energy of the first excited state (ratio 1.9-2.1).

Experimentally, the TPA values of the e’ states of the 3C and 7C systems have been reported as 10 GM and 557 GM, respectively.46 However, the computed values consistently overestimate the TPA cross-sections, in many cases by multiple orders of magnitude. This is due to the limitations of these methodologies in predicting the ratios of the excited-state energies and the transition dipole moments between excited states. Although ADC(3/2) still overestimates both of the experimentally reported values by roughly an order of magnitude, its reasonable performance for both excited-state energies and transition dipole moments makes it currently the best performer among the methods studied. In contrast, ADC(2)-s, ADC(2)-x, SAC-CI, and 28

ACS Paragon Plus Environment

Page 29 of 38

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

CC2 all yield unphysically large TPA cross-sections in the majority of cases due to their limitations discussed in the previous Section. Interestingly, INDO/MRDCI provides relatively small TPA cross-sections for all of the cyanines except for 5C; although this approach correctly predicts small transition dipole moments among the excited states, the Ege’/Ege ratio is unphysically larger than 2 for the longer streptocyanines. These results highlight once again the significant challenges in accurately modeling the polymethine TPA properties, as minor inaccuracies in the excited-state characteristics can lead to substantial changes in the accuracy of the NLO properties.

For states e’ and e”, results from the simplified three-state model and the response theory approach are consistently in excellent agreement, which demonstrates that the TPA into these states is predominantly due to coupling through the first excited state as a virtual state. Due to the excited-state spacing in the polymethines, which have a low lying one-photon allowed state and a strongly coupled higher-lying second excited state which is two-photon allowed, the SOS expression converges much more quickly with the number of excited states in the summation than has previously been observed for ethylene using ADC/ISR approaches.38 As shown in Tables S4 and S5, the excellent agreement between the three-state model and response theory approach holds even for excited states up to several eV higher in energy than states e’ and e”.

5. Conclusions and outlook

Accurate modeling of the polymethine excited-state properties, in particular, energies and transition dipole moments among the first several excited states, is critical to evaluate their third-order NLO properties. Here, we have focused on the first several excited states of four streptocyanine molecules with varying backbone length and the resulting description of the 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

Page 30 of 38

two-photon absorption and the third-order polarizability . In particular, we compared the strict and extended second-order Algebraic Diagrammatic Construction schemes along with the intermediate state representation to a few benchmark ADC(3/2) calculations to investigate the influence of doubly excited configurations. These results have also been compared to other computational methods like CC2, SAC-CI SD-R, and the semi-empirical INDO/MRDCI algorithm, which have previously been used to study polymethine NLO properties. For the simulation of two-photon absorption spectra, an essential three-state model has been compared with a methodology based on matrix inversion and the direct inversion of the iterative subspace (DIIS) algorithm; the results are found to be very comparable and validate the sum-over-states expression, which quickly converges.

The linear and non-linear optical properties of polymethines are notoriously hard to simulate as the energies as well as the characteristics of the lowest excited states need to be properly described. Highly accurate methods must be applied to obtain chemically relevant results. In our analysis, ADC(3/2) has been found to be the most appropriate method that is consistently and reasonably accurate across all properties. Although the Re() calculated by means of a limited essential-state formalism somewhat overestimates the experimental values, the ADC(3/2) TPA cross sections are found to be more reasonable than those from any of the other methodologies studied. However, its computational cost prohibits application to larger polymethines that are of experimental interest.

Comparison of the ADC(3/2) results to the ADC(2)-s and ADC(2)-x schemes reveals that each of the latter two methods has strengths and marked weaknesses. For an appropriate choice of the second-order ADC methodology, it is important to consider the relevant properties of interest. ADC(2)-x is known to underestimate excited-state energies; however, it has a clear 30

ACS Paragon Plus Environment

Page 31 of 38

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

advantage with respect to ADC(2)-s in the description of the transition dipole moments ee’ and

ee” because of its more complete treatment of doubly excited states. Thus, Re() is better described by ADC(2)-s, since its magnitude depends primarily on the first excited-state properties.

Among the non-ADC methodologies, the excited-state energies obtained by CC2 are consistent with ADC(2)-s and show an Ege’/Ege ratio very close to 1.7, but the transition dipole moments between excited states are greatly overestimated due to the limited treatment of double excitations; thus, CC2, which describes response properties only through first order, performs poorly for both Re() and TPA (i.e., Im()). The properties calculated by SAC-CI SD-R are consistently among the less accurate, while the excited-state energies decrease too rapidly with increasing polymethine length; to describe streptocyanines, other methods bearing a similar computational cost could therefore be selected. We note that ADC(2)-x, CCSD and SAC-CI SD-R scale like N6, while the computational effort for CC2 and ADC(2)-s follows a N5 scaling, with N the number of orbitals. INDO/MRDCI, which allows for a proper treatment of double excitations at a very low computational cost, bears reasonable transition dipole moments and a Re() which shows an acceptable trend; however, the excited-state energies are only moderately reliable.

The benchmarking performed here can be used to inform the choice of computational methodologies to understand the properties of new polymethines of experimental interest. In addition, since the methods described here are still limited in their accuracy, this study suggests that further work to accurately compute the nonlinear optical properties of polymethines is needed. One natural extension is an improved description of the third-order polarizability  for 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

Page 32 of 38

the various ADC schemes by taking into account transition dipole moments between the ground state and higher-lying excited states. Further on, complex response theory can be considered, as it has been done for applications in Hartree-Fock or DFT calculations with the conjugate residual and preconditioned iterative subspace algorithms,85 as well as in coupled cluster theory with the complex polarization propagator approaches.86

32

ACS Paragon Plus Environment

Page 33 of 38

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

Acknowledgments S.K. is grateful to the Georgia Institute of Technology for their hospitality in Fall 2012. The authors thank Dr. Paul Winget for computational assistance. J.L.B. acknowledges generous support of this work by KAUST.

Supplementary information The Supporting Information is available free of charge on the ACS Publications website at DOI:… Description of the calculation of excitations, TPA spectra and third-order polarizabilities by means of ADC; decomposition of the e, e’, and e” excited states for 3C, 7C, and 9C at various levels of theory; extended tables of the computed vertical excitation energies of the first excited states of the four cyanines, together with the transition dipole moments from the ground to the first, second, and third excited states; real parts of the static and dynamic third-order polarizabilities; results of the TPA calculations obtained using response theory and the threestate model for the four investigated compounds.

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

Page 34 of 38

References (1) Marder, S. R.; Gorman, C. B.; Meyers, F.; Perry, J. W.; Bourhill, G.; Brédas, J. L.; Pierce, B. M. Science 1994, 265, 632. (2) Meyers, F.; Marder, S. R.; Pierce, B. M.; Bredas, J. L. J. Am. Chem. Soc. 1994, 116, 10703. (3) Gieseking, R. L.; Mukhopadhyay, S.; Risko, C.; Marder, S. R.; Brédas, J.-L. Adv. Mater. 2014, 26, 68. (4) Hales, J. M.; Barlow, S.; Kim, H.; Mukhopadhyay, S.; Brédas, J.-L.; Perry, J. W.; Marder, S. R. Chem. Mater. 2014, 26, 549. (5) Hales, J. M.; Matichak, J.; Barlow, S.; Ohira, S.; Yesudas, K.; Bredas, J. L.; Perry, J. W.; Marder, S. R. Science 2010, 327, 1485. (6) Stegeman, G. I.; Stolen, R. H. J. Opt. Soc. Am. B 1989, 6, 652. (7) Le Guennic, B.; Jacquemin, D. Accounts of chemical research 2015, 48, 530. (8) Send, R.; Valsson, O.; Filippi, C. J. Chem. Theory Comput. 2011, 7, 444. (9) Fu, J.; Padilha, L. A.; Hagan, D. J.; Van Stryland, E. W.; Przhonska, O. V.; Bondar, M. V.; Slominsky, Y. L.; Kachkovskii, A. J. Opt. Soc. Am. B 2007, 24, 56. (10) Barlow, S.; Brédas, J. L.; Getmanenko, Y. A.; Gieseking, R. L.; Hales, J. M.; Kieu, K.; Kim, H.; Marder, S. R.; Norwood, R. A.; Perry, J. W.; Peyghambarian, N.; Risko, C.; Shahin, S.; Zhang, Y. Mater. Horiz. 2014, 1, 577. (11) Anderson, R. J. M.; Holtom, G. R.; McClain, W. M. J. Chem. Phys. 1979, 70, 4310. (12) Kohler, B. E.; Terpougov, V. J. Chem. Phys. 1996, 104, 9297. (13) Petek, H.; Bell, A. J.; Choi, Y. S.; Yoshihara, K.; Tounge, B. A.; Christensen, R. L. J. Chem. Phys. 1993, 98, 3777. (14) Salem, L. (15) Jacquemin, D. J. Phys. Chem. A 2011, 115, 2442. (16) Nikolajewski, H. E.; Dahne, S.; Leupold, D.; Hirsch, B. Tetrahedron 1968, 24, 6685. (17) Konig, W.; Regner, W. Ber. Dtsch. Chem. Ges. 1930, 63, 2823. (18) Angeli, C.; Pastore, M. J. Chem. Phys. 2011, 134, 184302. (19) Schreiber, M.; Buß, V.; Fülscher, M. P. Phys. Chem. Chem. Phys. 2001, 3, 3906. (20) Mukhopadhyay, S.; Risko, C.; Marder, S. R.; Brédas, J.-L. Chem. Sci. 2012, 3, 3103. (21) Gieseking, R. L.; Mukhopadhyay, S.; Risko, C.; Brédas, J.-L. ACS Photon. 2014, 1, 261. (22) Karaca, S.; Elmacı, N. J. Mol. Struct.: THEOCHEM 2009, 915, 149. (23) Jacquemin, D.; Zhao, Y.; Valero, R.; Adamo, C.; Ciofini, I.; Truhlar, D. G. J. Chem. Theory Comput. 2012, 8, 1255. (24) Dreuw, A.; Head-Gordon, M. Chem. Rev. 2005, 105, 4009. (25) Starcke, J. H.; Wormit, M.; Schirmer, J.; Dreuw, A. Chem. Phys. 2006, 329, 39. (26) Haettig, C. Adv. Quant. Chem. 2005, 50, 37. (27) Saitow, M.; Mochizuki, Y. Chem. Phys. Lett. 2012, 525, 144. (28) Wormit, M.; Rehn, D. R.; Harbach, P. H. P.; Wenzel, J.; Krauter, C. M.; Epifanovsky, E.; Dreuw, A. Mol. Phys. 2014, 112, 774. (29) Mertins, F.; Schirmer, J. Phys. Rev. A 1996, 53, 2140. (30) Mertins, F.; Schirmer, J.; Tarantelli, A. Phys. Rev. A 1996, 53, 2153. (31) Dreuw, A.; Wormit, M. WIREs Comput Mol Sci 2015, 5, 82. (32) Schirmer, J. Phys. Rev. A 1982, 26, 2395. (33) Trofimov, A. B.; Schirmer, J. J. Phys. B 1995, 28, 2299. (34) Trofimov, A. B.; Stelter, G.; Schirmer, J. J. Chem. Phys. 1999, 111, 9982. (35) Schirmer, J.; Trofimov, A. B. J. Chem. Phys. 2004, 120, 11449. 34

ACS Paragon Plus Environment

Page 35 of 38

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

(36) Harbach, P. H.; Wormit, M.; Dreuw, A. The Journal of chemical physics 2014, 141, 064113. (37) Knippenberg, S.; Starcke, J. H.; Wormit, M.; Dreuw, A. Molecular Physics 2010, 108, 2801. (38) Knippenberg, S.; Rehn, D. R.; Wormit, M.; Starcke, J. H.; Rusakova, I. L.; Trofimov, A. B.; Dreuw, A. The Journal of chemical physics 2012, 136, 064107. (39) Wenzel, J.; Wormit, M.; Dreuw, A. Journal of chemical theory and computation 2014, 10, 4583. (40) Wenzel, J.; Holzer, A.; Wormit, M.; Dreuw, A. The Journal of chemical physics 2015, 142, 214104. (41) Tolbert, L. M.; Zhao, X. J. Am. Chem. Soc. 1997, 119, 3253. (42) Lepkowicz, R. S.; Przhonska, O. V.; Hales, J. M.; Fu, J.; Hagan, D. J.; Van Stryland, E. W.; Bondar, M. V.; Slominsky, Y. L.; Kachkovski, A. D. Chem. Phys. 2004, 305, 259. (43) Hu, H.; Przhonska, O. V.; Terenziani, F.; Painelli, A.; Fishman, D.; Ensley, T. R.; Reichert, M.; Webster, S.; Bricks, J. L.; Kachkovski, A. D.; Hagan, D. J.; Van Stryland, E. W. Phys. Chem. Chem. Phys. 2013, 15, 7666. (44) Terenziani, F.; Przhonska, O. V.; Webster, S.; Padilha, L. A.; Slominsky, Y. L.; Davydenko, I. G.; Gerasov, A. O.; Kovtun, Y. P.; Shandura, M. P.; Kachkovski, A. D.; Hagan, D. J.; Van Stryland, E. W.; Painelli, A. J. Phys. Chem. Lett. 2010, 1, 1800. (45) Gieseking, R. L.; Ravva, M. K.; Coropceanu, V.; Brédas, J.-L. J. Phys. Chem. C 2016, DOI: 10.1021/acs.jpcc.6b02100. (46) Johr, T.; Wernke, W.; Pfeiffer, M.; Lau, A.; Dahne, L. Chem. Phys. Lett. 1995, 246, 521. (47) Wernke, W.; Pfeiffer, M.; Johr, T.; Lau, A.; Grahn, W.; Johannes, H.-H.; Dahne, L. Chem. Phys. 1997, 216, 337. (48) Ryabitsky, A. B.; Kachkovski, A. D.; Przhonska, O. V. J. Mol. Struct.: THEOCHEM 2007, 802, 75. (49) Orr, B. J.; Ward, J. F. Mol. Phys. 1971, 20, 513. (50) Garito, A. F.; Teng, C. C. Proc. SPIE 1986, 0613, 146. (51) Helfin, J. R.; Garito, A. F. In Electroresponsive Molecular and Polymeric Systems; Skotheim, T. A., Ed.; Marcel Dekker Inc.: New York, 1991, p 1. (52) Pierce, B. M. Proc. SPIE 1991, 1560, 148. (53) Dirk, C. W.; Cheng, L.-T.; Kuzyk, M. G. Int. J. Quantum Chem. 1992, 43, 27. (54) Brédas, J.-L.; Adant, C.; Tackx, P.; Persoons, A.; Pierce, B. M. Chem. Rev. 1994, 94, 243. (55) Wortmann, R.; Bishop, D. M. J. Chem. Phys. 1998, 108, 1001. (56) Chai, J.-D.; Head-Gordon, M. J. Chem. Phys. 2008, 128, 084106. (57) Hariharan, P. C.; Pople, J. A. Theor. Chim. Acta 1973, 28, 213. (58) Huzinaga, S. Comp. Phys. Rep. 1993, 2, 281. (59) Remya, K.; Suresh, C. H. Journal of computational chemistry 2013, 34, 1341. (60) Dunning, T. H. J. Chem. Phys. 1989, 90, 1007. (61) Szalay, P. G.; Aquino, A. J. A.; Barbatti, M.; Lischka, H. Chemical Physics 2011, 380, 9. (62) Knippenberg, S.; Deleuze, M. S. Journal of Electron Spectroscopy and Related Phenomena 2010, 178-179, 61. (63) Marian, C. M.; Gilka, N. J. Chem. Theor. Comput. 2008, 4, 1501. (64) Nakatsuji, H. Chem. Phys. Lett. 1978, 59, 362. (65) Nakatsuji, H.; Hirao, K. J. Chem. Phys. 1978, 68, 2053. (66) Ridley, J.; Zerner, M. Theor. Chim. Acta 1973, 32, 111. 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

Page 36 of 38

(67) Ridley, J. E.; Zerner, M. C. Theor. Chim. Acta 1976, 42, 223. (68) Gieseking, R. L.; Risko, C.; Bredas, J. L. The journal of physical chemistry letters 2015, 6, 2158. (69) Mataga, N.; Nichimoto, K. Z. Phys. Chem. 1957, 13, 140. (70) Ohno, K. Theor. Chim. Acta 1964, 2, 219. (71) Zoa v3.0, J. P. Calbert, Service de Chimie des Matriaux Nouveaux, Universit de Mons, Mons (Belgium) - http://zoa.freeservers.com. (72) Y. Shao, L. F. Molnar, Y. Jung, J. Kussmann, C. Ochsenfeld, S. T. Brown, A. T. Gilbert, L. V. Slipchenko, S. V. Levchenko, D. P. O’Neill, R. A. D. Jr, R. C. Lochan, T. Wang, G. J. Beran, N. A. Besley, J. M. Herbert, C. Y. Lin, T. V. Voorhis, S. H. Chien, A. Sodt, R. P. Steele, V. A. Rassolov, P. E. Maslen, P. P. Korambath, R. D. Adamson, B. Austin, J. Baker, E. F. C. Byrd, H. Dachsel, R. J. Doerksen, A. Dreuw, B. D. Dunietz, A. D. Dutoi, T. R. Furlani, S. R. Gwaltney, A. Heyden, S. Hirata, C.-P. Hsu, G. Kedziora, R. Z. Khalliulin, P. Klunzinger, A. M. Lee, M. S. Lee, W. Liang, I. Lotan, N. Nair, B. Peters, E. I. Proynov, P. A. Pieniazek, Y. M. Rhee, J. Ritchie, E. Rosta, C. D. Sherrill, A. C. Simmonett, J. E. Subotnik, H. L. W. III, W. Zhang, A. T. Bell, and A. K. Chakraborty, “Advances in methods and algorithms in a modern quantum chemistry program package,” Phys. Chem. Chem. Phys. 8, 3172–3191 (2006). (73) M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, . Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, and D. J. Fox, “Gaussian 09 Revision B.01,” Gaussian Inc. Wallingford CT 2009. (74) M. C. Zerner, G. H. Loew, R. F. Kichner, and U. T. Mueller-Westerhoff, J. Am. Chem. Soc. 102, 589 (1980). (75) TURBOMOLE V6.3.1 2011, a development of University of Karlsruhe and Forschungszentrum Karlsruhe GmbH, 1989-2007, TURBOMOLE GmbH, since 2007; available from http://www.turbomole.com. (76) K. Aidas, C. Angeli, K. L. Bak, V. Bakken, R. Bast, L. Boman, O. Christiansen, R. Cimiraglia, S. Coriani, P. Dahle, E. K. Dalskov, U. Ekström, T. Enevoldsen, J. J. Eriksen, P. Ettenhuber, B. Fernández, L. Ferrighi, H. Fliegl, L. Frediani, K. Hald, A. Halkier, C. Hättig, H. Heiberg, T. Helgaker, A. C. Hennum, H. Hettema, E. Hjertenæs, S. Høst, I.-M. Høyvik, M. F. Iozzi, B. Jansik, H. J. Aa. Jensen, D. Jonsson, P. Jørgensen, J. Kauczor, S. Kirpekar, T. Kjærgaard, W. Klopper, S. Knecht, R. Kobayashi, H. Koch, J. Kongsted, A. Krapp, K. Kristensen, A. Ligabue, O. B. Lutnæs, J. I. Melo, K. V. Mikkelsen, R. H. Myhre, C. Neiss, C. B. Nielsen, P. Norman, J. Olsen, J. M. H. Olsen, A. Osted, M. J. Packer, F. Pawlowski, T. B. Pedersen, P. F. Provasi, S. Reine, Z. Rinkevicius, T. A. Ruden, K. Ruud, V. Rybkin, P. Salek, C. C. M. Samson, A. Sánchez de Merás, T. Saue, S. P. A. Sauer, B. Schimmelpfennig, K. Sneskov, A. H. Steindal, K. O. Sylvester-Hvid, P. R. Taylor, A. M. Teale, E. I. Tellgren, D. P. Tew, A. J. Thorvaldsen, L. Thøgersen, O. Vahtras, M. A. Watson, D. J. D. Wilson, M. Ziolkowski, and H. Ågren, "The Dalton quantum chemistry program system", WIREs Comput. Mol. Sci. (doi: 10.1002/wcms.1172) 36

ACS Paragon Plus Environment

Page 37 of 38

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

(77) Dalton, a molecular electronic structure program, Release DALTON2013.4 (2013), see http://daltonprogram.org. (78) Pulay, P. Chem. Phys. Lett. 1980, 73, 393. (79) Gürtler, O.; Dähne, S. Z. Phys. Chem., Leipzig 1974, 255, 501. (80) Laurent, A. D.; Le Guennic, B.; Jacquemin, D. Theor. Chim. Acta 2016, 135, 173. (81) Azarias, C.; Jacquemin, D. J. Phys. Chem. A. 2016, 120, 2824. (82) Fihey, A.; Favennec, A.; Le Guennic, B.; Jacquemin, D. Phys. Chem. Chem. Phys. 2016, 18, 9358. (83) Mewes, J. M.; You, Z. Q.; Wormit, M.; Kriesche, T.; Herbert, J. M.; Dreuw, A. The journal of physical chemistry. A 2015, 119, 5446. (84) Terenziani, F.; Painelli, A.; Katan, C.; Charlot, M.; Blanchard-Desce, M. J. Am. Chem. Soc. 2006, 128, 15742. (85) Kauczor, J.; Jørgensen, P.; Norman, P. J. Chem. Theory Comput. 2011, 7, 1610. (86) Coriani, S.; Christiansen, O.; Fransson, T.; Norman, P. Phys. Rev. A 2012, 85, 022507.

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

TOC 221x155mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 38 of 38