Resolving Femtosecond Solvent Reorganization Dynamics in an Iron Complex by Nonadiabatic Dynamics Simulations

Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations

Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.

Water Solvent Reorganization upon Ultrafast Resonant Stimulated X-ray Raman Excitation of a Metalloporphyrin Dimer

We propose an X-ray Raman pump-X-ray diffraction probe scheme to follow solvation dynamics upon charge migration in a solute molecule. The X-ray Raman pump selectively prepares a valence electronic wavepacket in the solute, while the probe provides information about the entire molecular ensemble. A combination of molecular dynamics and ab initio quantum chemistry simulations is applied to a Zn-Ni porphyrin dimer in water. Using time-resolved X-ray diffraction and pair distribution functions, we extracted solvation shell dynamics.

Vibrational Coherences of the Photoinduced Mixed-Valent Creutz-Taube Ion Revealed by Excited State Dynamics

A recent study of photoinduced mixed-valency in the one-electron reduced form (μ-pz)[RuII(NH3)5]24+ of the Creutz-Taube ion used transient absorption spectroscopy with vis-NIR broadband detection to uncover a mixed-valent excited state with a typical intervalence charge transfer band and a nanosecond lifetime [Pieslinger et al. Angew. Chem., Int. Ed. 2022, 61, e202211747]. Herein, we use excited state dynamics simulations with implicit solvation to elucidate the electronic and vibrational evolution in the first 10 ps after the optical excitation. A manifold of excited states with weak interaction between the metal centers is populated already at time zero due to the breakdown of the Condon approximation and dominates the population of electronic states at short time scales (<0.5 ps). A long-lived vibrational wave packet mostly confined to oscillations of the metal center-bridge distances is observed. The oscillations are traced to the electronic structure properties of states with weak metal-metal coupling. The long-lived mixed-valent excited state of the Creutz-Taube ion analogue is formed vibrationally cold and has a more compact geometry. While experimentally, intersystem crossing and vibrational relaxation were deduced to be completed within 1 ps, our analysis indicates that both processes might persist at longer times.

Eliminating finite-size effects on the calculation of x-ray scattering from molecular dynamics simulations

Structural studies using x-ray scattering methods for investigating molecules in solution are shifting focus toward describing the role and effects of the surrounding solvent. However, forward models based on molecular dynamics (MD) simulations to simulate structure factors and x-ray scattering from interatomic distributions such as radial distribution functions (RDFs) face limitations imposed by simulations, particularly at low values of the scattering vector q. In this work, we show how the value of the structure factor at q = 0 calculated from RDFs sampled from finite MD simulations is effectively dependent on the size of the simulation cell. To eliminate this error, we derive a new scheme to renormalize the sampled RDFs based on a model of the excluded volume of the particle-pairs they were sampled from, to emulate sampling from an infinite system. We compare this new correction method to two previous RDF-correction methods, developed for Kirkwood-Buff theory applications. We present a quantitative test to assess the reliability of the simulated low-q scattering signal and show that our RDF-correction successfully recovers the correct q = 0 limit for neat water. We investigate the effect of MD-sampling time on the RDF-corrections, before advancing to a molecular example system, comprised of a transition metal complex solvated in a series of water cells with varying densities. We show that our correction recovers the correct q = 0 behavior for all densities. Furthermore, we employ a simple continuum scattering model to dissect the total scattering signal from the solvent-solvent structural correlations in a solute-solvent model system to find two distinct contributions: a non-local density-contribution from the finite, fixed cell size in NVT simulations, and a local contribution from the solvent shell. We show how the second contribution can be approximated without also including the finite-size contribution. Finally, we provide a "best-practices"-checklist for experimentalists planning to incorporate explicit solvation MD simulations in future work, offering guidance for improving the accuracy and reliability of structural studies using x-ray scattering methods in solution.

Nanosecond Metal-to-Ligand Charge-Transfer State in an Fe(II) Chromophore: Lifetime Enhancement via Nested Potentials

Examples of Fe complexes with long-lived (≥1 ns) charge-transfer states are limited to pseudo-octahedral geometries with strong σ-donor chelates. Alternative strategies based on varying both coordination motifs and ligand donicity are highly desirable. Reported herein is an air-stable, tetragonal Fe^II complex, Fe(HMTI)(CN)₂ (HMTI = 5,5,7,12,12,14-hexamethyl-1,4,8,11-tetraazacyclotetradeca-1,3,8,10-tetraene), with a 1.25 ns metal-to-ligand charge-transfer (MLCT) lifetime. The structure has been determined, and the photophysical properties have been examined in a variety of solvents. The HMTI ligand is highly π-acidic due to low-lying π*(C═N), which enhances Δ_Fe via stabilizing t_2g orbitals. The inflexible geometry of the macrocycle results in short Fe-N bonds, and density functional theory calculations show that this rigidity results in an unusual set of nested potential energy surfaces. Moreover, the lifetime and energy of the MLCT state depends strongly on the solvent environment. This dependence is caused by modulation of the axial ligand-field strength by Lewis acid-base interactions between the solvent and the cyano ligands. This work represents the first example of a long-lived charge transfer state in an Fe^II macrocyclic species.

Photoactive Metal-to-Ligand Charge Transfer Excited States in 3d6 Complexes with Cr0, MnI, FeII, and CoIII

Many coordination complexes and organometallic compounds with the 4d⁶ and 5d⁶ valence electron configurations have outstanding photophysical and photochemical properties, which stem from metal-to-ligand charge transfer (MLCT) excited states. This substance class makes extensive use of the most precious and least abundant metal elements, and consequently there has been a long-standing interest in first-row transition metal compounds with photoactive MLCT states. Semiprecious copper(I) with its completely filled 3d subshell is a relatively straightforward and well explored case, but in 3d⁶ complexes the partially filled d-orbitals lead to energetically low-lying metal-centered (MC) states that can cause undesirably fast MLCT excited state deactivation. Herein, we discuss recent advances made with isoelectronic Cr⁰, Mn^I, Fe^II, and Co^III compounds, for which long-lived MLCT states have become accessible over the past five years. Furthermore, we discuss possible future developments in the search for new first-row transition metal complexes with partially filled 3d subshells and photoactive MLCT states for next-generation applications in photophysics and photochemistry.

Solvent Dynamics of Aqueous Halides before and after Photoionization

Electron transfer reactions can be strongly influenced by solvent dynamics. We study the photoionization of halides in water as a model system for such reactions. There are no internal nuclear degrees of freedom in the solute, allowing the dynamics of the solvent to be uniquely identified. We simulate the equilibrium solvent dynamics for Cl^-, Br^-, I^-, and their respective neutral atoms in water, comparing quantum mechanical/molecular mechanical (QM/MM) and classical molecular dynamics (MD) methods. On the basis of the obtained configurations, we calculate the extended X-ray absorption fine structure (EXAFS) spectra rigorously based on the MD snapshots and compare them in detail with other theoretical and experimental results available in the literature. We find our EXAFS spectra based on QM/MM MD simulations in good agreement with their experimental counterparts for the ions. Classical MD simulations for the ions lead to EXAFS spectra that agree equally well with the experiment when it comes to the oscillatory period of the signal, even though they differ from the QM/MM radial distribution functions extracted from the MD. The amplitude is, however, considerably overestimated. This suggests that to judge the reliability of theoretical simulation methods or to elucidate fine details of the atomistic dynamics of the solvent based on EXAFS spectra, the amplitude as well as the oscillatory period need to be considered. If simulations fail qualitatively, as does the classical MD for the aqueous neutral halogen atoms, the resulting EXAFS will also be strongly affected in both oscillatory period and amplitude. The good reliability of QM/MM-based EXAFS simulations, together with clear qualitative differences in the EXAFS spectra found between halides and their atomic counterparts, suggests that a combined theory and experimental EXAFS approach is suitable for elucidating the nonequilibrium solvent dynamics in the photoionization of halides and possibly also for electron transfer reactions in more complex systems.

Triplet Emitting C^N^C Cyclometalated Dibenzo[c,h]Acridine Pt(II) Complexes

In a series of Pt(II) complexes [Pt(dba)(L)] containing the very rigid, dianionic, bis-cyclometalating, tridentate C^N^C2− heterocyclic ligand dba2− (H2dba = dibenzo[c,h]acridine), the coligand (ancillary ligand) L = dmso, PPh3, CNtBu and Me2Imd (N,N’-dimethylimidazolydene) was varied in order to improve its luminescence properties. Beginning with the previously reported dmso complex, we synthesized the PPh3, CNtBu and Me2Imd derivatives and characterized them by elemental analysis, 1H (and 31P) NMR spectroscopy and MS. Cyclic voltammetry showed partially reversible reduction waves ranging between −1.89 and −2.10 V and increasing along the series Me2Imd < dmso ≈ PPh3 < CNtBu. With irreversible oxidation waves ranging between 0.55 (L = Me2Imd) and 1.00 V (dmso), the electrochemical gaps range between 2.65 and 2.91 eV while increasing along the series Me2Imd < CNtBu < PPh3 < dmso. All four complexes show in part vibrationally structured long-wavelength absorption bands peaking at around 530 nm. TD-DFT calculated spectra agree quite well with the experimental spectra, with only a slight redshift. The photoluminescence spectra of all four compounds are very similar. In fluid solution at 298 K, they show broad, only partially structured bands, with maxima at around 590 nm, while in frozen glassy matrices at 77 K, slightly blue-shifted (~580 nm) bands with clear vibronic progressions were found. The photoluminescence quantum yields ΦL ranged between 0.04 and 0.24, at 298 K, and between 0.80 and 0.90 at 77 K. The lifetimes τ at 298 K ranged between 60 and 14040 ns in Ar-purged solutions and increased from 17 to 43 µs at 77 K. The TD-DFT calculated emission spectra are in excellent agreement with the experimental findings. In terms of high ΦL and long τ, the dmso and PPh3 complexes outperform the CNtBu and Me2Imd derivatives. This is remarkable in view of the higher ligand strength of Me2Imd, compared with all other coligands, as concluded from the electrochemical data.