Simulation of Ultrafast Excited-State Dynamics in Fe(II) Complexes: Assessment of Electronic Structure Descriptions

Mátyás Pápai
DOI: https://doi.org/10.26434/chemrxiv-2024-sfn8p
2024-09-18
Abstract:The assessment of electronic structure descriptions utilized in the simulation of ultrafast excited-state dynamics of Fe(II) complexes is presented. Herein, we evaluate the performance of the RPBE, OPBE, BLYP, B3LYP, B3LYP*, CAM-B3LYP, and LC-BLYP (time-dependent) density functional theory (DFT/TD-DFT) methods in full-dimensional trajectory surface hopping (TSH) simulations carried out on linear vibronic coupling (LVC) potentials. We exploit the existence of time-resolved X-ray emission spectroscopy (XES) data for the [Fe(bmip)2]2+ and [Fe(terpy)2]2+ prototypes for dynamics between metal-to-ligand charge transfer (MLCT) and metal-centered (MC) states, which serve as a reference to benchmark the calculations (bmip = 2,6-bis(3-methyl-imidazole-1-ylidine)-pyridine, terpy = 2,2’:6’,2''-terpyridine). The results show that the simulated ultrafast population dynamics between MLCT and MC states with various spin multiplicilities (singlet, triplet, quintet) highly depend on the utilized DFT/TD-DFT method with the percentage of exact (Hartree-Fock) exchange being the governing factor. Importantly, B3LYP* is the only DFT/TD-DFT method reproducing all important aspects of the experimentally resolved dynamics for both complexes, signalling an optimal balance in the description of MLCT-MC energetics. This work demonstrates the power of combining TSH/LVC dynamics simulations with time-resolved experimental reference data to benchmark full-dimensional potential energy surfaces.
Chemistry
What problem does this paper attempt to address?
The paper aims to address the challenges faced in the ultrafast excited-state dynamics simulation of iron (II) complexes. Specifically: 1. **Evaluation of Electronic Structure Description Methods**: The paper evaluates the performance of several Density Functional Theory (DFT) and its time-dependent version (TD-DFT) methods in the ultrafast excited-state dynamics simulation of iron (II) complexes. These methods include RPBE, OPBE, BLYP, B3LYP, B3LYP*, CAM-B3LYP, and LC-BLYP. 2. **Experimental Data Benchmark**: Using time-resolved X-ray emission spectroscopy (XES) data as a reference, the dynamics between metal-to-ligand charge transfer (MLCT) and metal-centered (MC) states of two prototype iron (II) complexes [Fe(bmip)₂]²⁺ and [Fe(terpy)₂]²⁺ were simulated. 3. **Method Comparison**: By comparing the results obtained from different methods, it was found that the B3LYP* method best reproduces the experimentally observed dynamic behavior, indicating its optimal balance in describing MLCT-MC energetics. 4. **Influence Factor Analysis**: The study found that the proportion of exact exchange (Hartree-Fock exchange) included is a key factor in determining the simulation results. Pure GGA methods (such as RPBE, OPBE, BLYP) do not include any exact exchange, leading to an over-stabilization of the MLCT state, thus failing to correctly simulate the experimentally observed phenomena. In summary, through detailed simulations and comparisons with experimental data, the paper validates the strengths and weaknesses of different DFT/TD-DFT methods in describing the ultrafast excited-state dynamics of iron (II) complexes, and identifies the B3LYP* method as the best performer in this field.