Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEnergies of ions in water and nanopores within Density Functional Theory
105
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Accurate calculations of electrostatic potentials and treatment of substrate polarizability are critical for predicting the permeation of ions inside water-filled nanopores. The {it ab initio} molecular dynamics method (AIMD), based on Density Functional Theory (DFT), accounts for the polarizability of materials, water, and solutes, and it should be the method of choice for predicting accurate electrostatic energies of ions. In practice, DFT coupled with the use of periodic boundary conditions in a charged system leads to large energy shifts. Results obtained using different DFT packages may vary because of the way pseudopotentials and long-range electrostatics are implemented. Using maximally localized Wannier functions, we apply robust corrections that yield relatively unambiguous ion energies in select molecular and aqueous systems and inside carbon nanotubes. Large binding energies are predicted for ions in metallic carbon nanotube arrays, while with consistent definitions Na$^+$ and Cl$^-$ energies are found to exhibit asymmetries comparable with those computed using non-polarizable water force fields.
rate research
Read More
Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is know about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single particle calculations tend to overestimate spin gaps while they underestimate charge gaps.
Localized basis sets in the projector augmented wave formalism allow for computationally efficient calculations within density functional theory (DFT). However, achieving high numerical accuracy requires an extensive basis set, which also poses a fundamental problem for the interpretation of the results. We present a way to obtain a reduced basis set of atomic orbitals through the subdiagonalization of each atomic block of the Hamiltonian. The resulting local orbitals (LOs) inherit the information of the local crystal field. In the LO basis, it becomes apparent that the Hamiltonian is nearly block-diagonal, and we demonstrate that it is possible to keep only a subset of relevant LOs which provide an accurate description of the physics around the Fermi level. This reduces to some extent the redundancy of the original basis set, and at the same time it allows one to perform post-processing of DFT calculations, ranging from the interpretation of electron transport to extracting effective tight-binding Hamiltonians, very efficiently and without sacrificing the accuracy of the results.
Excitons are electron-hole pairs appearing below the band gap in insulators and semiconductors. They are vital to photovoltaics, but are hard to obtain with time-dependent density-functional theory (TDDFT), since most standard exchange-correlation (xc) functionals lack the proper long-range behavior. Furthermore, optical spectra of bulk solids calculated with TDDFT often lack the required resolution to distinguish discrete, weakly bound excitons from the continuum. We adapt the Casida equation formalism for molecular excitations to periodic solids, which allows us to obtain exciton binding energies directly. We calculate exciton binding energies for both small- and large-gap semiconductors and insulators, study the recently proposed bootstrap xc kernel [S. Sharma et al., Phys. Rev. Lett. 107, 186401 (2011)], and extend the formalism to triplet excitons.
The accurate description of the optical spectra of insulators and semiconductors remains an important challenge for time-dependent density-functional theory (TDDFT). Evidence has been given in the literature that TDDFT can produce bound as well as continuum excitons for specific systems, but there are still many unresolved basic questions concerning the role of dynamical exchange and correlation (xc). In particular, the role of the long spatial range and the frequency dependence of the xc kernel $f_{rm xc}$ for excitonic binding are still not very well explored. We present a minimal model for excitons in TDDFT, consisting of two bands from a one-dimensional Kronig-Penney model and simple approximate xc kernels, which allows us to address these questions in a transparent manner. Depending on the system, it is found that adiabatic xc kernels can produce a single bound exciton, and sometimes two bound excitons, where the long spatial range of $f_{rm xc}$ is not a necessary condition. It is shown how the Wannier model, featuring an effective electron-hole interaction, emerges from TDDFT. The collective, many-body nature of excitons is explicitly demonstrated.
Due to the strongly nonlocal nature of $f_{xc}({bf r},{bf r},omega)$ the {em scalar} exchange and correlation (xc) kernel of the time-dependent density-functional theory (TDDFT), the formula for Q the friction coefficient of an interacting electron gas (EG) for ions tends to give a too large value of Q for heavy ions in the medium- and low-density EG, if we adopt the local-density approximation (LDA) to $f_{xc}({bf r},{bf r},omega)$, even though the formula itself is formally exact. We have rectified this unfavorable feature by reformulating the formula for Q in terms of the {em tensorial} xc kernel of the time dependent current-density functional theory, to which the LDA can be applied without intrinsic difficulty. Our numerical results find themselves in a considerably better agreement with the experimental stopping power of Al and Au for slow ions than those previously obtained within the LDA to the TDDFT.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
