Search

Core-electron excitations in solvated systems, influenced by solvent geometry and hydrogen bonding, make X-ray absorption spectroscopy (XAS) a valuable tool for assessing solvent-solute interactions. However, calculating XAS spectra with electronic-structure methods has proven challenging due to a delicate interplay between correlation and solvation effects. This study provides a computational procedure for XAS modeling in solvated systems, with water-solvated ammonia and ammonium systems serving as probes. Exploring methodological challenges, we investigate explicit embedding models, specifically the polarizable embedding family, including polarizable density embedding and extended polarizable density embedding. Our linear-response time-dependent density functional theory (LR-TDDFT) XAS calculations reveal the efficiency of this approach, with extended polarizable density embedding emerging as a robust improvement over polarizable density embedding. Contrary to some recent literature, our study challenges the belief that LR-TDDFT cannot accurately describe XAS spectra of ammonia and ammonium solvated in water.

We extend the polarizable density embedding (PDE) model to support the calculation of nuclear magnetic resonance (NMR) shielding constants using gauge-including atomic orbitals (GIAOs) within a density functional theory (DFT) framework. The PDE model divides the total system into fragments, describing some by quantum mechanics (QM) and the others through an embedding model. The PDE model uses anisotropic polarizabilities, inter-fragment two-electron Coulomb integrals, and a non-local repulsion operator to emulate the QM effects. The terms involving Coulomb integrals are straightforwardly extended with GIAOs. In contrast, we consider two approaches to handle the gauge dependency of the non-local operator, employing either simple symmetrization or a gauge transformation. We find the latter approach to be most stable with respect to increasing the basis set size of the QM region. We examine the accuracy of the PDE model for calculating NMR shielding constants on several solutes in a water solution. The performance is compared with the classical polarizable embedding (PE) model in addition to supermolecular reference calculations. Based on these systems, we address the basis set convergence characteristics and the QM region size requirements. Furthermore, we investigate the performance of the PDE model for a system with significant electron spill-out. In many cases, we find that the PDE model outperforms the PE model, especially regarding the accuracy of nuclear shielding constants when using small QM region sizes and in systems with significant electron spill-out.

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 2 O) 6 ] 3+ and [Ir(H 2 O) 6 ] 3+ 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.

Most chemistry, including chemistry where relativistic effects are important, occurs in an environment, and in many cases, this environment has a significant effect on the chemistry. In nonrelativistic quantum chemistry, a lot of progress has been achieved with respect to including environments such as a solvent or protein in the calculations, and now is the time to extend the possibilities for also doing this in relativistic quantum chemistry. The polarizable embedding (PE) model efficiently incorporates electrostatic effects of the environment by describing it as a collection of localized electric multipoles and polarizabilities obtained through quantum chemical calculations. In this article, we present the theory and implementation of four- and exact two-component Hamiltonians within a PE framework. We denote the methods the PE-4c-DFT and PE-X2C-DFT models. The models include a linear response formalism to calculate time-dependent (TD) properties: PE-TD-4c-DFT and PE-TD-X2C-DFT. With this first implementation, we calculate the PE-TD-4c-PBE0 excitation energies of the TcO 4 - and ReO 4 - ions in an explicit water solvent. This initial investigation focuses on the relative size of relativistic and solvent contributions to the excitation energies. The solvent effect is divided into an indirect solvent effect due to the structural perturbation of the XO 4 - ion and a direct electrostatic effect. The relativistic effects as well as both types of solvent effects are found to contribute to a shift in the excitation energies, but they do so to different extents depending on the ion and the electronic transition in question.

The polarizable embedding (PE) approach is a flexible embedding model where a preselected region out of a larger system is described quantum mechanically, while the interaction with the surrounding environment is modeled through an effective operator. This effective operator represents the environment by atom-centered multipoles and polarizabilities derived from quantum mechanical calculations on (fragments of) the environment. Thereby, the polarization of the environment is explicitly accounted for. Here, we present the coupling of the PE approach with the density matrix renormalization group (DMRG). This PE-DMRG method is particularly suitable for embedded subsystems that feature a dense manifold of frontier orbitals which requires large active spaces. Recovering such static electron-correlation effects in multiconfigurational electronic structure problems, while accounting for both electrostatics and polarization of a surrounding environment, allows us to describe strongly correlated electronic structures in complex molecular environments. We investigate various embedding potentials for the well-studied first excited state of water with active spaces that correspond to a full configuration-interaction treatment. Moreover, we study the environment effect on the first excited state of a retinylidene Schiff base within a channelrhodopsin protein. For this system, we also investigate the effect of dynamical correlation included through short-range density functional theory.

Computational methods that can accurately and effectively predict all types of electronic excitations for any molecular system are missing in the toolbox of the computational chemist. Although various Kohn-Sham density-functional methods (KS-DFT) fulfill this aim in some cases, they become inadequate when the molecule has near-degeneracies and/or low-lying double-excited states. To address these issues we have recently proposed multiconfiguration short-range density-functional theory-MC-srDFT-as a new tool in the toolbox. While initial applications for systems with multireference character and double excitations have been promising, it is nevertheless important that the accuracy of MC-srDFT is at least comparable to the best KS-DFT methods also for organic molecules that are typically of single-reference character. In this paper we therefore systematically investigate the performance of MC-srDFT for a selected benchmark set of electronic excitations of organic molecules, covering the most common types of organic chromophores. This investigation confirms the expectation that the MC-srDFT method is accurate for a broad range of excitations and comparable to accurate wave function methods such as CASPT2, NEVPT2, and the coupled cluster based CC2 and CC3.

We present a combination of the polarizable embedding (PE) scheme with the complex polarization propagator (CPP) method with the aim of calculating response properties including relaxation for large and complex systems. This new approach, termed PE-CPP, will benefit from the highly advanced description of the environmental electrostatic potential and polarization in the PE method as well as the treatment of near-resonant effects in the CPP approach. The PE-CPP model has been implemented in a Kohn-Sham density functional theory approach, and we present pilot calculations exemplifying the implementation for the UV/vis and carbon K-edge X-ray absorption spectra of the protein plastocyanin. Furthermore, technical details associated with a PE-CPP calculation are discussed.

Even for pure substances, the deduction of the absolute configuration is not always straightforward since there is no direct link between the magnitude and sign of the optical rotation and the absolute configuration. It would be very useful to use computations of the optical rotation to link experimentally measured optical rotations to an absolute configuration. Such electronic structure calculations of the optical rotation typically employ regular energy optimized basis sets from wave function theory, and especially the aug-cc-pVDZ basis set has been popular. Here, we have carried out extrapolation of the optical rotation to the basis set limits for nine small or medium sized molecules, using basis sets developed specifically for DFT and magnetic properties (aug-pcS-n series). We suggest that assignment of absolute configuration by comparisons between theoretical and experimental optical rotations may be improved by employing different wavelengths, and accordingly the optical rotation at two wavelengths (589.3 and 355.0 nm) has been investigated. Several fitting schemes were used to estimate the optical rotations at the basis set limit. It was found that use of the aug-cc-pVDZ basis set often leads to results that deviate significantly form the basis set limit results, especially at 355.0 nm but also at 589.3 nm. The double-ζ aug-pcS-1 basis set usually provides results which are closer to the limiting values. The basis set requirements are generally more severe at 355.0 nm, where also the aug-cc-pVTZ and 6-311++G(3fd,3dp) basis sets show significant deviations from the basis set limit results, while the aug-pcS-2 basis set always leads to results within an acceptable deviation.

The hyperfine coupling tensor of electron paramagnetic resonance (EPR), describing the interaction between an electron and a given nuclei, depends strongly on the electron density at the nucleus. With standard Gaussian-type orbital basis sets (GTOs), employed in most calculations, it is difficult to obtain converged results of the hyperfine coupling tensor, and basis sets with more flexible core regions have therefore been devised. To this class of core property basis sets belong the aug-cc-pVTZ-J basis sets developed for the s- and p-block atoms. Here, we extend the aug-cc-pVTZ-J basis sets to include the 3d elements Sc-Zn. The converged optimal basis sets are throughout the series described by a (25s17p10d3f2g)/[17s10p7d3f2g] contraction scheme, where four tight s-, one tight p-, and one tight d-type function have been added to the original aug-cc-pVTZ basis sets. The basis sets are generally contracted, and molecular orbital coefficients are used as contraction coefficients. By validation studies with different functionals and compounds, it is shown that the values of the contraction coefficient are effectively independent of the compound used in their generation and the exchange-correlation functional employed in the calculation.