Static polarizabilities and C6 dispersion coefficients using the algebraic-diagrammatic construction scheme for the complex polarization propagator

Efficient and Parallel Implementation of Real and Complex Response Functions Employing the Second-Order Algebraic-Diagrammatic Construction Scheme for the Polarization Propagator

We present the implementation of an efficient matrix-folded formalism for the evaluation of complex response functions and the calculation of transition properties at the level of the second-order algebraic-diagrammatic construction (ADC(2)) scheme. The underlying algorithms, in combination with the adopted hybrid MPI/OpenMP parallelization strategy, enabled calculations of the UV/vis spectra of a guanine oligomer series ranging up to 1032 contracted basis functions, thereby utilizing vast computational resources from up to 32,768 CPU cores. Further analysis of the convergence behavior of the involved iterative subspace algorithms revealed the superiority of a frequency-separated treatment of response equations even for a large spectral window, including 101 frequencies. We demonstrate the applicability to general quantum mechanical operators by the first reported electronic circular dichroism spectrum calculated with a complex polarization propagator approach at the ADC(2) level of theory.

Cotton-Mouton Effect Using Algebraic Diagrammatic Construction Schemes in the Intermediate State Representation Formalism

The Cotton-Mouton effect is theoretically investigated for a selected set of molecules by using a novel computational methodology based on algebraic diagrammatic construction (ADC) schemes in the intermediate state representation (ISR) formulation. Therefore, the electronic contributions to the frequency-dependent polarizabilities and, for the first time, to the magnetizabilities as well as mixed electric and magnetic hypermagnetizabilities have been computed in the ADC/ISR framework. In addition to calculation of the Cotton-Mouton constant and the birefringence, the gauge origin dependence of the computed tensors and the applied methodology are thoroughly investigated. The new ADC/ISR methodology, employing the recently presented responsefun package, is applied to a test set of Ne and small molecules (H2, HF, O2, CO2, and benzene) and compared to data from the experiment as well as other ab initio methods. The presented theoretical ab initio ADC/ISR approach is a substantial extension of the available computational methods for the investigation of complex nonlinear properties, however, with a gauge origin dependence inherent to the method that decreases with increasing perturbation order.

responsefun: Fun with Response Functions in the Algebraic Diagrammatic Construction Framework

We present the open-source responsefun package, which implements a universally applicable procedure for computing molecular response properties within the algebraic diagrammatic construction (ADC) framework, exploiting the intermediate state representation (ISR) approach. With symbolic mathematics, the user can simply enter textbook sum-over-states (SOS) expressions from time-dependent perturbation theory, which are then automatically translated into the corresponding symbolic ADC/ISR formulations. Using the data structures provided by the hybrid Python/C++ module adcc for calculating excited states with ADC, the specified response property is directly evaluated, and the result is returned to the user. Employing the novel responsefun package, we present the first ADC/ISR calculations of second-order hyperpolarizability tensors and three-photon-absorption matrix elements.

Efficient Protocol for Computing MCD Spectra in a Broad Frequency Range Combining Resonant and Damped CC2 Quadratic Response Theory

Coupled cluster response theory offers a path to high-accuracy calculations of spectroscopic properties, such as magnetic circular dichroism (MCD). However, divergence or slow convergence issues are often encountered for electronic transitions in high-energy regions with a high density of states. This is here addressed for MCD by an implementation of damped quadratic response theory for resolution-of-identity coupled cluster singles-and-approximate-doubles (RI-CC2), along with an implementation of the MCD A term from resonant response theory. Combined, damped and resonant response theory calculations provide an efficient strategy to obtain MCD spectra over a broad frequency range and for systems that include highly symmetric molecules with degenerate excited states. The protocol is illustrated by application to zinc tetrabenzoporphyrin in the energy region of 2-8 eV and comparison to experimental data. Timings are reported for the resonant and damped approaches, showing that a greater part of the calculation time is consumed by the construction of the building blocks for the final MCD ellipticity. A recommendation on how to use the procedure is outlined.

Solving response expressions in the ADC/ISR framework

We present an implementation for the calculation of molecular response properties using the algebraic-diagrammatic construction (ADC)/intermediate state representation approach. For the second-order ADC model [ADC(2)], a memory-efficient ansatz avoiding the storage of double excitation amplitudes is investigated. We compare the performance of different numerical algorithms for the solution of the underlying response equations for ADC(2) and show that our approach also strongly improves the convergence behavior for the investigated algorithms compared with the standard implementation. All routines are implemented in an open-source Python library.

Density Functional Response Calculations of Dispersion Coefficients C6 and C9 of Closed- and Open-Shell Systems

Dipole polarizabilities and C₆ and C₉ dispersion coefficients are computed for closed- and open-shell atoms and molecules, using dynamic (time-dependent) density functional (TD-DFT) linear response theory as implemented in the response module of the NWChem quantum chemistry package. The response module is capable of accurate calculations of these properties, based on spin-restricted and spin-unrestricted formalisms. The calculated static polarizabilities and dispersion coefficients are compared to available experimental and other theoretical data. The behavior of the dynamic polarizability at imaginary frequencies is analyzed for differently sized closed- and open-shell systems. An interpolation method enforcing the monotonic decrease of the polarizability with increasing imaginary frequency is beneficial for the integration used to obtain C₆ and C₉. Scaling of the TD-DFT data by ratios of the static polarizability, which can be calculated with a variety of methods, including highly accurate theories, may be used as a leading-order correction.

Accurate Spectral Properties within Double-Hybrid Density Functional Theory: A Spin-Scaled Range-Separated Second-Order Algebraic-Diagrammatic Construction-Based Approach

Our second-order algebraic-diagrammatic construction [ADC(2)]-based double-hybrid (DH) ansatz (J. Chem. Theory Comput. 2019, 15, 4440. DOI: 10.1021/acs.jctc.9b00391) is combined with range-separation techniques. In the present scheme, both the exchange and the correlation contributions are range-separated, while spin-scaling approaches are also applied. The new methods are thoroughly tested for the most popular benchmark sets including 250 singlet and 156 triplet excitations, as well as 80 oscillator strengths. It is demonstrated that the range separation for the correlation contributions is highly recommended for both the genuine and the ADC(2)-based DH approaches. Our results show that the latter scheme slightly but consistently outperforms the former one for single excitation dominated transitions. Furthermore, states with larger fractions of double excitations are assessed as well, and challenging charge-transfer excitations are also discussed, where the recently proposed spin-scaled long-range corrected DHs fail. The suggested iterative fourth-power scaling RS-PBE-P86/SOS-ADC(2) method, using only three adjustable parameters, provides the most robust and accurate excitation energies within the DH theory. In addition, the relative error of the oscillator strengths is reduced by 65% compared to the best genuine DH functionals.

Benchmarking Correlated Methods for Static and Dynamic Polarizabilities: The T145 Data Set Evaluated with RPA, RPA(D), HRPA, HRPA(D), SOPPA, SOPPA(CC2), SOPPA(CCSD), CC2, and CCSD

Due to the importance of predicting static and dynamic polarizabilities, the performance of various correlated linear response methods including random phase approximation (RPA), RPA(D), higher-order random phase approximation (HRPA), HRPA(D), second-order polarization propagator approximation (SOPPA), SOPPA(CC2), SOPPA(CCSD), CC2, and CCSD has been evaluated against CCSD(T) (static case) and CCSD (dynamic cases) for the T145 set of 145 organic molecules. The benchmark reveals that the HRPA(D) method has the best performance for both static and dynamic polarizabilities apart from CCSD. RPA(D) ranks second for the dynamic cases and third for the static case. Using coupled-cluster amplitudes in SOPPA(CCSD) and SOPPA(CC2), the SOPPA results are significantly improved. The HRPA method has the largest deviations from the reference values for both cases. In general, according to the performance and computational cost of the methods, the HRPA(D) and RPA(D) methods are proposed for calculations of static and dynamic polarizabilities of this and similar sets of molecules.

Damped (linear) response theory within the resolution-of-identity coupled cluster singles and approximate doubles (RI-CC2) method

An implementation of a complex solver for the solution of the linear equations required to compute the complex response functions of damped response theory is presented for the resolution-of-identity (RI) coupled cluster singles and approximate doubles (CC2) method. The implementation uses a partitioned formulation that avoids the storage of double excitation amplitudes to make it applicable to large molecules. The solver is the keystone element for the development of the damped coupled cluster response formalism for linear and nonlinear effects in resonant frequency regions at the RI-CC2 level of theory. Illustrative results are reported for the one-photon absorption cross section of C₆₀, the electronic circular dichroism of n-helicenes (n = 5, 6, 7), and the C₆ dispersion coefficients of a set of selected organic molecules and fullerenes.

Core–valence-separated coupled-cluster-singles-and-doubles complex-polarization-propagator approach to X-ray spectroscopies

The iterative subspace algorithm to solve the CCSD complex linear response equations has been modified to include a core–valence separation projection step to overcome convergence problems. Illustrative results are reported for XAS, XCD, XES and RIXS.

Combined Density Functional and Algebraic-Diagrammatic Construction Approach for Accurate Excitation Energies and Transition Moments

A composite of time-dependent density functional theory (TDDFT) and the second-order algebraic-diagrammatic construction [ADC(2)] approach is presented for efficient calculation of spectral properties of molecules. Our method can be regarded as a new excited-state double-hybrid (DH) approach or a dressed TDDFT scheme, but it can also be interpreted as an empirically tuned ADC(2) model. Several combinations of exchange-correlation functionals and spin-scaling schemes are explored. Our best-performing method includes the Perdew, Burke, and Ernzerhof exchange and Perdew's 1986 correlation functional and employs the scaled-opposite-spin approximation for the higher-order terms. The computation time of the new method scales as the fourth power of the system size, and an efficient cost-reduction approach is also presented, which further speeds up the calculations. Our benchmark calculations show that the proposed model outperforms not only the existing DH approaches and ADC(2) variants but also the considerably more expensive coupled-cluster methods.

Deriving Force-Field Parameters from First Principles Using a Polarizable and Higher Order Dispersion Model

In this work we propose a strategy based on quantum mechanical (QM) calculations to parametrize a polarizable force field for use in molecular dynamics (MD) simulations. We investigate the use of multiple atoms-in-molecules (AIM) strategies to partition QM determined molecular electron densities into atomic subregions. The partitioned atomic densities are subsequently used to compute atomic dispersion coefficients from effective exchange-hole-dipole moment (XDM) calculations. In order to derive values for the repulsive van der Waals parameters from first principles, we use a simple volume relation to scale effective atomic radii. Explicit inclusion of higher order dispersion coefficients was tested for a series of alkanes, and we show that combining C6 and C8 attractive terms together with a C11 repulsive potential yields satisfying models when used in combination with our van der Waals parameters and electrostatic and bonded parameters as directly obtained from quantum calculations as well. This result highlights that explicit inclusion of higher order dispersion terms could be viable in simulation, and it suggests that currently available QM analysis methods allow for first-principles parametrization of molecular mechanics models.

Polarizability of Helium, Neon, and Argon: New Perspectives for Gas Metrology

With dielectric-constant gas thermometry, the molar polarizability of helium, neon, and argon has been determined with relative standard uncertainties of about 2 parts per million. A series of isotherms measured with the three noble gases and two different experimental setups led to this unprecedented level of uncertainty. These data are crucial for scientists in the field of gas metrology, working on pressure and temperature standards. Furthermore, with the new benchmark values for neon and argon, theoretical calculations, today about 3 orders of magnitude larger in uncertainty, can be checked and improved.