Beyond Explored Functionals: A Computational Journey of Two-Photon Absorption

We present a thorough investigation into the efficacy of 19 density functional theory (DFT) functionals, relative to RI-CC2 results, for computing two-photon absorption (2PA) cross sections (σ2PA) and key dipole moments (|μ00|, |μ11|, |Δμ|, |μ01|) for a series of coumarin dyes in the gas-phase. The functionals include different categories, including local density approximation (LDA), generalized gradient approximation (GGA), hybrid-GGA (H-GGA), range-separated hybrid-GGA (RSH-GGA), meta-GGA (M-GGA), and hybrid M-GGA (HM-GGA), with 14 of them being subjected to analysis for the first time with respect to predicting σ2PA values. Analysis reveals that functionals integrating both short-range (SR) and long-range (LR) corrections, particularly those within the RSH-GGA and HM-GGA classes, outperform the others. Furthermore, the range-separation approach was found more impactful compared to the varying percentages of Hartree-Fock exchange (HF Ex) within different functionals. The functionals traditionally recommended for 2PA do not appear among the top 9 in our study, which is particularly interesting, as these top-performing functionals have not been previously investigated in this context. This list is dominated by M11, QTP variants, ωB97X, ωB97X-V, and M06-2X, surpassing the performance of other functionals, including the commonly used CAM-B3LYP.

Ultralight aerogels via supramolecular polymerization of a new chiral perfluoropyridine-based sulfonimidamide organogelator

Chiral and enantiopure perfluorinated sulfonimidamides act as low-molecular weight gelators at low critical gelation concentration (<1 mg mL-1) via supramolecular polymerization in nonpolar organic solvents and more heterogenic mixtures, such as biodiesel and oil. Freeze-drying of the organogel leads to ultralight aerogel with extremely low density (1 mg mL-1). The gelation is driven by hydrogen bonding resulting in a helical molecular ordering and unique fibre assemblies as confirmed by scanning electron microscopy, CD spectroscopy, and computational modeling of the supramolecular structure.

Dynamical Effects of Solvation on Norbornadiene/Quadricyclane Systems

Molecules that can undergo reversible chemical transformations following the absorption of light, the so-called molecular photoswitches, have attracted increasing attention in technologies, such as solar energy storage. Here, the optical and thermochemical properties of the photoswitch are central to its applicability, and these properties are influenced significantly by solvation. We investigate the effects of solvation on two norbornadiene/quadricyclane photoswitches. Emphasis is put on the energy difference between the two isomers and the optical absorption as these are central to the application of the systems in solar energy storage. Using a combined classical molecular dynamics and quantum mechanical/molecular mechanical computational scheme, we showcase that the dynamic effects of solvation are important. In particular, it is found that standard implicit solvation models generally underestimate the energy difference between the two isomers and overestimate the strength of the absorption, while the explicit solvation spectra are also less red-shifted than those obtained using implicit solvation models. We also find that the absorption spectra of the two systems are strongly correlated with specific dihedral angles. Altogether, this highlights the importance of including the dynamic effects of solvation.

Variational Pair-Density Functional Theory: Dealing with Strong Correlation at the Protein Scale

Multiconfigurational pair-density functional theory (MC-PDFT) offers a promising solution to the challenges faced by traditional density functional theory (DFT) in addressing molecular systems containing transition metals, open-shells, or strong correlations in general. By utilizing both the density and on-top pair-density, MC-PDFT can make use of a more flexible multiconfigurational wave function to capture the necessary static correlation, while the pair-density functional also includes the effect of dynamic correlation. So far, MC-PDFT has been used after a multiconfigurational self-consistent field (MCSCF) step, using the orbitals and configuration interaction coefficients from the converged MCSCF wave function to compute PDFT energies and properties. Here, instead, we propose to perform a direct optimization of the wave function using the pair-density functionals, resulting in a variational formulation of MC-PDFT. We derive the expressions for the wave function gradient and illustrate their similarity to standard MCSCF equations. Furthermore, we illustrate the accuracy on a set of singlet-triplet gaps as well as dissociation curves. Our findings highlight one of MC-PDFT's standout features: a reduced dependency on the active space size compared to conventional multiconfigurational wave function methodologies. Additionally, we show that the computational cost of MC-PDFT is potentially lower than MCSCF and often on-par with standard Kohn-Sham DFT, which is demonstrated by performing a MC-PDFT calculation of the entire ferredoxin protein with 1447 atoms and nearly 12 000 basis functions.

Tensor Hypercontraction of Cluster Perturbation Theory: Quartic Scaling Perturbation Series for the Coupled Cluster Singles and Doubles Ground-State Energies

Even though cluster perturbation theory has been shown to be a robust noniterative alternative to coupled cluster theory, it is still plagued by high order polynomial computational scaling and the storage of higher order tensors. We present a proof-of-concept strategy for implementing a cluster perturbation theory ground-state energy series for the coupled cluster singles and doubles energy with N4 computational scaling using tensor hypercontraction (THC). The reduction in computational scaling by two orders is achieved by decomposing two electron repulsion integrals, doubles amplitudes and multipliers, as well as selected double intermediates to the THC format. Using the outlined strategy, we showcase that the THC pilot implementations retain numerical accuracy to within 1 kcal/mol relative to corresponding conventional and density fitting implementations, and we empirically verify the N4 scaling.

Toward Accurate Two-Photon Absorption Spectrum Simulations: Exploring the Landscape beyond the Generalized Gradient Approximation

In this Letter, we present a pioneering analysis of the density functional approximations (DFAs) beyond the generalized gradient approximation (GGA) for predicting two-photon absorption (2PA) strengths of a set of push-pull π-conjugated molecules. In more detail, we have employed a variety of meta-generalized gradient approximation (meta-GGA) functionals, including SCAN, MN15, and M06-2X, to assess their accuracy in describing the 2PA properties of a chosen set of 48 organic molecules. Analytic quadratic response theory is employed for these functionals, and their performance is compared against the previously studied DFAs and reference data obtained at the coupled-cluster CC2 level combined with the resolution-of-identity approximation (RI-CC2). A detailed analysis of the meta-GGA functional performance is provided, demonstrating that they improve upon their predecessors in capturing the key electronic features of the π-conjugated two-photon absorbers. In particular, the Minnesota functional MN15 shows very promising results as it delivers pleasingly accurate chemical rankings for two-photon transition strengths and excited-state dipole moments.

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.

Distinct Heterocyclic Moieties Govern the Selectivity of Thiophene-Vinylene-Based Ligands Towards Aβ or Tau Pathology in Alzheime's Disease

Distinct aggregated proteins are correlated with numerous neurodegenerative diseases and the development of ligands that selectively detect these pathological hallmarks is vital. Recently, the synthesis of thiophene-based optical ligands, denoted bi-thiophene-vinyl-benzothiazoles (bTVBTs), that could be utilized for selective assignment of tau pathology in brain tissue with Alzheime's disease (AD) pathology, was reported. Herein, we investigate the ability of these ligands to selectively distinguish tau deposits from aggregated amyloid-β (Aβ), the second AD associated pathological hallmark, when replacing the terminal thiophene moiety with other heterocyclic motifs. The selectivity for tau pathology was reduced when introducing specific heterocyclic motifs, verifying that specific molecular interactions between the ligands and the aggregates are necessary for selective detection of tau deposits. In addition, ligands having certain heterocyclic moieties attached to the central thiophene-vinylene building block displayed selectivity to aggregated Aβ pathology. Our findings provide chemical insights for the development of ligands that can distinguish between aggregated proteinaceous species consisting of different proteins and might also aid in creating novel agents for clinical imaging of tau pathology in AD.

Heterocyclic [9]Helicenes Exhibiting Bright Circularly Polarized Luminescence

We describe a concise synthetic strategy for the preparation of heterocyclic [9]helicenes and a simple preparative-scale protocol for the optical resolution of the resulting M- and P-enantiomers. The helicenes were characterized by single-crystal X-ray diffraction along with a range of spectroscopic and computational techniques. A fluorescence quantum yield of up to 65 % was observed, and the chiroptical properties of both M- and P-helicenes revealed large dissymmetry factors. The circularly polarized luminescence brightness reaches up to 17 M-1  cm-1 , as measured experimentally and verified computationally, which makes this the highest circularly polarized luminescence brightness among heterocyclic helicenes. We describe how chiroptical properties (both circular dichroism and circularly polarized luminescence) can be described and predicted using quantum chemical calculations. The synthetic approach also reveals by-products that originate from internal oxidation reactions, presumably mediated by the close proximity of the π-surfaces in the helicene structure.

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.

Algebraic Diagrammatic Construction Schemes Employing the Intermediate State Formalism: Theory, Capabilities, and Interpretation

Algebraic diagrammatic construction (ADC) schemes represent a family of ab initio methods for the calculation of excited electronic states and electron-detached and -attached states. All ADC methods have been demonstrated to possess great potential for molecular applications, e.g., for the calculation of absorption or photoelectron spectra or electron attachment processes. ADC originates from Green's function or propagator theory; however, most recent ADC developments heavily rely on the intermediate state representation or effective Liouvillian formalisms, which comprise new ADC methods and computational schemes for high-order properties. The different approaches for the calculation of excitation energies, ionization potentials, and electron affinities are intimately related, and they provide a coherent description of these quantities at equivalent levels of theory and with comparable errors. Most quantum chemical program packages contain ADC methods; however, the most complete ADC suite of methods can be found in the recent release of Q-Chem.

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. This paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

Self-Assembly of Chiro-Optical Materials from Nonchiral Oligothiophene-Porphyrin Derivatives and Random Coil Synthetic Peptides

Biomimetic chiral optoelectronic materials can be utilized in electronic devices, biosensors and artificial enzymes. Herein, this work reports the chiro-optical properties and architectural arrangement of optoelectronic materials generated from self-assembly of initially nonchiral oligothiophene-porphyrin derivatives and random coil synthetic peptides. The photo-physical- and structural properties of the materials were assessed by absorption-, fluorescence- and circular dichroism spectroscopy, as well as dynamic light scattering, scanning electron microscopy and theoretical calculations. The materials display a three-dimensional ordered helical structure and optical activity that are observed due to an induced chirality of the optoelectronic element upon interaction with the peptide. Both these properties are influenced by the chemical composition of the oligothiophene-porphyrin derivative, as well as the peptide sequence. We foresee that our findings will aid in developing self-assembled optoelectronic materials with dynamic architectonical accuracies, as well as offer the possibility to generate the next generation of materials for a variety of bioelectronic applications.

Modeling Absorption and Emission Spectroscopies of Symmetric and Asymmetric Azaoxahelicenes in Vacuum and Solution

Helicenes are of general interest due to the significant chiral signals in both absorption- and emission-based spectroscopy. Herein, the spectroscopic properties of four recently synthesized azaoxahelicenes are studied using density functional theory methods. The azaoxahelicenes have 7, 9, 10, and 13 units and one to two complete turns of the structure. UV-vis absorption and electronic circular dichroism spectra are determined both in vacuum and in solution using explicit solvation through a combined molecular dynamics/polarizable embedding framework. Additionally, emission and circularly polarized luminescence spectra are determined based on vibronic calculations. The resulting spectra are in good agreement with the experimentally available data, highlighting that both absorption- and emission-based spectra of the systems can be modeled computationally such that reliable predictions can be made for systems that are yet to be synthesized.

Fast Overlapping Block Processing Algorithm for Feature Extraction

In many video and image processing applications, the frames are partitioned into blocks, which are extracted and processed sequentially. In this paper, we propose a fast algorithm for calculation of features of overlapping image blocks. We assume the features are projections of the block on separable 2D basis functions (usually orthogonal polynomials) where we benefit from the symmetry with respect to spatial variables. The main idea is based on a construction of auxiliary matrices that virtually extends the original image and makes it possible to avoid a time-consuming computation in loops. These matrices can be pre-calculated, stored and used repeatedly since they are independent of the image itself. We validated experimentally that the speed up of the proposed method compared with traditional approaches approximately reaches up to 20 times depending on the block parameters.

Analytical gradients for core-excited states in the algebraic diagrammatic construction (ADC) framework

Expressions for analytical molecular gradients of core-excited states have been derived and implemented for the hierarchy of algebraic diagrammatic construction (ADC) methods up to extended second-order within the core-valence separation (CVS) approximation. We illustrate the use of CVS-ADC gradients by determining relaxed core-excited state potential energy surfaces and optimized geometries for water, formic acid, and benzene. For water, our results show that in the dissociative lowest core-excited state, a linear configuration is preferred. For formic acid, we find that the O K-edge lowest core-excited state is non-planar, a fact that is not captured by the equivalent core approximation where the core-excited atom with its hole is replaced by the "Z + 1" neighboring atom in the periodic table. For benzene, the core-excited state gradients are presented along the Jahn-Teller distorted geometry of the 1s → π* excited state. Our development may pave a new path to studying the dynamics of molecules in their core-excited states.

Efficient Open-Source Implementations of Linear-Scaling Polarizable Embedding: Use Octrees to Save the Trees

We present open-source implementations of the linear-scaling fast multipole method (FMM) within the polarizable embedding (PE) model for efficient treatment of large polarizable environments in computational spectroscopy simulations. The implementations are tested for accuracy, efficiency, and usability on model systems as well as more realistic biomolecular systems. We explain how FMM parameters affect the calculation of molecular properties and show that PE calculations employing FMM can be carried out in a black-box manner. The efficiency of the linear-scaling approach is demonstrated by simulating the UV/vis spectrum of a chromophore in an environment of more than 1 million polarizable sites. Our implementations are interfaced to several open-source quantum chemistry programs, making computational spectroscopy simulations within the PE model and FMM available to a large variety of methods and a broad user base.

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 C₆ dispersion coefficients of carbon fullerenes are well established, the widespread reported scalings for the latter (ranging from N^2.2 to N^2.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 C₆ coefficients show scalings of N^1.2 and N^2.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.

Efficient implementation of isotropic cubic response functions for two-photon absorption cross sections within the self-consistent field approximation

Within the self-consistent field approximation, computationally tractable expressions for the isotropic second-order hyperpolarizability have been derived and implemented for the calculation of two-photon absorption cross sections. The novel tensor average formulation presented in this work allows for the evaluation of isotropic damped cubic response functions using only ∼3.3% (one-photon off-resonance regions) and ∼10% (one-photon resonance regions) of the number of auxiliary Fock matrices required when explicitly calculating all the needed individual tensor components. Numerical examples of the two-photon absorption cross section in the one-photon off-resonance and resonance regions are provided for alanine-tryptophan and 2,5-dibromo-1,4-bis(2-(4-diphenylaminophenyl)vinyl)-benzene. Furthermore, a benchmark set of 22 additional small- and medium-sized organic molecules is considered. In all these calculations, a quantitative assessment is made of the reduced and approximate forms of the cubic response function in the one-photon off-resonance regions and results demonstrate a relative error of less than ∼5% when using the reduced expression as compared to the full form of the isotropic cubic response function.

Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems

The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.

Psi4 1.4: Open-source software for high-throughput quantum chemistry

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.