1
|
Hansen BL, Markmann V, Pápai M, Lenzen P, Haubro ML, Mikeházi A, Németh Z, Vancza A, Levantino M, Zerdane S, Zederkof DB, Khakhulin D, Dohn AO, Nielsen MM, Haldrup K, Vankó G. Excited-state structural characterization of a series of nanosecond-lived [Fe(terpy)2]2+ derivatives using x-ray solution scattering. J Chem Phys 2025; 162:124308. [PMID: 40162540 DOI: 10.1063/5.0237705] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2024] [Accepted: 02/20/2025] [Indexed: 04/02/2025] Open
Abstract
[ F e ( t e r p y ) 2 ] 2 + (terpy = 2,2':6',2″-terpyridine) is a transition metal complex where the spin state is photoswitchable and where the properties of the metal-centered quintet excited state (5MC) can be tuned by substituting different electron withdrawing or electron donating groups on the 4' position of the terpyridine. To better understand the physics determining the photoswitching performance, a deeper insight into the positions of the relevant potential energy surfaces and the molecular structure of the 5MC state is needed. We present a structural investigation based on Time Resolved x-ray Solution Scattering (TR-XSS) by which we determine the average dFe-N bond-length elongation following population of the 5MC state as well as the lifetime of this state in a series of seven modified [Fe(terpy)2]2+ systems in aqueous solution following photo-excitation. The analysis of the TR-XSS data is supported by Density Functional Theory (DFT) and Molecular Dynamics calculations. The quintet state lifetime is determined to vary by more than a factor of 10 (from 1.5 to 16 ns) based on the electron withdrawing/donating properties of the substituting group. Both the DFT calculations and the structural analysis of the experimental data show that the main photo-induced change in metal-ligand bond lengths ΔdFe-N is ∼0.2 Å for all systems.
Collapse
Affiliation(s)
- Bianca L Hansen
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Verena Markmann
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Mátyás Pápai
- HUN-REN Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| | - Philipp Lenzen
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Morten Lunn Haubro
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Antal Mikeházi
- HUN-REN Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| | - Zoltán Németh
- HUN-REN Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| | - Andor Vancza
- HUN-REN Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| | - Matteo Levantino
- European Synchrotron Radiation Facility (ESRF), CS40220, 38043 Grenoble, Cedex 9, France
| | - Serhane Zerdane
- European Synchrotron Radiation Facility (ESRF), CS40220, 38043 Grenoble, Cedex 9, France
| | | | | | - Asmus Ougaard Dohn
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Martin Meedom Nielsen
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - Kristoffer Haldrup
- Department of Physics, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| | - György Vankó
- HUN-REN Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
2
|
Pápai M. Simulation of Ultrafast Excited-State Dynamics in Fe(II) Complexes: Assessment of Electronic Structure Descriptions. J Chem Theory Comput 2025; 21:560-574. [PMID: 39752586 PMCID: PMC11780750 DOI: 10.1021/acs.jctc.4c01331] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2024] [Revised: 12/02/2024] [Accepted: 12/09/2024] [Indexed: 01/29/2025]
Abstract
The assessment of electronic structure descriptions utilized in the simulation of the ultrafast excited-state dynamics of Fe(II) complexes is presented. Herein, we evaluate the performance of the RPBE, OPBE, BLYP, B3LYP, B3LYP*, PBE0, TPSSh, 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 multiplicities (singlet, triplet, and quintet) highly depend on the utilized DFT/TD-DFT method, with the percentage of exact (Hartree-Fock) exchange being the governing factor. Importantly, B3LYP* and TPSSh are the only DFT/TD-DFT methods with satisfactory performance, best reproducing the experimentally resolved dynamics for both complexes, signaling 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.
Collapse
Affiliation(s)
- Mátyás Pápai
- HUN-REN Wigner Research
Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
3
|
Branching mechanism of photoswitching in an Fe(II) polypyridyl complex explained by full singlet-triplet-quintet dynamics. Commun Chem 2023; 6:7. [PMID: 36697805 PMCID: PMC9829715 DOI: 10.1038/s42004-022-00796-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Accepted: 12/09/2022] [Indexed: 01/11/2023] Open
Abstract
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 interpretation of the photorelaxation remains controversial. 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 involving two sequential processes: a dominant 3MLCT→3MC(3T2g)→3MC(3T1g)→5MC, and a minor 3MLCT→3MC(3T2g)→5MC component. (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) polypyridines, and fully describe the photorelaxation dynamics.
Collapse
|
4
|
Pápai M. Toward Simulation of Fe(II) Low-Spin → High-Spin Photoswitching by Synergistic Spin-Vibronic Dynamics. J Chem Theory Comput 2022; 18:1329-1339. [PMID: 35199532 PMCID: PMC8908767 DOI: 10.1021/acs.jctc.1c01184] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
A new theoretical
approach is presented and applied for the simulation
of Fe(II) low-spin (LS, singlet, t2g6eg0) → high-spin (HS, quintet, t2g4eg2) photoswitching dynamics of the octahedral
model complex [Fe(NCH)6]2+. The utilized synergistic
methodology heavily exploits the strengths of complementary electronic
structure and spin-vibronic dynamics methods. Specifically, we perform
3D quantum dynamics (QD) and full-dimensional trajectory surface hopping
(TSH, in conjunction with a linear vibronic coupling model), with
the modes for QD selected by TSH. We follow a hybrid approach which
is based on the application of time-dependent density functional theory
(TD-DFT) excited-state potential energy surfaces (PESs) and multiconfigurational
second-order perturbation theory (CASPT2) spin–orbit couplings
(SOCs). Our method delivers accurate singlet–triplet–quintet
intersystem crossing (ISC) dynamics, as assessed by comparison to
our recent high-level ab initio simulations and related
time-resolved experimental data. Furthermore, we investigate the capability
of our simulations to identify the location of ISCs. Finally, we assess
the approximation of constant SOCs (calculated at the Franck–Condon
geometry), whose validity has central importance for the combination
of TD-DFT PESs and CASPT2 SOCs. This efficient methodology will have
a key role in simulating LS → HS dynamics for more complicated
cases, involving higher density of states and varying electronic character,
as well as the analysis of ultrafast experiments.
Collapse
Affiliation(s)
- Mátyás Pápai
- Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|