The Influence of the Exchange-Correlation Functional on the Non-Interacting Kinetic Energy and Its Implications for Orbital-Free Density Functional ApproximationsSCI北大核心CSCD

In this work it is shown that the kinetic energy and the exchange-correlation energy are mutual dependent on each other.This aspect is first derived in an orbital-free context.It is shown that the total Fermi potential depends on the density only,the individual parts,the Pauli kinetic energy and the exchange-correlation energy,however,are orbital dependent and as such mutually influence each other.The numerical investigation is performed for the orbital-based non-interacting Kohn-Sham system in order to avoid additional effects due to further approximations of the kinetic energy.The numerical influence of the exchange-correlation functional on the non-interacting kinetic energy is shown to be of the orderof a few Hartrees.For chemical purposes,however,the energetic performance as a function of the nuclear coordinates is much more important than total energies.Therefore,the effect on the bond dissociation curve was studied exemplarily for the carbon monoxide.The data reveals that,the mutual influence between the exchange-correlation functional and the kinetic energy has a significant influence on bond dissociation energies and bond distances.Therefore,the effect of the exchange-correlation treatment must be considered in the design of orbital-free density functional approximations for the kinetic energy.     

Effective homogeneity of the exchange-correlation and non-interacting kinetic energy functionals under density scaling

Correlated electron densities, experimental ionisation potentials, and experimental electron affinities are used to investigate the homogeneity of the exchange-correlation and non-interacting kinetic energy functionals of Kohn-Sham density functional theory under density scaling. Results are presented for atoms and small molecules, paying attention to the influence of the integer discontinuity and the choice of the electron affinity. For the exchange-correlation functional, effective homogeneities are highly system-dependent on either side of the integer discontinuity. By contrast, the average homogeneity-associated with the potential that averages over the discontinuity-is generally close to 4/3 when the discontinuity is computed using positive affinities for systems that do bind an excess electron and negative affinities for those that do not. The proximity to 4/3 becomes increasingly pronounced with increasing atomic number. Evaluating the discontinuity using a zero affinity in systems that do not bind an excess electron instead leads to effective homogeneities on the electron abundant side that are close to 4/3. For the non-interacting kinetic energy functional, the effective homogeneities are less system-dependent and the effect of the integer discontinuity is less pronounced. Average values are uniformly below 5/3. The study provides information that may aid the development of improved exchange-correlation and non-interacting kinetic energy functionals.     

On the evaluation of the non-interacting kinetic energy in density functional theory

The utility of both an orbital-free and a single-orbital expression for computing the non-interacting kinetic energy in density functional theory is investigated for simple atomic systems. The accuracy of both expressions is governed by the extent to which the Kohn-Sham equation is solved for the given exchange-correlation functional and so special attention is paid to the influence of finite Gaussian basis sets. The orbital-free expression is a statement of the virial theorem and its accuracy is quantified. The accuracy of the single-orbital expression is sensitive to the choice of Kohn-Sham orbital. The use of particularly compact orbitals is problematic because the failure to solve the Kohn-Sham equation exactly in regions where the orbital has decayed to near-zero leads to unphysical behaviour in regions that contribute to the kinetic energy, rendering it inaccurate. This problem is particularly severe for core orbitals, which would otherwise appear attractive due to their formally nodeless nature. The most accurate results from the single-orbital expression are obtained using the relatively diffuse, highest occupied orbitals, although special care is required at orbital nodes.     

Some formal properties of ensemble density functionals

Formal properties of ensemble density functionals are examined. Expressions for the difference between energy functionals where the particle number differs by one are constructed in terms of their first functional derivatives for the universal energy functional, the electron–electron repulsion energy functional, and the interacting kinetic energy functional. Equations that must be satisfied by second and higher order functional derivatives are derived. It is also shown that the shape of ${\delta V_{ee}[\rho]\over\delta\rho({\bf r})}$ and ${\delta K[\rho]\over\delta\rho({\bf r})}$ , the functional derivatives of the mutual electron–electron repulsion, and kinetic energy, respectively, are separately particle number independent for particle numbers between successive integers. © 2013 Wiley Periodicals, Inc.     

Phase Space View of Ensembles of Excited States

The density functional theory and its extension to ensembles of excited states can be formalized as thermodynamics. However, these theories are not unique because one of their key quantities, the kinetic energy density,can be defined in several ways. Usually, the everywhere positive gradient form is applied; however, other forms are also acceptable, provided they integrate to the true kinetic energy. Recently, a kinetic energy density of the ground-state theory maximizing the information entropy has been proposed. Here, ensemble kinetic energy density, leading to extremum information entropy, is derived via constrained search. The corresponding ensemble temperature is found to be constant.     

Transport of matter and energy in a mesoscopic thermo-hydrodynamic approach

We derive a thermo-hydrodynamic theory for particles and energy flow, based on a nonequilibrium grand-canonical ensemble formalism. The time-dependent kinetic coefficients are explicitly given in terms of microscopic mechanical quantities. The time evolution equations describing the coupled flow of energy and particles are derived. The second-rank tensorial fluxes of current of energy and particles present in the nonequilibrium ensemble are nondiagonal. We obtain a generalized Fick's Law, which presents the effect of the energy flow on the particle diffusion equation.     

On the single-electron local kinetic energy

Ｉｎｒｅｃｅｎｔｄｅｃａｄｅｓ,ｇｒｅａｔｐｒｏｇｒｅｓｓｈａｓｂｅｅｎｍａｄｅｉｎｔｈｅｄｅｎｓｉｔｙｆｕｎｃｔｉｏｎａｌｔｈｅｏｒｙａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｓ.ＨｏｈｅｎｂｅｒｇＫｏｈｎｔｈｅｏｒｅｍｉｓｉｔｓｂａｓｉｓ.Ｔａｋｉｎｇｔｈｅｏｎｅｅｌｅｃｔｒｏｎｄｅｎｓｉｔｙａｓａｆｕｎｄａｍｅｎｔａｌｖａｒｉａｂｌｅ,ｔｈｉｓｔｈｅｏｒｙｔｒｅａｔｓｔｈｅｔｏｔａｌｅｎｅｒｇｙＥ[ρ]ｏｆｔｈｅｅｌｅｃｔｒｏｎｉｃｇｒｏｕｎｄｓｔａｔｅｆｏｒａｍｏｌｅｃｕｌｅａｓａｆｕｎｃｔｉｏｎａｌ…     

Nuclear magnetic resonance shielding tensors calculated with kinetic energy density-dependent exchange-correlation functionals

Exchange-correlation functionals that depend on the non-interacting kinetic energy density (τ) break gauge invariance in the presence of a magnetic field. This yields incorrect results for molecular magnetic properties. We propose a simple generalization of the kinetic energy density that resolves this problem. Our modification is validated by computing NMR absolute isotropic shielding constants and shielding anisotropies with τ-dependent functionals for a representative set of molecules. The accuracy of τ-dependent approximations is found to surpass that of common generalized gradient approximations (GGA) and hybrid functionals for strongly deshielded nuclei.     

Self-consistent meta-generalized gradient approximation study of adsorption of aromatic molecules on noble metal surfaces

The adsorption of benzene, pyridine, and two nucleobases on the Au(111) surface has been investigated using a fully relaxed, self-consistent meta-generalized gradient approximation (meta-GGA) density functional theory setup with the M06-L functional. The meta-GGA based molecule-surface separations are shortened and the adsorption bond strengths of the molecules are greatly improved over the virtually non-interacting results obtained when using a plain GGA exchange-correlation functional. The nucleobases containing oxygen atoms show higher corrugation with adsorption site and orientation than the other aromatic molecules considered. The adsorption of pentacene is studied on Au, Ag, and Cu surfaces. In agreement with experiment, the adsorption energies are found to increase with decreasing nobleness, but the dependency is underestimated. We point out how the kinetic energy density can discriminate between covalent and non-covalent bonding regions of orbital overlap.     

Method of local-scaling transformations and density functional theory in quantum chemistry. III. The energy density functional: Spin-restricted approach

The rigorous derivation of the energy density functional is proposed within the framework of the spinfree, or spin-restricted formulation of the energy density functional theory. It is shown particularly that the kinetic energy density functional is given by a sum of the Weizsacker term and the so-called “modified” Thomas–Fermi one. The variational principle is formulated for the energy density functional theory in terms of the Euler–Lagrange equation, and the virial theorem is proposed.     

Study on kinetics of combustion of brick-shaped carbonaceous materials

Combustion of brick-shaped carbonaceous materials (carbon deposits from coke oven, coke and electrographite) was carried out in thermobalance in static air. Analysis of kinetics of the process was carried out using both classical (Arrhenius law) and newer (three-parametric equation) methods. In classical approach two types of kinetic equations were used in calculations: differential and integral. The results obtained show that, independently on kinetic variables (α – conversion degree or m – mass of sample) used in differential equations, kinetics of combustion of brick-shaped carbonaceous materials is characterized by only one pair of Arrhenius coefficients: activation energy (E) and pre-exponential constant (A). At the same time the integral equation demonstrates distinction in relation to methods based on differential equations, generating higher activation energies and separate isokinetic effect (IE). Parallel IE shows that kinetic analysis has to encompass activation energy in connection to second coefficient, pre-exponential constant A, depending on assumptions made for kinetic equations. On the other hand three-parametric equation allows describing kinetic of combustion in alternative way using only one experimental value – initial temperature in form of point of initial oxidation (PIO) – and also offers new methods of interpretation of the process.     

Differential geometry based solvation model. III. Quantum formulation

Solvation is of fundamental importance to biomolecular systems. Implicit solvent models, particularly those based on the Poisson-Boltzmann equation for electrostatic analysis, are established approaches for solvation analysis. However, ad hoc solvent-solute interfaces are commonly used in the implicit solvent theory. Recently, we have introduced differential geometry based solvation models which allow the solvent-solute interface to be determined by the variation of a total free energy functional. Atomic fixed partial charges (point charges) are used in our earlier models, which depends on existing molecular mechanical force field software packages for partial charge assignments. As most force field models are parameterized for a certain class of molecules or materials, the use of partial charges limits the accuracy and applicability of our earlier models. Moreover, fixed partial charges do not account for the charge rearrangement during the solvation process. The present work proposes a differential geometry based multiscale solvation model which makes use of the electron density computed directly from the quantum mechanical principle. To this end, we construct a new multiscale total energy functional which consists of not only polar and nonpolar solvation contributions, but also the electronic kinetic and potential energies. By using the Euler-Lagrange variation, we derive a system of three coupled governing equations, i.e., the generalized Poisson-Boltzmann equation for the electrostatic potential, the generalized Laplace-Beltrami equation for the solvent-solute boundary, and the Kohn-Sham equations for the electronic structure. We develop an iterative procedure to solve three coupled equations and to minimize the solvation free energy. The present multiscale model is numerically validated for its stability, consistency and accuracy, and is applied to a few sets of molecules, including a case which is difficult for existing solvation models. Comparison is made to many other classic and quantum models. By using experimental data, we show that the present quantum formulation of our differential geometry based multiscale solvation model improves the prediction of our earlier models, and outperforms some explicit solvation model.     

Integrals and derivatives for correlated Gaussian functions using matrix differential calculus

The matrix differential calculus is applied for the first time to a quantum chemical problem via new matrix derivations of integral formulas and gradients for Hamiltonian matrix elements in a basis of correlated Gaussian functions. Requisite mathematical background material on Kronecker products, Hadamard products, the vec and vech operators, linear structures, and matrix differential calculus is presented. New matrix forms for the kinetic and potential energy operators are presented. Integrals for overlap, kinetic energy, and potential energy matrix elements are derived in matrix form using matrix calculus. The gradient of the energy functional with respect to the correlated Gaussian exponent matrices is derived. Burdensome summation notation is entirely replaced with a compact matrix notation that is both theoretically and computationally insightful. © 1996 John Wiley & Sons, Inc.     

Non-local relation between kinetic and exchange energy densities in Hartree–Fock theory

A non-local generalization K( r, r' ) of the kinetic energy t( r ) such that t( r ) = ∫K( r, r' ) dr' is defined using the idempotency property of the Hartree–Fock first-order density matrix. This is, in turn, related by means of an explicit differential equation to the non-local exchange energy density X( r, r' ). The relationship is illustrated for a couple of examples: with the Fermi-hole in a uniform electron gas, of importance in the local density version of density functional theory, and with inhomogeneous electron systems.     

Determination of energetics and kinetics from single-particle intermittency and ensemble-averaged fluorescence intensity decay of quantum dots

Quantification of energetics and kinetics for the band-edge exciton states of quantum dots and the long-lived dark state is important for better understanding of the underlying mechanism for single-particle intermittency and ensemble fluorescence intensity decay. Based on a multistate diffusion-reaction model by extending our previous studies, we analyze experimental data from ensemble measurements and fluorescence intermittency of single quantum dots and determine important molecular-based quantities such as Stokes shift, free energy gap, activation energy, reorganization energy, and other kinetic parameters.     

Periodic boundary condition induced breakdown of the equipartition principle and other kinetic effects of finite sample size in classical hard-sphere molecular dynamics simulation

We examine consequences of the non-Boltzmann nature of probability distributions for one-particle kinetic energy, momentum, and velocity for finite systems of classical hard spheres with constant total energy and nonidentical masses. By comparing two cases, reflecting walls (NVE or microcanonical ensemble) and periodic boundaries (NVEPG or molecular dynamics ensemble), we describe three consequences of the center-of-mass constraint in periodic boundary conditions: the equipartition theorem no longer holds for unequal masses, the ratio of the average relative velocity to the average velocity is increased by a factor of [N/(N-1)]1/2, and the ratio of average collision energy to average kinetic energy is increased by a factor of N/(N-1). Simulations in one, two, and three dimensions confirm the analytic results for arbitrary dimension.     

Kinetic Analysis of Non-isothermal DSC Data. Computer-aided test of its applicability

The applicability of the kinetic analysis of data obtained by non-isothermal differential scanning calorimetry (DSC) is discussed. The Johnson-Mehl-Avrami (JMA) model was used for the computer simulation of DSC traces subsequently analysed by common methods of kinetic analysis of non-isothermal data. For the temperature-independent kinetic exponent n of the JMA equation, the kinetic analysis was shown to provide correct results, e.g. a correct kinetic model and apparent activation energy. On the other hand, for the temperature-dependent kinetic exponent, there is a great possibility of erroneous determination of the correct kinetic model and apparent activation energy, especially at higher heating rates. Since the temperature dependence of n cannot be determined on the basis of non-isothermal DSC experiments, conclusions must be drawn with appropriate caution. This revised version was published online in August 2006 with corrections to the Cover Date.     

3,7-二硝基亚氨基-2,4-二氮杂双环[3,3,0]辛烷的放热分解反应动力学

The kinetic parameters of the exothermic decomposition reaction of the title compound in a temperature-programmed mode have been studied by means of DSC. The DSC data obtained are fitted to the integral, differential and exothermic rate equations by linear least-squares, iterative, combined dichotomous and least-squares methods, respectively. After establishing the most probable general expression of differential and integral mechanism functions by the logical choice method, the corresponding values of the apparent activation energy (Ea), pre-exponential factor (A) and reaction order (n) will be obtained by the exothermic rate equation. The results show that the empirical kinetic model function in differential form and the values of Ea and A of this reaction are (1-α)0.44, 230.4 kJ/mol and 1018.16 s-1, respectively. With the help of the heating rate and obtained kinelic parameters, the kinetic equation of the exothermic decomposition reaction process of the title compound is proposed. The critical temperature of thermal explosion of the compound is 302.6℃. The above-mentioned kinetic parameters are quite useful for analyzing and evaluating the stability and thermal change rule of the title compound.     

Exact exchange-only density functional theory by means of local scaling transformations

The density-constrained variation of the kinetic energy of a non-interacting system carried out within the framework of local-scaling transformations is employed in the present work for the purpose of determining exchange-only Kohri Sham orbitals for the beryllium atom as well as the Kohn-Sham exchange-only correlation potential. The starting basis functions are those of Clementi-Roetti and Raffenetti-types. We have also performed an optimization via a densitydriven method. The resulting exchange-only potentials for the examples studied are almost indistinguishable from the Talman-Shadwick potential for the beryllium atom.