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.

Decomposition of molecular properties

We review recent work on property decomposition techniques using quantum chemical methods and discuss some topical applications in terms of quantum mechanics-molecular mechanics calculations and the constructing of properties of large molecules and clusters. Starting out from the so-called LoProp decomposition scheme [Gagliardi et al., J. Chem. Phys., 2004, 121, 4994] for extracting atomic and inter-atomic contributions to molecular properties we show how this method can be generalized to localized frequency-dependent polarizabilities, to localized hyperpolarizabilities and to localized dispersion coefficients. Some applications of the generalized decomposition technique are reviewed - calculations of frequency-dependent polarizabilities, Rayleigh scattering of large clusters, and calculations of hyperpolarizabilities of proteins.

Nonequilibrium Environment Dynamics in a Frequency-Dependent Polarizable Embedding Model

Hybrid quantum mechanical/molecular mechanical (QM/MM) models are some of the most powerful and computationally feasible approaches to account for solvent effects or more general environmental perturbations on quantum chemical systems. In their more recent formulations (known as polarizable embedding) they can account for electrostatic and mutual polarization effects between the QM and the MM subsystems. In this paper, a polarizable embedding scheme based on induced dipoles that is able both to describe electron evolution of the embedded QM system in an efficient manner as well as to capture the frequency dependent behavior of the solvent is proposed, namely, ωMMPol. The effects of this frequency-dependent solvent on a time-dependent model system-the Rabi oscillations of H₂⁺ in a resonant field-are considered. The solvent is shown to introduce only mild perturbations when the excitation frequencies of the solvent and the solute are off-resonant. However, the dynamics of the H₂⁺ are fundamentally changed in the presence of a near-resonant excitation solvent. The effectiveness of ωMMPol to simulating realistic chemical systems is demonstrated by capturing charge transfer dynamics within a solvated system.

A quantum-mechanical perspective on linear response theory within polarizable embedding

We present a derivation of linear response theory within polarizable embedding starting from a rigorous quantum-mechanical treatment of a composite system. To this aim, two different subsystem decompositions (symmetric and nonsymmetric) of the linear response function are introduced and the pole structures as well as residues of the individual terms are discussed. In addition to providing a thorough justification for the descriptions used in polarizable embedding models, this theoretical analysis clarifies which form of the response function to use and highlights complications in separating out subsystem contributions to molecular properties. The basic features of the presented expressions and various approximate forms are illustrated by their application to a composite model system.

Local decomposition of imaginary polarizabilities and dispersion coefficients

We present a new way to compute the two-body contribution to the dispersion energy using ab initio theory.

First Hyperpolarizability of Collagen Using the Point Dipole Approximation

The application of localized hyperpolarizabilities to predict a total protein hyperpolarizability is presented for the first time, using rat-tail collagen as a demonstration example. We employ a model comprising the quadratic Applequist point-dipole approach, the so-called LoProp transformation, and a procedure with molecular fractionation using conjugate caps to determine the atomic and bond contributions to the net β tensor of the collagen [(PPG)10]3 triple-helix. By using Tholes exponential damping modification to the dyadic tensor in the Applequist equations, a correct qualitative agreement with experiment is found. The intensity of the βHRS signal and the depolarization ratios are best reproduced by decomposing the LoProp properties into the atomic positions and using Tholes exponential damping with the original damping parameter. Some ramifications of the model for general protein property optimization are briefly discussed.

Hyperpolarizabilities of extended molecular mechanical systems

The equations for decomposition of the molecular hyperpolarizabilty into atomic sites is derived alongside with equations to calculate the cluster hyperpolarizability as resulting from the second order induced dipole moments in the point-dipole model.

Excited states in large molecular systems through polarizable embedding

Using the polarizable embedding model enables rational design of light-sensitive functional biological materials.