Relativistic Polarizable Embedding

A computational study of the aqueous pertechnetate anion: Elucidation of the hydration structure and spectroscopic properties

A computational study of the TcO4- anion in water has been carried out, the study includes an analysis of the structure of the solvation shells, the hydration energy and the electronic transitions that affect the shape of the UV-Vis absorption spectrum. The solvated system was characterized by combined QM/MM molecular dynamics simulations. The second-order perturbation theory restricted active space method and the novel self-consistent field density matrix renormalization group approach were employed to describe the static correlation. Two distinct solvation shells were found with 24 and 65 water molecules, respectively. The water molecules surrounding the anion have a general orientation in which one hydrogen is oriented away from the oxyanion while the remaining atoms are at a similar distance. However, there are specific water molecules that form hydrogen bonds with the oxygens of the oxyanion, while others are oriented with their oxygen atom towards the anion. The observed excitations in the UV-Vis spectrum are of a T2 character, with the main source of the observed behavior being the charge transfer from oxygen atoms to the central technetium atom. The calculations show that the most intense band of the spectrum is broader and has a blue tail with respect to the gas phase spectrum. This difference is due to the lower symmetry caused by the aqueous environment, which allows different states to mix and leads to broadening of the band.

Exact Two-Component Relativistic Polarizable Density Embedding

We have implemented the fragment-based polarizable density embedding (PDE) model within a relativistic framework building on the eXact 2-Component (X2C) relativistic Hamiltonian, thereby taking the PDE method to a relativistic framework. The PDE model provides a robust solution to the electron-leakage problem, and we show that this newly implemented model offers an accurate way to model solvated systems possessing significant relativistic effects. To demonstrate the model's performance, we perform comparative calculations of the K- and L2,3-edge spectra of water-solvated cysteine (both protonated and deprotonated) with the X2C Hamiltonian. Particularly, with counterions such as Na+ in the solvent, electron leakage clearly shows in the older polarizable embedding model through spurious peaks in the spectra. However, when the PDE model is employed, these spurious peaks disappear.

Accuracy of One- and Two-Photon Intensities with the Extended Polarizable Density Embedding Model

The recently developed extended polarizable density embedding (PDE-X) model is evaluated for the spectroscopic properties of organic chromophores solvated in water, including both one- and two-photon absorption properties. The PDE-X embedding model systematically improves vertical excitation energies over the preceding polarizable density embedding model (PDE). PDE-X shows more modest improvements over existing embedding models for oscillator strengths and two-photon absorption cross-sections, which are more sensitive properties. We argue that the origin of these discrepancies is related to the description of polarization effects, suggesting directions for future development of the embedding model.

A method to capture the large relativistic and solvent effects on the UV-vis spectra of photo-activated metal complexes

We have recently developed a method based on relativistic time-dependent density functional theory (TD-DFT) that allows the calculation of electronic spectra in solution (Creutzberg, Hedegård, J. Chem. Theory Comput.18, 2022, 3671). This method treats the solvent explicitly with a classical, polarizable embedding (PE) description. Furthermore, it employs the complex polarization propagator (CPP) formalism which allows calculations on complexes with a dense population of electronic states (such complexes are known to be problematic for conventional TD-DFT). Here, we employ this method to investigate both the dynamic and electronic effects of the solvent for the excited electronic states of trans-trans-trans-[Pt(N₃)₂(OH)₂(NH₃)₂] in aqueous solution. This complex decomposes into species harmful to cancer cells under light irradiation. Thus, understanding its photo-physical properties may lead to a more efficient method to battle cancer. We quantify the effect of the underlying structure and dynamics by classical molecular mechanics simulations, refined with a subsequent DFT or semi-empirical optimization on a cluster. Moreover, we quantify the effect of employing different methods to set up the solvated system, e.g., how sensitive the results are to the method used for the refinement, and how large a solvent shell that is required. The electronic solvent effect is always included through a PE potential.

Solvent effects in four-component relativistic electronic structure theory based on the reference interaction-site model

A combined method of the Dirac-Hartree-Fock (DHF) method and the reference interaction-site model (RISM) theory is reported; this is the initial implementation of the coupling of the four-component relativistic electronic structure theory and an integral equation theory of molecular liquids. In the method, the DHF and RISM equations are solved self-consistently, and therefore the electronic structure of the solute, including relativistic effects, and the solvation structure are determined simultaneously. The formulation is constructed based on the variational principle with respect to the Helmholtz energy, and analytic free energy gradients are also derived using the variational property. The method is applied to the iodine ion (I^- ), methyl iodide (CH₃ I), and hydrogen chalcogenide (H₂ X, where X = O-Po) in aqueous solutions, and the electronic structures of the solutes, as well as the solvation free energies and their component analysis, solvent distributions, and solute-solvent interactions, are discussed.

Polarizable Embedding Complex Polarization Propagator in Four- and Two-Component Frameworks

Explicit embedding methods combined with the complex polarization propagator (CPP) enable the modeling of spectroscopy for increasingly complex systems with a high density of states. We present the first derivation and implementation of the CPP in four- and exact-two-component (X2C) polarizable embedding (PE) frameworks. We denote the developed methods PE-4c-CPP and PE-X2C-CPP, respectively. We illustrate the methods by estimating the solvent effect on ultraviolet-visible (UV-vis) and X-ray atomic absorption (XAS) spectra of [Rh(H₂O)₆]³⁺ and [Ir(H₂O)₆]³⁺ immersed in aqueous solution. We moreover estimate solvent effects on UV-vis spectra of a platinum complex that can be photochemically activated (in water) to kill cancer cells. Our results clearly show that the inclusion of the environment is required: UV-vis and (to a lesser degree) XAS spectra can become qualitatively different from vacuum calculations. Comparison of PE-4c-CPP and PE-X2C-CPP methods shows that X2C essentially reproduces the solvent effect obtained with the 4c methods.

Investigating the influence of relativistic effects on absorption spectra for platinum complexes with light-activated activity against cancer cells

We report the first systematic investigation of relativistic effects on the UV-vis spectra of two prototype complexes for so-called photo-activated chemotherapy (PACT), trans-trans-trans-[Pt(N3)2(OH)2(NH3)2] and cis-trans-cis-[Pt(N3)2(OH)2(NH3)2]. In PACT, design of new drugs requires in-depth understanding of the photo-activation mechanisms. A first step is usually to rationalize their UV-vis spectra for which time-dependent density functional theory (TD-DFT) is an indispensable tool. We carried out TD-DFT calculations with a systematic series of non-relativistic (NR), scalar-relativistic (SR), and four-component (4c) Hamiltonians. As expected, large differences are found between spectra calculated within 4c and NR frameworks, while the most intense features (found at higher energies below 300 nm) can be reasonably well reproduced within a SR framework. It is also shown that effective core potentials (ECPs) yield essentially similar results as all-electron SR calculations. Yet the underlying transitions can be strongly influenced by spin-orbit coupling, which is only present in the 4c framework: while this can affect both intense and less intense transitions in the spectra, the effect is most pronounced for weaker transitions at lower energies, above 300 nm. Since the investigated complexes are activated with light of wavelengths above 300 nm, employing a method with explicit inclusion of spin-orbit coupling may be crucial to rationalize the activation mechanism.

The DIRAC code for relativistic molecular calculations

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.

CPPE: An Open-Source C++ and Python Library for Polarizable Embedding

We present a modular open-source library for polarizable embedding (PE) named CPPE. The library is implemented in C++, and it additionally provides a Python interface for rapid prototyping and experimentation in a high-level scripting language. Our library integrates seamlessly with existing quantum chemical program packages through an intuitive and minimal interface. Until now, CPPE has been interfaced to three packages, Q-Chem, Psi4, and PySCF. Furthermore, we show CPPE in action using all three program packages for a computational spectroscopy application. With CPPE, host program interfaces only require minor programming effort, paving the way for new combined methodologies and broader availability of the PE model.