Search

Investigating the atomic interactions that govern solvation structure of transition metal complexes.

The time-resolved x-ray absorption spectrum of the BT-1T cation (BT-1T+) is theoretically simulated in order to investigate the charge transfer reaction of the system. We employ both trajectory surface hopping and quantum dynamics to simulate the structural evolution over time and the changes in the state populations. To compute the static x-ray absorption spectra (XAS) of the ground and excited states, we apply both the time-dependent density functional theory and the coupled cluster singles and doubles method. The results obtained are in good agreement between the methods. It is, furthermore, found that the small structural changes that occur during the reaction have little effect on the static XAS. Hence, the tr-XAS can be computed based on the state populations determined from a nuclear dynamics simulation and one set of static XAS calculations, utilizing the ground state optimized geometry. This approach can save considerable computational resources, as the static spectra need not to be calculated for all geometries. As BT-1T is a relatively rigid molecule, the outlined approach should only be considered when investigating non-radiative decay processes in the vicinity of the Franck-Condon point.

We present a sub-picosecond resolved investigation of the structural solvent reorganization and geminate recombination dynamics following 400 nm two-photon excitation and photodetachment of a valence p electron from the aqueous atomic solute, I

It has long been known that irradiation with visible light converts Fe(II) polypyridines from their low-spin (singlet) to high-spin (quintet) state, yet mechanistic interpreation of the photorelaxation remains contraversal. Herein, we simulate the full singlet-triplet-quintet dynamics of the [Fe(terpy)2]2+ (terpy = 2,2’:6’,2’’-terpyridine) complex in full dimension, in order to clarify the complex photodynamics. Importantly, we report a branching mechanism with two sequential components: 3MLCT → 3MC(3T2g) → 5MC and 3MLCT → 3MC(3T2g) → 3MC(3T1g) → 5MC (MLCT = metal-to-ligand charge transfer, MC = metal-centered). While the direct 3MLCT → 5MC mechanism is considered as a relevant alternative, we show that it could only be operative and thus lead to competing pathways in the absence of 3MC states. The quintet state is populated on the sub-picosecond timescale involving non-exponential dynamics and coherent Fe-N breathing oscillations. The results are in agreement with the available time-resolved experimental data on Fe(II) polypiridines, and fully describe the photorelaxation dynamics.

When characterizing transition metal complexes and their functionalities, the importance of including the solvent as an active participant is becoming more and more apparent. Whereas many studies have evaluated long-range dispersion effects inside organic molecules and organometallics, less is known about their role in solvation. Here, we have analysed the components within solute-solvent and solvent-solvent interactions of one of the most studied iron-based photoswitch model systems, in two spin states. We find that long-range dispersion effects modulate the coordination significantly, and that this is accurately captured by density functional theory models including dispersion corrections. We furthermore correlate gas-phase relaxed complex-water clusters to thermally averaged molecular densities. This shows how the gas-phase interactions translate to solution structure, quantified through 3D molecular densities, angular distributions, and radial distribution functions. We show that finite-size simulation cells can cause the radial distribution functions to have artificially enlarged amplitudes. Finally, we quantify the effects of many-body interactions within the solvent shells, and find that almost a fifth of the total interaction energy of the solute-shell system in the high-spin state comes from many-body contributions, which cannot be captured by by pair-wise additive force field methods.

A new theoretical approach is presented and applied for the simulation of Fe(II) low-spin (LS, singlet, t

We present a sub-picosecond resolved investigation of the structural solvent reorganization and geminate recombination dynamics following 400 nm two-photon excitation and photodetachment of a valence p electron from the aqueous atomic solute, I−(aq). The measurements utilized time-resolved X-ray Absorption Near Edge Structure (TR-XANES) spectroscopy and X-ray Solution Scattering (TR-XSS) at the Linac Coherent Light Source x-ray free electron laser in a laser pump/x-ray probe experiment. The XANES measurements around the L1-edge of the generated nascent iodine atoms (I0) yield an average electron ejection distance from the iodine parent of 7.4 ± 1.5 Å with an excitation yield of about 1/3 of the 0.1M NaI aqueous solution. The kinetic traces of the XANES measurement are in agreement with a purely diffusion-driven geminate iodine–electron recombination model without the need for a long-lived (I0:e−) contact pair. Nonequilibrium classical molecular dynamics simulations indicate a delayed response of the caging H2O solvent shell and this is supported by the structural analysis of the XSS data: We identify a two-step process exhibiting a 0.1 ps delayed solvent shell reorganization time within the tight H-bond network and a 0.3 ps time constant for the mean iodine–oxygen distance changes. The results indicate that most of the reorganization can be explained classically by a transition from a hydrophilic cavity with a well-ordered first solvation shell (hydrogens pointing toward I−) to an expanded cavity around I0 with a more random orientation of the H2O molecules in a broadened first solvation shell.

The photoisomerization reaction of azobenzene in both directions have been investigated with a density functional theory based approach using the surface hopping procedure with forced jumps. While the cis-to-trans isomerization was found to be a stepwise reaction along the CNNC dihedral angle, the trans-to-cis isomerization was observed to be one smooth step. The further unbiased full-dimensional analysis of the cis-to-trans isomerization revealed that, while the CNNC dihedral angle is an important degree of freedom for describing the reaction, it is insufficient for describing all of the dynamics. For a fuller picture two coupled modes must be considered. The trans-to-cis isomerization on the other hand was found to be well described along only the CNNC dihedral angle, and its longer timescale could be ascribed to the slow oscillations of this degree of freedom rather than a potential energy barrier in the excited state. The timescales observed in this study was found to be in good agreement with experiment, and thus this work provides insights into the interpretation of experimental observations. Finally, investigations of the structures of the CIs for both reactive and non-reactive trajectories showed a heavy functional dependency.

Inthis work, we focus on macrocyclic structures comprised of two dihydroazulene (DHA) units and one azobenzene (AZB) unit and the possibility for photoisomerizing one unit selectively by tuning the excitation energies of each individual unit. An unfortunate overlap between the absorption bands of DHA and AZB as well as trans-and cis-AZB prevents us to have a full control onthese macrocyclic structures, and their absorption bands need to be separated. By means of time-dependent density-functional theory calculations, we investigate the effects of ortho substitutions of the AZB unitby fluorine and chlorine atoms on the absorption spectra of the DHA/AZB macrocycles. The calculations on the isolated AZB show that substitutions lead to distortion of the planar molecular structure because of the repulsive interactions between halogen atoms and a systematic blueshift of the ππ* bands between 25 and 50 nm. Moreover, separations between 10 and 48 nm, depending on the substituent, are observed in the nπ* bands. The results from the calculations on the substituted AZB–DHA–DHA macrocycles reveal significant separations of the DHA/trans-AZB and trans-/cis-AZB absorption bands by values of 46–73 and 15–52 nm, respectively, for different substitutions. We realize that ortho substitutions with mixed fluorine–chlorine-atoms can provide the best separations in both ππ* and nπ* bands of AZB–DHA–DHA photoisomers. The results of this work offer a guideline for designing and synthesizing new, efficient, and highly controllable materials applicable in devices for optical data storage and molecular electronics.

The Variational Multi-Configurational Gaussian (vMCG) approach offers a framework to perform exact trajectory-based quantum dynamics. Herein we use two model vibronic coupling Hamiltonians of pyrazine to explore, for the first time, the influence of the coupling between the external field and the Gaussian basis functions (GBFs) in vMCG on the dynamics. We show that when the excitation pulse is short compared to the nuclear dynamics, vertical projection without a field and explicit description of the external field converge. For longer pulses, a sizeable change is observed. We demonstrate that comparatively few GBFs are sufficient to provide qualitative agreement to MCTDH dynamics and a quantitative agreement can be achieved using ∼100 GBFs. Longer pulses require more GBFs due to the prolonged coupling between the ground and excited states. Throughout this work the single set formalism offers the fastest convergence.

We present a computational study on nonadiabatic excited-state dynamics initiated from the 3p Rydberg states of trimethylamine (TMA). We utilise a methodology based on full-dimensional (39 D) trajectory surface-hopping (TSH) simulations, in which propagation is carried out on on-the-fly density functional theory (DFT)/time-dependent DFT (TD-DFT) potentials. Both our electronic structure benchmarks to high-level ab initio methods (EOM-CCSD, CASPT2) and TSH simulations demonstrate high-accuracy of the applied CAM-B3LYP functional for the description of Rydberg excited states. Based on our excited-state simulations, we construct the following mechanistic picture: when pumped resonantly to the 3p Rydberg manifold, TMA coherently vibrates along the planarisation mode with a period of 104 fs and an exponential coherence decay time constant of 240 fs. Nonadiabatic dynamics occur on a faster (∼1 ps) and a slower (∼3 ps) timescale, along the N-C stretching mode by mixing with a dissociative σN-C* state. As a minor relaxation channel, 3p → 3s internal conversion occurs via branching at the σN-C*/3s intersection. We find that photodissociaton is hardly observable within 3 ps (1%), which is a failure of the range-separated hybrid CAM-B3LYP functional, as a consequence of its static electron correlation deficiency at long range. In contrast, pure DFT (GGA-BLYP) provides an accurate long-range description (19% dissociation yield), also supported by comparison to recent ultrafast experiments, even if the Rydberg state energies are significantly underestimated (>1 eV). Finally, we reveal the crucial role of vibrational coherence and energy transfer from the planarisation mode for N-C bond activation and resulting nonadiabatic dynamics. The present work illustrates the importance of nuclear-electronic coupling for excited-state dynamics and branching at conical intersections.

The photoinduced ring-opening reaction is a key process in the functioning of dihydroazulene/vinylheptafulvene (DHA/VHF) photoswitches. Over the years, the mechanism of this reaction has been extensively debated. Herein, by means of nonadiabatic trajectory dynamics simulations and quantum chemistry calculations, we present the first detailed and comprehensive investigation on the mechanism of the photoinduced ring-opening reaction of DHA. The results show the crucial role of the excited-state ring planarization process for the bond breaking. Our dynamics simulations show that the DHA ring opening is an ultrafast reaction that does not follow exponential kinetics but exhibits ballistic dynamics. Upon photoexcitation, the planarization occurs within 300-500 fs. This leads to the ring-opening reaction and concurrent decay of the molecule to the ground state within 100 fs through an S1 → S0 internal conversion process toward forming the VHF isomer. These results are consistent with previous ultrafast time-resolved experiments and lead to a thorough understanding of the DHA/VHF photoconversion.

We have employed a range of ultrafast X-ray spectroscopies in an effort to characterize the lowest energy excited state of [Fe(dcpp)2]2+ (where dcpp is 2,6-(dicarboxypyridyl)pyridine). This compound exhibits an unusually short excited-state lifetime for a low-spin Fe(II) polypyridyl complex of 270 ps in a room-temperature fluid solution, raising questions as to whether the ligand-field strength of dcpp had pushed this system beyond the 5T2/3T1 crossing point and stabilizing the latter as the lowest energy excited state. Kα and Kβ X-ray emission spectroscopies have been used to unambiguously determine the quintet spin multiplicity of the long-lived excited state, thereby establishing the 5T2 state as the lowest energy excited state of this compound. Geometric changes associated with the photoinduced ligand-field state conversion have also been monitored with extended X-ray absorption fine structure. The data show the typical average Fe-ligand bond length elongation of ∼0.18 Å for a 5T2 state and suggest a high anisotropy of the primary coordination sphere around the metal center in the excited 5T2 state, in stark contrast to the nearly perfect octahedral symmetry that characterizes the low-spin 1A1 ground state structure. This study illustrates how the application of time-resolved X-ray techniques can provide insights into the electronic structures of molecules-in particular, transition metal complexes-that are difficult if not impossible to obtain by other means.

We have employed a range of ultrafast X-ray spectroscopies in an effort to characterize the lowest energy excited state of [Fe(dcpp)

The solvent-mediatedexcited-state dynamics of the COOH-functionalized Fe-carbene photosensitizer [Fe(bmicp)2]2+ (bmicp= 2,6-bis(3-methyl-imidazole-1-ylidine)-4-carboxy-pyridine) is studied by time-dependent density functional theory, as well as classical and quantum dynamics simulations. We demonstrate the crucial role of the polar acetonitrile solvent in stabilizing the metal-to-ligand charge transfer (MLCT) states of the investigated molecule using the conductor polarizable continuum model. This leads to dynamics that avoid sub-ps back electron transfer to the metal and an exceptionally long-lived 1MLCT state that does not undergo sub-ps 1MLCT → 3MLCT intersystem crossing as it is energetically isolated. We identify two components of the excited-state solvent reorganization process: an initial rotation (∼300 fs) and diffusional dynamics within the local cage surrounding the rotated solvent molecule (∼2 ps). Finally, it is found that the relaxation of the solvent only slightly affects the excited-state population dynamics of [Fe(bmicp)2]2+.

Publisher's version (útgefin grein), Light-driven molecular reactions are dictated by the excited state potential energy landscape, depending critically on the location of conical intersections and intersystem crossing points between potential surfaces where non-adiabatic effects govern transition probabilities between distinct electronic states. While ultrafast studies have provided significant insight into electronic excited state reaction dynamics, experimental approaches for identifying and characterizing intersections and seams between electronic states remain highly system dependent. Here we show that for 3d transition metal systems simultaneously recorded X-ray diffuse scattering and X-ray emission spectroscopy at sub-70 femtosecond time-resolution provide a solid experimental foundation for determining the mechanistic details of excited state reactions. In modeling the mechanistic information retrieved from such experiments, it becomes possible to identify the dominant trajectory followed during the excited state cascade and to determine the relevant loci of intersections between states. We illustrate our approach by explicitly mapping parts of the potential energy landscape dictating the light driven low-to-high spin-state transition (spin crossover) of [Fe(2,2′-bipyridine)3]2+, where the strongly coupled nuclear and electronic dynamics have been a source of interest and controversy. We anticipate that simultaneous X-ray diffuse scattering and X-ray emission spectroscopy will provide a valuable approach for mapping the reactive trajectories of light-triggered molecular systems involving 3d transition metals., Use of the Linac Coherent Light Source (LCLS), SLAC National Accelerator Laboratory, is supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC02-76SF00515. Work by KK, MER, and KJG was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division. KL acknowledges a Melvin and Joan Lane Stanford Graduate Fellowship. KSK, KH, MMN, MGL, PV, EB, FBH, MC gratefully acknowledge DANSCATT support for the beamtime efforts. TCBH, EB, MGL, MIP, KBM, and MMN gratefully acknowledge support by the Danish Council for Independent Research under grant no. DFF-4002-00272B. MIP, MMN, and KBM acknowledge support by the Danish Council for Independent Research under grant no. DFF-8021-00347B. KSK gratefully acknowledge the support of the Carlsberg Foundation and the Danish Council for Independent Research. MIP, DSS, EB, and GV acknowledge support from the ‘Lendület’ (Momentum) Program of the Hungarian Academy of Sciences (LP2013-59), the Government of Hungary and the European Regional Development Fund under grant No. VEKOP-2.3.2-16-2017-00015, the European Research Council via contract ERC-StG-259709 (X-cited!), the Hungarian Scientific Research Fund (OTKA) under contract K 109257, and the National Research, Development and Innovation Fund (NKFIH FK 124460). ZN acknowledges support from the Bolyai Fellowship of the Hungarian Academy of Sciences. SC acknowledges the ELI-ALPS project (GINOP-2.3.6-15-2015-00001) which is supported by the European Union and co-financed by the European Regional Development Fund. VS, JU, and PP acknowledge support from the Knut and Alice Wallenberg Foundation (KAW). SK acknowledges the support from Knut & Alice Wallenberg foundation (KAW 2014.0370). AOD acknowledges support by the Icelandic Research Fund (Grant# 196279-051).

Simulation of the ultrafast excited-state dynamics and elastic X-ray scattering of the [Fe(bmip)2]2+ [bmip = 2,6-bis(3-methyl-imidazole-1-ylidine)-4-pyridine] complex is presented and analyzed. We employ quantum wavepacket dynamics simulations on a 5-dimensional potential energy surface (PES) calculated by time-dependent density functional theory with 26 coupled diabatic states. The simulations are initiated by explicit inclusion of a time-dependent electromagnetic field. In the case of resonant excitation into singlet metal-to-ligand charge transfer (1MLCT) states, kinetic (exponential) population dynamics are observed with small nuclear motion. In agreement with transient optical absorption spectroscopy experiments, we observe a subpicosecond 1MLCT → 3MLCT intersystem crossing and a subsequent decay into triplet metal-centered (3MC) states on a picosecond time scale. The simulated time-resolved difference scattering signal is dominated by the 3MC component, for which the structural distortions are significant. On the other hand, excitation into 1MC states leads to ballistic (nonexponential) population dynamics with strong nuclear motion. The reason for these ballistic dynamics is that in this case, the excitation occurs into a nonequilibrium region, i.e., far from the minimum of the 1MC PES. This results in wavepacket dynamics along the principal breathing mode, which is clearly visible in both the population dynamics and difference scattering. Finally, the importance of decomposing the difference scattering into components by electronic states is highlighted, information which is not accessible from elastic X-ray scattering experiments.

In Fe[N(SiPh2Me)2]2 molecular magnet, an extraordinary B ≈ 92 T hyperfine field was found in the 5 K 57Fe Mossbauer spectrum under an external magnetic field of 0.1 T. This evidences the presence of an unquenched orbital angular moment at the central iron atom. Fe[N(SiPh2Me)2]2 complex is thus shown to represent a further example of low-coordinate iron complexes where quasi free-ion magnetism visualizes itself through an unquenched orbital moment. Magnetization measurements and hysteresis in magnetization indicated exchange coupling and nanosized magnetic units.

Understanding and subsequently being able to manipulate the excited-state decay pathways of functional transition-metal complexes is of utmost importance in order to solve grand challenges in solar energy conversion and data storage. Herein, we perform quantum chemical calculations and spin-vibronic quantum dynamics simulations on the Fe-N-heterocyclic carbene complex, [Fe(btbip)2]2+ (btbip = 2,6-bis(3-tert-butyl-imidazole-1-ylidene)pyridine). The results demonstrate that a relatively minor structural change compared to its parent complex, [Fe(bmip)2]2+ (bmip = 2,6-bis(3-methyl-imidazole-1-ylidene)pyridine), completely alters the excited-state relaxation. Ultrafast deactivation of the initially excited metal-to-ligand charge transfer (1,3MLCT) states occurs within 350 fs. In contrast to the widely adopted mechanism of Fe(II) photophysics, these states decay into close-lying singlet metal-centered (1MC) states. This occurs because the tert-butyl functionalization stabilizes the 1MC states, enabling the 1,3MLCT → 1MC population transfer to occur close to the Franck-Condon geometry, making the conversion very efficient. Subsequently, a spin cascade occurs within the MC manifold, leading to the population of triplet and quintet MC states. These results will inspire highly involved ultrafast experiments performed at X-ray free electron lasers and shall pave the way for the design of novel high-efficiency transition-metal-based functional molecules.

In liquid phase chemistry dynamic solute-solvent interactions often govern the path, ultimate outcome, and efficiency of chemical reactions. These steps involve many-body movements on subpicosecond time scales and thus ultrafast structural tools capable of capturing both intramolecular electronic and structural changes, and local solvent structural changes are desired. We have studied the intra- and intermolecular dynamics of a model chromophore, aqueous [Fe(bpy)3](2+), with complementary X-ray tools in a single experiment exploiting intense XFEL radiation as a probe. We monitored the ultrafast structural rearrangement of the solute with X-ray emission spectroscopy, thus establishing time zero for the ensuing X-ray diffuse scattering analysis. The simultaneously recorded X-ray diffuse scattering patterns reveal slower subpicosecond dynamics triggered by the intramolecular structural dynamics of the photoexcited solute. By simultaneous combination of both methods only, we can extract new information about the solvation dynamic processes unfolding during the first picosecond (ps). The measured bulk solvent density increase of 0.2% indicates a dramatic change of the solvation shell around each photoexcited solute, confirming previous ab initio molecular dynamics simulations. Structural changes in the aqueous solvent associated with density and temperature changes occur with ∼1 ps time constants, characteristic for structural dynamics in water. This slower time scale of the solvent response allows us to directly observe the structure of the excited solute molecules well before the solvent contributions become dominant.

Ultrafast pump-probe spectroscopy within the X-ray regime is now possible owing to the development of X-ray Free Electrons Lasers (X-FELs) and is opening new opportunities for the direct probing of femtosecond evolution of the nuclei, the electronic and spin degrees of freedom. In this contribution we use wavepacket dynamics of the photoexcited decay of a new Fe(ii) complex, [Fe(bmip)2]2+ (bmip = 2,6-bis(3-methyl-imidazole-1-ylidine)pyridine), to simulate the experimental observables associated with femtosecond Fe K-edge X-ray Absorption Near-Edge Structure (XANES) and X-ray emission (XES) spectroscopy. We show how the evolution of the nuclear wavepacket is translated into the spectroscopic signal and the sensitivity of these approaches for following excited state dynamics.

Recent ultrafast experiments have unveiled the time scales of vibrational cooling and decoherence upon photoexcitation of the diplatinum complex [Pt2(P2O5H2)4]4– in solvents. Here, we contribute to the understanding of the structure and dynamics of the lowest lying singlet excited state of the model photocatalyst by performing potential energy surface calculations and Born–Oppenheimer molecular dynamics simulations in the gas phase and in water. Solvent effects were treated using a multiscale quantum mechanics/molecular mechanics approach. Fast sampling was achieved with a modified version of delta self-consistent field implemented in the grid-based projector-augmented wave density functional theory code. The known structural parameters and the PESs of the first singlet and triplet excited states are correctly reproduced. Besides, the simulations deliver clear evidence that pseudorotation of the ligands in the excited state leads to symmetry lowering of the Pt2P8 core. Coherence decay of Pt–Pt stretching vibrations in solution was found to be governed by vibrational cooling, which is in agreement with previous ultrafast experiments. We also show that the flow of excess Pt–Pt vibrational energy is first directed toward vibrational modes involving the ligands, with the solvent favoring intramolecular vibrational energy redistribution. The results are supported by thorough vibrational analysis in terms of generalized normal modes.

The excited-state dynamics of two functional Fe-carbene complexes, [Fe(bmip)2]2+ (bmip = 2,6-bis(3-methyl-imidazole-1-ylidene)-pyridine) and [Fe(btbip)2]2+ (btbip = 2,6-bis(3-tert-butyl-imidazole-1-ylidene)pyridine), are studied using the spin-vibronic model. In contrast to the usual projection of the ground state nuclear wavefunction onto an excited state surface, the dynamics are initiated by an explicit interaction term between the external time-dependent electric field (laser pulse) and the transition dipole moment of the molecule. The results show that the spin-vibronic model, as constructed directly from electronic structure calculations, exhibits erroneous, polarization-dependent relaxation dynamics stemming from artificial interference of coupled relaxation pathways. This is due to the lack of rotational invariance in the description of excitation into degenerate states. We introduce and discuss a correction using the spherical basis and complex transition dipole moments. This modification in the Hamiltonian leads to rotationally invariant excitation and produces polarization-independent population dynamics.

The journal of physical chemistry / C 121(5), 2620 - 2626 (2017). doi:10.1021/acs.jpcc.6b12940, We probe the dynamics of valence electrons in photoexcited [Fe(terpy)$_2$]$^{2+}$ in solution to gain deeper insight into the Fe−ligand bond changes. We use hard X-rayemission spectroscopy (XES), which combines element specificity and high penetration with sensitivity to orbital structure, making it a powerful technique for molecular studies in a wide variety of environments. A picosecond-time-resolved measurement of the complete 1s X-ray emission spectrum captures the transient photoinduced changes and includes the weak valence-to-core (vtc) emission lines that correspond to transitions from occupied valence orbitals to the nascent core-hole. Vtc-XES offers particular insight into the molecular orbitals directly involved in the light-driven dynamics; a change in the metal− ligand orbital overlap results in an intensity reduction and a blue energy shift in agreement with our theoretical calculations and more subtle features at the highest energies reflect changes in the frontier orbital populations., Published by ACS, Washington, DC

Magnetization and low temperature in-field 57Fe Mossbauer spectroscopy measurements have been performed on a Fe[C(SiMe3)3]2/ferrihydrite hetero-mixture. The results indicate the presence of ferromagnetic coupling of magnetic moments involving Fe[C(SiMe3)3]2 with a hyperfine magnetic field of about 151 T, attributable mainly to the non-frozen atomic orbital contribution. The present findings show the sensitivity of single-ion molecular magnets to local alterations of their lattice-environment and might explain and reconcile some of the differences found in the literature for the observed bulk magnetic properties of the title Fe[C(SiMe3)3]2 compound.

Earth-abundant transition-metal complexes are desirable for sensitizers in dye-sensitized solar cells or photocatalysts. Iron is an obvious choice, but the energy level structure of its typical polypyridyl complexes, featuring low-lying metal-centered states, has made such complexes useless as energy converters. Recently, we synthesized a novel iron-N-heterocyclic carbene complex exhibiting a remarkable 100-fold increase of the lifetime compared to previously known iron(II) complexes. Here, we rationalize the measured excited-state dynamics with DFT and TD-DFT calculations. The calculations show that the exceptionally long excited-state lifetime (∼9 ps) is achieved for this Fe complex through a significant destabilization of both triplet and quintet metal-centered scavenger states compared to other Fe(II) complexes. In addition, a shallow (3)MLCT potential energy surface with a low-energy transition path from the (3)MLCT to (3)MC and facile crossing from the (3)MC state to the ground state are identified as key features for the excited-state deactivation.

As a demonstration of the analysis of the electronic structure and the nuclear dynamics from time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS), we present the TR-NEXAFS spectra of pyrazine following the excitation to the 1B2u(ππ*) state. The spectra are calculated combining the frozen-core/core-valence separated equation-of-motion coupled cluster singles and doubles approach for the spectral signatures and the multiconfiguration time-dependent Hartree method for the wave packet propagation. The population decay from the 1B2u(ππ*) state to the 1B3u(nπ*) and 1Au(nπ*) states, followed by oscillatory flow of population between the 1B3u(nπ*) and 1Au(nπ*) states, is interpreted by means of visualization of the potential energy curves and the reduced nuclear densities. By examining the electronic structure of the three valence-excited states and the final core-excited states, we observe that the population dynamics is explicitly reflected in the TR-NEXAFS spectra, especially when the heteroatoms are selected as the X-ray absorption sites. This work illustrates the feasibility of extracting fine details of molecular photophysical processes from TR-NEXAFS spectra by using currently available theoretical methods.

We study the structural dynamics of photoexcited [Co(terpy)_{2}]^{2+} in an aqueous solution with ultrafast x-ray diffuse scattering experiments conducted at the Linac Coherent Light Source. Through direct comparisons with density functional theory calculations, our analysis shows that the photoexcitation event leads to elongation of the Co-N bonds, followed by coherent Co-N bond length oscillations arising from the impulsive excitation of a vibrational mode dominated by the symmetrical stretch of all six Co-N bonds. This mode has a period of 0.33 ps and decays on a subpicosecond time scale. We find that the equilibrium bond-elongated structure of the high spin state is established on a single-picosecond time scale and that this state has a lifetime of ∼7 ps.

Fe(II) complexes have long been assumed unsuitable as photosensitizers because of their low-lying nonemissive metal centered (MC) states, which inhibit electron transfer. Herein, we describe the excited-state relaxation of a novel Fe(II) complex that incorporates N-heterocyclic carbene ligands designed to destabilize the MC states. Using first-principles quantum nuclear wavepacket simulations we achieve a detailed understanding of the photoexcited decay mechanism, demonstrating that it is dominated by an ultrafast intersystem crossing from 1MLCT–3MLCT proceeded by slower kinetics associated with the conversion into the 3MC states. The slowest component of the 3MLCT decay, important in the context of photosensitizers, is much longer than related Fe(II) complexes because the population transfer to the 3MC states occurs in a region of the potential where the energy gap between the 3MLCT and 3MC states is large, making the population transfer inefficient.

The journal of physical chemistry / C 120(35), 19537 - 19546(2016). doi:10.1021/acs.jpcc.6b06005, Metalloporphyrins are prominent building blocks in the synthetic toolbox of advanced photodriven molecular devices. When the central ion is paramagnetic, the relaxation pathways within the manifold of excited states are highly intricate so that unravelling the intramolecular energy and electron transfer processes is usually a very complex task. This fact is critically hampering the development of applications based on the enhanced coupling offered by the electronic exchange interaction. In this work, the dynamics of charge separation in a copper porphyrin-fullerene are studied with several complementary spectroscopic tools across the electromagnetic spectrum (from near-infrared to X-ray wavelengths), each of them providing specific diagnostics. Correlating the various rates clearly demonstrates that the lifetime of the photoinduced charge-separated state exceeds by about 10-fold that of the isolated photoexcited CuII porphyrin. As revealed by the spectral modifications in the XANES region, this stabilization is accompanied by a transient change in covalency around the CuII center, which is induced by an enhanced interaction with the C60 moiety. This experimental finding is further confirmed by state-of-the art calculations using DFT and TD-DFT including dispersion effects that explain the electrostatic and structural origins of this interaction, as the CuIIP cation becomes ruffled and approaches closer to the fullerene in the charge-separated state. From a methodological point of view, these results exemplify the potential of multielectron excitation features in transient X-ray spectra as future diagnostics of subfemtosecond electronic dynamics. From a practical point of view, this work is paving the way for elucidating out-of-equilibrium electron transfer events coupled to magnetic interaction processes on their intrinsic time-scales., Published by Soc., Washington, DC

The photochemical reactions performed by transition metal complexes have been proposed as viable routes towards solar energy conversion and storage into other forms that can be conveniently used in our everyday applications. In order to develop efficient materials, it is necessary to identify, characterize and optimize the elementary steps of the entire process on the atomic scale. To this end, we have studied the photoinduced electronic and structural dynamics in two heterobimetallic ruthenium–cobalt dyads, which belong to the large family of donor–bridge–acceptor systems. Using a combination of ultrafast optical and X-ray absorption spectroscopies, we can clock the light-driven electron transfer processes with element and spin sensitivity. In addition, the changes in local structure around the two metal centers are monitored. These experiments show that the nature of the connecting bridge is decisive for controlling the forward and the backward electron transfer rates, a result supported by quantum chemistry calculations. More generally, this work illustrates how ultrafast optical and X-ray techniques can disentangle the influence of spin, electronic and nuclear factors on the intramolecular electron transfer process. Finally, some implications for further improving the design of bridged sensitizer-catalysts utilizing the presented methodology are outlined.

Building a detailed understanding of the structure function relationship is a crucial step in the optimization of molecular photocatalysts employed in water splitting schemes. The optically dark nature of their active sites usually prevents a complete mapping of the photoinduced dynamics. In this work, transient X-ray absorption spectroscopy highlights the electronic and geometric changes that affect such a center in a bimetallic model complex. Upon selective excitation of the ruthenium chromophore, the cobalt moiety is reduced through intramolecular electron transfer and undergoes a spin flip accompanied by an average bond elongation of 0.20 +/- 0.03 angstrom. The analysis is supported by simulations based on density functional theory structures (B3LYP*/TZVP) and FEFF 9.0 multiple scattering calculations. More generally, these results exemplify the large potential of the technique for tracking elusive intermediates that impart unique functionalities in photochemical devices.

Theoretical predictions show that depending on the populations of the Fe 3dxy, 3dxz, and 3dyz orbitals two possible quintet states can exist for the high-spin state of the photoswitchable model system [Fe(terpy)2]2+. The differences in the structure and molecular properties of these 5B2 and 5E quintets are very small and pose a substantial challenge for experiments to resolve them. Yet for a better understanding of the physics of this system, which can lead to the design of novel molecules with enhanced photoswitching performance, it is vital to determine which high-spin state is reached in the transitions that follow the light excitation. The quintet state can be prepared with a short laser pulse and can be studied with cutting-edge time-resolved X-ray techniques. Here we report on the application of an extended set of X-ray spectroscopy and scattering techniques applied to investigate the quintet state of [Fe(terpy)2]2+ 80 ps after light excitation. High-quality X-ray absorption, nonresonant emission, and resonant emission spectra as well as X-ray diffuse scattering data clearly reflect the formation of the high-spin state of the [Fe(terpy)2]2+ molecule; moreover, extended X-ray absorption fine structure spectroscopy resolves the Fe–ligand bond-length variations with unprecedented bond-length accuracy in time-resolved experiments. With ab initio calculations we determine why, in contrast to most related systems, one configurational mode is insufficient for the description of the low-spin (LS)–high-spin (HS) transition. We identify the electronic structure origin of the differences between the two possible quintet modes, and finally, we unambiguously identify the formed quintet state as 5E, in agreement with our theoretical expectations.

Theoretical predictions show that depending on the populations of the Fe 3dxy, 3dxz, and 3dyz orbitals two possible quintet states can exist for the high-spin state of the photoswitchable model system [Fe(terpy)2]2+. The differences in the structure and molecular properties of these 5B2 and 5E quintets are very small and pose a substantial challenge for experiments to resolve them. Yet for a better understanding of the physics of this system, which can lead to the design of novel molecules with enhanced photoswitching performance, it is vital to determine which high-spin state is reached in the transitions that follow the light excitation. The quintet state can be prepared with a short laser pulse and can be studied with cutting-edge time-resolved X-ray techniques. Here we report on the application of an extended set of X-ray spectroscopy and scattering techniques applied to investigate the quintet state of [Fe(terpy)2]2+ 80 ps after light excitation. High-quality X-ray absorption, nonresonant emission, and resonant emission spectra as well as X-ray diffuse scattering data clearly reflect the formation of the high-spin state of the [Fe(terpy)2]2+ molecule; moreover, extended X-ray absorption fine structure spectroscopy resolves the Fe–ligand bond-length variations with unprecedented bond-length accuracy in time-resolved experiments. With ab initio calculations we determine why, in contrast to most related systems, one configurational mode is insufficient for the description of the low-spin (LS)–high-spin (HS) transition. We identify the electronic structure origin of the differences between the two possible quintet modes, and finally, we unambiguously identify the formed quintet state as 5E, in agreement with our theoretical expectations.

X-ray spectroscopies, when combined in laser-pump, X-ray-probe measurement schemes, can be powerful tools for tracking the electronic and geometric structural changes that occur during the course of a photoinitiated chemical reaction. X-ray absorption spectroscopy (XAS) is considered an established technique for such measurements, and X-ray emission spectroscopy (XES) of the strongest core-to-core emission lines (Kα and Kβ) is now being utilized. Flux demanding valence-to-core XES promises to be an important addition to the time-resolved spectroscopic toolkit. In this paper we present measurements and density functional theory calculations on laser-excited, solution-phase ferrocyanide that demonstrate the feasibility of valence-to-core XES for time-resolved experiments. We discuss technical improvements that will make valence-to-core XES a practical pump-probe technique.

The electronic structure relevant to low spin (LS)↔high spin (HS) transitions in Fe(II) coordination compounds with a FeN6 core are studied. The selected [Fe(tz)6]2+ (1) (tz = 1H-tetrazole), [Fe(bipy)3]2+ (2) (bipy = 2,2′-bipyridine), and [Fe(terpy)2]2+ (3) (terpy = 2,2′:6′,2″-terpyridine) complexes have been actively studied experimentally, and with their respective mono-, bi-, and tridentate ligands, they constitute a comprehensive set for theoretical case studies. The methods in this work include density functional theory (DFT), time-dependent DFT (TD-DFT), and multiconfigurational second order perturbation theory (CASPT2). We determine the structural parameters as well as the energy splitting of the LS–HS states (ΔEHL) applying the above methods and comparing their performance. We also determine the potential energy curves representing the ground and low-energy excited singlet, triplet, and quintet d6 states along the mode(s) that connect the LS and HS states. The results indicate that while DFT is well suited for the prediction of structural parameters, an accurate multiconfigurational approach is essential for the quantitative determination of ΔEHL. In addition, a good qualitative agreement is found between the TD-DFT and CASPT2 potential energy curves. Although the TD-DFT results might differ in some respect (in our case, we found a discrepancy at the triplet states), our results suggest that this approach, with due care, is very promising as an alternative for the very expensive CASPT2 method. Finally, the two-dimensional (2D) potential energy surfaces above the plane spanned by the two relevant configuration coordinates in [Fe(terpy)2]2+ were computed at both the DFT and CASPT2 levels. These 2D surfaces indicate that the singlet–triplet and triplet–quintet states are separated along different coordinates, i.e., different vibration modes. Our results confirm that in contrast to the case of complexes with mono- and bidentate ligands, the singlet–quintet transitions in [Fe(terpy)2]2+ cannot be described using a single configuration coordinate.

The performance of six frequently used density functional theory (DFT) methods (RPBE, OLYP, TPSS, B3LYP, B3LYP*, and TPSSh) in the prediction of Mössbauer isomer shifts(δ) and quadrupole splittings (ΔEQ) is studied for an extended and diverse set of Fe complexes. In addition to the influence of the applied density functional and the type of the basis set, the effect of the environment of the molecule, approximated with the conducting-like screening solvation model (COSMO) on the computed Mössbauer parameters, is also investigated. For the isomer shifts the COSMO-B3LYP method is found to provide accurate δ values for all 66 investigated complexes, with a mean absolute error (MAE) of 0.05 mm s–1 and a maximum deviation of 0.12 mm s–1. Obtaining accurate ΔEQ values presents a bigger challenge; however, with the selection of an appropriate DFT method, a reasonable agreement can be achieved between experiment and theory. Identifying the various chemical classes of compounds that need different treatment allowed us to construct a recipe for ΔEQ calculations; the application of this approach yields a MAE of 0.12 mm s–1 (7% error) and a maximum deviation of 0.55 mm s–1 (17% error). This accuracy should be sufficient for most chemical problems that concern Fe complexes. Furthermore, the reliability of the DFT approach is verified by extending the investigation to chemically relevant case studies which include geometric isomerism, phase transitions induced by variations of the electronic structure (e.g., spin crossover and inversion of the orbital ground state), and the description of electronically degenerate triplet and quintet states. Finally, the immense and often unexploited potential of utilizing the sign of the ΔEQ in characterizing distortions or in identifying the appropriate electronic state at the assignment of the spectral lines is also shown.

The sensitivity of the 1s X-ray emission and high-energy-resolution fluorescence-detected X-ray absorption spectroscopies (XES and HERFD-XAS) to resolve the variations in the chemical state (electronic structure and local coordination) of Br has been investigated for a selected set of compounds including NaBrO3, NH4Br and C2H4Br2 (1,2-dibromoethane). For the Br K-edge XAS, employing the HERFD mode significantly increases the energy resolution, which demonstrates that with a crystal spectrometer used as a detector the absorption technique becomes a more powerful analytical tool. In the case of XES, the experimental results as well as the density functional theory (DFT) modeling both show that the chemical sensitivity of the main 1s diagram emission lines (Kα1,2 and Kβ1,3) is rather limited. However, the valence-to-core (Kβ2) region of XES displays significant shape and intensity variations, as expected for transitions having the same final states as those of photoemission spectroscopy. The spectra are in good agreement with the molecular orbital description delivered by DFT calculations. Calculations for an extended series of Br compounds confirm that valence-to-core XES can serve as a probe for chemical analysis, and, being a hard X-ray photon-in/photon-out technique, it is particularly well-suited for in situ investigations of molecular transformations, even on the ultrafast time scales down to femtosecond time resolution.

Characterizing the geometric and electronic structures of individual photoexcited dye molecules in solution is an important step towards understanding the interfacial properties of photo-active electrodes. The broad family of "red sensitizers" based on osmium(II) polypyridyl compounds often undergoes small photo-induced structural changes which are challenging to characterize. In this work, X-ray transient absorption spectroscopy with picosecond temporal resolution is employed to determine the geometric and electronic structures of the photoexcited triplet state of [Os(terpy)2]2+ (terpy: 2,2':6',2″-terpyridine) solvated in methanol. From the EXAFS analysis, the structural changes can be characterized by a slight overall expansion of the first coordination shell [OsN6]. DFT calculations supports the XTA results. They also provide additional information about the nature of the molecular orbitals that contribute to the optical spectrum (with TD-DFT) and the near-edge region of the X-ray spectra.