Differentiability of Lieb functional in electronic density functional theory
NASA Astrophysics Data System (ADS)
Lammert, Paul E.
A solid understanding of the Lieb functional FL is important because of its centrality in the foundations of electronic density functional theory. A basic question is whether directional derivatives of FL at an ensemble-V-representable density are given by (minus) the potential. A widely accepted purported proof that FL is Gâteaux differentiable at EV-representable densities would say, ?yes.? But that proof is fallacious, as shown here. FL is not Gâteaux differentiable in the normal sense, nor is it continuous. By means of a constructive approach, however, we are able to show that the derivative of FL at an EV-representable density ?0 in the direction of ?1 is given by the potential if ?0 and ?1 are everywhere strictly greater than zero, and they and the ground state wave function have square integrable derivatives through second order.
Gutzwiller density functional theory for correlated electron systems
Ho, K. M.; Schmalian, J.; Wang, C. Z.
2008-02-04
We develop a density functional theory (DFT) and formalism for correlated electron systems by taking as reference an interacting electron system that has a ground state wave function which exactly obeys the Gutzwiller approximation for all one-particle operators. The solution of the many-electron problem is mapped onto the self-consistent solution of a set of single-particle Schroedinger equations, analogously to standard DFT-local density approximation calculations.
Optimal-transport formulation of electronic density-functional theory
NASA Astrophysics Data System (ADS)
Buttazzo, Giuseppe; De Pascale, Luigi; Gori-Giorgi, Paola
2012-06-01
The most challenging scenario for Kohn-Sham density-functional theory, that is, when the electrons move relatively slowly trying to avoid each other as much as possible because of their repulsion (strong-interaction limit), is reformulated here as an optimal transport (or mass transportation theory) problem, a well-established field of mathematics and economics. In practice, we show that to solve the problem of finding the minimum possible internal repulsion energy for N electrons in a given density ρ(r) is equivalent to find the optimal way of transporting N-1 times the density ρ into itself, with the cost function given by the Coulomb repulsion. We use this link to set the strong-interaction limit of density-functional theory on firm ground and to discuss the potential practical aspects of this reformulation.
The exact density functional for two electrons in one dimension
NASA Astrophysics Data System (ADS)
Cohen, Aron; Mori-Sanchez, Paula
The exact universal density functional F [ ρ ] is calculated for real space two-electron densities in one dimension ρ (x) with a soft-Coulomb interaction. It is calculated by the Levy constrained search F [ ρ ] =minΨ-->ρ < Ψ | \\Tcirc +\\Vcircee | Ψ > over wavefunctions of a two-dimensional Hilbert space Ψ (x1 ,x2) --> ρ (x1) and can be directly visualized. We do an approximate constrained search via density matrices and a direct approximation to natural orbitals. This allows us to make an accurate approximation to the exact functional that is calculated using a search over potentials. We investigate the exact functional and the performance of many approximations on some of the most challenging electronic structure in two-electron systems, from strongly-correlated electron transfer to the description of a localized-delocalized transition. The exact Kohn-Sham potential, vs (x) , and exact Kohn-Sham eigenvalues, ɛi, are calculated and this allows us to discuss the band-gap problem versus the perspective of the exact density functional F [ ρ ] for all numbers of electrons. We calculate the derivative discontinuity of the exact functional in an example of a Mott-Insulator, one-dimensional stretched H2.
Excitations and benchmark ensemble density functional theory for two electrons
Pribram-Jones, Aurora; Burke, Kieron; Yang, Zeng-hui; Ullrich, Carsten A.; Trail, John R.; Needs, Richard J.
2014-05-14
A new method for extracting ensemble Kohn-Sham potentials from accurate excited state densities is applied to a variety of two-electron systems, exploring the behavior of exact ensemble density functional theory. The issue of separating the Hartree energy and the choice of degenerate eigenstates is explored. A new approximation, spin eigenstate Hartree-exchange, is derived. Exact conditions that are proven include the signs of the correlation energy components and the asymptotic behavior of the potential for small weights of the excited states. Many energy components are given as a function of the weights for two electrons in a one-dimensional flat box, in a box with a large barrier to create charge transfer excitations, in a three-dimensional harmonic well (Hooke's atom), and for the He atom singlet-triplet ensemble, singlet-triplet-singlet ensemble, and triplet bi-ensemble.
Excess electrons in ice: a density functional theory study.
Bhattacharya, Somesh Kr; Inam, Fakharul; Scandolo, Sandro
2014-02-21
We present a density functional theory study of the localization of excess electrons in the bulk and on the surface of crystalline and amorphous water ice. We analyze the initial stages of electron solvation in crystalline and amorphous ice. In the case of crystalline ice we find that excess electrons favor surface states over bulk states, even when the latter are localized at defect sites. In contrast, in amorphous ice excess electrons find it equally favorable to localize in bulk and in surface states which we attribute to the preexisting precursor states in the disordered structure. In all cases excess electrons are found to occupy the vacuum regions of the molecular network. The electron localization in the bulk of amorphous ice is assisted by its distorted hydrogen bonding network as opposed to the crystalline phase. Although qualitative, our results provide a simple interpretation of the large differences observed in the dynamics and localization of excess electrons in crystalline and amorphous ice films on metals. PMID:24401958
Momentum distribution function of the electron gas at metallic densities
NASA Astrophysics Data System (ADS)
Takada, Yasutami; Yasuhara, H.
1991-10-01
The momentum distribution function n(k) of the electron gas is calculated in the effective-potential-expansion method at metallic densities. The recently established self-consistency relation between n(k) and the correlation energy [Y. Takada and T. Kita, J. Phys. Soc. Jpn. 60, 25 (1991)] is employed to check the accuracy of our results. This check shows that the effective-potential-expansion method provides probably the exact and at least more accurate results of n(k) than all the other methods that have given n(k) thus far.
Ionospheric mapping functions based on electron density fields
NASA Astrophysics Data System (ADS)
Zus, Florian; Deng, Zhiguo; Heise, Stefan; Wickert, Jens
2016-04-01
We developed an ionospheric Mapping Function (MF) for the Global Navigation Satellite System (GNSS) which is based on the electron density field of the International Reference Ionosphere (IRI). The station specific MF utilizes a look-up table which contains a set of ray-traced ionospheric delays. Hence, unlike the simple MFs that are currently in use, the developed MF depends on the time, location, elevation and azimuth angle. Ray-bending is taken into account, which implies that the MF depends on the carrier frequency as well. This frequency dependency of the MF can be readily used to examine higher-order ionospheric effects due to ray-bending. We compare the proposed MF with the so-called single layer model MF and find significant differences in particular around the equatorial anomaly. In so-far as the proposed MF is based on a realistic electron density field (IRI) our comparison shows the potential error of the single-layer model MF in practice. We conclude that the developed MF concept might be valuable in the GNSS Total Electron Content estimation. The frequency dependency of the MF can be used to mitigate higher-order ionospheric effects.
Electronic properties of graphene nanoribbons: A density functional investigation
Kumar, Sandeep Sharma, Hitesh
2015-05-15
Density functional theory calculations have been performed on graphene nano ribbons (GNRs) to investigate the electronic properties as a function of chirality, size and hydrogenation on the edges. The calculations were performed on GNRs with armchair and zigzag configurations with 28, 34, 36, 40, 50, 56, 62, 66 carbon atoms. The structural stability of AGNR and ZGNR increases with the size of nanoribbon where as hydrogenation of GNR tends to lowers their structural stability. All GNRs considered have shown semiconducting behavior with HOMO-LUMO gap decreasing with the increase in the GNR size. The hydrogenation of GNR decreases its HOMO-LUMO gap significantly. The results are in agreement with the available experimental and theoretical results.
Equation satisfied by the energy-density functional for electron-electron mutual Coulomb repulsion
Joubert, Daniel P.
2011-10-15
It is shown that the electron-electron mutual Coulomb repulsion energy-density functional V{sub ee}{sup {gamma}}[{rho}] satisfies the equationV{sub ee}{sup {gamma}}[{rho}{sub N}{sup 1}]-V{sub ee}{sup {gamma}}[{rho}{sub N-1}{sup {gamma}}]={integral}d{sup 3}r({delta}V{sub ee}{sup {gamma}}[{rho}{sub N}{sup 1}]/{delta}{rho}{sub N}{sup 1}(r))[{rho}{sub N}{sup 1}(r)-{rho}{sub N-1}{sup {gamma}}(r)], where {rho}{sub N}{sup 1}(r) and {rho}{sub N-1}{sup {gamma}}(r) are N-electron and (N-1)-electron densities determined from the same adiabatic scaled external potential of the N-electron system at coupling strength {gamma}.
NASA Astrophysics Data System (ADS)
Putaja, A.; Eich, F. G.; Baldsiefen, T.; Räsänen, E.
2016-03-01
Physically valid and numerically efficient approximations for the exchange and correlation energy are critical for reduced-density-matrix-functional theory to become a widely used method in electronic structure calculations. Here we examine the physical limits of power functionals of the form f (n ,n') =(nn')α for the scaling function in the exchange-correlation energy. To this end we obtain numerically the minimizing momentum distributions for the three- and two-dimensional homogeneous electron gas, respectively. In particular, we examine the limiting values for the power α to yield physically sound solutions that satisfy the Lieb-Oxford lower bound for the exchange-correlation energy and exclude pinned states with the condition n (k )<1 for all wave vectors k . The results refine the constraints previously obtained from trial momentum distributions. We also compute the values for α that yield the exact correlation energy and its kinetic part for both the three- and two-dimensional electron gas. In both systems, narrow regimes of validity and accuracy are found at α ≳0.6 and at rs≳10 for the density parameter, corresponding to relatively low densities.
NASA Astrophysics Data System (ADS)
Tachikawa, Hiroto; Iyama, Tetsuji; Kawabata, Hiroshi
2016-05-01
Electronic structures and formation mechanism of hydrogen functionalized carbon nanotube (CNT) have been investigated by means of density functional theory (DFT) method. The mechanism of hydrogen addition reaction to the CNT surface was also investigated. Pure and boron-nitrogen (BN) substituted CNT (denoted by CNT and BN-CNT, respectively) were examined as the carbon nanotubes. It was found that the additions of hydrogen atom to B (boron atom) and C (carbon atom) sites of BN-CNT proceed without activation barrier, whereas the hydrogenation of N (nitrogen atom) site needs the activation energy. The electronic states of hydrogen functionalized CNT and BN-CNT were discussed on the basis of theoretical results.
De Proft, F; Van Alsenoy, C; Peeters, A; Langenaeker, W; Geerlings, P
2002-09-01
In the Hirshfeld partitioning of the electron density, the molecular electron density is decomposed in atomic contributions, proportional to the weight of the isolated atom density in the promolecule density, constructed by superimposing the isolated atom electron densities placed on the positions the atoms have in the molecule. A maximal conservation of the information of the isolated atoms in the atoms-in-molecules is thereby secured. Atomic charges, atomic dipole moments, and Fukui functions resulting from the Hirshfeld partitioning of the electron density are computed for a large series of molecules. In a representative set of organic and hypervalent molecules, they are compared with other commonly used population analysis methods. The expected bond polarities are recovered, but the charges are much smaller compared to other methods. Condensed Fukui functions for a large number of molecules, undergoing an electrophilic or a nucleophilic attack, are computed and compared with the HOMO and LUMO densities, integrated over the Hirshfeld atoms in molecules. PMID:12116389
Electronic states of alkyl-radical-functionalized C20 fullerene using density functional theory
NASA Astrophysics Data System (ADS)
Abe, Shigeaki; Kawano, Shimpei; Toida, Yu; Nakamura, Mariko; Inoue, Satoshi; Sano, Hidehiko; Yoshida, Yasuhiro; Kawabata, Hiroshi; Tachikawa, Hiroto
2016-03-01
The structures and electronic states of alkyl-radical-functionalized C20 fullerenes (denoted by C20-R) have been investigated using density functional theory (DFT). The different alkyl radicals investigated were methyl, ethyl, propyl, and butyl radicals. The DFT calculation indicated that the alkyl radical binds to the carbon atom of C20 in the on-top site, thus forming a strong C-C single bond. The binding energies of the alkyl radicals to C20 were calculated to be 83.9-86.6 kcal/mol at the CAM-B3LYP/6-311G(d,p) level. The electronic states of the C20-R complex are discussed on the basis of the theoretical results.
Electronic states of aryl radical functionalized graphenes: Density functional theory study
NASA Astrophysics Data System (ADS)
Tachikawa, Hiroto; Kawabata, Hiroshi
2016-06-01
Functionalized graphenes are known as a high-performance molecular device. In the present study, the structures and electronic states of the aryl radical functionalized graphene have been investigated by the density functional theory (DFT) method to elucidate the effects of functionalization on the electronic states of graphene (GR). Also, the mechanism of aryl radical reaction with GR was investigated. The benzene, biphenyl, p-terphenyl, and p-quaterphenyl radicals [denoted by (Bz) n (n = 1–4), where n means numbers of benzene rings in aryl radical] were examined as aryl radicals. The DFT calculation of GR–(Bz) n (n = 1–4) showed that the aryl radical binds to the carbon atom of GR, and a C–C single bond was formed. The binding energies of aryl radicals to GR were calculated to be ca. 6.0 kcal mol‑1 at the CAM-B3LYP/6-311G(d,p) level. It was found that the activation barrier exists in the aryl radical addition: the barrier heights were calculated to be 10.0 kcal mol‑1. The electronic states of GR–(Bz) n were examined on the basis of theoretical results.
Many-electron expansion: A density functional hierarchy for strongly correlated systems
NASA Astrophysics Data System (ADS)
Zhu, Tianyu; de Silva, Piotr; van Aggelen, Helen; Van Voorhis, Troy
2016-05-01
Density functional theory (DFT) is the de facto method for the electronic structure of weakly correlated systems. But for strongly correlated materials, common density functional approximations break down. Here, we derive a many-electron expansion (MEE) in DFT that accounts for successive one-, two-, three-, ... particle interactions within the system. To compute the correction terms, the density is first decomposed into a sum of localized, nodeless one-electron densities (ρi). These one-electron densities are used to construct relevant two- (ρi+ρj ), three- (ρi+ρj+ρk ), ... electron densities. Numerically exact results for these few-particle densities can then be used to correct an approximate density functional via any of several many-body expansions. We show that the resulting hierarchy gives accurate results for several important model systems: the Hubbard and Peierls-Hubbard models in 1D and the pure Hubbard model in 2D. We further show that the method is numerically convergent for strongly correlated systems: applying successively higher order corrections leads to systematic improvement of the results. MEE thus provides a hierarchy of density functional approximations that applies to both weakly and strongly correlated systems.
Reimers, Jeffrey R; Solomon, Gemma C; Gagliardi, Alessio; Bilić, Ante; Hush, Noel S; Frauenheim, Thomas; Di Carlo, Aldo; Pecchia, Alessandro
2007-07-01
A review is presented of the nonequilibrium Green's function (NEGF) method "gDFTB" for evaluating elastic and inelastic conduction through single molecules employing the density functional tight-binding (DFTB) electronic structure method. This focuses on the possible advantages that DFTB implementations of NEGF have over conventional methods based on density functional theory, including not only the ability to treat large irregular metal-molecule junctions with high nonequilibrium thermal distributions but perhaps also the ability to treat dispersive forces, bond breakage, and open-shell systems and to avoid large band lineup errors. New results are presented indicating that DFTB provides a useful depiction of simple gold-thiol interactions. Symmetry is implemented in DFTB, and the advantages it brings in terms of large savings of computational resources with significant increase in numerical stability are described. The power of DFTB is then harnessed to allow the use of gDFTB as a real-time tool to discover the nature of the forces that control inelastic charge transport through molecules and the role of molecular symmetry in determining both elastic and inelastic transport. Future directions for the development of the method are discussed. PMID:17530826
Putz, Mihai V.
2009-01-01
The density matrix theory, the ancestor of density functional theory, provides the immediate framework for Path Integral (PI) development, allowing the canonical density be extended for the many-electronic systems through the density functional closure relationship. Yet, the use of path integral formalism for electronic density prescription presents several advantages: assures the inner quantum mechanical description of the system by parameterized paths; averages the quantum fluctuations; behaves as the propagator for time-space evolution of quantum information; resembles Schrödinger equation; allows quantum statistical description of the system through partition function computing. In this framework, four levels of path integral formalism were presented: the Feynman quantum mechanical, the semiclassical, the Feynman-Kleinert effective classical, and the Fokker-Planck non-equilibrium ones. In each case the density matrix or/and the canonical density were rigorously defined and presented. The practical specializations for quantum free and harmonic motions, for statistical high and low temperature limits, the smearing justification for the Bohr’s quantum stability postulate with the paradigmatic Hydrogen atomic excursion, along the quantum chemical calculation of semiclassical electronegativity and hardness, of chemical action and Mulliken electronegativity, as well as by the Markovian generalizations of Becke-Edgecombe electronic focalization functions – all advocate for the reliability of assuming PI formalism of quantum mechanics as a versatile one, suited for analytically and/or computationally modeling of a variety of fundamental physical and chemical reactivity concepts characterizing the (density driving) many-electronic systems. PMID:20087467
NASA Astrophysics Data System (ADS)
Sun, Jianwei; Perdew, John P.; Yang, Zenghui; Peng, Haowei
2016-05-01
The uniform electron gas and the hydrogen atom play fundamental roles in condensed matter physics and quantum chemistry. The former has an infinite number of electrons uniformly distributed over the neutralizing positively charged background, and the latter only one electron bound to the proton. The uniform electron gas was used to derive the local spin density approximation to the exchange-correlation functional that undergirds the development of the Kohn-Sham density functional theory. We show here that the ground-state exchange-correlation energies of the hydrogen atom and many other 1- and 2-electron systems are modeled surprisingly well by a different local spin density approximation (LSDA0). LSDA0 is constructed to satisfy exact constraints but agrees surprisingly well with the exact results for a uniform two-electron density in a finite, curved three-dimensional space. We also apply LSDA0 to excited or noded 1-electron densities, where it works less well. Furthermore, we show that the localization of the exact exchange hole for a 1- or 2-electron ground state can be measured by the ratio of the exact exchange energy to its optimal lower bound.
Sun, Jianwei; Perdew, John P; Yang, Zenghui; Peng, Haowei
2016-05-21
The uniform electron gas and the hydrogen atom play fundamental roles in condensed matter physics and quantum chemistry. The former has an infinite number of electrons uniformly distributed over the neutralizing positively charged background, and the latter only one electron bound to the proton. The uniform electron gas was used to derive the local spin density approximation to the exchange-correlation functional that undergirds the development of the Kohn-Sham density functional theory. We show here that the ground-state exchange-correlation energies of the hydrogen atom and many other 1- and 2-electron systems are modeled surprisingly well by a different local spin density approximation (LSDA0). LSDA0 is constructed to satisfy exact constraints but agrees surprisingly well with the exact results for a uniform two-electron density in a finite, curved three-dimensional space. We also apply LSDA0 to excited or noded 1-electron densities, where it works less well. Furthermore, we show that the localization of the exact exchange hole for a 1- or 2-electron ground state can be measured by the ratio of the exact exchange energy to its optimal lower bound. PMID:27208927
Guido, Ciro A. Cortona, Pietro; Adamo, Carlo; Institut Universitaire de France, 103 Bd Saint-Michel, F-75005 Paris
2014-03-14
We extend our previous definition of the metric Δr for electronic excitations in the framework of the time-dependent density functional theory [C. A. Guido, P. Cortona, B. Mennucci, and C. Adamo, J. Chem. Theory Comput. 9, 3118 (2013)], by including a measure of the difference of electronic position variances in passing from occupied to virtual orbitals. This new definition, called Γ, permits applications in those situations where the Δr-index is not helpful: transitions in centrosymmetric systems and Rydberg excitations. The Γ-metric is then extended by using the Natural Transition Orbitals, thus providing an intuitive picture of how locally the electron density changes during the electronic transitions. Furthermore, the Γ values give insight about the functional performances in reproducing different type of transitions, and allow one to define a “confidence radius” for GGA and hybrid functionals.
Nishimoto, Yoshio
2015-09-01
We develop a formalism for the calculation of excitation energies and excited state gradients for the self-consistent-charge density-functional tight-binding method with the third-order contributions of a Taylor series of the density functional theory energy with respect to the fluctuation of electron density (time-dependent density-functional tight-binding (TD-DFTB3)). The formulation of the excitation energy is based on the existing time-dependent density functional theory and the older TD-DFTB2 formulae. The analytical gradient is computed by solving Z-vector equations, and it requires one to calculate the third-order derivative of the total energy with respect to density matrix elements due to the inclusion of the third-order contributions. The comparison of adiabatic excitation energies for selected small and medium-size molecules using the TD-DFTB2 and TD-DFTB3 methods shows that the inclusion of the third-order contributions does not affect excitation energies significantly. A different set of parameters, which are optimized for DFTB3, slightly improves the prediction of adiabatic excitation energies statistically. The application of TD-DFTB for the prediction of absorption and fluorescence energies of cresyl violet demonstrates that TD-DFTB3 reproduced the experimental fluorescence energy quite well. PMID:26342360
Nishimoto, Yoshio
2015-09-07
We develop a formalism for the calculation of excitation energies and excited state gradients for the self-consistent-charge density-functional tight-binding method with the third-order contributions of a Taylor series of the density functional theory energy with respect to the fluctuation of electron density (time-dependent density-functional tight-binding (TD-DFTB3)). The formulation of the excitation energy is based on the existing time-dependent density functional theory and the older TD-DFTB2 formulae. The analytical gradient is computed by solving Z-vector equations, and it requires one to calculate the third-order derivative of the total energy with respect to density matrix elements due to the inclusion of the third-order contributions. The comparison of adiabatic excitation energies for selected small and medium-size molecules using the TD-DFTB2 and TD-DFTB3 methods shows that the inclusion of the third-order contributions does not affect excitation energies significantly. A different set of parameters, which are optimized for DFTB3, slightly improves the prediction of adiabatic excitation energies statistically. The application of TD-DFTB for the prediction of absorption and fluorescence energies of cresyl violet demonstrates that TD-DFTB3 reproduced the experimental fluorescence energy quite well.
Density functional theory for d- and f-electron materials and compounds
Mattson, Ann E.; Wills, John M.
2016-02-12
Here, the fundamental requirements for a computationally tractable Density Functional Theory-based method for relativistic f- and (nonrelativistic) d-electron materials and compounds are presented. The need for basing the Kohn–Sham equations on the Dirac equation is discussed. The full Dirac scheme needs exchange-correlation functionals in terms of four-currents, but ordinary functionals, using charge density and spin-magnetization, can be used in an approximate Dirac treatment. The construction of a functional that includes the additional confinement physics needed for these materials is illustrated using the subsystem-functional scheme. If future studies show that a full Dirac, four-current based, exchange-correlation functional is needed, the subsystemmore » functional scheme is one of the few schemes that can still be used for constructing functional approximations.« less
Long-range density-matrix-functional theory: Application to a modified homogeneous electron gas
Pernal, Katarzyna
2010-05-15
We propose a method that employs functionals of the one-electron reduced density matrix (density matrix) to capture long-range effects of electron correlation. The complementary short-range regime is treated with density functionals. In an effort to find approximations for the long-range density-matrix functional, a modified power functional is applied to the homogeneous electron gas with Coulomb interactions replaced by their corresponding long-range counterparts. For the power {beta}=1/2 and the range-separation parameter {omega}=1/r{sub s}, the functional reproduces the correlation and the kinetic correlation energies with a remarkable accuracy for intermediate and large values of r{sub s}. Analysis of the Euler equation corresponding to this functional reveals correct r{sub s} expansion of the correlation energy in the limit of large r{sub s}. The first expansion coefficient is in very good agreement with that obtained from the modified Wigner-Seitz model.
Lehtovaara, Lauri; Havu, Ville; Puska, Martti
2011-10-21
We present an all-electron method for time-dependent density functional theory which employs hierarchical nonuniform finite-element bases and the time-propagation approach. The method is capable of treating linear and nonlinear response of valence and core electrons to an external field. We also introduce (i) a preconditioner for the propagation equation, (ii) a stable way to implement absorbing boundary conditions, and (iii) a new kind of absorbing boundary condition inspired by perfectly matched layers. PMID:22029294
Density functional theory for low-energy electron-molecule scattering
NASA Astrophysics Data System (ADS)
Burke, Kieron; Wasserman, Adam
2004-09-01
Time-dependent density functional theory (TDDFT) is becoming popular as an approach to time-dependent electronic problems[1]. In the weak field regime, TDDFT predicts electronic transition frequencies and optical spectra of atoms, molecules, clusters, and solids, with an accuracy comparable to high-level wavefunction calculations at a fraction of the computational cost[2]. For large systems, TDDFT is the method of choice. Given the importance of correlation effects in low-energy electron-molecule scattering, extracting scattering amplitudes from TDDFT appears desirable. I will review this background, and outline how this can be done[3]. Detailed results will be shown by Wasserman in another talk. [1] Time-Dependent Density Functional Theory, M.A.L. Marques and E.K.U. Gross, Annu. Rev. Phys. Chem. 55, 427 (2004). [2] Time-dependent density functional theory in quantum chemistry, F. Furche and K. Burke, to appear in 1st vol. of Annu. Rev. of Computational Chemistry (2004) [3] Electron-molecule scattering from time-dependent density functional theory A. Wasserman, N.T. Maitra, and K. Burke, submitted (see http:dft.rutgers.edu/pubs/publist.html).
Employing homogeneity properties of density functionals to determine the total electronic energy
NASA Astrophysics Data System (ADS)
Morrison, Robert C.; Parr, Robert G.
1996-05-01
The exchange-correlation energy is calculated from expressions for (Exc+Tc)/2 and (Exc-Tc)/2. The expression for Exc+Tc is known exactly. The expression for Exc-Tc is not known exactly and we approximate Gxc-Tc, where Gxc=Exc+J/N, by assuming that it is a functional that is homogeneous of degree 1 in the electron density and that Tc is homogeneous of degree 0 in the electron density. The resulting formula for the total energy reduces the average of the errors in the Parr and Ghosh
NASA Astrophysics Data System (ADS)
Farzanehpour, Mehdi; Tokatly, Ilya; Nano-Bio Spectroscopy Group; ETSF Scientific Development Centre Team
2015-03-01
We present a rigorous formulation of the time-dependent density functional theory for interacting lattice electrons strongly coupled to cavity photons. We start with an example of one particle on a Hubbard dimer coupled to a single photonic mode, which is equivalent to the single mode spin-boson model or the quantum Rabi model. For this system we prove that the electron-photon wave function is a unique functional of the electronic density and the expectation value of the photonic coordinate, provided the initial state and the density satisfy a set of well defined conditions. Then we generalize the formalism to many interacting electrons on a lattice coupled to multiple photonic modes and prove the general mapping theorem. We also show that for a system evolving from the ground state of a lattice Hamiltonian any density with a continuous second time derivative is locally v-representable. Spanish Ministry of Economy and Competitiveness (Grant No. FIS2013-46159-C3-1-P), Grupos Consolidados UPV/EHU del Gobierno Vasco (Grant No. IT578-13), COST Actions CM1204 (XLIC) and MP1306 (EUSpec).
Time-Dependent Electronic Populations in Fragment-Based Time-Dependent Density Functional Theory.
Mosquera, Martín A; Wasserman, Adam
2015-08-11
Conceiving a molecule as being composed of smaller molecular fragments, or subunits, is one of the pillars of the chemical and physical sciences and leads to productive methods in quantum chemistry. Using a fragmentation scheme, efficient algorithms can be proposed to address problems in the description of chemical bond formation and breaking. We present a formally exact time-dependent density functional theory for the electronic dynamics of molecular fragments with a variable number of electrons. This new formalism is an extension of previous work [Phys. Rev. Lett. 111, 023001 (2013)]. We also introduce a stable density-inversion method that is applicable to time-dependent and ground-state density functional theories and their extensions, including those discussed in this work. PMID:26574438
Sun, Shih-Jye; Lin, Ken-Huang; Li, Jia-Yun; Ju, Shin-Pon
2014-10-07
The simulated annealing basin-hopping method incorporating the penalty function was used to predict the lowest-energy structures for ultrathin tungsten nanowires and nanotubes of different sizes. These predicted structures indicate that tungsten one-dimensional structures at this small scale do not possess B.C.C. configuration as in bulk tungsten material. In order to analyze the relationship between multi-shell geometries and electronic transfer, the electronic and structural properties of tungsten wires and tubes including partial density of state and band structures which were determined and analyzed by quantum chemistry calculations. In addition, in order to understand the application feasibility of these nanowires and tubes on nano-devices such as field emitters or chemical catalysts, the electronic stability of these ultrathin tungsten nanowires was also investigated by density functional theory calculations.
Tait, E W; Ratcliff, L E; Payne, M C; Haynes, P D; Hine, N D M
2016-05-18
Experimental techniques for electron energy loss spectroscopy (EELS) combine high energy resolution with high spatial resolution. They are therefore powerful tools for investigating the local electronic structure of complex systems such as nanostructures, interfaces and even individual defects. Interpretation of experimental electron energy loss spectra is often challenging and can require theoretical modelling of candidate structures, which themselves may be large and complex, beyond the capabilities of traditional cubic-scaling density functional theory. In this work, we present functionality to compute electron energy loss spectra within the onetep linear-scaling density functional theory code. We first demonstrate that simulated spectra agree with those computed using conventional plane wave pseudopotential methods to a high degree of precision. The ability of onetep to tackle large problems is then exploited to investigate convergence of spectra with respect to supercell size. Finally, we apply the novel functionality to a study of the electron energy loss spectra of defects on the (1 0 1) surface of an anatase slab and determine concentrations of defects which might be experimentally detectable. PMID:27094207
NASA Astrophysics Data System (ADS)
Tait, E. W.; Ratcliff, L. E.; Payne, M. C.; Haynes, P. D.; Hine, N. D. M.
2016-05-01
Experimental techniques for electron energy loss spectroscopy (EELS) combine high energy resolution with high spatial resolution. They are therefore powerful tools for investigating the local electronic structure of complex systems such as nanostructures, interfaces and even individual defects. Interpretation of experimental electron energy loss spectra is often challenging and can require theoretical modelling of candidate structures, which themselves may be large and complex, beyond the capabilities of traditional cubic-scaling density functional theory. In this work, we present functionality to compute electron energy loss spectra within the onetep linear-scaling density functional theory code. We first demonstrate that simulated spectra agree with those computed using conventional plane wave pseudopotential methods to a high degree of precision. The ability of onetep to tackle large problems is then exploited to investigate convergence of spectra with respect to supercell size. Finally, we apply the novel functionality to a study of the electron energy loss spectra of defects on the (1 0 1) surface of an anatase slab and determine concentrations of defects which might be experimentally detectable.
Density functional theory description of electronic properties of wurtzite zinc oxide
NASA Astrophysics Data System (ADS)
Franklin, L.; Ekuma, C. E.; Zhao, G. L.; Bagayoko, D.
2013-05-01
We report calculated, electronic properties of wurtzite zinc oxide (w-ZnO). We solved self-consistently the two inherently coupled equations of density functional theory (DFT), following the Bagayoko, Zhao, and Williams (BZW) method as enhanced by the work of Ekuma and Franklin (BZW-EF). We employed a local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO). Most of the calculated, electronic properties of w-ZnO are in excellent agreement with experiment, including our zero temperature band gap of 3.39 eV and the electron effective mass. The doubly self-consistent approach utilized in this work points to the ability of theory to predict accurately key properties of semiconductors and hence to inform and to guide the design and fabrication of semiconductor-based devices.
Density functional calculation of the structural and electronic properties of germanium quantum dots
Anas, M. M.; Gopir, G.
2015-04-24
We apply first principles density functional computational methods to study the structures, densities of states (DOS), and higher occupied molecular orbital (HOMO) – lowest unoccupied molecular orbital (LUMO) gaps of selected free-standing Ge semiconductor quantum dots up to 1.8nm. Our calculations are performed using numerical atomic orbital approach where linear combination of atomic orbital was applied. The surfaces of the quantum dots was passivized by hydrogen atoms. We find that surface passivation does affect the electronic properties associated with the changes of surface state, electron localization, and the energy gaps of germanium nanocrystals as well as the confinement of electrons inside the quantum dots (QDs). Our study shows that the energy gaps of germanium quantum dots decreases with the increasing dot diameter. The size-dependent variations of the computed HOMO-LUMO gaps in our quantum dots model were found to be consistent with the effects of quantum confinement reported in others theoretical and experimental calculation.
Density Functional Theory of Structural and Electronic Properties of III-N Semiconductors
Guerel, H. Hakan; Akinci, Oezden; Uenlue, Hilmi
2010-11-01
In this wok, we present the density functional theory (DFT) calculations of cubic III-N based semiconductors by using the full potential linear augmented plane-wave method plus local orbitals as implemented in the WIEN2k code. Our aim is to predict the pressure effect on structural and electronic properties of III-N binaries and ternaries. Results are given for structural properties (e.g., lattice constant, elastic constants, bulk modulus, and its pressure derivative) and electronic properties (e.g., band structure, density of states, band gaps and band widths) of GaAs, GaN, AlN, and InN binaries and GaAsN ternaries. The proposed model uses GGA exchange-correlation potential to determine band gaps of semiconductors at {Gamma}, L and X high symmetry points of Brillouin zone. The results are found in good agreement with available experimental data for structural and electronic properties of these semiconductors.
Angyan, Janos G.; Gerber, Iann C.; Savin, Andreas; Toulouse, Julien
2005-07-15
Long-range exchange and correlation effects, responsible for the failure of currently used approximate density functionals in describing van der Waals forces, are taken into account explicitly after a separation of the electron-electron interaction in the Hamiltonian into short- and long-range components. We propose a 'range-separated hybrid' functional based on a local density approximation for the short-range exchange-correlation energy, combined with a long-range exact exchange energy. Long-range correlation effects are added by a second-order perturbational treatment. The resulting scheme is general and is particularly well adapted to describe van der Waals complexes, such as rare gas dimers.
Density functional study on electronic properties of P-doped spinel silicon carbon nitride
NASA Astrophysics Data System (ADS)
Zhang, Yufen; Zhao, Xian; Cheng, Xiufeng; Mu, Yuguang
2008-08-01
We performed density functional calculations on the electronic properties of P-doped spinel silicon carbon nitride. When Si is replaced by C at the tetrahedral sites of P-doped c-Si 3N 4, the band gap can be adjusted, and an insulator-to-metal transition is predicted to occur at the C-to-Si ratio of 0.27. Finally, some possible examinations and potential applications for the large band-gap reduction are discussed.
NASA Astrophysics Data System (ADS)
Ruggenthaler, Michael; Flick, Johannes; Pellegrini, Camilla; Appel, Heiko; Tokatly, Ilya V.; Rubio, Angel
2014-07-01
In this work, we give a comprehensive derivation of an exact and numerically feasible method to perform ab initio calculations of quantum particles interacting with a quantized electromagnetic field. We present a hierarchy of density-functional-type theories that describe the interaction of charged particles with photons and introduce the appropriate Kohn-Sham schemes. We show how the evolution of a system described by quantum electrodynamics in Coulomb gauge is uniquely determined by its initial state and two reduced quantities. These two fundamental observables, the polarization of the Dirac field and the vector potential of the photon field, can be calculated by solving two coupled, nonlinear evolution equations without the need to explicitly determine the (numerically infeasible) many-body wave function of the coupled quantum system. To find reliable approximations to the implicit functionals, we present the appropriate Kohn-Sham construction. In the nonrelativistic limit, this density-functional-type theory of quantum electrodynamics reduces to the density-functional reformulation of the Pauli-Fierz Hamiltonian, which is based on the current density of the electrons and the vector potential of the photon field. By making further approximations, e.g., restricting the allowed modes of the photon field, we derive further density-functional-type theories of coupled matter-photon systems for the corresponding approximate Hamiltonians. In the limit of only two sites and one mode we deduce the appropriate effective theory for the two-site Hubbard model coupled to one photonic mode. This model system is used to illustrate the basic ideas of a density-functional reformulation in great detail and we present the exact Kohn-Sham potentials for our coupled matter-photon model system.
NASA Astrophysics Data System (ADS)
Gong, Sai; Liu, Bang-Gui
2012-05-01
TiO2 has been recently used to realize high-temperature ferromagnetic semiconductors. In fact, it has been widely used for a long time as white pigment and sunscreen because of its whiteness, high refractive index, and excellent optical properties. However, its electronic structures and the related properties have not been satisfactorily understood. Here, we use Tran and Blaha's modified Becke-Johnson (TB-mBJ) exchange potential (plus a local density approximation correlation potential) within the density functional theory to investigate electronic structures and optical properties of rutile and anatase TiO2. Our comparative calculations show that the energy gaps obtained from mBJ method agree better with the experimental results than that obtained from local density approximation (LDA) and generalized gradient approximation (GGA), in contrast with substantially overestimated values from many-body perturbation (GW) calculations. As for optical dielectric functions (both real and imaginary parts), refractive index, and extinction coefficients as functions of photon energy, our mBJ calculated results are in excellent agreement with the experimental curves. Our further analysis reveals that these excellent improvements are achieved because mBJ potential describes accurately the energy levels of Ti 3d states. These results should be helpful to understand the high temperature ferromagnetism in doped TiO2. This approach can be used as a standard to understand electronic structures and the related properties of such materials as TiO2.
Electron affinities for rare gases and some actinides from local-spin-density-functional theory
Guo, Y.; Wrinn, M.C.; Whitehead, M.A. )
1989-12-01
The negative ions of the rare gases (He, Ne, Ar, Kr, Xe, and Rn) and some actinides (Pu, Am, Bk, Cf, and Es) have been calculated self-consistently by the generalized exchange local-spin-density-functional theory with self-interaction correction and correlation. The electron affinities were obtained as the differences between the statistical total energies of the negative ions and neutral atoms; the electron affinities were positive around several millirydbergs. Consequently, the negative ions are predicted stable for the rare gases and actinides.
A hybrid density functional study on the electron and hole trap states in anatase titanium dioxide.
Yamamoto, Takenori; Ohno, Takahisa
2012-01-14
We present a theoretical study on electron and hole trap states in the bulk and (001) surface of anatase titanium dioxide using screened hybrid density functional calculations. In both the bulk and surface, calculations suggest that the neutral and ionized oxygen vacancies are possible electron traps. The doubly ionized oxygen vacancy is the most stable in the bulk, and is a candidate for a shallow donor in colorless anatase crystals. The hole trap states are localized at oxygen anions in both the bulk and surface. The self-trapped electron centered at a titanium cation cannot be produced in the bulk, but can be formed at the surface. The electron trap level at the surface oxygen vacancy is consistent with observations by photoelectron spectroscopy. The optical absorptions and luminescence in UV-irradiated anatase nanoparticles are found to come from the surface self-trapped hole and the surface oxygen vacancy. PMID:22127526
Corsini, Niccolò R. C. Greco, Andrea; Haynes, Peter D.; Hine, Nicholas D. M.; Molteni, Carla
2013-08-28
We present an implementation in a linear-scaling density-functional theory code of an electronic enthalpy method, which has been found to be natural and efficient for the ab initio calculation of finite systems under hydrostatic pressure. Based on a definition of the system volume as that enclosed within an electronic density isosurface [M. Cococcioni, F. Mauri, G. Ceder, and N. Marzari, Phys. Rev. Lett.94, 145501 (2005)], it supports both geometry optimizations and molecular dynamics simulations. We introduce an approach for calibrating the parameters defining the volume in the context of geometry optimizations and discuss their significance. Results in good agreement with simulations using explicit solvents are obtained, validating our approach. Size-dependent pressure-induced structural transformations and variations in the energy gap of hydrogenated silicon nanocrystals are investigated, including one comparable in size to recent experiments. A detailed analysis of the polyamorphic transformations reveals three types of amorphous structures and their persistence on depressurization is assessed.
NASA Astrophysics Data System (ADS)
Dimakis, Nicholas; Navarro, Nestor E.; Velazquez, Julian; Salgado, Andres
2015-04-01
Periodic density functional calculations on graphene monolayers with and without an iron adatom have been used to elucidate iron-graphene adsorption and its effects on graphene electronic and vibrational properties. Density-of-states calculations and charge density contour plots reveal charge transfer from the iron s orbitals to the d orbitals, in agreement with past reports. Adsorbed iron atoms covalently bind to the graphene substrate, verified by the strong hybridization of iron d-states with the graphene bands in the energy region just below the Fermi level. This adsorption is weak and compared to the well-analyzed CO adsorption on Pt: It is indicated by its small adsorption energy and the minimal change of the substrate geometry due to the presence of the iron adatoms. Graphene vibrational spectra are analyzed though a systematic variation of the graphene supercell size. The shifts of graphene most prominent infrared active vibrational modes due to iron adsorption are explored using normal mode eigenvectors.
Extracting the density profile of an electronic wave function in a quantum dot
NASA Astrophysics Data System (ADS)
Boyd, Erin E.; Westervelt, Robert M.
2011-11-01
We use a model of a one-dimensional nanowire quantum dot to demonstrate the feasibility of a scanning probe microscope (SPM) imaging technique that can extract both the energy of an electron state and the amplitude of its wave function using a single instrument. This imaging technique can probe electrons that are buried beneath the surface of a low-dimensional semiconductor structure and provide valuable information for the design of quantum devices. A conducting SPM tip, acting as a movable gate, measures the energy of an electron state using Coulomb blockade spectroscopy. When the tip is close to the nanowire dot, it dents the wave function Ψ(x) of the quantum state, changing the electron's energy by an amount proportional to |Ψ(x)|2. By recording the change in energy as the SPM tip is moved along the length of the dot, the density profile of the electronic wave function can be found along the length of the quantum dot.
NASA Astrophysics Data System (ADS)
Basurto, Luis; Zope, Rajendra R.; Baruah, Tunna
2016-05-01
We report an electronic structure study of a multichromophoric molecular complex containing two of each borondipyrromethane dye, Zn-tetraphenyl-porphyrin, bisphenyl anthracene and a fullerene. The snowflake shaped molecule behaves like an antenna capturing photon at different frequencies and transferring the photon energy to the porphyrin where electron transfer occurs from the porphyrin to the fullerene. The study is performed within density functional formalism using large polarized Guassian basis sets (12,478 basis functions in total). The energies of the HOMO and LUMO states in the complex, as adjudged by the ionization potential and the electron affinity values, show significant differences with respect to their values in participating subunits in isolation. These differences are also larger than the variations of the ionization potential and electron affinity values observed in non-bonded C60-ZnTPP complexes in co-facial arrangement or end-on orientations. An understanding of the origin of these differences is obtained by a systematic study of the effect of structural strain, the presence of ligands, the effect of orbital delocalization on the ionization energy and the electron affinity. Finally, a few lowest charge transfer energies involving electronic transitions from the porphyrin component to the fullerene subunit of the complex are predicted.
Applications of methods beyond density functional theory to the study of correlated electron systems
NASA Astrophysics Data System (ADS)
Sims, Hunter Robert
The difficulty in accurately treating systems in which electron-electron interactions are the dominant physics has plagued condensed matter physics for decades. Currently, there exist many different computational techniques designed to improve upon density functional theory to varying degrees of accuracy. To date, no unified, parameter-free method exists that is guaranteed to yield the correct answer for all materials. Consequently, proper treatment of such systems often requires a combination of several methods, allowing one to check them against one another when their regions of validity overlap and to expand one's reach when a single method cannot reliably describe all of the physics at work. In this dissertation, I present discussion and, when appropriate, brief derivations of several of the most prominent electronic structure methods currently in use---from the local density approximation through LDA+DMFT. I then present several investigations into the electronic and magnetic structure of materials of potential interest for information technology that also illustrate the current state of affairs in computational condensed matter physics. I explore the intersite exchange interactions in CrO2 within density functional theory (with and without Hubbard "+U" corrections) and evaluate these results through analytic and numerical means. I study the dependence of the mysterious magnetization of Fe16N2 on crystal and electronic structure and employ a wide range of techniques in an attempt to bring greater rigor and deeper understanding to the widely-varying reports on this material. In conjunction with others' careful experimental analysis, I provide a picture of the band structure of the magnetic insulator NiFe2O4 that reveals a novel hierarchy in its band gaps and suggests applications in spintronics and possibly other areas. Finally, I employ dynamical mean-field theory to study the behavior of impurity states in elemental semiconductors, using H impurities in Ge as
Density functional study on electronic properties of P-doped spinel silicon carbon nitride
Zhang Yufen; Zhao Xian Cheng Xiufeng; Mu Yuguang
2008-08-15
We performed density functional calculations on the electronic properties of P-doped spinel silicon carbon nitride. When Si is replaced by C at the tetrahedral sites of P-doped c-Si{sub 3}N{sub 4}, the band gap can be adjusted, and an insulator-to-metal transition is predicted to occur at the C-to-Si ratio of 0.27. Finally, some possible examinations and potential applications for the large band-gap reduction are discussed. - Graphical abstract: We performed density functional calculations to predict the insulator-to-metal transition by replacing Si by C at the tetrahedral sites of P-doped c-Si{sub 3}N{sub 4}.
NASA Astrophysics Data System (ADS)
Prasad, O.; Sinha, L.; Misra, N.; Narayan, V.; Kumar, N.; Kumar, A.
2010-09-01
The present work deals with the structural, electronic, and vibrational analysis of rivastigmine. Rivastigmine, an antidementia medicament, is credited with significant therapeutic effects on the cognitive, functional, and behavioural problems that are commonly associated with Alzheimer’s dementia. For rivastigmine, a number of minimum energy conformations are possible. The geometry of twelve possible conformers has been analyzed and the most stable conformer was further optimized at a higher basis set. The electronic properties and vibrational frequencies were then calculated using a density functional theory at the B3LYP level with the 6-311+G(d, p) basis set. The different molecular surfaces have also been drawn to understand the activity of the molecule. A narrower frontier orbital energy gap in rivastigmine makes it softer and more reactive than water and dimethylfuran. The calculated value of the dipole moment is 2.58 debye.
Li, Zun-Yun; Wang, Hai-Long; Lu, Tong-Tong; He, Tian-Jing; Liu, Fan-Chen; Chen, Dong-Ming
2007-08-01
The ground-state structure and electronic and vibrational spectra of octaethylporphyrin diacid (H4OEP2+) have been studied with the density functional theory. The geometrical parameters computed with B3LYP, PBE1PBE and mPW1PW91 functionals and 6-31G* basis sets are well consistent with the experimental values. Electronic absorption spectrum of H4OEP2+ has been studied with the time-dependent DFT method, and the calculated excitation energies and oscillator strengths are compared with the experimental results. The Raman and IR spectra of H4OEP2+ and the Raman spectrum of its N-deuterated analogue (D4OEP2+) were measured. The observed Raman and IR bands have been assigned based on the frequency calculations at the B3LYP/6-31G* level of theory. PMID:17142087
Hao, Feng Mattsson, Ann E.; Armiento, Rickard
2014-05-14
We have previously proposed that further improved functionals for density functional theory can be constructed based on the Armiento-Mattsson subsystem functional scheme if, in addition to the uniform electron gas and surface models used in the Armiento-Mattsson 2005 functional, a model for the strongly confined electron gas is also added. However, of central importance for this scheme is an index that identifies regions in space where the correction provided by the confined electron gas should be applied. The electron localization function (ELF) is a well-known indicator of strongly localized electrons. We use a model of a confined electron gas based on the harmonic oscillator to show that regions with high ELF directly coincide with regions where common exchange energy functionals have large errors. This suggests that the harmonic oscillator model together with an index based on the ELF provides the crucial ingredients for future improved semi-local functionals. For a practical illustration of how the proposed scheme is intended to work for a physical system we discuss monoclinic cupric oxide, CuO. A thorough discussion of this system leads us to promote the cell geometry of CuO as a useful benchmark for future semi-local functionals. Very high ELF values are found in a shell around the O ions, and take its maximum value along the Cu–O directions. An estimate of the exchange functional error from the effect of electron confinement in these regions suggests a magnitude and sign that could account for the error in cell geometry.
Actinide electronic structure based on the Dirac equation and density functional theory
NASA Astrophysics Data System (ADS)
Wills, John M.; Mattsson, Ann E.
2013-03-01
Density functional theory (DFT) provides a formally predictive basis for predicting the structural properties of actinides. Although available approximations to the exchange/correlation functional provide accurate predictions for many materials, they fail qualitatively and sometimes quantitatively when applied to actinides. Major contributors to this deficiency are an inadequate treatment of confinement physics and an incomplete treatment of relativity in the underlying equations. The development of a functional correctly incorporating confinement physics with a proper treatment of relativity would provide definitive, internally consistent predictions of actinide properties. To enable the development of such a functional and quantify the predictions of currently available functionals, we have developed an efficient first-principles electronic structure method based on the Dirac equation. Results are compared with current methods, and the implications for relativistic density functionals discussed. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.
Dale, Stephen G.; Johnson, Erin R.
2015-11-14
Exploration of the solvated electron phenomena using density-functional theory (DFT) generally results in prediction of a localised electron within an induced solvent cavity. However, it is well known that DFT favours highly delocalised charges, rendering the localisation of a solvated electron unexpected. We explore the origins of this counterintuitive behaviour using a model Kevan-structure system. When a polarisable-continuum solvent model is included, it forces electron localisation by introducing a strong energetic bias that favours integer charges. This results in the formation of a large energetic barrier for charge-hopping and can cause the self-consistent field to become trapped in local minima thus converging to stable solutions that are higher in energy than the ground electronic state. Finally, since the bias towards integer charges is caused by the polarisable continuum, these findings will also apply to other classical polarisation corrections, as in combined quantum mechanics and molecular mechanics (QM/MM) methods. The implications for systems beyond the solvated electron, including cationic DNA bases, are discussed.
Two-electron Rabi oscillations in real-time time-dependent density-functional theory.
Habenicht, Bradley F; Tani, Noriyuki P; Provorse, Makenzie R; Isborn, Christine M
2014-11-14
We investigate the Rabi oscillations of electrons excited by an applied electric field in several simple molecular systems using time-dependent configuration interaction (TDCI) and real-time time-dependent density-functional theory (RT-TDDFT) dynamics. While the TDCI simulations exhibit the expected single-electron Rabi oscillations at a single resonant electric field frequency, Rabi oscillations in the RT-TDDFT simulations are a two-electron process. The existence of two-electron Rabi oscillations is determined both by full population inversion between field-free molecular orbitals and the behavior of the instantaneous dipole moment during the simulations. Furthermore, the Rabi oscillations in RT-TDDFT are subject to an intensity threshold of the electric field, below which Rabi oscillations do not occur and above which the two-electron Rabi oscillations occur at a broad range of frequencies. It is also shown that at field intensities near the threshold intensity, the field frequency predicted to induce Rabi oscillations by linear response TDDFT only produces detuned Rabi oscillations. Instead, the field frequency that yields the full two-electron population inversion and Rabi oscillation behavior is shown to be the average of single-electron transition frequencies from the ground S0 state and the doubly-excited S2 state. The behavior of the two-electron Rabi oscillations is rationalized via two possible models. The first model is a multi-photon process that results from the electric field interacting with the three level system such that three level Rabi oscillations may occur. The second model suggests that the mean-field nature of RT-TDDFT induces paired electron propagation. PMID:25399137
Two-electron Rabi oscillations in real-time time-dependent density-functional theory
Habenicht, Bradley F.; Tani, Noriyuki P.; Provorse, Makenzie R.; Isborn, Christine M.
2014-11-14
We investigate the Rabi oscillations of electrons excited by an applied electric field in several simple molecular systems using time-dependent configuration interaction (TDCI) and real-time time-dependent density-functional theory (RT-TDDFT) dynamics. While the TDCI simulations exhibit the expected single-electron Rabi oscillations at a single resonant electric field frequency, Rabi oscillations in the RT-TDDFT simulations are a two-electron process. The existence of two-electron Rabi oscillations is determined both by full population inversion between field-free molecular orbitals and the behavior of the instantaneous dipole moment during the simulations. Furthermore, the Rabi oscillations in RT-TDDFT are subject to an intensity threshold of the electric field, below which Rabi oscillations do not occur and above which the two-electron Rabi oscillations occur at a broad range of frequencies. It is also shown that at field intensities near the threshold intensity, the field frequency predicted to induce Rabi oscillations by linear response TDDFT only produces detuned Rabi oscillations. Instead, the field frequency that yields the full two-electron population inversion and Rabi oscillation behavior is shown to be the average of single-electron transition frequencies from the ground S{sub 0} state and the doubly-excited S{sub 2} state. The behavior of the two-electron Rabi oscillations is rationalized via two possible models. The first model is a multi-photon process that results from the electric field interacting with the three level system such that three level Rabi oscillations may occur. The second model suggests that the mean-field nature of RT-TDDFT induces paired electron propagation.
Hubert, Mickaël; Hedegård, Erik D; Jensen, Hans Jørgen Aa
2016-05-10
Computational methods that can accurately and effectively predict all types of electronic excitations for any molecular system are missing in the toolbox of the computational chemist. Although various Kohn-Sham density-functional methods (KS-DFT) fulfill this aim in some cases, they become inadequate when the molecule has near-degeneracies and/or low-lying double-excited states. To address these issues we have recently proposed multiconfiguration short-range density-functional theory-MC-srDFT-as a new tool in the toolbox. While initial applications for systems with multireference character and double excitations have been promising, it is nevertheless important that the accuracy of MC-srDFT is at least comparable to the best KS-DFT methods also for organic molecules that are typically of single-reference character. In this paper we therefore systematically investigate the performance of MC-srDFT for a selected benchmark set of electronic excitations of organic molecules, covering the most common types of organic chromophores. This investigation confirms the expectation that the MC-srDFT method is accurate for a broad range of excitations and comparable to accurate wave function methods such as CASPT2, NEVPT2, and the coupled cluster based CC2 and CC3. PMID:27058733
Double-hybrid density functional theory for excited electronic states of molecules
NASA Astrophysics Data System (ADS)
Grimme, Stefan; Neese, Frank
2007-10-01
Double-hybrid density functionals are based on a mixing of standard generalized gradient approximations (GGAs) for exchange and correlation with Hartree-Fock (HF) exchange and a perturbative second-order correlation part (PT2) that is obtained from the Kohn-Sham (GGA) orbitals and eigenvalues. This virtual orbital-dependent functional (dubbed B2PLYP) contains only two empirical parameters that describe the mixture of HF and GGA exchange (ax) and of the PT2 and GGA correlation (ac), respectively. Extensive testing has recently demonstrated the outstanding accuracy of this approach for various ground state problems in general chemistry applications. The method is extended here without any further empirical adjustments to electronically excited states in the framework of time-dependent density functional theory (TD-DFT) or the closely related Tamm-Dancoff approximation (TDA-DFT). In complete analogy to the ground state treatment, a scaled second-order perturbation correction to configuration interaction with singles (CIS(D)) wave functions developed some years ago by Head-Gordon et al. [Chem. Phys. Lett. 219, 21 (1994)] is computed on the basis of density functional data and added to the TD(A)-DFT/GGA excitation energy. The method is implemented by applying the resolution of the identity approximation and the efficiency of the code is discussed. Extensive tests for a wide variety of molecules and excited states (of singlet, triplet, and doublet multiplicities) including electronic spectra are presented. In general, rather accurate excitation energies (deviations from reference data typically <0.2eV) are obtained that are mostly better than those from standard functionals. Still, systematic errors are obtained for Rydberg (too low on average by about 0.3eV) and charge-transfer transitions but due to the relatively large ax parameter (0.53), B2PLYP outperforms most other functionals in this respect. Compared to conventional HF-based CIS(D), the method is more robust in
NASA Astrophysics Data System (ADS)
Khaira, Jobanpreet S.; Jain, Richa N.; Chakraborty, Brahmananda; Ramaniah, Lavanya M.
2015-06-01
The electronic structure of yttrium-doped Silicon Carbide Nanotubes has been theoretically investigated using first principles density functional theory (DFT). Yttrium atom is bonded strongly on the surface of the nanotube with a binding energy of 2.37 eV and prefers to stay on the hollow site at a distance of around 2.25 Å from the tube. The semi-conducting nanotube with chirality (4, 4) becomes half mettalic with a magnetic moment of 1.0 µB due to influence of Y atom on the surface. There is strong hybridization between d orbital of Y with p orbital of Si and C causing a charge transfer from d orbital of the Y atom to the tube. The Fermi level is shifted towards higher energy with finite Density of States for only upspin channel making the system half metallic and magnetic which may have application in spintronic devices.
Khaira, Jobanpreet S.; Jain, Richa N.; Chakraborty, Brahmananda; Ramaniah, Lavanya M.
2015-06-24
The electronic structure of yttrium-doped Silicon Carbide Nanotubes has been theoretically investigated using first principles density functional theory (DFT). Yttrium atom is bonded strongly on the surface of the nanotube with a binding energy of 2.37 eV and prefers to stay on the hollow site at a distance of around 2.25 Å from the tube. The semi-conducting nanotube with chirality (4, 4) becomes half mettalic with a magnetic moment of 1.0 µ{sub B} due to influence of Y atom on the surface. There is strong hybridization between d orbital of Y with p orbital of Si and C causing a charge transfer from d orbital of the Y atom to the tube. The Fermi level is shifted towards higher energy with finite Density of States for only upspin channel making the system half metallic and magnetic which may have application in spintronic devices.
Conformation and electronic structure of polyethylene: A density-functional approach
NASA Astrophysics Data System (ADS)
Miao, M. S.; van Camp, P. E.; van Doren, V. E.; Ladik, J. J.; Mintmire, J. W.
1996-10-01
Two different local-density approximations, the Gáspár-Kohn-Sham and the Perdew-Zunger approximations, of the density-functional method have been used to calculate structural and electronic properties of polyethylene systems with several different dihedral angles. For each system, the CC bond lengths and the CCC and HCH bond angles are optimized simultaneously. All the parameters appear to be strongly coupled with torsional freedom and vary with the change in dihedral angle in a pattern similar to that of the total energy. The total energy has an absolute minimum for the planar zigzag conformation but a distinct local minimum for the quasistable helical conformation. Another minimum occurs in the energy curve close to this gauche minimum. The calculated valence and conduction bands are discussed and compared with other theoretical calculations and experiment.
Sjostrom, Travis; Daligault, Jérôme
2015-12-01
We validate the application of our recent orbital-free density functional theory (DFT) approach [Phys. Rev. Lett. 113, 155006 (2014);] for the calculation of ionic and electronic transport properties of dense plasmas. To this end, we calculate the self-diffusion coefficient, the viscosity coefficient, the electrical and thermal conductivities, and the reflectivity coefficient of hydrogen and aluminum plasmas. Very good agreement is found with orbital-based Kohn-Sham DFT calculations at lower temperatures. Because the computational costs of the method do not increase with temperature, we can produce results at much higher temperatures than is accessible by the Kohn-Sham method. Our results for warm dense aluminum at solid density are inconsistent with the recent experimental results reported by Sperling et al. [Phys. Rev. Lett. 115, 115001 (2015)]. PMID:26764850
Ionic and electronic transport properties in dense plasmas by orbital-free density functional theory
NASA Astrophysics Data System (ADS)
Sjostrom, Travis; Daligault, Jérôme
2015-12-01
We validate the application of our recent orbital-free density functional theory (DFT) approach [Phys. Rev. Lett. 113, 155006 (2014), 10.1103/PhysRevLett.113.155006;] for the calculation of ionic and electronic transport properties of dense plasmas. To this end, we calculate the self-diffusion coefficient, the viscosity coefficient, the electrical and thermal conductivities, and the reflectivity coefficient of hydrogen and aluminum plasmas. Very good agreement is found with orbital-based Kohn-Sham DFT calculations at lower temperatures. Because the computational costs of the method do not increase with temperature, we can produce results at much higher temperatures than is accessible by the Kohn-Sham method. Our results for warm dense aluminum at solid density are inconsistent with the recent experimental results reported by Sperling et al. [Phys. Rev. Lett. 115, 115001 (2015), 10.1103/PhysRevLett.115.115001].
NASA Astrophysics Data System (ADS)
Kulwinder, Kaur; Ranjan, Kumar
2016-05-01
We study the effect of pressure on electronic and thermoelectric properties of Mg2Si using the density functional theory and Boltzmann transport equations. The variation of lattice constant, band gap, bulk modulus with pressure is also analyzed. Further, the thermoelectric properties (Seebeck coefficient, electrical conductivity, electronic thermal conductivity) have been studied as a function of temperature and pressure up to 1200 K. The results show that Mg2Si is an n-type semiconductor with a band gap of 0.21 eV. The negative value of the Seebeck coefficient at all pressures indicates that the conduction is due to electrons. With the increase in pressure, the Seebeck coefficient decreases and electrical conductivity increases. It is also seen that, there is practically no effect of pressure on the electronic contribution of thermal conductivity. The paper describes the calculation of the lattice thermal conductivity and figure of merit of Mg2Si at zero pressure. The maximum value of figure of merit is attained 1.83×10‑3 at 1000 K. The obtained results are in good agreement with the available experimental and theoretical results. Project supported by the Council of Scientific & Industrial Research (CSIR), India.
Ionic thermal effects on photo-electron emission within time-dependent density-functional theory
NASA Astrophysics Data System (ADS)
Gao, Cong-Zhang; Dinh, Phuong Mai; Reinhard, Paul-Gerhard; Suraud, Eric
2016-02-01
We study the impact of thermal fluctuations of cluster/molecule shape on photo-electron spectra (PES) and photo-electron angular distributions (PAD) using a detailed time-dependent simulation of the emission dynamics and thermal ionic motion. Basis of the description is time-dependent density-functional theory (TDDFT) coupled to molecular dynamics for ionic motion. Test cases are small Na clusters and the C3 molecule. For Na clusters, we find that PES signals are rather robust for one-photon processes while large smearing of the pattern are observed at lower frequencies in multi-photon processes. This effect can be related to the typical spectral response of the metal clusters. PAD are generally much more robust than PES. The C3 molecule produces a greater variety of thermal response. This happens because this molecule has eigenmodes with much different softness.
Electronic and magnetic properties of substituted BN sheets: A density functional theory study
Zhou, Yungang; Yang, Ping; Wang, Zhiguo; Zu, Xiaotao T.; Xiao, Hai Yan; Sun, Xin; Khaleel, Mohammad A.; Gao, Fei
2011-04-15
Using density functional calculations, we investigate the geometries, electronic structures and magnetic properties of hexagonal BN sheets with 3d transition metal (TM) and nonmetal atoms embedded in three types of vacancies: VB, VN, and VB+N. We show that some embedded configurations, except TM atoms in VN vacancy, are stable in BN sheet and yield interesting phenomena. For instance, the band gaps and magnetic moments of BN sheet can be tuned depending on the embedded dopant species and vacancy type. In particular, embedment such as Cr in VB+N, Co in VB, and Ni in VB leads to half-metallic BN sheets interesting for spin filter applications. From the investigation of Mn-chain (CMn) embedments, a regular 1D structure can be formed in BN sheet as an electron waveguide, a metal nanometer wire with a single atom thickness.
NASA Astrophysics Data System (ADS)
Xiao, Ling-Ping; Zeng, Zhi; Chen, Xiao-Jia
2016-06-01
The pressure effect on the geometrical and electronic structures of crystalline naphthalene is calculated up to 30 GPa by performing density functional calculations. The lattice parameters a, b, and c, decrease by 1.77 Å (-20.4%), 0.85 Å (-14.1%), and 0.91 Å (-8.2%), respectively, while the monoclinic angle β increases by 3.95° in this pressure region. At the highest pressure of 30 GPa the unit cell volume decreases by 62.7%. The detailed analysis of the molecular arrangement within crystal structure reveals that the molecular motion becomes more and more localized, and hints towards the evolution of intermolecular interaction with pressure. Moreover, the electronic structure of naphthalene under high pressure is also discussed. A pressure induced decrease of the band gap is observed.
NASA Astrophysics Data System (ADS)
Almeida, Nuno M. S.; McKinlay, Russell G.; Paterson, Martin J.
2015-01-01
The lowest excited electronic states of the permanganate ion MnO4- are calculated using a hierarchy of coupled cluster response approaches, as well as time-dependent density functional theory. It is shown that while full linear response coupled cluster with singles and doubles (or higher) performs well, that permanganate represents a stern test for approximate coupled cluster response models, and that problems can be traced to very large orbital relaxation effects. TD-DFT is reasonably robust although errors around 0.6 eV are still observed. In order to further investigate the strong correlations prevalent in the electronic ground state large-scale RASSCF calculations were also performed. Again very large orbital relaxation in the correlated wavefunction is observed. Although the system can qualitatively be described by a single configuration, multi-reference diagnostic values show that care must be taken in this and similar metal complexes.
Quantum chemistry without wave functions: two-electron reduced density matrices.
Mazziotti, David A
2006-03-01
For 50 years, progress toward the direct calculation of the ground-state two-electron reduced density matrix (2-RDM) was stymied from an inability to constrain the 2-RDM to represent an N-electron wave function. Recent advances in theory and optimization have realized the direct calculation of the 2-RDM. A variational 2-RDM procedure, using first-order semidefinite programming, has been shown to capture multireference correlation effects important at nonequilibrium geometries [Mazziotti, Phys. Rev. Lett. 2004, 93, No. 213001]. This method emerged from research on a nonvariational calculation of the 2-RDM by the contracted Schrödinger equation. Both approaches will be discussed and illustrated. PMID:16548509
Electron bubbles in liquid helium: Density functional calculations of infrared absorption spectra
Grau, Victor; Barranco, Manuel; Mayol, Ricardo; Pi, Marti
2006-02-01
Within density functional theory, we have calculated the energy of the transitions from the ground state to the first two excited states in the electron bubbles in liquid helium at pressures from zero to about the solidification pressure. For {sup 4}He at low temperatures, our results are in very good agreement with infrared absorption experiments. Above a temperature of {approx}2 K, we overestimate the energy of the 1s-1p transition. We attribute this to the break down of the Franck-Condon principle due to the presence of helium vapor inside the bubble. Our results indicate that the 1s-2p transition energies are sensitive not only to the size of the electron bubble, but also to its surface thickness. We also present results for the infrared transitions in the case of liquid {sup 3}He, for which we lack experimental data.
A revised electronic Hessian for approximate time-dependent density functional theory.
Ziegler, Tom; Seth, Michael; Krykunov, Mykhaylo; Autschbach, Jochen
2008-11-14
Time-dependent density functional theory (TD-DFT) at the generalized gradient level of approximation (GGA) has shown systematic errors in the calculated excitation energies. This is especially the case for energies representing electron transitions between two separated regions of space or between orbitals of different spatial extents. It will be shown that these limitations can be attributed to the electronic ground state Hessian G(GGA). Specifically, we shall demonstrate that the Hessian G(GGA) can be used to describe changes in energy due to small perturbations of the electron density (Deltarho), but it should not be applied to one-electron excitations involving the density rearrangement (Deltarho) of a full electron charge. This is in contrast to Hartree-Fock theory where G(HF) has a trust region that is accurate for both small perturbations and one-electron excitations. The large trust radius of G(HF) can be traced back to the complete cancellation of Coulomb and exchange terms in Hartree-Fock (HF) theory representing self-interaction (complete self-interaction cancellation, CSIC). On the other hand, it is shown that the small trust radius for G(GGA) can be attributed to the fact that CSIC is assumed for GGA in the derivation of G(GGA) although GGA (and many other approximate DFT schemes) exhibits incomplete self-interaction cancellation (ISIC). It is further shown that one can derive a new matrix G(R-DFT) with the same trust region as G(HF) by taking terms due to ISIC properly into account. Further, with TD-DFT based on G(R-DFT), energies for state-to-state transitions represented by a one-electron excitation (psi(i)-->psi(a)) are approximately calculated as DeltaE(ai). Here DeltaE(ai) is the energy difference between the ground state Kohn-Sham Slater determinant and the energy of a Kohn-Sham Slater determinant where psi(i) has been replaced by psi(a). We make use of the new Hessian in two numerical applications involving charge-transfer excitations. It is
Mezey, Paul G
2014-09-16
Conspectus Just as complete molecules have no boundaries and have "fuzzy" electron density clouds approaching zero density exponentially at large distances from the nearest nucleus, a physically justified choice for electron density fragments exhibits similar behavior. Whereas fuzzy electron densities, just as any fuzzy object, such as a thicker cloud on a foggy day, do not lend themselves to easy visualization, one may partially overcome this by using isocontours. Whereas a faithful representation of the complete fuzzy density would need infinitely many such isocontours, nevertheless, by choosing a selected few, one can still obtain a limited pictorial representation. Clearly, such images are of limited value, and one better relies on more complete mathematical representations, using, for example, density matrices of fuzzy fragment densities. A fuzzy density fragmentation can be obtained in an exactly additive way, using the output from any of the common quantum chemical computational techniques, such as Hartree-Fock, MP2, and various density functional approaches. Such "fuzzy" electron density fragments properly represented have proven to be useful in a rather wide range of applications, for example, (a) using them as additive building blocks leading to efficient linear scaling macromolecular quantum chemistry computational techniques, (b) the study of quantum chemical functional groups, (c) using approximate fuzzy fragment information as allowed by the holographic electron density theorem, (d) the study of correlations between local shape and activity, including through-bond and through-space components of interactions between parts of molecules and relations between local molecular shape and substituent effects, (e) using them as tools of density matrix extrapolation in conformational changes, (f) physically valid averaging and statistical distribution of several local electron densities of common stoichiometry, useful in electron density databank mining, for
Technology Transfer Automated Retrieval System (TEKTRAN)
Density functional theory/B3LYP has been employed to optimize the conformations of selected 4-arylflavan-3-ols and their phenolic methyl ether 3-O-acetates. The electronic circular dichroism spectra of the major conformers have been calculated using time-dependent density functional theory to valida...
NASA Astrophysics Data System (ADS)
Kweon, Kyoung E.; Hwang, Gyeong S.
2012-10-01
The structure and property prediction of metal oxides can significantly be improved by incorporating exact Hartree-Fock (HF) exchange into density functional theory (DFT), which is the so-called hybrid DFT. We explored the impact of HF exchange inclusion on the predicted structural, bonding, and electronic properties of bismuth vanadate (BiVO4), with particular attention to the difference between its monoclinic and tetragonal scheelite phases. The applied exchange-correlation (xc) functionals include the gradient corrected Perdew-Burke-Ernzerhof (PBE) and the PBE-HF hybrid functionals with HF exchange amounts of 10%, 25%, and 50%. We find that the PBE-HF25% yields a monoclinic structure in very close agreement with the experimentally determined structure, while the PBE-HF50% tends to overestimate the monoclinic distortion and the PBE/PBE-HF10% can hardly identify a distinct monoclinic configuration at ambient conditions. Electronic structure analysis reveals that the increasing monoclinic distortion with the amount of HF exchange is related to the enhancement of hybridization between Bi 6s-O 2p antibonding states and unoccupied Bi 6p states. The bonding mechanisms and band structures of the monoclinic and tetragonal phases of BiVO4 were also investigated, and we discuss how the predictions are sensitive to the xc functional choice.
A density functional theory investigation of the electronic structure and spin moments of magnetite
NASA Astrophysics Data System (ADS)
Noh, Junghyun; Osman, Osman I.; Aziz, Saadullah G.; Winget, Paul; Brédas, Jean-Luc
2014-08-01
We present the results of density functional theory (DFT) calculations on magnetite, Fe3O4, which has been recently considered as electrode in the emerging field of organic spintronics. Given the nature of the potential applications, we evaluated the magnetite room-temperature cubic Fd\\bar{3}m phase in terms of structural, electronic, and magnetic properties. We considered GGA (PBE), GGA + U (PBE + U), and range-separated hybrid (HSE06 and HSE(15%)) functionals. Calculations using HSE06 and HSE(15%) functionals underline the impact that inclusion of exact exchange has on the electronic structure. While the modulation of the band gap with exact exchange has been seen in numerous situations, the dramatic change in the valence band nature and states near the Fermi level has major implications for even a qualitative interpretation of the DFT results. We find that HSE06 leads to highly localized states below the Fermi level while HSE(15%) and PBE + U result in delocalized states around the Fermi level. The significant differences in local magnetic moments and atomic charges indicate that describing room-temperature bulk materials, surfaces and interfaces may require different functionals than their low-temperature counterparts.
Lowering the density of electronic defects on organic-functionalized Si(100) surfaces
Peng, Weina DeBenedetti, William J. I.; Kim, Seonjae; Chabal, Yves J.; Hines, Melissa A.
2014-06-16
The electrical quality of functionalized, oxide-free silicon surfaces is critical for chemical sensing, photovoltaics, and molecular electronics applications. In contrast to Si/SiO{sub 2} interfaces, the density of interface states (D{sub it}) cannot be reduced by high temperature annealing because organic layers decompose above 300 °C. While a reasonable D{sub it} is achieved on functionalized atomically flat Si(111) surfaces, it has been challenging to develop successful chemical treatments for the technologically relevant Si(100) surfaces. We demonstrate here that recent advances in the chemical preparation of quasi-atomically-flat, H-terminated Si(100) surfaces lead to a marked suppression of electronic states of functionalized surfaces. Using a non-invasive conductance-voltage method to study functionalized Si(100) surfaces with varying roughness, a D{sub it} as low as 2.5 × 10{sup 11} cm{sup −2}eV{sup −1} is obtained for the quasi-atomically-flat surfaces, in contrast to >7 × 10{sup 11} cm{sup −2}eV{sup −1} on atomically rough Si(100) surfaces. The interfacial quality of the organic/quasi-atomically-flat Si(100) interface is very close to that obtained on organic/atomically flat Si(111) surfaces, opening the door to applications previously thought to be restricted to Si(111)
Lowering the density of electronic defects on organic-functionalized Si(100) surfaces
NASA Astrophysics Data System (ADS)
Peng, Weina; DeBenedetti, William J. I.; Kim, Seonjae; Hines, Melissa A.; Chabal, Yves J.
2014-06-01
The electrical quality of functionalized, oxide-free silicon surfaces is critical for chemical sensing, photovoltaics, and molecular electronics applications. In contrast to Si/SiO2 interfaces, the density of interface states (Dit) cannot be reduced by high temperature annealing because organic layers decompose above 300 °C. While a reasonable Dit is achieved on functionalized atomically flat Si(111) surfaces, it has been challenging to develop successful chemical treatments for the technologically relevant Si(100) surfaces. We demonstrate here that recent advances in the chemical preparation of quasi-atomically-flat, H-terminated Si(100) surfaces lead to a marked suppression of electronic states of functionalized surfaces. Using a non-invasive conductance-voltage method to study functionalized Si(100) surfaces with varying roughness, a Dit as low as 2.5 × 1011 cm-2eV-1 is obtained for the quasi-atomically-flat surfaces, in contrast to >7 × 1011 cm-2eV-1 on atomically rough Si(100) surfaces. The interfacial quality of the organic/quasi-atomically-flat Si(100) interface is very close to that obtained on organic/atomically flat Si(111) surfaces, opening the door to applications previously thought to be restricted to Si(111).
Using time-dependent density functional theory in real time for calculating electronic transport
NASA Astrophysics Data System (ADS)
Schaffhauser, Philipp; Kümmel, Stephan
2016-01-01
We present a scheme for calculating electronic transport within the propagation approach to time-dependent density functional theory. Our scheme is based on solving the time-dependent Kohn-Sham equations on grids in real space and real time for a finite system. We use absorbing and antiabsorbing boundaries for simulating the coupling to a source and a drain. The boundaries are designed to minimize the effects of quantum-mechanical reflections and electrical polarization build-up, which are the major obstacles when calculating transport by applying an external bias to a finite system. We show that the scheme can readily be applied to real molecules by calculating the current through a conjugated molecule as a function of time. By comparing to literature results for the conjugated molecule and to analytic results for a one-dimensional model system we demonstrate the reliability of the concept.
NASA Astrophysics Data System (ADS)
Sun, Haitao; Tang, Ke; Li, Yanmin; Su, Chunfang; Zhou, Zhengyu; Wang, Zhizhong
The effect of hydrogen bond interactions on ionization potentials (IPs) and electron affinities (EAs) of thymine-formamide complexes (T-F) have been investigated employing the density functional theory B3LYP at 6-311++G(d, p) basis set level. All complexes experience a geometrical change on either electron detachment or attachment, and the change might be facilitated or hindered according to the strength of the hydrogen-bonding interaction involved. The strength of hydrogen bonds presents an opposite changing trend on the two processes. A more important role that H-bonding interaction plays in the process of electron attachment than in the process of electron detachment can be seen by a comparison of the IPs and EAs of complexes with that of isolated thymine. Futhermore, the EAs of isolated thymine are in good agreement with the experimental values (AEA is 0.79 eV, VEA is -0.29 eV [Wetmore et al., Chem Phys Lett 2000, 322, 129]). The calculated total NPA charge distributions reveal that nearly all the negative charges locate on thymine monomer in the anions and even in the cationic states, there are a few negative charges on thymine monomer. An analysis of dissociation energies predicts the processes T-F+→ T++ F and T-F- → T- + F to be the most energetically favorable for T-F+ and T-F-, respectively. Content:text/plain; charset="UTF-8"
Electronic transport properties of one dimensional lithium nanowire using density functional theory
Thakur, Anil; Kumar, Arun; Chandel, Surjeet; Ahluwalia, P. K.
2015-05-15
Single nanowire electrode devices are a unique platform for studying as energy storage devices. Lithium nanowire is of much importance in lithium ion batteries and therefore has received a great deal of attention in past few years. In this paper we investigated structural and electronic transport properties of Li nanowire using density functional theory (DFT) with SIESTA code. Electronic transport properties of Li nanowire are investigated theoretically. The calculations are performed in two steps: first an optimized geometry for Li nanowire is obtained using DFT calculations, and then the transport relations are obtained using NEGF approach. SIESTA and TranSIESTA simulation codes are used in the calculations correspondingly. The electrodes are chosen to be the same as the central region where transport is studied, eliminating current quantization effects due to contacts and focusing the electronic transport study to the intrinsic structure of the material. By varying chemical potential in the electrode regions, an I-V curve is traced which is in agreement with the predicted behavior. Agreement of bulk properties of Li with experimental values make the study of electronic and transport properties in lithium nanowires interesting because they are promising candidates as bridging pieces in nanoelectronics. Transmission coefficient and V-I characteristic of Li nano wire indicates that Li nanowire can be used as an electrode device.
Electronic transport properties of one dimensional lithium nanowire using density functional theory
NASA Astrophysics Data System (ADS)
Thakur, Anil; Kumar, Arun; Chandel, Surjeet; Ahluwalia, P. K.
2015-05-01
Single nanowire electrode devices are a unique platform for studying as energy storage devices. Lithium nanowire is of much importance in lithium ion batteries and therefore has received a great deal of attention in past few years. In this paper we investigated structural and electronic transport properties of Li nanowire using density functional theory (DFT) with SIESTA code. Electronic transport properties of Li nanowire are investigated theoretically. The calculations are performed in two steps: first an optimized geometry for Li nanowire is obtained using DFT calculations, and then the transport relations are obtained using NEGF approach. SIESTA and TranSIESTA simulation codes are used in the calculations correspondingly. The electrodes are chosen to be the same as the central region where transport is studied, eliminating current quantization effects due to contacts and focusing the electronic transport study to the intrinsic structure of the material. By varying chemical potential in the electrode regions, an I-V curve is traced which is in agreement with the predicted behavior. Agreement of bulk properties of Li with experimental values make the study of electronic and transport properties in lithium nanowires interesting because they are promising candidates as bridging pieces in nanoelectronics. Transmission coefficient and V-I characteristic of Li nano wire indicates that Li nanowire can be used as an electrode device.
NASA Astrophysics Data System (ADS)
Nawa, Kenji; Nakamura, Kohji; Akiyama, Toru; Ito, Tomonori; Weinert, Michael
2015-03-01
Interest in single organometallic molecule and that adsorbed on solid surfaces has rapidly increased because of possible novel applications. For molecules with transition metals (TMs), the d-electron configuration is an essential aspect of their electronic and magnetic properties, and correlation effects can not be excluded. Here, we investigate systematically the electron configuration and correlation effects for prototypical organometallic molecules of tridimensional metallocene (TMCp2) and planer phthalocyanine (TMPc). Calculations were carried out based on the constraint density functional theory (DFT) by using the full-potential linearized augmented plane wave method that incorporates an on-site Coulomb interaction correction + U . We find that these correlation effects play a key role in determining the ground state of the electron configuration for the organometallic molecules. The calculated ground states of TMCp2, where TM =Cr, Mn, Fe, Co, and Ni, obtained by constraint DFT with +U reproduce the experimentally determined structures of 3E2 g , 6A1 g , 1A1 g , 2E1 g , and 3A2 g , respectively. Results for the TMPc will be also presented.
Vassilev, Peter; Louwerse, Manuel J; Baerends, Evert Jan
2005-12-15
Ab initio density functional theory molecular dynamics simulations of the solvated states of the hydroxyl radical and hydroxide ion are performed using the Becke-Lee-Yang-Parr (BLYP) exchange-correlation functional (Becke, A. D. Phys. Rev. A 1988, 38, 3098. Lee, C.; Yang, W.; Parr, R. G. Phys. Rev. B 1988, 37, 785). The structures of the solvation shells of the two species are examined. It is found that the OH radical forms a relatively well-defined solvation complex with four neighboring water molecules. Three of these molecules are hydrogen bonded to the OH, while the fourth is hemibonded via a three-electron two-centered bond between the oxygen atoms of the OH and water. The activity and the diffusion mechanism of the OH radical in water is discussed in comparison with the OH- ion. Although the results are partially influenced by the tendency of the BLYP density functional to overestimate hemibonded structure, the present simulations suggest that the widely accepted picture of rapid diffusion of OH radical in water through hydrogen exchange reaction may need to be reconsidered. PMID:16375337
NASA Astrophysics Data System (ADS)
Xiao, Bing; Sun, Jianwei; Ruzsinszky, Adrienn; Perdew, John P.
2014-08-01
We employ semilocal density functionals [local spin-density approximation (LSDA), Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and meta-GGAs)], LSDA plus Hubbard U (LSDA+U) theory, a nonlocal range-separated Heyd-Scuseria-Ernzerhof hybrid functional (HSE06), and the random-phase approximation (RPA) to assess their performances for the ground-state magnetism and electronic structure of a strongly correlated metal, rutile VO2. Using recent quantum Monte Carlo results as the benchmark, all tested semilocal and hybrid functionals as well as the RPA (with PBE inputs) predict the correct magnetic ground states for rutile VO2. The observed paramagnetism could arise from temperature-disordered local spin moments or from the thermal destruction of these moments. All semilocal functionals also give the correct ground-state metallicity for rutile VO2. However, in the ferromagnetic (FM) and antiferromagnetic (AFM) phases, LSDA+U and HSE06 incorrectly predict rutile VO2 to be a Mott-Hubbard insulator. For the computed electronic structures of FM and AFM phases, we find that the Tao-Perdew-Staroverov-Scuseria (TPSS) and revised TPSS (revTPSS) meta-GGAs give strong 2p-3d hybridizations, resulting in a depopulation of the 2p bands of O atoms, in comparison with other tested meta-GGAs. The regularized TPSS (regTPSS) and meta-GGAs made simple, i.e., MGGA_MS0 and MGGA_MS2, which are free of the spurious order-of-limits problem of TPSS and revTPSS, give electronic states close to those of the PBE GGA and LSDA. In comparison to experiment, semilocal functionals predict better equilibrium cell volumes for rutile VO2 in FM and AFM states than in the spin-unpolarized state. For meta-GGAs, a monotonic decrease of the exchange enhancement factor Fx(s,α) with α for small s, as in the MGGA_MS functionals, leads to large (probably too large) local magnetic moments in spin-polarized states.
NASA Astrophysics Data System (ADS)
Ma, Liang-Cai; Zhang, Jian-Min; Xu, Ke-Wei
2013-02-01
By using first-principles calculations based on density-functional theory, we have systematically investigated the equilibrium structure and electronic properties of 22 free-standing copper (Cu) nanowires having different cross-sections with 1-14 Cu atoms per unit cell. The structural properties of the studied Cu nanowires were greatly different from those of face centered cubic bulk Cu. For each wire the equilibrium lattice constant was obtained. The binding energy increases with increasing atom number per unit cell in different structures. As for the polygonal structures of a fixed cross-section, the preferred structures should be the staggered ones which contain a linear chain along the wire axis that passes through the center of the polygons, where each chain atom is just located at a point equidistant from the planes of polygons. All the nanowires are metallic. The numbers of conduction channels were determined from the electronic band structures. In general, the number of conduction channels increases when the nanowire becomes thicker. The density of charge revealed delocalized metallic bonding for all studied Cu nanowires.
NASA Astrophysics Data System (ADS)
Terrier, Cyril; Vitorge, Pierre; Gaigeot, Marie-Pierre; Spezia, Riccardo; Vuilleumier, Rodolphe
2010-07-01
Structural and electronic properties of La3+ immersed in bulk water have been assessed by means of density functional theory (DFT)-based Car-Parrinello molecular dynamics (CPMD) simulations. Correct structural properties, i.e., La(III)-water distances and La(III) coordination number, can be obtained within the framework of Car-Parrinello simulations providing that both the La pseudopotential and conditions of the dynamics (fictitious mass and time step) are carefully set up. DFT-MD explicitly treats electronic densities and is shown here to provide a theoretical justification to the necessity of including polarization when studying highly charged cations such as lanthanoids(III) with classical MD. La3+ was found to strongly polarize the water molecules located in the first shell, giving rise to dipole moments about 0.5 D larger than those of bulk water molecules. Finally, analyzing Kohn-Sham orbitals, we found La3+ empty 4f orbitals extremely compact and to a great extent uncoupled from the water conduction band, while the 5d empty orbitals exhibit mixing with unoccupied states of water.
NASA Astrophysics Data System (ADS)
Corsini, Niccolò R. C.; Greco, Andrea; Hine, Nicholas D. M.; Molteni, Carla; Haynes, Peter D.
2013-08-01
We present an implementation in a linear-scaling density-functional theory code of an electronic enthalpy method, which has been found to be natural and efficient for the ab initio calculation of finite systems under hydrostatic pressure. Based on a definition of the system volume as that enclosed within an electronic density isosurface [M. Cococcioni, F. Mauri, G. Ceder, and N. Marzari, Phys. Rev. Lett. 94, 145501 (2005)], 10.1103/PhysRevLett.94.145501, it supports both geometry optimizations and molecular dynamics simulations. We introduce an approach for calibrating the parameters defining the volume in the context of geometry optimizations and discuss their significance. Results in good agreement with simulations using explicit solvents are obtained, validating our approach. Size-dependent pressure-induced structural transformations and variations in the energy gap of hydrogenated silicon nanocrystals are investigated, including one comparable in size to recent experiments. A detailed analysis of the polyamorphic transformations reveals three types of amorphous structures and their persistence on depressurization is assessed.
Corsini, Niccolò R C; Greco, Andrea; Hine, Nicholas D M; Molteni, Carla; Haynes, Peter D
2013-08-28
We present an implementation in a linear-scaling density-functional theory code of an electronic enthalpy method, which has been found to be natural and efficient for the ab initio calculation of finite systems under hydrostatic pressure. Based on a definition of the system volume as that enclosed within an electronic density isosurface [M. Cococcioni, F. Mauri, G. Ceder, and N. Marzari, Phys. Rev. Lett. 94, 145501 (2005)], it supports both geometry optimizations and molecular dynamics simulations. We introduce an approach for calibrating the parameters defining the volume in the context of geometry optimizations and discuss their significance. Results in good agreement with simulations using explicit solvents are obtained, validating our approach. Size-dependent pressure-induced structural transformations and variations in the energy gap of hydrogenated silicon nanocrystals are investigated, including one comparable in size to recent experiments. A detailed analysis of the polyamorphic transformations reveals three types of amorphous structures and their persistence on depressurization is assessed. PMID:24006984
NASA Astrophysics Data System (ADS)
Gunasekaran, Sethu; Rajalakshmi, K.; Kumaresan, Subramanian
2013-08-01
The Fourier transform (FT-IR) spectrum of Levofloxacin was recorded in the region 4000-400 cm-1 and a complete vibrational assignment of fundamental vibrational modes of the molecule was carried out using density functional method. The observed fundamental modes have been compared with the harmonic vibrational frequencies computed using DFT (B3LYP) method by employing 6-31 G (d, p) basis sets. The most stable geometry of the molecule under investigation has been determined from the potential energy scan. The first-order hyperpolarizability (βo) and other related properties (μ, αo) of Levofloxacin are calculated using density functional theory (DFT) on a finite field approach. UV-vis spectrum of the molecule was recorded and the electronic properties, such as HOMO and LUMO energies were performed by DFT using 6-31 G (d, p) basis sets. Stability of the molecule arising from hyperconjugative interactions, charge delocalization have been analyzed using natural bond orbital analysis (NBO). The calculated HOMO and LUMO energies show that, the charge transfer occurs within the molecule. The other molecular properties like molecular electrostatic potential (MESP), Mulliken population analysis and thermodynamic properties of the title molecule have been calculated.
Gunasekaran, Sethu; Rajalakshmi, K; Kumaresan, Subramanian
2013-08-01
The Fourier transform (FT-IR) spectrum of Levofloxacin was recorded in the region 4000-400 cm(-1) and a complete vibrational assignment of fundamental vibrational modes of the molecule was carried out using density functional method. The observed fundamental modes have been compared with the harmonic vibrational frequencies computed using DFT (B3LYP) method by employing 6-31 G (d, p) basis sets. The most stable geometry of the molecule under investigation has been determined from the potential energy scan. The first-order hyperpolarizability (βo) and other related properties (μ, αo) of Levofloxacin are calculated using density functional theory (DFT) on a finite field approach. UV-vis spectrum of the molecule was recorded and the electronic properties, such as HOMO and LUMO energies were performed by DFT using 6-31 G (d, p) basis sets. Stability of the molecule arising from hyperconjugative interactions, charge delocalization have been analyzed using natural bond orbital analysis (NBO). The calculated HOMO and LUMO energies show that, the charge transfer occurs within the molecule. The other molecular properties like molecular electrostatic potential (MESP), Mulliken population analysis and thermodynamic properties of the title molecule have been calculated. PMID:23685802
Electronic structure and defect properties of Tl6SeI4: Density functional calculations
NASA Astrophysics Data System (ADS)
Biswas, Koushik; Du, Mao-Hua; Singh, David J.
2012-10-01
We report density functional calculations of electronic structure, phase diagram, and dielectric, optical, and defect properties of Tl6SeI4. We discuss how electronic structure and defect properties affect resistivity and carrier mobility-lifetime (μτ) products in Tl6SeI4. We find large Born effective charges due to covalency involving Tl-6p states. High Born charges generally enhance the static dielectric constant. This provides a mechanism for effective screening of charged defects and impurities. We find that high resistivity can be obtained under near-stoichiometric growth conditions via Fermi level pinning near the middle of the band gap by shallow donors and acceptors, as opposed to deep traps that can give high resistivity, but at the expense of short carrier drift lengths. Defect calculations also reveal the presence of deep native donors that may cause electron trapping. The experimentally observed good μτ products may be explained by a combination of small effective masses and effective screening of charged defects. High resistivity and good μτ products make Tl6SeI4 a promising room-temperature radiation detector material. We also show the calculated defect diffusion barriers, which affect defect migration under external bias in a detector.
Liu, Xia; Tan, Yingzi; Li, Xiuling; Wu, Xiaojun; Pei, Yong
2015-08-28
The electronic and magnetic properties of transition metal (TM = Sc, Ti, V, Cr and Mn) atom incorporated single and double one-dimensional (1D) styrene molecular wires confined on the hydrogen-terminated Si(100) surface are explored for the first time by means of spin-polarized density functional theory, denoted as Si-[TM(styrene)]. It is unveiled that TM atoms bind asymmetrically to the adjacent phenyl rings, which leads to novel electronic and magnetic properties in stark contrast to the well-studied gas phase TM-benzene molecular wires. Si-[Mn(styrene)]∞ and Si-[Cr(styrene)]∞ single molecular wires (SMWs) are a ferromagnetic semiconductor and half metal, respectively. Creation of H-atom defects on the silicon surface can introduce an impurity metallic band, which leads to novel half-metallic magnetism of a Si-[Mn(styrene)]∞ system. Moreover, double molecular wires (DMWs) containing two identical or hetero SMWs are theoretically designed. The [Mn(styrene)]∞-[Cr(styrene)]∞ DMW exhibits half-metallic magnetism where the spin-up and spin-down channels are contributed by two single molecular wires. Finally, we demonstrate that introducing a TM-defect may significantly affect the electronic structure and magnetic properties of molecular wires. These studies provide new insights into the structure and properties of surface supported 1-D sandwiched molecular wires and may inspire the future experimental synthesis of substrate confined organometallic sandwiched molecular wires. PMID:26219748
NASA Astrophysics Data System (ADS)
Endo, Kazunaka
2016-02-01
In the Auger electron spectra (AES) simulations, we define theoretical modified kinetic energies of AES in the density functional theory (DFT) calculations. The modified kinetic energies correspond to two final-state holes at the ground state and at the transition-state in DFT calculations, respectively. This method is applied to simulate Auger electron spectra (AES) of 2nd periodic atom (Li, Be, B, C, N, O, F)-involving substances (LiF, beryllium, boron, graphite, GaN, SiO2, PTFE) by deMon DFT calculations using the model molecules of the unit cell. Experimental KVV (valence band electrons can fill K-shell core holes or be emitted during KVV-type transitions) AES of the (Li, O) atoms in the substances agree considerably well with simulation of AES obtained with the maximum kinetic energies of the atoms, while, for AES of LiF, and PTFE substance, the experimental F KVV AES is almost in accordance with the spectra from the transitionstate kinetic energy calculations.
NASA Astrophysics Data System (ADS)
Ferri, Nicola; Distasio, Robert A., Jr.; Ambrosetti, Alberto; Car, Roberto; Scheffler, Matthias; Tkatchenko, Alexandre
2015-03-01
Ubiquitous long-range van der Waals (vdW) interactions play a fundamental role in the structure and stability of a wide range of systems. Within the DFT framework, the vdW energy represents a crucial, but tiny part of the total energy, hence its influence on the electronic density, n (r) , and electronic properties is typically assumed to be rather small. Here, we address this question via a fully self-consistent (SC) implementation of the interatomic Tkatchenko-Scheffler vdW functional and its extension to surfaces. Self-consistency leads to large changes in the binding energies and electrostatic moments of highly polarizable alkali metal dimers. For some metal surfaces, vdW interactions increase dipole moments and induce non-trivial charge rearrangements, leading to visible changes in the metal workfunctions. Similar behavior is observed for molecules adsorbed on metals. Our study reveals a non-trivial connection between electrostatics and long-range electron correlation effects.
Yuan, H. K.; Chen, H. Tian, C. L.; Kuang, A. L.; Wang, J. Z.
2014-04-21
Gadolinium-oxide clusters in various sizes and stoichiometries have been systematically studied by employing the density functional theory with the generalized gradient approximation. The clusters in bulk stoichiometry are relatively more stable and their binding energies increase with the increasing size. Stoichiometric (Gd{sub 2}O{sub 3}){sub n} clusters of n = 1–3 prefer cage-like structures, whereas the clusters of n = 4–30 prefer compact structures layered by wedge-like units and exhibit a rough feature toward the bulk-like arrangement with small disorders of atomic positions. The polyhedral-cages analogous to carbon-fullerenes are stable isomers yet not the minimum energy configurations. Their stabilities can be improved by embedding one oxygen atom or a suitable cage to form core-shell configurations. The mostly favored antiferromagnetic couplings between adjacent Gd atoms are nearly degenerated in energy with their ferromagnetic couplings, resulting in super-paramagnetic characters of gadolinium-oxide clusters. The Ruderman-Kittel-Kasuya-Yosida (RKKY)-type mechanism together with the superexchange-type mechanism plays cooperation role for the magnetic interactions in clusters. We present, as a function of n, calculated binding energies, ionization potential, electron affinity, and electronic dipole moment.
NASA Astrophysics Data System (ADS)
Yuan, H. K.; Chen, H.; Tian, C. L.; Kuang, A. L.; Wang, J. Z.
2014-04-01
Gadolinium-oxide clusters in various sizes and stoichiometries have been systematically studied by employing the density functional theory with the generalized gradient approximation. The clusters in bulk stoichiometry are relatively more stable and their binding energies increase with the increasing size. Stoichiometric (Gd2O3)n clusters of n = 1-3 prefer cage-like structures, whereas the clusters of n = 4-30 prefer compact structures layered by wedge-like units and exhibit a rough feature toward the bulk-like arrangement with small disorders of atomic positions. The polyhedral-cages analogous to carbon-fullerenes are stable isomers yet not the minimum energy configurations. Their stabilities can be improved by embedding one oxygen atom or a suitable cage to form core-shell configurations. The mostly favored antiferromagnetic couplings between adjacent Gd atoms are nearly degenerated in energy with their ferromagnetic couplings, resulting in super-paramagnetic characters of gadolinium-oxide clusters. The Ruderman-Kittel-Kasuya-Yosida (RKKY)-type mechanism together with the superexchange-type mechanism plays cooperation role for the magnetic interactions in clusters. We present, as a function of n, calculated binding energies, ionization potential, electron affinity, and electronic dipole moment.
NASA Astrophysics Data System (ADS)
Prabhaharan, M.; Prabakaran, A. R.; Srinivasan, S.; Gunasekaran, S.
2015-03-01
The present work has been carried out a combined experimental and theoretical study on molecular structure, vibrational spectra and NBO analysis of cyanuric acid. The FT-IR (100-4000 cm-1) and FT-Raman spectra (400-4000 cm-1) of cyanuric acid were recorded. In DFT methods, Becke's three parameter exchange-functional (B3) combined with gradient-corrected correlation functional of Lee, Yang and Parr (LYP) by implementing the split-valence polarized 6-31G(d,p) and 6-31++G(d,p) basis sets have been considered for the computation of the molecular structure optimization, vibrational frequencies, thermodynamic properties and energies of the optimized structures. The density functional theory (DFT) result complements the experimental findings. The electronic properties, such as HOMO-LUMO energies and molecular electrostatic potential (MESP) are also performed. Mulliken population analysis on atomic charges is also calculated. The first order hyperpolarizability (βtotal) of this molecular system and related properties (β, μ and Δα) are calculated using DFT/B3LYP/6-31G (d,p) and B3LYP/6-311++G(d,p) methods. The thermodynamic functions (heat capacity, entropy and enthalpy) from spectroscopic data by statistical methods were also obtained for the range of temperature 50-1000 K.
NASA Astrophysics Data System (ADS)
Romero, Nichols Anthony
Two decades after the discovery of C60, fullerenes continue to be intensely studied for their diverse properties, not the least of which is their unusually high superconducting transition temperature T c (up to 40 K in Cs3C60). Ab initio electronic structure methods are able to contribute to our understanding of these materials by providing electron densities, band structures, density-of-states, binding energies, and even the electron-phonon coupling. Because of the increasing computational feasibility of these large-scale calculations, these methods now play a prominent role in verifying and guiding experimental investigation of new materials. This dissertation presents results of a theoretical investigation of several scenarios where Tc enhancement exceeding those found in the alkali-doped fullerides could be exhibited: (1) field-effect doping of C60 layers, (2) C28-derived molecular solids and (3) covalently bonded C28 solids. The method employed is the Kohn-Sham formulation of density functional theory. Simpler tight-binding calculations are also used when appropriate. The study of field-effect doping was stimulated by the reports of large Tc enhancements by Schon et al., which were later retracted and found to be falsified. Even before the legitimacy of these reports came into question, we concluded that our calculations did not substantiate their claims. We present our electronic structure calculations and conclusions which are independent and potentially useful for future work on field-effect devices. The main part of this dissertation is a separate study on C28-based solids motivated by theoretical arguments suggesting that solids based on smaller fullerenes could exhibit a Tc enhancement. Among many molecular solids formed from closed-shell C28-derived molecules, we have found that solid C28H4 binds weakly and exhibit many of the salient features of solid C60, with an estimated Tc of 58 K. In this same spirit, we also study covalent solids based on
NASA Astrophysics Data System (ADS)
Di Valentin, M.; Salvadori, E.; Barone, V.; Carbonera, D.
2013-10-01
Advanced electron paramagnetic resonance (EPR) techniques, in combination with Density Functional theory (DFT), have been applied to the comparative study of carotenoid triplet states in two major photosynthetic antenna complexes, the Peridinin-chlorophyll a-protein of dinoflagellates and the light-harvesting complex II of higher plants. Carotenoid triplet states are populated by triplet-triplet energy transfer (TTET) from chlorophyll molecules to photoprotect the system from singlet oxygen formation under light-stress conditions. The TTET process is strongly dependent on the relative arrangement and on the electronic properties of the triplet states involved. The proposed spectroscopic approach exploits the concept of spin conservation during TTET, which leads to recognisable spin polarisation effects in the time-resolved and field-swept echo-detected EPR spectra. The electron spin polarisation produced at the carotenoid acceptor site depends on the initial polarisation of the chlorophyll donor and on the relative geometrical arrangement of the donor-acceptor zero-field splitting axes. We have demonstrated that a proper analysis of the spectra in the framework of spin angular momentum conservation allows to derive the pathways of TTET and to gain insight into the structural requirements of this mechanism for those antenna complexes, whose X-ray structure is available. We have further proved that this method, developed for natural antenna complexes of known X-ray structure, can be extended to systems lacking structural information in order to derive the relative arrangement of the partners in the energy transfer process. The structural requirements for efficient TTET, obtained from time-resolved and pulse EPR, have been complemented by a detailed description of the electronic structure of the carotenoid triplet state, provided by pulse Electron-Nuclear DOuble Resonance (ENDOR) experiments. Triplet-state hyperfine couplings of the α- and β-protons of the
Drew, Simon C; Young, Charles G; Hanson, Graeme R
2007-04-01
The electron paramagnetic resonance spin Hamiltonian parameters of mononuclear thiomolybdenyl complexes based upon the tris(pyrazolyl)borate ligand, together with their molybdenyl analogues, are calculated using density functional theory. The electronic g matrix and 95Mo hyperfine matrix are calculated as second-order response properties from the coupled-perturbed Kohn-Sham equations. The scalar relativistic zero-order regular approximation (ZORA) is used with an all-electron basis and an accurate mean-field spin-orbit operator which includes all one- and two-electron terms. The principal values and relative orientations of the g and A interaction matrices obtained from the experimental spectra in a previous EPR study are compared with those obtained from unrestricted Kohn-Sham calculations at the BP86 and B3LYP level, and the latter are found to be in good quantitative agreement. A quasi-restricted approach is used to analyze the influence of the various molecular orbitals on g and A. In all complexes the ground state magnetic orbital is dX2-Y2-based and the orientation of the A matrix is directly related to the orientation of this orbital. The largest single contribution to the orientation of the g matrix arises from the spin-orbit coupling of the dYZ-based lowest-unoccupied molecular orbital into the ground state. A number of smaller, cumulative charge-transfer contributions augment the d-d contributions. A comparison of the theoretical EPR parameters obtained using both crystallographic and gas-phase geometry-optimized structures of Tp*MoO(bdt) (Tp* = hydrotris(3,5-dimethylpyrazol-1-yl)borate, bdt = 1,2-benzenedithiolate) suggests a correspondence between the metal-dithiolate fold angle and the angle of noncoincidence between g and A. PMID:17305330
NASA Astrophysics Data System (ADS)
Kaewmaraya, T.; Ramzan, M.; Löfâs, H.; Ahuja, Rajeev
2013-01-01
In this article, we use hybrid density functional (HSE06) to study the crystal and electronic structures and optical properties of well known phase change memory material Ge2Sb2Te5. We calculate the structural parameters, band gaps, and dielectric functions of three stable structures of this material. We also analyze the electron charge distribution using the Bader's theory of charge analysis. We find that hybrid density functional slightly overestimates the value of "c" parameter. However, overall, our results calculated with the use of hybrid density functional (HSE06) are very close to available experimental values than calculated with the use of Perdew Burke-Ernzerhof functional. Specifically, the electronic band gap values of this material calculated with HSE06 are in good agreement with the available experimental data in the literature. Furthermore, we perform the charge analysis and find that naive ionic model fails to explain the charge distribution between the constituent atoms, showing the complex nature of this compound.
NASA Astrophysics Data System (ADS)
Liu, Xia; Tan, Yingzi; Li, Xiuling; Wu, Xiaojun; Pei, Yong
2015-08-01
The electronic and magnetic properties of transition metal (TM = Sc, Ti, V, Cr and Mn) atom incorporated single and double one-dimensional (1D) styrene molecular wires confined on the hydrogen-terminated Si(100) surface are explored for the first time by means of spin-polarized density functional theory, denoted as Si-[TM(styrene)]. It is unveiled that TM atoms bind asymmetrically to the adjacent phenyl rings, which leads to novel electronic and magnetic properties in stark contrast to the well-studied gas phase TM-benzene molecular wires. Si-[Mn(styrene)]∞ and Si-[Cr(styrene)]∞ single molecular wires (SMWs) are a ferromagnetic semiconductor and half metal, respectively. Creation of H-atom defects on the silicon surface can introduce an impurity metallic band, which leads to novel half-metallic magnetism of a Si-[Mn(styrene)]∞ system. Moreover, double molecular wires (DMWs) containing two identical or hetero SMWs are theoretically designed. The [Mn(styrene)]∞-[Cr(styrene)]∞ DMW exhibits half-metallic magnetism where the spin-up and spin-down channels are contributed by two single molecular wires. Finally, we demonstrate that introducing a TM-defect may significantly affect the electronic structure and magnetic properties of molecular wires. These studies provide new insights into the structure and properties of surface supported 1-D sandwiched molecular wires and may inspire the future experimental synthesis of substrate confined organometallic sandwiched molecular wires.The electronic and magnetic properties of transition metal (TM = Sc, Ti, V, Cr and Mn) atom incorporated single and double one-dimensional (1D) styrene molecular wires confined on the hydrogen-terminated Si(100) surface are explored for the first time by means of spin-polarized density functional theory, denoted as Si-[TM(styrene)]. It is unveiled that TM atoms bind asymmetrically to the adjacent phenyl rings, which leads to novel electronic and magnetic properties in stark contrast to
Uzunova, Ellie L; Mikosch, Hans
2012-03-29
The dimers of cobalt oxide (CoO)(2) with cyclic and open bent structure are studied with the B1LYP density functional; the ordering of states is validated by the CCSD(T) method. The D(2h)-symmetry rhombic dioxide Co(2)O(2) with antiferromagnetically ordered electrons on cobalt centers is the global minimum. The cyclic peroxide Co(2)(O(2)) with side-on-bonded dioxygen in (7)B(2) ground state is separated from the global minimum by an energy gap of 3.15 eV. The dioxide is highly reactive as indicated by the high value of proton affinity and chemical reactivity indices. The four-member ring structures are more stable than those with three-member ring or chain configuration. The thermodynamic stability toward dissociation to CoO increases upon carbonylation, whereas proton affinity and reactivity with release of molecular oxygen also increase. The global minimum of Co(2)O(2)(CO)(6) corresponds to a triplet state (3)A" with oxygen atoms shifted above the molecular plane of the rhombic dioxide Co(2)O(2). The SOMO-LUMO gap in the ground-state carbonylated dioxide is wider, compared to the same gap in the bare dicobalt dioxide. The peroxo-isomer Co(2)(O(2))(CO)(6) retains the planar Co(2)(O(2)) ring and is only stable in a high-spin state (7)A". The carbonylated clusters have increased reactivity in both redox and nucleophilic reactions, as a result of the increased electron density in the Co(2)O(2)-ring area. PMID:22397598
Domain boundaries in silicene: Density functional theory calculations on electronic properties
NASA Astrophysics Data System (ADS)
Xiao, Hong-Jun; Zhang, Li-Zhi; Du, Shi-Xuan; Gao, Hong-Jun
2015-08-01
By using density functional theory (DFT)-based first-principles calculations, the structural stability and electronic properties for two kinds of silicene domain boundaries, forming along armchair edge and zigzag edge, have been investigated. The results indicate that a linkage of tetragonal and octagonal rings (4|8) appears along the armchair edge, while a linkage of paired pentagonal and octagonal rings (5|5|8) appears along the zigzag edge. Different from graphene, the buckling properties of silicene lead to two mirror symmetrical edges of silicene line-defect. The formation energies indicate that the 5|5|8 domain boundary is more stable than the 4|8 domain boundary. Similar to graphene, the calculated electronic properties show that the 5|5|8 domain boundaries exhibit metallic properties and the 4|8 domain boundaries are half-metal. Both domain boundaries create the perfect one-dimensional (1D) metallic wires. Due to the metallic properties, these two kinds of nanowires can be used to build the silicene-based devices. Project supported by the National Natural Science Foundation of China (Grant Nos. 61390501 and 51325204), the National Basic Research Program of China (Grant Nos. 2011CB808401 and 2011CB921702), and the Tainjin Supercomputing Center, Chinese Academy of Sciences.
NASA Astrophysics Data System (ADS)
Andrzejak, Marcin; Sterzel, Mariusz; Pawlikowski, Marek T.
2005-07-01
The absorption spectra of the N-(2,5-di- tert-butylphenyl) phthalimide ( 1-), N-(2,5-di- tert-butylphenyl)-1,8-naphthalimide ( 2-) and N-(2,5-di- tert-butylphenyl)-perylene-3,4-dicarboximide ( 3-) anion radicals are studied in terms of time dependent density functional theory (TDDFT). For these anion radicals a large number electronic states (from 30 to 60) was found in the visible and near-IR regions (5000-45000 cm -1). In these regions the TD/B3LYP treatment at the 6-1+G* level is shown to reproduce satisfactorily the empirical absorption spectra of all three anion radicals studied. The most apparent discrepancies between purely electronic theory and the experiment could be found in the excitation region corresponding to D0→ D1 transitions in the 2- and 3- molecules. For these species we argue that the structures seen in the lowest energy part of the absorptions of the 2- and 3- species are very likely due to Franck-Condon (FC) activity of the totally symmetric vibrations not studied in this Letter.
NASA Astrophysics Data System (ADS)
Wang, Jun-Fei; Fu, Xiao-Nan; Zhang, Xiao-Dong; Wang, Jun-Tao; Li, Xiao-Dong; Jiang, Zhen-Yi
2016-08-01
The structural, elastic, electronic, and thermodynamic properties of thermoelectric material MgAgSb in γ,β,α phases are studied with first-principles calculations based on density functional theory. The optimized lattice constants accord well with the experimental data. According to the calculated total energy of the three phases, the phase transition order is determined from α to γ phase with cooling, which is in agreement with the experimental result. The physical properties such as elastic constants, bulk modulus, shear modulus, Young’s modulus, Poisson’s ratio, and anisotropy factor are also discussed and analyzed, which indicates that the three structures are mechanically stable and each has a ductile feature. The Debye temperature is deduced from the elastic properties. The total density of states (TDOS) and partial density of states (PDOS) of the three phases are investigated. The TDOS results show that the γ phase is most stable with a pseudogap near the Fermi level, and the PDOS analysis indicates that the conduction band of the three phases is composed mostly of Mg-3s, Ag-4d, and Sb-5p. In addition, the changes of the free energy, entropy, specific heat, thermal expansion of γ-MgAgSb with temperature are obtained successfully. The obtained results above are important parameters for further experimental and theoretical tuning of doped MgAgSb as a thermoelectric material at high temperature. Project supported by the National Natural Science Foundation of China (Grant No. 11504088), the Fund from Henan University of Technology, China (Grant Nos. 2014YWQN08 and 2013JCYJ12), the Natural Science Fund from the Henan Provincial Education Department, China (Grant No. 16A140027), the Natural Science Foundation of Shaanxi Province of China (Grant Nos. 2013JQ1018 and 15JK1759), and the Science Foundation of Northwest University of China (Grant No. 14NW23).
Density-functional theory study of gramicidin A ion channel geometry and electronic properties
Todorović, Milica; Bowler, David R.; Gillan, Michael J.; Miyazaki, Tsuyoshi
2013-01-01
Understanding the mechanisms underlying ion channel function from the atomic-scale requires accurate ab initio modelling as well as careful experiments. Here, we present a density functional theory (DFT) study of the ion channel gramicidin A (gA), whose inner pore conducts only monovalent cations and whose conductance has been shown to depend on the side chains of the amino acids in the channel. We investigate the ground state geometry and electronic properties of the channel in vacuum, focusing on their dependence on the side chains of the amino acids. We find that the side chains affect the ground state geometry, while the electrostatic potential of the pore is independent of the side chains. This study is also in preparation for a full, linear scaling DFT study of gA in a lipid bilayer with surrounding water. We demonstrate that linear scaling DFT methods can accurately model the system with reasonable computational cost. Linear scaling DFT allows ab initio calculations with 10 000–100 000 atoms and beyond, and will be an important new tool for biomolecular simulations. PMID:24068174
Nagy, A; Amovilli, C
2008-03-21
In the ground state, the pair density n can be determined by solving a single auxiliary equation of a two-particle problem. Electron-electron cusp condition and asymptotic behavior for the Pauli potential of the effective potential of the two-particle equation are presented. PMID:18361562
Density-functional study of atomic and electronic structures of multivacancies in silicon carbide
NASA Astrophysics Data System (ADS)
Iwata, Jun-Ichi; Shinei, Chikara; Oshiyama, Atsushi
2016-03-01
We report the density-functional calculations that provide a firm theoretical framework to identify the multivacancies and unravel the underlying physics in the most stable silicon carbide polytype 4H-SiC. The calculations with the generalized gradient approximation (GGA) for the Si and C monovacancy, VSi and VC, have clarified the significantly lower formation energy of the C vacancy accompanied by the efficient pairing relaxation of the surrounding Si dangling bonds. Our GGA calculation also predicts a stable next-neighbor VCVC divacancy which is lower in the formation energy than the usual nearest-neighbor divacancy VCVSi discussed in the past. We also perform the calculations with the hybrid functional and confirm the stability of the next-neighbor VCVC divacancy. Our calculations indeed clarify that it is possible to detect the VCVC with its peculiar hyperfine coupling constants by the electron paramagnetic resonance (EPR) measurements. Based on the structural characteristics and the energetics for the monovacancy, we further propose an extended dangling-bond-counting (EDBC) model to choose the energetically favorable topological network of the vacant sites for the multivacancy. The GGA calculations combined with the EDBC model reveal that V3 and V6 are energetically favorable. The stable V3 is a nearest-neighbor complex of VC-VSi-VC , whereas the V6 is the high-symmetry V5 (the central VSi surrounded by four VC) plus a symmetry-breaking next-neighbor VC. We perform the GGA calculations for the electronic structure of such V3 and V6 and discuss the possibility of detecting these multivacancies. In particular, the EPR-detected ANN1 center is provisionally identified as the doubly positive trivacancy VC-VSi-VC .
Assessment of the ΔSCF density functional theory approach for electronic excitations in organic dyes
Kowalczyk, T.; Yost, S. R.; Van Voorhis, T.
2010-01-01
This paper assesses the accuracy of the ΔSCF method for computing low-lying HOMO→LUMO transitions in organic dye molecules. For a test set of vertical excitation energies of 16 chromophores, surprisingly similar accuracy is observed for time-dependent density functional theory and for ΔSCF density functional theory. In light of this performance, we reconsider the ad hoc ΔSCF prescription and demonstrate that it formally obtains the exact stationary density within the adiabatic approximation, partially justifying its use. The relative merits and future prospects of ΔSCF for simulating individual excited states are discussed.
NASA Astrophysics Data System (ADS)
Chen, Hai-Hua; Bi, Yan; Cheng, Yan; Ji, Guangfu; Cai, Lingcang
2012-10-01
The elastic properties, electronic structure and thermodynamic behavior of the TaB have been investigated for the first time in this work. Using first-principles plane-wave ultrasoft-pseudopotential density functional theory (DFT), the ground state properties and equation of state of TaB have been obtained. The average zero-pressure bulk modulus of TaB is 302 GPa. By analyzing the elastically anisotropic behavior and the relative structure parameters of TaB, we found that the crystal cell along the b-axis was more compressible than along the a and c axes. The calculated ratio of bulk modulus and shear modulus (B/G) for TaB is 1.58, demonstrating that TaB is rather brittle. From the elastic stiffness constants, we found that TaB in the Cmcm phase is mechanically stable. The calculated hardness of TaB is 28.6 GPa which is close to the previous data. Moreover, using the Gibbs 2 model, the thermodynamic properties such as the thermal expansion and Debye temperature of TaB have been obtained firstly. At the ambient temperature, the Debye temperatures of TaB are 792 K and 845 K from GGA calculation and LDA calculation, respectively.
Superheavy Element Chemistry by Relativistic Density Functional Theory Electronic Structure Modeling
NASA Astrophysics Data System (ADS)
Zaitsevskii, A. V.; Polyaev, A. V.; Demidov, Yu. A.; Mosyagin, N. S.; Lomachuk, Yu. V.; Titov, A. V.
2015-06-01
Two-component density functional theory in its non-collinear formulation combined with the accurate relativistic electronic structure model defined by shape-consistent small-core pseudopotentials (PP/RDFT) provides a robust basis of efficient computational schemes for predicting energetic and structural properties of complex polyatomic systems including superheavy elements (SHEs). Because of the exceptional role of thermochromatography in the experiments on the "chemical" identification of SHEs with atomic numbers Z ≥ 112, we focus on the description of the adsorption of single SHE atoms on the surfaces of solids through cluster modeling of adsorption complexes. In some cases our results differ significantly from those of previous theoretical studies. The results of systematic comparative studies on chemical bonding in simple molecules of binary compounds of SHEs and their nearest homologs with most common light elements, obtained at the PP/RDFT level and visualized through the "chemical graphs", provide the understanding of the general chemistry of SHEs which at present cannot be derived from the experimental data. These results are used to discuss the main trends in changing chemical properties of the elements in the given group of the periodic table and demonstrate the specificity of SHEs.
NASA Astrophysics Data System (ADS)
Turel, N.; Arikan, F.
2010-12-01
Ionospheric channel characterization is an important task for both HF and satellite communications. The inherent space-time variability of the ionosphere can be observed through total electron content (TEC) that can be obtained using GPS receivers. In this study, within-the-hour variability of the ionosphere over high-latitude, midlatitude, and equatorial regions is investigated by estimating a parametric model for the probability density function (PDF) of GPS-TEC. PDF is a useful tool in defining the statistical structure of communication channels. For this study, a half solar cycle data is collected for 18 GPS stations. Histograms of TEC, corresponding to experimental probability distributions, are used to estimate the parameters of five different PDFs. The best fitting distribution to the TEC data is obtained using the maximum likelihood ratio of the estimated parametric distributions. It is observed that all of the midlatitude stations and most of the high-latitude and equatorial stations are distributed as lognormal. A representative distribution can easily be obtained for stations that are located in midlatitude using solar zenith normalization. The stations located in very high latitudes or in equatorial regions cannot be described using only one PDF distribution. Due to significant seasonal variability, different distributions are required for summer and winter.
NASA Astrophysics Data System (ADS)
Abdulsattar, Mudar Ahmed
2016-04-01
Energetically stable Ge diamondoids are used as building blocks to investigate Ge nanocrystal properties using density functional theory (DFT). Energy gap, HOMO and LUMO of Ge diamondoids vary according to confinement theory with shape fluctuations. Ge diamondoid vibrational force constants reach 2.18 mDyne/Å which is between that of bulk silicon and tin. Ge-Ge vibrational frequencies and reduced masses reach 357.47 cm‑1 and 41.93 amu, respectively, which are higher than the values of bulk Ge. Size variation of UV-Vis shows that the maximum optical peak moves from 163 nm to nearly 290 nm as the size of Ge diamondoids and molecules increases. The higher peak value approaches that of experimental Ge quantum dots at 300 nm. NMR spectra of Ge diamondoids are analyzed as a function of diamondoid sizes. Our results show that all investigated spectroscopic tools are sensitive to molecular or nanocrystals size. NMR is particularly good size indicator. Natural bond orbital (NBO) population analysis shows present diamondoid bondings differ from ideal sp3 bonding. The bonding for Ge valence orbitals is in the range (4s1.244p2.60) ‑ (4s1.274p2.73) depending on distance between Ge atom and diamondoid surface. Highest Ge diamondoids vibrational longitudinal optical (LO) mode is blue shifted with respect to experimental bulk value which is the opposite case for C and Si. H surface atom effects on electronic and vibrational properties are discussed.
Gidopoulos, Nikitas I.; Gross, E. K. U.
2014-01-01
A novel treatment of non-adiabatic couplings is proposed. The derivation is based on a theorem by Hunter stating that the wave function of the complete system of electrons and nuclei can be written, without approximation, as a Born–Oppenheimer (BO)-type product of a nuclear wave function, X(R), and an electronic one, ΦR(r), which depends parametrically on the nuclear configuration R. From the variational principle, we deduce formally exact equations for ΦR(r) and X(R). The algebraic structure of the exact nuclear equation coincides with the corresponding one in the adiabatic approximation. The electronic equation, however, contains terms not appearing in the adiabatic case, which couple the electronic and the nuclear wave functions and account for the electron–nuclear correlation beyond the BO level. It is proposed that these terms can be incorporated using an optimized local effective potential. PMID:24516183
Partition density functional theory
NASA Astrophysics Data System (ADS)
Nafziger, Jonathan
Partition density functional theory (PDFT) is a method for dividing a molecular electronic structure calculation into fragment calculations. The molecular density and energy corresponding to Kohn Sham density-functional theory (KS-DFT) may be exactly recovered from these fragments. Each fragment acts as an isolated system except for the influence of a global one-body 'partition' potential which deforms the fragment densities. In this work, the developments of PDFT are put into the context of other fragment-based density functional methods. We developed three numerical implementations of PDFT: One within the NWChem computational chemistry package using basis sets, and the other two developed from scratch using real-space grids. It is shown that all three of these programs can exactly reproduce a KS-DFT calculation via fragment calculations. The first of our in-house codes handles non-interacting electrons in arbitrary one-dimensional potentials with any number of fragments. This code is used to explore how the exact partition potential changes for different partitionings of the same system and also to study features which determine which systems yield non-integer PDFT occupations and which systems are locked into integer PDFT occupations. The second in-house code, CADMium, performs real-space calculations of diatomic molecules. Features of the exact partition potential are studied for a variety of cases and an analytical formula determining singularities in the partition potential is derived. We introduce an approximation for the non-additive kinetic energy and show how this quantity can be computed exactly. Finally a PDFT functional is developed to address the issues of static correlation and delocalization errors in approximations within DFT. The functional is applied to the dissociation of H2 + and H2.
Visualization of electronic density
Grosso, Bastien; Cooper, Valentino R.; Pine, Polina; Hashibon, Adham; Yaish, Yuval; Adler, Joan
2015-04-22
An atom’s volume depends on its electronic density. Although this density can only be evaluated exactly for hydrogen-like atoms, there are many excellent numerical algorithms and packages to calculate it for other materials. 3D visualization of charge density is challenging, especially when several molecular/atomic levels are intertwined in space. We explore several approaches to 3D charge density visualization, including the extension of an anaglyphic stereo visualization application based on the AViz package to larger structures such as nanotubes. We will describe motivations and potential applications of these tools for answering interesting questions about nanotube properties.
Lykissa, Iliana; Li, Shu-Yi; Granqvist, Claes G.; Niklasson, Gunnar A.; Ramzan, Muhammad; Chakraborty, Sudip; Ahuja, Rajeev
2014-05-14
Thin films of V{sub 2}O{sub 5} were prepared by sputter deposition onto transparent and electrically conducting substrates and were found to be X-ray amorphous. Their electrochemical density of states was determined by chronopotentiometry and displayed a pronounced low-energy peak followed by an almost featureless contribution at higher energies. These results were compared with density functional theory calculations for amorphous V{sub 2}O{sub 5}. Significant similarities were found between measured data and computations; specifically, the experimental low-energy peak corresponds to a split-off part of the conduction band apparent in the computations. Furthermore, the calculations approximately reproduce the experimental band gap observed in optical measurements.
Structural and electronic properties of poly(vinyl alcohol) using density functional theory
Dabhi, Shweta Jha, Prafulla K.
2014-04-24
The first principles calculations have been carried out to investigate the structural, electronic band structure density of states along with the projected density of states for poly(vinyl alcohol). Our structural calculation suggests that the poly(vinyl alcohol) exhibits monoclinic structure. The calculated structural lattice parameters are in excellent agreement with available experimental values. The band structure calculations reveal that the direct and indirect band gaps are 5.55 eV and 5.363 eV respectively in accordance with experimental values.
Oshiyama, Atsushi Iwata, Jun-Ichi; Uchida, Kazuyuki; Matsushita, Yu-Ichiro
2015-03-21
We show that our real-space finite-difference scheme allows us to perform density-functional calculations for nanometer-scale targets containing more than 100 000 atoms. This real-space scheme is applied to twisted bilayer graphene, clarifying that Moiré pattern induced in the slightly twisted bilayer graphene drastically modifies the atomic and electronic structures.
Svane, A.; Trygg, J.; Johansson, B.; Eriksson, O. |
1997-09-01
Electronic-structure calculations of elemental praseodymium are presented. Several approximations are used to describe the Pr f electrons. It is found that the low-pressure, trivalent phase is well described using either the self-interaction corrected (SIC) local-spin-density (LSD) approximation or the generalized-gradient approximation (GGA) with spin and orbital polarization (OP). In the SIC-LSD approach the Pr f electrons are treated explicitly as localized with a localization energy given by the self-interaction of the f orbital. In the GGA+OP scheme the f-electron localization is described by the onset of spin and orbital polarization, the energetics of which is described by spin-moment formation energy and a term proportional to the total orbital moment, L{sub z}{sup 2}. The high-pressure phase is well described with the f electrons treated as band electrons, in either the LSD or the GGA approximations, of which the latter describes more accurately the experimental equation of state. The calculated pressure of the transition from localized to delocalized behavior is 280 kbar in the SIC-LSD approximation and 156 kbar in the GGA+OP approach, both comparing favorably with the experimentally observed transition pressure of 210 kbar. {copyright} {ital 1997} {ital The American Physical Society}
Visualization of electronic density
NASA Astrophysics Data System (ADS)
Grosso, Bastien; Cooper, Valentino R.; Pine, Polina; Hashibon, Adham; Yaish, Yuval; Adler, Joan
2015-10-01
The spatial volume occupied by an atom depends on its electronic density. Although this density can only be evaluated exactly for hydrogen-like atoms, there are many excellent algorithms and packages to calculate it numerically for other materials. Three-dimensional visualization of charge density is challenging, especially when several molecular/atomic levels are intertwined in space. In this paper, we explore several approaches to this, including the extension of an anaglyphic stereo visualization application based on the AViz package for hydrogen atoms and simple molecules to larger structures such as nanotubes. We will describe motivations and potential applications of these tools for answering interesting physical questions about nanotube properties.
Mentel, Ł M; van Meer, R; Gritsenko, O V; Baerends, E J
2014-06-01
For chemistry an accurate description of bond weakening and breaking is vital. The great advantage of density matrix functionals, as opposed to density functionals, is their ability to describe such processes since they naturally cover both nondynamical and dynamical correlation. This is obvious in the Löwdin-Shull functional, the exact natural orbital functional for two-electron systems. We present in this paper extensions of this functional for the breaking of a single electron pair bond in N-electron molecules, using LiH, BeH(+), and Li2 molecules as prototypes. Attention is given to the proper formulation of the functional in terms of not just J and K integrals but also the two-electron L integrals (K integrals with a different distribution of the complex conjugation of the orbitals), which is crucial for the calculation of response functions. Accurate energy curves are obtained with extended Löwdin-Shull functionals along the complete dissociation coordinate using full CI calculations as benchmark. PMID:24907988
NASA Astrophysics Data System (ADS)
Mentel, Ł. M.; van Meer, R.; Gritsenko, O. V.; Baerends, E. J.
2014-06-01
For chemistry an accurate description of bond weakening and breaking is vital. The great advantage of density matrix functionals, as opposed to density functionals, is their ability to describe such processes since they naturally cover both nondynamical and dynamical correlation. This is obvious in the Löwdin-Shull functional, the exact natural orbital functional for two-electron systems. We present in this paper extensions of this functional for the breaking of a single electron pair bond in N-electron molecules, using LiH, BeH+, and Li2 molecules as prototypes. Attention is given to the proper formulation of the functional in terms of not just J and K integrals but also the two-electron L integrals (K integrals with a different distribution of the complex conjugation of the orbitals), which is crucial for the calculation of response functions. Accurate energy curves are obtained with extended Löwdin-Shull functionals along the complete dissociation coordinate using full CI calculations as benchmark.
Mentel, Ł. M.; Meer, R. van; Gritsenko, O. V.; Baerends, E. J.
2014-06-07
For chemistry an accurate description of bond weakening and breaking is vital. The great advantage of density matrix functionals, as opposed to density functionals, is their ability to describe such processes since they naturally cover both nondynamical and dynamical correlation. This is obvious in the Löwdin-Shull functional, the exact natural orbital functional for two-electron systems. We present in this paper extensions of this functional for the breaking of a single electron pair bond in N-electron molecules, using LiH, BeH{sup +}, and Li{sub 2} molecules as prototypes. Attention is given to the proper formulation of the functional in terms of not just J and K integrals but also the two-electron L integrals (K integrals with a different distribution of the complex conjugation of the orbitals), which is crucial for the calculation of response functions. Accurate energy curves are obtained with extended Löwdin-Shull functionals along the complete dissociation coordinate using full CI calculations as benchmark.
van Meer, R; Gritsenko, O V; Baerends, E J
2014-01-14
Time dependent density matrix functional theory in its adiabatic linear response formulation delivers exact excitation energies ωα and oscillator strengths fα for two-electron systems if extended to the so-called phase including natural orbital (PINO) theory. The Löwdin-Shull expression for the energy of two-electron systems in terms of the natural orbitals and their phases affords in this case an exact phase-including natural orbital functional (PILS), which is non-primitive (contains other than just J and K integrals). In this paper, the extension of the PILS functional to N-electron systems is investigated. With the example of an elementary primitive NO functional (BBC1) it is shown that current density matrix functional theory ground state functionals, which were designed to produce decent approximations to the total energy, fail to deliver a qualitatively correct structure of the (inverse) response function, due to essential deficiencies in the reconstruction of the two-body reduced density matrix (2RDM). We now deduce essential features of an N-electron functional from a wavefunction Ansatz: The extension of the two-electron Löwdin-Shull wavefunction to the N-electron case informs about the phase information. In this paper, applications of this extended Löwdin-Shull (ELS) functional are considered for the simplest case, ELS(1): one (dissociating) two-electron bond in the field of occupied (including core) orbitals. ELS(1) produces high quality ωα(R) curves along the bond dissociation coordinate R for the molecules LiH, Li2, and BH with the two outer valence electrons correlated. All of these results indicate that response properties are much more sensitive to deficiencies in the reconstruction of the 2RDM than the ground state energy, since derivatives of the functional with respect to both the NOs and the occupation numbers need to be accurate. PMID:24437859
NASA Astrophysics Data System (ADS)
van Meer, R.; Gritsenko, O. V.; Baerends, E. J.
2014-01-01
Time dependent density matrix functional theory in its adiabatic linear response formulation delivers exact excitation energies ωα and oscillator strengths fα for two-electron systems if extended to the so-called phase including natural orbital (PINO) theory. The Löwdin-Shull expression for the energy of two-electron systems in terms of the natural orbitals and their phases affords in this case an exact phase-including natural orbital functional (PILS), which is non-primitive (contains other than just J and K integrals). In this paper, the extension of the PILS functional to N-electron systems is investigated. With the example of an elementary primitive NO functional (BBC1) it is shown that current density matrix functional theory ground state functionals, which were designed to produce decent approximations to the total energy, fail to deliver a qualitatively correct structure of the (inverse) response function, due to essential deficiencies in the reconstruction of the two-body reduced density matrix (2RDM). We now deduce essential features of an N-electron functional from a wavefunction Ansatz: The extension of the two-electron Löwdin-Shull wavefunction to the N-electron case informs about the phase information. In this paper, applications of this extended Löwdin-Shull (ELS) functional are considered for the simplest case, ELS(1): one (dissociating) two-electron bond in the field of occupied (including core) orbitals. ELS(1) produces high quality ωα(R) curves along the bond dissociation coordinate R for the molecules LiH, Li2, and BH with the two outer valence electrons correlated. All of these results indicate that response properties are much more sensitive to deficiencies in the reconstruction of the 2RDM than the ground state energy, since derivatives of the functional with respect to both the NOs and the occupation numbers need to be accurate.
Meer, R. van; Gritsenko, O. V.; Baerends, E. J.
2014-01-14
Time dependent density matrix functional theory in its adiabatic linear response formulation delivers exact excitation energies ω{sub α} and oscillator strengths f{sub α} for two-electron systems if extended to the so-called phase including natural orbital (PINO) theory. The Löwdin-Shull expression for the energy of two-electron systems in terms of the natural orbitals and their phases affords in this case an exact phase-including natural orbital functional (PILS), which is non-primitive (contains other than just J and K integrals). In this paper, the extension of the PILS functional to N-electron systems is investigated. With the example of an elementary primitive NO functional (BBC1) it is shown that current density matrix functional theory ground state functionals, which were designed to produce decent approximations to the total energy, fail to deliver a qualitatively correct structure of the (inverse) response function, due to essential deficiencies in the reconstruction of the two-body reduced density matrix (2RDM). We now deduce essential features of an N-electron functional from a wavefunction Ansatz: The extension of the two-electron Löwdin-Shull wavefunction to the N-electron case informs about the phase information. In this paper, applications of this extended Löwdin-Shull (ELS) functional are considered for the simplest case, ELS(1): one (dissociating) two-electron bond in the field of occupied (including core) orbitals. ELS(1) produces high quality ω{sub α}(R) curves along the bond dissociation coordinate R for the molecules LiH, Li{sub 2}, and BH with the two outer valence electrons correlated. All of these results indicate that response properties are much more sensitive to deficiencies in the reconstruction of the 2RDM than the ground state energy, since derivatives of the functional with respect to both the NOs and the occupation numbers need to be accurate.
NASA Astrophysics Data System (ADS)
Mata-Pinzon, Zaahel; Valladares, Ariel Alberto; Valladares, Alexander; Valladares, Renela Maria
2014-03-01
The properties of materials are strongly related to their atomic topology, especially when we compare properties related to the ordered and disordered phases. Using Density Functional Theory methods on 64-atom supercells we obtain the structure and calculate the electronic density of states (eDOS) as a function of the concentration of lead, thallium or antimony in an amorphous bismuth supercell. This is done to investigate how the eDOS affects the superconducting transition temperature (Tc), taking into account the measurements made by Shier and Ginsberg[2] on contaminated amorphous bismuth thin films. Supported by CONACYT and DGAPA (UNAM).
NASA Astrophysics Data System (ADS)
Kurashige, Yuki
2014-06-01
Recent advances in quantum chemical density matrix renormalisation group (DMRG) theory are presented. The DMRG, originally devised as an alternative to the exact diagonalisation in condensed matter physics, has become a powerful quantum chemical method for molecular systems that exhibit a multireference character, e.g., excited states, π-conjugated systems, transition metal complexes, and in particular for large systems by combining it with conventional multireference electron correlation methods. The capability of the current quantum chemical DMRG is demonstrated for an application involving the potential energy curve of the chromium dimer, which is one of the most demanding multireference systems and thus requires the best electronic structure treatment for non-dynamical and dynamical correlation as well as large basis sets.
Predicting a new photocatalyst and its electronic properties by density functional theory
NASA Astrophysics Data System (ADS)
Sarker, Pranab; Prasher, Dixit; Gaillard, Nicolas; Huda, Muhammad N.
2013-10-01
A new material CuBiW2O8 is reported here which is suitable for photocatalysts for solar-to-hydrogen generation by splitting water through photoelectrochemical approach. By density functional theory total energy calculations along with extensive mineral database search of relevant oxides, the crystal structures of CuBiW2O8 has been determined, which agrees well with the experimental result. We have analyzed the thermodynamical stability of this material. Its stability was found to be comparable to other well-known oxides, such as CuWO4. The band structure calculation reveals that it has a suitable band gap. In addition to this, density of states and optical absorption calculations show favorable features of a photocatalyst.
NASA Astrophysics Data System (ADS)
Ferri, Nicola; DiStasio, Robert A.; Ambrosetti, Alberto; Car, Roberto; Tkatchenko, Alexandre
2015-05-01
How strong is the effect of van der Waals (vdW) interactions on the electronic properties of molecules and extended systems? To answer this question, we derived a fully self-consistent implementation of the density-dependent interatomic vdW functional of Tkatchenko and Scheffler [Phys. Rev. Lett. 102, 073005 (2009)]. Not surprisingly, vdW self-consistency leads to tiny modifications of the structure, stability, and electronic properties of molecular dimers and crystals. However, unexpectedly large effects were found in the binding energies, distances, and electrostatic moments of highly polarizable alkali-metal dimers. Most importantly, vdW interactions induced complex and sizable electronic charge redistribution in the vicinity of metallic surfaces and at organic-metal interfaces. As a result, a substantial influence on the computed work functions was found, revealing a nontrivial connection between electrostatics and long-range electron correlation effects.
Andrés, Juan; Berski, Sławomir; Silvi, Bernard
2016-07-01
Probing the electron density transfers during a chemical reaction can provide important insights, making possible to understand and control chemical reactions. This aim has required extensions of the relationships between the traditional chemical concepts and the quantum mechanical ones. The present work examines the detailed chemical insights that have been generated through 100 years of work worldwide on G. N. Lewis's ground breaking paper on The Atom and the Molecule (Lewis, G. N. The Atom and the Molecule, J. Am. Chem. Soc. 1916, 38, 762-785), with a focus on how the determination of reaction mechanisms can be reached applying the bonding evolution theory (BET), emphasizing how curly arrows meet electron density transfers in chemical reaction mechanisms and how the Lewis structure can be recovered. BET that combines the topological analysis of the electron localization function (ELF) and Thom's catastrophe theory (CT) provides a powerful tool providing insight into molecular mechanisms of chemical rearrangements. In agreement with physical laws and quantum theoretical insights, BET can be considered as an appropriate tool to tackle chemical reactivity with a wide range of possible applications. Likewise, the present approach retrieves the classical curly arrows used to describe the rearrangements of chemical bonds for a given reaction mechanism, providing detailed physical grounds for this type of representation. The ideas underlying the valence-shell-electron pair-repulsion (VSEPR) model applied to non-equilibrium geometries provide simple chemical explanations of density transfers. For a given geometry around a central atom, the arrangement of the electronic domain may comply or not with the VSEPR rules according with the valence shell population of the considered atom. A deformation yields arrangements which are either VSEPR defective (at least a domain is missing to match the VSEPR arrangement corresponding to the geometry of the ligands), VSEPR compliant
Magnetic properties of f-electron systems in spin-polarized relativistic density functional theory
NASA Astrophysics Data System (ADS)
Yamagami, H.; Mavromaras, A.; Kübler, J.
1997-12-01
The magnetic ground state of the series of lanthanide and actinide trivalent ions is investigated by means of spin-polarized relativistic spin-density functional theory. In the local density functional approximation (LDA) an internal effective magnetic field due to exchange and correlation couples to the spin degrees of freedom. The resulting set of coupled Dirac equations yields ground-state multiplets that obey the well-known Hund's rules. This remarkable result comes about by the coupling of the j = l + 1/2 with the j = l - 1/2 states due to the exchange - correlation potential that is, as usual, the functional derivative of the exchange - correlation energy with respect to the spin magnetic moment. The effect of the coupling is shown to depend on the varying relative strengths of spin - orbit coupling and exchange splitting within the f series. Since in the f levels the internal exchange splitting dominates rather than the spin - orbit splitting, the energy level scheme is that of the Paschen - Back effect, and thus features of the Russell - Saunders coupling persist in spite of relativistic effects.
Itinerant Electron Magnets: Curie Temperature and Susceptibility in Density-Functional Theory
NASA Astrophysics Data System (ADS)
Mohn, Peter; Khmelevskyi, Sergei
Groundstate properties of solids are astonishingly well described by the local density functional approximation (LDA) [1]. This is also true of metallic magnets for which the situation was recently described by this author [2]. Excited-state properties of magnets (and other systems), however, are still a great challenge and it was believed until recently that the band-picture, based in the LDA, fails entirely in describing magnetism at elevated temperatures. We emphasize here that this is not so, attempting first to expose the reason why it was thought that the band-picture fails. Since historically the underlying physical picture was developed by Stoner and Wohlfarth, we begin with a discussion of their theory using, however, an approach that reveals the essential assumptions. This is Mermin's [3] finite-temperature density functional theory. Two points emerge: one is the essential noncollinearity of the magnetic moments at finite temperatures, the other is the form of the exchange-corr elation contribution to the thermodynamic potential. We know how to deal with noncollinear order in the LDA and we explain how we might use this knowledge to advance the issue. Exchange and correlation at finite temperatures are, however, at the present state not well understood. This statement not only applies to density functional theory but also to many-body treatments addressed at this workshop. We will show in particular that the theory of magnons in the band-LDA-picture at low temperatures is in good shape. At high temperatures we opt for a theory involving spin fluctuations and argue that, although broad features of the magnetic phase transition are described satisfactorily, many details await further improvements.
NASA Astrophysics Data System (ADS)
Skone, Jonathan; Govoni, Marco; Galli, Giulia
Dielectric-dependent hybrid [DDH] functionals have recently been shown to yield highly accurate energy gaps and dielectric constants for a wide variety of solids, at a computational cost considerably less than standard GW calculations. The fraction of exact exchange included in the definition of DDH functionals depends (self-consistently) on the dielectric constant of the material. In the present talk we introduce a range-separated (RS) version of DDH functionals where short and long-range components are matched using material dependent, non-empirical parameters. Comparing with state of the art GW calculations and experiment, we show that such RS hybrids yield accurate electronic properties of both molecules and solids, including energy gaps, photoelectron spectra and absolute ionization potentials. This work was supported by NSF-CCI Grant Number NSF-CHE-0802907 and DOE-BES.
NASA Astrophysics Data System (ADS)
Rathi, Somilkumar; Ray, Asok
2008-03-01
Ab initio calculations within the framework of hybrid density functional theory and finite cluster approximation have been performed for the electronic and geometric structures of three different types of armchair germanium carbide nanotubes from (3, 3) to (11, 11). Full geometry and spin optimizations with unrestricted symmetry have been performed. A detailed comparison of the structures and stabilities of the three types of nanotubes will be presented. The dependence of the electronic band gaps on the respective tube diameters, energy density of states, dipole moments as well as Mulliken charge distributions have been investigated. Radial buckling of tube along with bond length variations is also studied. All armchair GeC nanotubes investigated so far are semiconducting in nature. Applications in the field of nano-optoelectronic devices, molecular electronics and band gap engineering are envisioned for GeC nanotubes.
NASA Astrophysics Data System (ADS)
Nolan, Michael; Grigoleit, Sonja; Sayle, Dean C.; Parker, Stephen C.; Watson, Graeme W.
2005-02-01
We present periodic density functional theory (DFT) calculations of bulk ceria and its low index surfaces (1 1 1), (1 1 0) and (1 0 0). We find that the surface energies increase in the order (1 1 1) > (1 1 0) > (1 0 0), while the magnitude of the surface relaxations follows the inverse order. The electronic properties of the bulk and surfaces are analysed by means of the electronic density of states and the electron density. We demonstrate that the bonding in pure ceria is partially covalent and analysis of the resulting electronic states confirms the presence of localised Ce 4f states above the Fermi level. The surface atoms show only a small change in the charge distribution in comparison to the bulk and from the DOS the main differences are due to the changes in the oxygen 2p and cerium 5 d states. Investigation of the atomic and electronic structure of an oxygen vacancy on the (1 0 0) surface shows the problems DFT can have with the description of strongly localised systems, wrongly predicting electron delocalisation over all of the cerium atoms in the simulation cell. We demonstrate an improvement in the description of the strongly correlated cerium 4f states in partially reduced ceria by applying the DFT+U methodology, which leads to the appearance of a new gap state between the valence band and the empty Ce 4f band. Analysis of the partial charge density shows that these states are localised on the Ce III ions neighbouring the oxygen vacancy. In terms of classical defect chemistry, the vacancy is bound by two neighbouring Ce III ions, which have been reduced from Ce IV, i.e. VO··+2CeCe'. The remaining Ce ions are in the Ce IV oxidation state. The localisation of Ce 4f electrons modifies the predicted structure of the defective surface.
Ensling, David; Thissen, Andreas; Laubach, Stefan; Schmidt, Peter C.; Jaegermann, Wolfram
2010-11-15
The electronic properties of LiCoO{sub 2} have been studied by theoretical band-structure calculations (using density functional theory) and experimental methods (photoemission). Synchrotron-induced photoelectron spectroscopy, resonant photoemission spectroscopy (ResPES), and soft x-ray absorption (XAS) have been applied to investigate the electronic structure of both occupied and unoccupied states. High-quality PES spectra were obtained from stoichiometric and highly crystalline LiCoO{sub 2} thin films deposited ''in situ'' by rf magnetron sputtering. An experimental approach of separating oxygen- and cobalt-derived (final) states by ResPES in the valence-band region is presented. The procedure takes advantage of an antiresonant behavior of cobalt-derived states at the 3p-3d excitation threshold. Information about the unoccupied density of states has been obtained by O K XAS. The structure of the Co L absorption edge is compared to semiempirical charge-transfer multiplet calculations. The experimental results are furthermore compared with band-structure calculations considering three different exchange potentials [generalized gradient approximation (GGA), using a nonlocal Hubbard U (GGA+U) and using a hybrid functional (Becke, three-parameter, Lee-Yang-Parr [B3LYP])]. For these different approaches total density of states and partial valence-band density of states have been investigated. The best qualitative agreement with experimental results has been obtained by using a GGA+U functional with U=2.9 eV.
Effects of d-electrons in pseudopotential screened-exchange density functional calculations
NASA Astrophysics Data System (ADS)
Lee, Byounghak; Wang, Lin-Wang; Canning, Andrew
2008-06-01
We report a theoretical study on the role of shallow d states in the screened-exchange local density approximation (sX-LDA) band structure of binary semiconductor systems. We found that inaccurate pseudo-wave functions can lead to (1) an overestimation of the screened-exchange interaction between the localized d states and the delocalized higher energy s and p states, and (2) an underestimation of the screened-exchange interaction between the d states. The resulting sX-LDA band structures have substantially smaller band gaps compared with experiments. We correct the pseudo-wave functions of d states by including the semicore s and p states of the same shell in the valence states. The correction of pseudo-wave functions yields band gaps and d-state binding energies in good agreement with experiments and the full potential linearized augmented plane wave sX-LDA calculations. Compared with the quasiparticle GW method, our sX-LDA results shows not only similar quality on the band gaps but also much better d-state binding energies. Combined with its capability of ground-state structure calculation, the sX-LDA is expected to be a valuable theoretical tool for the II-VI and III-V (especially the III-N) bulk semiconductors and nanostructure studies.
Density Functionals of Chemical Bonding
Putz, Mihai V.
2008-01-01
The behavior of electrons in general many-electronic systems throughout the density functionals of energy is reviewed. The basic physico-chemical concepts of density functional theory are employed to highlight the energy role in chemical structure while its extended influence in electronic localization function helps in chemical bonding understanding. In this context the energy functionals accompanied by electronic localization functions may provide a comprehensive description of the global-local levels electronic structures in general and of chemical bonds in special. Becke-Edgecombe and author’s Markovian electronic localization functions are discussed at atomic, molecular and solid state levels. Then, the analytical survey of the main workable kinetic, exchange, and correlation density functionals within local and gradient density approximations is undertaken. The hierarchy of various energy functionals is formulated by employing both the parabolic and statistical correlation degree of them with the electronegativity and chemical hardness indices by means of quantitative structure-property relationship (QSPR) analysis for basic atomic and molecular systems. PMID:19325846
NASA Astrophysics Data System (ADS)
Çatıkkaş, Berna; Aktan, Ebru; Yalçın, Ergin
2016-08-01
This work deals with the optimized molecular structure, vibrational spectra, nonlinear optic (NLO) and frontier molecule orbital (FMO) properties of 1-Methyl-2-phenyl-3-(1,3,4-thiadiazol-2-yldiazenyl)-1H-indole (MPI) by quantum chemical calculations. The Fourier transform infrared (FT-MIR and FT-FIR) and Raman spectra of 1-Methyl-2-phenyl-3-(1,3,4-thiadiazol-2-yldiazenyl)-1H-indole (MPI) were recorded in the region (4000-400 cm-1 and 400-30 cm-1) and (3200-92 cm-1), respectively. The analysis and complete vibrational assignments of the fundamental modes of the MPI molecule were carried out by using the observed FT-IR and FT-Raman data and calculated Total Energy Distribution (TED) according to Scaled Quantum Mechanics procedure. The calculated geometrical parameters of the MPI molecule are in agreement with the obtained values from XRD studies. On the other hand, the difference between the scaled and observed wavenumber values of the most of the fundamentals are very small. 1H NMR and 13C NMR chemical shift values, and energy gap between LUMO-HOMO and molecular electrostatic potential (MEP) were investigated by using density functional theory (B3LYP) methods. UV/Visible spectra and λ maximum absorption values, the oscillator strengths in the chloroform, methanol and DMSO solvation in combination with different basis sets were calculated by using the time-dependent density functional theory (TD-DFT). Additionally, the predicted nonlinear optical (NLO) properties of the MPI are quite greater than that of urea at the B3LYP/6-31++G(d,p) level.
NASA Astrophysics Data System (ADS)
Hung, L.; Guedj, C.; Bernier, N.; Blaise, P.; Olevano, V.; Sottile, F.
2016-04-01
We present the valence electron energy-loss spectrum and the dielectric function of monoclinic hafnia (m -HfO2) obtained from time-dependent density-functional theory (TDDFT) predictions and compared to energy-filtered spectroscopic imaging measurements in a high-resolution transmission-electron microscope. Fermi's golden rule density-functional theory (DFT) calculations can capture the qualitative features of the energy-loss spectrum, but we find that TDDFT, which accounts for local-field effects, provides nearly quantitative agreement with experiment. Using the DFT density of states and TDDFT dielectric functions, we characterize the excitations that result in the m -HfO2 energy-loss spectrum. The sole plasmon occurs between 13 and 16 eV, although the peaks ˜28 and above 40 eV are also due to collective excitations. We furthermore elaborate on the first-principles techniques used, their accuracy, and remaining discrepancies among spectra. More specifically, we assess the influence of Hf semicore electrons (5 p and 4 f ) on the energy-loss spectrum, and find that the inclusion of transitions from the 4 f band damps the energy-loss intensity in the region above 13 eV. We study the impact of many-body effects in a DFT framework using the adiabatic local-density approximation (ALDA) exchange-correlation kernel, as well as from a many-body perspective using "scissors operators" matched to an ab initio G W calculation to account for self-energy corrections. These results demonstrate some cancellation of errors between self-energy and excitonic effects, even for excitations from the Hf 4 f shell. We also simulate the dispersion with increasing momentum transfer for plasmon and collective excitation peaks.
Borges, P. D. E-mail: lscolfaro@txstate.edu; Scolfaro, L. E-mail: lscolfaro@txstate.edu
2014-12-14
The thermoelectric properties of indium nitride in the most stable wurtzite phase (w-InN) as a function of electron and hole concentrations and temperature were studied by solving the semiclassical Boltzmann transport equations in conjunction with ab initio electronic structure calculations, within Density Functional Theory. Based on maximally localized Wannier function basis set and the ab initio band energies, results for the Seebeck coefficient are presented and compared with available experimental data for n-type as well as p-type systems. Also, theoretical results for electric conductivity and power factor are presented. Most cases showed good agreement between the calculated properties and experimental data for w-InN unintentionally and p-type doped with magnesium. Our predictions for temperature and concentration dependences of electrical conductivity and power factor revealed a promising use of InN for intermediate and high temperature thermoelectric applications. The rigid band approach and constant scattering time approximation were utilized in the calculations.
Johansson, Malin B; Baldissera, Gustavo; Valyukh, Iryna; Persson, Clas; Arwin, Hans; Niklasson, Gunnar A; Osterlund, Lars
2013-05-22
The optical and electronic properties of nanocrystalline WO3 thin films prepared by reactive dc magnetron sputtering at different total pressures (Ptot) were studied by optical spectroscopy and density functional theory (DFT) calculations. Monoclinic films prepared at low Ptot show absorption in the near infrared due to polarons, which is attributed to a strained film structure. Analysis of the optical data yields band-gap energies Eg ≈ 3.1 eV, which increase with increasing Ptot by 0.1 eV, and correlate with the structural modifications of the films. The electronic structures of triclinic δ-WO3, and monoclinic γ- and ε-WO3 were calculated using the Green function with screened Coulomb interaction (GW approach), and the local density approximation. The δ-WO3 and γ-WO3 phases are found to have very similar electronic properties, with weak dispersion of the valence and conduction bands, consistent with a direct band-gap. Analysis of the joint density of states shows that the optical absorption around the band edge is composed of contributions from forbidden transitions (>3 eV) and allowed transitions (>3.8 eV). The calculations show that Eg in ε-WO3 is higher than in the δ-WO3 and γ-WO3 phases, which provides an explanation for the Ptot dependence of the optical data. PMID:23614973
Maaloul, L.; Gangwar, R. K.; Morel, S.; Stafford, L.
2015-11-15
Langmuir probe and trace rare gases optical emission spectroscopy were used to analyze the spatial structure of the electron density and electron energy distribution function (EEDF) in a cylindrical Ar magnetron plasma reactor used for sputter-deposition of ZnO-based thin films. While a typical Bessel (zero order) diffusion profile was observed along the radial direction for the number density of charged particles at 21 cm from the ZnO target, a significant rise of these populations with respect to the Bessel function was seen in the center of the reactor at 4 cm from the magnetron surface. As for the EEDF, it was found to transform from a more or less Maxwellian far from the target to a two-temperature Maxwellian with a depletion of high-energy electrons where magnetic field confinement effects become important. No significant change in the behavior of the electron density and EEDF across a wide range of pressures (5–100 mTorr) and self-bias voltages (115–300 V) was observed during magnetron sputtering of Zn, ZnO, and In{sub 2}O{sub 3} targets. This indicates that sputtering of Zn, In, and O atoms do not play a very significant role on the electron particle balance and electron heating dynamics, at least over the range of experimental conditions investigated.
NASA Astrophysics Data System (ADS)
Polin, Daniel; Ziegler, Joshua; Malozovsky, Yuriy; Bagayoko, Diola
We present the findings of ab-initio calculations of electronic, transport, and structural properties of cubic sodium oxide (Na2O). These results were obtained using density functional theory (DFT), specifically a local density approximation (LDA) potential, and the linear combination of Gaussian orbitals (LCGO). Our implementation of LCGO followed the Bagayoko, Zhao, and Williams method as enhanced by the work of Ekuma and Franklin (BZW-EF). We describe the electronic band structure of Na2O with a direct band gap of 2.22 eV. Our results include predicted values for the electronic band structure and associated energy eigenvalues, the total and partial density of states (DOS and pDOS), the equilibrium lattice constant of Na2O, and the bulk modulus. We have also calculated the electron and holes effective masses in the Γ to L, Γ to X, and Γ to K directions. Acknowledgments: This work was funded in part by the National Science Foundation (NSF) and the Louisiana Board of Regents, through LASiGMA [Award Nos. EPS- 1003897, NSF (2010-15)-RII-SUBR] and NSF HRD-1002541, the US Department of Energy - National, Nuclear Security Administration (NNSA) (Award No. DE- NA0002630), LaSPACE, and LONI-SUBR.
Filatov, Michael; Huix-Rotllant, Miquel; Burghardt, Irene
2015-05-14
State-averaged (SA) variants of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, SA-REKS and state-interaction (SI)-SA-REKS, implement ensemble density functional theory for variationally obtaining excitation energies of molecular systems. In this work, the currently existing version of the SA-REKS method, which included only one excited state into the ensemble averaging, is extended by adding more excited states to the averaged energy functional. A general strategy for extension of the REKS-type methods to larger ensembles of ground and excited states is outlined and implemented in extended versions of the SA-REKS and SI-SA-REKS methods. The newly developed methods are tested in the calculation of several excited states of ground-state multi-reference systems, such as dissociating hydrogen molecule, and excited states of donor–acceptor molecular systems. For hydrogen molecule, the new method correctly reproduces the distance dependence of the lowest excited state energies and describes an avoided crossing between the doubly excited and singly excited states. For bithiophene–perylenediimide stacked complex, the SI-SA-REKS method correctly describes crossing between the locally excited state and the charge transfer excited state and yields vertical excitation energies in good agreement with the ab initio wavefunction methods.
Yan, Xun-Wang; Faculty of Physics and Electronic Technology, Hubei University, Wuhan 430062; State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Science, Beijing 100190, China and School of Physics and Electrical Engineering, Anyang Normal University, Henan 455000 ; Huang, Zhongbing; Faculty of Physics and Electronic Technology, Hubei University, Wuhan 430062 ; Lin, Hai-Qing
2013-11-28
By the first principle calculations based on the van der Waals density functional theory, we study the crystal structures and electronic properties of La-doped phenanthrene. Two stable atomic geometries of La{sub 1}phenanthrene are obtained by relaxation of atomic positions from various initial structures. The structure-I is a metal with two energy bands crossing the Fermi level, while the structure-II displays a semiconducting state with an energy gap of 0.15 eV, which has an energy gain of 0.42 eV per unit cell compared to the structure-I. The most striking feature of La{sub 1}phenanthrene is that La 5d electrons make a significant contribution to the total density of state around the Fermi level, which is distinct from potassium doped phenanthrene and picene. Our findings provide an important foundation for the understanding of superconductivity in La-doped phenanthrene.
Ramesh Babu, K; Vaitheeswaran, G
2014-06-11
We report a systematic first principles density functional study on the electronic structure, elastic and optical properties of nitrogen based solid hydrogen storage materials LiNH2, NaNH2, KNH2, and RbNH2. The ground state structural properties are calculated by using standard density functional theory, and also dispersion corrected density functional theory. We find that van der Waals interactions are dominant in LiNH2 whereas they are relatively weak in other alkali metal amides. The calculated elastic constants show that all the compounds are mechanically stable and LiNH2 is found to be a stiffer material among the alkali metal amides. The melting temperatures are calculated and follow the order RbNH2 < KNH2 < NaNH2 < LiNH2. The electronic band structure is calculated by using the Tran–Blaha modified Becke–Johnson potential and found that all the compounds are insulators, with a considerable band gap. The [NH2]− derived states completely dominate in the entire valence band region while the metal atom states occupy the conduction band. The calculated band structure is used to analyze the different interband optical transitions occurring between valence and conduction bands. Our calculations show that these materials have considerable optical anisotropy. PMID:25932472
NASA Astrophysics Data System (ADS)
Gori-Giorgi, Paola; Perdew, John P.
2002-10-01
We construct analytic formulas that represent the coupling-constant-averaged pair distribution function gxc(rs,ζ,kFu) of a three-dimensional nonrelativistic ground-state electron gas constrained to a uniform density with density parameter rs=(9π/4)1/3/kF and relative spin polarization ζ over the whole range 0
The magnetic and electronic structure of vanadyl pyrophosphate from density functional theory.
Cheng, Mu-Jeng; Nielsen, Robert J; Tahir-Kheli, Jamil; Goddard, William A
2011-05-28
We have studied the magnetic structure of the high symmetry vanadyl pyrophosphate ((VO)(2)P(2)O(7), VOPO), focusing on the spin exchange couplings, using density functional theory (B3LYP) with the full three-dimensional periodicity. VOPO involves four distinct spin couplings: two larger couplings exist along the chain direction (a-axis), which we predict to be antiferromagnetic, J(OPO) = -156.8 K and J(O) = -68.6 K, and two weaker couplings appear along the c (between two layers) and b directions (between two chains in the same layer), which we calculate to be ferromagnetic, J(layer) = 19.2 K and J(chain) = 2.8 K. Based on the local density of states and the response of spin couplings to varying the cell parameter a, we found that J(OPO) originates from a super-exchange interaction through the bridging -O-P-O- unit. In contrast, J(O) results from a direct overlap of 3d(x(2)-y(2)) orbitals on two vanadium atoms in the same V(2)O(8) motif, making it very sensitive to structural fluctuations. Based on the variations in V-O bond length as a function of strain along a, we found that the V-O bonds of V-(OPO)(2)-V are covalent and rigid, whereas the bonds of V-(O)(2)-V are fragile and dative. These distinctions suggest that compression along the a-axis would have a dramatic impact on J(O), changing the magnetic structure and spin gap of VOPO. This result also suggests that assuming J(O) to be a constant over the range of 2-300 K whilst fitting couplings to the experimental magnetic susceptibility is an invalid method. Regarding its role as a catalyst, the bonding pattern suggests that O(2) can penetrate beyond the top layers of the VOPO surface, converting multiple V atoms from the +4 to +5 oxidation state, which seems crucial to explain the deep oxidation of n-butane to maleic anhydride. PMID:21503330
A van der Waals density functional built upon the electron-gas foundation
NASA Astrophysics Data System (ADS)
Hyldgaard, Per; Berland, Kristian; Schröder, Elsebeth
2015-03-01
The vdW-DF method is designed to be a systematic extension of the constraint-based generalized-gradient approximation (GGA) and can therefore serve as general purpose density functional [PRB 90, 075148 (2014)]. Yet the early versions can have issues both with bulk systems and with weak chemisorption. We present a recent nonempirical version, vdW-DF-cx [J. Chem. Phys. 140, 18A539 (2014), PRB 89, 035412 (2014)], that resolves these issues. The version is designed to have a consistent combination of exchange and correlation. We show that it performs well for inter-molecular binding and that it can even be better than PBE for describing cohesion and structure of molecules and solids. These results validate the robustness of the vdW-DF plasmon-pole model, which we show is closed linked to the exchange correlation hole of constraint-based GGA. The work was supported by the Swedish Research Council (VR), by the Chalmers Areas of Advance: Materials, and by the Swedish National Infrastructure for Computing.
NASA Astrophysics Data System (ADS)
Suhasini, M.; Sailatha, E.; Gunasekaran, S.; Ramkumaar, G. R.
2015-11-01
The Fourier transform infrared (FT-IR) and FT-Raman spectra of Lornoxicam were recorded in the region 4000-450 cm-1 and 4000-50 cm-1 respectively. Density functional theory (DFT) has been used to calculate the optimized geometrical parameters, atomic charges, and vibrational wavenumbers and intensity of the vibrational bands. The computed vibrational wave numbers were compared with the FT-IR and FT-Raman experimental data. The computational calculations at DFT/B3LYP level with 6-31G(d,p) and 6-31++G(d,p) basis sets. The complete vibrational assignments were performed on the basis of the potential energy distribution (PED) of the Vibrational modes calculated using Vibrational Energy Distribution Analysis (VEDA 4) program. The oscillator's strength calculated by TD-DFT and Lornoxicam is approach complement with the experimental findings. The NMR chemical shifts 13C and 1H were recorded and calculated using the gauge independent atomic orbital (GIAO) method. The Natural charges and intermolecular contacts have been interpreted using Natural Bond orbital (NBO) analysis and the HOMO-LUMO energy gap has been calculated. The thermodynamic properties like Entropy, Enthalpy, Specific heat capacity and zero vibrational energy have been calculated. Besides, molecular electrostatic potential (MEP) was investigated using theoretical calculations.
A hybrid-exchange density functional study of the bonding and electronic structure in bulk CuFeS2.
Martínez-Casado, Ruth; Chen, Vincent H-Y; Mallia, Giuseppe; Harrison, Nicholas M
2016-05-14
The geometric, electronic, and magnetic properties of bulk chalcopyrite CuFeS2 have been investigated using hybrid-exchange density functional theory calculations. The results are compared with available theoretical and experimental data. The theoretical description of the bonding and electronic structure in CuFeS2 is analyzed in detail and compared to those computed for chalcocite (CuS2) and greigite (Fe3S4). In particular, the behavior of the 3d electrons of Fe(3+) is discussed in terms of the Hubbard-Anderson model in the strongly correlated regime and found to be similarly described in both materials by an on-site Coulomb repulsion (U) of ∼8.9 eV and a transfer integral (t) of ∼0.3 eV. PMID:27179495
A hybrid-exchange density functional study of the bonding and electronic structure in bulk CuFeS2
NASA Astrophysics Data System (ADS)
Martínez-Casado, Ruth; Chen, Vincent H.-Y.; Mallia, Giuseppe; Harrison, Nicholas M.
2016-05-01
The geometric, electronic, and magnetic properties of bulk chalcopyrite CuFeS2 have been investigated using hybrid-exchange density functional theory calculations. The results are compared with available theoretical and experimental data. The theoretical description of the bonding and electronic structure in CuFeS2 is analyzed in detail and compared to those computed for chalcocite (CuS2) and greigite (Fe3S4). In particular, the behavior of the 3d electrons of Fe3+ is discussed in terms of the Hubbard-Anderson model in the strongly correlated regime and found to be similarly described in both materials by an on-site Coulomb repulsion (U) of ˜8.9 eV and a transfer integral (t) of ˜0.3 eV.
NASA Astrophysics Data System (ADS)
Nawa, Kenji; Kitaoka, Yukie; Nakamura, Kohji; Imamura, Hiroshi; Akiyama, Toru; Ito, Tomonori; Weinert, M.
2016-07-01
The ground-state electronic configurations of the correlated organometallic metallocenes, M Cp2,M =V , Cr, Mn, Fe, Co, and Ni, are investigated using constraint density functional theory combined with nonempirical Ueff parameters determined from linear-response theory. The relative stability of the various d -orbital electronic configurations of these organometallic molecules is found to be sensitive to the amount of correlation. Using nonempirical values of Ueff, the calculated electronic configurations are in agreement with the experiments: 4A2 g ,3E2 g ,6A1 g ,1A1 g ,2E1 g , and 3A2 g for the VCp2,CrCp2,MnCp2,FeCp2,CoCp2 , and NiCp2, respectively.
NASA Astrophysics Data System (ADS)
Mo, Chongjie; Yang, Yu; Kang, Wei; Zhang, Ping
2016-04-01
The electronic structure and optical properties for the (U,Th)O2 compound are systematically studied by employing the Heyd-Scuseria-Ernzerh method (HSE) of screened hybrid density functional. The electronic band gap of (U,Th)O2 is predicted to be 3.06 eV, in the middle of the values of UO2 and ThO2. Based on wavefunction analysis, we conclude (U,Th)O2 to be a Mott insulator in its ground state. The frequency dependent dielectric functions and optical properties are then calculated and compared with those of ThO2 and UO2. At the visible light frequency range, the adsorption coefficients for ThO2, UO2 and (U,Th)O2 are totally different, which gives an accessible method to predict the proportion of U atoms in an arbitrary unknown (U,Th)O2 compounds from the adsorption spectrum of visible lights.
Interstellar Electron Density Spectra
NASA Astrophysics Data System (ADS)
Lambert, Hendrick Clark
This study concerns the investigation of the form of the wavenumber spectrum of the Galactic electron density fluctuations through an examination of the scattering of the radio pulses emitted by pulsars as they propagate through the diffuse ionized interstellar gas. A widely used model for the electron density spectrum is based on the simple power-law: Pne(q)∝ q-β, where β = 11/3 is usually assumed, corresponding to Kolmogorov's turbulence spectrum. The simple Kolmogorov model provides satisfactory agreement for observations along many lines of sight; however, major inconsistencies remain. The inconsistencies suggest that an increase in the ratio of the power between the high (10-8[ m]-1≤ q<=10-7[ m]-1) and low (10-13[ m]-1≤ q<=10-12[ m]-1) wavenumbers is needed. This enhancement in the ratio can in turn be achieved by either including an inner scale, corresponding to a dissipation scale for the turbulent cascade, in the Kolmogorov spectrum or by considering steeper spectra. Spectra with spectral exponents β > 4 have been in general rejected based on observations of pulsar refractive scintillations. The special case of β = 4 has been given little attention and is analyzed in detail. Physically, this 'β = 4' model corresponds to the random distribution, both in location and orientation, of discrete objects with relatively sharp boundaries across the line of sight. An outer scale is included in the model to account for the average size of such objects. We compare the predictions of the inner-scale and β = 4 models both with published observations and observations we made as part of this investigation. We conclude that the form of the wavenumber spectrum is dependent on the line of sight. We propose a composite spectrum featuring a uniform background turbulence in presence of randomly distributed discrete objects, as modeled by the β = model.
NASA Astrophysics Data System (ADS)
Grimme, Stefan
2013-06-01
Two approximations in the Tamm-Dancoff density functional theory approach (TDA-DFT) to electronically excited states are proposed which allow routine computations for electronic ultraviolet (UV)- or circular dichroism (CD) spectra of molecules with 500-1000 atoms. Speed-ups compared to conventional time-dependent DFT (TD-DFT) treatments of about two to three orders of magnitude in the excited state part at only minor loss of accuracy are obtained. The method termed sTDA ("s" for simplified) employs atom-centered Löwdin-monopole based two-electron repulsion integrals with the asymptotically correct 1/R behavior and perturbative single excitation configuration selection. It is formulated generally for any standard global hybrid density functional with given Fock-exchange mixing parameter ax. The method performs well for two standard benchmark sets of vertical singlet-singlet excitations for values of ax in the range 0.2-0.6. The mean absolute deviations from reference data are only 0.2-0.3 eV and similar to those from standard TD-DFT. In three cases (two dyes and one polypeptide), good mutual agreement between the electronic spectra (up to 10-11 eV excitation energy) from the sTDA method and those from TD(A)-DFT is obtained. The computed UV- and CD-spectra of a few typical systems (e.g., C60, two transition metal complexes, [7]helicene, polyalanine, a supramolecular aggregate with 483 atoms and about 7000 basis functions) compare well with corresponding experimental data. The method is proposed together with medium-sized double- or triple-zeta type atomic-orbital basis sets as a quantum chemical tool to investigate the spectra of huge molecular systems at a reliable DFT level.
Thakkar, Ajit J. Wu, Taozhe
2015-10-14
Static electronic dipole polarizabilities for 135 molecules are calculated using second-order Møller-Plesset perturbation theory and six density functionals recently recommended for polarizabilities. Comparison is made with the best gas-phase experimental data. The lowest mean absolute percent deviations from the best experimental values for all 135 molecules are 3.03% and 3.08% for the LC-τHCTH and M11 functionals, respectively. Excluding the eight extreme outliers for which the experimental values are almost certainly in error, the mean absolute percent deviation for the remaining 127 molecules drops to 2.42% and 2.48% for the LC-τHCTH and M11 functionals, respectively. Detailed comparison enables us to identify 32 molecules for which the discrepancy between the calculated and experimental values warrants further investigation.
Sato, Shunsuke A; Taniguchi, Yasutaka; Shinohara, Yasushi; Yabana, Kazuhiro
2015-12-14
We develop methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta-generalized gradient approximation was proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional was proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt the predictor-corrector step for a stable time evolution. We have developed a method to evaluate electronic excitation energy without referring to the energy functional which is unknown for the TB-mBJ potential. For the HSE functional, we have developed a method for the operation of the Fock-like term in Fourier space to facilitate efficient use of massive parallel computers equipped with graphic processing units. We compare electronic excitations in silicon and germanium induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: they are close to each other using TB-mBJ and HSE and are much smaller in LDA. At high laser intensities close to the damage threshold, electronic excitation energies do not differ much among the three cases. PMID:26671367
NASA Astrophysics Data System (ADS)
Sato, Shunsuke A.; Taniguchi, Yasutaka; Shinohara, Yasushi; Yabana, Kazuhiro
2015-12-01
We develop methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta-generalized gradient approximation was proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional was proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt the predictor-corrector step for a stable time evolution. We have developed a method to evaluate electronic excitation energy without referring to the energy functional which is unknown for the TB-mBJ potential. For the HSE functional, we have developed a method for the operation of the Fock-like term in Fourier space to facilitate efficient use of massive parallel computers equipped with graphic processing units. We compare electronic excitations in silicon and germanium induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: they are close to each other using TB-mBJ and HSE and are much smaller in LDA. At high laser intensities close to the damage threshold, electronic excitation energies do not differ much among the three cases.
Sato, Shunsuke A.; Taniguchi, Yasutaka; Shinohara, Yasushi; Yabana, Kazuhiro
2015-12-14
We develop methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta-generalized gradient approximation was proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional was proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt the predictor-corrector step for a stable time evolution. We have developed a method to evaluate electronic excitation energy without referring to the energy functional which is unknown for the TB-mBJ potential. For the HSE functional, we have developed a method for the operation of the Fock-like term in Fourier space to facilitate efficient use of massive parallel computers equipped with graphic processing units. We compare electronic excitations in silicon and germanium induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: they are close to each other using TB-mBJ and HSE and are much smaller in LDA. At high laser intensities close to the damage threshold, electronic excitation energies do not differ much among the three cases.
Baruah, Tunna; Garnica, Amanda; Paggen, Marina; Basurto, Luis; Zope, Rajendra R
2016-04-14
We study the electronic structure of C60 fullerenes functionalized with a thiophene-diketo-pyrrolopyrrole-thiophene based chromophore using density functional theory combined with large polarized basis sets. As the attached chromophore has electron donor character, the functionalization of the fullerene leads to a donor-acceptor (DA) system. We examine in detail the effect of the linker and the addition site on the electronic structure of the functionalized fullerenes. We further study the electronic structure of these DA complexes with a focus on the charge transfer excitations. Finally, we examine the interface of the functionalized fullerenes with the widely used poly(3-hexylthiophene-2,5-diyl) (P3HT) donor. Our results show that all functionalized fullerenes with an exception of the C60-pyrrolidine [6,6], where the pyrrolidine is attached at a [6,6] site, have larger electron affinities relative to the pristine C60 fullerene. We also estimate the quasi-particle gap, lowest charge transfer excitation energy, and the exciton binding energies of the functionalized fullerene-P3MT model systems. Results show that the exciton binding energies in these model complexes are slightly smaller compared to a similarly prepared phenyl-C61-butyric acid methyl ester (PCBM)-P3MT complex. PMID:27083718
NASA Astrophysics Data System (ADS)
Baruah, Tunna; Garnica, Amanda; Paggen, Marina; Basurto, Luis; Zope, Rajendra R.
2016-04-01
We study the electronic structure of C60 fullerenes functionalized with a thiophene-diketo-pyrrolopyrrole-thiophene based chromophore using density functional theory combined with large polarized basis sets. As the attached chromophore has electron donor character, the functionalization of the fullerene leads to a donor-acceptor (DA) system. We examine in detail the effect of the linker and the addition site on the electronic structure of the functionalized fullerenes. We further study the electronic structure of these DA complexes with a focus on the charge transfer excitations. Finally, we examine the interface of the functionalized fullerenes with the widely used poly(3-hexylthiophene-2,5-diyl) (P3HT) donor. Our results show that all functionalized fullerenes with an exception of the C60-pyrrolidine [6,6], where the pyrrolidine is attached at a [6,6] site, have larger electron affinities relative to the pristine C60 fullerene. We also estimate the quasi-particle gap, lowest charge transfer excitation energy, and the exciton binding energies of the functionalized fullerene-P3MT model systems. Results show that the exciton binding energies in these model complexes are slightly smaller compared to a similarly prepared phenyl-C61-butyric acid methyl ester (PCBM)-P3MT complex.
Uzunova, Ellie L.; Mikosch, Hans
2014-01-14
The iron oxide dimers (FeO){sub 2} and their peroxide isomers are studied with the B3LYP density functional as bare clusters and as hexacarbonyls. Among the bare clusters the planar four-member ring structures are more stable than the non-planar ones and the rhombic dioxide Fe{sub 2}O{sub 2} with antiferromagnetically ordered electrons on iron centers is the global minimum. Water adsorption on the bare diiron dioxide is exothermic, but dissociation does not occur. Carbonylation favors a non-planar Fe{sub 2}O{sub 2} ring for both the dioxides and the peroxides and high electron density at the Fe centers is induced, evidenced by the natural charge distribution, the high proton affinity, and the values of global electronegativity and hardness. The iron dioxide hexacarbonyl Fe{sub 2}O{sub 2}(CO){sub 6} is diamagnetic in the state of the global minimum. It is separated from the next low-lying triplet state by a small energy gap of 0.22 eV. Time-dependent density functional theory methods were applied to examine electron excitations from the ground state to the low-lying triplet states in the hexacarbonyls and their adsorption complexes with water. Singlet-to-triplet state excitations occur via ligand-to-metal charge transfer in the hexacarbonyls; in the adsorption complexes excitations from the oxygen lone pairs to the adsorption center also occur and they appear in the IR-visible region. The lowest energy singlet and triplet state reaction paths for water splitting were followed. On the singlet potential energy surface (PES), water splitting is spontaneous, while for the triplet PES an activation barrier of 14.1 kJ mol{sup −1} was determined.
NASA Astrophysics Data System (ADS)
Nawa, Kenji; Nakamura, Kohji; Akiyama, Toru; Ito, Tomonori; Weinert, Michael
Effective on-site Coulomb interactions (Ueff) and electron configurations in the localized d and f orbitals of metal complexes in transition-metal oxides and organometallic molecules, play a key role in the first-principles search for the true ground-state. However, wide ranges of values in the Ueff parameter of a material, even in the same ionic state, are often reported. Here, we revisit this issue from constraint density functional theory (DFT) by using the full-potential linearized augmented plane wave method. The Ueff parameters for prototypical transition-metal oxides, TMO (TM =Mn, Fe, Co, Ni), were calculated by the second derivative of the total energy functional with respect to the d occupation numbers inside the muffin-tin (MT) spheres as a function of the sphere radius. We find that the calculated Ueff values depend significantly on the MT radius, with a variation of more than 3 eV when the MT radius changes from 2.0 to 2.7 a.u., but importantly an identical valence band structure can be produced in all the cases, with an approximate scaling of Ueff. This indicates that a simple transferability of the Ueff value among different calculation methods is not allowed. We further extend the constraint DFT to treat various electron configurations of the localized d-orbitals in organometallic molecules, TMCp2 (TM =Cr, Mn, Fe, Co, Ni), and find that the calculated Ueff values can reproduce the experimentally determined ground-state electron configurations.
NASA Astrophysics Data System (ADS)
Marana, Naiara L.; Albuquerque, Anderson R.; La Porta, Felipe A.; Longo, Elson; Sambrano, Julio R.
2016-05-01
Periodic density functional theory calculations with the B3LYP hybrid functional and all-electron Gaussian basis set were performed to simulate the structural and electronic properties as well as the strain and formation energies of single-walled ZnO nanotubes (SWZnONTs) and Carbon nanotubes (SWCNTs) with different chiralities as functions of their diameters. For all SWZnONTs, the band gap, strain energy, and formation energy converge to ~4.5 eV, 0.0 eV/atom, and 0.40 eV/atom, respectively. This result suggests that the nanotubes are formed more easily from the surface than from the bulk. For SWCNTs, the strain energy is always positive, while the formation energy is negative for armchair and zigzag nanotubes, therefore suggesting that these types of nanotubes can be preferentially formed from the bulk. The electronic properties of SWCNTs depend on the chirality; all armchair nanotubes are metallic, while zigzag and chiral nanotubes can be metallic or semiconducting, depending on the n and m vectors.
NASA Astrophysics Data System (ADS)
Rivelino, Roberto; Malaspina, Thaciana; Fileti, Eudes E.
2009-01-01
We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C60(OH)n] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke’s three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C60(OH)36 . Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1 , 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C60(OH)24 isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4 , 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.
NASA Astrophysics Data System (ADS)
Kastlunger, Georg; Stadler, Robert
2015-03-01
To define the conductance of single-molecule junctions with a redox functionality in an electrochemical cell, two conceptually different electron transport mechanisms, namely, coherent tunneling and vibrationally induced hopping, compete with each other, where implicit parameters of the setup such as the length of the molecule and the applied gate voltage decide which mechanism is the dominant one. Although coherent tunneling is most efficiently described within Landauer theory and the common theoretical treatment of electron hopping is based on Marcus theory, both theories are adequate for the processes they describe without introducing accuracy-limiting approximations. For a direct comparison, however, it has to be ensured that the crucial quantities obtained from electronic structure calculations, i.e., the transmission function T (E ) in Landauer theory and the transfer integral V , the reorganization energy λ , and the driving force Δ G0 in Marcus theory, are derived from similar grounds, as pointed out by Nitzan and coworkers in a series of publications. In this paper our framework is a single-particle picture, for which we perform density functional theory calculations for the conductance corresponding to both transport mechanisms for junctions with the central molecule containing one, two, or three Ruthenium centers, from which we extrapolate our results in order to define the critical length of the transition point of the two regimes which we identify at 5.76nm for this type of molecular wire. We also discuss trends in the dependence on an electrochemically induced gate potential.
Density Functional Study of Structures and Electron Affinities of BrO4F/BrO4F−
Gong, Liangfa; Xiong, Jieming; Wu, Xinmin; Qi, Chuansong; Li, Wei; Guo, Wenli
2009-01-01
The structures, electron affinities and bond dissociation energies of BrO4F/BrO4F− species have been investigated with five density functional theory (DFT) methods with DZP++ basis sets. The planar F-Br…O2…O2 complexes possess 3A′ electronic state for neutral molecule and 4A′ state for the corresponding anion. Three types of the neutral-anion energy separations are the adiabatic electron affinity (EAad), the vertical electron affinity (EAvert), and the vertical detachment energy (VDE). The EAad value predicted by B3LYP method is 4.52 eV. The bond dissociation energies De (BrO4F → BrO4-mF + Om) (m = 1–4) and De− (BrO4F− → BrO4-mF− + Om and BrO4F− → BrO4-mF + Om−) are predicted. The adiabatic electron affinities (EAad) were predicted to be 4.52 eV for F-Br…O2…O2 (3A′←4A′) (B3LYP method). PMID:19742128
NASA Astrophysics Data System (ADS)
Wang, Lu; Li, Qing-Fang; Yang, Cui-Hong; Wei, Yue-Ling; Zhu, Xing-Feng; Rao, Wei-Feng
2016-05-01
The geometric and electronic structures of the ethylene oxide (EO) molecule adsorbed on Si(100)-(2 × 1) surface were investigated by using the density-functional theory calculations. All possible adsorbed structures were considered and it was found that only four adsorption structures are stable. The calculations of the formation energy revealed the most stable conformation and demonstrated that the nature of Si-O bond significantly affects the stability of adsorption systems. The analysis of corresponding electronic structures showed that two adsorbed structures are still semiconductor compounds but the other two are not. In particular, the EO after adsorbing was found to be connected via a ring-opening reaction where the molecule forms a five-membered ring together with the surface of dimer silicon atoms, and the produced five-membered ring is almost perpendicular to the silicon surface.
Structural and electronic properties of Sb{sub n} (n=2-10) clusters using density-functional theory
Zhou Xiaolan; Zhao Jijun; Chen Xiaoshuang; Lu Wei
2005-11-15
Using a density-functional calculation within the generalized gradient approximation, we have investigated the structural and electronic properties of antimony clusters. The lowest-energy structures of neutral and anionic Sb{sub n} clusters up to n=10 were determined and compared with those of light-element group-V clusters (As{sub n},P{sub n}). The size dependence of binding energies, gaps between lowest-unoccupied and highest-occupied molecular orbitals, second-order difference of total energies, and adiabatic electron affinities (EAs) of Sb{sub n} clusters were discussed. Even-odd oscillation as well as magic-number behavior at Sb{sub 4} were observed. Qualitative agreement was found between theoretical and experimental EAs.
Li, Xifeng; Sevilla, Michael D; Sanche, Léon
2003-11-12
The discovery of DNA strand breaks induced by low energy secondary electrons sparks a necessity to elucidate the mechanism. Through theoretical studies based on a sugar-phosphate-sugar model that mimics a backbone section of the DNA strand, it is found that bond cleavages at 3' or 5'C-O sites after addition of an electron are possible with a ca. 10 kcal/mol activation barrier. Moreover, the potential energy surfaces show that dissociation at both sites is highly favorable thermodynamically. Although the phosphate group in DNA is not a favored site for electron attachment because of competitive electron transfer to the bases, any electrons which attach to phosphates on first encounter may induce strand breaks even when the electron energy is near zero eV. These findings have profound implication as low energy secondary electrons are abundantly generated in all types of ionization radiation. PMID:14599198
Prociuk, Alexander; Van Kuiken, Ben; Dunietz, Barry D
2006-11-28
Electronic transmission through a metal-molecule-metal system is calculated by employing a Green's function formalism in the scattering based scheme. Self-energy models representing the bulk and the potential bias are used to describe electron transport through the molecular system. Different self-energies can be defined by varying the partition between device and bulk regions of the metal-molecule-metal model system. In addition, the self-energies are calculated with different representations of the bulk through its Green's function. In this work, the dependence of the calculated transmission on varying the self-energy subspaces is benchmarked. The calculated transmission is monitored with respect to the different choices defining the self-energy model. In this report, we focus on one-dimensional model systems with electronic structures calculated at the density functional level of theory. PMID:17144733
Kuwahara, Riichi; Tadokoro, Yoichi; Ohno, Kaoru
2014-08-28
In this paper, we calculate kinetic and potential energy contributions to the electronic ground-state total energy of several isolated atoms (He, Be, Ne, Mg, Ar, and Ca) by using the local density approximation (LDA) in density functional theory, the Hartree-Fock approximation (HFA), and the self-consistent GW approximation (GWA). To this end, we have implemented self-consistent HFA and GWA routines in our all-electron mixed basis code, TOMBO. We confirm that virial theorem is fairly well satisfied in all of these approximations, although the resulting eigenvalue of the highest occupied molecular orbital level, i.e., the negative of the ionization potential, is in excellent agreement only in the case of the GWA. We find that the wave function of the lowest unoccupied molecular orbital level of noble gas atoms is a resonating virtual bound state, and that of the GWA spreads wider than that of the LDA and thinner than that of the HFA. PMID:25173006
Nordlund, Dennis; Odelius, Michael; Bluhm, Hendrik; Ogasawara, Hirohito; Pettersson, Lars G.M.; Nilsson, Anders
2008-04-29
We present valence photoelectron emission spectra of liquid water in comparison with gas-phase water, ice close to the melting point, low temperature amorphous and crystalline ice. All aggregation states have major electronic structure changes relative to the free molecule, with rehybridization and development of bonding and anti-bonding states accompanying the hydrogen bond formation. Sensitivity to the local structural order, most prominent in the shape and splitting of the occupied 3a{sub 1} orbital, is understood from the electronic structure averaging over various geometrical structures, and reflects the local nature of the orbital interaction.
Density Functional Theory Studies of the Electronic Structure of Solid State Actinide Oxides
Wen, Xiaodong; Martin, Richard L.; Henderson, Thomas M.; Scuseria, Gustavo E.
2013-02-13
The actinide oxides have been extensively studied in the context of the nuclear fuel cycle. They are also of fundamental interest as members of a class of strongly correlated materials, the Mott insulators. Their complex physical and chemical properties make them challenging systems to characterize, both experimentally and theoretically. Chiefly, this is because actinide oxides can exhibit both electronic localization and electronic delocalization and have partially occupied f orbitals, which can lead to multiple possibilities for ground states. Of particular concern for theoretical work is that the large number of competing states display strong correlations which are dffcult to capture with computationally tractable methods.
2014-01-01
In systems in atomic scale and nanoscale such as clusters or agglomerates constituted by particles from a few to less than 100 atoms, quantum confinement effects are very important. Their optical and electronic properties are often dependent on the size of the systems and the way in which the atoms in these clusters are bonded. Generally, these nanostructures display optical and electronic properties significantly different to those found in corresponding bulk materials. Silicon agglomerates embedded in silicon rich oxide (SRO) films have optical properties, which have been reported to be directly dependent on silicon nanocrystal size. Furthermore, the room temperature photoluminescence (PL) of SRO has repeatedly generated a huge interest due to its possible applications in optoelectronic devices. However, a plausible emission mechanism has not been widely accepted in the scientific community. In this work, we present a short review about the experimental results on silicon nanoclusters in SRO considering different techniques of growth. We focus mainly on their size, Raman spectra, and photoluminescence spectra. With this as background, we employed the density functional theory with a functional B3LYP and a basis set 6-31G* to calculate the optical and electronic properties of clusters of silicon (constituted by 15 to 20 silicon atoms). With the theoretical calculation of the structural and optical properties of silicon clusters, it is possible to evaluate the contribution of silicon agglomerates in the luminescent emission mechanism, experimentally found in thin SRO films. PMID:25276105
Espinosa-Torres, Néstor D; la Luz, David Hernández-de; Flores-Gracia, José Francisco J; Luna-López, José A; Martínez-Juárez, Javier; Vázquez-Valerdi, Diana E
2014-01-01
In systems in atomic scale and nanoscale such as clusters or agglomerates constituted by particles from a few to less than 100 atoms, quantum confinement effects are very important. Their optical and electronic properties are often dependent on the size of the systems and the way in which the atoms in these clusters are bonded. Generally, these nanostructures display optical and electronic properties significantly different to those found in corresponding bulk materials. Silicon agglomerates embedded in silicon rich oxide (SRO) films have optical properties, which have been reported to be directly dependent on silicon nanocrystal size. Furthermore, the room temperature photoluminescence (PL) of SRO has repeatedly generated a huge interest due to its possible applications in optoelectronic devices. However, a plausible emission mechanism has not been widely accepted in the scientific community. In this work, we present a short review about the experimental results on silicon nanoclusters in SRO considering different techniques of growth. We focus mainly on their size, Raman spectra, and photoluminescence spectra. With this as background, we employed the density functional theory with a functional B3LYP and a basis set 6-31G* to calculate the optical and electronic properties of clusters of silicon (constituted by 15 to 20 silicon atoms). With the theoretical calculation of the structural and optical properties of silicon clusters, it is possible to evaluate the contribution of silicon agglomerates in the luminescent emission mechanism, experimentally found in thin SRO films. PMID:25276105
NASA Astrophysics Data System (ADS)
Lacombe, Lionel; Dinh, P. Huong Mai; Reinhard, Paul-Gerhard; Suraud, Eric; Sanche, Leon
2015-08-01
We present an extension of standard time-dependent density functional theory (TDDFT) to include the evaluation of rare reaction channels, taking as an example of application the theoretical modelling of electron attachment to molecules. The latter process is of great importance in radiation-induced damage of biological tissue for which dissociative electron attachment plays a decisive role. As the attachment probability is very low, it cannot be extracted from the TDDFT propagation whose mean field provides an average over various reaction channels. To extract rare events, we augment TDDFT by a perturbative treatment to account for the occasional jumps, namely electron capture in our test case. We apply the modelling to electron attachment to H2O, H3O+, and (H2O)2. Dynamical calculations have been done at low energy (3-16 eV). We explore, in particular, how core-excited states of the targets show up as resonances in the attachment probability. Contribution to the Topical Issue "COST Action Nano-IBCT: Nano-scale Processes Behind Ion-Beam Cancer Therapy", edited by Andrey Solov'yov, Nigel Mason, Gustavo García, Eugene Surdutovich.
Morzan, Uriel N.; Ramírez, Francisco F.; Scherlis, Damián A. E-mail: mcgl@qb.ffyb.uba.ar; Lebrero, Mariano C. González E-mail: mcgl@qb.ffyb.uba.ar
2014-04-28
This article presents a time dependent density functional theory (TDDFT) implementation to propagate the Kohn-Sham equations in real time, including the effects of a molecular environment through a Quantum-Mechanics Molecular-Mechanics (QM-MM) hamiltonian. The code delivers an all-electron description employing Gaussian basis functions, and incorporates the Amber force-field in the QM-MM treatment. The most expensive parts of the computation, comprising the commutators between the hamiltonian and the density matrix—required to propagate the electron dynamics—, and the evaluation of the exchange-correlation energy, were migrated to the CUDA platform to run on graphics processing units, which remarkably accelerates the performance of the code. The method was validated by reproducing linear-response TDDFT results for the absorption spectra of several molecular species. Two different schemes were tested to propagate the quantum dynamics: (i) a leap-frog Verlet algorithm, and (ii) the Magnus expansion to first-order. These two approaches were confronted, to find that the Magnus scheme is more efficient by a factor of six in small molecules. Interestingly, the presence of iron was found to seriously limitate the length of the integration time step, due to the high frequencies associated with the core-electrons. This highlights the importance of pseudopotentials to alleviate the cost of the propagation of the inner states when heavy nuclei are present. Finally, the methodology was applied to investigate the shifts induced by the chemical environment on the most intense UV absorption bands of two model systems of general relevance: the formamide molecule in water solution, and the carboxy-heme group in Flavohemoglobin. In both cases, shifts of several nanometers are observed, consistently with the available experimental data.
NASA Astrophysics Data System (ADS)
Morzan, Uriel N.; Ramírez, Francisco F.; Oviedo, M. Belén; Sánchez, Cristián G.; Scherlis, Damián A.; Lebrero, Mariano C. González
2014-04-01
This article presents a time dependent density functional theory (TDDFT) implementation to propagate the Kohn-Sham equations in real time, including the effects of a molecular environment through a Quantum-Mechanics Molecular-Mechanics (QM-MM) hamiltonian. The code delivers an all-electron description employing Gaussian basis functions, and incorporates the Amber force-field in the QM-MM treatment. The most expensive parts of the computation, comprising the commutators between the hamiltonian and the density matrix—required to propagate the electron dynamics—, and the evaluation of the exchange-correlation energy, were migrated to the CUDA platform to run on graphics processing units, which remarkably accelerates the performance of the code. The method was validated by reproducing linear-response TDDFT results for the absorption spectra of several molecular species. Two different schemes were tested to propagate the quantum dynamics: (i) a leap-frog Verlet algorithm, and (ii) the Magnus expansion to first-order. These two approaches were confronted, to find that the Magnus scheme is more efficient by a factor of six in small molecules. Interestingly, the presence of iron was found to seriously limitate the length of the integration time step, due to the high frequencies associated with the core-electrons. This highlights the importance of pseudopotentials to alleviate the cost of the propagation of the inner states when heavy nuclei are present. Finally, the methodology was applied to investigate the shifts induced by the chemical environment on the most intense UV absorption bands of two model systems of general relevance: the formamide molecule in water solution, and the carboxy-heme group in Flavohemoglobin. In both cases, shifts of several nanometers are observed, consistently with the available experimental data.
Huang, Ruiqi; Wang, Qingxia; Cai, Xiaolin; Li, Chong; Jia, Yu; Wang, Fei; Wang, Sanjun
2015-07-15
Employing the first-principles combined with hybrid functional calculations, the electronic and magnetic properties of GaAs doped with a N{sub 2} molecule are investigated in this work. We find that in Ga{sub 32}As{sub 31}(N{sub 2}){sub As} the N-N split is able to saturate the dangling bond of Ga atom ,form sp{sup 3}-like hybridization, and simultaneously supply an extra localized electron, leading to a magnetic ground state with a magnetic moment of ∼1μ{sub B}. This magnetic ground state is different from previously nonmagnetic results predicted by PBE functional, which results from the self-interaction error inherent in semi-local density functional theory. Moreover, the band gap of magnetic ground state of Ga{sub 32}As{sub 31}(N{sub 2}){sub As} alloy decreases, which is relative to GaAs . Finally we discuss and explain why the magnetism is not discovered in previous experiments and theories.
Panek, Jaroslaw; Latajaka, Zdzislaw
2000-12-26
Results of DFT and MP4 calculations on AlNO2 and GaNO2 molecules are presented. One Cs and two C2v structures (two minima and one TS) are found and their energies and vibrational frequencies are reported and discussed. The minima are close in energy and lie ca. 70 kcal mol-1 below reactants (M+NO2). More insight is obtained via topological analysis of electron density and electron localization function (ELF). It is shown that the molecules are bound mainly via electrostatic interactions, and there is a significant charge transfer from metal atom to the NO2 moiety. Detailed analysis of the ELF shows that the loss of stability of gallium complexes with respect to aluminium structures is best explained by (antibonding) influence of gallium semi-cored electrons.
NASA Astrophysics Data System (ADS)
Fortunato, Leandro F.; Zubieta, Carolina E.; Fuente, Silvia A.; Belelli, Patricia G.; Ferullo, Ricardo M.
2016-11-01
We report a density functional theory (DFT) investigation on the interaction of tiny Aun (n = 1-5) clusters with the bare and hydroxylated (110) surfaces of goethite (α-FeOOH). Both adsorption and atom-by-atom nucleation processes were modeled. The adsorption is shown to be strong on the bare surface and takes place preferentially through the interaction of Au atoms with unsaturated surface oxygen anions, accompanied with an electronic charge transfer from the metal to the support. Au3, Au4 and Au5 planar structures resulted to be particularly stable due to polarization effects; indeed, Coulombic repulsion between basal Au atoms and surface oxygen anions promotes the displacement of the electronic density toward terminal Au atoms producing a Au+δ(basal)/Au-δ(terminal) polarization. On the hydroxylated surface, Au clusters adsorb more weakly with respect to the bare surface, mainly through monocoordinated surface hydroxyl groups and tricoordinated oxygen ions. Concerning the nucleation mechanism, while on the hydroxylated surface the nucleation energy is governed by the spin of the interacting systems, on the bare surface polarization effects seems to play a predominant role.
Effects of semicore d-electrons in screened-exchange density functional methods
NASA Astrophysics Data System (ADS)
Lee, Byounghak; Wang, Lin-Wang
2007-03-01
We report a theoretical study on the role of shallow d states in the screened-exchange local density approximation (sX-LDA) band structure of binary semiconductor systems. We found that the inaccurate pseudo-wavefunctions can lead to 1) an overestimation of the screened-exchange interaction between the localized d states and the delocalized higher energy s and p states and 2) an underestimation of the screened-exchange interaction between the d states. The resulting sX-LDA band structures have substantially smaller band gaps compared with experiments. We correct the pseudo-wavefunctions of d states by including the s and p states of the same shell in the valence states. The correction of pseudo-wavefunctions yields band gaps and the d state binding energy in good agreement with experiments. Compared with the quasi-particle GW method, our sX-LDA results shows not only similar quality band gaps but also much better d state binding energy. As an example, we present sX-LDA results of s-d coupling in zinc-blende semiconductors and compare them with LDA+U results. We also present an efficient method to correct the pseudo-wavefunction exchange-integral error by using projection of wavefunctions onto atomic orbitals.
Effects of d-electrons in pseudopotential screened-exchange density functional calculations
Lee, Byounghak; Canning, Andrew; Wang, Lin-Wang
2007-09-12
We report a theoretical study on the role of shallow d states in the screened-exchange local density approximation (sX-LDA) band structure of binary semiconductor systems.We found that inaccurate pseudo-wavefunctions can lead to 1) an overestimation of the screened-exchange interaction betweenthe localized d states and the delocalized higher energy s and p states and 2) an underestimation of the screened-exchange interaction between the d states. The resulting sX-LDA band structures have substantially smaller band gaps compared with experiments. We correct the pseudo-wavefunctions of d states by including the semicore s and p states of the same shell in the valence states. The correction of pseudo-wavefunctions yields band gaps and d state binding energies in good agreement with experiments and the full potential linearized augmented plane wave sX-LDA calculations. Compared with the quasi-particle GW method, our sX-LDA results shows not only similar quality on the band gaps but also much better d state binding energies. Combined with its capability of ground state structure calculation, the sX-LDA is expected to be a valuable theoretical tool for the II-VI and III-V (especially the III-N) bulk semiconductors and nanostructure studies.
Effects of d-electrons in pseudopotential screened-exchange density functional calculations
Wang, Lin-Wang; Lee, Byounghak; Canning, Andrew; Wang, Lin-Wang
2008-08-11
We report a theoretical study on the role of shallow d states in the screened-exchange local density approximation (sX-LDA) band structure of binary semiconductor systems. We found that the inaccurate pseudo-wavefunctions can lead to (1) an overestimation of the screened-exchange interaction between the localized d states and the delocalized higher energy s and p states and (2) an underestimation of the screened-exchange interaction between the d states. The resulting sX-LDA band structures have substantially smaller band gaps compared with experiments. We correct the pseudo-wavefunctions of d states by including the semicore s and p states of the same shell in the valence states. The correction of pseudo-wavefunctions yields band gaps and the d state binding energy with good agreements with experiments and the full potential linearized augmented planewave (FLAPW) calculations. Compared with the quasi-particle GW method, our sX-LDA results shows not only similar quality on the band gaps but also much better d state binding energy. Combined with its capability of ground state structure calculation, the sX-LDA is expected to be a valuable theoretical tool for the II-VI and III-V (especially the III-N) bulk semiconductors and nanostructure studies.
Liu, Jie; Liang, WanZhen
2011-07-01
We present the analytical expression and computer implementation for the second-order energy derivatives of the electronic excited state with respect to the nuclear coordinates in the time-dependent density functional theory (TDDFT) with Gaussian atomic orbital basis sets. Here, the Tamm-Dancoff approximation to the full TDDFT is adopted, and therefore the formulation process of TDDFT excited-state Hessian is similar to that of configuration interaction singles (CIS) Hessian. However, due to the replacement of the Hartree-Fock exchange integrals in CIS with the exchange-correlation kernels in TDDFT, many quantitative changes in the derived equations are arisen. The replacement also causes additional technical difficulties associated with the calculation of a large number of multiple-order functional derivatives with respect to the density variables and the nuclear coordinates. Numerical tests on a set of test molecules are performed. The simulated excited-state vibrational frequencies by the analytical Hessian approach are compared with those computed by CIS and the finite-difference method. It is found that the analytical Hessian method is superior to the finite-difference method in terms of the computational accuracy and efficiency. The numerical differentiation can be difficult due to root flipping for excited states that are close in energy. TDDFT yields more exact excited-state vibrational frequencies than CIS, which usually overestimates the values. PMID:21744894
Ahangari, Morteza Ghorbanzadeh; Fereidoon, Abdolhosein; Ganji, Masoud Darvish
2013-08-01
We performed first principles calculations based on density functional theory (DFT) to investigate the effect of epoxy monomer content on the electronic and mechanical properties of single-walled carbon nanotubes (SWCNTs). Our calculation results reveal that interfacial interaction increases with increasing numbers of epoxy monomers on the surface of SWCNTs. Furthermore, density of states (DOS) results showed no orbital hybridization between the epoxy monomers and nanotubes. Mulliken charge analysis shows that the epoxy polymer carries a positive charge that is directly proportional to the number of monomers. The Young's modulus of the nanotubes was also studied as a function of monomer content. It was found that, with increasing number of monomers on the nanotubes, the Young's modulus first decreases and then approaches a constant value. The results of a SWCNT pullout simulation suggest that the interfacial shear stress of the epoxy/SWCNT complex is approximately 68 MPa. These results agreed well with experimental results, thus proving that the simulation methods used in this study are viable. PMID:23609226
Proynov, Emil; Liu, Fenglai; Shao, Yihan; Kong, Jing
2012-01-01
In a recent letter [E. Proynov, Y. Shao, and J. Kong, Chem. Phys. Lett. 493, 381 (2010)10.1016/j.cplett.2010.05.029], Becke's B05 model of nondynamic electron correlation in density functional theory was implemented self-consistently with computational efficiency (the “SCF-RI-B05” scheme). Important modifications of the algorithm were done in order to make the self-consistency feasible. In the present work, we give a complete account of the SCF-RI-B05 algorithm, including all the formulae for the analytical representation of the B05 functional and for its self-consistent field (SCF) potential. The average performance of the SCF-RI-B05 method reported in the above letter was somewhat less accurate, compared to the original B05 implementation, mainly because the parameters of the original B05 model were optimized with post-local-spin-density calculations. In this work, we report improved atomization energies with SCF-RI-B05, based on a SCF re-optimization of its four linear parameters. The re-optimized SCF-RI-B05 scheme is validated also on reaction barriers, and on the subtle energetics of NO dimer, an exemplary system of strong nondynamic correlation. It yields both the binding energy and the singlet-triplet splitting of the NO dimer correctly, and close to the benchmarks reported in the literature. PMID:22280739
All-electron Kohn–Sham density functional theory on hierarchic finite element spaces
Schauer, Volker; Linder, Christian
2013-10-01
In this work, a real space formulation of the Kohn–Sham equations is developed, making use of the hierarchy of finite element spaces from different polynomial order. The focus is laid on all-electron calculations, having the highest requirement onto the basis set, which must be able to represent the orthogonal eigenfunctions as well as the electrostatic potential. A careful numerical analysis is performed, which points out the numerical intricacies originating from the singularity of the nuclei and the necessity for approximations in the numerical setting, with the ambition to enable solutions within a predefined accuracy. In this context the influence of counter-charges in the Poisson equation, the requirement of a finite domain size, numerical quadratures and the mesh refinement are examined as well as the representation of the electrostatic potential in a high order finite element space. The performance and accuracy of the method is demonstrated in computations on noble gases. In addition the finite element basis proves its flexibility in the calculation of the bond-length as well as the dipole moment of the carbon monoxide molecule.
Yang, Zhi; Xiong, Shi-Jie
2008-09-28
The geometries stability, electronic properties, and magnetism of Y(n)O clusters up to n=14 are systematically studied with density functional theory. In the lowest-energy structures of Y(n)O clusters, the equilibrium site of the oxygen atom gradually moves from an outer site of the cluster, via a surface site, and finally, to an interior site as the number of the Y atoms increases from 2 to 14. Starting from n=12, the O atom falls into the center of the cluster with the Y atoms forming the outer frame. The results show that clusters with n=2, 4, 8, and 12 are more stable than their respective neighbors, and that the total magnetic moments of Y(n)O clusters are all quite small except Y(12)O cluster. The lowest-energy structure of Y(12)O cluster is a perfect icosahedron with a large magnetic moment 6mu(B). In addition, we find that the total magnetic moments are quenched for n=2, 6, and 8 due to the closed-shell electronic configuration. The calculated ionization potentials and electron affinities are in good agreement with the experimental results, which imply that the present theoretical treatments are satisfactory. PMID:19045024
NASA Astrophysics Data System (ADS)
Zhong, X.; Rungger, I.; Zapol, P.; Heinonen, O.
2015-03-01
Understanding electronic properties of substoichiometric phases of titanium oxide such as Magnéli phase T i4O7 is crucial in designing and modeling resistive switching devices. Here we present our study on Magnéli phase T i4O7 together with rutile Ti O2 and T i2O3 using density functional theory methods with atomic-orbital-based self-interaction correction (ASIC). We predict a new antiferromagnetic (AF) ground state in the low temperature (LT) phase, and we explain energy difference with a competing AF state using a Heisenberg model. The predicted energy ordering of these states in the LT phase is calculated to be robust in a wide range of modeled isotropic strain. We have also investigated the dependence of the electronic structures of the Ti-O phases on stoichiometry. The splitting of titanium t2 g orbitals is enhanced with increasing oxygen deficiency as Ti-O is reduced. The electronic properties of all these phases can be reasonably well described by applying ASIC with a "standard" value for transition metal oxides of the empirical parameter α of 0.5 representing the magnitude of the applied self-interaction correction.
Regeta, Khrystyna; Bannwarth, Christoph; Grimme, Stefan; Allan, Michael
2015-06-28
The technique of low energy (0-30 eV) electron impact spectroscopy, originally developed for gas phase molecules, is applied to room temperature ionic liquids (IL). Electron energy loss (EEL) spectra recorded near threshold, by collecting 0-2 eV electrons, are largely continuous, assigned to excitation of a quasi-continuum of high overtones and combination vibrations of low-frequency modes. EEL spectra recorded by collecting 10 eV electrons show predominantly discrete vibrational and electronic bands. The vibrational energy-loss spectra correspond well to IR spectra except for a broadening (∼0.04 eV) caused by the liquid surroundings, and enhanced overtone activity indicating a contribution from resonant excitation mechanism. The spectra of four representative ILs were recorded in the energy range of electronic excitations and compared to density functional theory multireference configuration interaction (DFT/MRCI) calculations, with good agreement. The spectra up to about 8 eV are dominated by π-π* transitions of the aromatic cations. The lowest bands were identified as triplet states. The spectral region 2-8 eV was empty in the case of a cation without π orbitals. The EEL spectrum of a saturated solution of methylene green in an IL band showed the methylene green EEL band at 2 eV, indicating that ILs may be used as a host to study nonvolatile compounds by this technique in the future. PMID:26018044
NASA Astrophysics Data System (ADS)
Hahn, Torsten; Rückerl, Florian; Liebing, Simon; Pederson, Mark
We present our experimental and theoretical results on novel Picene/F4TCNQ and Manganese-Phthalocyanine/F4TCNQ donor / acceptor systems. We apply the recently developed Fermi-orbital based approach for self-interaction corrected density functional theory (FO-SIC DFT) to these materials and compare the results to standard DFT calculations and to experimental data obtained by photoemission spectroscopy. We focus our analysis on the description of the magnitude of the ground state charge transfer and on the details of the formed hybrid orbitals. Further, we show that for weakly bound donor / acceptor systems the FO-SIC approach delivers a more realistic description of the electronic structure compared to standard DFT calculations Support by DFG FOR1154 is greatly acknowledged.
NASA Astrophysics Data System (ADS)
Huang, Y.-T.; Diehl, R. D.; Pulkkinen, A.; Pussi, K.
2015-09-01
Alkali metal adsorption systems provide important models for chemisorption. Low-energy electron diffraction experiments and density functional theory calculations were carried out for the adsorption of potassium on Pb(1 0 0). The stable structure for all submonolayer coverages was found to be the commensurate c(2 × 2) structure, with potassium atoms located in substitutional sites in the top substrate layer. This structure is temperature activated and occurs for adsorption or annealing of the film above 200 K. This finding is consistent with an earlier theory that proposed that for substrates with low energies of vacancy formation, substitutional structures can be the most stable. The structural and vibrational parameters deduced from the experiment are in agreement with the calculated values, and these values fit well into and add to the database of alkali metal adsorption properties.
NASA Astrophysics Data System (ADS)
Zhang, Xiao; Guo, Ling; Guo, Jian; Ren, Ning-Ning
2014-03-01
The binding energy, dissociation energy, ionization potentials, electron affinities, gap and stability of small Al n Pt ( n = 1-15) clusters, in comparison with pure aluminum clusters have been systematically investigated by means of density functional calculations at the B3LYP level. The growth patten for Al n Pt clusters is that the Pt atom substituted the surface atom of the Al n + 1 clusters for n < 13. Starting from n = 13, the Pt atom completely falls into the center of the Al-frame. The Pt atom substituted the center atom of the Al n + 1 clusters to form the Pt-encapsulated Aln geometries for n > 13. We also find that the impurity Pt atom causes local structural distortion due to different atomic radii and different bonding characteristics. The clusters with total atom numbers of 2, 7, and 11 exhibit high stability.
Huang, Y-T; Diehl, R D; Pulkkinen, A; Pussi, K
2015-09-01
Alkali metal adsorption systems provide important models for chemisorption. Low-energy electron diffraction experiments and density functional theory calculations were carried out for the adsorption of potassium on Pb(1 0 0). The stable structure for all submonolayer coverages was found to be the commensurate c(2 × 2) structure, with potassium atoms located in substitutional sites in the top substrate layer. This structure is temperature activated and occurs for adsorption or annealing of the film above 200 K. This finding is consistent with an earlier theory that proposed that for substrates with low energies of vacancy formation, substitutional structures can be the most stable. The structural and vibrational parameters deduced from the experiment are in agreement with the calculated values, and these values fit well into and add to the database of alkali metal adsorption properties. PMID:26202693
Density-functional study of the electronic and optical properties of the spinel compound CuIr2S4
NASA Astrophysics Data System (ADS)
Sarkar, Soumyajit; de Raychaudhury, Molly; Saha-Dasgupta, T.
2009-03-01
Using first-principles density-functional calculations we have computed the electronic and optical properties of spinel compound CuIr2S4 , which undergoes a structural phase transition, accompanied by a metal-insulator transition at a temperature of about 230 K. The nature of this transition has been discussed in the literature in terms of both the correlated singlet formation picture as well as the orbitally driven Peierls transition picture. Our first-principles calculations find little role of correlation. Our calculated reflectivity and conductivity data for both the high-temperature and low-temperature phases are found to be in good agreement with measured data by Wang [Phys. Rev. B 69, 153104 (2004)].
NASA Astrophysics Data System (ADS)
Zheng, Y.; Neville, J. J.; Brion, C. E.; Wang, Y.; Davidson, E. R.
1994-11-01
The binding energy spectra and momentum distributions of all valence orbitals of acetone have been studied by electron momentum spectroscopy (EMS) and SCF, MRSD-CI, and density functional theory (DFT) calculations. The experiment was performed using a multichannel EMS spectrometer at a total energy of 1200 eV. Binding energy spectra measured in the energy range of 6-60 eV are compared with the results of OVGF and 2ph-TDA many-body Green's function calculations. In the inner valence region strong splitting of the 5a 1 and 4a 1 orbitals due to final state electron correlation is observed. The distribution of energies and pole strengths predicted by the Green's function calculations deviates considerably from the measured ionization energies and strengths in the innervalence region. The measured momentum distributions are compared with calculations at the level of the target Hartree-Fock approximation (THFA) using the SCF method and the target Kohn-Sham approximation (TKSA) using DFT and the local-density approximation. Basis sets used for the SCF calculations ranged from the simplest (STO-3G) to large (204-GTO) and for the DFT calculations very large atomic natural orbital (ANO) basis sets were used. The effects of electron correlation and relaxation are also investigated in MRSD-CI calculations of the full ion-neutral overlap amplitude using large and saturated basis sets. In general, the THFA model with an intermediate basis set and very diffuse functions (6-311 + +G ∗∗) and with a near Hartree-Fock limit SCF wavefunction (204-GTO), and the TKSA-DFT model with an ANO basis set all provide reasonable predictions of momentum distributions for most orbitals. However, none of these calculations gives a completely satisfactory description of the momentum distribution of the HOMO (5b 2) orbital.
The mapping of electronic energy distributions using experimental electron density.
Tsirelson, Vladimir G
2002-08-01
It is demonstrated that the approximate kinetic energy density calculated using the second-order gradient expansion with parameters of the multipole model fitted to experimental structure factors reproduces the main features of this quantity in a molecular or crystal position space. The use of the local virial theorem provides an appropriate derivation of approximate potential energy density and electronic energy density from the experimental (model) electron density and its derivatives. Consideration of these functions is not restricted by the critical points in the electron density and provides a comprehensive characterization of bonding in molecules and crystals. PMID:12149553
Gulans, Andris; Kontur, Stefan; Meisenbichler, Christian; Nabok, Dmitrii; Pavone, Pasquale; Rigamonti, Santiago; Sagmeister, Stephan; Werner, Ute; Draxl, Claudia
2014-09-10
Linearized augmented planewave methods are known as the most precise numerical schemes for solving the Kohn-Sham equations of density-functional theory (DFT). In this review, we describe how this method is realized in the all-electron full-potential computer package, exciting. We emphasize the variety of different related basis sets, subsumed as (linearized) augmented planewave plus local orbital methods, discussing their pros and cons and we show that extremely high accuracy (microhartrees) can be achieved if the basis is chosen carefully. As the name of the code suggests, exciting is not restricted to ground-state calculations, but has a major focus on excited-state properties. It includes time-dependent DFT in the linear-response regime with various static and dynamical exchange-correlation kernels. These are preferably used to compute optical and electron-loss spectra for metals, molecules and semiconductors with weak electron-hole interactions. exciting makes use of many-body perturbation theory for charged and neutral excitations. To obtain the quasi-particle band structure, the GW approach is implemented in the single-shot approximation, known as G(0)W(0). Optical absorption spectra for valence and core excitations are handled by the solution of the Bethe-Salpeter equation, which allows for the description of strongly bound excitons. Besides these aspects concerning methodology, we demonstrate the broad range of possible applications by prototypical examples, comprising elastic properties, phonons, thermal-expansion coefficients, dielectric tensors and loss functions, magneto-optical Kerr effect, core-level spectra and more. PMID:25135665
Bykov, Dmytro; Neese, Frank
2015-10-01
In this Forum Article, an extensive discussion of the mechanism of six-electron, seven-proton nitrite reduction by the cytochrome c nitrite reductase enzyme is presented. On the basis of previous studies, the entire mechanism is summarized and a unified picture of the most plausible sequence of elementary steps is presented. According to this scheme, the mechanism can be divided into five functional stages. The first phase of the reaction consists of substrate binding and N-O bond cleavage. Here His277 plays a crucial role as a proton donor. In this step, the N-O bond is cleaved heterolytically through double protonation of the substrate. The second phase of the mechanism consists of two proton-coupled electron-transfer events, leading to an HNO intermediate. The third phase involves the formation of hydroxylamine, where Arg114 provides the necessary proton for the reaction. The second N-O bond is cleaved in the fourth phase of the mechanism, again triggered by proton transfer from His277. The Tyr218 side chain governs the fifth and last phase of the mechanism. It consists of radical transfer and ammonia formation. Thus, this mechanism implies that all conserved active-site side chains work in a concerted way in order to achieve this complex chemical transformation from nitrite to ammonia. The Forum Article also provides a detailed discussion of the density functional theory based cluster model approach to bioinorganic reactivity. A variety of questions are considered: the resting state of enzyme and substrate binding modes, interaction with the metal site and with active-site side chains, electron- and proton-transfer events, substrate dissociation, etc. PMID:26237518
π vs σ-Radical States of One-Electron Oxidized DNA/RNA Bases: A Density Functional Theory Study
Kumar, Anil; Sevilla, Michael D.
2013-01-01
As a result of their inherent planarity, DNA base radicals generated by one electron oxidation/reduction or bond cleavage form π- or σ-radicals. While most DNA base systems form π-radicals there are a number of nucleobase analogs such as one-electron oxidized 6-azauraci1, 6-azacytosine, and 2-thiothymine or one-electron reduced 5-bromouracil that form more reactive σ-radicals. Elucidating the availability of these states within DNA, base radical electronic structure is important to the understanding of the reactivity of DNA base radicals in different environments. In this work, we address this question by the calculation of the relative energies of π- and σ-radical states in DNA/RNA bases and their analogs. We used density functional theory B3LYP/6-31++G** method to optimize the geometries of π- and σ-radicals in Cs symmetry (i.e., planar) in the gas phase and in solution using the polarized continuum model (PCM). The calculations predict that σ- and π-radical states in one electron oxidized bases of thymine, T(N3-H)•, and uracil, U(N3-H)• are very close in energy, i.e., the π-radical is only ca. 4 kcal/mol more stable than the σ-radical. For the one electron oxidized radicals of cytosine, C•+, C(N4-H)•, adenine, A•+, A(N6-H)•, and guanine, G•+, G(N2-H)•, G(N1-H)• the π-radicals are ca. 16 to 41 kcal/mol more stable than their corresponding σ-radicals. Inclusion of solvent (PCM) is found to stabilize the π- over σ-radical of each of the systems. U(N3-H)• with three discrete water molecules in the gas phase, is found to form a three-electron σ bond between N3 atom of uracil and O atom of a water molecule but on inclusion of full solvation and discrete hydration the π-radical remains most stable.. PMID:24000793
NASA Astrophysics Data System (ADS)
Goings, Joshua J.; Li, Xiaosong
2016-06-01
One of the challenges of interpreting electronic circular dichroism (ECD) band spectra is that different states may have different rotatory strength signs, determined by their absolute configuration. If the states are closely spaced and opposite in sign, observed transitions may be washed out by nearby states, unlike absorption spectra where transitions are always positive additive. To accurately compute ECD bands, it is necessary to compute a large number of excited states, which may be prohibitively costly if one uses the linear-response time-dependent density functional theory (TDDFT) framework. Here we implement a real-time, atomic-orbital based TDDFT method for computing the entire ECD spectrum simultaneously. The method is advantageous for large systems with a high density of states. In contrast to previous implementations based on real-space grids, the method is variational, independent of nuclear orientation, and does not rely on pseudopotential approximations, making it suitable for computation of chiroptical properties well into the X-ray regime.
NASA Astrophysics Data System (ADS)
Chae, Kisung; Kim, Hanchul
2013-02-01
We perform spin-polarized density functional theory calculations for a Zn vacancy on the ZnO(10bar 10) surface. Two stable configurations of the surface Zn vacancy are found, and the activation energy barrier is estimated to be ˜0.01 eV. The lower energy configuration has a newly formed surface Zn-O bond to restore the bulk-like structure on the surface. Due to the newly formed bond, the vacancy state in the band gap is characterized by a complicated hybridization of neighboring surface and subsurface atoms and by a more extended electron density. Despite such a hybridization, the surface Zn vacancy is found to have a robust magnetic moment of 1 μ B , implying that surface Zn vacancies may be responsible for the ferromagnetism observed in ZnO thin films and nanoparticles. Simulated scanning tunneling microscope images show that the two structures of the surface Zn vacancy can be distinguished in the filled-state images.
NASA Astrophysics Data System (ADS)
Govindasamy, P.; Gunasekaran, S.; Ramkumaar, G. R.
2014-09-01
The Fourier transform infrared (FT-IR) and FT-Raman spectra of N-(4-hydroxy phenyl) acetamide (N4HPA) of painkiller agent were recorded in the region 4000-450 cm-1 and 4000-50 cm-1 respectively. Density functional theory (DFT) has been used to calculate the optimized geometrical parameter, atomic charges, and vibrational wavenumbers and intensity of the vibrational bands. The computed vibrational wave numbers were compared with the FT-IR and FT-Raman experimental data. The computational calculations at DFT/B3LYP level with 6-31G(d,p), 6-31++G(d,p), 6-311G(d,p) and 6-311++G(d,p) basis sets. The complete vibrational assignments were performed on the basis of the potential energy distribution (PED) of the vibrational modes calculated using Vibrational energy distribution analysis (VEDA 4) program. The oscillator’s strength calculated by TD-DFT and N4HPA is approach complement with the experimental findings. The NMR chemical shifts 13C and 1H were recorded and calculated using the gauge independent atomic orbital (GIAO) method. The molecular electrostatic potential (MESP) and electron density surfaces of the molecule were constructed. The Natural charges and intermolecular contacts have been interpreted using Natural Bond orbital (NBO) analysis the HOMO-LUMO energy gap has been calculated. The thermodynamic properties like entropy, heat capacity and zero vibrational energy have been calculated.
NASA Astrophysics Data System (ADS)
Rathi, Somilkumar J.; Ray, Asok K.
2008-08-01
Ab initio calculations within the framework of hybrid density functional theory and the finite cluster approximation have been performed for the electronic and geometric structures of three different types of armchair germanium carbide nanotube, from (3, 3) to (11, 11). Full geometry and spin optimizations with unrestricted symmetry have been performed. Physically pertinent quantities of interest such as the cohesive energies, band gaps, radial buckling, density of states, dipole moments, and Mulliken charge distributions have been investigated in detail for all nanotubes. For type I nanotubes, the largest cohesive energy obtained is 4.092 eV/atom, whereas for type II and type III nanotubes, the values are 3.987 eV/atom and 3.968 eV/atom, respectively. For optimized type I nanotubes, Ge atoms moved toward the tube axis and C atoms moved in the opposite direction after relaxation, opposite to the trends observed in types II and III. The band gaps for type I nanotubes are larger than the bulk 3C-GeC gap, varying between 2.666 and 3.016 eV, while type II and type III nanotubes have significantly lower band gaps, with all nanotubes being semiconducting in nature. Mulliken charge analysis indicates primarily ionic behavior for type I GeC nanotubes and a mixed ionic with covalent behavior for the other two types. None of the tubes appear to be magnetic. Applications in the field of nano-optoelectronic devices, molecular electronics, and band gap engineering are envisioned for GeC nanotubes.
Rathi, Somilkumar J; Ray, Asok K
2008-08-20
Ab initio calculations within the framework of hybrid density functional theory and the finite cluster approximation have been performed for the electronic and geometric structures of three different types of armchair germanium carbide nanotube, from (3, 3) to (11, 11). Full geometry and spin optimizations with unrestricted symmetry have been performed. Physically pertinent quantities of interest such as the cohesive energies, band gaps, radial buckling, density of states, dipole moments, and Mulliken charge distributions have been investigated in detail for all nanotubes. For type I nanotubes, the largest cohesive energy obtained is 4.092 eV/atom, whereas for type II and type III nanotubes, the values are 3.987 eV/atom and 3.968 eV/atom, respectively. For optimized type I nanotubes, Ge atoms moved toward the tube axis and C atoms moved in the opposite direction after relaxation, opposite to the trends observed in types II and III. The band gaps for type I nanotubes are larger than the bulk 3C-GeC gap, varying between 2.666 and 3.016 eV, while type II and type III nanotubes have significantly lower band gaps, with all nanotubes being semiconducting in nature. Mulliken charge analysis indicates primarily ionic behavior for type I GeC nanotubes and a mixed ionic with covalent behavior for the other two types. None of the tubes appear to be magnetic. Applications in the field of nano-optoelectronic devices, molecular electronics, and band gap engineering are envisioned for GeC nanotubes. PMID:21730632
Density functional study of structural and electronic properties of Al{sub n}@C{sub 60}
Dhiman, Shobhna; Kumar, Ranjan; Dharamvir, Keya
2014-04-24
Fullerene derivatives have been shown to make contributions in many types of applications. Ab initio investigation of structural and electronic properties of aluminum doped endohedral fullerene has been performed using numerical atomic orbital density functional theory. We have obtained ground state structures for Al{sub n}@C{sub 60} (n=1–10). Which shows that C{sub 60} molecule can accommodate maximum of nine aluminum atoms, for n > 9 the cage eventually break. Encapsulated large number of aluminum atoms leads to deformation of cage with diameter varies from 7.16Å to 7.95Å. Binding energy/Al atom is found to increase till n = 4 and after that it decreases with the number of Al atoms with a sudden increase for n=10 due to breakage of C{sub 60} cage and electronic affinity first increases till n=4 then it decreases up to n=9 with a sharp increase for n=10. Ionization potential also first increases and then decreases. Homo-Lumo gap decreases till n=3 with a sharp increase for n=4, after that it shows an oscillatory nature. The results obtained are consistent with available theoretical and experimental results. The ab-initio calculations were performed using SIESTA code with generalized gradient approximation (GGA)
NASA Astrophysics Data System (ADS)
Sun, Xu; Gu, You-song; Wang, Xue-qiang; Zhang, Yue
2012-06-01
The electronic properties and stability of Li-doped ZnO with various defects have been studied by calculating the electronic structures and defect formation energies via first-principles calculations using hybrid Hartree-Fock and density functional methods. The results from formation energy calculations show that Li pair complexes have the lowest formation energy in most circumstances and they consume most of the Li content in Li doped ZnO, which make the p-type conductance hard to obtain. The formation of Li pair complexes is the main obstacle to realize p-type conductance in Li doped ZnO. However, the formation energy of LiZn decreases as environment changes from Zn-rich to O-rich and becomes more stable than that of Li-pair complexes at highly O-rich environment. Therefore, p-type conductance can be obtained by Li doped ZnO grown or post annealed in oxygen rich atmosphere.
Walsh, Aron; Wei, S.-H.; Yan Yanfa; Al-Jassim, M. M.; Turner, John A.; Woodhouse, Michael; Parkinson, B. A.
2007-10-15
A systematic study of nine binary and ternary spinel oxides formed from Co, Al, and Fe is presented by means of density functional theory. Analysis of the structural, magnetic, and electronic properties through the series of materials is carried out. Preference for the octahedral spinel sites are found in the order Fe
A finite-temperature density functional study of electron self-trapping in 3He and 4He.
Jin, Dafei; Guo, Wei
2012-06-28
We introduce a compact finite-temperature density functional model to study electron self-trapping in both liquid and vapor (3)He and (4)He. This model can quantitatively reproduce the most essential thermodynamic properties of (3)He and (4)He along their liquid-vapor coexistence lines. The structures and energetics of self-trapped electron bubbles on the 1S ground state and 1P excited state are particularly investigated. Our results show that 1S and 1P bubbles exist in liquid at any temperature, whereas 1S bubbles exist in vapor only above 1.6 K in (3)He and above 2.8 K in (4)He, 1P bubbles exist in vapor only above 2.5 K in (3)He and 4.0 K in (4)He. An initially spherical 1P bubble is unstable against deformation towards a peanut shape. In liquid, a peanut-shaped 1P bubble is held from fission by surface tension until reaching the liquid-vapor critical point, whereas in vapor it always splits into two smaller bubbles. The existence of 1P bubbles in finite-temperature liquid helium and their fission instability in helium vapor reveal interesting physics in this system. PMID:22755590
Gedanken densities and exact constraints in density functional theory
Perdew, John P.; Ruzsinszky, Adrienn; Sun, Jianwei; Burke, Kieron
2014-05-14
Approximations to the exact density functional for the exchange-correlation energy of a many-electron ground state can be constructed by satisfying constraints that are universal, i.e., valid for all electron densities. Gedanken densities are designed for the purpose of this construction, but need not be realistic. The uniform electron gas is an old gedanken density. Here, we propose a spherical two-electron gedanken density in which the dimensionless density gradient can be an arbitrary positive constant wherever the density is non-zero. The Lieb-Oxford lower bound on the exchange energy can be satisfied within a generalized gradient approximation (GGA) by bounding its enhancement factor or simplest GGA exchange-energy density. This enhancement-factor bound is well known to be sufficient, but our gedanken density shows that it is also necessary. The conventional exact exchange-energy density satisfies no such local bound, but energy densities are not unique, and the simplest GGA exchange-energy density is not an approximation to it. We further derive a strongly and optimally tightened bound on the exchange enhancement factor of a two-electron density, which is satisfied by the local density approximation but is violated by all published GGA's or meta-GGA’s. Finally, some consequences of the non-uniform density-scaling behavior for the asymptotics of the exchange enhancement factor of a GGA or meta-GGA are given.
Eljarrat, Alberto; Sastre, Xavier; Peiró, Francesca; Estradé, Sónia
2016-06-01
In the present work, the dielectric response of III-nitride semiconductors is studied using density functional theory (DFT) band structure calculations. The aim of this study is to improve our understanding of the features in the low-loss electron energy-loss spectra of ternary alloys, but the results are also relevant to optical and UV spectroscopy results. In addition, the dependence of the most remarkable features with composition is tested, i.e. applying Vegard's law to band gap and plasmon energy. For this purpose, three wurtzite ternary alloys, from the combination of binaries AlN, GaN, and InN, were simulated through a wide compositional range (i.e., Al x Ga1-x N, In x Al1-x N, and In x Ga1-x N, with x=[0,1]). For this DFT calculations, the standard tools found in Wien2k software were used. In order to improve the band structure description of these semiconductor compounds, the modified Becke-Johnson exchange-correlation potential was also used. Results from these calculations are presented, including band structure, density of states, and complex dielectric function for the whole compositional range. Larger, closer to experimental values, band gap energies are predicted using the novel potential, when compared with standard generalized gradient approximation. Moreover, a detailed analysis of the collective excitation features in the dielectric response reveals their compositional dependence, which sometimes departs from a linear behavior (bowing). Finally, an advantageous method for measuring the plasmon energy dependence from these calculations is explained. PMID:26868876
Pinjari, Rahul V; Gejji, Shridhar P
2009-02-19
Inverted cucurbit[n]uril (i(x)CB[n], x = 1, 2; n = 6-8), the enantiomers of cucurbit[n]uril (CB[n]) comprising one or more inverted glycouril units, show distinct selectivity in recognition toward the guest by the virtue of shape and dimensions of its cavity. The iCB[n] (x = 1 and n = 6, 7) are isolated as intermediates during the synthesis of CB[n]. In this work, density functional theory using the hybrid B3LYP functional has been employed to derive the electronic structure and the NMR chemical shifts in the i(x)CB[n] hosts. The present calculations have shown that the inversion of the glycouril unit of CB[6] and CB[7] engenders a destabilization by 4.2 and 5.7 kJ mol(-1), respectively, and, as opposed to this, the iCB[8] is favored by 18.6 kJ mol(-1) over the corresponding CB[8] host. Likewise, i2CB[7] possessing two inverted glycourils are highly destabilized over CB[7]. A large separation of the inverted glycouril units reduces the repulsion between methine protons inside the cavity, rendering the 1,4-i2CB[n] (n = 7 or 8) to be of lowest energy. Stabilization energies from the self-consistent reaction field (SCRF) theory are calculated with water, ethanol, and tetrahydrofuran (THF) as solvents. Unlike in gas phase and other solvents, the stabilization hierarchy iCB[6] < iCB[7] < iCB[8] has been predicted in THF. Molecular electrostatic potential (MESP) was used to gauge the cavity shape of these hosts. Consequently the iCB[6] reveals a half-sprocket-like cavity; an additional tooth for each glycouril in the succeeding iCB[n] homologue was noticed. In the case of the 1,5-i2CB[8] enantiomer, the cavity turns out to be rectangular. The deeper MESP minima near the ureido oxygens suggest strong electrostatic interactions with the guest at the iCB[6] portals. The electron-rich region within the cavity explains the large affinity of CB[n] toward the electron deficient guests. The electronic distribution and shape and size of the cavity thus derived provide insights
NASA Astrophysics Data System (ADS)
Suhai, Sándor
1995-06-01
Structural and energetic aspects of the Peierls-type lattice dimerization were investigated in infinite, one-dimensional, periodic trans-polyacetylene (t-PA) using many-body perturbation theory (MBPT) and density-functional theory (DFT). Cohesive properties and dimerization parameters were obtained first for the classical Coulomb potential in the Hartree approximation and then by gradually turning on exchange and correlation potentials. Besides the nonlocal Hartree-Fock exchange, several other exchange functionals were used incorporating gradient corrections as well. For MBPT, electron correlation was included up to the fourth order of the Mo/ller-Plesset scheme and the behavior of lattice sums for different PT terms was analyzed in detail. The electrostatic part of the infinite lattice sums was computed by the multipole expansion technique. In solving the polymer Kohn-Sham equations, the performance of several different correlation potentials was studied again including different gradient corrections. Atomic basis sets of systematically increasing size, in the range of double-zeta to triple-zeta (TZ) up to TZ (3df,3p2d), were used in all calculations to construct the symmetry-adapted (Bloch-type) polymer wave functions, to fully optimize the structures, and to extrapolate different physical quantities to the limit of a hypothetical infinite basis set. Comparison of the different DFT results with MBPT and with experiments demonstrated the importance of gradient terms both for exchange and correlation. On the other hand, the best DFT functional, using a medium-size atomic basis set, excellently reproduced the cohesive and dimerization energies obtained for infinite t-PA at the MP4/TZ(3d2f,3p2d) level and provided dimerization parameters close to experiment. The experimentally observed lattice spacing of 2.46+/-0.01 Å will be correctly predicted both at the MBPT and DFT levels with 2.48 and 2.44 Å, respectively.
Ching, W. Y.; Aryal, Sitram; Rulis, Paul; Schnick, Wolfgang
2011-04-15
Using density-functional-theory-based ab initio methods, the electronic structure and physical properties of the newly synthesized nitride BeP{sub 2}N{sub 4} with a phenakite-type structure and the predicted high-pressure spinel phase of BeP{sub 2}N{sub 4} are studied in detail. It is shown that both polymorphs are wide band-gap semiconductors with relatively small electron effective masses at the conduction-band minima. The spinel-type phase is more covalently bonded due to the increased number of P-N bonds for P at the octahedral sites. Calculations of mechanical properties indicate that the spinel-type polymorph is a promising superhard material with notably large bulk, shear, and Young's moduli. Also calculated are the Be K, P K, P L{sub 3}, and N K edges of the electron energy-loss near-edge structure for both phases. They show marked differences because of the different local environments of the atoms in the two crystalline polymorphs. These differences will be very useful for the experimental identification of the products of high-pressure syntheses targeting the predicted spinel-type phase of BeP{sub 2}N{sub 4}.
NASA Astrophysics Data System (ADS)
Lee, Ming-Hsien; Liu, Po-Liang; Hong, Yung-An; Chou, Yen-Ting; Hong, Jia-Yang; Siao, Yu-Jin
2013-02-01
We conduct first-principles total-energy density functional calculations to study the band structures in Ge1-xSnx infrared semiconductor alloys. The norm-conserving optimized pseudopotentials of Ge and Sn have been constructed for electronic structure calculations. The composition-bandgap relationships in Ge1-xSnx lattices are evaluated by a detailed comparison of structural models and their electronic band structures. The critical Sn composition related to the transition from indirect- to direct-gap in Ge1-xSnx alloys is estimated to be as low as x ˜ 0.016 determined from the parametric fit. Our results show that the crossover Sn concentration occurs at a lower critical Sn concentration than the values predicted from the absorption measurements. However, early results indicate that the reliability of the critical Sn concentration from such measurements is hard to establish, since the indirect gap absorption is much weaker than the direct gap absorption. We find that the direct band gap decreases exponentially with the Sn composition over the range 0
NASA Astrophysics Data System (ADS)
Nagae, Yuki; Kurosawa, Masashi; Shibayama, Shigehisa; Araidai, Masaaki; Sakashita, Mitsuo; Nakatsuka, Osamu; Shiraishi, Kenji; Zaima, Shigeaki
2016-08-01
We have carried out density functional theory (DFT) calculation for Si1‑ x Sn x alloy and investigated the effect of the displacement of Si and Sn atoms with strain relaxation on the lattice constant and E– k dispersion. We calculated the formation probabilities for all atomic configurations of Si1‑ x Sn x according to the Boltzmann distribution. The average lattice constant and E– k dispersion were weighted by the formation probability of each configuration of Si1‑ x Sn x . We estimated the displacement of Si and Sn atoms from the initial tetrahedral site in the Si1‑ x Sn x unit cell considering structural relaxation under hydrostatic pressure, and we found that the breaking of the degenerated electronic levels of the valence band edge could be caused by the breaking of the tetrahedral symmetry. We also calculated the E– k dispersion of the Si1‑ x Sn x alloy by the DFT+U method and found that a Sn content above 50% would be required for the indirect–direct transition.
NASA Astrophysics Data System (ADS)
Akiyama, Toru; Nakamura, Kohji; Ito, Tomonori; Freeman, Arthur J.
2014-03-01
The electronic band structures and excited states of III-V semiconductors such as GaP, AlP, AlAs, and AlSb for various polytypes are determined employing the screened-exchange density functional calculations implemented in the full-potential linearized augmented plane-wave methods. We demonstrate that GaP and AlSb in the wurtzite (WZ) structure have direct gap while III-V semiconductors in the zinc blende, 4H, and 6H structures considered in this study exhibit an indirect gap. Furthermore, we find that inclusion of Al atoms less than 17% and 83% in the hexagonal AlxGa1-xP and AlxGa1-xAs alloys, respectively, leads to a direct transition with a gap energy of ˜2.3 eV. The feasibility of III-V semiconductors with a direct gap in WZ structure offers a possible crystal structure engineering to tune the optical properties of semiconductor materials.
Akiyama, Toru; Nakamura, Kohji; Ito, Tomonori; Freeman, Arthur J.
2014-03-31
The electronic band structures and excited states of III-V semiconductors such as GaP, AlP, AlAs, and AlSb for various polytypes are determined employing the screened-exchange density functional calculations implemented in the full-potential linearized augmented plane-wave methods. We demonstrate that GaP and AlSb in the wurtzite (WZ) structure have direct gap while III-V semiconductors in the zinc blende, 4H, and 6H structures considered in this study exhibit an indirect gap. Furthermore, we find that inclusion of Al atoms less than 17% and 83% in the hexagonal Al{sub x}Ga{sub 1−x}P and Al{sub x}Ga{sub 1−x}As alloys, respectively, leads to a direct transition with a gap energy of ∼2.3 eV. The feasibility of III-V semiconductors with a direct gap in WZ structure offers a possible crystal structure engineering to tune the optical properties of semiconductor materials.
Severino, Joyce Ferreira; Goodman, Bernard A; Kay, Christopher W M; Stolze, Klaus; Tunega, Daniel; Reichenauer, Thomas G; Pirker, Katharina F
2009-04-15
Electron paramagnetic resonance spectroscopy and density functional theory calculations have been used to investigate the redox properties of the green tea polyphenols (GTPs) (-)-epigallocatechin gallate (EGCG), (-)-epigallocatechin (EGC), and (-)-epicatechin gallate (ECG). Aqueous extracts of green tea and these individual phenols were autoxidized at alkaline pH and oxidized by superoxide anion (O(2)(-)) radicals in dimethyl sulfoxide. Several new aspects of the free radical chemistry of GTPs were revealed. EGCG can be oxidized on both the B and the D ring. The B ring was the main oxidation site during autoxidation, but the D ring was the preferred site for O(2)(-) oxidation. Oxidation of the D ring was followed by structural degradation, leading to generation of a radical identical to that of oxidized gallic acid. Alkaline autoxidation of green tea extracts produced four radicals that were related to products of the oxidation of EGCG, EGC, ECG, and gallic acid, whereas the spectra from O(2)(-) oxidation could be explained solely by radicals generated from EGCG. Assignments of hyperfine coupling constants were made by DFT calculations, allowing the identities of the radicals observed to be confirmed. PMID:19439236
Park, Ji-Sang; Kang, Joongoo; Yang, Ji-Hui; Metzger, Wyatt; Wei, Su-Huai
2015-01-15
Using first-principles density functional calculations, we investigate the relative stability and electronic structure of the grain boundaries (GBs) in zinc-blende CdTe. Among the low-Σ-value symmetric tilt Σ3 (111), Σ3 (112), Σ5 (120), and Σ5 (130) GBs, we show that the Σ3 (111)GB is always the most stable due to the absence of dangling bonds and wrong bonds. The Σ5 (120) GBs, however, are shown to be more stable than the Σ3 (112) GBs, even though the former has a higher Σ value, and the latter is often used as a model system to study GB effects in zinc-blende semiconductors. Furthermore,more » we find that although containing wrong bonds, the Σ5 (120) GBs are electrically benign due to the short wrong bond lengths, and thus are not as harmful as the Σ3 (112) GBs also having wrong bonds but with longer bond lengths.« less
Qian, Zekan; Li, Rui; Hou, Shimin; Xue, Zengquan; Sanvito, Stefano
2007-11-21
An efficient self-consistent approach combining the nonequilibrium Green's function formalism with density functional theory is developed to calculate electron transport properties of molecular devices with quasi-one-dimensional (1D) electrodes. Two problems associated with the low dimensionality of the 1D electrodes, i.e., the nonequilibrium state and the uncertain boundary conditions for the electrostatic potential, are circumvented by introducing the reflectionless boundary conditions at the electrode-contact interfaces and the zero electric field boundary conditions at the electrode-molecule interfaces. Three prototypical systems, respectively, an ideal ballistic conductor, a high resistance tunnel junction, and a molecular device, are investigated to illustrate the accuracy and efficiency of our approach. PMID:18035901
NASA Astrophysics Data System (ADS)
Wopperer, P.; Dinh, P. M.; Reinhard, P.-G.; Suraud, E.
2015-02-01
There are various ways to analyze the dynamical response of clusters and molecules to electromagnetic perturbations. Particularly rich information can be obtained from measuring the properties of electrons emitted in the course of the excitation dynamics. Such an analysis of electron signals covers observables such as total ionization, Photo-Electron Spectra (PES), Photoelectron Angular Distributions (PAD), and ideally combined PES/PAD. It has a long history in molecular physics and was increasingly used in cluster physics as well. Recent progress in the design of new light sources (high intensity, high frequency, ultra short pulses) opens new possibilities for measurements and thus has renewed the interest on these observables, especially for the analysis of various dynamical scenarios, well beyond a simple access to electronic density of states. This, in turn, has motivated many theoretical investigations of the dynamics of electronic emission for molecules and clusters up to such a complex and interesting system as C60. A theoretical tool of choice is here Time-Dependent Density Functional Theory (TDDFT) propagated in real time and on a spatial grid, and augmented by a Self-Interaction Correction (SIC). This provides a pertinent, robust, and efficient description of electronic emission including the detailed pattern of PES and PAD. A direct comparison between experiments and well founded elaborate microscopic theories is thus readily possible, at variance with more demanding observables such as for example fragmentation or dissociation cross sections. The purpose of this paper is to describe the theoretical tools developed on the basis of real-time and real-space TDDFT and to address in a realistic manner the analysis of electronic emission following irradiation of clusters and molecules by various laser pulses. After a general introduction, we shall present in a second part the available experimental results motivating such studies, starting from the simplest
NASA Astrophysics Data System (ADS)
Badaeva, Ekaterina; Feng, Yong; Gamelin, Daniel R.; Li, Xiaosong
2008-05-01
The electronic structures of pure and Co2+-doped ZnO quantum dots (QDs) with sizes up to 300 atoms were investigated using three different density functional theory approximations: local spin density approximation (LSDA), gradient-corrected Perdew-Burke-Ernzerhof (PBE) and the hybrid PBE1 functionals with LANL2DZ pseudo-potential and associated basis set. Qualitative agreement among the three methods is found for the pure ZnO nanostructures, but only the hybrid functional reproduces the correct bandgap energies quantitatively. For Co2+-doped ZnO QDs, both LSDA and PBE incorrectly model interactions between Co2+ d levels and the valence band of ZnO, which will strongly impair predictions of dopant-carrier magnetic exchange interactions based on such computational results. Experimental observations are reproduced well in calculations at the hybrid PBE1 level of theory, making this the method of choice for exploring the magnetism of transition metal ions in ZnO QDs computationally. The qualitative features of the Co2+ 3d levels do not change appreciably with changes in cluster size over the range examined, leading to size-dependent dopant-band edge energy differences. The results presented here thus provide an experimentally calibrated framework for future ab initio descriptions of dopant-carrier and dopant-dopant magnetic exchange interactions in diluted magnetic semiconductors (DMS) nanocrystals.
Wavelet analysis of electron-density maps.
Main, P; Wilson, J
2000-05-01
The wavelet transform is a powerful technique in signal processing and image analysis and it is shown here that wavelet analysis of low-resolution electron-density maps has the potential to increase their resolution. Like Fourier analysis, wavelet analysis expresses the image (electron density) in terms of a set of orthogonal functions. In the case of the Fourier transform, these functions are sines and cosines and each one contributes to the whole of the image. In contrast, the wavelet functions (simply called wavelets) can be quite localized and may only contribute to a small part of the image. This gives control over the amount of detail added to the map as the resolution increases. The mathematical details are outlined and an algorithm which achieves a resolution increase from 10 to 7 A using a knowledge of the wavelet-coefficient histograms, electron-density histogram and the observed structure amplitudes is described. These histograms are calculated from the electron density of known structures, but it seems likely that the histograms can be predicted, just as electron-density histograms are at high resolution. The results show that the wavelet coefficients contain the information necessary to increase the resolution of electron-density maps. PMID:10771431
Berland, Kristian; Arter, Calvin A; Cooper, Valentino R; Lee, Dr. Kyuho; Lundqvist, Prof. Bengt I.; Schroder, Prof. Elsebeth; Thonhauser, Prof. Timo; Hyldgaard, Per
2014-01-01
The theoretical description of sparse matter attracts much interest, in particular for those groundstate properties that can be described by density functional theory (DFT). One proposed approach, the van der Waals density functional (vdW-DF) method, rests on strong physical foundations and offers simple yet accurate and robust functionals. A very recent functional within this method called vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B, in print] stands out in its attempt to use an exchange energy derived from the same plasmon-based theory from which the nonlocal correlation energy was derived. Encouraged by its good performance for solids, layered materials, and aromatic molecules, we apply it to several systems that are characterized by competing interactions. These include the ferroelectric response in PbTiO3, the adsorption of small molecules within metal-organic frameworks (MOFs), the graphite/diamond phase transition, and the adsorption of an aromaticmolecule on the Ag(111) surface. Our results indicate that vdW-DF-cx is overall well suited to tackle these challenging systems. In addition to being a competitive density functional for sparse matter, the vdW-DF-cx construction presents a more robust general purpose functional that could be applied to a range of materials problems with a variety of competing interactions.
Berland, Kristian; Molecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, California 94720 ; Arter, Calvin A.; Thonhauser, T.; Cooper, Valentino R.; Lee, Kyuho; Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720 ; Lundqvist, Bengt I.; Schröder, Elsebeth; Hyldgaard, Per
2014-05-14
The theoretical description of sparse matter attracts much interest, in particular for those ground-state properties that can be described by density functional theory. One proposed approach, the van der Waals density functional (vdW-DF) method, rests on strong physical foundations and offers simple yet accurate and robust functionals. A very recent functional within this method called vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412 (2014)] stands out in its attempt to use an exchange energy derived from the same plasmon-based theory from which the nonlocal correlation energy was derived. Encouraged by its good performance for solids, layered materials, and aromatic molecules, we apply it to several systems that are characterized by competing interactions. These include the ferroelectric response in PbTiO{sub 3}, the adsorption of small molecules within metal-organic frameworks, the graphite/diamond phase transition, and the adsorption of an aromatic-molecule on the Ag(111) surface. Our results indicate that vdW-DF-cx is overall well suited to tackle these challenging systems. In addition to being a competitive density functional for sparse matter, the vdW-DF-cx construction presents a more robust general-purpose functional that could be applied to a range of materials problems with a variety of competing interactions.
Berland, Kristian; Arter, Calvin A; Cooper, Valentino R; Lee, Kyuho; Lundqvist, Bengt I; Schröder, Elsebeth; Thonhauser, T; Hyldgaard, Per
2014-05-14
The theoretical description of sparse matter attracts much interest, in particular for those ground-state properties that can be described by density functional theory. One proposed approach, the van der Waals density functional (vdW-DF) method, rests on strong physical foundations and offers simple yet accurate and robust functionals. A very recent functional within this method called vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412 (2014)] stands out in its attempt to use an exchange energy derived from the same plasmon-based theory from which the nonlocal correlation energy was derived. Encouraged by its good performance for solids, layered materials, and aromatic molecules, we apply it to several systems that are characterized by competing interactions. These include the ferroelectric response in PbTiO3, the adsorption of small molecules within metal-organic frameworks, the graphite/diamond phase transition, and the adsorption of an aromatic-molecule on the Ag(111) surface. Our results indicate that vdW-DF-cx is overall well suited to tackle these challenging systems. In addition to being a competitive density functional for sparse matter, the vdW-DF-cx construction presents a more robust general-purpose functional that could be applied to a range of materials problems with a variety of competing interactions. PMID:24832347
Otsuka, Takao; Sumita, Masato; Izawa, Hironori; Morihashi, Kenji
2016-07-21
Electron transfer (ET) states of 1-methyl-3-(N-(1,8-naphthalimidyl)ethyl)imidazolium iodide are responsible for its photophysics. Investigation of an ET state based on constrained density functional theory (CDFT) revealed that nonradiative decay from the ET excited state is mediated by the interaction of the iodine atom with the 1,8-naphthalimide or the imidazolium group. PMID:27222312
Bettinger, Holger F; Tönshoff, Christina; Doerr, Markus; Sanchez-Garcia, Elsa
2016-01-12
While the optical spectra of the acene series up to pentacene provide textbook examples for the annulation principle, the spectra of the larger members are much less understood. The present work provides an investigation of the optically allowed excited states of the acene series from pentacene to nonacene, the largest acene observed experimentally, using the density functional based multireference configuration method (DFT/MRCI). For this purpose, the ten lowest energy states of the B2u and B3u irreducible representations were computed. In agreement with previous computational investigations, the electronic wave functions of the acenes acquire significant multireference character with increasing acene size. The HOMO → LUMO excitation is the major contributor to the (1)La state (p band, B2u) also for the larger acenes. The oscillator strength decreases with increasing length. The (1)Lb state (α band, B3u), so far difficult to assign for the larger acenes due to overlap with photoprecursor bands, becomes almost insensitive to acene length. The (1)Bb state (β band, B3u) also moves only moderately to lower energy with increasing acene size. Excited states of B3u symmetry that formally result from double excitations involving HOMO, HOMO-1, LUMO, and LUMO+1 decrease in energy much faster with system size. One of them (D1) has very small oscillator strength but becomes almost isoenergetic with the (1)La state for nonacene. The other (D2) also has low oscillator strength as long as it is higher in energy than (1)Bb. Once it is lower in energy than the (1)Bb state, both states interact strongly resulting in two states with large oscillator strengths. The emergence of two strongly absorbing states is in agreement with experimental observations. The DFT/MRCI computations reproduce experimental excitation energies very well for pentacene and hexacene (within 0.1 eV). For the larger acenes deviations are larger (up to 0.2 eV), but qualitative agreement is observed. PMID
NASA Astrophysics Data System (ADS)
Kristyan, Sandor
2013-11-01
Using orbital-free framework, a simple numerical optimization of the density functional for ground state electronic energy is described for any type of functional approximation, demonstrated via the example of linear combinations of homogeneous functionals of the density. The numerical recipe is given and analyzed for solution: Originating from the linear dependence of nuclear-electron attraction functional on one-electron density ( V ne[ ρ 0( r 1)] = -ΣA = 1,…,MZA∫ ρ 0( r 1)rA1 -1d r 1), and a quadratic LCAO approximation for ρ 0, the optimization can be done with iterative use of lin-solver. This quadratic approximation, as simplest educated choice for ρ 0, is compared and analyzed algebraically to the HF-SCF one in the Appendices. We call the attention that the introduction of a self-consistent field optimization of non-linear density functional is a new element in this part of the related, general theory.
Vargas, Jorge; Springborg, Michael; Kirtman, Bernard
2014-02-07
The response to an electrostatic field is determined through simple model calculations, within both the restricted Hartree-Fock and density functional theory methods, for long, finite as well as infinite, periodic chains. The permanent dipole moment, μ{sub 0}, the polarizability, α, and the hyperpolarizabilities β and γ, calculated using a finite-field approach, are extensively analyzed. Our simple model allows for treatment of large systems and for separation of the properties into atomic and unit-cell contributions. That part of the response properties attributable to the terminations of the finite system change into delocalized current contributions in the corresponding infinite periodic system. Special emphasis is placed on analyzing the reasons behind the dramatic overestimation of the response properties found with density functional theory methods presently in common use.
Unexpected properties of a density functional
Karwowski, J.; Stanke, M.
2005-02-01
An observation on a pathological behavior of an exact density functional derived from either relativistic (Dirac) or nonrelativistic (Levy-Leblond) quantum-mechanical equation is reported. As expected, in the case of a one-electron atom the variational minimum of this functional is equal to the exact ground-state energy. However, apart from the correct density, this minimum is reached also by an infinite set of densities which do not correspond to the exact wave function. This paradoxical property of the functional is related to the multicomponent structure of both Dirac and Levy-Leblond wave functions. In particular, imposing the correct boundary conditions upon the trial densities removes only a part of the fake solutions. The results of this study demonstrate that in density-functional theories derived from models based on multicomponent wave functions, one should not expect any simple relation between the accuracy of the energy and the correctness of the corresponding density.
Mouesca, Jean-Marie
2014-01-01
The goal of this "how to" chapter is to present in a way as simple and practical as possible some of the concepts, key issues, and practices behind the so-called broken symmetry (BS) state which is widely used within the density functional theory (DFT) (for a very nice but thoughtful introduction to DFT (without equations!), read Perdew et al. (J Chem Theory Comput 5:902-908, 2009)) community to compute energetic as well as spectroscopic properties pertaining to (poly-)radicals, bioinorganic clusters (especially those containing transition metal ions), etc. Such properties encompass exchange coupling constants J (molecular magnetism) but also (among other things) g-tensors and hyperfine coupling tensors A (from electron paramagnetic resonance), isomer shifts δ and quadrupolar tensors ΔE Q (from Mössbauer), etc.Hopefully, this chapter will appeal to those DFT practitioners who would like to understand the basics behind the BS state and help them "demystify" some of the issues involved with them. More technical issues will only be alluded to, and appropriate references will be given for those interested to go beyond this mere introduction. This chapter is however not a review of the field. Consequently, it will be primarily based on my own experience. The goal here (in the spirit of a "how to" chapter) is to accompany the readers' thoughts in a progressive way along increasingly complex issues rather than encumbering the same thoughts with too complicate mathematical details (the few derivations which are given will therefore be explicit). Moreover, I will emphasize in this chapter the interplay between the computation of BS states on the one hand, and the derivation of phenomenological models on the other hand, whose parameters can be supplied from appropriate BS states. Finally, this chapter is dedicated to Louis Noodleman (Scripps Research Institute, CA, USA), pioneer (Noodleman, J Chem Phys 74:5737-5743, 1981; Noodleman, Chem Phys 109:131-143, 1986) and
Density-dependent covariant energy density functionals
Lalazissis, G. A.
2012-10-20
Relativistic nuclear energy density functionals are applied to the description of a variety of nuclear structure phenomena at and away fromstability line. Isoscalar monopole, isovector dipole and isoscalar quadrupole giant resonances are calculated using fully self-consistent relativistic quasiparticle randomphase approximation, based on the relativistic Hartree-Bogoliubovmodel. The impact of pairing correlations on the fission barriers in heavy and superheavy nuclei is examined. The role of pion in constructing desnity functionals is also investigated.
NASA Astrophysics Data System (ADS)
Zhang, Qian; Liang, Xiao; Chen, Bor-Yann; Chang, Chang-Tang
2015-12-01
This study tended to decipher the mechanism of photo degradation of azo dyes, which bond was favorable to be broken for application of wastewater decolorization. That is, from chemical structure perspective, the critical substituents to affect electron donor/acceptor for dye degradation would be identified in this research. The model reactive blacks (RB5), reactive blue 171 (RB171) and reactive red 198 (RR198) were degraded by graphene loaded TiO2, indicating how the electron withdrawing and releasing groups affect azo dye degradability. The byproducts and intermediate products were analyzed by ultraviolet-visible spectroscopy (UV-vis), gas chromatography-mass spectrometry (GC-MS) and ion chromatography (IC). Furthermore, the radicals involved in the reaction were found by electron paramagnetic resonance (ESR) to confirm the main oxidized species of hydroxyl radicals rather than the light generated positive holes. The finding revealed that the breakages of the bonds were due to the electron density changes around the bonds. This principle can be applicable not only for RB5 degradation, but also for reactive blue 171 (RB171), reactive red 198 (RR198) and some other textile dyes.
Density functionals from deep learning
NASA Astrophysics Data System (ADS)
McMahon, Jeffrey
Density-functional theory is a formally exact description of a many-body quantum system in terms of its density; in practice, however, approximations to the universal density functional (DF) are necessary. Machine learning has recently been proposed as a novel approach to discover such a DF (or components of it). Conventional machine learning algorithms, however, are limited in their ability to process data in their raw form, leading to invariance and/or sensitivity issues. In this presentation, an alternative approach based on deep learning will be demonstrated. Deep learning allows computational models that are capable of discovering intricate structure in large and/or high-dimensional data sets with multiple levels of abstraction, and do not suffer from the aforementioned issues. Results from the application of this approach to the prediction of the kinetic-energy DF of noninteracting electrons will be presented. Using theoretical results from computer science, a connection between the underlying model and the theorems of Hohenberg and Kohn will also be suggested.
Mazziotti, David A.
2010-06-15
An improved parametrization of the two-electron reduced density matrix (2-RDM) [D. A. Mazziotti, Phys. Rev. Lett. 101, 253002 (2008)] was recently shown to yield energies and properties that are markedly better than those calculated by traditional ab initio methods of similar computational scaling. In this paper a family of such energy functionals, generalizing the ones obtained previously, is derived through the use of (i) p-particle contraction relations based on the contraction of the cumulant expansions of p-particle RDMs and (ii) Cauchy-Schwarz relations that arise from an important set of N-representability constraints known as the two-positivity conditions. The 2-RDMs are explicitly parameterized in terms of the first-order part of the cumulant 2-RDM and, for the inclusion of single excitations, a second-order part of the 1-RDM. In contrast to earlier formulations based on the coefficients from configuration interaction with single and double excitations (CISD), the cumulant-based parametric 2-RDM methods, from the properties of cumulants, are rigorously size extensive. We also show that writing the energy functionals in terms of correlated 1-RDMs and cumulant 2-RDMs reduces the computational cost of the parametric 2-RDM methods to that of CISD. Applications are made to ground-state energies of several molecules, equilibrium bond distances, and frequencies of HF, F{sub 2}, and CO, the relative energy of the cis and trans isomers of HO{sub 3}{sup -}, and the HCN-HNC isomerization reaction. For bond breaking in hydrogen fluoride the improved and more efficient parametric 2-RDM methods yield energies with similar accuracies at both equilibrium and nonequilibrium geometries in 6-31G** and polarized valence quadruple-{zeta} basis sets. Computed 2-RDMs very nearly satisfy well-known N-representability conditions.
Arora, Priya; Moudgil, R. K.; Bhukal, Nisha
2015-05-15
Static density-density correlation function has been calculated for a spin-polarized two-dimensional quantum electron fluid by including the first-order exchange and self-energy corrections to the random-phase approximation (RPA). This is achieved by determining these corrections to the RPA linear density-density response function, obtained by solving the equation of motion for the single-particle Green’s function. Resulting infinite hierarchy of equations (involving higher-order Green’s functions) is truncated by factorizing the two-particle Green’s function as a product of the single-particle Green’s function and one-particle distribution function. Numerical results of correlation function are compared directly against the quantum Monte Carlo simulation data due to Tanatar and Ceperley for different coupling parameter (r{sub s}) values. We find almost exact agreement for r{sub s}=1, with a noticeable improvement over the RPA. Its quality, however, deteriorates with increasing r{sub s}, but correction to RPA is quite significant.
NASA Astrophysics Data System (ADS)
Pamuk, Betül; Baima, Jacopo; Dovesi, Roberto; Calandra, Matteo; Mauri, Francesco
2016-07-01
We investigate the capability of density functional theory (DFT) to appropriately describe the spin susceptibility, χs, and the intervalley electron-phonon coupling in LixZrNCl . At low doping, LixZrNCl behaves as a two-dimensional two-valley electron gas, with parabolic bands. In such a system, χs increases with decreasing doping because of the electron-electron interaction. We show that DFT with local functionals (LDA/GGA) is not capable of reproducing this behavior. The use of exact exchange in Hartree-Fock (HF) or in DFT hybrid functionals enhances χs. HF, B3LYP, and PBE0 approaches overestimate χs, whereas the range-separated HSE06 functional leads to results similar to those obtained in the random phase approximation (RPA) applied to a two-valley two-spin electron gas. Within HF, LixZrNCl is even unstable towards a ferromagnetic state for x <0.16 . The intervalley phonons induce an imbalance in the valley occupation that can be viewed as the effect of a pseudomagnetic field. Thus, similarly to what happens for χs, the electron-phonon coupling of intervalley phonons is enhanced by the electron-electron interaction. Only hybrid DFT functionals capture such an enhancement and the HSE06 functional reproduces the RPA results presented in M. Calandra et al. [Phys. Rev. Lett. 114, 077001 (2015), 10.1103/PhysRevLett.114.077001]. These results imply that the description of the susceptibility and electron-phonon coupling with a range-separated hybrid functional would be important also in other two-dimensional weakly doped semiconductors, such as transition-metal dichalcogenides and graphene.
NASA Astrophysics Data System (ADS)
Rüger, Robert; van Lenthe, Erik; Heine, Thomas; Visscher, Lucas
2016-05-01
We propose a new method of calculating electronically excited states that combines a density functional theory based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed time-dependent density functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB.
Rüger, Robert; van Lenthe, Erik; Heine, Thomas; Visscher, Lucas
2016-05-14
We propose a new method of calculating electronically excited states that combines a density functional theory based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed time-dependent density functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB. PMID:27179467
Arcelus, Oier; Li, Chunmei; Rojo, Teófilo; Carrasco, Javier
2015-06-01
Clarifying the electronic structure of sodium superoxide (NaO2) is a key step in understanding the electrochemical behavior of Na-O2 batteries. Here we report a density functional theory study to explore the effect of atomic structure and morphology on the electronic properties of different model systems: NaO2 bulk, (100) surface, and small (NaO2)n clusters (n = 3-8). We found that a correct description of the open-shell 2p electrons of O2(-) requires the use of a hybrid functional, which reveals a clear insulating nature of all of the investigated systems. This sheds light onto the capacity limitations of pure NaO2 as a discharge product and highlights the need for developing new strategies to enhance its electron transport in the optimization of Na-O2 cells. PMID:26266497
NASA Astrophysics Data System (ADS)
Bhandary, Sumanta; Schüler, Malte; Thunström, Patrik; di Marco, Igor; Brena, Barbara; Eriksson, Olle; Wehling, Tim; Sanyal, Biplab
2016-04-01
A proper theoretical description of the electronic structure of the 3 d orbitals in the metal centers of functional metalorganics is a challenging problem. We apply density functional theory and an exact diagonalization method in a many-body approach to study the ground-state electronic configuration of an iron porphyrin (FeP) molecule. Our study reveals that the consideration of multiple Slater determinants is important, and FeP is a potential candidate for realizing a spin crossover due to a subtle balance of crystal-field effects, on-site Coulomb repulsion, and hybridization between the Fe-d orbitals and ligand N-p states. The mechanism of switching between two close-lying electronic configurations of Fe-d orbitals is shown. We discuss the generality of the suggested approach and the possibility to properly describe the electronic structure and related low-energy physics of the whole class of correlated metal-centered organometallic molecules.
NASA Astrophysics Data System (ADS)
Perron, H.; Mellier, T.; Domain, C.; Roques, J.; Simoni, E.; Drot, R.; Catalette, H.
2007-08-01
Periodic density functional theory (DFT) calculations using plane-wave basis sets were performed in order to study the bulk of nickel ferrite NiFe2O4. The local spin density approximation (LSDA) and the generalized gradient approximation (GGA) formalism were used, and it appeared that the LSDA failed to describe the magnetic structure of this compound. However, the GGA formalism gave reliable results in good agreement with experimental data for the lattice parameters, the electronic properties and the bulk modulus. In addition, the calculated density of states of the metallic species d block as well as their local magnetic moments were correlated to the crystal-field theory. Then, a charge deformation map was computed and, as expected from the electronegativity scale, the electron excess is localized around oxygen atoms along the bond axes. The formation energies of metallic vacancies are in good agreement with the inverse spinel structure experimentally observed.