Nature Communications

Nat Commun. 9: 4502

Ultrafast carbon monoxide photolysis and heme spin-crossover in myoglobin via nonadiabatic quantum dynamics

http://orcid.org/0000-0003-1897-1916Konstantin Falahati¹, Hiroyuki Tamura², http://orcid.org/0000-0002-9727-9049Irene Burghardt

burghardt@chemie.uni-frankfurt.de

¹, http://orcid.org/0000-0002-2131-7328Miquel Huix-Rotllant

miquel.huixrotllant@univ-amu.fr

¹³

1. 0000 0004 1936 9721grid.7839.5Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Frankfurt, 60438 Germany

2. 0000 0001 2151 536Xgrid.26999.3dDepartment of Chemical System Engineering, The University of Tokyo, Tokyo, 113-8656 Japan

3. 0000 0004 4902 8637grid.462456.7Aix Marseille Univ, CNRS, ICR, Marseille, France

© The Author(s) 2018

DOI: 10.1038/s41467-018-06615-1

Published online: 29 October 2018

Abstract

Light absorption of myoglobin triggers diatomic ligand photolysis and a spin crossover transition of iron(II) that initiate protein conformational change. The photolysis and spin crossover reactions happen concurrently on a femtosecond timescale. The microscopic origin of these reactions remains controversial. Here, we apply quantum wavepacket dynamics to elucidate the ultrafast photochemical mechanism for a heme–carbon monoxide (heme–CO) complex. We observe coherent oscillations of the Fe–CO bond distance with a period of 42 fs and an amplitude of ∼1 Å. These nuclear motions induce pronounced geometric reorganization, which makes the CO dissociation irreversible. The reaction is initially dominated by symmetry breaking vibrations inducing an electron transfer from porphyrin to iron. Subsequently, the wavepacket relaxes to the triplet manifold in ∼75 fs and to the quintet manifold in ∼430 fs. Our results highlight the central role of nuclear vibrations at the origin of the ultrafast photodynamics of organometallic complexes.

Abstract

Myoglobin bound to carbon monoxide undergoes an ultrafast light-induced reaction, which ends up in a photolyzed carbon monoxide and a spin transition of the iron center. Here, the authors employ quantum wavepacket dynamics to show that photolysis precedes the spin transition, a mechanism dominated by strong electron-nuclear couplings.

Introduction

Hemeproteins containing a porphyrin-iron complex (heme) play a major role in the storage and transport of diatomic molecules¹. The diatomic ligand dissociation occurs concurrently to a spin-crossover (SCO) transition from low-spin (LS) to high-spin (HS) of the FeII center². The heme complex with carbon monoxide (CO) in myoglobin is one of the most studied hemeproteins of this kind³–¹⁸. Heme–CO is initially in a singlet LS state (S = 0), which transforms to a quintet HS state (S = 2) upon CO dissociation¹⁰. The unbound CO in the myoglobin cavity along with the motions of the remaining heme initiate a “protein quake” that opens a channel for the diatomic molecule to escape³. Over a century ago, Halden and Lorrain discovered that the CO dissociation can be initiated photochemically¹⁹, by a mechanism which is still under debate. The photolysis can be initiated upon absorption to the lowest singlet porphyrin band (¹Q) or the second porphyrin band (¹B)²⁰. Here, we focus on the photochemistry of the heme–CO complex upon excitation to the lowest Q-band of porphyrin.

The heme-CO photolysis is an ultrafast process. Recent pump-probe X-ray experiments of myoglobin, with an initial pump pulse exciting the ¹Q state, seem to agree on a two-step kinetic reaction: (i) a first step taking <50–70 fs, attributed to both CO photolysis and partial SCO, and (ii) a completion of the spin transition to the HS state in ~300–400 fs⁵,⁶,²⁰. Despite the numerous studies on heme-CO photolysis, the kinetics and mechanism of dissociation are still under debate, notably regarding the ultrafast nature of the reaction, and the spin and character of the photolytic state. The rate of photolysis has never been experimentally reported, apart from the upper bound of 50–70 fs. As far as the photolytic state is concerned, the most widely accepted hypothesis is that dissociation occurs from a metal-to-ligand charge-transfer (MLCT) state¹⁵–¹⁸,²⁰. Still, experiments and theory do not provide a unified picture to date. On the one hand, Franzen, Martin et al., based on time-resolved absorption and Raman experiments, described a “rapid spin state change that must precede photolysis”²⁰. They consider that photolysis is occurring from a triplet metal→porphyrin ring transfer (³MLCT dissociation)²⁰,²¹. In the model by Franzen et al., the ultrafast reaction is due to valence tautomerism. This mechanism implies a rapid interconversion of several quasi-degenerate electronic states involving d-transitions, some of which are dissociative for the CO bond. On the other hand, Head-Gordon et al., employing time-dependent density-functional theory, proposed a dissociation in a singlet metal → σCO* state (¹MLCT dissociation)¹⁵–¹⁸. In the model of Head-Gordon et al., CO photolysis occurs in a Marcus-like process, in which the initial ¹Q population is transferred to the photolytic state after crossing a barrier of <0.2 eV. However, this hypothesis seems in contradiction with experimental observations of absence of fluorescence emission, an ultrafast reaction and a unity quantum yield for the heme-CO complex²². The theoretical method used provided an incomplete description of ultrafast heme-CO photolysis, since the coupled electron-nuclear motion was neglected. Recent theoretical studies of similar iron complexes indicate the fundamental role of nuclear motions in explaining the ultrafast nature of intersystem crossing (ISC)²³–²⁵. Such electron-nuclear strong coupling has been experimentally observed in time-resolved X-ray absorption in a Fe(II) tris(2,2′-bipyridine) complex²⁶,²⁷ and also in the ferricyanide ion²⁸.

Here, we resolve the mechanism of ultrafast heme–CO photolysis by means of quantum wavepacket dynamics, which accounts for nuclear and electronic coherent motion²⁹. We find that CO photolysis occurs in around 20 fs in the ¹MLCT band, prior to the spin transition. The ultrafast nature of the photolysis is due to strong vibronic couplings and a band of quasi-continuous states. Upon dissociation, the SCO process induces a sequential transition of the remaining heme first into the triplet ³MLCT and second to the HS quintet manifold ⁵MLCT.

Results

Heme–CO model and electronic structure

The model taken as representative of the active center of the myoglobin protein consists of an iron(II) encapsulated in a porphyin ring, and axially ligated to a CO and an imidazole. Hereafter this model is referred to as heme–CO for the complex with carbon monoxide and simply heme for the remaining ligands. The LS state of the heme–CO complex has an electronic configuration (a₂)²(a₁)² (ex*,ey*)0 for the porphyrin and (dxy)²(dxz, dyz)⁴(dz2)0 (dx2-y2)0 for Fe. In its minimum energy structure, the iron and porphyrin are coplanar and the CO is bound to Fe in an upright position perpendicular to the plane. The Fe–C(O) bond is 1.80 Å, around 0.26 Å shorter than the Fe–N bonds with porphyrin or imidazole, which are 2.03 Å and 2.06 Å, respectively. Upon photolysis, the Fe changes into a quintet state, with electronic configuration (dxy)¹(dxz, dyz)³dz21 (dx2-y2)1. In the HS, the Fe–CO bond is unstable. The iron is displaced from the porphyrin plane, distorted to a square-pyramidal structure. The Fe–N bond length is 2.10 Å for the porphyrin nitrogens and 2.17 Å for imidazole.

The full electronic spectrum of heme–CO in the Franck-Condon region is shown in Table 1. This spectrum has been computed in the gas phase, at the CASSCF(10,9)/CASPT2/ANO-RCC-VDZP level of theory for the minimum energy structure of the S₀ ground state (the active space is shown in Supplementary Fig. 1). More details can be found on the effect on the spectrum of the basis set (Supplementary Table 1), the level of theory (Supplementary Table 2) and the model (Supplementary Tables 3–5). At this geometry, the HS state is found 0.7 eV above the LS state. The bright states are the ¹Q states, which correspond to a π → π* transition localized on the porphyrin ring. The ¹Qx and ¹Qy states have an oscillator strength of 1.37⋅10⁻² and 9.13 ⋅ 10⁻³ au for Qx and Qy, respectively. In our model, these transitions are found at 2.73–2.79 eV, while experimentally these bands are found at 2.14–2.30 eV³⁰. These difference are not unexpected, since on the one hand we do not take into account the electrostatic environment of the protein, and on the other hand the statistical average of several protein conformations. Quantum Mechanics/Molecular Mechanics (QM/MM) calculations decrease the overall spectrum by around 0.84 eV for the lowest part of the spectrum. More importantly, the energetic order and the gap between the states in our model as compared with QM/MM is essentially conserved (Supplementary Table 6). This is essential to guarantee a realistic dynamical treatment.

See full tableTable 1. Electronic spectrum of the heme–CO model in the gas phase

The lowest valence states (up to ca. 2 eV) correspond to dark doubly-degenerate metal-centered (MC) states, involving transitions from the occupied dxy, dxz and dyz orbitals to the unoccupied dz2 and dx2-y2 orbitals. The crystal field splitting of the iron orbitals corresponds to a square pyramidal complex. The MLCT band is formed by states of mainly MLCT character, although some MC and LMCT character is also observed (Supplementary Table 7 for a detailed analysis of the state character). These bands of states have strong multi-configurational character, and were not correctly represented in previous single-reference studies based on density-functional theory¹⁵,¹⁸.

Quantum dynamics

Quantum wavepacket dynamics has been performed using a vibronic model containing the main vibrational coordinates (heme doming, symmetry-breaking, rotational and dissociative modes, see Supplementary Fig. 2 for a plot of the vibrations and Methods section for further details.) The evolution of diabatic state populations during the first 500 fs is shown in Fig. 1a. The results have been obtained by sampling 10 independent quantum dynamics simulations with different initial conditions. The initial conditions are obtained by projecting the heme–CO geometries extracted from molecular dynamics snapshots of a myoglobin protein onto our vibronic model Hamiltonian (Supplementary Table 8.) In the Figure, only the states which are mainly populated are shown (Supplementary Fig. 3 for the dynamics of all states.) Following the dynamics of each spin manifold, a sequential transfer is found to occur, first from singlet to triplet (S → T) and then from triplet to quintet (T → Q). The time constant for S → T transfer is estimated as 76 ± 15 fs, while the T → Q decay time is 429 ± 70 fs, in perfect agreement with the experimental rates of Cammarata et al. and Franzen et al.⁵,⁶,²⁰ (for the timescales for each initial condition, see Supplementary Table 9 and Supplementary Note 1.) In the singlet manifold, we observe an abrupt decay of the ¹Q population, which is mainly transferred to the singlet ¹MLCT manifold. The rate for this transfer is estimated to 26 ± 7 fs. The overall reaction is thus sequential ¹Q → ¹MLCT → ³MLCT → ⁵MLCT (see Fig. 1b). The initial step is a complete transfer from the ¹Q state to the ¹MLCT state, resulting in a negligible ¹Q population after ~100 fs. During the remaining transfer steps the three spin manifolds coexist. This is because the three manifolds are close in energy and strongly mixed through spin-orbit coupling. However, a clearly dominant spin state prevails during different intervals: singlet for t < 75 fs, triplet for 75 < t < 425 fs and quintet at later times.

Fig. 1View larger versionFig. 1. Photodynamics of photolysis and spin-crossover. Quantum photodynamics of heme–CO complex during the first 0.5 ps, with initial conditions averaged over 10 molecular dynamics snapshots. a Evolution of diabatic populations for states ¹Q (magenta), ¹MLCT (orange), triplet band (green) and quintet band (blue). The ¹Q population rapidly decays giving rise to ¹MLCT population dominating by 75 fs, at which point the triplet population increases. The quintet population builds up more slowly, and evolves into the dominant state at around 350 fs. b Schematic representation of the reaction mechanism and interpretation in terms of time constants. Upon initial excitation to the Q-band, the metal-to-ligand charge transfer (MLCT) state is populated in ∼25 fs. In a second step, the system relaxes to the triplet (∼75 fs) and to the lowest quintet state (∼430 fs). Black arrows indicate the direction of the electron transfer and the main nuclear motions. c Evolution of the Fe–C(O) distance (magenta, left axis) and the Fe out-of-plane distance (black, right axis). Large amplitude motions are observed with a period of oscillation of 40 fs. The amplitude of oscillation is initially 0.9 Å and converges towards a value of 2.2 Å. At this distance, the CO is essentially photolyzed. The standard deviation of these geometric values is shown as a shaded area. In the inset, the Fourier transform of the Fe–C(O) oscillations is shown (in cm⁻¹)

In Fig. 1c, we show the evolution of the Fe–C(O) and the Fe out-of-plane distances for the first 500 fs (see Supplementary Figs. 4–6 for the corresponding evolution per set of initial conditions). Clearly, we observe damped large amplitude stretching motions for the Fe–C distance, which converges to a quasi-stationary value of 2.2 Å. These motions are coherent, with a period of 42 fs. The Fourier transform of this signal shows a frequency of 800 cm⁻¹, much faster than the Fe–C(O) stretching frequency of 488 cm⁻¹. In 21 fs, the Fe–CO distance oscillates between the equilibrium distance (1.7–1.8 Å) and 2.5 Å. These are the typical distances at which the ground-state dissociation occurs⁷. After the first oscillation, the wavepacket is essentially in the MLCT bands, which are repulsive for CO such that the CO does not recombine and undergoes continued oscillations. Relaxation of the structure damps these oscillations until an equilibrium value of 2.2 Å is reached. Both the equilibrium value and the coherent frequency might be disrupted by the presence of the protein, an effect not included in the present model. This would introduce a faster decoherence of the Fe and CO interactions than what we observe in our model. As for the evolution of the Fe out-of-plane motion, we observe oscillations of ±0.2 Å around the initial position, in good agreement with recent time-resolved X-ray crystallography data⁴. As expected, the iron center reacts to the elongation of the Fe–CO distance by following the CO. Very rapidly, it oscillates back to the porphyrin plane, where part of the energy is dissipated to histidine. In the second and subsequent oscillations, the oscillations of Fe and CO are out of phase, indicating that the bond is photolyzed. A large dispersion of the Fe–C distances, which is appearing after 0.2 ps, is due to the delocalization of the wavepacket over the 1,3,5MLCT bands, indicating an unbound CO  molecule. This dispersion is less marked for the Fe out-of-plane distance due to the more localized nature of the Fe atom, bounded strongly to the heme and the proximal histidine.

Experimental evidence shows that the CO is photolyzed within 70 fs⁵,⁶,²⁰. In order to determine whether dissociation occurs from the singlet or the triplet state, we have considered a reduced model with only the singlet manifold, in which we eliminate the coupling parameters to the triplet and quintet states (Supplementary Figs. 7–9). In this case, an ultrafast transfer from Q → ¹MLCT is still observed. The Fe–C(O) distance again exhibits oscillations with an amplitude of 0.7 Å. The fact that the amplitude is smaller than in the full model can be ascribed to the absence spin-orbit couplings which introduce strong mixing of the ¹MLCT surfaces with ³MLCT states, that extend the region of oscillations. From these results, we can infer that the transfer to the ¹MLCT manifold is sufficient to dissociate the CO, although the presence of ³MLCT favors further the photolysis. The dissociation is thus happening between 0.5 and 1.5 periods of the Fe–C(O) oscillation (20–60 fs), when the wavepacket is mainly in the ¹MLCT state.

After photolysis, the SCO mechanism brings the system sequentially to a HS state. After ~75 fs, the ³MLCT band is clearly dominant, although a residual ³MC participation is also observed (Supplementary Fig. 1). From the triplet manifold, a population transfer to the HS states is slowly building up, until it becomes dominant at around 400 fs. At this time, the wavepacket is distributed across the ⁵MLCT and ⁵MC excited states, staying trapped in these states. This is consistent with the slow time constant of ~3 ps observed in myoglobin, which has been attributed to a vibrational cooling due to the protein²⁰,³¹. Such excited state trapping has been observed recently in similar iron complexes²⁷. This trapping can be understood in terms of the potential energy surfaces depicted in Fig. 2. Upon excitation to the ¹Q band, a sequential ultrafast transfer, first to ¹MLCT and then to ³MLCT occurs. These transfers are energetically favored and therefore extremely fast. In the ³MLCT manifold, the wavepacket relaxes in a distribution of ³MLCT and ³MC states. The ³MC state can undergo direct transfer to the ⁵MLCT state, which is, however, energetically less favorable (and thus slower) than the singlet → triplet MLCT. The transition from the lowest ⁵MLCT to the ⁵MC state is energetically unfavoured. On the other hand, the ³MC states are energetically separated from the rest of the states, such that a transfer to the ⁵MC manifold becomes inefficient.

Fig. 2View larger versionFig. 2. Potential energy surface along a dissociative mode. Plot of the potential energy surfaces along the Fe–C(O) distance. The LS state S₀, and the ¹Q states are shown in black and magenta, respectively. The singlet manifold is shown in orange, the triplet manifold in green and the quintet manifold in blue. The metal-centered (MC) states are shown explicitly, whereas the metal-to-ligand charge transfer (MLCT) states are depicted as a quasi-continuous band. Further, a schematic evolution of the wavepacket dynamics is shown. After absorption to the Q-band, the wavepacket undergoes large-amplitude oscillations in the Fe–CO coordinate on the ¹MLCT state. Then, it cascades down acquiring more ³MLCT and ⁵MLCT character as the spin crossover transitions occur. Finally, the wavepacket disperses, and the Fe–C(O) distance oscillations decrease to a value of 2.2 Å

Vibronic effects

The initial structure is close to a C4v symmetric geometry, in which the CO is in upright position and the porphyrin ring is close to the square planar conformation. The origin of the heme–CO bond is the typical σ-donation π-back-donation mechanism. A partial electron transfer from an occupied σ orbital of CO to an empty dz2 orbital of Fe occurs simultaneously to a partial back transfer of electron density from the occupied dyz and dxz of Fe to the empty πy,CO* and πx,CO* orbital, respectively. At this geometry, the π orbitals of porphyrin do not overlap with the d-orbitals.

Different types of vibrations can weaken the heme–CO bond, thus liberating the diatomic molecule. In Fig. 3, a schematic representation of the most representative vibrations is shown (for a detailed analysis of all vibrations, see Supplementary Table 10). The most important vibration is the dissociative coordinate (Fig. 3a), which corresponds to a symmetric Fe–CO stretching. This is the main reaction coordinate for CO photolysis. The Fe–CO stretching decreases the overlap between the iron and CO orbitals, thus decreasing both the strength of σ-donation and the π-back-donation. Once the CO is released, the remaining complex is then stabilized in the quintet state. This brings the Fe center out-of-plane, stabilized by π-interactions with the porphyrin antibonding orbitals.

Fig. 3View larger versionFig. 3. Vibrational motions and effect on the Fe–CO interaction. Representation of the nuclear motions (left) and the effect on the Fe–CO orbital interaction (right). The σCO and dxz/dyz orbitals (in blue) are occupied, while ex/y, πx∕y,CO* and dz2 orbitals (in red) are unoccupied. Blue arrows indicate the direction of the partial electron transfer, due to the overlap between an occupied and an unoccupied orbital (solid and dashed arrows indicate strong and weak overlap). Magenta arrows indicate the vibrational motion of the specified atoms. a In the initial structure, the Fe–CO bond is strongest due to partial electron transfer dxz∕dyz→πx*∕πy* (π-back donation) and dz2←σCO (σ-donation). b Dissociative vibration, which corresponds to a stretching of the Fe–CO bond. This structure destabilizes the Fe–CO bond and stabilizes by a partial electron transfer dxz∕dyz→ex*∕ey* when the Fe atom is out of plane. c Rotational vibrations, which correspond to the rotation of the CO moiety. This vibration decreases the overlap between Fe and CO orbitals. d Symmetry-breaking vibration, in which the porphyrin is distorted inducing Jahn-Teller and pseudo Jahn-Teller vibronic couplings in the excited state. This vibration has no effect on the Fe–CO bond, but activates the coupling between the porphyrin and the Fe densities

Activation of the dissociative vibration is fundamental for releasing the CO from the heme. This vibration is of a₁ symmetry, as it is a vibration along the principal symmetry axis. However, the initial photon is absorbed exclusively in the porphyrin moiety. The porphyrin optical Q-state is E-symmetric, thus localized in the perpendicular plane with respect to the principal axis. This implies that the dissociative vibration is not activated in the quasi-C4v initial geometry. The dissociation of CO thus requires an initial energy transfer from the porphyrin plane to the Fe center, which can only be possible by populating symmetry breaking vibrations. In the excited states, strong symmetry breaking of Jahn-Teller (JT) and pseudo Jahn-Teller (pJT) type occurs³²,³³. The JT effect is the sudden increase of inter-state couplings between an E-degenerate state along symmetry breaking vibrational coordinates. These JT couplings induce a double-well minimum in the potential energy surfaces at elongated coordinates, while the energy is maximum at the initially symmetric structure. The pJT effect is analogous to the JT effect but applies to non-degenerate states, usually requiring more than a single vibration and more than one electronic state. Both mechanisms are characterized by strong vibronic couplings. In the particular case of the heme–CO complex, the E-symmetry of the Q-band can be broken by b₁ and b₂ vibrations, in the so-called E⊗b JT mechanism³². A representative vibration of this type is depicted in Fig. 3d. Such a vibration is not directly coupled to the Fe–CO bond, but is fundamental to induce the interstate coupling between the Q-band and the MLCT band. The relaxation energy induced by such vibrations can be as large as 0.5 eV, which rapidly leads the wavepacket far from the Franck-Condon region. Large amplitude motions have been experimentally observed in similar iron complexes²⁶.

The rotation of CO (see Fig. 3c) has been discussed at length in the literature, generating a controversy on whether it has an effect on the photochemistry⁸,⁹,¹⁶,³⁴. The rotation of CO slightly weakens the Fe–CO bond due to a less efficient overlap between the densities of each moiety. Furthermore, this rotation induces a Renner-Teller type of coupling (rovibronic coupling), in which e-symmetry bending modes of CO and a₁-type vibrations couple E and A₁ states, in what is known as the (E⊕A₁)⊗(e⊕a₁) effect. This mechanism can also couple the Q-band (E symmetry) and the MLCT band (A₁ symmetry) to the CO bending (e symmetry) and the dissociative iron-CO stretching (a₁ symmetry). The population evolution of the singlet manifold excluding the JT modes— and thus the rotational modes that are responsible of the Q-MLCT coupling—shows a similar population and Fe–C kinetics. This indicates that any symmetry breaking coordinate, related to protein fluctuations, rotations or symmetry breaking Jahn-Teller effects, will accelerate the coupling to the MLCT band.

The porphyrin doming (a₁ vibration of the porphyrin ring, see modes 8,12 in Supplementary Fig. 3) has been pointed out as an important motion for the photodynamics⁵,⁶,³⁵. The doming motion is out-of-plane, and has a slow period of 371 fs. It is of the same symmetry as the dissociative mode, to which it is strongly coupled. The potential energy surface cuts along the doming mode indicate a weak reorganization energy and a small vibronic coupling for all potential energy surfaces. Therefore, even though this motion is necessary to stabilize the high-spin state, it plays a passive role and does not significantly affect the short-time dynamics of photolysis (Supplementary Figs. 10–12). Therefore, we conclude that the doming motion is in fact a consequence of the CO dissociation, that is strongly coupled to the Fe–CO dissociation coordinate, and is activated during the second kinetic step of 280 fs⁵,⁶,²⁰.

Spin crossover

The photochemical spin crossover mechanism is an important reaction in many organometallic complexes, in which ligand chromophores allow a usually fast low-to-high spin transition of the complex. Recent theoretical and experimental studies for similar iron complexes have shown that ligands play a fundamental role in modulating the time constants and the electronic pathway during the photoreaction⁴–⁶. The spin crossover mechanism is mediated by the relativistic spin-orbit coupling (SOC). Although this is a purely electronic transition, the vibrational effects on the SOC seem fundamental to explain the ultrafast nature of spin crossover in metallic complexes. In the particular case of heme–CO, the main SOC strengths are about 200–300 cm⁻¹, which are of similar size as the JT vibronic couplings. Therefore, internal conversion and intersystem crossing are competing mechanisms in the dynamical relaxation mechanism.

In the heme–CO complex, the initially populated Q-band of porphyrin exhibits a negligible SOC with any other state. In general, SOC is strong when metal density is present in the electronic state. Therefore, the wavepacket is initially dominated by the symmetry breaking vibronic couplings in the singlet manifold, leading to the transition ¹Q → ¹MLCT and ¹Q → ¹MC. Both the ¹MLCT and ¹MC have strong spin-orbit coupling strength with ³MLCT and ³MC states, respectively. In general, the SOC is strong when the selection rules ΔS = 1 and ΔL = 1 are obeyed. The latter condition implies that one of the d orbitals changes orientation between the triplet and the singlet state. For example, a strong coupling occurs between the singlet state of character dxy→dx2-y2 with the triplet state of character dxz→dx2-y2. We observe that singlet-triplet transitions occur at a faster rate in the MLCT band than within the MC states (Supplementary Fig. 4). At around 50 fs, the ³MLCT population becomes dominant, at the expense of the singlet manifold populations. This is because there are a larger number of crossing points between MLCT states of different multiplicity than between MC states (see Fig. 2 and Supplementary Fig. 4). This results in a dominant population of ³MLCT and ⁵MLCT states during the spin-crossover transitions. The ³Q band triplets, which were pointed out as playing a role in the photophysics of heme–CO, contribute negligibly to the dynamics²⁰. The ¹MC state exhibits a residual population after 100 fs, following prevalent transfer to the ³MC state, with a population which is two times smaller than the ³MLCT population. Finally, the ³MLCT → ⁵MLCT transfer builds up population in the quintet band. A ³MC → ⁵MC transfer, however, is not observed, due to the fact that the ³MC and ⁵MC states do not cross in an accessible region. Our results show rather that the wavepacket remains trapped in the ⁵MLCT band. The trapping occurs in regions of the energy surface where relaxation to lower states (in this case MC states) is slow due to large gaps and small vibronic couplings, whereas a band of electronic states acts as an energy dissipation force, thus localizing the wavepacket in an excited state. Such observation for heme–CO is in line with the valence tautomerism model of Franzen, Martin et al.²⁰ In summary, the existence of quasi-continuous MLCT bands degenerate to the ¹Q state and strong vibronic couplings are at the origin of the ultrafast photolysis and spin-crossover in the heme–CO complex.

Discussion

The mechanism of heme–CO photolysis and spin crossover has been elucidated by means of high-dimensional quantum dynamical calculations. The results indicate that the vibronic mechanism can be summarized in terms of a sequential ¹Q → ¹MLCT/MC → ³MLCT → ⁵MLCT transfer. During the first step, which is completed within 20–60 fs, the CO is dissociated, prior to the spin-crossover mechanism. The dissociation occurs in the singlet MLCT excited states. For these states, the σ − π bond between Fe and CO bond is weakened. Subsequently, the wavepacket relaxes stepwise to the quintet band under the effect of spin-orbit coupling. We assign the three experimental rates observed for the heme–CO complex in myoglobin, recently reported in Ref. ²⁰ and Refs. ⁵,⁶ as follows: (i) the 50–70 fs time constant corresponds to the singlet-triplet spin crossover, in which the population is mainly in the ³MLCT excited state; (ii) the 300–400 fs time constant corresponds to the triplet-quintet spin crossover, in which the population is mainly in the ⁵MLCT excited state; (iii) the 2.4 ps time constant we tentatively assign to the relaxation of the ⁵MLCT excited state to the final HS state of the heme.

Methods

Electronic structure

Unconstrained structure optimization and normal-mode calculations of S₀ in gas phase have been performed with density-functional theory (DFT) using the B3LYP exchange-correlation functional and the LAN2DZ basis set³⁶–³⁹. These computations have been performed with Gaussian09⁴⁰. Excited states have been computed with CASSCF(10,9)/CASPT2⁴¹,⁴². The active orbitals are the Fe(II) 3d-orbitals and the Gouterman’s porphyrin π-orbitals (see Supplementary Fig. 2 and Supplementary Note 2). State-averaged calculations have been performed for 35 singlet states, 50 triplet states, and 30 quintet states. In all calculations, we employed relativistic Douglas-Kroll-Hess (DKH) Hamiltonian with the ANO-RCC-VDZ basis set⁴³–⁴⁵. Larger basis sets do not significantly affect the energetics (Supplementary Table 2). Multi-state CASPT2 has been performed with frozen core and 116 frozen virtual orbitals. Similarly to coupled-cluster⁴⁶, deleting virtual orbitals accelerates computations and keeps the essential physics (Supplementary Table 2). In particular, the spin crossover transition is well represented and the DFT minima are also CASPT2 minima. Spin-orbit couplings have been obtained by expanding the DKH Hamiltonian in the basis of 20 singlets, 20 triplets and 20 quintets⁴⁵. All multi-configuration results have been obtained with Molcas⁴⁷.

Vibronic Hamiltonian

The model Hamiltonian for representing the dissociation and the spin crossover is of the form:

1. Ĥ(q)= ∑S=02ĤS(q)+ ∑S,S′=12V^SS′SO,[Image: 41467_2018_6615_Article_Equ1.gif]

in which S = 0, 1, 2 is the total spin quantum number, q is the mass-frequency-weighted vibrational coordinate vector, V^SO is the spin-orbit coupling. The vibronic Hamiltonian, represented in the basis of ground-state normal modes, is given by

2. ĤS(q)=T^S(q)+V^S(q),[Image: 41467_2018_6615_Article_Equ2.gif]

containing the harmonic kinetic energy and potential operator. The potential is given by

3. V^S(q)= ∑ms=-∣S∣∣S∣∑I,J=120δIJEIS+VIJS(q)SΨImsSΨJms+c.c.,[Image: 41467_2018_6615_Article_Equ3.gif]

in which ms is the azimuthal spin and SΨIms is the state wavefunction. The EIS energies are diabatic at the S₀ minimum geometry; at the quasi-C4v geometry these can be considered quasi-diabatic. The VIJS(q) potential is derived from one-dimensional potential energy surface cuts which were diabatized by fitting to the following functional form,

4. VIIS(q)=∑i=151i!∂iEIS(q)∂qiq=0qi+αISeβISq-γIS+DIS1-eεISq-qIS2,[Image: 41467_2018_6615_Article_Equ4.gif]

5. VIJS(q)=vIJS(q)q.[Image: 41467_2018_6615_Article_Equ5.gif]

Here, the inter-state couplings VIJS(q) are approximated within a linear vibronic model. The set of parameters are determined by a least square fitting

6. fyS=1N∑IσEImod,SyIS,qσ-EIPT2,S(qσ)2,[Image: 41467_2018_6615_Article_Equ6.gif]

where yS=∂iEIS∂qi,αIS,βIS,γIS,εIS,qIS,vIJS, qσ is the discretized value of a normal mode coordinate and EImod,S(qσ) and EIPT2,S(qσ) are the spin-free adiabatic energies from the model Hamiltonian and the CASPT2 calculation, respectively. The latter are obtained by one-dimensional potential energy surface cuts along each ground-state normal mode direction. The passive normal modes have been kept frozen along the cuts.

The fitting of quasi-continuous bands requires an accurate initial guess. We employ the following protocol: (i) quasi-diabatic surfaces are obtained by maximum overlap of CI coefficients; (ii) pre-fitting is done on the quasi-diabatic surfaces by employing first and second order Taylor expansion or the exponential function, (iii) the adiabatic PT2 surfaces are fitted with the full diagonal Hamiltonian, and (iv) the full Hamiltonian is fitted until the total error is <0.10 eV.

Vibrational modes

Vibrational coordinates have been chosen by performing one-dimensional potential energy cuts and selecting the modes with largest reorganization energy. There are a total of 15 explicit vibrational modes in the model: 1 symmetric and 1 asymmetric Fe-CO stretchings, 2 doming modes, 1 CO stretching, 1 in-plane, and 1 out-of-plane porphyrin vibrations, and 4 doubly degenerate modes (Jahn-Teller modes of porphyrin out-of-plane, a porphyrin in-plane mode and a pure CO bending.) These modes have been selected in order to reproduce the symmetry breaking that allows the Q → MLCT ultrafast transition, the dissociation and the relaxation of the quintet heme. A schematic plot of these vibrations can be found in Supplementary Fig. 3. Duschinsky rotation effects and anharmonic couplings between vibrational modes have been neglected in our model; in particular, the latter are negligibly small compared to the vibronic couplings (see Supplementary Table 11 and Supplementary Note 3).

Quantum dynamics

The Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method⁴⁸,⁴⁹ (Heidelberg package, version 8.5.5) was used to perform wavepacket propagation up to 1 picosecond for a vibronic coupling model comprising 179 electronic states of singlet, triplet and quintet multiplicities, and 15 vibrational modes exhibiting strong vibronic coupling. A four-layer representation of the wavefunction was employed, where 8 single-particle functions (spf) were chosen per layer, with a primitive harmonic-oscillator discrete variable representation (DVR) comprising up to 80 DVR points (see Supplementary Note 4 and Supplementary Fig. 13). Initial conditions were constructed in accordance with the relative oscillator strengths of the singlet ligand and metal-centered states at the Franck-Condon geometry.

Supplementary Material

Open In Web Browser

Footnotes

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Supplementary Information accompanies this paper at 10.1038/s41467-018-06615-1.

Acknowledgements

MHR acknowledges financial support from Agence Nationale de la Recherche (grant ANR-16-CE29-0008, BIOMAGNET) and the Alexander von Humboldt-Foundation. K.F. and I.B. gratefully acknowledge funding by the Deutsche Forschungsgemeinschaft via RTG 1986 “Complex Scenarios of Light Control”. I.B. thanks Prof. Gerhard Hummer (MPI Biophysics, Frankfurt) for making MD snapshots available for the present study. All authors thank Prof. Nicolas Ferré, Dr. Padmabati Mondal, Dr. Anthony Kermagoret and the reviewers for useful comments. This work was granted access to the HPC resources of Aix-Marseille Université financed by the project Equip@Meso (ANR-10-EQPX-29-01) of the program “Investissements d’Avenir” supervised by the Agence Nationale de la Recherche.

Author contributions

M.H.R. designed the research in collaboration with I.B. M.H.R., I.B., and H.T. developed the vibronic Hamiltonian. M.H.R. and K.F. performed the electronic structure calculations. M.H.R., K.F., and I.B. performed the quantum dynamics. M.H.R. wrote the paper. All authors contributed to the discussions and gave comments on the manuscript.

Data availability

The authors declare that the data supplementary the findings of this study are available within the paper and its supplementary information files. Parameters for the model Hamiltonian and potential energy surfaces along each vibrational mode are available from the corresponding authors upon request.

Notes

Competing interests

The authors declare no competing interests.

Articles from Nature Communications are provided here courtesy of Nature Publishing Group

PMC Copyright Notice

The articles available from the PMC site are protected by copyright, even though access is free. Copyright is held by the respective authors or publishers who provide these articles to PMC. Users of PMC are responsible for complying with the terms and conditions defined by the copyright holder.

Users should assume that standard copyright protection applies to articles in PMC, unless an article contains an explicit license statement that gives a user additional reuse or redistribution rights. PMC does not allow automated/bulk downloading of articles that have standard copyright protection.

See the copyright notice on the PMC site, https://www.ncbi.nlm.nih.gov/pmc/about/copyright/, for further details and specific exceptions.

References

Messerschmidt, A., Huber, R., Poulos, T. & Wieghardt, K. (eds.) Handbooks of Metalloproteins, Vol. 1 (John Wiley and Sons, 2001).

Chung, L. W., Li, X. & Morokuma, K. Modeling Enzymatic Reactions in Metalloenzymes and Photobiology by Quantum Mechanics (QM) and Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations. Chap. 3. Quantum Biochemistry (Wiley–VCH Verlag GmbH, 2010).

Ansari A, et al., authors. Protein states and proteinquakes. Proc. Natl Acad. Sci. USA. 1985;82:5000–5004. [CrossRef] [PubMed]

Barends TR, et al., authors. Direct observation of ultrafast collective motions in CO myoglobin upon ligand dissociation. Science. 2015;350:445–450. [CrossRef] [PubMed]

Levantino M, Schiro G, Lemke HT, Cottone G, authors. Ultrafast myoglobin structural dynamics observed with an X–ray free–electron laser. Nat. Commun. 2015;6:6772 [CrossRef] [PubMed]

Levantino M, et al., authors. Observing heme doming in myoglobin with Femtosecond X–ray absorption spectroscopy. Struc. Dynam. 2015;2:041713 [CrossRef]

Harvey JN, author. DFT computation of the intrinsic barrier to CO geminate recombination with heme compounds. J. Am. Chem. Soc. 2000;122:12401–12402. [CrossRef]

Spiro TG, Kozlowski PM, authors. Is the CO adduct of myoglobin bent, and does it matter? Acc. Chem. Res. 2001;34:137–144. [CrossRef] [PubMed]

Rovira C, Schulze B, Eichinger M, Evanseck JD, Parrinello M, authors. Influence of the heme pocket conformation on the structure and vibrations of the Fe–CO bond in myoglobin: a QM/MM density functional study. Biophys. J. 2001;81:435–445. [CrossRef] [PubMed]

Harvey JN, Spin–forbidden CO, authors. Ligand recombination in myoglobin. Faraday Discuss. 2004;127:165–177. [CrossRef] [PubMed]

Meier C, author. Laser control of vibrational excitation in carboxyhemoglobin: a quantum wavepacket study. J. Chem. Phys. 2005;123:044504 [CrossRef] [PubMed]

Debnath A, Falvo C, Meier C, authors. State-selective excitation of the CO stretch in carboxyhemoglobin by mid-IR laser pulse shaping: a theoretical investigation. J. Phys. Chem. A. 2013;117:12884–12888. [CrossRef] [PubMed]

Kühn O, author. Multidimensional vibrational quantum dynamics of CO-heme compounds: ultrafast IVR mediated Fe-CO bond–breaking after CO excitation? Chem. Phys. Lett. 2005;402:48–53. [CrossRef]

Waleh A, Loew GH, authors. Quantum mechanical studies of the photodissociation of carbonylheme complexes. J. Am. Chem. Soc. 1982;104:2346–2351. [CrossRef]

Dreuw A, Dunletz BD, Head-Gordon M, authors. Characterization of the relevant excited states in the photodissociation of CO–ligated hemoglobin and myoglobin. J. Am. Chem. Soc. 2002;124:12070–12071. [CrossRef] [PubMed]

De Angelis F, Car R, Spiro TG, authors. Role of ligand bending in the photodissociation of O₂ vs CO–heme: a time–dependent density functional study. J. Am. Chem. Soc. 2003;125:15710–15711. [CrossRef] [PubMed]

Ohta T, Pal B, Kitagawa T, authors. Excited state property of hardly photodissociable heme–CO adduct studied by time–dependent density functional theory. J. Phys. Chem. B. 2005;109:21110–21117. [CrossRef] [PubMed]

Dunietz BD, Dreuw A, Head-Gordon M, authors. Initial steps of the photodissociation of the CO ligated heme group. J. Phys. Chem. B. 2003;107:5623–5629. [CrossRef]

Halden J, Lorrain SJ, authors. The oxygen tension of arterial blood. J. Physiol. 1896;20:497–520. [CrossRef]

Franzen S, Kiger L, Poyart C, Martin J, authors. Heme photolysis occurs by ultrafast excited state metal–to–ring charge transfer. Biophys. J. 2001;80:2372 [CrossRef] [PubMed]

Petrich JW, Poyart C, Martin JL, authors. Photophysics and reactivity of heme proteins: a femtosecond absorption study of hemoglobin, myoglobin, and protoheme. Biochem. 1988;27:4049–4060. [CrossRef] [PubMed]

Adar F, Gouterman M, Aronowitz S, authors. Fluorescence, resonance Raman, and raditionless decay in several hemoproteins. J. Phys. Chem. 1976;80:2184–2191. [CrossRef]

Sousa C, et al., authors. Ultrafast deactivation mechanism of the excited singlet in the light-induced spin crossover of [Fe(2,2′-bipyridine)₃]²⁺. Chem. Eur. J. 2013;19:17541–17551. [CrossRef] [PubMed]

Eng J, Gourlaouen C, Gindensperger E, Daniel C, authors. Spin-vibronic quantum dynamics for ultrafast excited-state processes. Acc. Chem. Res. 2015;48:809–817. [CrossRef] [PubMed]

Pápai M, Penfold TJ, Moller KB, authors. Effect of tert-butyl functionalization on the photoexcited decay of a Fe(II)-N-heterocyclic carbene complex. J. Phys. Chem. C. 2016;120:17234–17241. [CrossRef]

Auböck G, Chergui M, authors. Sub–50–fs photoinduced spin crossover in [Fe(bpy)₃]²⁺. Nat. Chem. 2015;7:629–633. [CrossRef] [PubMed]

Lemke HT, et al., authors. Coherent structural trapping through wave packet dispersion during photoinduced spin state switching. Nat. Comm. 2017;6:6772

Engel N, et al., authors. Light-induced relaxation dynamics of the ferricyanide ion revisited by ultrafast XUV photoelectron spectroscopy. Phys. Chem. Chem. Phys. 2017;19:14248–14255. [CrossRef] [PubMed]

Manthe U, author. Wavepacket dynamics and the multi-configurational time-dependent Hartree approach. J. Phys. Condens. Matter. 2017;29:253001 [CrossRef] [PubMed]

Bowen WJ, author. The absorption spectra and extinction coefficients of myoglobin. J. Biol. Chem. 1949;179:235–245. [PubMed]

Mizutani Y, Kitagawa T, authors. Direct observation of cooling of heme upon photodissociation of carbonmonoxy myoglobin. Science. 1997;278:443–446. [CrossRef] [PubMed]

Bersuker, I. The Jahn-Teller Effect and Vibronic Interactions in Modern Chemistry. Modern Inorganic Chemistry (Springer Verlag, 1984).

Bersuker I, author. Pseudo-Jahn-Teller effect–a two-state paradigm in formation, deformation, and transformation of molecular systems and solids. Chem. Rev. 2013;113:1351–1390. [CrossRef] [PubMed]

Carme R, Parrinello M, authors. Harmonic and anharmonic dynamics of Fe–CO and Fe–O₂ in heme models. Biophys. J. 2000;78:93–100. [CrossRef] [PubMed]

Klug DD, et al., authors. Doming modes and dynamics of model heme compounds. Proc. Natl Acad. Sci. USA. 2002;99:12526–12530. [CrossRef] [PubMed]

Becke AD, author. Density–functional exchange–energy approximation with correct asymptotic behavior. Phys. Rev. A. 1988;38:3098 [CrossRef]

Lee C, Yang W, Parr RG, authors. Development of the Colle–Salvetti correlation–energy formula into a functional of the electron density. Phys. Rev. B. 1988;37:785 [CrossRef]

Vosko SH, Wilk L, Nusair M, authors. Accurate spin–dependent electron liquid correlation energies for local spin density calculations: a critical analysis. Can. J. Phys. 1980;58:1200–1211. [CrossRef]

Hay PJ, Wadt WR, authors. Ab initio effective core potentials for molecular calculations. potentials for K to Au including the outermost core orbitals. J. Chem. Phys. 1985;82:299 [CrossRef]

Frisch, M. J. et al. Gaussian 09, Revision D.01, Gaussian, Inc., Wallingford CT (2009).

Roos, B. J. The complete active space self–consistent field method and its applications in electronic structure calculations., Vol. 69 In Adv. Chem. Phys. 399 (John Wiley and Sons, Ltd., 1987).

Andersson K, Malmqvist P, Roos BO, authors. Second–order perturbation theory with a complete active space self–consistent field reference function. J. Chem. Phys. 1992;96:1218–1226. [CrossRef]

Roos BO, Lindh R, Malmqvist PA, Veryazov V, Widmark PO, authors. Main group atoms and dimers studied with a new relativistic ANO basis set. J. Phys. Chem. A. 2005;108:2851 [CrossRef]

Roos BO, Lindh R, Malmqvist PA, Veryazvov V, Widmark PO, authors. New relativistic ANO basis sets for transition metal atoms. J. Phys. Chem. A. 2005;109:6575–6579. [CrossRef] [PubMed]

Roos BO, Malmqvist PA, authors. Relativistic quantum chemistry: the multiconfigurational approach. Phys. Chem. Chem. Phys. 2004;6:2919–2927. [CrossRef]

Send R, Kaila VRI, Sundholm D, authors. Reduction of the virtual space for coupled–cluster excitation energies of large molecules and embedded systems. J. Chem. Phys. 2011;134:214114 [CrossRef] [PubMed]

Aquilante F, et al., authors. MOLCAS 8: new capabilities for multiconfigurational quantum chemical calculations across the periodic table. J. Comput. Chem. 2016;37:506–541. [CrossRef] [PubMed]

Meyer HD, Manthe U, Cederbaum LS, authors. The multi–configurational time–dependent Hartree approach. Chem. Phys. Lett. 1990;165:73–78. [CrossRef]

Wang H, Thoss M, authors. Multilayer formulation of the multiconfiguration time–dependent Hartree theory. J. Chem. Phys. 2003;119:1289 [CrossRef]

[Back]

Table 1.

Electronic spectrum of the heme–CO model in the gas phase

State | Energy | Character | Sym.
1⁵MC | 0.70 | dyz,dxy→dx2-y2,dz2 | E
2⁵MC | 0.74 | dxz,dxy→dx2-y2,dz2 | E
1³MC | 0.89 | dxy→dx2-y2 | B₂
3⁵MC | 1.40 | dxz,dyz→dx2-y2,dz2 |
2³MC | 1.43 | dxz→dz2 | E
3³MC | 1.51 | dyz→dz2 | E
4³MC | 1.62 | dxy→dz2 | A₂
5³MC | 1.86 | dyz→dx2-y2 | E
6³MC | 1.90 | dxz→dx2-y2 | E
1¹MC | 1.90 | dxy→dx2-y2 | B₂
1³MLCT | 2.15 | a1,dyz,dxy→e,dx2-y2,dz2 |
2¹MC | 2.35 | dyz→dz2 | E
3¹MC | 2.36 | dxy→dz2 | E
1³Q₁ | 2.43 | a₂ → e | E
⁵MLCT | 2.48 | a2,dyz,dxy→e,dx2-y2,dz2 |
2³Q₂ | 2.49 | a₂ → e | E
1³MLCT | 2.56 | a1,dyz,dxy→e,dx2-y2,dz2 |
¹Qx | 2.73 | a₂ → e/a₁ → e | E
¹Qy | 2.79 | a₂ → e/a₁ → e | E
¹MLCT | 2.83 | a2,dyz→e,dz2 |

Electronic states up to the bright ¹Qx,y state are shown. The excitation energies (in eV) are computed at the CASSCF(10,9)/CASPT2/ANO-RCC-VDZP level of theory at the S₀ minimum energy structure obtained at the B3LYP/LANL2DZ level of theory. In addition, the state nature, the character of the dominant transitions and the symmetry label is shown

[Back]

[Back]

Fig. 1.

Photodynamics of photolysis and spin-crossover. Quantum photodynamics of heme–CO complex during the first 0.5 ps, with initial conditions averaged over 10 molecular dynamics snapshots. a Evolution of diabatic populations for states ¹Q (magenta), ¹MLCT (orange), triplet band (green) and quintet band (blue). The ¹Q population rapidly decays giving rise to ¹MLCT population dominating by 75 fs, at which point the triplet population increases. The quintet population builds up more slowly, and evolves into the dominant state at around 350 fs. b Schematic representation of the reaction mechanism and interpretation in terms of time constants. Upon initial excitation to the Q-band, the metal-to-ligand charge transfer (MLCT) state is populated in ∼25 fs. In a second step, the system relaxes to the triplet (∼75 fs) and to the lowest quintet state (∼430 fs). Black arrows indicate the direction of the electron transfer and the main nuclear motions. c Evolution of the Fe–C(O) distance (magenta, left axis) and the Fe out-of-plane distance (black, right axis). Large amplitude motions are observed with a period of oscillation of 40 fs. The amplitude of oscillation is initially 0.9 Å and converges towards a value of 2.2 Å. At this distance, the CO is essentially photolyzed. The standard deviation of these geometric values is shown as a shaded area. In the inset, the Fourier transform of the Fe–C(O) oscillations is shown (in cm⁻¹)

[Back]

[Back]

Fig. 2.

Potential energy surface along a dissociative mode. Plot of the potential energy surfaces along the Fe–C(O) distance. The LS state S₀, and the ¹Q states are shown in black and magenta, respectively. The singlet manifold is shown in orange, the triplet manifold in green and the quintet manifold in blue. The metal-centered (MC) states are shown explicitly, whereas the metal-to-ligand charge transfer (MLCT) states are depicted as a quasi-continuous band. Further, a schematic evolution of the wavepacket dynamics is shown. After absorption to the Q-band, the wavepacket undergoes large-amplitude oscillations in the Fe–CO coordinate on the ¹MLCT state. Then, it cascades down acquiring more ³MLCT and ⁵MLCT character as the spin crossover transitions occur. Finally, the wavepacket disperses, and the Fe–C(O) distance oscillations decrease to a value of 2.2 Å

[Back]

[Back]

Fig. 3.

Vibrational motions and effect on the Fe–CO interaction. Representation of the nuclear motions (left) and the effect on the Fe–CO orbital interaction (right). The σCO and dxz/dyz orbitals (in blue) are occupied, while ex/y, πx∕y,CO* and dz2 orbitals (in red) are unoccupied. Blue arrows indicate the direction of the partial electron transfer, due to the overlap between an occupied and an unoccupied orbital (solid and dashed arrows indicate strong and weak overlap). Magenta arrows indicate the vibrational motion of the specified atoms. a In the initial structure, the Fe–CO bond is strongest due to partial electron transfer dxz∕dyz→πx*∕πy* (π-back donation) and dz2←σCO (σ-donation). b Dissociative vibration, which corresponds to a stretching of the Fe–CO bond. This structure destabilizes the Fe–CO bond and stabilizes by a partial electron transfer dxz∕dyz→ex*∕ey* when the Fe atom is out of plane. c Rotational vibrations, which correspond to the rotation of the CO moiety. This vibration decreases the overlap between Fe and CO orbitals. d Symmetry-breaking vibration, in which the porphyrin is distorted inducing Jahn-Teller and pseudo Jahn-Teller vibronic couplings in the excited state. This vibration has no effect on the Fe–CO bond, but activates the coupling between the porphyrin and the Fe densities

[Back]