RASPT2/RASSCF vs Range-Separated/Hybrid DFT Methods: Assessing the Excited States of a Ru(II)bipyridyl Complex

Computational Exploration of the Mechanism of Action of a Sorafenib-Containing Ruthenium Complex as an Anticancer Agent for Photoactivated Chemotherapy

Ruthenium(II) polypyridyl complexes are being tested as potential anticancer agents in different therapies, which include conventional chemotherapy and light-activated approaches. A mechanistic study on a recently synthesized dual-action Ru(II) complex [Ru(bpy)2(sora)Cl]+ is described here. It is characterized by two mono-dentate leaving ligands, namely, chloride and sorafenib ligands, which make it possible to form a di-aquo complex able to bind DNA. At the same time, while the released sorafenib can induce ferroptosis, the complex is also able to act as a photosensitizer according to type II photodynamic therapy processes, thus generating one of the most harmful cytotoxic species, 1O2. In order to clarify the mechanism of action of the drug, computational strategies based on density functional theory are exploited. The photophysical properties of the complex, which include the absorption spectrum, the kinetics of ISC, and the character of all the excited states potentially involved in 1O2 generation, as well as the pathway providing the di-aquo complex, are fully explored. Interestingly, the outcomes show that light is needed to form the mono-aquo complex, after releasing both chloride and sorafenib ligands, while the second solvent molecule enters the coordination sphere of the metal once the system has come back to the ground-state potential energy surface. In order to simulate the interaction with canonical DNA, the di-aquo complex interaction with a guanine nucleobase as a model has also been studied. The whole study aims to elucidate the intricate details of the photodissociation process, which could help with designing tailored metal complexes as potential anticancer agents.

Benchmarks for transition metal spin-state energetics: why and how to employ experimental reference data?

Accurate prediction of energy differences between alternative spin states of transition metal complexes is essential in computational (bio)inorganic chemistry-for example, in characterization of spin crossover materials and in the theoretical modeling of open-shell reaction mechanisms-but it remains one of the most compelling problems for quantum chemistry methods. A part of this challenge is to obtain reliable reference data for benchmark studies, as even the highest-level applicable methods are known to give divergent results. This Perspective discusses two possible approaches to method benchmarking for spin-state energetics: using either theoretically computed or experiment-derived reference data. With the focus on the latter approach, an extensive general review is provided for the available experimental data of spin-state energetics and their interpretations in the context of benchmark studies, targeting the possibility of back-correcting the vibrational effects and the influence of solvents or crystalline environments. With a growing amount of experience, these effects can be now not only qualitatively understood, but also quantitatively modeled, providing the way to derive nearly chemically accurate estimates of the electronic spin-state gaps to be used as benchmarks and advancing our understanding of the phenomena related to spin states in condensed phases.

Computational Assessment of a Dual-Action Ru(II)-Based Complex: Photosensitizer in Photodynamic Therapy and Intercalating Agent for Inducing DNA Damage

A combined quantum-mechanical and classical molecular dynamics study of a recent Ru(II) complex with potential dual anticancer action is reported here. The main basis for the multiple action relies on the merocyanine ligand, whose electronic structure allows the drug to be able to absorb within the therapeutic window and in turn efficiently generate ¹O₂ for photodynamic therapy application and to intercalate within two nucleobases couples establishing reversible electrostatic interactions with DNA. TDDFT outcomes, which include the absorption spectrum, triplet states energy, and spin-orbit matrix elements, evidence that the photosensitizing activity is ensured by an MLCT state at around 660 nm, involving the merocyanine-based ligand, and by an efficient ISC from such state to triplet states with different characters. On the other hand, the MD exploration of all the possible intercalation sites within the dodecamer B-DNA evidences the ability of the complex to establish several electrostatic interactions with the nucleobases, thus potentially inducing DNA damage, though the simulation of the absorption spectra for models extracted by each MD trajectory shows that the photosensitizing properties of the complex remain unaltered. The computational results support that the anti-tumor effect may be related to multiple mechanisms of action.

Theoretical investigation of [Ru(bpy)2(HAT)]2+ (HAT = 1,4,5,8,9,12-hexaazatriphenylene; bpy = 2,2'-bipyridine): Photophysics and reactions in excited state

In this article, Density Functional Theory based calculations, including dispersion corrections, PBE0(D3BJ)/Def2-TZVP(-f), were performed to elucidate the photophysics of the [Ru(bpy)₂(HAT)]²⁺ complex in water. In addition, the thermodynamics of the charge and electron transfer excited state reactions of this complex with oxygen, nitric oxide and Guanosine-5'-monophosphate nucleotide (GMP) were investigated. The first singlet excite state, S₁, strongly couples with the second and third triplet excited states (T₂ and T₃) giving rise to a high intersystem crossing rate of 6.26 × 10¹¹ s^-1 which is ∼10⁶ greater than the fluorescence rate decay. The thermodynamics of the excited reactions revealed that all electron transfer reactions investigated are highly favorable, due mainly to the high stability of the triply charged radical cation ²PS^•3+ species formed after the electron has been transferred. Excited state electron transfer from the GMP nucleotide to the complex is also highly favorable (ΔG_sol = -92.6 kcal/mol), showing that this complex can be involved in the photooxidation of DNA, in line with experimental findings. Therefore, the calculations allow to conclude that the [Ru(bpy)₂(HAT)]²⁺ complex can act in Photodynamic therapy through both mechanisms type I and II, through electron transfer from and to the complex and triplet-triplet energy transfer, generating ROS, RNOS and through DNA photooxidation. In addition, the work also opens a perspective of using this complex for the in-situ generation of the singlet nitroxyl (¹NO^-) species, which can have important applications for the generation of HNO and may have, therefore, important impact for physiological studies involving HNO.

Unraveling the Symmetry Effects on the Second-Order Nonlinear Optical Responses of Molecular Switches: The Case of Ruthenium Complexes

Owing to their odd order, second-order nonlinear optical (NLO) responses are very sensitive to symmetry. Therefore, within hyper-Rayleigh scattering (HRS) technique, the symmetry impacts the amplitude of the molecular responses, the HRS first hyperpolarizability (β_HRS), and the depolarization ratio (DR). Starting from a challenging octupolar structure bearing six ruthenium(II) ammine centers π-conjugated via quaterpyridyl moieties to a tris-chelated zinc(II) core, together with its Λ shape and one-dimensional analogues built by replacing one or two Ru-quaterpyridyl moieties with bipyridine moieties, (time-dependent) density functional theory calculations have been enacted to unravel the symmetry-NLO response relationships as well as their Ru^II/III redox-triggered switching effects. The one-dimensional and Λ-shaped NLOphores present β_HRS values ∼3 times larger than those of the octupolar system, for both Ru oxidation states. However, using the few-state valence bond-charge transfer models demonstrates that the β_HRS response of the octupolar complex cation can become larger than those of its one-dimensional and Λ-shaped analogues provided stronger donor-acceptor groups are employed. In parallel, the DRs decrease from a strong dipolar character (DR ≈ 6) for the one-dimensional chromophore to a weaker dipolar character (DR ≈ 5) for the Λ-shaped one and to a clear octupolar character (DR ≈ 1.7) for the last one. In all cases, the β responses originate mostly from metal-to-ligand charge transfer excited states, as revealed using a new scheme for analyzing the variations in electron density upon excitation. The Ru^II/III oxidations lead to a strong decrease in the β_HRS responses, which is attributed to the loss of the donor character of the Ru centers and therefore to the reduction of the push-pull π-conjugation. These results demonstrate that the NLO contrast and the NLO switching behavior of these Ru cations are maintained for the different molecular symmetries. Finally, the character of the β responses of the oxidized species, as revealed by the DR values, further evidences a clear evolution from dipolar to octupolar NLOphores.

Metal-to-Ligand Charge-Transfer Spectrum of a Ru-Bipyridine-Sensitized TiO2 Cluster from Embedded Multiconfigurational Excited-State Theory

Understanding optical properties of the dye molecule in dye-sensitized solar cells (DSSCs) from first-principles quantum mechanics can contribute to improving the efficiency of such devices. While density functional theory (DFT) and time-dependent DFT have been pivotal in simulating optoelectronic properties of photoanodes used in DSSCs at the atomic scale, questions remain regarding DFT's adequacy and accuracy to furnish critical information needed to understand the various excited-state processes involved. Here, we simulate the absorption spectra of a dye-sensitized solar cell analogue, comprised of a Ru-bipyridine (Ru-bpy) dye molecule and a small TiO₂ cluster via DFT and via an accurate embedded correlated wavefunction (CW) theory. We generated CW spectra for the adsorbed Ru-bpy dye via a recently introduced capped density functional embedding theory or capped-DFET (to generate the embedding potential that accounts for the interaction of the molecule and the TiO₂ cluster). We then combined capped-DFET with the accurate but expensive multiconfigurational complete active space second-order perturbation theory (CASPT2)-embedded CASPT2. Because the CW theory is conducted on only a portion of the total system in the presence of an embedding potential that describes that portion's interaction with its environment, we efficiently obtain CW-quality predictions that reflect local properties of the entire system. Specifically, for example, with capped-DFET and embedded CW theory, we can simulate accurately a plethora of metal-to-ligand charge-transfer excited properties at a manageable computational cost. Here, we predict detailed electronic spectra within the visible region, featuring the lowest three singlet and triplet excited states, along with predictions of the singlets' lifetimes. We illustrated these results using a Jablonski diagram that show the relative energy position of the singlet and longer-lived triplet excited states and analyzed and proposed relaxation paths for the excited state corresponding to the most intense but short-lived absorption (interconversion, intersystem crossing, fluorescence, and phosphorescence) that may lead to longer-lived excited states necessary for efficient charge separation required to generate current in solar cells.

The Electronic Structures of CoGe n-/0 ( n = 1-3) Clusters from Multiconfigurational CASSCF/CASPT2 and RASSCF/RASPT2 Calculations

Density functional theory and multiconfigurational CASPT2 and RASPT2 methods are employed to investigate the low-lying states of CoGe _n^-/0 ( n = 1-3) clusters. With the RASPT2 approach, the active space is extended to 14 orbitals for CoGe^-/0, 17 orbitals for CoGe₂^-/0, and 20 orbitals for CoGe₃^-/0. These active spaces include the 3d, 4s, and 4d of Co and 4p of Ge. The 4d of Co is incorporated into these active spaces in order to account for the important double-shell effect of Co. The structural parameters, vibrational frequencies, and relative energies of the low-lying states of CoGe _n^-/0 ( n = 1-3) are reported. The ground states of CoGe _n^- ( n = 1-3) are computed to be ³Φ of linear CoGe^-, ³B₁ of cyclic CoGe₂^-, and ³B₁ of cyclic CoGe₃^- isomer. The ground states of the neutral clusters are calculated to be ²Δ of linear CoGe, ⁴B₁ of cyclic CoGe₂, and ⁴A″ of tetrahedral CoGe₃ isomer. The calculated adiabatic and vertical detachment energies of the anionic ground states are in agreement with the experimental values as observed in the 266 nm anion photoelectron spectra.

Excited-State Proton Transfer Controls Irreversibility of Photoisomerization in Mononuclear Ruthenium(II) Monoaquo Complexes: A DFT Study

The detailed DFT investigation clears the working mechanism of the irreversible photoisomerization of trans-[Ru(tpy)(pynp)(OH2)](2+) (TA) and cis-[Ru(tpy)(pynp)(OH2)](2+) (CA) complexes. Both TA and CA complexes present two types of low lying triplet states, one resulting from a triplet metal-ligand charge-transfer (TMLCT) and the other from a triplet metal-centered d-d transition (TMC). The vertical excitation of the singlet ground state of the complexes leads to a singlet excited state, which undergoes ultrafast decay to the corresponding TMLCT. For TA, this TMLCT transforms with a low barrier to a TMC state. The dissociative nature of the TMC state leads to easy water removal to produce a five-coordinate intermediate that can isomerize via rotation of a pynp ligand and proceed towards the CA product. For CA, however, during this excitation and intersystem crossing process, an excited-state proton transfer (ESPT) occurs and the resultant TMLCT is very much stabilized with a very strong Ru(II)-OH bond; the high barrier from this TMLCT blocks conversion to a TMC state and thus prevents isomerization from the cis to the trans isomer. This high barrier also prevents the possibility of the isomerization process from TA to CA solely on the adiabatic triplet pathway. Instead, crossing points (XMC-CB, XMC-CA) near the minimum of the triplet metal-centered state of the cis isomer provide nonadiabatic decay channels to the ground-state S0--CA, which completes the photoisomerization pathway from TA to CA.

Spectroscopic and second-order nonlinear optical properties of Ruthenium(ii) complexes: a DFT/MRCI and ADC(2) study

We present an assessment of correlated electronic structure methods for the nonlinear optical properties of Ru(ii) dyes.

Charge-transfer excited states in phosphorescent organo-transition metal compounds: a difficult case for time dependent density functional theory?

A series of 17 platinum(ii) and iridium(iii) complexes have been investigated theoretically and experimentally to elucidate the charge-transfer character in emission from the lowest triplet state. TDDFT is found to be surprisingly accurate.

Computational insights into the photodeactivation dynamics of phosphors for OLEDs: a perspective

In this perspective we highlight recent computational efforts to unravel competing photodeactivation mechanisms of radiative and non-radiative nature of phosphors.

Pyridine imines as ligands in luminescent iridium complexes

Biscyclometallated iridium complexes [Ir(ppz)₂(X^Y)][PF₆] (X^Y = pyridine imine) have been synthesised; N-alkyl complexes are emissive whilst N-aryl ones are only weakly emissive, DFT calculations shed light on the reasons for this difference.

Theoretical modelling of photoswitching of hyperpolarisabilities in ruthenium complexes

Static excited-state polarisabilities and hyperpolarisabilities of three Ru^II ammine complexes are computed at the density functional theory (DFT) and several correlated ab initio levels. Most accurate modelling of the low energy electronic absorption spectrum is obtained with the hybrid functionals B3LYP, B3P86 or M06 for the complex [Ru^II(NH₃)₅(MeQ⁺)]³⁺ (MeQ⁺=N-methyl-4,4′-bipyridinium, 3) in acetonitrile. The match with experimental data is less good for [Ru^II(NH₃)₅L]³⁺ (L=N-methylpyrazinium, 2; N-methyl-4-{E,E-4-(4-pyridyl)buta-1,3-dienyl}pyridinium, 4). These calculations confirm that the first dipole- allowed excited state (FDAES) has metal-to-ligand charge-transfer (MLCT) character. Both the solution and gas-phase results obtained for 3 by using B3LYP, B3P86 or M06 are very similar to those from restricted active-space SCF second-order perturbation theory (RASPT2) with a very large basis set and large active space. However, the time-dependent DFT λ_max predictions from the long-range corrected functionals CAM-B3LYP, LC-ωPBE and wB97XB and also the fully ab initio resolution of identity approximate coupled-cluster method (gas-phase only) are less accurate for all three complexes. The ground state (GS) two-state approximation first hyperpolarisability β_2SA for 3 from RASPT2 is very close to that derived experimentally via hyper-Rayleigh scattering, whereas the corresponding DFT-based values are considerably larger. The β responses calculated by using B3LYP, B3P86 or M06 increase markedly as the π-conjugation extends on moving along the series 2→4, for both the GS and FDAES species. All three functionals predict substantial FDAES β enhancements for each complex, increasing with the π-conjugation, up to about sevenfold for 4. Also, the computed second hyperpolarisabilities γ generally increase in the FDAES, but the results vary between the different functionals.

Theoretical Investigation of the Electronic Structure of Fe(II) Complexes at Spin-State Transitions

The electronic structure relevant to low spin (LS)↔high spin (HS) transitions in Fe(II) coordination compounds with a FeN₆ core are studied. The selected [Fe(tz)₆]²⁺ (1) (tz = 1H-tetrazole), [Fe(bipy)₃]²⁺ (2) (bipy = 2,2′-bipyridine), and [Fe(terpy)₂]²⁺ (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 (ΔE_HL) 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 d⁶ 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 ΔE_HL. 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)₂]²⁺ 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)₂]²⁺ cannot be described using a single configuration coordinate.