Electronically Excited Solute Described by RISM ... - ACS Publications

9 Oct 2018 - with Multireference Perturbation Theory: Vertical Excitation. Energies of Bioimaging Probes. Ryosuke Y. Shimizu,. †. Takeshi Yanai,. â€...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Quantum Electronic Structure

Electronically Excited Solute Described by RISM Approach Coupled with Multi-Reference Perturbation Theory: Vertical Excitation Energies of Bio-Imaging Probes Ryosuke Yale Shimizu, Takeshi Yanai, Yuki Kurashige, and Daisuke Yokogawa J. Chem. Theory Comput., Just Accepted Manuscript • DOI: 10.1021/acs.jctc.8b00599 • Publication Date (Web): 09 Oct 2018 Downloaded from http://pubs.acs.org on October 12, 2018

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

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

Electronically Excited Solute Described by RISM Approach Coupled with Multi-Reference Perturbation Theory: Vertical Excitation Energies of Bio-Imaging Probes Ryosuke Y. Shimizu,† Takeshi Yanai,†,‡,¶ Yuki Kurashige,§ and Daisuke Yokogawa∗,†,‡ †Department of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan ‡Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan ¶Japan Science and Technology Agency, PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan §Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan E-mail: [email protected] Phone: +81 52 789 2851

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 For theoretically studying molecules with fluorescence in the near-infrared region, high-accuracy determination of state energy level is required for meaningful analyses since the spectra of interest are of very narrow energy range. In particular, these molecules are in many cases handled in solution; therefore, consideration of the solvation effect is essential upon calculation together with the electronic structure of the excited state. Earlier studies showed that they cannot be described with conventional methods such as PCM-TD-DFT, yielding results far from experimental data. Here, we have developed a new method by combining a solvation theory based on statistical mechanics (RISM) and a multi-reference perturbation theory (CASPT2) with the extension of the density matrix renormalization group (DMRG) reference states for calculating the photochemical properties of near-infrared molecules and have obtained higher-accuracy prediction.

1

Introduction Development of new fluorescent molecules have enabled the visualization of various bio-

logical events and have supported advances in the field of bio-imaging. 1–3 Ideal fluorescence for these probes have been within the visible to near-infrared (NIR) region since biological tissue can be efficiently penetrated, allowing deep tissue imaging in many higher organisms. 4 While conventional NIR molecules have had fluorescence in the so-called “first NIR window,” approximately 700-1000 nm, recently, molecules with fluorescence in the “second NIR window,” approximately 1000-1700 nm, have captured attention; these molecules are expected to further improve the performance of bio-imaging, yet are underexplored. 5,6 Considering the difficulty of measuring transition energies experimentally, theoretical calculations hold a great advantage in accelerating the study of these molecules by providing complementary or direct description of the excitation. Nonetheless, since these molecules possess photochemical properties in a narrow energy range, highly accurate calculations are required for 2

ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27 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

meaningful analyses. Accomplishing this task requires achieving high accuracy in the calculation of both the electronic structure of the excited solute and the solvation effect. One of the most widely used methods for calculating excited states is the time-dependent density functional theory (TD-DFT). 7,8 Although TD-DFT can evaluate large systems at an affordable cost, it often has uncontrollable difficulties in correctly assessing the electronic structure of long π-conjugated systems, 9–11 a general characteristic of the fluorescent molecules in focus. In terms of solvation, the polarizable continuum model (PCM) 12,13 is widely employed. Again, its computational cost is reasonable; however, this method is known to be unable to illustrate local solvation structures such as hydrogen bonding well. 14 Although the combination of these methods, PCM-TD-DFT, 15 is widely employed for calculating excited states in solution, it is reportedly incapable of providing the desired accuracy of the aforementioned type of molecules. 16,17 Hence, there has been a marked challenge in developing a method to rectify the limitations of the conventional approaches and acquire higher accuracy with reasonable computational cost. In this article, we present the development of a new method by combining the complete active space second-order perturbation theory (CASPT2) 18,19 and the reference interaction site model (RISM) 20 for the models of excitation and solvation, respectively. The CASPT2 method is a multi-reference perturbation theory, which is capable of accurately describing correlated many-electron wave functions at the second-order perturbation level using the multi-configurational states determined by a complete active space self-consistent field method (CASSCF) 21 calculation as references. Limitations in choosing the active space for CAS calculations have recently been greatly eased with the use of the density matrix renormalization group (DMRG) in DMRG-CASSCF 22–24 and DMRG-CASPT2. 25,26 The RISM method is a solvation theory of a statistical mechanical approach, in which the solvation structure is represented by correlation functions. Since the correlation functions are obtained by solving the RISM equations analytically, the calculation is efficient compared to methods

3

ACS Paragon Plus Environment

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

Page 4 of 27

such as QM/MM-MD which demands large computational cost for sampling based on MD simulations. Furthermore, information of the structure of the solvent molecule is preserved in RISM, facilitating the description of local interactions such as hydrogen bonding. 27,28 These advantages have prompted the combination of RISM with various QM approaches (RISM-SCF 29–33 and RISM-SCF-SEDD 34–36 ).

2

Method In combining RISM and CASPT2, we first defined the Helmholtz free energy of the solute

in solution calculated with the CASPT2 method as, 31,32

ˆ |Ψ⟩ + ∆µ + E2 ACASPT2 ≡ ⟨Ψ| H = ACASSCF + E2 ,

(1)

ˆ is the Hamiltonian constructed in the gas phase, Ψ is the CAS reference where the operator H function, ∆µ is the excess chemical potential, and E2 is the second-order energy. Singer and Chandler’s hypernetted-chain (HNC) closure relation 37 was employed for the excess chemical potential term given as follows: [ ] ∫ 1 2 1 1∑ ρs dr cαs (r) − hαs (r) + hαs (r)cαs (r) . ∆µ = − β α,s 2 2

(2)

Greek subscripts refer to the interaction sites of the solute and the Roman subscripts refer to the sites of the solvent, ρs is the number density of the solvent, and β is 1/kB T , where kB is the Boltzmann constant and T is temperature. The functions hαs and cαs are the total and direct correlation functions, respectively, and are obtained by solving the RISM equations. 20 The Helmholtz free energy ACASSCF defined in eq 1 is regarded as a functional of the correlation functions hαs , cαs , and tαs = hαs − cαs , and the CAS reference function Ψ.

4

ACS Paragon Plus Environment

Page 5 of 27 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

This free energy is minimized by defining the following Lagrangian: 31

L = A [c, h, t, Ψ] − E (⟨Ψ|Ψ⟩ − 1) .

(3)

Taking the variation of the Lagrangian with respect to the functions yields, [ ] ∫ 1∑ ∂∆µ ∂∆µ ∂∆µ ρs dr δL = − δtαs (r) + δhαs (r) + δcαs (r) β α,s ∂tαs ∂hαs ∂cαs ˆ + Vt ∂p − E |Ψ⟩ = 0, + ⟨δΨ| H ∂ ⟨Ψ| where Vi = −

∑ s

and

∫ ∫ qs ρs

hαs (r)fi (r′ ) drdr′ ′ |r − r |

(4)

(i ∈ α),

∑ ∂p [Ξ + (m − 1)Γ]−1 R′pq a†p aq . = ∂ ⟨Ψ| pq

(5)

(6)

The operators a†p and aq are creation and annihilation operators respectively, f is the auxiliary basis set (ABS), and the remaining matrices are as defined in previous work. 17,38 From the Ornstein-Zernike equation, HNC closure relation, and the RISM equation, the terms in the square brackets in eq 4 become 0. Hence, eq 4 simplifies to,

ˆ solv − E |Ψ⟩ = 0, ⟨δΨ| H

(7)

ˆ solv as, where we have defined the solvated Hamiltonian H ˆ solv = H ˆ + Vt ∂p . H ∂ ⟨Ψ|

5

ACS Paragon Plus Environment

(8)

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 6 of 27

RISM calculations are performed with the electron density distribution on each atomic site. The distributions are obtained from the one-particle reduced density matrix defined as, γpq = ⟨Ψ| a†p aq |Ψ⟩ .

(9)

For molecules with large active spaces which require the DMRG, the reference wave function is simply replaced by the active-space DMRG wave function |ΨDMRG ⟩. 25,26 The absorption energy ω a was then calculated as the difference of the excited state energy [E]

[G]

ECASPT2 and the free energy at the ground state ACASPT2 : 17 [E]

[G]

ω a = ECASPT2 − ACASPT2 ,

(10)

where [G] ˆ [G] |Ψ⟩ + E2[G] + ∆µ[G] , ACASPT2 = ⟨Ψ| H

(11)

[E] ˆ [E] |Ψ⟩ + E2[E] + ∆µ[G] + (V[G] )t (p[E] − p[G] ). ECASPT2 = ⟨Ψ| H

(12) [G]

The superscripts G and E represent the ground and excited states, respectively. Vi

is

the electrostatic potential on the i-th solute site induced by the ground state solvation structure at the ground state, and p[G] and p[E] are the matrices including the expansion coefficients of the ABS for approximating the electron density of the ground and excited states, respectively. 17,38 For the absorption process, it is assumed that the relaxation of the solvation structure does not occur at the excited state following the Franck-Condon principle. 17 This can be observed in eq 12, where the sum of the first and second term is the energy of the solute at the excited state in the potential induced by the ground state solvation structure. Although it is well known that consideration of vibronic effects is important in computing absorption maxima, 39,40 we plan to implement these effects into our method in the future, and limit the present study to computing only vertical excitation energies.

6

ACS Paragon Plus Environment

Page 7 of 27 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 sketch of our implementation of the method is depicted in Scheme 1. For the ground state calculation, first Ψ is determined from ground state-specific (DMRG-)CASCI and orbital optimization steps. Then, the one-particle reduced density matrix γpq , formed from Ψ, is passed to the RISM code for computing the electrostatic potential Vi on each solute site induced by the solvation structure. These electrostatic potentials are used to construct the ˆ solv for the successive iteration. Once Ψ converges, the obtained solsolvated Hamiltonian H vated Hamiltonian is used in ground state-specific (DMRG-)CASPT2 to apply perturbation in solution. In computing the absorption energy, the potentials obtained from the ground state calculation, V[G] , are used to construct the solvated Hamiltonian for determining Ψ at the excited state; this corresponds to the calculation of the excited state with the solvation structure fixed at the ground state. Finally, perturbation is applied in solution for the excited state via state-specific (DMRG-)CASPT2 with the obtained solvated Hamiltonian. Since it is unnatural to include solvation effects on state-averaged electron density, all calculations are conducted state-specifically. START

START V[G]

,V

(DMRG-)CASCI

(DMRG-)CASCI

Orbital Optimization

Orbital Optimization

RISM Calculation

Update V ---

converged?

converged?

ˆ solv H

ˆ solv H

(DMRG-)CASPT2

(DMRG-)CASPT2

END

END

Scheme 1: Flowchart of the RISM-CASPT2 method for the ground (left) and excited (right) state in calculating absorption energies.

7

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

3

Computational Details The Lennard-Jones parameters used for the RISM calculations were taken from OPLS-

AA 41 and are provided in the supporting information together with solvent parameters. The solvent data required for the RISM calculations were computed using ex-RISM. 42,43 The basis sets employed were aug-cc-pVDZ 44 for nitrogen and oxygen, and cc-pVDZ 45 for the other atoms. The optimized geometries used for calculations involving RISM were computed using the CAM-B3LYP level of theory with symmetry employed at 300 K in respective solutions with RISM using the gamess program package 46,47 in which our method has been implemented in the past. The same program package was used for calculating absorption energies with RISM-TD-DFT. Calculations involving PCM employed optimized geometries computed with the CAM-B3LYP level of theory using PCM in respective solutions. The PCM-TD-DFT calculations were conducted with the gaussian program package 48 and PCM-MRMP2 calculations were conducted with the gamess program package. As written in the previous section, the reference wave functions for RISM-CASPT2 were determined with DMRG-CASSCF approach when using large active spaces. In this study, the active-space DMRG calculations were performed with 256 spin-adapted renormalized basis states. This size of renormalized basis produces near convegent DMRG results when accounting for quasi-one dimensional π conjugations, such as the systems considered in this study. The fourth-order cumulants were neglected in an approximate manner in the solution of the CASPT2 equation, for which the reduced density matrices of up to the third order were calculated based on the DMRG algorithm. 26 All calculations involved in the RISM-CASPT2 method and its extension with the active-space DMRG reference states were conducted with the orz program package 49 in which we have newly implemented our method into. Symmetry was not employed for all absorption energy calculations.

8

ACS Paragon Plus Environment

Page 8 of 27

Page 9 of 27 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

Results and Discussion As an assessment of the newly developed method, we have calculated the absorption

energies of the molecules given in Scheme 2. To evaluate the reliability of the solvation theory of choice and the accuracy of the electronic structure, we have selected MQZ for a molecule with high solvatochromism in absorption, 4AP as a molecule with absorption in the near UV region, and CY3 and ICG were chosen for long π-conjugated systems, related to bio-imaging probes. 3,5 All π and π ∗ orbitals were included in the active space. The size of the active space was (12e,11o) for MQZ and (14e,12o) for 4AP. For CY3 and ICG, the active spaces were set to (20e,19o) and (32e,31o), respectively; their configuration space exceeds the capability of the conventional CASSCF/CASPT2 method, but were allowed to be handled by its DMRG extension. The localized form of 31 active orbitals used as DMRG sites for ICG is displayed in the supporting information. All experimental absorption energies given refer to the absorption maxima of the molecules which are results under the presence of molecular vibrations; however, these effects are neglected in our calculations. Note that when using DMRG reference wave functions, RISM-CASPT2 calculations are denoted as RISM-DMRG-CASPT2.

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 27

O -O

NH N+

H2N O

N-Methyl-6-oxyquinolone (MQZ)

N+

4-Aminophthalimide (4AP)

N

Modeled indocarbocyanine (CY3)

N+

N

Modeled indocyanine green (ICG)

Scheme 2: Target molecules.

Calculated absorption energies of MQZ in water (WAT), methanol (MET), acetonitrile (ACN), and dimethyl sulfoxide (DMSO) are summarized in Table 1. For comparison, the results of PCM-TD-DFT calculations using the CAM-B3LYP 50 and PBE0 51 functionals are included. Comparing the TD-DFT level results with and without PCM, the solvation effect in the absorption energies was shown to be insensitive to the solvent polarity, constantly shifting gas-phase energies by an energy of approximately 1.0 eV and 0.9 eV for CAM-B3LYP and PBE0, respectively. Hence, the PCM-TD-DFT predictions showed no solvatochromism in absorption and were generally quantitatively away from the experimental absorption energies. The effect of the solvation model can be checked by using RISM in place of PCM at the TD-DFT level; this level of theory is referred to as RISM-TD-DFT. 17 We confirmed that the solvation energies obtained by RISM-TD-DFT vary depending on the solvents, showing an appreciable solvatochromic behavior in absorption. The use of RISM greatly improved the absorption energies with an error of less than 0.3 eV relative to the experimental data. 10

ACS Paragon Plus Environment

Page 11 of 27 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 RISM-TD-DFT results with the PBE0 functional especially showed good agreement with the experimental data; however, this is not always the case, and we will address on this later. In terms of the results from calculations employing multi-reference perturbation theories, using PCM as the solvation model again showed constant shifting of the gas-phase energies, not reproducing the solvatochromism in absorption. On the other hand, the RISM-CASPT2 method successfully illustrated solvatochromism in the absorption energies of MQZ, which were accurately predicted with an error of less than 0.2 eV relative to the experimental data. It was thus shown that the CASPT2 treatment generally further reduced the error in prediction, compared to those of the RISM-based TD-DFT results. Table 1: Calculated absorption energies of MQZ in various solutions (eV)

PCM-TD-DFT CAM-B3LYP PBE0 RISM-TD-DFT CAM-B3LYP PBE0 PCM-MRMP2 RISM-CASPT2 Exp.a a

Gas

ACN DMSO MET WAT

2.30 2.15

3.34 3.02

3.33 3.01

3.34 3.02

3.38 3.06

2.30 2.15 2.14 2.04 –

2.68 2.44 2.44 2.38 2.45

2.70 2.46 2.44 2.38 2.44

2.94 2.66 2.43 2.73 2.88

3.23 2.91 2.45 3.21 3.03

References 52–54.

Indeed, vibronic effects hold an important role in predicting absorption maxima. 39,40 However, these effects can be ignored to some extent for molecules with absorption in the near UV region where the vibronic contribution becomes relatively smaller when compared to the total absorption energy. Calculating the absorption energies of 4AP with the RISMCASPT2 method in the same solvents as MQZ showed similar accuracy with the error within approximately 0.2 eV relative to the experimental results (Table 2). These results provided further evidence of the accuracy of the method and proved that agreement with

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 27

the experimental results was not due to compensation of errors. Table 2: Calculated absorption energies of 4AP with RISM-CASPT2 in various solutions (eV)

RISM-CASPT2 Exp.a

Gas ACN DMSO MET WAT 3.80 3.62 3.55 3.54 3.40 3.67 3.49 3.32 3.36 3.36

a

Reference 55. The experimental value in hexane is used as a comparison to the calculations in the gas-phase.

One of the advantages of RISM is its ability to describe local solvation structures. To illuminate this point, we have computed the charge density ρq of the solvent around the solute site α by the following equation:

ρqα (r) = 4πr2



qs ρs gαs (r),

(13)

s

where gαs (r) is the radial distribution function representing the solvation structure. The charge density of the solvents around O- resulting from RISM-CASPT2 calculations of MQZ have been summarized in Figure 1. This solute site was chosen to focus on the degree of hydrogen bonding of each solvent. Results showed high positive charges present in the proximity of the O- site of MQZ in water. This indicates that there was strong hydrogen bonding between the solute site and solvent. Similarly, methanol also showed positive charges around the solute site, but was less than that of water. This corresponds to the solvation energy in methanol being less than that of water. Acetonitrile and DMSO showed nearly identical charge densities, and almost no positive charges existed around the solute site, which resulted in the small and identical solvation energies for the calculations in these solvents.

12

ACS Paragon Plus Environment

Page 13 of 27 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

WAT MET ACN DMSO

Figure 1: The charge density of the solvents, water, methanol, acetonitrile, and DMSO, around O- from the absorption energy calculations of MQZ by RISM-CASPT2.

As mentioned earlier, calculating absorption energies within the NIR region requires highly-accurate theoretical treatment. In order to highlight this point, the error in the calculated absorption energies compared to the experimental values 56 of CY3 and ICG, in DMSO and water, respectively, have been summarized in Figure 2. The absorption energies calculated by both PCM-TD-DFT and RISM-TD-DFT with the PBE0 functional had large errors, approximately 0.65 eV. They exceed acceptable errors for reliably studying the molecules of interest in the “second NIR window” (1000-1700 nm), 5 which spans for 0.5 eV. On the other hand, the RISM-DMRG-CASPT2 method predicted the absorption energies with the error reduced to within 0.2 eV, which seems to be satisfactory. These results illustrated the high applicability of RISM-CASPT2 in combination with the DMRG technique describing highly correlated π electronic states in solution.

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

Figure 2: Error in calculated absorption energies of CY3 and ICG by PCM-TD-DFT, RISMTD-DFT, and RISM-DMRG-CASPT2. The PBE0 functional was employed for TD-DFT calculations. Experimental values from Reference 56. The consideration of the multi-reference character seemingly played a key role in delivering the high accuracy. The natural orbital occupation number analysis based on the RISM-CASSCF wave function showed that the occupations of HOMO and LUMO digress from a simple high-spin open-shell character as the molecules become larger in size (Table 3). This means that the single-reference (or mean-field) treatment of excited states was acceptable for small systems such as MQZ, but inadequate in accounting for the energy levels of the larger systems such as CY3 and ICG. Table 3: Natural occupation numbers of HOMO and LUMO from RISM-CASSCF molecule HOMO LUMO MQZ 0.94 1.07 CY3 1.04 0.98 ICG 1.24 0.78

5

Conclusions To summarize, we have presented a new method for calculating excited states in solu-

tion by combining RISM and (DMRG-)CASPT2. Absorption energy calculations of MQZ, 14

ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27 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

CY3, and ICG with our new method successfully reproduced experimental measurements with satisfactory accuracy which suffices analysis of these bio-imaging related systems. The improvement was made by ameliorating the treatments of the excited electronic structure and the solvation effect simultaneously. The newly developed RISM-CASPT2 method and its extension with the DMRG reference are highly applicable for the theoretical study of the fluorescent molecules of interest within the “second NIR window.”

Acknowledgement This work was partly supported by JSPS KAKENHI Grant Numbers JP15K05385, JP16H06353, for Daisuke Yokogawa, JSPS KAKENHI Grant Numbers JP16H04101, JP15H01097 and JST PRESTO Grant Number JP17937609, for Takeshi Yanai, and JSPS KAKENHI Grant Number JP16H00855, for Yuki Kurashige.

Supporting Information Available The following files are available free of charge. • Lennard-Jones parameters and solvent parameters. • Localized active orbitals used as DMRG sites for ICG.

References (1) Wang, C.; Fukazawa, A.; Taki, M.; Sato, Y.; Higashiyama, T.; Yamaguchi, S. A Phosphole Oxide Based Fluorescent Dye with Exceptional Resistance to Photobleaching: A Practical Tool for Continuous Imaging in STED Microscopy. Angew. Chem. Int. Ed. 2015, 54, 15213–15217.

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

(2) Tsuchiya, Y.; Yoshimura, M.; Sato, Y.; Kuwata, K.; Toh, S.; Holbrook-Smith, D.; Zhang, H.; McCourt, P.; Itami, K.; Kinoshita, T.; Hagihara, S. Probing Strigolactone Receptors in Striga Hermonthica with Fluorescence. Science 2015, 349, 864–868. (3) Terai, T.; Nagano, T. Small-Molecule Fluorophores and Fluorescent Probes for Bioimaging. Pflügers Arch. 2013, 465, 347–359. (4) Valeur, B.; Berberan-Santos, M. N. Molecular Fluorescence, 2nd ed.; Wiley-VCH: Weinheim, 2012. (5) Hong, G.; Antaris, A. L.; Dai, H. Near-Infrared Fluorophores for Biomedical Imaging. Nat. Biomed. Eng. 2017, 1, 0010. (6) Tao, Z.; Hong, G.; Shinji, C.; Chen, C.; Diao, S.; Antaris, A. L.; Zhang, B.; Zou, Y.; Dai, H. Biological Imaging Using Nanoparticles of Small Organic Molecules with Fluorescence Emission at Wavelengths Longer than 1000 nm. Angew. Chem. Int. Ed. 2013, 52, 13002–13006. (7) Runge, E.; Gross, E. K. U. Density-Functional Theory for Time-Dependent Systems. Phys. Rev. Lett. 1984, 52, 997–1000. (8) Gross, E.; Kohn, W. In Density Functional Theory of Many-Fermion Systems; Löwdin, P.-O., Ed.; Advances in Quantum Chemistry; Academic Press, 1990; Vol. 21; pp 255–291. (9) Grimme, S.; Parac, M. Substantial Errors from Time-Dependent Density Functional Theory for the Calculation of Excited States of Large π Systems. ChemPhysChem 2003, 4, 292–295. (10) Fabian, J. TDDFT-Calculations of Vis/NIR Absorbing Compounds. Dyes Pigm. 2010, 84, 36–53.

16

ACS Paragon Plus Environment

Page 16 of 27

Page 17 of 27 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

(11) Kurashige, Y.; Yanai, T. Theoretical Study of the π → π ∗ Excited States of Oligoacenes: A Full π-Valence DMRG-CASPT2 Study. Bull. Chem. Soc. Jpn. 2014, 87, 1071–1073. (12) Cossi, M.; Barone, V.; Cammi, R.; Tomasi, J. Ab Initio Study of Solvated Molecules: A New Implementation of the Polarizable Continuum Model. Chem. Phys. Lett. 1996, 255, 327–335. (13) Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105, 2999–3094. (14) Nagy, P. I.; Alagona, G.; Ghio, C. Theoretical Investigation of Tautomeric Equilibria for Isonicotinic Acid, 4-Pyridone, and Acetylacetone in Vacuo and in Solution. J. Chem. Theory Comput. 2007, 3, 1249–1266. (15) Scalmani, G.; Frisch, M. J.; Mennucci, B.; Tomasi, J.; Cammi, R.; Barone, V. Geometries and Properties of Excited States in the Gas Phase and in Solution: Theory and Application of a Time-Dependent Density Functional Theory Polarizable Continuum Model. J. Chem. Phys. 2006, 124, 094107. (16) Charaf-Eddin, A.; Le Guennic, B.; Jacquemin, D. Excited-States of BODIPY-Cyanines: Ultimate TD-DFT Challenges? RSC Adv. 2014, 4, 49449–49456. (17) Yokogawa, D. Time-Dependent Density Functional Theory (TD-DFT) Coupled with Reference Interaction Site Model Self-Consistent Field Explicitly Including Spatial Electron Density Distribution (RISM-SCF-SEDD). J. Chem. Phys. 2016, 145, 094101. (18) Andersson, K.; Malmqvist, P. A.; Roos, B. O.; Sadlej, A. J.; Wolinski, K. Second-Order Perturbation Theory with a CASSCF Reference Function. J. Phys. Chem. 1990, 94, 5483–5488. (19) Andersson, K.; Malmqvist, P.; Roos, B. O. Second-Order Perturbation Theory with a

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

Complete Active Space Self-Consistent Field Reference Function. J. Chem. Phys. 1992, 96, 1218–1226. (20) Hansen, J.-P.; McDonald, I. R. Theory of Simple Liquids, 4th ed.; Academic Press: Oxford, 2013. (21) Roos, B. O.; Taylor, P. R.; Siegbahn, P. E. A Complete Active Space SCF Method (CASSCF) Using a Density Matrix Formulated Super-CI Approach. Chem. Phys. 1980, 48, 157–173. (22) Zgid, D.; Nooijen, M. The Density Matrix Renormalization Group Self-Consistent Field Method: Orbital Optimization with the Density Matrix Renormalization Group Method in the Active Space. J. Chem. Phys. 2008, 128, 144116. (23) Ghosh, D.; Hachmann, J.; Yanai, T.; Chan, G. K.-L. Orbital Optimization in the Density Matrix Renormalization Group, with Applications to Polyenes and β-Carotene. J. Chem. Phys. 2008, 128, 144117. (24) Kurashige, Y.; Yanai, T. High-Performance Ab Initio Density Matrix Renormalization Group Method: Applicability to Large-Scale Multireference Problems for Metal Compounds. J. Chem. Phys. 2009, 130, 234114. (25) Kurashige, Y.; Yanai, T. Second-Order Perturbation Theory with a Density Matrix Renormalization Group Self-Consistent Field Reference Function: Theory and Application to the Study of Chromium Dimer. J. Chem. Phys. 2011, 135, 094104. (26) Kurashige, Y.; Chalupský, J.; Lan, T. N.; Yanai, T. Complete Active Space SecondOrder Perturbation Theory with Cumulant Approximation for Extended Active-Space Wavefunction from Density Matrix Renormalization Group. J. Chem. Phys. 2014, 141, 174111.

18

ACS Paragon Plus Environment

Page 18 of 27

Page 19 of 27 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

(27) Yokogawa, D.; Ono, K.; Sato, H.; Sakaki, S. Theoretical Study on Aquation Reaction of cis-Platin Complex: RISM-SCF-SEDD, a Hybrid Approach of Accurate Quantum Chemical Method and Statistical Mechanics. Dalton Trans. 2011, 40, 11125–11130. (28) Sato, H.; Sakaki, S. Comparison of Electronic Structure Theories for Solvated Molecules: RISM-SCF Versus PCM. J. Phys. Chem. A 2004, 108, 1629–1634. (29) Ten-no, S.; Hirata, F.; Kato, S. A Hybrid Approach for the Solvent Effect on the Electronic Structure of a Solute Based on the RISM and Hartree-Fock Equations. Chem. Phys. Lett. 1993, 214, 391–396. (30) Ten-no, S.; Hirata, F.; Kato, S. Reference Interaction Site Model Self-Consistent Field Study for Solvation Effect on Carbonyl Compounds in Aqueous Solution. J. Chem. Phys. 1994, 100, 7443–7453. (31) Sato, H.; Hirata, F.; Kato, S. Analytical Energy Gradient for the Reference Interaction Site Model Multiconfigurational Self-Consistent-Field Method: Application to 1,2-Difluoroethylene in Aqueous Solution. J. Chem. Phys. 1996, 105, 1546–1551. (32) Higashi, M.; Kato, S. Theoretical Study on Electronic and Spin Structures of [Fe2 S2 ]2+,+ Cluster: Reference Interaction Site Model Self-Consistent Field (RISM-SCF) and Multireference Second-Order Møller-Plesset Perturbation Theory (MRMP) Approach. J. Phys. Chem. A 2005, 109, 9867–9874. (33) Yamazaki, S.; Kato, S. Solvent Effect on Conical Intersections in Excited-State 9HAdenine: Radiationless Decay Mechanism in Polar Solvent. J. Am. Chem. Soc. 2007, 129, 2901–2909. (34) Yokogawa, D.; Sato, H.; Sakaki, S. New Generation of the Reference Interaction Site Model Self-Consistent Field Method: Introduction of Spatial Electron Density Distribution to the Solvation Theory. J. Chem. Phys. 2007, 126 .

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

(35) Yokogawa, D.; Sato, H.; Sakaki, S. Analytical Energy Gradient for Reference Interaction Site Model Self-Consistent Field Explicitly Including Spatial Electron Density Distribution. J. Chem. Phys. 2009, 131, 214504. (36) Yokogawa, D. Coupled Cluster Theory Combined with Reference Interaction Site Model Self-Consistent Field Explicitly Including Spatial Electron Density Distribution. J. Chem. Theory Comput. 2018, 14, 2661–2666. (37) Singer, S. J.; Chandler, D. Free Energy Functions in the Extended RISM Approximation. Mol. Phys. 1985, 55, 621–625. (38) Yokogawa, D. New Fitting Approach of Electrostatic Potential for Stable Quantum Mechanical Calculations Using the Reference Interaction Site Model. Chem. Phys. Lett. 2013, 587, 113–117. (39) Avila Ferrer, F. J.; Cerezo, J.; Stendardo, E.; Improta, R.; Santoro, F. Insights for an Accurate Comparison of Computational Data to Experimental Absorption and Emission Spectra: Beyond the Vertical Transition Approximation. J. Chem. Theory Comput. 2013, 9, 2072–2082. (40) Santoro, F.; Jacquemin, D. Going Beyond the Vertical Approximation with TimeDependent Density Functional Theory. WIREs Comput Mol Sci 2016, 6, 460–486. (41) Jorgensen, W. L.; Maxwell, D. S.; Tirado-Rives, J. Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. (42) Andersen, H. C.; Chandler, D. Optimized Cluster Expansions for Classical Fluids. I. General Theory and Variational Formulation of the Mean Spherical Model and Hard Sphere Percus-Yevick Equations. J. Chem. Phys. 1972, 57, 1918–1929.

20

ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27 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

(43) Hirata, F.; Rossky, P. J. An Extended RISM Equation for Molecular Polar Fluids. Chem. Phys. Lett. 1981, 83, 329–334. (44) Kendall, R. A.; Dunning, Jr., T. H.; Harrison, R. J. Electron Affinities of the First-Row Atoms Revisited. Systematic Basis Sets and Wave Functions. J. Chem. Phys. 1992, 96, 6796–6806. (45) Dunning, Jr., T. H. Gaussian Basis Sets for use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. (46) Schmidt, M. W.; Baldridge, K. K.; Boatz, J. A.; Elbert, S. T.; Gordon, M. S.; Jensen, J. H.; Koseki, S.; Matsunaga, N.; Nguyen, K. A.; Su, S.; Windus, T. L.; Dupuis, M.; Montgomery, J. A. General Atomic and Molecular Electronic Structure System. J. Comput. Chem. 1993, 14, 1347–1363. (47) Gordon, M. S.; Schmidt, M. W. In Theory and Applications of Computational Chemistry; Dykstra, C. E., Frenking, G., Kim, K. S., Scuseria, G. E., Eds.; Elsevier: Amsterdam, 2005; pp 1167–1189. (48) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannen-

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

berg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09 Revision D.01. Gaussian Inc. Wallingford CT 2009. (49) Yanai, T.; Kurashige, Y.; Mizukami, W.; Chalupský, J.; Lan, T. N.; Saitow, M. Density Matrix Renormalization Group for Ab Initio Calculations and Associated Dynamic Correlation Methods: A Review of Theory and Applications. Int. J. Quantum Chem. 2015, 115, 283–299. (50) Yanai, T.; Tew, D. P.; Handy, N. C. A New Hybrid Exchange-Correlation Functional using the Coulomb-Attenuating Method (CAM-B3LYP). Chem. Phys. Lett. 2004, 393, 51–57. (51) Adamo, C.; Barone, V. Toward Reliable Density Functional Methods Without Adjustable Parameters: The PBE0 Model. J. Chem. Phys. 1999, 110, 6158–6170. (52) Schmode, S.; Ludwig, R. Utilization of the Dye N-Methyl-6-oxyquinolone as an Optical Acidometer in Molecular Solvents and Protic Ionic Liquids. Chem. Commun. 2017, 53, 10761–10764. (53) Kim, T. G.; Topp, M. R. Ultrafast Excited-State Deprotonation and Electron Transfer in Hydroxyquinoline Derivatives. J. Phys. Chem. A 2004, 108, 10060–10065. (54) Mason, S. F.; Philp, J.; Smith, B. E. Prototropic Equilibria of Electronically Excited Molecules. Part II. 3-, 6-, and 7-Hydroxyquinoline. J. Chem. Soc. A 1968, 3051–3056. (55) Mukherjee, S.; Sahu, K.; Roy, D.; Mondal, S. K.; Bhattacharyya, K. Solvation Dynamics of 4-Aminophthalimide in Dioxane-Water Mixture. Chem. Phys. Lett. 2004, 384, 128– 133. (56) Benson, R. C.; Kues, H. A. Absorption and Fluorescence Properties of Cyanine Dyes. J. Chem. Eng. Data 1977, 22, 379–383. 22

ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27 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

Graphical TOC Entry

23

ACS Paragon Plus Environment

START 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

Journal of Chemical Theory and Computation

START

Page 24 of 27

V[G]

,V

(DMRG-)CASCI

(DMRG-)CASCI

Orbital Optimization

Orbital Optimization

RISM Calculation

Update V ---

converged?

converged?

ˆ solv H

ˆ solv H

(DMRG-)CASPT2

END

ACS Paragon Plus Environment

(DMRG-)CASPT2

END

-

OChemical Theory and Computation Journal of

Page 25 of 27 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

N+

N-Methyl-6-oxyquinolone (MQZ)

N+

N

Modeled indocarbocyanine (CY3)

N+

N ACS Paragon Plus Environment

Modeled indocyanine green (ICG)

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

WAT MET ACN DMSO

ACS Paragon Plus Environment

Page 26 of 27

Page 27 of 27 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

Journal of Chemical Theory and Computation

ACS Paragon Plus Environment