Complex Polarization Propagator Research Articles

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.

Complex Linear Response Functions for a Multiconfigurational Self-Consistent Field Wave Function in a High Performance Computing Environment.

We present novel developments for the highly efficient evaluation of complex linear response functions of a multiconfigurational self-consistent field (MCSCF) wave function as implemented in MultiPsi. Specifically, expressions for the direct evaluation of linear response properties at given frequencies using the complex polarization propagator (CPP) approach have been implemented, within both the Tamm-Dancoff approximation (TDA) and the random phase approximation (RPA). Purely real algebra with symmetric and antisymmetric trial vectors in a shared subspace is used wherein the linear response equations are solved. Two bottlenecks of large scale MC-CPP calculations, namely, the memory footprint and computational time, are addressed. The former is addressed by limiting the size of the subspace of trial vectors by using singular value decomposition (SVD) on either orbital or CI subspaces. The latter is addressed using an efficient parallel implementation as well as the strategy of dynamically adding linear response equations at near-convergence to neighboring roots. Furthermore, a novel methodology for decomposing MC-CPP spectra in terms of intuitive orbital excitations in an approximate fashion is presented. The performance of the code is illustrated with several numerical examples, including the X-ray spectrum of a molecule with nearly one hundred atoms. Additionally, for X-ray spectroscopy, the effect of including or excluding the core orbital in the active space on small covalent metal complexes is 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(H2O)6]3+ and [Ir(H2O)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.

Size-dependent polarizabilities and van der Waals dispersion coefficients of fullerenes from large-scale complex polarization propagator calculations.

While the anomalous non-additive size-dependencies of static dipole polarizabilities and van der Waals C6 dispersion coefficients of carbon fullerenes are well established, the widespread reported scalings for the latter (ranging from N2.2 to N2.8) call for a comprehensive first-principles investigation. With a highly efficient implementation of the linear complex polarization propagator, we have performed Hartree-Fock and Kohn-Sham density functional theory calculations of the frequency-dependent polarizabilities for fullerenes consisting of up to 540 carbon atoms. Our results for the static polarizabilities and C6 coefficients show scalings of N1.2 and N2.2, respectively, thereby deviating significantly from the previously reported values obtained with the use of semi-classical/empirical methods. Arguably, our reported values are the most accurate to date as they represent the first ab initio or first-principles treatment of fullerenes up to a convincing system size.

Characterization of tetrakis(thiadiazole)porphyrazine metal complexes by magnetic circular dichroism and magnetic circularly polarized luminescence.

The magnetic circular dichroism (MCD) spectra of metal complexes of tetrakis(thiadiazole)porphyrazines ([TTDPzM] with M = 2HI , ZnII , MgII (H2 O), and CdII ) have been recorded in dimethyl formamide solution. Together with the UV-Vis spectra, the MCD spectra provide useful information about the structure and electronic properties of the complexes. The experimental UV-Vis and MCD spectra compare pretty well with DFT calculations of two sorts, based either on the sum-over-states (SOS) approach or on the complex polarization propagator approach. They further corroborate the findings and interpretation of MCD spectra of porphyrazines based on the model of Michl for peripheral molecular orbitals. Magnetic circularly polarized luminescence (MCPL) spectra, quite uncommon in the literature, have been recorded for [TTDPzM] (M = 2HI , ZnII , MgII (H2 O)).

Dynamic Polarizability and Higher-Order Electric Properties of Fluorene, Carbazole, and Dibenzofuran

Static electric properties, from the dipole moment to the second-hyperpolarizability tensor γ, of the 3-membered, isoelectronic ring molecules, fluorene (FL), carbazole (CR), and dibenzofuran (DBF), have been calculated at various levels of approximation. The electron correlation effects have been included at the coupled-cluster (CC) level, using CCSD and CC2 versions of the method. DFT calculations with the CAM-B3LYP functional have also been performed, and the results are compared to the CC values. The electric property-tailored Pol basis set and its more compact Z3Pol version have been employed in all static calculations. Differences between dipole polarizability values computed at the Pol and Z3Pol bases have been found to be almost negligible. Therefore, all components of the frequency-dependent dipole polarizability tensor α(-ω;ω) have been determined at the CAM-B3LYP/Z3Pol level. Divergence occurring at electronic resonances has been eliminated using the complex polarization propagator (CPP) formalism, explicitly introducing an imaginary iΓ parameter to account approximately for the finite lifetime of the excited state. The imaginary part of the dipole polarizability Im α(-ω;ω) has been calculated for a wide range of external radiation energies up to 10 eV, and its maxima have been compared to the calculated vertical electronic excitation energies.

X-ray absorption spectrum simulations of hexagonal ice.

We calibrate basis sets and performance of two theoretical approaches to compute X-ray absorption spectra (XAS) of condensed water by comparison to experiments on hexagonal ice Ih. We apply both the transition-potential half-core-hole approach and the complex polarization propagator using four different models of the crystal with increasing oxygen and proton disorder but find poor agreement with experiments. We note that there are large variations in experimental spectra depending on detection mode and how the ice samples were prepared, which leads us to critically investigate what structures were actually prepared and measured in each case. This is done by using a Monte Carlo-based fitting technique which fits the spectra based on a library of precomputed spectra and assigns weights to contributions from different model structures. These are then used to generate O-O and O-H radial distribution functions and tetrahedrality parameters associated with each of the measured spectra. We find that all spectra are associated with sharp peaks at the oxygen positions in the perfect lattice, but with significant disorder around these positions. We suggest that presently available XAS of hexagonal ice are not fully representative of the perfect crystalline lattice, but contain varying amounts of defects and possible contributions from low-density amorphous ice.

One-Photon Absorption Properties from a Hybrid Polarizable Density Embedding/Complex Polarization Propagator Approach for Polarizable Solutions.

We present a formulation of the polarizable density embedding (PDE) method in combination with the complex polarization propagator (CPP) method for the calculation of absorption spectra of molecules in solutions. The method is particularly useful for the calculation of near-edge X-ray absorption fine structure (NEXAFS) spectra. We compare the performance of PDE-CPP with the previously formulated polarizable embedding (PE)-CPP model for the calculation of the NEXAFS spectra of adenine, formamide, glycine, and adenosine triphosphate (ATP) in water at the carbon and nitrogen K-edges, as well as of formamide and glycine at the oxygen K-edge. In general, we find only minor differences between the performance of PDE and PE for the targeted parts of the spectra, except in the case of transitions involving Rydberg states, for which nonelectrostatic effects are found to be important.

Resonant Inelastic X-ray Scattering Amplitudes and Cross Sections in the Algebraic Diagrammatic Construction/Intermediate State Representation (ADC/ISR) Approach.

A scheme has been derived and implemented to gain computational access to the full electronic part of the Kramers-Heisenberg-Dirac (KHD) expression for resonant and nonresonant inelastic scattering amplitudes. Our implementation of this scheme is based on the complex polarization propagator in the algebraic diagrammatic construction (ADC) framework and within its intermediate state representation (ISR). The hierarchy of the second- and third-order ADC/ISR computational schemes known as ADC(2), ADC(2)-x, and ADC(3/2) is considered, and the calculated resonant inelastic X-ray scattering (RIXS) amplitudes and transition strengths for water are in excellent agreement with recent experimental data.

Model Order Reduction Algorithm for Estimating the Absorption Spectrum.

The ab initio description of the spectral interior of the absorption spectrum poses both a theoretical and computational challenge for modern electronic structure theory. Due to the often spectrally dense character of this domain in the quantum propagator's eigenspectrum for medium-to-large sized systems, traditional approaches based on the partial diagonalization of the propagator often encounter oscillatory and stagnating convergence. Electronic structure methods which solve the molecular response problem through the solution of spectrally shifted linear systems, such as the complex polarization propagator, offer an alternative approach which is agnostic to the underlying spectral density or domain location. This generality comes at a seemingly high computational cost associated with solving a large linear system for each spectral shift in some discretization of the spectral domain of interest. In this work, we present a novel, adaptive solution to this high computational overhead based on model order reduction techniques via interpolation. Model order reduction reduces the computational complexity of mathematical models and is ubiquitous in the simulation of dynamical systems and control theory. The efficiency and effectiveness of the proposed algorithm in the ab initio prediction of X-ray absorption spectra is demonstrated using a test set of challenging water clusters which are spectrally dense in the neighborhood of the oxygen K-edge. On the basis of a single, user defined tolerance we automatically determine the order of the reduced models and approximate the absorption spectrum up to the given tolerance. We also illustrate that, for the systems studied, the automatically determined model order increases logarithmically with the problem dimension, compared to a linear increase of the number of eigenvalues within the energy window. Furthermore, we observed that the computational cost of the proposed algorithm only scales quadratically with respect to the problem dimension.

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

An implementation of the damped linear response function, or complex polarization propagator, using the algebraic-diagrammatic construction (ADC) scheme has been developed and utilized for the calculation of electric-dipole polarizabilities and C6 dispersion coefficients. Four noble gases (He, Ne, Ar, and Kr), five n-alkanes (methane, ethane, propane, butane, and pentane), three carbonyls (formaldehyde, acetaldehyde, and acetone), and three unsaturated hydrocarbons (ethene, acetylene, and benzene) have been treated with the hierarchical set of models ADC(2), ADC(2)-x, and ADC(3/2), and comparison has been made to results obtained with damped linear response Hartree–Fock (HF) and coupled cluster singles and doubles (CCSD) theory as well as high-quality experimental estimates via the dipole oscillator strength distribution approach. This study marks the first ADC calculations of C6 dispersion coefficients and the first ADC(3/2) calculations of static polarizabilities. Results at CCSD and ADC(3/2) levels of theory are shown to be of similar quality, with electron correlation effects increasing the molecular property values for all calculations except CCSD considerations of ethene and acetylene (attributed to an overestimation of bond electron density at HF level of theory). The discrepancies between CCSD and ADC(3/2) are partially due to ADC overestimating anisotropies, and discrepancies with respect to experimental values are partially due to the lack of zero-point vibrational effects in the present study.

A Polarization Propagator for Nonlinear X-ray Spectroscopies.

A complex polarization propagator approach has been developed to third order and implemented in density functional theory (DFT), allowing for the direct calculation of nonlinear molecular properties in the X-ray wavelength regime without explicitly addressing the excited-state manifold. We demonstrate the utility of this propagator method for the modeling of coherent near-edge X-ray two-photon absorption using, as an example, DFT as the underlying electronic structure model. Results are compared with the corresponding near-edge X-ray absorption fine structure spectra, illuminating the differences in the role of symmetry, localization, and correlation between the two spectroscopies. The ramifications of this new technique for nonlinear X-ray research are briefly discussed.

Surface-Enhanced Raman Scattering Due to Charge-Transfer Resonances: A Time-Dependent Density Functional Theory Study of Ag13-4-Mercaptopyridine

We have used time-dependent density functional theory in conjunction with the CAM-B3LYP functional and MWB28/aug-cc-pVDZ basis set to determine non-, near-, and on-resonance Raman spectra for a complex formed by 4-mercaptopyridine (4-Mpy) binding with a Ag13 cluster via the thiolate Ag–S bond. Geometry optimizations of the Ag13-4-Mpy complex showed an on-top structure directly bound to one Ag atom with the ring of the molecule almost flat with respect to two Ag atoms of the complex. The corresponding B3LYP/MWB28/aug-cc-pVDZ geometry is also an on-top structure directly bound to one Ag atom, but the molecule is directed away from the surface. The near-resonance Raman calculations were carried out in the infinite lifetime approximation, while the on-resonant Raman excitation profiles were calculated with the complex polarization propagator (CPP) approach, introducing a half width at half-maximum spectral broadening of 0.2 eV. Calculation of the UV–vis spectra of the isolated 4-Mpy and of the Ag13-4-Mpy comp...

Four-Component Damped Density Functional Response Theory Study of UV/Vis Absorption Spectra and Phosphorescence Parameters of Group 12 Metal-Substituted Porphyrins.

The influences of group 12 (Zn, Cd, Hg) metal-substitution on the valence spectra and phosphorescence parameters of porphyrins (P) have been investigated in a relativistic setting. In order to obtain valence spectra, this study reports the first application of the damped linear response function, or complex polarization propagator, in the four-component density functional theory framework [as formulated in Villaume et al. J. Chem. Phys. 2010 , 133 , 064105 ]. It is shown that the steep increase in the density of states as due to the inclusion of spin-orbit coupling yields only minor changes in overall computational costs involved with the solution of the set of linear response equations. Comparing single-frequency to multifrequency spectral calculations, it is noted that the number of iterations in the iterative linear equation solver per frequency grid-point decreases monotonously from 30 to 0.74 as the number of frequency points goes from one to 19. The main heavy-atom effect on the UV/vis-absorption spectra is indirect and attributed to the change of point group symmetry due to metal-substitution, and it is noted that substitutions using heavier atoms yield small red-shifts of the intense Soret-band. Concerning phosphorescence parameters, the adoption of a four-component relativistic setting enables the calculation of such properties at a linear order of response theory, and any higher-order response functions do not need to be considered-a real, conventional, form of linear response theory has been used for the calculation of these parameters. For the substituted porphyrins, electronic coupling between the lowest triplet states is strong and results in theoretical estimates of lifetimes that are sensitive to the wave function and electron density parametrization. With this in mind, we report our best estimates of the phosphorescence lifetimes to be 460, 13.8, 11.2, and 0.00155 s for H2P, ZnP, CdP, and HgP, respectively, with the corresponding transition energies being equal to 1.46, 1.50, 1.38, and 0.89 eV.

4-Component relativistic calculations of L3 ionization and excitations for the isoelectronic species UO2(2+), OUN(+) and UN2.

We present a 4-component relativistic study of uranium 2p3/2 ionization and excitation in the isoelectronic series UO2(2+), OUN(+) and UN2. We calculate ionization energies by ΔSCF at the Hartree-Fock (HF) and Kohn-Sham (KS) level of theory. At the ΔHF level we observe a perfectly linear chemical shift of ionization energies with respect to uranium atomic charges obtained from projection analysis. We have also developed a non-canonical 2nd-order Møller-Plesset code for wave function based correlation studies. We observe the well-known failure of Koopmans' theorem for core ionization due to the dominance of orbital relaxation over electron correlation effects. More unexpectedly, we find that the correlation contribution has the same sign as the relaxation contribution and show that this is due to a strong coupling of relaxation and correlation. We simulate uranium L3 XANES spectra, dominated by 2p3/2 → U6d transitions, by restricted excitation window time-dependent density functional theory (REW-TDDFT) and the complex polarization propagator (CPP) approach and demonstrate that they give identical spectra when the same Lorentz broadening is chosen. We also simulate XANES spectra by the Hartree-Fock based static exchange (STEX) method and show how STEX excitation energies can be reproduced by time-dependent Hartree-Fock calculations within the Tamm-Dancoff approximation. We furthermore show that Koopmans' theorem provide a correct approximation of ionization energies in the linear response regime and use this observation to align REW-TDDFT and CPP spectra with STEX ones. We point out that the STEX method affords the most detailed assignment of spectra since it employs virtual orbitals optimized for the selected core ionization. The calculated XANES spectra reflect the loss of bound virtual orbitals as the molecular charge is reduced along the isoelectronic series.

Beyond the electric-dipole approximation: A formulation and implementation of molecular response theory for the description of absorption of electromagnetic field radiation.

We present a formulation of molecular response theory for the description of a quantum mechanical molecular system in the presence of a weak, monochromatic, linearly polarized electromagnetic field without introducing truncated multipolar expansions. The presentation focuses on a description of linear absorption by adopting the energy-loss approach in combination with the complex polarization propagator formulation of response theory. Going beyond the electric-dipole approximation is essential whenever studying electric-dipole-forbidden transitions, and in general, non-dipolar effects become increasingly important when addressing spectroscopies involving higher-energy photons. These two aspects are examined by our study of the near K-edge X-ray absorption fine structure of the alkaline earth metals (Mg, Ca, Sr, Ba, and Ra) as well as the trans-polyenes. In following the series of alkaline earth metals, the sizes of non-dipolar effects are probed with respect to increasing photon energies and a detailed assessment of results is made in terms of studying the pertinent transition electron densities and in particular their spatial extension in comparison with the photon wavelength. Along the series of trans-polyenes, the sizes of non-dipolar effects are probed for X-ray spectroscopies on organic molecules with respect to the spatial extension of the chromophore.

Enhancement of Internal Motions of Lysozyme through Interaction with Gold Nanoclusters and its Optical Imaging

Understanding the interaction of gold nanoclusters with proteins has important ramifications in various fields. We present a study of the interaction between gold nanoclusters and lysozyme investigated using classical molecular dynamics and center-of-mass pulling simulations. The results reveal that the gold nanoclusters induce significant structural changes in lysozyme. Because the internal motions of lysozyme are related to its function, the changes in these internal motions have been quantified using principal component analysis of the molecular dynamics trajectories. The internal motions of lysozyme that are important for its function have been altered because of the interaction with the gold nanocluster. We have also explored how these induced changes in the lysozyme structure affect specific optical properties of the gold nanocluster using the complex polarization propagator method within the time-dependent density functional theory framework, which is of relevance for studies of the optical imaging...

On the interplay between chirality and exciton coupling: a DFT calculation of the circular dichroism in π-stacked ethylene.

The chirality of stacked weakly interacting π-systems was interpreted in terms of Frenkel exciton states and the formation of excitonic circular dichroism (CD) bands was monitored for ethylene stacks of varying sizes. Convergence of CD bands with respect to the system size was observed for stacks involving around 10 molecules. By means of rotation around the C-C double bond in ethylene, chirality was induced in the monomeric system and which was shown to dominate the spectral responses, even for polymer aggregates. In helical assemblies of chiral entities, there will always be a mix of excitonic and monomeric contributions to the CD signal and it is demonstrated that the complex polarization propagator approach in combination with Density Functional Theory is a suitable method to address this situation.

Efficient Calculations of Molecular Linear Response Properties for Spectral Regions.

Molecular spectra can be determined from molecular response functions, by solving the so-called damped response equations using the complex polarization propagator approach. The overall structure of response equations is identical for variational wave functions such as the Hartree-Fock, multiconfiguration self-consistent field, and Kohn-Sham density functional theory, and the key program module is the linear response equation solver. We present an implementation of the solver using the algorithm with symmetrized vectors, optimized for addressing spectral regions of a width of some 5-10 eV and a resolution below 0.1 eV. The work is illustrated by the consideration of UV-vis as well as near carbon K -edge absorption spectra of the C60 fullerene. We demonstrate that it is possible to converge tightly response equations for hundreds of optical frequencies in resonance regions of the spectrum at a cost not much exceeding the solution of a single response equation in the nonresonant region. Our work is implemented in the molecular orbital based module of the Dalton program and serves as a documentation of the code distributed in the Dalton2013 release version.

Nuclear spin circular dichroism

Recent years have witnessed a growing interest in magneto-optic spectroscopy techniques that use nuclear magnetization as the source of the magnetic field. Here we present a formulation of magnetic circular dichroism (CD) due to magnetically polarized nuclei, nuclear spin-induced CD (NSCD), in molecules. The NSCD ellipticity and nuclear spin-induced optical rotation (NSOR) angle correspond to the real and imaginary parts, respectively, of (complex) quadratic response functions involving the dynamic second-order interaction of the electron system with the linearly polarized light beam, as well as the static magnetic hyperfine interaction. Using the complex polarization propagator framework, NSCD and NSOR signals are obtained at frequencies in the vicinity of optical excitations. Hartree-Fock and density-functional theory calculations on relatively small model systems, ethene, benzene, and 1,4-benzoquinone, demonstrate the feasibility of the method for obtaining relatively strong nuclear spin-induced ellipticity and optical rotation signals. Comparison of the proton and carbon-13 signals of ethanol reveals that these resonant phenomena facilitate chemical resolution between non-equivalent nuclei in magneto-optic spectra.