Adiabatic connection interaction strength interpolation method made accurate for the uniform electron gas.

The adiabatic connection interaction strength interpolation (ISI)-like method provides a high-level expression for the correlation energy, being, in principle, exact not only in the weak-interaction limit, where it recovers the second-order Görling-Levy perturbation term, but also in the strong-interaction limit that is described by the strictly correlated electron approach. In this work, we construct a genISI functional made accurate for the uniform electron gas, a solid-state physics paradigm that is a very difficult test for ISI-like correlation functionals. We assess the genISI functional for various jellium spheres with the number of electrons Z ≤ 912 and for the non-relativistic noble atoms with Z ≤ 290.

Regularized and Opposite Spin-Scaled Functionals from Møller-Plesset Adiabatic Connection─Higher Accuracy at Lower Cost.

Noncovalent interactions (NCIs) play a crucial role in biology, chemistry, material science, and everything in between. To improve pure quantum-chemical simulations of NCIs, we propose a methodology for constructing approximate correlation energies by combining an interpolation along the Møller-Plesset adiabatic connection (MP AC) with a regularization and spin-scaling strategy applied to MP2 correlation energies. This combination yields κ-SPL2, which exhibits superior accuracy for NCIs compared to any of the individual strategies.

TURBOMOLE: Today and Tomorrow.

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches.

Minimal Auxiliary Basis Set Approach for the Electronic Excitation Spectra of Organic Molecules.

We report a minimal auxiliary basis model for time-dependent density functional theory (TDDFT) with hybrid density functionals that can accurately reproduce excitation energies and absorption spectra from TDDFT while reducing cost by about 2 orders of magnitude. Our method, dubbed TDDFT-ris, employs the resolution-of-the-identity technique with just one -type auxiliary basis function per atom for the linear response operator, where the Gaussian exponents are parametrized across the periodic table using tabulated atomic radii with a single global scaling factor. By tuning on a small test set, we determine a single functional-independent scale factor that balances errors in excitation energies and absorption spectra.

Orbital-free methods for plasmonics: Linear response.

Plasmonic systems, such as metal nanoparticles, are widely used in different areas of application, going from biology to photovoltaics. The modeling of the optical response of such systems is of fundamental importance to analyze their behavior and to design new systems with required properties. When the characteristic sizes/distances reach a few nanometers, nonlocal and spill-out effects become relevant and conventional classical electrodynamics models are no more appropriate.

Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit.

Adiabatic connection models (ACMs), which interpolate between the limits of weak and strong interaction, are powerful tools to build accurate exchange-correlation functionals. If the exact weak-interaction expansion from the second-order perturbation theory is included, a self-consistent implementation of these functionals is challenging and still absent in the literature. In this work, we fill this gap by presenting a fully self-consistent-field (SCF) implementation of some popular ACM functionals.

Plasmon Couplings from Subsystem Time-Dependent Density Functional Theory.

Many applications in plasmonics are related to the coupling between metallic nanoparticles (MNPs) or between an emitter and a MNP. The theoretical analysis of such a coupling is thus of fundamental importance to analyze the plasmonic behavior and to design new systems. While classical methods neglect quantum and spill-out effects, time-dependent density functional theory (TD-DFT) considers all of them and with Kohn-Sham orbitals delocalized over the whole system.

Noncovalent Interactions from Models for the Møller-Plesset Adiabatic Connection.

Given the omnipresence of noncovalent interactions (NCIs), their accurate simulations are of crucial importance across various scientific disciplines. Here we construct accurate models for the description of NCIs by an interpolation along the Møller-Plesset adiabatic connection (MP AC). Our interpolation approximates the correlation energy, by recovering MP2 at small coupling strengths and the correct large-coupling strength expansion of the MP AC, recently shown to be a functional of the Hartree-Fock density.

Minimal auxiliary basis set for time-dependent density functional theory and comparison with tight-binding approximations: Application to silver nanoparticles.

The modeling of optical spectra of plasmonic nanoparticles via first-principles approaches is computationally expensive; thus, methods with high accuracy/computational cost ratio are required. Here, we show that the Time-Dependent Density Functional Theory (TDDFT) approach can be strongly simplified if only one s-type function per atom is employed in the auxiliary basis set, with a properly optimized exponent. This approach (named TDDFT-as, for auxiliary s-type) predicts excitation energies for silver nanoparticles with different sizes and shapes with an average error of only 12 meV compared to reference TDDFT calculations.

Performance of Semilocal Kinetic Energy Functionals for Orbital-Free Density Functional Theory.

We assess several generalized gradient approximations (GGAs) and Laplacian-level meta-GGAs (LL-MGGA) kinetic energy (KE) functionals for orbital-free density functional theory calculations of bulk metals and semiconductors, considering equilibrium distances, bulk moduli, total and kinetic energies, and the electron densities. We also considered the effects of the pseudopotentials, the vacancy formation energies, and the bond lengths of molecular dimers. We found that LL-MGGA KE functionals are distinctively superior to GGA functionals, showing the importance of the Laplacian of the density in the functional construction.

Investigation of the Exchange-Correlation Potentials of Functionals Based on the Adiabatic Connection Interpolation.

We have studied the correlation potentials produced by various adiabatic connection models (ACMs) for several atoms and molecules. The results have been compared to accurate reference potentials (coupled cluster and quantum Monte Carlo results) as well as to state-of-the-art ab initio DFT approaches. We have found that all the ACMs yield correlation potentials that exhibit a correct behavior, quite resembling scaled second-order Görling-Levy (GL2) potentials and including most of the physically meaningful features of the accurate reference data.

Optical properties of plasmonic core-shell nanomatryoshkas: a quantum hydrodynamic analysis.

Plasmonic response of the metallic structure characterized by sub-nanometer dielectric gaps can be strongly affected by nonlocal or quantum effects. In this paper, we investigate these effects in spherical Na and Au nanomatryoshka structures with sub-nanometer core-shell separation. We use the state-of-the-art quantum hydrodynamic theory (QHT) to study both near-field and far-field optical properties of these systems: results are compared with the classical local response approximation (LRA), Thomas-Fermi hydrodynamic theory (TF-HT), and the reference time-dependent density functional theory (TD-DFT).

Restoring Size Consistency of Approximate Functionals Constructed from the Adiabatic Connection.

Approximate exchange-correlation functionals built by modeling in a nonlinear way the adiabatic connection (AC) integrand of density functional theory have many attractive features, being virtually parameter-free and satisfying different exact properties, but they also have a fundamental flaw: they violate the size-consistency condition, crucial to evaluate interaction energies of molecular systems. We show that size consistency in the AC-based functionals can be restored in a very simple way at no extra computational cost. Results on a large set of benchmark molecular interaction energies show that functionals based on the interaction strength interpolation approximations are significantly more accurate than second-order perturbation theory.

Assessment of interaction-strength interpolation formulas for gold and silver clusters.

The performance of functionals based on the idea of interpolating between the weak- and the strong-interaction limits the global adiabatic-connection integrand is carefully studied for the challenging case of noble-metal clusters. Different interpolation formulas are considered and various features of this approach are analyzed. It is found that these functionals, when used as a correlation correction to Hartree-Fock, are quite robust for the description of atomization energies, while performing less well for ionization potentials.

Modified Fourth-Order Kinetic Energy Gradient Expansion with Hartree Potential-Dependent Coefficients.

Using the semiclassical neutral atom theory, we developed a modified fourth-order kinetic energy (KE) gradient expansion (GE4m) that keeps unchanged all the linear-response terms of the uniform electron gas and gives a significant improvement with respect to the known semilocal functionals for both large atoms and jellium surfaces. On the other hand, GE4m is not accurate for light atoms; thus, we modified the GE4m coefficients making them dependent on a novel ingredient, the reduced Hartree potential, recently introduced in the Journal of Chemical Physics 2016, 145, 084110, in the context of exchange functionals. The resulting KE gradient expansion functional, named uGE4m, belongs to the novel class of u-meta-generalized-gradient-approximations (uMGGA) whose members depend on the conventional ingredients (i.

Electrostatic Mechanophores in Tuneable Light-Emitting Piezopolymer Nanowires.

Electromechanical coupling through piezoelectric polymer chains allows the emission of organic molecules in active nanowires to be tuned. This effect is evidenced by highly bendable arrays of counter-ion dye-doped nanowires made of a poly(vinylidenefluoride) copolymer. A reversible redshift of the dye emission is found upon the application of dynamic stress during highly accurate bending experiments.

Laplacian-dependent models of the kinetic energy density: Applications in subsystem density functional theory with meta-generalized gradient approximation functionals.

The development of semilocal models for the kinetic energy density (KED) is an important topic in density functional theory (DFT). This is especially true for subsystem DFT, where these models are necessary to construct the required non-additive embedding contributions. In particular, these models can also be efficiently employed to replace the exact KED in meta-Generalized Gradient Approximation (meta-GGA) exchange-correlation functionals allowing to extend the subsystem DFT applicability to the meta-GGA level of theory.

Interaction-Strength Interpolation Method for Main-Group Chemistry: Benchmarking, Limitations, and Perspectives.

We have tested the original interaction-strength-interpolation (ISI) exchange-correlation functional for main group chemistry. The ISI functional is based on an interpolation between the weak and strong coupling limits and includes exact-exchange as well as the Görling-Levy second-order energy. We have analyzed in detail the basis-set dependence of the ISI functional, its dependence on the ground-state orbitals, and the influence of the size-consistency problem.

Hartree potential dependent exchange functional.

We introduce a novel non-local ingredient for the construction of exchange density functionals: the reduced Hartree parameter, which is invariant under the uniform scaling of the density and represents the exact exchange enhancement factor for one- and two-electron systems. The reduced Hartree parameter is used together with the conventional meta-generalized gradient approximation (meta-GGA) semilocal ingredients (i.e.

Accurate Kohn-Sham ionization potentials from scaled-opposite-spin second-order optimized effective potential methods.

One important property of Kohn-Sham (KS) density functional theory is the exact equality of the energy of the highest occupied KS orbital (HOMO) with the negative ionization potential of the system. This exact feature is out of reach for standard density-dependent semilocal functionals. Conversely, accurate results can be obtained using orbital-dependent functionals in the optimized effective potential (OEP) approach.