Recent Advances in Wave Function-Based Methods of Molecular-Property Calculations

N-Electron Valence Perturbation Theory with Reference Wave Functions from Quantum Computing: Application to the Relative Stability of Hydroxide Anion and Hydroxyl Radical

Quantum simulations of the hydroxide anion and hydroxyl radical are reported, employing variational quantum algorithms for near-term quantum devices. The energy of each species is calculated along the dissociation curve, to obtain information about the stability of the molecular species being investigated. It is shown that simulations restricted to valence spaces incorrectly predict the hydroxyl radical to be more stable than the hydroxide anion. Inclusion of dynamical electron correlation from nonvalence orbitals is demonstrated, through the integration of the variational quantum eigensolver and quantum subspace expansion methods in the workflow of N-electron valence perturbation theory, and shown to correctly predict the hydroxide anion to be more stable than the hydroxyl radical, provided that basis sets with diffuse orbitals are also employed. Finally, we calculate the electron affinity of the hydroxyl radical using an aug-cc-pVQZ basis on IBM's quantum devices.

Relativistic Density Functional NMR Tensors Analyzed with Spin-free Localized Molecular Orbitals

The implementation of fast relativistic methods based on density functional theory, in conjunction with localized molecular orbital (LMO) based analysis, allows straightforward interpretations of NMR parameters in terms of contributions from core shells, lone pairs, and bonds, for compounds containing elements from across the periodic table. We present a conceptual review of a frequently used LMO analysis of NMR parameters calculated in the presence of spin-orbit interactions and other relativistic effects. An accompanying example focuses on the ¹⁵ N shielding in a heavy metal complex.

Real-time electronic energy current and quantum energy flux in molecules

Intra- and inter-molecular electronic energy current is formulated by defining the probability current of electronic energy, called the energy flux. Among vast possible applications to electronic energy transfer phenomena, including chemical reaction dynamics, here we present a first numerical example from highly excited nonadiabatic electron wavepacket dynamics of a boron cluster B₁₂.

Predicting the Second-Order Nonlinear Optical Responses of Organic Materials: The Role of Dynamics

The last 30 years have witnessed an ever-growing application of computational chemistry for rationalizing the nonlinear optical (NLO) responses of organic chromophores. More specifically, quantum chemical calculations proved highly helpful in gaining fundamental insights into the factors governing the magnitude and character of molecular first hyperpolarizabilities (β), be they either intrinsic to the chromophore molecular structure and arising from symmetry, chemical substitution, or π-electron delocalization, or induced by external contributions such as the laser probe or solvation and polarization effects. Most theoretical reports assumed a rigid picture of the investigated systems, the NLO responses being computed solely at the most stable geometry of the chromophores. Yet, recent developments combining classical molecular dynamics (MD) simulations and DFT calculations have evidenced the significant role of structural fluctuations, which may induce broad distributions of NLO responses, and even generate them in some instances.This Account presents recent case studies in which theoretical simulations have highlighted these effects. The discussion specifically focuses on the simulation of the second-order NLO properties that can be measured experimentally either from Hyper-Rayleigh Scattering (HRS) or Electric-Field Induced Second Harmonic Generation (EFISHG). More general but technical topics concerning several aspects of the calculations of hyperpolarizabilities are instead discussed in the Supporting Information.Selected examples include organic chromophores, photochromic systems, and ionic complexes in the liquid phase, for which the effects of explicit solvation, concentration, and chromophore aggregation are emphasized, as well as large flexible systems such as peptide chains and pyrimidine-based helical polymers, in which the relative variations of the responses were shown to be several times larger than their average values. The impact of geometrical fluctuations is also illustrated for supramolecular architectures with the examples of nanoparticles formed by organic dipolar dyes in water solution, whose soft nature allows for large shape variations translating into huge fluctuations in time of their NLO response, and of self-assembled monolayers (SAMs) based on indolino-oxazolidine or azobenzene switches, in which the geometrical distortions of the photochromic molecules, as well as their orientational and positional disorder within the SAMs, highly impact their NLO response and contrast upon switching. Finally, the effects of the rigidity and fluidity of the surrounding are evidenced for NLO dyes inserted in phospholipid bilayers.

Power of the Sine Hamiltonian Operator for Estimating the Eigenstate Energies on Quantum Computers

Quantum computers have been shown to have tremendous potential in solving difficult problems in quantum chemistry. In this paper, we propose a new classical-quantum hybrid method, named as power of sine Hamiltonian operator (PSHO), to evaluate the eigenvalues of a given Hamiltonian (Ĥ). In PSHO, for any reference state |φ0⟩, the normalized energy of the sinn(H^τ)|φ0⟩ state can be determined. With the increase of the power, the initial reference state can converge to the eigenstate with the largest |sin(Eiτ)| value in the coefficients of the expansion of |φ0⟩, and the normalized energy of the sinn(H^τ)|φ0⟩ state converges to E_i. The ground- and excited-state energies of a Hamiltonian can be determined by taking different τ values. The performance of the PSHO method is demonstrated by numerical calculations of the H₄ and LiH molecules. Compared with the current popular variational quantum eigensolver method, PSHO does not need to design the ansatz circuits and avoids the complex nonlinear optimization problems. PSHO has great application potential in near-term quantum devices.

Ultrafast Evaluation of Two-Photon Absorption with Simplified Time-Dependent Density Functional Theory

This work presents the theoretical background to evaluate two-photon absorption (2PA) cross-sections in the framework of simplified time-dependent density functional theory (sTD-DFT). Our new implementation allows the ultrafast evaluation of 2PA cross-sections for large molecules based on a regular DFT ground-state determinant as well as a variant employing our tight-binding sTD-DFT-xTX flavor for very large systems. The method is benchmarked against higher-level calculations for trans-stilbene and typical fluorescent protein chromophores. For eGFP, a quadrupolar chromophore and its branched version, the flavine mono-nucleotide, and the iLOV protein, we compare sTD-DFT 2PA spectra to experimental ones. This includes extension and testing of our all-atom quantum chemistry methodology for the evaluation of 2PA for a system of ∼2000 atoms, providing striking agreement with the experimental spectrum.

Variational Quantum Computation of Molecular Linear Response Properties on a Superconducting Quantum Processor

Simulating response properties of molecules is crucial for interpreting experimental spectroscopies and accelerating materials design. However, it remains a long-standing computational challenge for electronic structure methods on classical computers. While quantum computers hold the promise of solving this problem more efficiently in the long run, existing quantum algorithms requiring deep quantum circuits are infeasible for near-term noisy quantum processors. Herein, we introduce a pragmatic variational quantum response (VQR) algorithm for response properties, which circumvents the need for deep quantum circuits. Using this algorithm, we report the first simulation of linear response properties of molecules including dynamic polarizabilities and absorption spectra on a superconducting quantum processor. Our results indicate that a large class of important dynamical properties, such as Green's functions, are within the reach of near-term quantum hardware using this algorithm in combination with suitable error mitigation techniques.

New pecJ-n (n = 1, 2) Basis Sets for High-Quality Calculations of Indirect Nuclear Spin–Spin Coupling Constants Involving 31P and 29Si: The Advanced PEC Method

In this paper, we presented new J-oriented basis sets, pecJ-n (n = 1, 2), for phosphorus and silicon, purposed for the high-quality correlated calculations of the NMR spin–spin coupling constants involving these nuclei. The pecJ-n basis sets were generated using the modified version of the property-energy consistent (PEC) method, which was introduced in our earlier paper. The modifications applied to the original PEC procedure increased the overall accuracy and robustness of the generated basis sets in relation to the diversity of electronic systems. Our new basis sets were successfully tested on a great number of spin–spin coupling constants, involving phosphorus or/and silicon, calculated within the SOPPA(CCSD) method. In general, it was found that our new pecJ-1 and pecJ-2 basis sets are very efficient, providing the overall accuracy that can be characterized by MAEs of about 3.80 and 1.98 Hz, respectively, against the benchmark data obtained with a large dyall.aae4z⁺ basis set of quadruple-ζ quality.

Cherry-Picking Resolvents: Recovering the Valence Contribution in X-ray Two-Photon Absorption within the Core-Valence-Separated Equation-of-Motion Coupled-Cluster Response Theory

Calculations of first-order response wave functions in the X-ray regime often diverge within correlated frameworks such as equation-of-motion coupled-cluster singles and doubles (EOM-CCSD), a consequence of the coupling with the valence ionization continuum. Here, we extend our strategy of introducing a hierarchy of approximations to the EOM-EE-CCSD resolvent (or, inversely, the model Hamiltonian) involved in the response equations for the calculation of X-ray two-photon absorption (X2PA) cross sections. We exploit the frozen-core core-valence separation (fc-CVS) scheme to first decouple the core and valence Fock spaces, followed by a separate approximate treatment of the valence resolvent. We demonstrate the robust convergence of X-ray response calculations within this framework and compare X2PA spectra of small benchmark molecules with the previously reported density functional theory results.

Fluorine spin-spin coupling constants of pentafluorobenzene revisited at the ab initio correlated levels

All possible spin-spin coupling constants, ¹⁹ F-¹⁹ F, ¹⁹ F-¹³ C, and ¹⁹ F-¹ H, of pentafluorobenzene were calculated at five different levels of theory, HF, DFT, SOPPA (CCSD), CCSD, and the SOPPA (CCSD)-based composite scheme with taking into account solvent, vibrational, relativistic, and correlation corrections. Most corrections were next to negligible for the long-range couplings but quite essential for the one-bond carbon-fluorine coupling constants. Hartree-Fock calculations were found to be entirely unreliable, while DFT results were comparable in accuracy with the data obtained using the wave function-based methods.

Stereochemical dependence of substituent γ-effects in the 19 F NMR shielding constants

The substituent α-, β-, and γ-effects of the elements of the second and third periods on ¹⁹ F NMR chemical shifts are evaluated including the establishment of stereochemical dependence of γ-effect, the latter particularly important in stereochemical studies of fluorine-containing compounds. Benchmark calculations performed for a series of 32 simple inorganic fluorine-containing molecules demonstrated a markedly good correlation between calculated and experimental fluorine chemical shifts characterized by a mean absolute error of 22.5 ppm in the range of about 900 ppm, which corresponds to a 2.5% error in the percentage terms.

Oscillator Strengths in the Framework of Equation of Motion Multilevel CC3

We present an efficient implementation of the equation of motion oscillator strengths for the closed-shell multilevel coupled cluster singles and doubles with perturbative triples method (MLCC3) in the electronic structure program e^T. The orbital space is split into an active part treated with CC3 and an inactive part computed at the coupled cluster singles and doubles (CCSD) level of theory. Asymptotically, the CC3 contribution scales as floating-point operations, where n_V is the total number of virtual orbitals while n_v and n_o are the number of active virtual and occupied orbitals, respectively. The CC3 contribution, thus, only scales linearly with the full system size and can become negligible compared to the cost of CCSD. We demonstrate the capabilities of our implementation by calculating the ultraviolet–visible spectrum of azobenzene and a core excited state of betaine 30 with more than 1000 molecular orbitals.

Recent Advances in Cartesian-Grid DFT in Atoms and Molecules

In the past several decades, density functional theory (DFT) has evolved as a leading player across a dazzling variety of fields, from organic chemistry to condensed matter physics. The simple conceptual framework and computational elegance are the underlying driver for this. This article reviews some of the recent developments that have taken place in our laboratory in the past 5 years. Efforts are made to validate a viable alternative for DFT calculations for small to medium systems through a Cartesian coordinate grid- (CCG-) based pseudopotential Kohn–Sham (KS) DFT framework using LCAO-MO ansatz. In order to legitimize its suitability and efficacy, at first, electric response properties, such as dipole moment (μ), static dipole polarizability (α), and first hyperpolarizability (β), are calculated. Next, we present a purely numerical approach in CCG for proficient computation of exact exchange density contribution in certain types of orbital-dependent density functionals. A Fourier convolution theorem combined with a range-separated Coulomb interaction kernel is invoked. This takes motivation from a semi-numerical algorithm, where the rate-deciding factor is the evaluation of electrostatic potential. Its success further leads to a systematic self-consistent approach from first principles, which is desirable in the development of optimally tuned range-separated hybrid and hyper functionals. Next, we discuss a simple, alternative time-independent DFT procedure, for computation of single-particle excitation energies, by means of “adiabatic connection theorem” and virial theorem. Optical gaps in organic chromophores, dyes, linear/non-linear PAHs, and charge transfer complexes are faithfully reproduced. In short, CCG-DFT is shown to be a successful route for various practical applications in electronic systems.

Efficient implementation of molecular CCSD gradients with Cholesky-decomposed electron repulsion integrals

We present an efficient implementation of ground and excited state coupled cluster singles and doubles (CCSD) gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition and density fitting are both inner projection methods, and, thus, similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3 O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. Furthermore, the computational time per optimization cycle is significantly lower compared to other implementations based on an inner projection method. We showcase the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory.

Nature of chemical bond and potential barrier in an invariant energy-orbital picture

Physical nature of the chemical bond and potential barrier is studied in terms of energy natural orbitals (ENOs), which are extracted from highly correlated electronic wavefunctions. ENO provides an objective one-electron picture about energy distribution in a molecule, just as the natural orbitals (NOs) represent one electron view about electronic charge distribution. ENO is invariant in the same sense as NO is, that is, ENOs converge to the exact ones as a series of approximate wavefunctions approach the exact one, no matter how the methods of approximation are adopted. Energy distribution analysis based on ENO can give novel insights about the nature of chemical bonding and formation of potential barriers, besides information based on the charge distribution alone. With ENOs extracted from full configuration interaction wavefunctions in a finite yet large enough basis set, we analyze the nature of chemical bonding of three low-lying electronic states of a hydrogen molecule, all being in different classes of the so-called covalent bond. The mechanism of energy lowering in bond formation, which gives a binding energy, is important, yet not the only concern for this small molecule. Another key notion in chemical bonding is whether a potential basin is well generated stiff enough to support a vibrational state(s) on it. Based on the virial theorem in the adiabatic approximation and in terms of the energy and force analyses with ENOs, we study the roles of the electronic kinetic energy and its nuclear derivative(s) on how they determine the curvature (or the force constant) of the potential basins. The same idea is applied to the potential barrier and, thereby, the transition states. The rate constant within the transition-state theory is formally shown to be described in terms of the electronic kinetic energy and the nuclear derivatives only. Thus, the chemical bonding and rate process are interconnected behind the scenes. Besides this aspect, we pay attention to (1) the effects of electron correlation that manifests itself not only in the orbital energy but also in the population of ENOs and (2) the role of nonadiabaticity (diabatic state mixing), resulting in double basins and a barrier on a single potential curve in bond formation. These factors differentiate a covalent bond into subclasses.

Cluster perturbation theory. VIII. First order properties for a coupled cluster state

We have extended cluster perturbation (CP) theory to comprehend the calculation of first order properties (FOPs). We have determined CP FOP series where FOPs are determined as a first energy derivative and also where the FOPs are determined as a generalized expectation value of the external perturbation operator over the coupled cluster state and its biorthonormal multiplier state. For S(D) orbital excitation spaces, we find that the CP series for FOPs that are determined as a first derivative in general in second order have errors of a few per cent in the singles and doubles correlation contribution relative to the targeted coupled cluster (CC) results. For a SD(T) orbital excitation space, we find that the CP series for FOPs determined as a generalized expectation value in second order have errors of about ten percent in the triples correlation contribution relative to the targeted CC results. These second order models therefore constitute viable alternatives for determining high quality FOPs.

Analyzing the Behavior of Spin Phases in External Magnetic Fields by Means of Spin-Constrained States

During molecular dissociation in the presence of an external uniform magnetic field, electrons flip their spin antiparallel to the magnetic field because of the stabilizing influence of the spin Zeeman operator. Although generalized Hartree-Fock descriptions furnish the optimal mean-field energetic description of such bond-breaking processes, they are allowed to break Ŝ_z symmetry, leading to intricate and unexpected spin phases and phase transitions. In this work, we show that the behavior of these molecular spin phases can be interpreted in terms of spin phase diagrams constructed by constraining states to target expectation values of projected spin. The underlying constrained states offer a complete electronic characterization of the spin phases and spin phase transitions, as they can be analyzed using standard quantum chemical tools. Because the constrained states effectively span the entire phase space, they could provide an excellent starting point for post-Hartree-Fock methods aimed at gaining more electron correlation or regaining spin symmetry.

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.

Quantum Chemical Approaches to the Calculation of NMR Parameters: From Fundamentals to Recent Advances

Quantum chemical methods for the calculation of indirect NMR spin–spin coupling constants and chemical shifts are always in progress. They never stay the same due to permanently developing computational facilities, which open new perspectives and create new challenges every now and then. This review starts from the fundamentals of the nonrelativistic and relativistic theory of nuclear magnetic resonance parameters, and gradually moves towards the discussion of the most popular common and newly developed methodologies for quantum chemical modeling of NMR spectra.

The second derivative of the electronic energy with respect to the compression scaling factor in the XP-PCM model: Theory and applications to compression response functions of atoms

We present the analytical theory for the second derivative of the electronic energy with respect to the scaling factor of the compression cavity within the eXtreme pressure polarizable continuum model (XP-PCM) for the study of compressed atomic and molecular systems. The theory has been exploited to study compression response functions describing how the atomic/molecular properties are effected by an external pressure. The response functions considered include the atomic compressibility and the pressure coefficients of the ionization energy (IE) and electron affinity (EA). The theory has been validated by numerical application to compressed neon, argon, and krypton atoms.

Cluster Perturbation Theory. VII. The convergence of Cluster Perturbation Expansions

The convergence of the recently developed cluster perturbation CP expansions (Pawlowski et al, J. Chem. Phys. 150 134108(2019)) is analyzed with the double purpose of developing the mathematical tools and concepts needed to describe these expansions at general order and to identify the factors that define the rate of convergence of CP series. To this end, the CP energy, amplitude, and Lagrangian multiplier equations as functions of the perturbation strength are developed. By determining the critical points, defined as the perturbation strengths for which the Jacobian become singular, the rate of convergence as well as the intruder and critical states are determined for five simple molecules: BH, CO, H2O, NH3, and HF. To describe the patterns of convergence for these expansions at orders lower than the high-order asymptotic limit, a model is developed, where the perturbation corrections arise from two critical points. It is shown that this model allows rationalization of the behavior of the perturbation corrections at much lower order than required for the onset of the asymptotic convergence. For the H2O, CO, and HF molecules, the pattern and rate of convergence is defined by critical states where the Fock-operator underestimates the excitation energies, whereas the pattern and rate of convergence for BH is defined by critical states where the Fock-operator overestimates the excitation energy. For the NH3 molecule, both forms of critical points are required to describe the convergence behavior up to at least order 25.

Theory of electronic resonances: fundamental aspects and recent advances

Electronic resonances are states that are unstable towards loss of electrons. They play critical roles in high-energy environments across chemistry, physics, and biology but are also relevant to processes under ambient conditions that involve unbound electrons. This feature article focuses on complex-variable techniques such as complex scaling and complex absorbing potentials that afford a treatment of electronic resonances in terms of discrete square-integrable eigenstates of non-Hermitian Hamiltonians with complex energy. Fundamental aspects of these techniques as well as their integration into molecular electronic-structure theory are discussed and an overview of some recent developments is given: analytic gradient theory for electronic resonances, the application of rank-reduction techniques and quantum embedding to them, as well as approaches for evaluating partial decay widths.

Extending NMR Quantum Computation Systems by Employing Compounds with Several Heavy Metals as Qubits

Nuclear magnetic resonance (NMR) is a spectroscopic method that can be applied to several areas. Currently, this technique is also being used as an experimental quantum simulator, where nuclear spins are employed as quantum bits or qubits. The present work is devoted to studying heavy metal complexes as possible candidates to act as qubit molecules. Nuclei such 113Cd, 199Hg, 125Te, and 77Se assembled with the most common employed nuclei in NMR-QIP implementations (1H, 13C, 19F, 29Si, and 31P) could potentially be used in heteronuclear systems for NMR-QIP implementations. Hence, aiming to contribute to the development of future scalable heteronuclear spin systems, we specially designed four complexes, based on the auspicious qubit systems proposed in our previous work, which will be explored by quantum chemical calculations of their NMR parameters and proposed as suitable qubit molecules. Chemical shifts and spin–spin coupling constants in four complexes were examined using the spin–orbit zeroth-order regular approximation (ZORA) at the density functional theory (DFT) level, as well as the relaxation parameters (T1 and T2). Examining the required spectral properties of NMR-QIP, all the designed complexes were found to be promising candidates for qubit molecules.

Theoretical analysis of the long-distance limit of NMR chemical shieldings

After some years of controversy, it was recently demonstrated how to obtain the correct long-distance limit [point-dipole approximation (PDA)] of pseudo-contact nuclear magnetic resonance chemical shifts from rigorous first-principles quantum mechanics [Lang et al., J. Phys. Chem. Lett. 11, 8735 (2020)]. This result confirmed the classical Kurland-McGarvey theory. In the present contribution, we elaborate on these results. In particular, we provide a detailed derivation of the PDA both from the Van den Heuvel-Soncini equation for the chemical shielding tensor and from a spin Hamiltonian approximation. Furthermore, we discuss in detail the PDA within the approximate density functional theory and Hartree-Fock theories. In our previous work, we assumed a relatively crude effective nuclear charge approximation for the spin-orbit coupling operator. Here, we overcome this assumption by demonstrating that the derivation is also possible within the fully relativistic Dirac equation and even without the assumption of a specific form for the Hamiltonian. Crucial ingredients for the general derivation are a Hamiltonian that respects gauge invariance, the multipolar gauge, and functional derivatives of the Hamiltonian, where it is possible to identify the first functional derivative with the electron number current density operator. The present work forms an important foundation for future extensions of the Kurland-McGarvey theory beyond the PDA, including induced magnetic quadrupole and higher moments to describe the magnetic hyperfine field.

Cluster perturbation theory. VI. Ground-state energy series using the Lagrangian

In this paper, we derive alternative cluster perturbation series to the energy cluster perturbation (ECP) series of Paper I. The ECP series were derived using the standard coupled cluster energy framework. Here, we use the coupled cluster Lagrangian framework to derive the Lagrangian cluster perturbation (LCP) series and show that a slightly modified order concept means that the ECP and the LCP series become identical. Using the Lagrangian, we also derive a perturbation series where total cluster amplitudes and multipliers are determined through the same orders as dictated by the 2n+1/2n+2 rule, the VCP series. The VCP energies have errors that are bilinear in the errors of the total cluster amplitudes and multipliers. Test calculations have been performed for S(D) and SD(T) orbital excitation spaces. The convergence of the test calculations can be divided into two groups. The first group contains molecules that exhibit a slow monotonic geometric convergence pattern in both orbital spaces. The second group, containing the rest of the molecules, exhibits a fast low order convergence. For the S(D) calculations, the deviations in fourth order from the coupled cluster singles and doubles energy are below 1.3 per cent for the LCP series and 0.3 per cent for the VCP series, respectively. For the SD(T) calculations, the second group hardly shows any difference between the low order convergence of the three series ECP, LCP, and LCP and at fourth order the deviations from the triples correlation energy are less than 1 per cent.

Probing Molecular Chirality of Ground and Electronically Excited States in the UV-vis and X-ray Regimes: An EOM-CCSD Study

We present several strategies for computing electronic circular dichroism (CD) spectra across different frequency ranges at the equation-of-motion coupled-cluster singles and doubles level of theory. CD spectra of both ground and electronically excited states are discussed. For selected cases, the approach is compared with coupled-cluster linear response results as well as time-dependent density functional theory. The extension of the theory to include the effect of spin-orbit coupling is presented and illustrated by calculations of X-ray CD spectra at the L-edge.

Revisiting the Fundamental Nature of Metal‐Ligand Bonding: An Impartial and Automated Fitting Procedure for Angular Overlap Model Parameters

The properties and reactivities of transition metal complexes are often discussed in terms of Ligand Field Theory (LFT), and with ab initio LFT a direct connection to quantum chemical wavefunctions was recently established. The Angular Overlap Model (AOM) is a widely used, ligand‐specific parameterization scheme of the ligand field splitting that has, however, been restricted by the availability and resolution of experimental data. Using ab initio LFT, we present here a generalised, symmetry‐independent and automated fitting procedure for AOM parameters that is even applicable to formally underdetermined or experimentally inaccessible systems. This method allows quantitative evaluations of assumptions commonly made in AOM applications, for example, transferability or the relative magnitudes of AOM parameters, and the response of the ligand field to structural or electronic changes. A two‐dimensional spectrochemical series of tetrahedral halido metalates ([M^IIX₄]²⁻, M=Mn−Cu) served as a case study. A previously unknown linear relationship between the halide ligands’ chemical hardness and their AOM parameters was found. The impartial and automated procedure for identifying AOM parameters introduced here can be used to systematically improve our understanding of ligand–metal interactions in coordination complexes.

Computational NMR of charged systems

This review covers NMR computational aspects of charged systems-carbocations, heterocations, and heteroanions, which were extensively studied in a number of laboratories worldwide, first of all, at the Loker Hydrocarbon Research Institute in California directed for several decades by a distinguished scientist, the Nobel laureate George Andrew Olah. The first part of the review briefly outlines computational background of the modern theoretical methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the historical results, advances, and perspectives of the computational NMR of classical carbocations like methyl cation, CH₃⁺ , and protonated methane, CH₅⁺ , together with their numerous homologs and derivatives. The third and the forth parts of this survey are focused on the NMR computational aspects of accordingly, heterocations and heteroanions, the organic and inorganic ions with a charge localized mainly on heteroatoms like boron, oxygen, nitrogen, and heavier elements.

Extended Benchmark Set of Main-Group Nuclear Shielding Constants and NMR Chemical Shifts and Its Use to Evaluate Modern DFT Methods

An extended theoretical benchmark set, NS372, for light main-group nuclear shieldings and NMR shifts has been constructed based on high-level GIAO-CCSD(T)/pcSseg-3//CCSD(T)/cc-pVQZ reference data. After removal of the large static-correlation cases O₃, F₃^-, and BH from the statistical evaluations for the ¹⁷O, ¹⁹F, and ¹¹B subsets, the benchmark comprises overall 372 shielding values in 117 molecules with a wide range of electronic-structure situations, containing 124 ¹H, 14 ¹¹B, 93 ¹³C, 43 ¹⁵N, 31 ¹⁷O, 47 ¹⁹F, 14 ³¹P, and 6 ³³S shielding constants. The CCSD(T)/pcSseg-3 data are shown to be close to the basis-set and method limit and thus provide an excellent benchmark to evaluate more approximate methods, such as density functional approaches. This dataset has been used to evaluate Hartree-Fock (HF) and MP2, and a wide range of exchange-correlation functionals from local density approximation (LDA) to generalized gradient approximations (GGAs) and meta-GGAs (focusing on their current-density functional implementations), as well as global hybrid, range-separated hybrid, local hybrid, and double-hybrid functionals. Starting with absolute shielding constants, the DSD-PBEP86 double hybrid is confirmed to provide the highest accuracy, with an aggregate relative mean absolute error (rel. MAE) of only 0.9%, followed by MP2 (1.1%). MP2 and double hybrids only show larger errors for a few systems with the largest static-correlation effects. The double-hybrid B2GP-PLYP, the two local hybrids cLH12ct-SsirPW92 and cLH12ct-SsifPW92, and the current-density functional meta-GGA cB97M-V follow closely behind (all 1.5%), as do some further functionals, cLH20t and cMN15-L (both 1.6%), as well as B2PLYP and KT3 (both 2.0%). Functionals on the lower rungs of the usual ladder offer the advantage of lower computational cost and access to larger molecules. Closer examination also reveals the best-performing methods for individual nuclei in the test set. Different ways of treating τ-dependent functionals are evaluated. When moving from absolute shielding constants to chemical shifts, some of the methods can benefit from systematic error compensation, and the overall error range somewhat narrows. Further methods now achieve the 2% threshold of relative MAEs, including functionals based on TPSS (TPSSh, cmPSTS).

Accelerating molecular property calculations with semiempirical preconditioning

Computing ab initio molecular linear response properties, e.g., electronic excitation energies and transition dipole moments, requires the solution of large eigenvalue problems or large systems of equations. These large eigenvalue problems or large systems of equations are commonly solved iteratively using Krylov space algorithms, such as the Davidson algorithm for eigenvalue problems. A critical ingredient in Krylov space algorithms is the preconditioner, which is used to generate optimal update vectors in each iteration. We propose to use semiempirical approximations as preconditioners to accelerate the calculation of ab initio properties. The crucial advantage to improving the preconditioner is that the converged result is unchanged, so there is no trade-off between accuracy and speedup. We demonstrate our approach by accelerating the calculation of electronic excitation energies and electric polarizabilities from linear response time-dependent density functional theory using the simplified time-dependent density functional theory semiempirical model. For excitation energies, the semiempirical preconditioner reduces the number of iterations on average by 37% and up to 70%. The semiempirical preconditioner reduces the number of iterations for computing the polarizability by 15% on average and up to 33%. Moreover, we show that the preconditioner can be further improved by tuning the empirical parameters that define the semiempirical model, leading to an additional reduction in the number of iterations by about 20%. Our approach bridges the gap between semiempirical models and ab initio methods and charts a path toward combining the speed of semiempirical models with the accuracy of ab initio methods.

Current density and molecular magnetic properties

We give an overview of the molecular response to an external magnetic field perturbing quantum mechanical systems. We present state-of-the-art methods for calculating magnetically-induced current-density susceptibilities. We discuss the essence and properties of current-density susceptibilities and how molecular magnetic properties can be calculated from them. We also review the theory of spin-current densities, how relativity affects current densities and magnetic properties. An overview of the magnetic ring-current criterion for aromaticity is given, which has implications on theoretical and experimental research. The recently reported theory of antiaromaticity and how molecular symmetry affects the magnetic response are discussed and applied to closed-shell paramagnetic molecules. The topology of magnetically induced current densities and its consequences for molecular magnetic properties are also presented with twisted and toroidal molecules as examples.

Measurement and Theory of Resonance Raman Optical Activity for Gases, Liquids, and Aggregates. What It Tells about Molecules

Resonance Raman Optical Activity (RROA) appeared as a natural extension of the nonresonance branch. It combines the structural sensitivity of chiroptical spectroscopy with the signal enhancement coming from the resonance of molecular electronic transitions with the excitation laser light. However, the idea has been hampered by many technical and theoretical problems that are being clarified only in recent years. We provide the theoretical basis and several examples documenting the problems, achievements, and potential of RROA, in particular in biomolecular studies.

Theoretical Investigation of Glycine Micro-Solvated. Energy and NMR Spin Spin Coupling Constants Calculations

Glycine in its neutral form can exist in the gas phase while its zwitterion form is more stable in water solution, but how many waters are actually necessary to stabilize the zwitterionic structure in the gas phase? Are the intramolecular isotropic spin spin coupling constants sensitive enough to accuse the change in the environment? or the conformer observed? These and related questions have been investigated by a computational study at the level of density functional theory employing the B3LYP functional and the 6-31++G**-J basis set. We found that at least two water molecules explicitly accounted for in the super-molecule structure are necessary to stabilize both conformers of glycine within a water polarizable continuum model. At least half of the SSCCs of both conformers are very stable to changes in the environment and at least four of them differ significantly between Neutral and Zwitterion conformation.

Self-assembling, structure and nonlinear optical properties of fluorescent organic nanoparticles in water

Owing to their intense emission, low toxicity and solubility in aqueous medium, fluorescent organic nanoparticles (FONs) have emerged as promising alternatives to inorganic ones for the realization of exogenous probes for bioimaging applications. However, the intimate structure of FONs in solution, as well as the role played by intermolecular interactions on their optical properties, remains challenging to study. Following a recent Second-Harmonic Scattering (SHS) investigation led by two of us [Daniel et al., ACS Photonics, 2015, 2, 1209], we report herein a computational study of the structural organization and second-order nonlinear optical (NLO) properties of FONs based on dipolar chromophores incorporating a hydrophobic triphenylamine electron-donating unit and a slightly hydrophilic aldehyde electron-withdrawing unit at their extremities. Molecular dynamics simulations of the FON formation in water are associated with quantum chemical calculations, to provide insight into the molecular aggregation process, the molecular orientation of the dipolar dyes within the nanoparticles, and the dynamical behavior of their NLO properties. Moreover, the impact of intermolecular interactions on the NLO responses of the FONs is investigated by employing the tight-binding version of the recently developed simplified time-dependent density functional theory (sTD-DFT) approach, allowing the all-atom quantum mechanics treatment of nanoparticles.

Relativistic study of parity-violating nuclear spin-rotation tensors

We present a four-component relativistic approach to describe the effects of the nuclear spin-dependent parity-violating (PV) weak nuclear forces on nuclear spin-rotation (NSR) tensors. The formalism is derived within the four-component polarization propagator theory based on the Dirac-Coulomb Hamiltonian. Such calculations are important for planning and interpretation of possible future experiments aimed at stringent tests of the standard model through the observation of PV effects in NSR spectroscopy. An exploratory application of this theory to the chiral molecules H₂X₂ (X = ¹⁷O, ³³S, ⁷⁷Se, ¹²⁵Te, and ²⁰⁹Po) illustrates the dramatic effect of relativity on these contributions. In particular, spin-free and spin-orbit effects are even of opposite signs for some dihedral angles, and the latter fully dominate for the heavier nuclei. Relativistic four-component calculations of isotropic nuclear spin-rotation constants, including parity-violating electroweak interactions, give frequency differences of up to 4.2 mHz between the H₂Po₂ enantiomers; on the nonrelativistic level of theory, this energy difference is 0.1 mHz only.

Information Entropy in Chemistry: An Overview

Basic applications of the information entropy concept to chemical objects are reviewed. These applications deal with quantifying chemical and electronic structures of molecules, signal processing, structural studies on crystals, and molecular ensembles. Recent advances in the mentioned areas make information entropy a central concept in interdisciplinary studies on digitalizing chemical reactions, chemico-information synthesis, crystal engineering, as well as digitally rethinking basic notions of structural chemistry in terms of informatics.

Second harmonic generation theory for a helical macromolecule with high sensitivity to structural disorder

Microscopic theory for the second harmonic generation in a helical molecular system is developed in the minimal coupling representation including non-local interaction effects. At the second order to the field we find a compact expression which combines dipolar, quadrupolar and magnetic contributions. A detailed derivation of the response is performed to specifically isolate the quadratic coupling terms, which we denote as the K coupling. Applying the theory to a helical macromolecule we find that the dipolar and quadrupolar contributions reflect the symmetry properties of the system and its homogeneity, while the K coupling contribution reveals inhomogeneities of the system.

Nitrogen K-Edge X-ray Absorption Spectra of Ammonium and Ammonia in Water Solution: Assessing the Performance of Polarizable Embedding Coupled Cluster Methods

The recent development of liquid jet and liquid leaf sample delivery systems allows for accurate measurements of soft X-ray absorption spectra in transmission mode of solutes in a liquid environment. As this type of measurement becomes increasingly accessible, there is a strong need for reliable theoretical methods for assisting in the interpretation of the experimental data. Coupled cluster methods have been extensively developed over the past decade to simulate X-ray absorption in the gas phase. Their performance for solvated species, on the contrary, remains largely unexplored. Here, we investigate the current state of the art of coupled cluster modeling of nitrogen K-edge X-ray absorption of aqueous ammonia and ammonium based on quantum mechanics/molecular mechanics, where both the level of coupled cluster calculations and polarizable embedding are scrutinized. The results are compared to existing experimental data as well as simulations based on transition potential density functional theory.

Density Functional Theory Study of Substitution Effects on the Second-Order Nonlinear Optical Properties of Lindquist-Type Organo-Imido Polyoxometalates

Density functional theory and time-dependent density functional theory have been enacted to investigate the effects of donor and acceptor on the first hyperpolarizability of Lindquist-type organo-imido polyoxometalates (POMs). These calculations employ a range-separated hybrid exchange-correlation functional (ωB97X-D), account for solvent effects using the implicit polarizable continuum model, and analyze the first hyperpolarizabilities by using the two-state approximation. They highlight the beneficial role of strong donors as well as of π-conjugated spacers (CH=CH rather than C≡C) on the first hyperpolarizabilities. Analysis based on the unit sphere representation confirms the one-dimensional push-pull π-conjugated character of the POMs substituted by donor groups and the corresponding value of the depolarization ratios close to 5. Furthermore, the use of the two-state approximation is demonstrated to be suitable for explaining the origin of the variations of the first hyperpolarizabilities as a function of the characteristics of a unique low-energy charge-transfer excited state and to attribute most of the first hyperpolarizability changes to the difference of dipole moment between the ground and that charge-transfer excited state.

Spatial Contributions to 1H NMR Chemical Shifts of Free-Base Porphyrinoids

A recently developed methodology for calculating, analyzing, and visualizing nuclear magnetic shielding densities is used for studying spatial contributions including ring-current contributions to 1H nuclear magnetic resonance (NMR) chemical shifts of aromatic and anti-aromatic free-base porphyrinoids. Our approach allows a visual inspection of the spatial origin of the positive (shielding) and negative (deshielding) contributions to the nuclear magnetic shielding constants. Diatropic and paratropic current-density fluxes yield both shielding and deshielding contributions implying that not merely the tropicity of the current density determines whether the contribution has a shielding or deshielding character. Instead the shielding or deshielding contribution is determined by the direction of the current-density flux with respect to the studied nucleus.