Vibronic Theory of Ultrafast Intersystem Crossing Dynamics in a Single

Feb 2, 2016 - Miguel Paez-Espejo , Mouhamadou Sy , and Kamel Boukheddaden. Journal of the American Chemical Society 2016 138 (9), 3202-3210...
2 downloads 0 Views 1MB Size
Letter pubs.acs.org/JPCL

Vibronic Theory of Ultrafast Intersystem Crossing Dynamics in a Single Spin-Crossover Molecule at Finite Temperature beyond the Born−Oppenheimer Approximation Nikolay Klinduhov†,‡ and Kamel Boukheddaden*,† †

Groupe d’Etude de la Matière Condensée, Université de Versailles, CNRS UMR 8635, 45 Avenue des Etats-Unis, 78035 Versailles Cedex, France ‡ Institute of Technical Acoustics, National Academy of Sciences of Belarus, 13 Lyudnikova st, 210023 Vitebsk, Belarus S Supporting Information *

ABSTRACT: Quantum density matrix theory is carried out to study the ultrafast dynamics of the photoinduced state in a spin-crossover (SC) molecule interacting with a heat bath. The investigations are realized at finite temperature and beyond the usual Born−Oppenheimer (BO) approach. We found that the SC molecule experiences in the photoexcited state (PES) a huge internal pressure, estimated at several gigapascals, partly released in an “explosive” way within ∼100 fs, causing large bond length oscillations, which dampen in the picosecond time scale because of internal conversion processes. During this regime, the BO approximation is not valid. Depending on the tunneling strength, the ultrafast relaxation may proceed through the thermodynamic metastable high-spin state or prevent it. Interestingly, we demonstrate that final relaxation toward the low-spin state always follows a local equilibrium pathway, where the BO approach is valid. Our formulation reconciles the nonequilibrium and the equilibrium properties of this fascinating phenomenon and opens the way to quantum studies on cluster molecules.

P

(change in the metal−ligand distances and angular distortions of the ligand, essentially) accompanying the electronic relaxation create an elastic energy barrier which stabilizes the photoinduced HS metastable state at low temperature. This problem benefited recently of nice time-resolved14−17 femtosecond X-ray absorption, combined with optical spectroscopy measurements, which revealed that the photoswitching involves first molecular breathing vibrations of the ligand which dampen rapidly to the benefit of molecular bending vibrations. Theoretical investigations18,19 have also considered the quantum description of this phenomenon at very low temperatures within the BO approximation. Here, we provide a detailed study of the quantum dynamics of the multichannel relaxation of the photoinduced state, described in the frame of vibronic theory,20 beyond the BO approximation, whose validity is revisited when going from nonequilibrium properties to asymptotic thermodynamic equilibrium. We restrict our study to the isolated molecule in contact with a heat bath. One may legitimately ask about the benefit of going beyond the BO approximation in the description of the quantum dynamics of a vibronic system. First, it is important to recall that the BO approximation does not allow the accurate reproduction of the

ulsed ultrafast photoinduced processes are observed in a wide range of physical phenomena both for gases and at the condensed phase. Among them, one can quote (not an exhaustive list) the photoinduced ultrafast electron-transfer,1,2 femtosecond photomagnetism,3,4 ultrafast photochemical reactions,5,6 etc. Some of them are single-site processes, whereas others belong to cooperative ultrafast phenomena also called domino effect,7 such as photoinduced structural phase transitions.8 Fe(II)-based spin-crossover (SC) materials converting thermally from the diamagnetic low-spin (LS, S = 0, t62ge0g) to the paramagnetic high-spin (HS, S = 2, t42ge2g) states are model systems of photoswitching9−12 between these two states. A crucial issue for this work is the preparation of the initial excited state of the molecule. According to experimental data, the very short (femtoseconds) pump pulse promotes the spin-crossover molecule from the ground LS state into an excited metal-toligand charge-transfer (quintet MLCT) state via a Franck− Condon transition, from which the system rapidly decays in a nonadiabatic way through multiple excited-state channels (including ligand-field states)13 coupled with coherent vibrational motions and internal conversion. During this complex process, electronic and vibrational states are strongly interlinked, and in most of the cases, BO approximation fails. We simulate here this complex mechanism with the simple Franck−Condon transition in a double-well potential (see Figure 1). The local structural molecular reorganization © 2016 American Chemical Society

Received: January 4, 2016 Accepted: February 2, 2016 Published: February 2, 2016 722

DOI: 10.1021/acs.jpclett.6b00014 J. Phys. Chem. Lett. 2016, 7, 722−727

Letter

The Journal of Physical Chemistry Letters

Figure 1. (a) Adiabatic potential energies versus the molecular deformation with the exact eigenvalues (dotted lines) and their corresponding nuclear probability densities at T = 3 K for J = 20 meV. The vertical arrow indicates the schematic Franck−Condon process. (b) Thermal variation of the average magnetization, ⟨σ⟩, showing a gradual spin-crossover transition.

eigenvalues, denoted Ei, are displayed in Figure 1a (dotted lines). The associated vibronic eigenstates, which have the form

excited eigenvalues and eigenstates of a quantum system when there is an active vibronic coupling between these states, as for SC systems. This implies that a change in the nuclei configuration from Q to Q + ΔQ (Q is a configuration coordinate) causes electronic transitions. Depending on the type of treatment (BO or exact quantum resolution), the system will explore a different pathway in the space coordinates, leading to different internal dynamics including in the intersystem crossing region. Thus, strong vibronic interaction and/or strong tunneling mixing lead to the breakdown of the BO approach. Here, we refer the reader to relevant reviews treating the problem of the dynamics of the photoinduced cis− trans transition21 beyond the BO limit and to useful fundamental works based on the application of stochastic Monte Carlo methods for solving nonadiabatic quantum correlated systems at finite temperature.22,23 In the present work, the SC molecule is described as a vibronic oscillator whose Hamiltonian is24−26 Hmol =

kQ 2 ω 2P 2 + + ( Δ − SQ ) σ z + Jσ x 2k 2

|φi⟩ =

∑ (cniLS|χnLS ⟩|LS⟩ + cniHS|χnHS ⟩|HS⟩)

(2)

n

lead to the corresponding computed (see Supporting Information) nuclear probability densities depicted in Figure 1a. This allows us to derive (see Supporting Information) the thermal dependence of the fictitious magnetization, ⟨σ⟩, displayed in Figure 1b for J = 20 meV, showing a gradual 2δ spin conversion at the transition temperature, Teq = k ln g ∼ B

120 K. The thermal dependence of ⟨Q⟩ is given in Figure S1a of the Supporting Information. To help the reader to obtain a comprehensive view and deeper insights into the quantum resolution of this problem, we provide, in the Supporting Information, the BO solution of Hamiltonian 1 that we compare in Figure S1b to the exact quantum resolution. Next, we define the reservoir (thermal bath) as a set of independent harmonic oscillators, whose Hamiltonian is

(1)

Hph =

∑ ℏwqbq+bq (3)

q

where ω2p2/2k is the kinetic energy, Q the molecular distortion, 1 and k the harmonic elastic constant; Δ = δ − 2 kBT ln g is the effective ligand field energy splitting, accounting for the energy gap, δ, and entropic contribution of the HS/LS degeneracy ratio, g.27,28 The parameter S is the vibronic coupling at the S LS = −k origin of the existence of two equilibrium positions, Q eq

where bq and b+q are the usual annihilation and creation phonon operators and wq is the frequency of the qth normal mode. Restricting ourselves to the bilinear contributions in different coordinates, the molecule−bath Hamiltonian is Hmol−ph =

∑ (gqbq+a + gq*bqa+)

(4)

q

S

HS and Q eq = k , for the bond length distances (see Figure 1a), and J is the tunneling contribution originating from high-order spin−orbit coupling. σz and σx are the usual Pauli matrices of a two-states fictitious spin, with eigenvalues (respectively eigenstates) for σz, −1 and +1 (respectively, |LS⟩ and |HS⟩). The calculations are performed using realistic parameter values: LS k = 7 eV/Å2 and S = 0.8eV/Å leading to ΔQeq = QHS eq − Qeq ∼ 0.2 Å, δ ≅ 20 meV (∼200 K), ln g = 4.6, J ∼ 10−30 meV, and ℏω ∼ 30 meV. Hamiltonian 1, whose matrix elements are given in section 1 of the Supporting Information, is solved exactly beyond the BO HS approximation in the 2 × ∞ basis, (|χLS n (Q)⟩⊗|LS⟩, |χn (Q)⟩⊗| LS HS HS⟩), where |χn (Q)⟩ and |χn (Q)⟩ are the respective wave functions of the LS and HS displaced harmonic oscillators. The

where (a + a+) ∼ Q. The Hamiltonian Hmol−ph is rewritten in the diagonal basis (eq 2) of the isolated molecule, and treating the molecular−bath interaction as a perturbation, we could apply the usual Redfield theory for the density matrix.29,30 The density matrix of the system, ρ (already traced over bath q variables), is written in the interaction representation ⎛ iH t ⎞ ⎛ iH t ⎞ ρI (t ) = exp⎜ 0 ⎟ ρ(t )exp⎜ − 0 ⎟ ⎝ ℏ ⎠ ⎝ ℏ ⎠ ⎛ it = exp⎜⎜ ⎝ℏ

⎛ ⎞ ⎟ ∑ Ej|φj⟩⟨φj|⎟ ρ(t )exp⎜⎜− it j ⎝ ℏ ⎠

⎞ ∑ Ej|φj⟩⟨φj|⎟⎟ j ⎠ (5)

723

DOI: 10.1021/acs.jpclett.6b00014 J. Phys. Chem. Lett. 2016, 7, 722−727

Letter

The Journal of Physical Chemistry Letters

Figure 2. Time dependence of the quantum averages of molecular distortion, ⟨Q⟩ (left panel), and fictitious spin, ⟨σ⟩ (middle panel), for T = 3 K for two values of the tunneling constant: (a) J = 2 meV and (b) J = 20 meV. The horizontal blue lines plotted in the ⟨Q⟩ and ⟨σ⟩ panels are the thermal equilibrium values. The insets sketch the enlarged views of the dynamics at short time scale. Right panel: corresponding phase portraits ⟨σ⟩ versus ⟨Q⟩ showing the behavior of the system in its phase space.

the excited HS state.16,32 This fast first process is omitted in this work. Concretely, the initial excited state, |Ψ0⟩, is built from the vibronic ground state, |φ0⟩ (see eq 2), as

leading after mathematical developments (detailed in section 2 of the Supporting Information) to the system of differential equations for the molecular density matrix elements31 dρmm dt dρmj dt

=

∑ (γmkρkk

|Ψ0⟩ =

− γkmρmm )

n

k≠i

= = −(iωmj + Γmj)ρmj ,

in which, ωmk =

Em − Ek ℏ

m≠j

∑ (cnLS0 |χnLS ⟩ + cnHS0 |χnHS ⟩)|HS⟩ ∑ (∑ cniHScmLS0S(m , n) + cniHScnHS0 )|φi⟩ in

(6)

m

where S(n,m) are the usual Franck−Condon integrals, given by HS S(n,m) = ⟨χLS n |χm ⟩. At very low temperature, the initial density matrix is written as ρ0 = |Ψ0⟩⟨Ψ0|, and its expression at finite temperature is given in section 4 of the Supporting Information. We have solved the set of nonlinear differential equations of motion (eq 6) for γ ∼ 1.5 ps−1, using Euler’s method (see section 5 of the Supporting Information for the technical details) for different tunneling interaction, J, in the low- and relatively high-temperature regimes. At each time, we compute the quantum averages of ⟨σ⟩ and ⟨Q⟩ and that of the elastic energy. Figure 2 illustrates the case of the lowtemperature limit, corresponding to kBT ≪ ℏω. Because of the initial density matrix, the system starts at point A, (⟨Q⟩ ∼ QLS eq ≈ − 0.1 Å and ⟨σ⟩ ∼ +1), from which it evolves, in a very short period (less than 100 fs), toward point B, expanding suddenly so as to exceed the equilibrium value of the HS distortion (QHS eq ≈ 0.1 Å), while the electronic state, (i.e., ⟨σ⟩ value), experiences a very weak change of less than 1%

(Em, is the mth eigenvalue of 1

Hamiltonian 1) and Γmj = 2 ∑k (γkm + γkj). The expressions of γmk, the damping factors, (detailed in section 3 of the Supporting Information) depend on the bath spectral function and are given by γmk = γ[akm2N (ωmk) + amk2(N (ωkm) + 1)], in which γ is the relaxation (damping) rate part and 1 N(ω) = ℏω/kBT . e −1 A crucial issue of this work is the assumption that the optical pulse duration (some femtoseconds) is short compared to the time scales of vibrations (picoseconds), which supports the vertical transition picture of Figure 1a. Moreover, the Franck− Condon transition is considered here between the electronic LS and HS states of the SC molecule, while in real cases it occurs between the ground state and the metal-to-ligand chargetransfer state (both LS), from which it decays very quickly to 724

DOI: 10.1021/acs.jpclett.6b00014 J. Phys. Chem. Lett. 2016, 7, 722−727

Letter

The Journal of Physical Chemistry Letters

Figure 3. Time dependence of the internal pressure, P (in GPa), for T = 3 K for two values of the tunneling constant: J = 2 meV (left panel) and J = 20 meV (right panel).

Figure 4. Blue curves: instantaneous elastic energy versus the instantaneous average displacement ⟨Q⟩ for T = 3 K showing the relaxation pathway in the adiabatic potential energy diagrams (black curves) for different values of the tunneling term, J (from left to right): 2, 20, and 50 meV.

(respectively 25%) for J = 2 meV (respectively J = 20 meV). The case J = 50 meV is presented in Figure S2. This first event is shorter for stronger tunneling constant J because it involves the strength of the coupling between the two electronic states. Next, between B and C, coherent vibrational motions take place because of internal energy transfer between the electronic and the vibrational degrees of freedom so as to reach the point C, located in the intersystem crossing region (see Figure 3). In this regime, whose duration is around 1 ps, the molecular vibrations dampen as a result of the coupling with the heat bath. This coupling, which was silent in the first regime, is activated by the interference between the vibrations and thermal exchange rate time windows. From point C, both “order parameters” ⟨Q⟩ and ⟨σz⟩ decay on time following stretched exponentials toward their thermal equilibrium values (blue horizontal line). The phase portraits of Figure 2, representing the system’s behavior in its phase space, also reveal the above-mentioned three dynamical regimes: (i) a very rapid molecular transition to the HS state for ⟨Q⟩ between A to B followed by (ii) a nonlinear dependence of ⟨σz⟩ versus ⟨Q⟩ between B and C. Here, ⟨Q⟩ oscillates with damping, while ⟨σz⟩ may oscillate or not, depending on the strength of the tunneling term. Finally from C to D, a linear dependence is obtained between ⟨σz⟩(t) and ⟨Q⟩(t), interpreted as being due to the occurrence of the BO regime where both quantities track

each other. This overall multiscale relaxation process can be approached through the time dependence of the initial internal stress (internal pressure) generated by the Franck−Condon process on the bond length, arising from the frustration caused by the incompatibility between the initial electronic state (HS, ⟨σz⟩ ∼ 1) and the vibrational state (⟨Q⟩ = QLS eq ). For a weak tunneling, the instantaneous stress is estimated as,

(

S

)

F = −k ⟨Q ⟩ − k ⟨σ z⟩

and vanishes at equilibrium when

S ⟨σ z⟩ k

⟨Q ⟩ = (see Figure 1a). Assuming the SC molecule is a hard sphere, its radius changes as R(t) = R0 + ⟨Q⟩(t), where R0 (= 2 Å) is the average radius between the LS and the HS states. F The internal pressure, evaluated as P(t ) = 4πR2 (given in

Figure 3 for J = 2 and 20 meV), reaches nominal values of 6−8 GPa in the first regime. This value is in excellent agreement with experimental findings of Brillouin scattering, atomic force microscopy, and spectroscopic ellipsometry studies.33−35 Interestingly, when the system enters the BO regime in point C, the stress vanishes (i.e., P(t) ∼ 0), whatever the value of the tunneling coupling, J, as displayed in Figure 3. To better illustrate the intersystem crossing and its intimate relation with the strength of the tunneling coupling, we computed the quantum average of the total potential energy, ⟨Epot⟩(t) during the low-temperature relaxation process, drawn 725

DOI: 10.1021/acs.jpclett.6b00014 J. Phys. Chem. Lett. 2016, 7, 722−727

The Journal of Physical Chemistry Letters



as a phase portrait, ⟨Epot⟩(t) versus ⟨Q⟩(t), on the energy potential diagrams (compare Figure 4). The corresponding high-temperature case is presented in Figures S3 and S4. We see clearly from these figures that the previously admitted36−38 idea that the relaxation path goes through the bottom of the HS state as a long-lived metastable state is true only for weak tunneling coupling and low temperature, where the BO regime takes place in a short time scale. Indeed, for J = 2 meV, the relaxation process in the Epot−Q portrait shows coherent nuclear motions around a fixed average value of ⟨Q⟩ = QHS eq . Although it is not shown here, the system ends its excursion in the LS state after very long time simulations. For J = 20 and 50 meV, the onset of the BO regime is longer because the molecular deformation and the spin state are hardly separated (see Figure S2) and the BO description fails. As a result, the system skips the pure HS state during its evolution. In summary, we proposed a vibronic microscopic description for the nonadiabatic relaxation of the photoexcited spincrossover molecule in contact with a heat bath that we solved exactly in the full quantum treatment, beyond the usually admitted BO approximation. We performed this study as a function of temperature and tunneling mixing, which plays a key role in the validity of the adiabatic description. We found in the low-temperature region the existence of three relaxation regimes during the intersystem crossing process, which are governed by the strength of the tunneling coupling. Interesting extensions of the present model relate with the existence of two active deformation coordinates in the course of the relaxation process,14 whose interplay will open new relaxation channels in a higher-dimension phase space and the consideration of a high-energy MLCT state in the PES process. The case of two or more interacting molecules is also a challenging problem in view of the study of photoinduced phase transitions and collective domino effect. These extensions can be carried out using new promising methods enabling the study of nonadiabatic effects at finite temperature based on quantum Monte Carlo calculations at the variational level.22,23 These stochastic techniques are indeed more appropriate to the large Hilbert space arising from the study of cooperative systems.



REFERENCES

(1) Kimel, A. V.; Kirilyuk, A.; Usachev, P. A.; Pisarev, R. V.; Balbashov, A. M.; Rasing, Th. Ultrafast Non-Thermal Control of Magnetization by Instantaneous Photomagnetic Pulses. Nature 2005, 435, 655−657. (2) Stanciu, C. D.; Hansteen, F.; Kimel, A. V.; Kirilyuk, A.; Tsukamoto, A.; Itoh, A.; Rasing, Th. All-Optical Magnetic Recording with Circularly Polarized Light. Phys. Rev. Lett. 2007, 99, 047601. See also Vahaplar, K.; Kalashnikova, A. M.; Kimel, A. V.; Hinzke, D.; Nowak, U.; Chantrell, R.; Tsukamoto, A.; Itoh, A.; Kirilyuk, A.; Rasing, Th. Ultrafast Path for Optical Magnetization Reversal via a Strongly Nonequilibrium State. Phys. Rev. Lett. 2009, 103, 117201. (3) Nasu, K. Photoinduced Phase Transitions; World Scientific: Singapore, 2004. (4) Yonemitsu, K.; Nasu, K. Theory of Photoinduced Phase Transitions in Itinerant Electron Systems. Phys. Rep. 2008, 465, 1−60. (5) Ishida, K.; Nasu, K. Nonlinearity in the Dynamics of Photoinduced Nucleation Process. Phys. Rev. Lett. 2008, 100, 116403. (6) Koshihara, S.; Tokura, Y.; Takeda, K.; Koda, T. Reversible Photoinduced Phase Transitions in Single Crystals of Polydiacetylenes. Phys. Rev. Lett. 1992, 68, 1148. Koshihara, S.; Tokura, Y.; Takeda, K.; Koda, T. Dynamics of a Photoinduced Phase Transition in Polydiacetylene Crystals. Phys. Rev. B: Condens. Matter Mater. Phys. 1995, 52, 6265−6272. (7) Cailleau, H.; Luty, T.; Koshihara, S.; Servol, M.; Lorenc, M.; Buron-Le Cointe, M.; Collet, E. PIPT from the Beginning to Future. Acta Phys. Pol., A 2012, 121, 297−306. (8) Luty, T.; Cailleau, H.; Koshihara, S.; Collet, E.; Takesada, M.; Lemée-Cailleau, M. H.; Buron-Le Cointe, M.; Nagaosa, N.; Tokura, Y.; Zienkiewicz, E.; Ouladdiaf, B. Static and Dynamic Order of Cooperative Multi-Electron Transfer. Europhys. Lett. 2002, 59, 619− 625. (9) Halcrow, M. A. Spin-Crossover Materials: Properties and Applications; Wiley: Oxford, U.K., 2013. (10) Varret, F.; Boukheddaden, K.; Chong, C.; Goujon, A.; Gillon, B.; Jeftic, J.; Hauser, A. Light-Induced Phase Separation in the [Fe(ptz)6] (BF4)2 Spin-Crossover Single Crystal. Eur. Phys. Lett. 2007, 77, 30007. (11) Chong, C.; Varret, F.; Boukheddaden, K. Evolution of SelfOrganized Spin Domains Under Light in Single-Crystalline [Fe(ptz)6] (BF4)2. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 81, 014104. (12) Desaix, A.; Roubeau, O.; Jeftic, J.; Haasnoot, J. G.; Boukheddaden, K.; Codjovi, E.; Linarès, J.; Noguès, M.; Varret, F. Light-Induced Bistability in Spin Transition Solids Leading to Thermal and Optical Hysteresis. Eur. Phys. J. B 1998, 6, 183−193. (13) Marino, A.; Chakraborty, P.; Servol, M.; Lorenc, M.; Collet, E.; Hauser, A. The Role of Ligand-Field States in the Ultrafast Photophysical Cycle of the Prototypical Iron(II) Spin-Crossover Compound [Fe(ptz)6] (BF4)2. Angew. Chem., Int. Ed. 2014, 53, 3863− 3867. (14) Cammarata, M.; Bertoni, R.; Lorenc, M.; Cailleau, H.; Di Matteo, S.; Mauriac, S.; Matar, F.; Lemke, H.; Chollet, M.; Ravy, S.; et al. Sequential Activation of Molecular Breathing and Bending during Spin-Crossover Photoswitching Revealed by Femtosecond Optical and X-Ray Absorption Spectroscopy. Phys. Rev. Lett. 2014, 113, 227402. (15) Bertoni, R.; Lorenc, M.; Tissot, A.; Boillot, M.; Collet, E. Femtosecond Photoswitching Dynamics and Microsecond Thermal Conversion Driven by Laser Heating in FeIII Spin-Crossover Solids. Coord. Chem. Rev. 2015, 282−283, 66−76. (16) Bertoni, R.; Lorenc, M.; Tissot, A.; Servol, M.; Boillot, M.-L.; Collet, E. Femtosecond Spin-State Photoswitching of Molecular Nanocrystals Evidenced by Optical Spectroscopy. Angew. Chem., Int. Ed. 2012, 51, 7485−7489. (17) Guérin, L.; Hébert, J.; Buron-Le Cointe, M.; Adachi, S.; Koshihara, S.; Cailleau, H.; Collet, E. Capturing One-Dimensional Precursors of a Photoinduced Transformation in a Material. Phys. Rev. Lett. 2010, 105, 246101.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.6b00014. Additional computational methods and data (PDF)



Letter

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by Université de Versailles SaintQuentin, Université Paris-Saclay, Centre National de la Recherche Scientifique and the “Agence Nationale de la Recherche” (ANR Project BISTA-MAT: ANR-12-BS07-003001). Financial support was provided by CNRS, Versailles University, and ANR-BISTAMAT program (ANR-12-BS070030). N.K. thanks Versailles University and CNRS for financial support as an invited researcher. 726

DOI: 10.1021/acs.jpclett.6b00014 J. Phys. Chem. Lett. 2016, 7, 722−727

Letter

The Journal of Physical Chemistry Letters (18) Chang, J.; Fedro, A. J.; van Veenendaal, M. Ultrafast Cascading Theory of Intersystem Crossings in Transition-Metal Complexes. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 82, 075124. (19) van Veenendaal, M.; Chang, J.; Fedro, A. J. Model of Ultrafast Intersystem Crossing in Photoexcited Transition-Metal Organic Compounds. Phys. Rev. Lett. 2010, 104, 067401. (20) Piepho, S. B.; Krausz, E. R.; Schatz, P. N. Vibronic Coupling Model for Calculation of Mixed Valence Absorption Profiles. J. Am. Chem. Soc. 1978, 100, 2996−3005. (21) Levine, B. G.; Martinez, T. J. Isomerization Through Conical Intersections. Annu. Rev. Phys. Chem. 2007, 58, 613−634. (22) Mazzola, G.; Sorella, S.; Zen, A. Finite-temperature Electronic Simulations Without the Born-Oppenheimer Constraint. J. Chem. Phys. 2012, 137, 134112. (23) Tubman, N. M.; Kylänpäa,̈ I.; Hammes-Schiffer, S.; Ceperley, D. M. Beyond the Born-Oppenheimer Approximation with Quantum Monte Carlo Methods. Phys. Rev. A: At., Mol., Opt. Phys. 2014, 90, 042507. (24) Boukheddaden, K.; Nishino, M.; Miyashita, S. Molecular Dynamics and Transfer Integral Investigations of an Elastic Anharmonic Model for Phonon-Induced Spin Crossover. Phys. Rev. Lett. 2008, 100, 177206. (25) Boukheddaden, K. Anharmonic Model for Phonon-Induced First-Order Transition in 1-D Spin-Crossover Solids. Prog. Theor. Phys. 2004, 112, 205−217. (26) D’Avino, G.; Painelli, A.; Boukheddaden, K. Vibronic Model for Spin Crossover Complexes. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, 104119. (27) Boukheddaden, K.; Shteto, I.; Hôo, B.; Varret, F. Dynamical Model for Spin-crossover Solids. I. Relaxation Effects in the Mean-field Approach. Phys. Rev. B: Condens. Matter Mater. Phys. 2000, 62, 14796; Dynamical Model for Spin-crossover Solids. II. Static and Dynamic Effects of Light in the Mean-field Approach. Phys. Rev. B: Condens. Matter Mater. Phys. 2000, 62, 14806. (28) Wajnflasz, J.; Pick, R. Transitions ≪Low Spin≫-≪High Spin≫ Dans Les Complexes De Fe2+. J. Phys. Colloq. France 1971, 32, C1-91− C1-92. (29) Carmichael, H. An Open Systems Approach to Quantum Optics; Springer-Verlag: Berlin, Germany, 1991. (30) Suárez, A.; Silbey, R. Hydrogen Tunneling in Condensed Media. J. Chem. Phys. 1991, 94, 4809−4816. (31) Agarwal, G. S. Quantum Optics; Cambridge University Press: Cambridge, U.K., 2012. (32) Gawelda, W.; Pham, V.-T.; Benfatto, M.; Zaushitsyn, Y.; Kaiser, M.; Grolimund, D.; Johnson, S. L.; Abela, R.; Hauser, A.; Bressler, C.; et al. Structural Determination of a Short-Lived Excited Iron(II) Complex by Picosecond X-Ray Absorption Spectroscopy. Phys. Rev. Lett. 2007, 98, 057401. (33) Jung, J.; Bruchhauser, F.; Feile, R.; Spiering, H.; Gütlich, P. The Cooperative Spin Transition in [FexZn1‑x(ptz)6](BF4)2: I. Elastic Properties − an Oriented Sample Rotation Study by Brillouin Spectroscopy. Z. Phys. B: Condens. Matter 1996, 100, 517−522. (34) Félix, G.; Mikolasek, M.; Peng, H.; Nicolazzi, W.; Molnár, G.; Chumakov, A. I.; Salmon, L.; Bousseksou, A. Lattice Dynamics in Spin-crossover Nanoparticles Through Nuclear Inelastic Scattering. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 91, 024422. (35) Boukheddaden, K.; Loutete-Dangui, E. D.; Codjovi, E.; Castro, M.; Rodriguéz-Velamazán, J. A.; Ohkoshi, S.; Tokoro, H.; Koubaa, M.; Abid, Y.; Varret, F. Experimental Access to Elastic and Thermodynamic Properties of RbMnFe(CN)6. J. Appl. Phys. 2011, 109, 013520. (36) Cannizzo, A.; Milne, C. J.; Consani, C.; Gawelda, W.; Bressler, Ch.; van Mourik, F.; Chergui, M. Light-induced Spin Crossover in Fe(II)-based Complexes: the Full Photocycle Unraveled by Ultrafast Optical and X-ray Spectroscopies. Coord. Chem. Rev. 2010, 254, 2677− 2686. (37) Gawelda, W.; Cannizzo, A.; Pham, V.-T.; El Nahhas, A.; Milne, C. J.; van der Veen, R.; Bressler, C.; Chergui, M. Light-Induced Spin Crossover Probed by Ultrafast Optical and X-ray Spectroscopies. Chimia 2007, 61 (4), 179−183.

(38) Ohkoshi, S.; Imoto, K.; Tsunobuchi, Y.; Takano, S.; Tokoro, H. Light-Induced Spin-crossover Magnet. Nat. Chem. 2011, 3, 564−569.

727

DOI: 10.1021/acs.jpclett.6b00014 J. Phys. Chem. Lett. 2016, 7, 722−727