Density functional study of double ionization energies

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 V xc known as statistical average of orbital potentials, after which the energy is computed from that density with 59 different exchange-correlation energy functionals E xc. For the 24 atoms, the two best E xc 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 E xc 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.3(2) eV. Many other hybrid functionals perform almost as well.     

Total energy evaluation in the Strutinsky shell correction method

We analyze the total energy evaluation in the Strutinsky shell correction method (SCM) of Ullmo et al. [Phys. Rev. B 63, 125339 (2001)], where a series expansion of the total energy is developed based on perturbation theory. In agreement with Yannouleas and Landman [Phys. Rev. B 48, 8376 (1993)], we also identify the first-order SCM result to be the Harris functional [Phys. Rev. B 31, 1770 (1985)]. Further, we find that the second-order correction of the SCM turns out to be the second-order error of the Harris functional, which involves the a priori unknown exact Kohn-Sham (KS) density, rho(KS)(r). Interestingly, the approximation of rho(KS)(r) by rho(out)(r), the output density of the SCM calculation, in the evaluation of the second-order correction leads to the Hohenberg-Kohn-Sham functional. By invoking an auxiliary system in the framework of orbital-free density functional theory, Ullmo et al. designed a scheme to approximate rho(KS)(r), but with several drawbacks. An alternative is designed to utilize the optimal density from a high-quality density mixing method to approximate rho(KS)(r). Our new scheme allows more accurate and complex kinetic energy density functionals and nonlocal pseudopotentials to be employed in the SCM. The efficiency of our new scheme is demonstrated in atomistic calculations on the cubic diamond Si and face-centered-cubic Ag systems.     

Dipole moments from atomic-number-dependent potentials in analytic density-functional theory

Molecular dipole moments of analytic density-functional theory are investigated. The effect of element-dependent exchange potentials on these moments are examined by comparison with conventional quantum-chemical methods and experiment for the subset of the extended G2 set of molecules that have nonzero dipole moment. Fitting the Kohn-Sham [Phys. Rev. 140, A1133 (1965)] potential itself makes a mean absolute error of less than 0.1 D. Variation of alpha (Slater's [Phys. Rev. 81, 385 (1951)] exchange parameter) values has far less effect on dipole moments than on energies. It is argued that in variable alpha methods one should choose the smaller of the two rather than the geometric mean of the two alpha values for the heteroatomic part of the linear-combination-atomic-orbital density. Calculations on the dipole moment of NH(2)(CH)(24)NO(2) are consistent with earlier calculations and show that varying the differences between alpha values for atoms with different atomic numbers has only short-ranged electrostatic effects.     

Spherically and system-averaged pair density functional theory

In a couple of recent papers Gori-Giorgi and Savin [Phys. Rev. A 71, 032513 (2005)] proposed a theory that provides simple radial equations to generate the spherically and system averaged pair density. In a recent density matrix functional theory [A. Nagy, Phys. Rev. A 66, 022505 (2002)] it was shown that the problem of an arbitrary system can be reduced to a two-particle problem. Based on this theory, via a double adiabatic connection, it is rigorously derived that the square root of the ground-state spherically and system averaged pair density is the solution of a simple radial equation, that is, contrary to the theory of Gori-Giorgi and Savin only a single equation has to be considered.     

The role of exchange-correlation functionals in the potential energy surface and dynamics of N(2) dissociation on W surfaces

We study the dissociative adsorption of N(2) on W(100) and W(110) by means of density functional theory and classical dynamics. Working with a full six-dimensional adiabatic potential energy surface (PES), we find that the theoretical results of the dynamical problem strongly depend on the choice of approximate exchange-correlation functional for the determination of the PES. We consider the Perdew-Wang-91 [Perdew et al., Phys. Rev. B 46, 6671 (1992)] and Perdew-Burke-Ernzerhof (RPBE) [Hammer et al., Phys. Rev. B 59, 7413 (1999)] functionals and carry out a systematic comparison between the dynamics determined by the respective PESs. Even though it has been shown in earlier works that the RPBE may provide better values for the chemisorption energies, our study brings evidence that it gives rise to a PES with excessive repulsion far from the surface.     

Density functional theory of fluids in nanopores: analysis of the fundamental measures theory in extreme dimensional-crossover situations

Two density functional theories, the fundamental measures theory of Rosenfeld [Phys. Rev. Lett. 63, 980 (1989)] and a subsequent approximation by Tarazona [Phys. Rev. Lett. 84, 694 (2000)] are applied to the study of the hard-sphere fluid in two situations: the cylindrical pore and the spherical cavity. The results are compared with those obtained with grand canonical ensemble Monte Carlo simulations. The differences between both theories are evaluated and interpreted in the terms of the dimensional crossover from three to one and zero dimensions.     

Density functional theory of liquid crystals and surface anchoring: hard Gaussian overlap-sphere and hard Gaussian overlap-surface potentials

This article applies the density functional theory to confined liquid crystals, comprised of ellipsoidal shaped particles interacting through the hard Gaussian overlap (HGO) potential. The extended restricted orientation model proposed by Moradi and co-workers [J. Phys.: Condens. Matter 17, 5625 (2005)] is used to study the surface anchoring. The excess free energy is calculated as a functional expansion of density around a reference homogeneous fluid. The pair direct correlation function (DCF) of a homogeneous HGO fluid is approximated, based on the optimized sum of Percus-Yevick and Roth DCF for hard spheres; the anisotropy introduced by means of the closest approach parameter, the expression proposed by Marko [Physica B 392, 242 (2007)] for DCF of HGO, and hard ellipsoids were used. In this study we extend an our previous work [Phys. Rev. E 72, 061706 (2005)] on the anchoring behavior of hard particle liquid crystal model, by studying the effect of changing the particle-substrate contact function instead of hard needle-wall potentials. We use the two particle-surface potentials: the HGO-sphere and the HGO-surface potentials. The average number density and order parameter profiles of a confined HGO fluid are obtained using the two particle-wall potentials. For bulk isotropic liquid, the results are in agreement with the Monte Carlo simulation of Barmes and Cleaver [Phys. Rev. E 71, 021705 (2005)]. Also, for the bulk nematic phase, the theory gives the correct density profile and order parameter between the walls.     

Rapid, accurate calculation of the s-wave scattering length

Transformation of the conventional radial Schro?dinger equation defined on the interval r ∈ [0, ∞) into an equivalent form defined on the finite domain y(r) ∈ [a, b] allows the s-wave scattering length a(s) to be exactly expressed in terms of a logarithmic derivative of the transformed wave function φ(y) at the outer boundary point y = b, which corresponds to r = ∞. In particular, for an arbitrary interaction potential that dies off as fast as 1/r(n) for n ≥ 4, the modified wave function φ(y) obtained by using the two-parameter mapping function r(y; ?r,β) = ?r[1 + 1/β tan(πy/2)] has no singularities, and a(s) = ?r[1 + 2/πβ 1/φ(1) dφ(1)/dy]. For a well bound potential with equilibrium distance r(e), the optimal mapping parameters are ?r ≈ r(e) and β ≈ n/2 - 1. An outward integration procedure based on Johnson's log-derivative algorithm [J. Comp. Phys. 13, 445 (1973)] combined with a Richardson extrapolation procedure is shown to readily yield high precision a(s)-values both for model Lennard-Jones (2n, n) potentials and for realistic published potentials for the Xe-e(-), Cs(2)(aΣ(u)(+)(3)), and (3, 4)He(2)(XΣ(g)(+)(1)) systems. Use of this same transformed Schro?dinger equation was previously shown [V. V. Meshkov et al., Phys. Rev. A 78, 052510 (2008)] to ensure the efficient calculation of all bound levels supported by a potential, including those lying extremely close to dissociation.     

Universal mathematical identities in density functional theory: results from three different spin-resolved representations

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(alpha),N(beta),nu(alpha)(r),nu(beta)(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" ([mu(alpha),mu(beta),nu(alpha)(r),nu(beta)(r)]) are derived, including the spin-resolved Berkowitz-Parr identity. The Legendre transform to the "density representation" ([rho(alpha)(r),rho(beta)(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(alpha),N(beta)] representation, and the [N=N(alpha)+N(beta),N(S)=N(alpha)-N(beta)] 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(alpha),N(beta)] representation may be more useful than the more common [N,N(S)] representation.     

Accelerating the convergence of the total energy evaluation in density functional theory calculations

A special feature of the Strutinsky shell correction method (SCM) [D. Ullmo et al., Phys. Rev. B 63, 125339 (2001)] and the recently proposed orbital-corrected orbital-free density functional theory (OO-DFT) [B. Zhou and Y. A. Wang, J. Chem. Phys. 124, 081107 (2006)] is that the second-order corrections are incorporated in the total energy evaluation. In the SCM, the series expansion of the total electronic energy is essentially the Harris functional with its second-order correction. Unfortunately, a serious technical problem for the SCM is the lack of the exact Kohn-Sham (KS) density rho KS(r) required for the evaluation of the second-order correction. To overcome this obstacle, we design a scheme that utilizes the optimal density from a high-quality density mixing scheme to approximate rho KS(r). Recently, we proposed two total energy density functionals, i.e., the Zhou-Wang-lambda (ZW lambda) and the Wang-Zhou-alpha (WZ alpha) functionals, for use in the OO-DFT method. If the two interpolation parameters, lambda and alpha, are chosen to allow the second-order errors of the ZW lambda and the WZ alpha functionals to vanish, these two functionals reduce to the Hohenberg-Kohn-Sham functional with its second-order correction. Again, the optimal density from a high-quality density mixing scheme is used to approximate rho KS(r) in the evaluation of lambda and alpha. This approach is tested in iterative KS-DFT calculations on systems with different chemical environments and can also be generalized for use in other iterative first-principles quantum chemistry methods.     

Single particle force distributions in simple fluids

The distribution function, W(F), of the magnitude of the net force, F, on particles in simple fluids is considered, which follows on from our previous publication [A. C. Bran?ka, D. M. Heyes, and G. Rickayzen, J. Chem. Phys. 135, 164507 (2011)] concerning the pair force, f, distribution function, P(f), which is expressible in terms of the radial distribution function. We begin by discussing the force on an impurity particle in an otherwise pure fluid but later specialize to the pure fluid, which is studied in more detail. An approximate formula, expected to be valid asymptotically, for W(F) referred to as, W(1)(F) is derived by taking into account only binary spatial correlations in the fluid. It is found that W(1)(F) = P(f). Molecular dynamics simulations of W for the inverse power (IP) and Lennard-Jones potential fluids show that, as expected, W(F) and P(f) agree well in the large force limit for a wide range of densities and potential forms. The force at which the maximum in W(F) occurs for the IP fluids follows a different algebraic dependence with density in low and high density domains of the equilibrium fluid. Other characteristic features in the force distribution functions also exhibit the same trends. An exact formula is derived relating W(F) to P(x)(F(x)), the distribution function of the x-cartesian components of the net force, F(x), on a particle. W(F) and P(x)(F(x)) have the same analytical forms (apart from constants) in the low and high force limits.     

Variationally stable calculations for molecular systems: polarizabilities and two-photon ionization cross section for the hydrogen molecule

The variationally stable method of Gao and Starace [B. Gao and A. F. Starace, Phys. Rev. Lett. 61, 404 (1988); Phys. Rev. A 39, 4550 (1989)] has been applied for the first time to the study of multiphoton processes in molecular systems. The generalization in theory is presented, as well as the calculation of properties such as the static and dynamic polarizabilities of the hydrogen molecule and the generalized two-photon ionization cross section. The Schwinger variational iterative method [R. R. Lucchese and V. McKoy, Phys. Rev. A 21, 112 (1980)] has been applied in the achievement of the photoelectron wave function, while a Hartree-Fock representation has been used for the target. This research has been motivated by the scarceness of ab initio calculations of molecular multiphoton ionization cross sections in the literature.     

Second-order Kohn-Sham perturbation theory: correlation potential for atoms in a cavity

Second-order perturbation theory based on the Kohn-Sham Hamiltonian leads to an implicit density functional for the correlation energy E(c) (MP2), which is explicitly dependent on both occupied and unoccupied Kohn-Sham single-particle orbitals and energies. The corresponding correlation potential v(c) (MP2), which has to be evaluated by the optimized potential method, was found to be divergent in the asymptotic region of atoms, if positive-energy continuum states are included in the calculation [Facco Bonetti et al., Phys. Rev. Lett. 86, 2241 (2001)]. On the other hand, Niquet et al., [J. Chem. Phys. 118, 9504 (2003)] showed that v(c) (MP2) has the same asymptotic -alpha(2r(4)) behavior as the exact correlation potential, if the system under study has a discrete spectrum only. In this work we study v(c) (MP2) for atoms in a spherical cavity within a basis-set-free finite differences approach, ensuring a completely discrete spectrum by requiring hard-wall boundary conditions at the cavity radius. Choosing this radius sufficiently large, one can devise a numerical continuation procedure which allows to normalize v(c) (MP2) consistent with the standard choice v(c)(r-->infinity)=0 for free atoms, without modifying the potential in the chemically relevant region. An important prerequisite for the success of this scheme is the inclusion of very high-energy virtual states. Using this technique, we have calculated v(c) (MP2) for all closed-shell and spherical open-shell atoms up to argon. One finds that v(c) (MP2) reproduces the shell structure of the exact correlation potential very well but consistently overestimates the corresponding shell oscillations. In the case of spin-polarized atoms one observes a strong interrelation between the correlation potentials of the two spin channels, which is completely absent for standard density functionals. However, our results also demonstrate that E(c) (MP2) can only serve as a first step towards the construction of a suitable implicit correlation functional: The fundamental variational instability of this functional is recovered for beryllium, for which a breakdown of the self-consistent Kohn-Sham iteration is observed. Moreover, even for those atoms for which the self-consistent iteration is stable, the results indicate that the inclusion of v(c) (MP2) in the total Kohn-Sham potential does not lead to an improvement compared to the complete neglect of the correlation potential.     

In search of definitive signatures of the elusive NCCO radical

Previous experimental assignments of the fundamental vibrational frequencies of NCCO have been brought into question by subsequent unsuccessful attempts to observe IR signatures of this radical at these frequencies. Here we compute the fundamental vibrational frequencies by applying second-order vibrational perturbation theory to the complete quartic force field computed at the all-electron (AE) coupled cluster singles, doubles, and perturbative triples level [CCSD(T)] with the correlation-consistent, polarized core-valence quadruple-zeta (cc-pCVQZ) basis set, which has tight functions to correctly describe core correlation. The AE-CCSD(T)/cc-pCVQZ geometric parameters are r(e)(N-C)=1.1623 A, r(e)(C-C)=1.4370 A, r(e)(C-O)=1.1758 A, theta(e)(N-C-C)=168.55 degrees , and theta(e)(C-C-O)=132.22 degrees . Our CCSD(T)/cc-pCVQZ values of the characteristic stretching frequencies nu(1) and nu(2) are 2171 and 1898 cm(-1), respectively, in stark contrast to the experimentally derived values of 2093 and 1774 cm(-1). Finally, focal-point extrapolations using correlation-consistent basis sets cc-pVXZ (X=D,T,Q,5,6) and electron correlation treatments as extensive as full coupled cluster singles, doubles, and triples (CCSDT) with perturbative accounting of quadruple excitations [CCSDT(Q)] determine the vibrationless barrier to linearity of NCCO and the dissociation energy (D(0)) of NCCO-->NC+CO to be 8.4 and 26.5 kcal mol(-1), respectively. Using our precisely determined dissociation energy, we recommend a new 0 K enthalpy of formation for NCCO of 50.9+/-0.3 kcal mol(-1).     

Density-functional theory with effective potential expressed as a direct mapping of the external potential: applications to atomization energies and ionization potentials

In this paper the authors further develop and apply the direct-mapping density functional theory to calculations of the atomization energies and ionization potentials. Single-particle orbitals are determined by solving the Kohn-Sham [Phys. Rev. A. 140, 1133 (1965)] equations with a local effective potential expressed in terms of the external potential. A two-parametric form of the effective potential for molecules is proposed and equations for optimization of the parameters are derived using the exchange-only approximation. Orbital-dependent correlation functional is derived from the second-order perturbation theory in its Moller-Plesset-type zeroth-order approximation based on the Kohn-Sham orbitals and orbital energies. The total atomization energies and ionization potentials computed with the second-order perturbation theory were found to be in agreement with experimental values and benchmark results obtained with ab initio wave mechanics methods.     

Efficient exact exchange approximations in density-functional theory

Two approaches to approximate the Slater potential component of local exact exchange of density-functional theory are investigated. The first approach employs density fitting of the electrostatic potential integrals over two occupied orbitals and the other approach approximates the "exact" Slater potential with the potential derived from the Becke-Roussel [Phys. Rev. A. 39, 3761 (1989)] model of the exchange hole. In both cases significant time savings can be achieved for larger systems compared to the calculation of the numerical Slater potential. It is then analyzed how well the orbitals obtained from the various total exchange potentials reproduce Hartree-Fock energies and molecular properties. A large range of atoms and small molecules has been utilized, including the three DNA bases adenine, thymine, and cytosine.     

Cross sections for electron impact excitation of the vibrationally resolved A 1Pi electronic state of carbon monoxide

The authors report new differential cross section measurements for electron impact excitation of the A (1)Pi(v(')) states of carbon monoxide. The energy range is 20-200 eV. They also reanalyze the A (1)Pi(v(')) manifold cross sections of Middleton et al. [J. Phys. B 26, 1743 (1993)] in order to provide a basis for comparison with our new vibrationally resolved differential cross sections. Excellent agreement is found between the two sets of measurements at all common energies. From 20 to 200 eV the present differential cross sections are extrapolated and integrated, and the corresponding integral excitation cross sections determined. New scaled Born integral cross sections, calculated as a part of the present study, are compared against these experimental integral cross sections, with excellent agreement being found for all the A (1)Pi(v(')=0-7)<--X (1)Sigma(g) (+)(v(")=0) transitions. In addition our scaled Born integral cross sections are found to be in excellent agreement between 300 and 1500 eV with those derived from the previous experiments of Lassettre and Skerbele [J. Chem. Phys. 54, 1597 (1971)] and of Zhong et al. [Phys. Rev. A 55, 1799 (1997)] and from near threshold to 15 eV with those derived from Zobel et al. [J. Phys. B 29, 813 (1996)] and Zetner et al. (J. Phys. B 31, 2395 (1998)].     

Investigation of excess adsorption, solvation force, and plate-fluid interfacial tension for Lennard-Jones fluid confined in slit pores

The excess Helmholtz free energy functional is formulated in terms of a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)] for a short ranged repulsion and a first-order mean-spherical approximation theory [Y. P. Tang, J. Chem. Phys. 118, 4140 (2003)] for a long ranged attraction. Within the framework of the density functional theory, the density profile, excess adsorption, solvation force, and plate-fluid interfacial tension of a Lennard-Jones fluid confined in slit pores are predicted, and the results agree well with the simulation data. The phase equilibria inside the slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal, and the plate-fluid interfacial tensions at equilibrium states are predicted consequently.     

B80 and B101-103 clusters: remarkable stability of the core-shell structures established by validated density functionals

Prompted by the very recent claim that the volleyball-shaped B(80) fullerene [X. Wang, Phys. Rev. B 82, 153409 (2010)] is lower in energy than the B(80) buckyball [N. G. Szwacki, A. Sadrzadeh, and B. I. Yakobson, Phys. Rev. Lett. 98, 166804 (2007)] and core-shell structure [J. Zhao, L. Wang, F. Li, and Z. Chen, J. Phys. Chem. A 114, 9969 (2010)], and inspired by the most recent finding of another core-shell isomer as the lowest energy B(80) isomer [S. De, A. Willand, M. Amsler, P. Pochet, L. Genovese, and S. Goedecher, Phys. Rev. Lett. 106, 225502 (2011)], we carefully evaluated the performance of the density functional methods in the energetics of boron clusters and confirmed that the core-shell construction (stuffed fullerene) is thermodynamically the most favorable structural pattern for B(80). Our global minimum search showed that both B(101) and B(103) also prefer a core-shell structure and that B(103) can reach the complete core-shell configuration. We called for great attention to the theoretical community when using density functionals to investigate boron-related nanomaterials.     

Solid/solid phase transitions in confined thin films: a zero temperature approach

We report a density functional theory study of confinement induced solid/solid phase transitions in a thin film (modeled as methane) at T=0. The solid film is confined by two graphite surfaces represented by a mean-field potential. As the wall separation is varied the structure of the confined film changes, which influences its density and the solvation force. Using the directly accessible grand canonical potential density we determine the stable phases and calculate the exact location of the phase transitions. We observe a series of phases having square and triangular symmetry. At low wall separations we find zig-zag buckling and an asymmetric buckled phase, whose structure is consistent with the strongest buckling instability of a triangular monolayer predicted by Chou and Nelson [Phys. Rev. E 48, 4611 (1993)] but, to our knowledge, has not been observed as a stable phase before. We find that the two-dimensional order parameters Psi(4) (square symmetry) and Psi(6) (triangular symmetry) show unphysical behavior in the transition region between square and triangular symmetry. Thus, in the present model they fail to predict the right location of the phase transitions.