skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Direct and Simple Experimental Crystallographic Method to Calculate Partial Charges with Atomistic Correspondence to ab-initio Methods
Derivation of partial charges in small and large scale molecular systems is important for modeling of various experimental and theoretical properties like dipole moments, auto-correlation functions, charge disparity, understanding of dispersion, benchmark of classical MD simulations and electrostatic potential energy surface mapping. A correspondence between theoretical calculations (based on single/small number of molecules) is usually established with macroscopic IR/Raman spectra or dipole moment measurements. Such comparisons are indirect and lack a fine mapping of electrostatic potential from theory to experiment. In a new approach developed as the experimental part of this work, partial charges are calculated from crystallographic model refinement. The experimental method exhibits a satisfactory correspondence with partial charges obtained using quantum chemistry calculations. Particularly, gas phase partial charges from CHELPG method and condensed phase Lowdin charges correlate well and validate this experimental method.  more » « less
Award ID(s):
1800105
PAR ID:
10187701
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
American Crystallographic Association National Meeting 2020
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract The restrained electrostatic potential (RESP) approach is a highly regarded and widely used method of assigning partial charges to molecules for simulations. RESP uses a quantum-mechanical method that yields fortuitous overpolarization and thereby accounts only approximately for self-polarization of molecules in the condensed phase. Here we present RESP2, a next generation of this approach, where the polarity of the charges is tuned by a parameter, δ, which scales the contributions from gas- and aqueous-phase calculations. When the complete non-bonded force field model, including Lennard-Jones parameters, is optimized to liquid properties, improved accuracy is achieved, even with this reduced set of five Lennard-Jones types. We argue that RESP2 with δ  ≈ 0.6 (60% aqueous, 40% gas-phase charges) is an accurate and robust method of generating partial charges, and that a small set of Lennard-Jones types is a good starting point for a systematic re-optimization of this important non-bonded term. 
    more » « less
  2. Of interest in understanding electronic structure, bulk physical properties, enthalpies of phase changes, dipole moments, and numerous other properties of molecules, is the determination of realistic partial atomic charges on atoms. Atoms may take on partial positive or negative charges due to polar covalent bonds, coordinate covalent bonds, or due to formal charges imposed by Lewis structure constraints. Traditional crystallographic refinement treats each atom as a neutral, spherical atom however. We present a ongoing developments of a mode of crystallographic model refinement that permits refinement of electron density at individual atoms in order to arrive at partial atomic charges of atoms in a crystallographic model. Comparison to calculated partial charges (CHELPG, NBO, Mulliken) from quantum calculations (DFT, MP2) in both the gas phase and crystalline state will be presented. 
    more » « less
  3. Modern crystallographic refinement methods treat each atom in a molecule as neutral with spherical electron density. Atoms, however, exhibit partial atomic charges arising from intramolecular forces via bonding. These partial charges are crucial for understanding electronic structure and bulk physical properties of molecules. Typically the polarity and polarizability of molecules are calculated using IR and Raman spectroscopy, respectively. While these techniques can be used on small molecules, fine elucidation of partial charges on individual atoms is still unrealized. Here we present crystallographic refinement developments that allow us to refine electron density around individual atoms to experimentally calculate partial atomic charges. Comparison between these experimentally calculated charges to theoretical quantum calculated charges will also be presented. 
    more » « less
  4. Abstract We report the application of our fragment‐based quantum chemistry model MIM (Molecules‐In‐Molecules) with electrostatic embedding. The method is termed “EE‐MIM (ElectrostaticallyEmbeddedMolecules‐In‐Molecules)” and accounts for the missing electrostatic interactions in the subsystems resulting from fragmentation. Point charges placed at the atomic positions are used to represent the interaction of each subsystem with the rest of the molecule with minimal increase in the computational cost. We have carefully calibrated this model on a range of different sizes of clusters containing up to 57 water molecules. The fragmentation methods have been applied with the goal of reproducing the unfragmented total energy at the MP2/6‐311G(d,p) level. Comparative analysis has been carried out between MIM and EE‐MIM to gauge the impact of electrostatic embedding. Performance of several different parameters such as the type of charge and levels of fragmentation are analyzed for the prediction of absolute energies. The use of background charges in subsystem calculations improves the performance of both one‐ and two‐layer MIM while it is noticeably important in the case of one‐layer MIM. Embedded charges for two‐layer MIM are obtained from a full system calculation at the low‐level. For one‐layer MIM, in the absence of a full system calculation, two different types of embedded charges, namely, Geometry dependent (GD) and geometry independent (GI) charges, are used. A self‐consistent procedure is employed to obtain GD charges. We have further tested our method on challenging charged systems with stronger intermolecular interactions, namely, protonated ammonia clusters (containing up to 30 ammonia molecules). The observations are similar to water clusters with improved performance using embedded charges. Overall, the performance of NPA charges as embedded charges is found to be the best. 
    more » « less
  5. Two-dimensional layered transition metal dichalcogenides are potential thermoelectric candidates with application in on-chip integrated nanoscale cooling and power generation. Here, we report a comprehensive experimental and theoretical study on the in-plane thermoelectric transport properties of thin 2H-MoTe2 flakes prepared in field-effect transistor geometry to enable electrostatic gating and modulation of the electronic properties. The thermoelectric power factor is enhanced by up to 45% using electrostatic modulation. The in-plane thermal conductivity of 9.8 ± 3.7 W m−1 K−1 is measured using the heat diffusion imaging method in a 25 nm thick flake. First-principles calculations are used to obtain the electronic band structure, phonon band dispersion, and electron–phonon scattering rates. The experimental electronic properties are in agreement with theoretical results obtained within energy-dependent relaxation time approximation. The thermal conductivity is evaluated using both the relaxation time approximation and the full iterative solution to the phonon Boltzmann transport equation. This study establishes a framework to quantitively compare first-principle-based calculations with experiments in 2D layered materials. 
    more » « less