Calculating and Characterizing the Charge Distributions in Solids
Indrani Choudhuri1, Donald G Truhlar1
1Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States.
Abstract:
Accurate estimation of the partial atomic charges on metal centers is useful for understanding electronic and catalytic properties of materials. However, different methods of calculating these charges may give quite different results; this issue has been more widely studied for molecules than for solids. Here we study the charges on the metal centers of a test set of 18 solids containing transition metals by using density functional theory with several density functionals (PBE, PBE+U, TPSS, revTPSS, HLE17, revM06-L, B3LYP, B3LYP*, and other exchange-modified B3LYP functionals) and four charge models (Bader, Hirshfeld, CM5, and DDEC6). The test set contains 12 systems with nonmagnetic metal centers (eight metal oxides (MO2), two metal sulfides (MS2), and two metal selenides (MSe2)) and six ferromagnetic transition metal complexes. Our study shows that, among the four types of charges, Bader charges are the highest and Hirshfeld charges are the lowest for all the systems, regardless of the functional being used. The CM5 charges are bigger than DDEC6 charges for MX2 with M = Ti or Mo and X = S or Se, but for the other 14 cases they are lower. We found that the most of the systems are sensitive to the Hubbard U parameters in PBE+U and to the percentage X of Hartree-Fock exchange in exchange-modified B3LYP; as we increase U or X, the charges on the metal atoms in MX2 increase steadily. Testing different density functionals shows charges calculated with higher Hubbard U parameters in PBE+U are comparable to B3LYP (with 20% Hartree-Fock exchange). Among four meta-GGA functionals studied, the charges with HLE17 have the closest agreement with B3LYP. The variation of charges with choice of charge model is greater than the variation with choice of density functional.
Related Concept Videos
Continuous Charge Distributions
The electric charge can also be subjected to an analogical...
Electric Field of a Charged Disk
The system's symmetry is in the cylindrical directions across the plane of the charge. As a result, the electric fields created by various surface charge elements nullify each other in the direction parallel to the surface. Thereby, the resulting electric field is perpendicular to the plane. Since the disk is...
Charge on a Conductor
Energy Associated With a Charge Distribution
Calculations of Electric Potential I
The ring is divided into infinitesimal small arcs such that point M is equidistant from all the arcs. Here, the cylindrical coordinate system is used to calculate the electric potential at point M. A general element of the arc between angles θ and θ + dθ is of the length Rdθ and has a charge of...
Potential Due to a Polarized Object


