Fractional Occupation Numbers Research Articles

Slope of the Delocalization Function Is Proportional to Analytical Hardness.

Conceptual Density Functional Theory (CDFT) has been extended beyond its traditional role in elucidating chemical reactivity to the development of density functional theory methods, e.g., the investigation of the delocalization error. This delocalization error causes the dependence of the energy on the number of electrons (N) to deviate from its exact piecewise linear behavior, an error which is the basis of many well-known limitations of commonly used density-functional approximations (DFAs). Following our previous work on the analytical hardness η± for pure functionals, we extend its application to hybrid and range-separated functionals. A comparison is made between the analytical hardness and the slope of the delocalization function introduced by Hait and Head-Gordon. Our results show that there is a linear relationship between its slope and the analytical hardness. An approximate scheme is presented to construct the energy vs N curve without fractional occupation number calculations. The extension to densities is discussed.

Multireference calculations on bond dissociation and biradical polycyclic aromatic hydrocarbons as guidance for fractional occupation number weighted density analysis in DFT calculations

This study explores open-shell biradical and polyradical molecular compounds based on extended multireference (MR) methods (MR-configuration interaction with singles and doubles (CISD) and MR-averaged quadratic coupled cluster (AQCC) approach) using the numbers of unpaired densities NU. These results were used to guide the analysis of the fractional occupation number weighted density (FOD) calculated within the finite temperature (FT) density functional theory (DFT) approach. As critical test examples, the dissociation of carbon–carbon (CC) single, double and triple bonds and a benchmark set of polycyclic aromatic hydrocarbons (PAHs) have been chosen. By examining single, double, and triple bond dissociations, we demonstrate the utility and accuracy but also limitations of the FOD analysis for describing these dissociation processes. In significant extension of previous work (Phys Chem Chem Phys 25: 27380–27393), the assessment of FOD applications for different classes of DFT functionals was performed examining the range-separated functionals ωB97XD, ωB97M-V, CAM-B3LYP, LC-ωPBE, and MN12-SX, the hybrid (M06-2X) functional and the double hybrid (B2P-LYP) functional. In all cases, strong correlations between NFOD and NU values are found. The major task was to develop a new linear regression formula for range-separated functionals allowing a convenient determination of the optimal electronic temperature Tel for the FT-DFT calculation. We also established an optimal temperature for the semiempirical extended tight-binding GFN2-xTB method. These findings significantly broaden the applicability of FOD analysis across various DFT functionals and semiempirical methods.

Experimental molecular structures in the gas phase at the upper size limit: The case of Si6Tip6.

Currently, the largest (ramax = 19.9Å) and by far the most complicated (234 atoms, C1 symmetry, 696 independent geometrical parameters, and 27 261 interatomic terms) experimental molecular structure of a cage-type Si6Tip6 (Tip = 2,4,6-iPr3C6H2) isomer has been investigated in the gas phase by the electron diffraction method (Tav = 645K) supplemented with theoretical simulations. A detailed analysis of the current possibilities for experimentally investigating large molecular structures is performed. A series of density functional theory approximations and the role of dispersion interactions have been benchmarked using the obtained data. Based on the refined geometry of Si6Tip6, various quantum-chemical methods have been applied for the investigation of the electronic structure of its Si6 core. In particular, natural bond orbital, quantum theory of atoms in molecules, interacting quantum atoms, fractional occupation number weighted density, and complete active space self-consistent field (CASSCF) methods were utilized. The diradical character of the molecule has been assessed by the UHF and CASSCF approximations. The problem of bonding between the hemispheroidal silicon atoms has been investigated.

Extending the Scope of Conceptual Density Functional Theory with Second Order Analytical Methodologies.

In the context of the growing impact of conceptual density functional theory (DFT) as one of the most successful chemical reactivity theories, response functions up to second order have now been widely applied; in recent years, among others, particular attention has been focused on the linear response function and also extensions to higher order have been put forward. As the larger part of these studies have been carried using a finite difference approach to compute these concepts, we now embarked on (an extension of) an analytical approach to conceptual DFT. With the ultimate aim of providing a complete set of analytically computable second order properties, including the softness and hardness kernels, the hardness as the simplest second order response function is scrutinized again with numerical results highlighting the difference in nature between the analytical hardness (referred to as hardness condition) and the Parr-Pearson absolute chemical hardness. The hardness condition is investigated for its capability to gauge the (de)localization error of density functional approximations (DFAs). The analytical Fukui function, besides overcoming the difficulties in the finite difference approach in treating negatively charged systems, also showcases the errors of deviating from the straight-line behavior using fractional occupation number calculations. Subsequently, the softness kernel and its atom-condensed inverse, the hardness matrix, are accessed through the Berkowitz-Parr relation. Revisiting the softness kernel confirms and extends previous discussions on how Kohn's Nearsightedness of Electronic Matter principle can be retrieved and identified as the physicist's version of the chemist's "transferability of functional groups" concept. The accurate, analytical hardness matrix evaluation on the other hand provides further support for the basics of Nalewajski's charge sensitivity analysis. Based on Parr and Liu's functional expansion of the energy functional, a new energy decomposition is introduced with an order of magnitude analysis of the different terms for a series of simple molecules both at their equilibrium geometry and upon variation in bond length and dihedral angle. Finally, for the first time, the perturbation expansion of the energy functional is studied numerically up to second order now that all response functions and integration techniques are at hand. The perturbation expansion energies are in excellent agreement with those obtained directly from DFA calculations giving confidence in the convergence of the perturbation series and its use in judging the importance of the different terms in reactivity investigations.

Evaluation of picture change effects on fractional occupation number states in noble gas atoms

In the two-component relativistic density functional theory, the picture change error (PCE), which originates from insufficient transformation of operators, should be corrected. In this study, we examine the PCE in the fractional occupation number (FON) state based on the spin-free infinite-order two-component Hamiltonian. The PCE for the total and orbital energy changes is estimated with respect to the FON electrons of the highest occupied molecular orbital and 1s core orbital in noble gas atoms. PCE is significant in core orbitals in heavy elements but relatively small in light elements and valence orbitals. The delocalization error, which can be represented by the total energy deviation from the behavior of the exact energy, is overestimated by the lack of picture change correction of the two-electron operator and underestimated by that of the density operator. Corresponding to these results, the PCE influences the value of orbital energies and slope of orbital energy change to FON.

Graph-based quantum response theory and shadow Born-Oppenheimer molecular dynamics.

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born-Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. The proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

One-body reduced density-matrix functional theory for the canonical ensemble

We establish one-body reduced density-matrix functional theory for the canonical ensemble in a finite basis set at an elevated temperature. Including temperature guarantees differentiability of the universal functional by occupying all states and additionally not fully occupying the states in a fermionic system. We use convexity of the universal functional and invertibility of the potential-to-1RDM map to show that the subgradient contains only one element which is equivalent to differentiability. This allows us to show that all 1RDMs with a purely fractional occupation number spectrum ($0 < n_i < 1 \; \forall_i$) are uniquely $v$-representable up to a constant.

Polyradical character assessment using multireference calculations and comparison with density-functional derived fractional occupation number weighted density analysis.

The biradicaloid character of different types of polycyclic aromatic hydrocarbons (PAHs) based on small band gaps is an important descriptor to assess their opto-electronic properties. In this work, the unpaired electron densities and numbers of unpaired electrons (NU values) calculated at the high-level multireference averaged quadratic coupled-cluster (MR-AQCC) method are used to develop a test set to assess the capabilities of different biradical descriptors based on density functional theory. A benchmark collection of 29 different compounds has been selected. The DFT descriptors contain primarily the fractional occupation number weighted electron density (FOD) based on simplified thermally-assisted-occupation density functional theory (TAO-DFT) calculations, but the singlet-triplet energy difference and other descriptors denoted as y0 and nLUNO have been considered as well. After adjustment of the literature-recommended finite temperatures, a very good, detailed agreement between unpaired density and FOD analysis is observed which is also manifested in excellent statistical correlations. The other two descriptors also show good correlations even though the absolute scaling is not satisfactory. A new linear fit of FOD data to the MR-AQCC reference values leads to an improved regression relation for determining the recommended finite temperature value in dependence of the Hartree-Fock exchange. This provides the basis for fast and reliable assessment of the biradical character of many classes of PAHs without the need for performing computationally extended MR calculations.

Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer.

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B2, C2, O2, and F2 in the local spin density approximation (LSDA), the Perdew-Burke-Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C2. This work shows that spin symmetry breaking in singlet C2, which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

DFT-Based Multireference Diagnostics in the Solid State: Application to Metal-Organic Frameworks.

When a many-body wave function of a system cannot be captured by a single determinant, high-level multireference (MR) methods are required to properly explain its electronic structure. MR diagnostics to estimate the magnitude of such static correlation have been primarily developed for molecular systems and range from low in computational cost to as costly as the full MR calculation itself. We report the first application of low-cost MR diagnostics based on the fractional occupation number calculated with finite-temperature DFT to solid-state systems. To compare the behavior of the diagnostics on solids and molecules, we select metal-organic frameworks (MOFs) as model materials because their reticular nature provides an intuitive way to identify molecular derivatives. On a series of closed-shell MOFs, we demonstrate that the DFT-based MR diagnostics are equally applicable to solids as to their molecular derivatives. The magnitude and spatial distribution of the MR character of a MOF are found to have a good correlation with those of its molecular derivatives, which can be calculated much more affordably in comparison to those of the full MOF. The additivity of MR character discussed here suggests the set of molecular derivatives to be a good representation of a MOF for both MR detection and ultimately for MR corrections, facilitating accurate and efficient high-throughput screening of MOFs and other porous solids.

Transition-potential coupled cluster II: optimisation of the core orbital occupation number

The issue of orbital relaxation in computational core-hole spectroscopy, specifically x-ray absorption, has been a major problem for methods such as equation-of-motion coupled cluster with singles and doubles (EOM-CCSD). The transition-potential coupled cluster (TP-CC) method is utilised to address this problem by including an explicit treatment of orbital relaxation via the use of reference orbitals with a fractional core occupation number. The value of the fractional occupation parameter λ was optimised for both TP-CCSD and XTP-CCSD methods in an element-specific manner due to the differences in atomic charge and energy scale. Additionally, TP-CCSD calculations using the optimised parameters were performed for the K-edge absorption spectra of gas-phase adenine and thymine. TP-CCSD reproduces the valence region well and requires smaller overall energy shifts in comparison to EOM-CCSD, while also improving on the relative position and intensities of several absorption peaks.

Core-hole delocalization for modeling x-ray spectroscopies: A cautionary tale.

The influence of core-hole delocalization for x-ray photoelectron, x-ray absorption, and x-ray emission spectrum calculations is investigated in detail using approaches including response theory, transition-potential methods, and ground state schemes. The question of a localized/delocalized vacancy is relevant for systems with symmetrically equivalent atoms, as well as near-degeneracies that can distribute the core orbitals over several atoms. We show that the issues relating to core-hole delocalization are present for calculations considering explicit core-hole states, e.g., when using a core-excited or core-ionized reference state or for fractional occupation numbers. As electron correlation eventually alleviates the issues, but even when using coupled-cluster single-double and perturbative triple, there is a notable discrepancy between core-ionization energies obtained with localized and delocalized core-holes (0.5eV for the carbon K-edge). Within density functional theory, the discrepancy correlates with the exchange interaction involving the core orbitals of the same spin symmetry as the delocalized core-hole. The use of a localized core-hole allows for a reasonably good inclusion of relaxation at a lower level of theory, whereas the proper symmetry solution involving a delocalized core-hole requires higher levels of theory to account for the correlation effects involved in orbital relaxation. For linear response methods, we further show that if x-ray absorption spectra are modeled by considering symmetry-unique sets of atoms, care has to be taken such that there are no delocalizations of the core orbitals, which would otherwise introduce shifts in absolute energies and relative features.

Air Stable Chalcogen-Doped Rubicenes with Diradical Character

Air Stable Chalcogen-Doped Rubicenes with Diradical Character

Reformulation of thermally assisted-occupation density functional theory in the Kohn–Sham framework

We reformulate the thermally assisted-occupation density functional theory (TAO-DFT) into the Kohn-Sham single-determinant framework and construct two new post-self-consistent field (post-SCF) static correlation correction schemes, named rTAO and rTAO-1. In contrast to the original TAO-DFT with the density in an ensemble form, in which each orbital density is weighted with a fractional occupation number, the ground-state density is given by a single-determinant wavefunction, a regular Kohn-Sham (KS) density, and total ground state energy is expressed in the normal KS form with a static correlation energy formulated in terms of the KS orbitals. In post-SCF calculations with rTAO functionals, an efficient energy scanning to quantitatively determine θ is also proposed. The rTAOs provide a promising method to simulate systems with strong static correlation as original TAO, but simpler and more efficient. We show that both rTAO and rTAO-1 is capable of reproducing most results from TAO-DFT without the additional functional Eθ used in TAO-DFT. Furthermore, our numerical results support that, without the functional Eθ, both rTAO and rTAO-1 can capture correct static correlation profiles in various systems.

Assessing Alkene Reactivity toward Cytochrome P450-Mediated Epoxidation through Localized Descriptors and Regression Modeling.

The prediction of sites of epoxidation by cytochrome P450s during metabolism is particularly important in drug design, as epoxides are capable of alkylating biological macromolecules. Reliable methods are needed to quantitatively predict P450-mediated epoxidation barriers for inclusion in high-throughput screening campaigns alongside protein-ligand docking. Utilizing the fractional occupation number weighted density (FOD) and orbital-weighted Fukui index (fw+) as descriptors of local reactivity and a data set of 36 alkene epoxidation barriers computed with density functional theory (DFT), we developed and validated a multiple linear regression model for the reliable estimation of epoxidation barriers using only substrate structures as input. Using our recommended level of theory (GFN2-xTB//GFN-FF), mean absolute errors in the training and test sets were found to be 0.66 and 0.70 kcal/mol, respectively, with coefficients of determination of ca. 0.80. We demonstrate the utility of this approach on three known substrates of CYP101A1 and further show that this approach is inappropriate for particularly electron-rich alkenes. By employing a modern semiempirical method on force-field-generated geometries, the required descriptors can be calculated on the millisecond timescale per structure, making the approach well suited for incorporation into high-throughput methodologies alongside docking.

Exact Second-Order Corrections and Accurate Quasiparticle Energy Calculations in Density Functional Theory.

We develop a second-order correction to commonly used density functional approximations (DFAs) to eliminate the systematic delocalization error. The method, based on the previously developed global scaling correction (GSC), is an exact quadratic correction to the DFA for the fractional charge behavior and uses the analytical second derivatives of the total energy with respect to fractional occupation numbers of the canonical molecular orbitals. For small and medium-size molecules, this correction leads to ground-state orbital energies that are a highly accurate approximation to the corresponding quasiparticle energies. It provides excellent predictions of ionization potentials, electron affinities, photoemission spectrum, and photoexcitation energies beyond previous approximate second-order approaches, thus showing potential for broad applications in computational spectroscopy.

Mixed Precision Fermi-Operator Expansion on Tensor Cores from a Machine Learning Perspective.

We present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia's A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent the fractional occupation numbers of the electronic states at finite temperatures.

Variational Density Functional Calculations of Excited States via Direct Optimization

The development of variational density functional theory approaches to excited electronic states is impeded by limitations of the commonly used self-consistent field (SCF) procedure. A method based on a direct optimization approach as well as the maximum overlap method is presented, and the performance is compared with previously proposed SCF strategies. Excited-state solutions correspond to saddle points of the energy as a function of the electronic degrees of freedom. The approach presented here makes use of a preconditioner determined with the help of the maximum overlap method to guide the convergence on a target nth-order saddle point. This method is found to be more robust and converge faster than previously proposed SCF approaches for a set of 89 excited states of molecules. A limited-memory formulation of the symmetric rank-one method for updating the inverse Hessian is found to give the best performance. A conical intersection for the carbon monoxide molecule is calculated without resorting to fractional occupation numbers. Calculations on the excited states of the hydrogen atom and a doubly excited state of the dihydrogen molecule using a self-interaction corrected functional are presented. For these systems, the self-interaction correction is found to improve the accuracy of density functional calculations of excited states.

Rapid Detection of Strong Correlation with Machine Learning for Transition-Metal Complex High-Throughput Screening

Despite its widespread use in chemical discovery, approximate density functional theory (DFT) is poorly suited to many targets, such as those containing open-shell, 3d transition metals that can be expected to have strong multireference (MR) character. For discovery workflows to be predictive, we need automated, low-cost methods that can distinguish the regions of chemical space where DFT should be applied from those where it should not. We curate more than 4800 open-shell transition-metal complexes up to hundreds of atoms in size from prior high-throughput DFT studies and evaluate affordable, finite-temperature DFT fractional occupation number (FON)-based MR diagnostics. We show that intuitive measures of strong correlation (i.e., the HOMO-LUMO gap) are not predictive of MR character as judged by FON-based diagnostics. Analysis of independently trained machine learning (ML) models to predict HOMO-LUMO gaps and FON-based diagnostics reveals differences in the metal and ligand sensitivity of the two quantities. We use our trained ML models to rapidly evaluate MR character over a space of ∼187000 theoretical complexes, identifying large-scale trends in spin-state-dependent MR character and finding small HOMO-LUMO gap complexes while ensuring low MR character.

DFT, DLPNO-CCSD(T), and NEVPT2 benchmark study of the reaction between ferrocenium and trimethylphosphine.

The reaction between ferrocenium and trimethylphosphine was studied using density functional theory (DFT), domain-based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)), and N-electron valence state perturbation theory (NEVPT2). The accuracy of the DFT functionals decreases compared to the DLPNO-CCSD(T) level in the following order: M06-L > TPSS > M06, BLYP > PBE, PBE0, B3LYP > > PWPB95 > > DSD-BLYP. The roles of thermochemical, continuum solvation (SMD), and counterpoise corrections were evaluated. Grimme's D3 empirical dispersion correction is essential for all functionals studied except M06 and M06-L. The reliability of the frequency calculations performed directly within the SMD was confirmed. The systems showed no significant multireference character according to T1 and T2 diagnostics and the fractional occupation number (FOD) weighted electron density analysis. The multireference NEVPT2 calculations gave qualitatively valid conclusions about the reaction mechanism. However, a multireference approach is generally not recommended because it requires arbitrary chosen active spaces.