Application of A Distributed Nucleus Approximation In Grid Based Minimization of the Kohn-Sham Energy Functional

In the distributed nucleus approximation we represent the singular nucleus as smeared over a smallportion of a Cartesian grid. Delocalizing the nucleus allows us to solve the Poisson equation for theoverall electrostatic potential using a linear scaling multigrid algorithm.This work is done in the context of minimizing the Kohn-Sham energy functionaldirectly in real space with a multiscale approach. The efficacy of the approximation is illustrated bylocating the ground state density of simple one electron atoms and moleculesand more complicated multiorbital systems.Comment: Submitted to JCP (July 1, 1995 Issue), latex, 27pages, 2figure

Timesaving Double-Grid Method for Real-Space Electronic-Structure Calculations

We present a simple and efficient technique in ab initio electronic-structure calculation utilizing real-space double-grid with a high density of grid points in the vicinity of nuclei. This technique promises to greatly reduce the overhead for performing the integrals that involves non-local parts of pseudopotentials, with keeping a high degree of accuracy. Our procedure gives rise to no Pulay forces, unlike other real-space methods using adaptive coordinates. Moreover, we demonstrate the potential power of the method by calculating several properties of atoms and molecules.Comment: 4 pages, 5 figure

Properties of short-range and long-range correlation energy density functionals from electron-electron coalescence

The combination of density functional theory with other approaches to the many-electron problem through the separation of the electron-electron interaction into a short-range and a long-range contribution is a promising method, which is raising more and more interest in recent years. In this work some properties of the corresponding correlation energy functionals are derived by studying the electron-electron coalescence condition for a modified (long-range-only) interaction. A general relation for the on-top (zero electron-electron distance) pair density is derived, and its usefulness is discussed with some examples. For the special case of the uniform electron gas, a simple parameterization of the on-top pair density for a long-range only interaction is presented and supported by calculations within the ``extended Overhauser model''. The results of this work can be used to build self-interaction corrected short-range correlation energy functionals.Comment: revised version, to appear in Phys. Rev.

Performance of the modified Becke-Johnson potential

Very recently, in the 2011 version of the Wien2K code, the long standing shortcome of the codes based on Density Functional Theory, namely, its impossibility to account for the experimental band gap value of semiconductors, was overcome. The novelty is the introduction of a new exchange and correlation potential, the modified Becke-Johnson potential (mBJLDA). In this paper, we report our detailed analysis of this recent work. We calculated using this code, the band structure of forty one semiconductors and found an important improvement in the overall agreement with experiment as Tran and Blaha [{\em Phys. Rev. Lett.} 102, 226401 (2009)] did before for a more reduced set of semiconductors. We find, nevertheless, within this enhanced set, that the deviation from the experimental gap value can reach even much more than 20%, in some cases. Furthermore, since there is no exchange and correlation energy term from which the mBJLDA potential can be deduced, a direct optimization procedure to get the lattice parameter in a consistent way is not possible as in the usual theory. These authors suggest that a LDA or a GGA optimization procedure is used previous to a band structure calculation and the resulting lattice parameter introduced into the 2011 code. This choice is important since small percentage differences in the lattice parameter can give rise to quite higher percentage deviations from experiment in the predicted band gap value.Comment: 10 pages, 2 figures, 5 Table

Exchange energy in the local Airy gas approximation

The Airy gas model of the edge electron gas is used to construct an exchange-energy functional which is an alternative to those obtained in the local density and generalized gradient approximations. Test calculations for rare gas atoms, molecules, solids and surfaces show that the Airy gas functional performs better than the local density approximation in all cases and better than the generalized gradient approximation for solids and surfaces. Typeset using REVTEX 1 Since the pioneering papers on density functional theory (DFT) [1,2] there has been a constant search for exchange-correlation functionals of chemical accuracy. This includes the works on the generalized gradient approximation (GGA) [3–7] which are dedicated efforts to construct local functionals for inhomogeneous systems ranging from atoms to solids based on the uniform electron gas, i.e., the local density approximation (LDA), and density gradient corrections, as well as the development of a number gradien

Rearrangement of cluster structure during fission processes

Results of molecular dynamics simulations of fission reactions $Na_{10}^{2+} \to Na_7^+ + Na_3^+$ and $Na_{18}^{2+} \to 2 Na_9^+$ are presented. Dependence of the fission barriers on isomer structure of the parent cluster is analyzed. It is demonstrated that the energy necessary for removing homothetic groups of atoms from the parent cluster is largely independent of the isomer form of the parent cluster. Importance of rearrangement of the cluster structure during the fission process is elucidated. This rearrangement may include transition to another isomer state of the parent cluster before actual separation of the daughter fragments begins and/or forming a "neck" between the separating fragments

Quasi-chemical study of Be2+^{2+}(aq) speciation

Be$^{2+}$(aq) hydrolysis can to lead to the formation of multi-beryllium clusters, but the thermodynamics of this process has not been resolved theoretically. We study the hydration state of an isolated Be$^{2+}$ ion using both the quasi-chemical theory of solutions and ab initio molecular dynamics. These studies confirm that Be$^{2+}$(aq) is tetra-hydrated. The quasi-chemical approach is then applied to then the deprotonation of Be(H_2O)_4^{2+}} to give BeOH(H_2O)_3{}^{+}}. The calculated pK$_a$ of 3.8 is in good agreement with the experimentally suggested value around 3.5. The calculated energetics for the formation of BeOHBe$^{3+}$ are then obtained in fair agreement with experiments.Comment: 11 pages, 3 figure

Diagnostics for specific PAHs in the far-IR: searching neutral naphthalene and anthracene in the Red Rectangle

Context. In the framework of the interstellar polycyclic aromatic hydrocarbons (PAHs) hypothesis, far-IR skeletal bands are expected to be a fingerprint of single species in this class. Aims. We address the question of detectability of low energy PAH vibrational bands, with respect to spectral contrast and intensity ratio with ``classical'' Aromatic Infrared Bands (AIBs). Methods. We extend our extablished Monte-Carlo model of the photophysics of specific PAHs in astronomical environments, to include rotational and anharmonic band structure. The required molecular parameters were calculated in the framework of the Density Functional Theory. Results. We calculate the detailed spectral profiles of three low-energy vibrational bands of neutral naphthalene, and four low-energy vibrational bands of neutral anthracene. They are used to establish detectability constraints based on intensity ratios with ``classical'' AIBs. A general procedure is suggested to select promising diagnostics, and tested on available Infrared Space Observatory data for the Red Rectangle nebula. Conclusions. The search for single, specific PAHs in the far-IR is a challenging, but promising task, especially in view of the forthcoming launch of the Herschel Space Observatory.Comment: 13 pages, 13 figures, accepted for publication in A&

Analysis of OPM potentials for multiplet states of 3d transition metal atoms

We apply the optimized effective potential method (OPM) to the multiplet energies of the 3d$^n$ transition metal atoms, where the orbital dependence of the energy functional with respect to orbital wave function is the single-configuration HF form. We find that the calculated OPM exchange potential can be represented by the following two forms. Firstly, the difference between OPM exchange potentials of the multiplet states can be approximated by the linear combination of the potentials derived from the Slater integrals $F^2({\rm 3d,3d})$ and $F^4({\rm 3d,3d})$ for the average energy of the configuration. Secondly, the OPM exchange potential can be expressed as the linear combination of the OPM exchange potentials of the single determinants.Comment: 15 pages, 6 figures, to be published in J. Phys.