Search

Artificial intelligence-enhanced quantum mechanical method 1 (AIQM1) is a general-purpose method that was shown to achieve high accuracy for many applications with a speed close to its baseline semiempirical quantum mechanical (SQM) method ODM2*. Here, we evaluate the hitherto unknown performance of out-of-the-box AIQM1 without any refitting for reaction barrier heights on eight datasets, including a total of ∼24 thousand reactions. This evaluation shows that AIQM1's accuracy strongly depends on the type of transition state and ranges from excellent for rotation barriers to poor for, e.g., pericyclic reactions. AIQM1 clearly outperforms its baseline ODM2* method and, even more so, a popular universal potential, ANI-1ccx. Overall, however, AIQM1 accuracy largely remains similar to SQM methods (and B3LYP/6-31G* for most reaction types) suggesting that it is desirable to focus on improving AIQM1 performance for barrier heights in the future. We also show that the built-in uncertainty quantification helps in identifying confident predictions. The accuracy of confident AIQM1 predictions is approaching the level of popular density functional theory methods for most reaction types. Encouragingly, AIQM1 is rather robust for transition state optimizations, even for the type of reactions it struggles with the most. Single-point calculations with high-level methods on AIQM1-optimized geometries can be used to significantly improve barrier heights, which cannot be said for its baseline ODM2* method. [ABSTRACT FROM AUTHOR]

A computational protocol is developed for efficient studies of partially reduced redox-active oxides using the self-consistent charge density functional tight-binding method. The protocol is demonstrated for ceria, which is a prototypical reducible oxide material. The underlying idea is to achieve a consistent (and harmonized) set of Slater–Koster (SK) tables with connected repulsive potentials that enable switching on and off the in-valence description of the Ce 4f states without serious loss of accuracy in structure and energetics. The implicit treatment of the Ce 4f states, with the use of f-in-core SK-tables, is found to lead to a significant decrease in computational time. More importantly, it allows for explicit control of the oxidation states of individual Ce atoms. This makes it possible to "freeze" the electronic configuration, thereby allowing the exploration of the energetics for various meta-stable configurations. We anticipate that the outlined strategy can help to shed light on the interplay between the size, shape, and redox activity for nanoceria and other related materials. [ABSTRACT FROM AUTHOR]

Here, extensions to quantum chemical nanoreactor molecular dynamics simulations for discovering complex reactive events are presented. The species-selective algorithm, where the nanoreactor effectively works for the selected desired reactants, was introduced to the original scheme. Moreover, for efficient simulations of large model systems with the modified approach, the divide-and-conquer linear-scaling density functional tight-binding method was exploited. Two illustrative applications of the polymerization of propylene and cyclopropane mixtures and the aggregation of sodium chloride from aqueous solutions indicate that species-selective quantum chemical nanoreactor molecular dynamics is a promising method to accelerate the sampling of multicomponent chemical processes proceeding under relatively mild conditions. [ABSTRACT FROM AUTHOR]

Herein, we present a new method to efficiently calculate electronically excited states in large molecular assemblies, consisting of hundreds of molecules. For this purpose, we combine the long-range corrected tight-binding density functional fragment molecular orbital method (FMO-LC-DFTB) with an excitonic Hamiltonian, which is constructed in the basis of locally excited and charge-transfer configuration state functions calculated for embedded monomers and dimers and accounts explicitly for the electronic coupling between all types of excitons. We first evaluate both the accuracy and efficiency of our fragmentation approach for molecular dimers and aggregates by comparing it with the full LC-TD-DFTB method. The comparison of the calculated spectra of an anthracene cluster shows a very good agreement between our method and the LC-TD-DFTB reference. The effective computational scaling of our method has been explored for anthracene clusters and for perylene bisimide aggregates. We demonstrate the applicability of our method by the calculation of the excited state properties of pentacene crystal models consisting of up to 319 molecules. Furthermore, the participation ratio of the monomer fragments to the excited states is analyzed by the calculation of natural transition orbital participation numbers, which are verified by the hole and particle density for a chosen pentacene cluster. The use of our FMO-LC-TDDFTB method will allow for future studies of excitonic dynamics and charge transport to be performed on complex molecular systems consisting of thousands of atoms. [ABSTRACT FROM AUTHOR]

Kernel-based methods, including Gaussian process regression (GPR) and generally kernel ridge regression, have been finding increasing use in computational chemistry, including the fitting of potential energy surfaces and density functionals in high-dimensional feature spaces. Kernels of the Matern family, such as Gaussian-like kernels (basis functions), are often used which allow imparting to them the meaning of covariance functions and formulating GPR as an estimator of the mean of a Gaussian distribution. The notion of locality of the kernel is critical for this interpretation. It is also critical to the formulation of multi-zeta type basis functions widely used in computational chemistry. We show, on the example of fitting of molecular potential energy surfaces of increasing dimensionality, the practical disappearance of the property of locality of a Gaussian-like kernel in high dimensionality. We also formulate a multi-zeta approach to the kernel and show that it significantly improves the quality of regression in low dimensionality but loses any advantage in high dimensionality, which is attributed to the loss of the property of locality. [ABSTRACT FROM AUTHOR]

We present analytic gradients for local density fitting Hartree–Fock (HF) and hybrid Kohn–Sham (KS) density functional methods. Due to the non-variational nature of the local fitting algorithm, the method of Lagrange multipliers is used to avoid the solution of the coupled perturbed HF and KS equations. We propose efficient algorithms for the solution of the arising Z-vector equations and the gradient calculation that preserve the third-order scaling and low memory requirement of the original local fitting algorithm. In order to demonstrate the speed and accuracy of our implementation, gradient calculations and geometry optimizations are presented for various molecular systems. Our results show that significant speedups can be achieved compared to conventional density fitting calculations without sacrificing accuracy. [ABSTRACT FROM AUTHOR]

Surface platinum hydride structures may exist and play a potentially important role during electrocatalysis and cathodic corrosion of Pt(111). Earlier work on platinum hydrides suggests that Pt may form clusters with multiple equivalents of hydrogen. Here, using thermodynamic methods and density functional theory, we compared several surface hydride structures on Pt(111). The structures contain multiple monolayers of hydrogen in or near the surface Pt layer. The hydrogen in these structures may bind the subsurface or reconstruct the surface both in the set of initial configurations and in the resulting (meta)stable structures. Multilayer stable configurations share one monolayer of subsurface H stacking between the top two Pt layers. The structure containing two monolayers (MLs) of H is formed at −0.29 V vs normal hydrogen electrode, is locally stable with respect to configurations with similar H densities, and binds H neutrally. Structures with 3 and 4 ML H form at −0.36 and −0.44 V, respectively, which correspond reasonably well to the experimental onset potential of cathodic corrosion on Pt(111). For the 3 ML configuration, the top Pt layer is reconstructed by interstitial H atoms to form a well-ordered structure with Pt atoms surrounded by four, five, or six H atoms in roughly square-planar and octahedral coordination patterns. Our work provides insight into the operando surface state during low-potential reduction reactions on Pt(111) and shows a plausible precursor for cathodic corrosion. [ABSTRACT FROM AUTHOR]

A systematic study is made of the accuracy and efficiency of a number of existing quadrature schemes for molecular Kohn–Sham Density-Functional Theory (DFT) using 408 molecules and 254 chemical reactions. Included are the fixed SG-x (x = 0–3) grids of Gill et al., Dasgupta, and Herbert, the 3-zone grids of Treutler and Ahlrichs, a fixed five-zone grid implemented in Molpro, and a new adaptive grid scheme. While all methods provide a systematic reduction of errors upon extension of the grid sizes, significant differences are observed in the accuracies for similar grid sizes with various approaches. For the tests in this work, the SG-x fixed grids are less suitable to achieve high accuracies in the DFT integration, while our new adaptive grid performed best among the schemes studied in this work. The extra computational time to generate the adaptive grid scales linearly with molecular size and is negligible compared with the time needed for the self-consistent field iterations for large molecules. A comparison of the grid accuracies using various density functionals shows that meta-GGA functionals need larger integration grids than GGA functionals to reach the same degree of accuracy, confirming previous investigations of the numerical stability of meta-GGA functionals. On the other hand, the grid integration errors are almost independent of the basis set, and the basis set errors are mostly much larger than the errors caused by the numerical integrations, even when using the smallest grids tested in this work. [ABSTRACT FROM AUTHOR]

We present a quasiparticle self-consistent GW (QSGW) implementation for periodic systems based on crystalline Gaussian basis sets. Our QSGW approach is based on a full-frequency analytic continuation GW scheme with Brillouin zone sampling and employs the Gaussian density fitting technique. We benchmark our QSGW implementation on a set of weakly correlated semiconductors and insulators as well as strongly correlated transition metal oxides, including MnO, FeO, CoO, and NiO. The band gap, band structure, and density of states are evaluated using finite size corrected QSGW. We find that although QSGW systematically overestimates the bandgaps of the tested semiconductors and transition metal oxides, it completely removes the dependence on the choice of density functionals and provides a more consistent prediction of spectral properties than G0W0 across a wide range of solids. This work paves the way for utilizing QSGW in ab initio quantum embedding for solids. [ABSTRACT FROM AUTHOR]

Revealing the effect of external applied potential on the reaction mechanism and product selectivity is of great significance in electrochemical studies. In this work, the grand canonical density functional theory method was applied to simulate the explicit electrocatalytic process of oxygen evolution reaction and electrochemical ozone production due to the O3 product sensitivity toward the applied potential. Over the Pt/Pd single atom embedded on B/N co-doped graphene (Pt/Pd-BNC) surface, crossover points of O2/O3 selectivity inversion were predicted to be 1.33 and 0.89 V vs standard hydrogen electrode, which were also consistent with the previous experimental results. An in-depth analysis of the energetic terms in the reaction free energies also found the considerable impact of the applied potential on the Helmholtz free energy term, with optimal potential predicted for the key elementary steps, and linear correlations between electrode potential (U) and reaction free energy were found for each elementary step. This study offers extensive knowledge on the potential effect on the O2/O3 selective formation on two-dimensional anode surfaces and provides new insights for investigating the reactivity/selectivity on electrode surfaces in real reaction conditions. [ABSTRACT FROM AUTHOR]

Polaron formation following optical absorption is a key process that defines the photophysical properties of many semiconducting transition metal oxides, which comprise an important class of materials with potential optoelectronic and photocatalytic applications. In this work, we use hematite (α-Fe2O3) as a model transition metal oxide semiconductor to demonstrate the feasibility of direct optical population of band edge polaronic states. We employ first-principles electron–phonon computations within the framework of the density functional theory+U+J method to reveal the presence of these states within a thermal distribution of phonon displacements and model their evolution with temperature. Our computations reproduce the temperature dependence of the optical dielectric function of hematite with remarkable accuracy and indicate that the band edge optical absorption and second-order resonance Raman spectra arise from polaronic optical transitions involving coupling to longitudinal optical phonons with energies greater than 50 meV. Additionally, we find that the resulting polaron comprises an electron localized to two adjacent Fe atoms with distortions that lie primarily along the coordinates of phonons with energies of 31 and 81 meV. [ABSTRACT FROM AUTHOR]

Most computational studies in chemistry and materials science are based on the use of density functional theory. Although the exact density functional is unknown, several density functional approximations (DFAs) offer a good balance of affordable computational cost and semi-quantitative accuracy for applications. The development of DFAs still continues on many fronts, and several new DFAs aiming for improved accuracy are published every year. However, the numerical behavior of these DFAs is an often-overlooked problem. In this work, we look at all 592 DFAs for three-dimensional systems available in Libxc 5.2.2 and examine the convergence of the density functional total energy based on tabulated atomic Hartree–Fock wave functions. We show that several recent DFAs, including the celebrated SCAN family of functionals, show impractically slow convergence with typically used numerical quadrature schemes, making these functionals unsuitable both for routine applications and high-precision studies, as thousands of radial quadrature points may be required to achieve sub-μEh accurate total energies for these functionals, while standard quadrature grids like the SG-3 grid only contain O (100) radial quadrature points. These results are both a warning to users to always check the sufficiency of the quadrature grid when adopting novel functionals, as well as a guideline to the theory community to develop better-behaved density functionals. [ABSTRACT FROM AUTHOR]

We study the performance of spin-component-scaled second-order Møller–Plesset perturbation theory (SCS-MP2) for the prediction of the lattice constant, bulk modulus, and cohesive energy of 12 simple, three-dimensional covalent and ionic semiconductors and insulators. We find that SCS-MP2 and the simpler scaled opposite-spin MP2 (SOS-MP2) yield predictions that are significantly improved over the already good performance of MP2. Specifically, when compared to experimental values with zero-point vibrational corrections, SCS-MP2 (SOS-MP2) yields mean absolute errors of 0.015 (0.017) Å for the lattice constant, 3.8 (3.7) GPa for the bulk modulus, and 0.06 (0.08) eV for the cohesive energy, which are smaller than those of leading density functionals by about a factor of two or more. We consider a reparameterization of the spin-scaling parameters and find that the optimal parameters for these solids are very similar to those already in common use in molecular quantum chemistry, suggesting good transferability and reliable future applications to surface chemistry on insulators. [ABSTRACT FROM AUTHOR]

Sun et al. [J. Chem. Phys. 144, 191101 (2016)] suggested that common density-functional approximations (DFAs) should exhibit large energy errors for excited states as a necessary consequence of orbital nodality. Motivated by self-interaction corrected density-functional calculations on many-electron systems, we continue their study with the exactly solvable 1s, 2p, and 3d states of 36 hydrogenic one-electron ions (H–Kr35+) and demonstrate with self-consistent calculations that state-of-the-art DFAs indeed exhibit large errors for the 2p and 3d excited states. We consider 56 functionals at the local density approximation (LDA), generalized gradient approximation (GGA) as well as meta-GGA levels, and several hybrid functionals such as the recently proposed machine-learned DM21 local hybrid functional. The best non-hybrid functional for the 1s ground state is revTPSS. As predicted by Sun et al., the 2p and 3d excited states are more difficult for DFAs, and LDA functionals turn out to yield the most systematic accuracy for these states among non-hybrid functionals. The best performance for the three states overall is observed with the BHandH global hybrid GGA functional, which contains 50% Hartree–Fock exchange and 50% LDA exchange. The performance of DM21 is found to be inconsistent, yielding good accuracy for some states and systems and poor accuracy for others. Based on these results, we recommend including a variety of one-electron cations in future training of machine-learned density functionals. [ABSTRACT FROM AUTHOR]

A new vacuum ultraviolet absorption (VUV) spectrum of azulene vapor has been obtained by using a synchrotron radiation source. The onset of the ultraviolet spectrum, previously reported by Sidman et al., has been analyzed in detail by Franck–Condon (FC) and Herzberg–Teller (HT) methods. The photoelectron spectral profile identifies the 3px-Rydberg state 00 band to be 131 cm−1 from the VUV maximum. Excited state energy levels were calculated by three independent methods: the wide scan VUV spectrum was correlated with symmetry adapted cluster configuration interaction calculations. The low energy portion of the spectrum was studied by both time dependent density functional theoretical methods (TDDFT) and multi-reference multi-root CI (MRD-CI). Equilibrium structures were determined for valence states at the TDDFT level. Rydberg states were determined by both TDDFT and MRD-CI. The FC + HT analyses were performed on the TDDFT wave-functions. The HT intensity profiles are generally low in intensity, relative to the FC ones; however, HT is dominant in the second singlet state (S2, 11A1). As a result, numerous non-symmetric modes, their overtones, and combination bands show considerable intensity in that band. Energies obtained from use of extremely diffuse s-, p-, d-, or f-character functions enabled realistic extrapolation to the IE1 for many Rydberg states (RS). The lowest RS (3b13s) based on IE2 lies at 4.804 eV with a quantum defect of 0.714. Differentiation between valence and RS is readily made using the second moments of the charge distribution. [ABSTRACT FROM AUTHOR]

Glutaredoxins are small enzymes that catalyze the oxidation and reduction of protein disulfide bonds by the thiol–disulfide exchange mechanism. They have either one or two cysteines in their active site, resulting in different catalytic reaction cycles that have been investigated in many experimental studies. However, the exact mechanisms are not yet fully known, and to our knowledge, no theoretical studies have been performed to elucidate the underlying mechanism. In this study, we investigated a proposed mechanism for the reduction of the disulfide bond in the protein HMA4n by a mutated monothiol Homo sapiens glutaredoxin and the co-substrate glutathione. The catalytic cycle involves three successive thiol–disulfide exchanges that occur between the molecules. To estimate the regioselectivity of the different attacks, classical molecular dynamics simulations were performed and the trajectories analyzed regarding the sulfur–sulfur distances and the attack angles between the sulfurs. The free energy profile of each reaction was obtained with hybrid quantum mechanical/molecular mechanical metadynamics simulations. Since this required extensive phase space sampling, the semi-empirical density functional tight-binding method was used to describe the reactive cysteines. For an accurate description, we used specific reaction parameters fitted to B3LYP energies of the thiol–disulfide exchange and a machine learned energy correction that was trained on coupled-cluster single double perturbative triple [CCSD(T)] energies of thiol–disulfide exchanges. Our calculations show the same regiospecificity as observed in the experiment, and the obtained barrier heights are about 12 and 20 kcal/mol for the different reaction steps, which confirms the proposed pathway. [ABSTRACT FROM AUTHOR]

The core excitation energies and related principal ionization energies are obtained for selected molecules using several density functionals and compared with benchmark equation-of-motion coupled cluster (EOM-CC) results. Both time-dependent and time-independent formulations of excitation spectra in the time-dependent density functional theory and the EOM-CC are employed to obtain excited states that are not always easily accessible with the time-independent method. Among those functionals, we find that the QTP(00) functional, which is only parameterized to reproduce the five IPs of water, provides excellent core IPs and core excitation energies, consistently yielding better excitation and ionization energies. We show that orbital eigenvalues of KS density functional theory play an important role in determining the accuracy of the excitation and photoelectron spectra. [ABSTRACT FROM AUTHOR]

Nonadiabatic excited state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, and charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamics (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. The nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Furthermore, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well to those obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials. [ABSTRACT FROM AUTHOR]

The semilocal form of an exchange hole is highly useful in developing non-local range-separated hybrid density functionals for finite and extended systems. The way to construct the conventional exact exchange hole model is based on either the Taylor series expansion or the reverse engineering technique from the corresponding exchange energy functional. Although the latter technique is quite popular in the context of generalized gradient approximation (GGA) functionals, the same for the meta-GGA functionals is not so much explored. Thus, in this study, we propose a reverse-engineered semilocal exchange hole of a meta-GGA functional, which only depends on the meta-GGA ingredient α (also known as the Pauli kinetic energy enhancement factor). The model is subsequently used to design the short-range-separated meta-GGA hybrid density functional. We show that the present method can be successfully applied for several challenging problems in the context of solids, especially for which the GGA based hybrid fails drastically. This assessment proves that the present functional is quite useful for materials sciences. Finally, we also use this method for several molecular test cases, where the results are also as comparative as its base semilocal functional. [ABSTRACT FROM AUTHOR]

The clustering, collision, and relaxation dynamics of pristine and doped helium nanodroplets is theoretically investigated in cases of pickup and clustering of heliophilic argon, collision of heliophobic cesium atoms, and coalescence of two droplets brought into contact by their mutual long-range van der Waals interaction. Three approaches are used and compared with each other. The He time-dependent density functional theory method considers the droplet as a continuous medium and accounts for its superfluid character. The ring-polymer molecular dynamics method uses a path-integral description of nuclear motion and incorporates zero-point delocalization while bosonic exchange effects are ignored. Finally, the zero-point averaged dynamics approach is a mixed quantum–classical method in which quantum delocalization is described by attaching a frozen wavefunction to each He atom, equivalent to classical dynamics with effective interaction potentials. All three methods predict that the growth of argon clusters is significantly hindered by the helium host droplet due to the impeding shell structure around the dopants and kinematic effects freezing the growing cluster in metastable configurations. The effects of superfluidity are qualitatively manifested by different collision dynamics of the heliophilic atom at high velocities, as well as quadrupole oscillations that are not seen with particle-based methods, for droplets experiencing a collision with cesium atoms or merging with each other. [ABSTRACT FROM AUTHOR]

In this Communication, we assess a panel of 18 double-hybrid density functionals for the modeling of the thermochemical and kinetic properties of an extended dataset of 449 organic chemistry reactions belonging to the BH9 database. We show that most of DHs provide a statistically robust performance to model barrier height and reaction energies in reaching the "chemical accuracy." In particular, we show that nonempirical DHs, such as PBE0-DH and PBE-QIDH, or minimally parameterized alternatives, such as ωB2PLYP and B2K-PLYP, succeed to accurately model both properties in a balanced fashion. We demonstrate, however, that parameterized approaches, such as ωB97X-2 or DSD-like DHs, are more biased to only one of both properties. [ABSTRACT FROM AUTHOR]

Diamine-functionalized M2(dobpdc) (M = Mg, Mn, Fe, Co, Zn) metal–organic frameworks (MOFs) are among a growing class of crystalline solids currently being intensively investigated for carbon capture as they exhibit a novel cooperative and selective CO2 adsorption mechanism and a step-shaped isotherm. To understand their CO2 adsorption behavior, ab initio calculations with near-chemical accuracy (∼6 kJ/mol, an average experimental error) are required. Here, we present density functional theory (DFT) calculations of CO2 adsorption in m-2-m–Zn2(dobpdc) (m-2-m = N,N′-dimethylethyle-nediamine and dobpdc4− = 4,4′-dioxidobiphenyl-3,3′-dicarboxylate) with different exchange–correlation functionals, including semilocal functionals [Perdew–Burke–Ernzerhof (PBE) and two revised PBE functionals], semiempirical pairwise corrections (D3 and Tkatchenko–Scheffler), nonlocal van der Waals (vdW) correlation functionals—vdW-optB88 (or vdW-DF-optB88), vdW-DF1, vdW-DF2, vdW-DF2-B86R (or rev-vdW-DF2), vdW-DF-cx (and vdW-DF-cx0), and revised VV10—and the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (GGA). Overall, we find that revPBE+D3 and RPBE+D3 show the best balance of performance for both the lattice parameters and the CO2 binding enthalpy of m-2-m–Zn2(dobpdc). revPBE+D3 and RPBE+D3 predict the m-2-m–Zn2(dobpdc) lattice parameters to within 1.4% of experiment and predict CO2 binding enthalpies of −68 kJ/mol, which compare reasonably well with the experiment (−57 kJ/mol). Although PBE (−57.7 kJ/mol), vdW-DF1 (−49.6 kJ/mol), and vdW-DF2 (−44.3 kJ/mol) are also found to predict the CO2 binding enthalpy with good accuracy, they overestimate lattice parameters and bond lengths. The other functionals considered predict the lattice parameters with the same accuracy as revPBE+D3 and RPBE+D3, but they overbind CO2 by around 26–50 kJ/mol. We find that the superior performance of revPBE+D3 and RPBE+D3 is sustained for the formation enthalpy and the lattice parameters of ammonium carbamate, a primary product of the cooperative CO2 insertion in diamine-functionalized M2(dobpdc) MOFs. Moreover, we find that their performance is derived from their larger repulsive exchange contributions to the CO2 binding enthalpy than the other functionals at the relevant range of the reduced density gradient value for the energetics of CO2 adsorption in the m-2-m–Zn2(dobpdc) MOF. A broader examination of the performance of RPBE+D3 for the structural parameters and CO2 binding enthalpies of 13 diamine-functionalized Mg2(dobpdc) MOFs further demonstrates that RPBE+D3 successfully reproduces experimental CO2 binding enthalpies and reveals a logarithmic relationship between the step pressure and the CO2 binding enthalpy of the diamine-functionalized Mg2(dobpdc) MOFs, consistent with experiments where available. The results of our benchmarking study can help guide the further development of versatile vdW-corrected DFT methods with predictive accuracy. [ABSTRACT FROM AUTHOR]

An accurate description of the electron correlation energy in closed- and open-shell molecules is shown to be obtained by a second-order perturbation theory (PT) termed REMP. REMP is a hybrid of the Retaining the Excitation degree (RE) and the Møller–Plesset (MP) PTs. It performs particularly encouragingly in an orbital-optimized variant (OO-REMP) where the reference wavefunction is given by an unrestricted Slater determinant whose spin orbitals are varied such that the total energy becomes a minimum. While the approach generally behaves less satisfactorily with unrestricted Hartree–Fock references, reasonable performance is observed for restricted Hartree–Fock and restricted open-shell Hartree–Fock references. Inclusion of single excitations to OO-REMP is investigated and found—as in similar investigations—to be dissatisfying as it deteriorates performance. For the non-multireference subset of the accurate W4-11 benchmark set of Karton et al. [Chem. Phys. Lett. 510, 165–178 (2011)], OO-REMP predicts most atomization and reaction energies with chemical accuracy (1 kcal mol−1) if complete-basis-set extrapolation with augmented and core-polarized basis sets is used. For the W4-11 related test-sets, the error estimates obtained with the OO-REMP method approach those of coupled-cluster with singles, doubles and perturbative triples [CCSD(T)] within 20%–35%. The best performance of OO-REMP is found for a mixing ratio of 20%:80% MP:RE, which is essentially independent of whether radical stabilization energies, barrier heights, or reaction energies are investigated. Orbital optimization is shown to improve the REMP approach for both closed and open shell cases and outperforms coupled-cluster theory with singles and doubles (CCSD), spin-component scaled Møller-Plesset theory at second order (SCS-MP2), and density functionals, including double hybrids in all the cases considered. [ABSTRACT FROM AUTHOR]

Electrons in zero external magnetic field can be studied with the Kohn–Sham (KS) scheme of either density functional theory (DFT) or spin-DFT (SDFT). The latter is normally used for open-shell systems because its approximations appear to model better the exchange and correlation (xc) functional, but also because, so far the application of DFT implied a closed-shell-like approximation. In the first part of this Communication, we show that correcting this error for open shells allows the approximate DFT xc functionals to become as accurate as those in SDFT. In the second part, we consider the behavior of SDFT for zero magnetic field. We show that the KS equations of SDFT emerge as the generalized KS equations of DFT in this limit, thus establishing a so far unknown link between the two theories. [ABSTRACT FROM AUTHOR]

Conformational polymorphs of organic molecular crystals represent a challenging test for quantum chemistry because they require careful balancing of the intra- and intermolecular interactions. This study examines 54 molecular conformations from 20 sets of conformational polymorphs, along with the relative lattice energies and 173 dimer interactions taken from six of the polymorph sets. These systems are studied with a variety of van der Waals-inclusive density functionals theory models; dispersion-corrected spin-component-scaled second-order Møller–Plesset perturbation theory (SCS-MP2D); and domain local pair natural orbital coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)]. We investigate how delocalization error in conventional density functionals impacts monomer conformational energies, systematic errors in the intermolecular interactions, and the nature of error cancellation that occurs in the overall crystal. The density functionals B86bPBE-XDM, PBE-D4, PBE-MBD, PBE0-D4, and PBE0-MBD are found to exhibit sizable one-body and two-body errors vs DLPNO-CCSD(T) benchmarks, and the level of success in predicting the relative polymorph energies relies heavily on error cancellation between different types of intermolecular interactions or between intra- and intermolecular interactions. The SCS-MP2D and, to a lesser extent, ωB97M-V models exhibit smaller errors and rely less on error cancellation. Implications for crystal structure prediction of flexible compounds are discussed. Finally, the one-body and two-body DLPNO-CCSD(T) energies taken from these conformational polymorphs establish the CP1b and CP2b benchmark datasets that could be useful for testing quantum chemistry models in challenging real-world systems with complex interplay between intra- and intermolecular interactions, a number of which are significantly impacted by delocalization error. [ABSTRACT FROM AUTHOR]

Range-separated hybrid approximations to the exchange–correlation density functional mix exact and semi-local exchange in a position-dependent manner. In their conventional form, the range separation is controlled by a constant parameter. Turning this constant into a density functional leads to a locally space-dependent range-separation function and thus a more powerful and flexible range-separation approach. In this work, we explore the self-consistent implementation of a local range-separated hybrid, taking into account a one-electron self-interaction correction and the behavior under uniform density scaling. We discuss different forms of the local range-separation function that depend on the electron density, its gradient, and the kinetic energy density. For test sets of atomization energies, reaction barrier heights, and total energies of atoms, we demonstrate that our best model is a clear improvement over common global range-separated hybrid functionals and can compete with density functionals that contain multiple empirical parameters. Promising results for equilibrium bond lengths, harmonic vibrational frequencies, and vertical ionization potentials further underline the potential and flexibility of our approach. [ABSTRACT FROM AUTHOR]

The structure of the double-layer formed at the surface of carbon electrodes is governed by the interactions between the electrode and the electrolyte species. However, carbon is notoriously difficult to simulate accurately, even with well-established methods such as electronic density functional theory and molecular dynamics. Here, we focus on the important case of a lithium ion in contact with the surface of graphite, and we perform a series of reference quantum Monte Carlo calculations that allow us to benchmark various electronic density functional theory functionals. We then fit an accurate carbon–lithium pair potential, which is used in molecular density functional theory calculations to determine the free energy of the adsorption of the ion on the surface in the presence of water. The adsorption profile in aqueous solution differs markedly from the gas phase results, which emphasize the role of the solvent on the properties of the double-layer. [ABSTRACT FROM AUTHOR]

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. We present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods. [ABSTRACT FROM AUTHOR]

Due to their improved accuracy, double-hybrid density functionals emerged as an important method for molecular electronic-structure calculations. The high computational costs of double-hybrid calculations in the condensed phase and the lack of efficient gradient implementations thereof inhibit a wide applicability for periodic systems. We present an implementation of forces and stress tensors for double-hybrid density functionals within the Gaussian and plane-waves electronic structure framework. The auxiliary density matrix method is used to reduce the overhead of the Hartree–Fock kernel providing an efficient and accurate methodology to tackle condensed phase systems. First applications to water systems of different densities and molecular crystals show the efficiency of the implementation and pave the way for advanced studies. Finally, we present large benchmark systems to discuss the performance of our implementation on modern large-scale computers. [ABSTRACT FROM AUTHOR]

We re-examine the recently introduced basis-set correction theory based on density-functional theory, which consists of correcting the basis-set incompleteness error of wave-function methods using a density functional. We use a one-dimensional model Hamiltonian with delta-potential interactions, which has the advantage of making easier to perform a more systematic analysis than for three-dimensional Coulombic systems while keeping the essence of the slow basis convergence problem of wave-function methods. We provide some mathematical details about the theory and propose a new variant of basis-set correction, which has the advantage of being suited to the development of an adapted local-density approximation. We show, indeed, how to develop a local-density approximation for the basis-set correction functional, which is automatically adapted to the basis set employed, without resorting to range-separated density-functional theory as in previous studies, but using instead a finite uniform electron gas whose electron–electron interaction is projected on the basis set. The work puts the basis-set correction theory on firmer ground and provides an interesting strategy for the improvement of this approach. [ABSTRACT FROM AUTHOR]

Benzvalyne (C6H4) is a bicyclic structural isomer of o-benzyne that some typically trusted levels of theory do not report as a minimum on the potential energy surface (PES). The structure was found to be a C2v minimum at the MCSCF, MP2, coupled-cluster single double, coupled-cluster single double triple (CCSDT)-1b, and CCSDT-2 levels of theory. Density functionals at the B3LYP-D3, B2PLYP-D3, and M06-D3 levels also produced a minimum structure. On the other hand, the CCSD(T), CCSD(T)-F12, and CCSDT-1a methods produced a single imaginary frequency for benzvalyne. However, the increase in the correlation for the CCSDT-1b and CCSDT-2 methods implies that benzvalyne is, in fact, a true, if highly strained, minimum on the PES. The C–C≡C bond angle was found to be only 108°; this angle is 180° for an unstrained C–C≡C triple bond moiety. As a result, the strain energy is notably high at 159 kcal mol−1. Comparing the strain energy of the rest of the molecule gives a strain energy of 92 kcal mol−1 for this triple bond region alone. The computed harmonic frequencies contain normal modes consisting of two hindered rotations of the C≡C diatomic part of the molecule, suggesting that the dissociation of this diatomic from the bicylobutane moiety may be important in the chemistry of this molecule. Because the putative C2v minimum is predicted to have a significant dipole moment (2.6 D), benzvalyne may be detectable in TMC-1, where the related o-benzyne molecule has recently been observed by radio astronomy. [ABSTRACT FROM AUTHOR]

The absence of a reliable formulation of the kinetic energy density functional has hindered the development of orbital free density functional theory. Using the data-aided learning paradigm, we propose a simple prescription to accurately model the kinetic energy density of any system. Our method relies on a dictionary of functional forms for local and nonlocal contributions, which have been proposed in the literature, and the appropriate coefficients are calculated via a linear regression framework. To model the nonlocal contributions, we explore two new nonlocal functionals—a functional that captures fluctuations in electronic density and a functional that incorporates gradient information. Since the analytical functional forms of the kernels present in these nonlocal terms are not known from theory, we propose a basis function expansion to model these seemingly difficult nonlocal quantities. This allows us to easily reconstruct kernels for any system using only a few structures. The proposed method is able to learn kinetic energy densities and total kinetic energies of molecular and periodic systems, such as H2, LiH, LiF, and a one-dimensional chain of eight hydrogens using data from Kohn–Sham density functional theory calculations for only a few structures. [ABSTRACT FROM AUTHOR]

The experimental and theoretical determination of the mean excitation energy, I(0), and the stopping power, S(v), of a material is of great interest in particle and material physics and radiation therapy. For calculations of I(0), the complete set of electronic transitions in a given basis set is required, effectively limiting such calculations to systems with a small number of electrons, even at the random-phase approximation (RPA)/time-dependent Hartree–Fock (TDHF) or time-dependent density-functional theory level. To overcome such limitations, we present here the implementation of a Lanczos algorithm adapted for the paired RPA/TDHF eigenvalue problem in the Dalton program and show that it provides good approximation of the entire RPA eigenspectra in a reduced space. We observe rapid convergence of I(0) with the number of Lanczos vectors as the algorithm favors the transitions with large contributions. In most cases, the algorithm recovers RPA I(0) values of up to 0.5% accuracy at less than a quarter of the full space size. The algorithm not only exploits the RPA paired structure to save computational resources but also preserves certain sum-over-states properties, as first demonstrated by Johnson et al. [Comput. Phys. Commun. 120, 155 (1999)]. The block Lanczos RPA solver, as presented here, thus shows promise for computing mean excitation energies for systems larger than what was computationally feasible before. [ABSTRACT FROM AUTHOR]

The principles of density-functional theory are studied for finite lattice systems represented by graphs. Surprisingly, the fundamental Hohenberg–Kohn theorem is found void, in general, while many insights into the topological structure of the density-potential mapping can be won. We give precise conditions for a ground state to be uniquely v-representable and are able to prove that this property holds for almost all densities. A set of examples illustrates the theory and demonstrates the non-convexity of the pure-state constrained-search functional. [ABSTRACT FROM AUTHOR]

Derivatives of BODIPY are popular fluorophores due to their synthetic feasibility, structural rigidity, high quantum yield, and tunable spectroscopic properties. While the characteristic absorption maximum of BODIPY is at 2.5 eV, combinations of functional groups and substitution sites can shift the peak position by ±1 eV. Time-dependent long-range corrected hybrid density functional methods can model the lowest excitation energies offering a semi-quantitative precision of ±0.3 eV. Alas, the chemical space of BODIPYs stemming from combinatorial introduction of—even a few dozen—substituents is too large for brute-force high-throughput modeling. To navigate this vast space, we select 77 412 molecules and train a kernel-based quantum machine learning model providing < 2 % hold-out error. Further reuse of the results presented here to navigate the entire BODIPY universe comprising over 253 giga (253 × 109) molecules is demonstrated by inverse-designing candidates with desired target excitation energies. [ABSTRACT FROM AUTHOR]

Simulating solids with quantum chemistry methods and Gaussian-type orbitals (GTOs) has been gaining popularity. Nonetheless, there are few systematic studies that assess the basis set incompleteness error (BSIE) in these GTO-based simulations over a variety of solids. In this work, we report a GTO-based implementation for solids and apply it to address the basis set convergence issue. We employ a simple strategy to generate large uncontracted (unc) GTO basis sets that we call the unc-def2-GTH sets. These basis sets exhibit systematic improvement toward the basis set limit as well as good transferability based on application to a total of 43 simple semiconductors. Most notably, we found the BSIE of unc-def2-QZVP-GTH to be smaller than 0.7 mEh per atom in total energies and 20 meV in bandgaps for all systems considered here. Using unc-def2-QZVP-GTH, we report bandgap benchmarks of a combinatorially designed meta-generalized gradient approximation (mGGA) functional, B97M-rV, and show that B97M-rV performs similarly (a root-mean-square-deviation of 1.18 eV) to other modern mGGA functionals, M06-L (1.26 eV), MN15-L (1.29 eV), and Strongly Constrained and Appropriately Normed (SCAN) (1.20 eV). This represents a clear improvement over older pure functionals such as local density approximation (1.71 eV) and Perdew-Burke-Ernzerhof (PBE) (1.49 eV), although all these mGGAs are still far from being quantitatively accurate. We also provide several cautionary notes on the use of our uncontracted bases and on future research on GTO basis set development for solids. [ABSTRACT FROM AUTHOR]

The electric g-tensor is a central quantity for the interpretation of electron paramagnetic resonance spectra. In this paper, a detailed derivation of the 1-electron contributions to the g-tensor is presented in the framework of linear response theory and the second-order Douglas-Kroll-Hess (DKH) transformation. Importantly, the DKH transformation in the presence of a magnetic field is not unique. Whether or not the magnetic field is included in the required Foldy-Wouthuysen transformation, different transformation matrices and, consequently, Hamiltonians result. In this paper, a detailed comparison of both approaches is presented, paying particular attention to the mathematical properties of the resulting Hamiltonians. In contrast to previous studies that address the g-tensor in the framework of DKH theory, the resulting terms are compared to those of the conventional Pauli theory and are given a physical interpretation. Based on these mathematical and physical arguments, we establish that the proper DKH transformation for systems with constant magnetic fields is based on a gauge-invariant Foldy-Wouthuysen transformation, i.e., a Foldy-Wouthuysen transformation including the magnetic field. Calculations using density functional theory (DFT) are carried out on a set of heavy, diatomic molecules, and a set of transition-metal complexes. Based on these calculations, the performance of the relativistic calculation with and without inclusion of picture-change effects is compared. Additionally, the g-tensor is calculated for the Lanthanide dihydrides. Together with the results from the other two molecular test sets, these calculations serve to quantify the magnitude of picture-change effects and elucidate trends across the periodic table. [ABSTRACT FROM AUTHOR]

We present a full density-functional-theory study taking into account the van der Waals interactions of a 2D supramolecular network adsorbed on the Si(111)[Square_Root]3x[Square_Root]3R30°-boron surface denoted SiB. We show that, contrarily to the previous calculations [B. Baris, V. Luzet, E. Duverger, Ph. Sonnet, F. Palmino, and F. Chérioux, Angew. Chem., Int. Ed. 50, 4094 (2011)] molecule-molecule interactions are attractive, thanks to van der Waals corrections which are essential to describe such systems. We confirm the importance of the substrate effect to achieve the molecular network on the boron doped silicon surface without covalent bond. Our simulated STM images, calculated in the framework of the bSKAN code, give better agreement with the experimental STM images than those obtained by the integrated LDOS calculations within the Tersoff-Hamann approximation. The tungsten tip presence is essential to retrieve three paired lobes as observed experimentally. The observed protrusions arise from the phenyl arms located above silicon adatoms. [ABSTRACT FROM AUTHOR]

The influence of a water solvent on the optical absorption properties of alizarin is investigated using time-dependent density functional theory (TDDFT). The solvent is modeled at two different levels of theory: a structureless dielectric medium, using a polarizable continuum model, and the explicit inclusion of water molecules which are treated at the same level of theory as the solute. Thermal effects on the photoabsorption spectra are included by combining TDDFT with first principle molecular dynamics. The effect of molecular distortions on the excitation energies and oscillator strengths is analyzed, and a strong correlation between particular structural and optical properties is found. [ABSTRACT FROM AUTHOR]

This paper supersedes previous theoretical approaches to conceptual DFT because it provides a unified and systematic approach to all of the commonly considered formulations of conceptual DFT, and even provides the essential mathematical framework for new formulations. Global, local, and nonlocal chemical reactivity indicators associated with the “closed-system representation” ([Nα,Nβ,να(r),νβ(r)]) of spin-polarized density functional theory (SP-DFT) are derived. The links between these indicators and the ones associated with the “open-system representation” ([μα,μβ,να(r),νβ(r)]) are derived, including the spin-resolved Berkowitz–Parr identity. The Legendre transform to the “density representation” ([ρα(r),ρβ(r)]) is performed, and the spin-resolved Harbola–Chattaraj–Cedillo–Parr identities linking the density representation to the closed-system and open-system representations are derived. Taken together, these results provide the framework for understanding chemical reactions from both the electron-following perspective (using either the closed-system or the open-system representation) and electron-preceding perspective (density representation). A powerful matrix-vector notation is developed; with this notation, identities in conceptual DFT become universal. Specifically, this notation allows the fundamental identities in conventional (spin-free) conceptual DFT, the [Nα,Nβ] representation, and the [N=Nα+Nβ,NS=Nα-Nβ] representation to be written in exactly the same forms. In cases where spin transfer and electron transfer are coupled (e.g., radical+molecule reactions), we believe that the [Nα,Nβ] representation may be more useful than the more common [N,NS] representation. [ABSTRACT FROM AUTHOR]

Finding the underlying mechanism from the statistical properties of an experimental two-state trajectory generated from dynamics in a complex on-off multisubstate kinetic scheme (KS) is the aim of many experiments. Since the data explicitly shows only transitions between substates of different states, information about the KS is lost, resulting in equivalence of KSs, i.e., the occurrence of different KSs that lead to the same data, in a statistical sense. In order to deal with this phenomenon, a canonical (unique) form of reduced dimensions (RD) is built from the data. RD forms are on-off networks with connections only between substates of different states, where the connections usually have nonexponential waiting time probability density functions. In this paper, we give a list of (about 50) relationships between properties of the data, the topology of reduced dimension forms, and features of KSs. Many of these relationships involve symmetries in RD forms, KSs, and the data and irreversible transitions in KSs. These relationships are useful both in theoretical analysis of on-off KSs and in the analysis of the data. [ABSTRACT FROM AUTHOR]

In this paper, double ionization energies (DIEs) of gas-phase atoms and molecules are calculated by energy difference method with density functional theory. To determine the best functional for double ionization energies, we first study 24 main group atoms in the second, third, and fourth periods. An approximation is used in which the electron density is first obtained from a density functional computation with the exchange-correlation potential Vxc known as statistical average of orbital potentials, after which the energy is computed from that density with 59 different exchange-correlation energy functionals Exc. For the 24 atoms, the two best Exc functional providing DIEs with average absolute deviation (AAD) of only 0.25 eV are the Perdew–Burke–Ernzerhof functional modified by Hammer et al. [Phys. Rev. B 59, 6413 (1999)] and one known as the Krieger–Chen–Iafrate–Savin functional modified by Krieger et al. (unpublished). Surprisingly, none of the 20 available hybrid functionals is among the top 15 functionals for the DIEs of the 24 atoms. A similar procedure is then applied to molecules, with opposite results: Only hybrid functionals are among the top 15 functionals for a selection of 29 molecules. The best Exc functional for the 29 molecules is found to be the Becke 1997 functional modified by Wilson et al. [J. Chem. Phys. 115, 9233 (2001)]. With that functional, the AAD from experiment for DIEs of 29 molecules is just under 0.5 eV. If the two suspected values for C2H2 and Fe(CO)5 are excluded, the AAD improves to 0.32 eV. Many other hybrid functionals perform almost as well. [ABSTRACT FROM AUTHOR]

Morphology control is important for practical applications of composite materials that consist of functional polymers and nanoparticles. Toward that end, block copolymers provide useful templates to arrange nanoparticles in the scaffold of self-organized polymer microdomains. This paper reports theoretical predictions for the distribution of nanoparticles in the lamellar structures of symmetric diblock copolymers on the basis of a polymer density functional theory (DFT) and the potential distribution theorem (PDT). The DFT predicts periodic spacing of lamellar structures in good agreement with molecular dynamics simulations. With the polymer structure from DFT as the input, the PDT is used to examine the effects of particle size, surface energy, polymer chain length, and compressibility on the distribution of nanoparticles in the limit of low particle density. It is found that the nanoparticle distribution depends not only on the particle size and surface energy but also on the local structure of the microdomain interface, polymer chain length, and compressibility. The theoretical predictions are compared well with experiments and simulations. [ABSTRACT FROM AUTHOR]

We propose a free-space density functional theory for polymer adsorption. The derivation within the framework of density functional theory leads to the splitting of the intrinsic free energy into an ideal-gas term and a residual term responsible for the intrinsic energy and the nonbonded interactions between monomers, respectively. A more reasonable treatment is adopted for the residual free energy to count for the monomer-monomer correlation underestimated by the local density approximation. An approach using propagators is proposed to calculate the single-chain partition function and the segment-density distributions, the three adsorption conformations as trains, loops, and tails are further described by propagators. Dirac's bra-ket notation used makes the derivation simpler and provides clearer physical meanings. The theoretical calculations for the adsorption of hard-sphere chains onto a nonadsorbing and an adsorbing hard wall show that the structure of the adsorption layer is strongly affected by the packing effect which has been underestimated by the previous lattice adsorption theory. [ABSTRACT FROM AUTHOR]

The ground state of odd-N polyynic oligomers CNH2 features kink solitons in carbon-carbon bond-length alternation (BLA) patterns. We perform a systematic first-principles computational study of neutral and singly charged kinks in long oligomers addressing relationships between BLA patterns, electron energy gaps, and accompanying distributions of spin and charge densities, both in vacuum and in the screening solvent environment. A quantitative comparison is made of the results derived with four different ab initio methods: from pure density-functional theory to pure Hartree–Fock (HF) and including two popular hybrid density functionals, B3LYP and BHandHLYP. A clear correlation is demonstrated between the derived spatial extent of kinks and the amount of HF exchange used in the functionals. For charged kinks, we find a substantial difference in the behavior of charge and spin densities. [ABSTRACT FROM AUTHOR]

The recently proposed constrained-pairing mean-field theory (CPMFT) is here extended to include exchange and correlation effects from density functional theory (DFT) via alternative densities. We transform from α and β spin densities to alternatives based on the total and on-top pair densities. This transformation is needed because CPMFT produces correct spin-symmetry and space-symmetry adapted densities that traditional DFT functionals are not designed to work with. The inclusion of DFT exchange and correlation effects in CPMFT is well founded both on practical and methodological reasons. We present multiple benchmarks showing that in many cases our model accurately reproduces unrestricted hybrid functional energies (both regular and range separated) and does so on the correct space-symmetry and spin-symmetry surface. Our approach affords efficient inclusion of dynamical correlation effects absent in CPMFT. [ABSTRACT FROM AUTHOR]

Hybrid functionals are responsible for much of the utility of modern Kohn-Sham density functional theory. When rigorously applied to solid-state metallic and small band gap systems, however, the slow decay of their nonlocal Hartree-Fock-type exchange makes hybrids computationally challenging and introduces unphysical effects. This can be remedied by using a range-separated hybrid which only keeps short-range nonlocal exchange, as in the functional of Heyd et al. [J. Chem. Phys. 118, 8207 (2003)]. On the other hand, many molecular properties require full long-range nonlocal exchange, which can also be included by means of a range-separated hybrid such as the recently introduced LC-ωPBE functional [O. A. Vydrov and G. E. Scuseria, J. Chem. Phys. 125, 234109 (2006)]. In this paper, we show that a three-range hybrid which mainly includes middle-range Hartree-Fock-type exchange and neglects long- and short-range Hartree-Fock-type exchange yields excellent accuracy for thermochemistry, barrier heights, and band gaps, emphasizing that the middle-range part of the 1/r potential seems crucial to accurately model these properties. [ABSTRACT FROM AUTHOR]

In Part I [R. E. Beckham and M. A. Bevan, J. Chem. Phys. 127, 164708 (2007)], results were presented for the sedimentation equilibrium of concentrated colloidal dispersions using confocal scanning laser microscopy experiments, Monte Carlo (MC) simulations, and a local density approximation perturbation theory. In this paper, we extended the modeling effort on those systems to include nonlocal density functional theory (DFT), which is capable of predicting the microstructure of the sediment at length scales comparable to the colloidal particle dimension. Specifically, we use a closure-based DFT formulation to predict interfacial colloidal sedimentation equilibrium density profiles. The colloid-colloid and colloid-surface interactions were modeled with DLVO screened electrostatic potentials using parameters taken directly from the experimental work. The DFT profiles were compared to the experimental and MC results from Part I. Good agreement was found for relatively dilute interfacial colloidal fluids, but agreement was less satisfactory as interfacial layering became more pronounced for conditions approaching the onset of interfacial crystallization. We also applied DFT in an inverse sense using the measured colloid density profile to extract the underlying colloid-surface potential; this can be thought of as a microscopic analog to the well-known procedure of using the macroscopic (coarse-grained) density profile to extract the osmotic equation of state. For the dilute interfacial fluid, the inverse DFT calculations reproduced the true colloid-surface potential to within 0.5kT at all elevations. [ABSTRACT FROM AUTHOR]

In this paper we employ all-electron ab initio time-dependent density functional theory based method to calculate the long range dipole-dipole dispersion coefficient (van der Waals coefficient) C6 of sodium atom clusters containing even number of atoms ranging from 2 to 20 atoms. The dispersion coefficients are obtained via Casimir-Polder relation [Phys. Rev. 3, 360 (1948)]. The calculations are carried out with two different exchange-correlation potentials: (i) the asymptotically correct statistical average of orbital potential (SAOP) and (ii) Vosko-Wilk-Nusair representation [Can. J. Phys. 58, 1200 (1980)] of exchange-correlation potential within local density approximation. A comparison with the other theoretical results has been performed. We also present the results for the static polarizabilities of sodium clusters and also compare them with other theoretical and experimental results. These comparisons reveal that the SAOP results for C6 and static polarizability are quite accurate and very close to the experimental results. We examine the relationship between volume of the cluster and van der Waals coefficient, and find that to a very high degree of correlation C6 scales as the square of the volume. We also present the results for van der Waals coefficient corresponding to cluster–Ar atom and cluster–N2 molecule interactions. [ABSTRACT FROM AUTHOR]

This paper presents response theory based on a spin-restricted Kohn-Sham formalism for computation of time-dependent and time-independent nonlinear properties of molecules with a high spin ground state. The developed approach is capable to handle arbitrary perturbations and constitutes an efficient procedure for evaluation of electric, magnetic, and mixed properties. Apart from presenting the derivation of the proposed approach, we show results from illustrating calculations of static and dynamic hyperpolarizabilities of small Si(3n+1)H(6n+3) (n=0,1,2) clusters which mimic Si(111) surfaces with dangling bond defects. The results indicate that the first hyperpolarizability tensor components of Si(3n+1)H(6n+3) have an ordering compatible with the measurements of second harmonic generation in SiO2/Si(111) interfaces and, therefore, support the hypothesis that silicon surface defects with dangling bonds are responsible for this phenomenon. The results exhibit a strong dependence on the quality of basis set and exchange-correlation functional, showing that an appropriate set of diffuse functions is required for reliable predictions of the first hyperpolarizability of open-shell compounds. [ABSTRACT FROM AUTHOR]