Analysis and modelling of atomic and molecular kohn-sham potentials. In: Nalewajski RF, editor. Density Functional Theory I. Berlin: Springer Berlin Heidelberg; 1996. pp. 107-67. [DOI: 10.1007/3-540-61091-x_4]

Chemical potential, derivative discontinuity, fractional electrons, jump of the Kohn-Sham potential, atoms as thermodynamic open systems, and other (mis)conceptions of the density functional theory of electrons in molecules

Many references exist in the density functional theory (DFT) literature to the chemical potential of the electrons in an atom or a molecule. The origin of this notion has been the identification of the Lagrange multiplier μ = ∂E/∂N in the Euler-Lagrange variational equation for the ground state density as the chemical potential of the electrons. We first discuss why the Lagrange multiplier in this case is an arbitrary constant and therefore cannot be a physical characteristic of an atom or molecule. The switching of the energy derivative ("chemical potential") from -I to -A when the electron number crosses the integer, called integer discontinuity or derivative discontinuity, is not physical but only occurs when the nonphysical noninteger electron systems and the corresponding energy and derivative ∂E/∂N are chosen in a specific discontinuous way. The question is discussed whether in fact the thermodynamical concept of a chemical potential can be defined for the electrons in such few-electron systems as atoms and molecules. The conclusion is that such systems lack important characteristics of thermodynamic systems and do not afford the definition of a chemical potential. They also cannot be considered as analogues of the open systems of thermodynamics that can exchange particles with an environment (a particle bath or other members of a Gibbsian ensemble). Thermodynamical (statistical mechanical) concepts like chemical potential, open systems, grand canonical ensemble etc. are not applicable to a few electron system like an atom or molecule. A number of topics in DFT are critically reviewed in light of these findings: jumps in the Kohn-Sham potential when crossing an integer number of electrons, the band gap problem, the deviation-from-straight-lines error, and the role of ensembles in DFT.

X-ray Raman scattering: a new in situ probe of molecular structure during nucleation and crystallization from liquid solutions

X-ray Raman scattering (XRS) has been used for in situ probing of solute molecule speciation in solution during cooling crystallization.

Accurate and efficient approximation to the optimized effective potential for exchange

We devise an efficient practical method for computing the Kohn-Sham exchange-correlation potential corresponding to a Hartree-Fock electron density. This potential is almost indistinguishable from the exact-exchange optimized effective potential (OEP) and, when used as an approximation to the OEP, is vastly better than all existing models. Using our method one can obtain unambiguous, nearly exact OEPs for any reasonable finite one-electron basis set at the same low cost as the Krieger-Li-Iafrate and Becke-Johnson potentials. For all practical purposes, this solves the long-standing problem of black-box construction of OEPs in exact-exchange calculations.

Prediction of iron K-edge absorption spectra using time-dependent density functional theory

Iron K-edge X-ray absorption pre-edge features have been calculated using a time-dependent density functional approach. The influence of functional, solvation, and relativistic effects on the calculated energies and intensities has been examined by correlation of the calculated parameters to experimental data on a series of 10 iron model complexes, which span a range of high-spin and low-spin ferrous and ferric complexes in O(h) to T(d) geometries. Both quadrupole and dipole contributions to the spectra have been calculated. We find that good agreement between theory and experiment is obtained by using the BP86 functional with the CP(PPP) basis set on the Fe and TZVP one of the remaining atoms. Inclusion of solvation yields a small improvement in the calculated energies. However, the inclusion of scalar relativistic effects did not yield any improved correlation with experiment. The use of these methods to uniquely assign individual spectral transitions and to examine experimental contributions to backbonding is discussed.

Describing static correlation in bond dissociation by Kohn–Sham density functional theory

We show that density functional theory within the RPA (random phase approximation for the exchange-correlation energy) provides a correct description of bond dissociation in H(2) in a spin-restricted Kohn-Sham formalism, i.e., without artificial symmetry breaking. We present accurate adiabatic connection curves both at equilibrium and beyond the Coulson-Fisher point. The strong curvature at large bond length implies important static (left-right) correlation, justifying modern hybrid functional constructions but also demonstrating their limitations. Although exact at infinite separation and accurate near the equilibrium bond length, the RPA dissociation curve displays unphysical repulsion at larger but finite bond lengths. Going beyond the RPA by including the exact exchange kernel (RPA+X), we find a similar repulsion. We argue that this deficiency is due to the absence of double excitations in adiabatic linear response theory. Further analyzing the H(2) dissociation limit we show that the RPA+X is not size consistent, in contrast to the RPA.