Electrostatic Potential Fitting Method Using Constrained Spatial Electron Density Expanded with Preorthogonal Natural Atomic Orbitals

2020 ◽  
Vol 124 (46) ◽  
pp. 9665-9673
Author(s):  
Daisuke Yokogawa ◽  
Kayo Suda

Computation ◽  
2019 ◽  
Vol 7 (4) ◽  
pp. 64
Author(s):  
Rafael Lopez ◽  
Frank Martinez ◽  
Ignacio Ema ◽  
Jose Manuel Garcia de la Vega ◽  
Guillermo Ramirez

An algorithm for the efficient computation of molecular electrostatic potential is reported. It is based on the partition/expansion of density into (pseudo) atomic fragments with the method of Deformed Atoms in Molecules, which allows to compute the potential as a sum of atomic contributions. These contributions are expressed as a series of irregular spherical harmonics times effective multipole moments and inverse multipole moments, including short-range terms. The problem is split into two steps. The first one consists of the partition/expansion of density accompanied by the computation of multipole moments, and its cost depends on the size of the basis set used in the computation of electron density within the Linear Combination of Atomic Orbitals framework. The second one is the actual computation of the electrostatic potential from the quantities calculated in the first step, and its cost depends on the number of computation points. For a precision in the electrostatic potential of six decimal figures, the algorithm leads to a dramatic reduction of the computation time with respect to the calculation from electron density matrix and integrals involving basis set functions.



Molecules ◽  
2021 ◽  
Vol 26 (14) ◽  
pp. 4227
Author(s):  
Alessandro Cossard ◽  
Silvia Casassa ◽  
Carlo Gatti ◽  
Jacques K. Desmarais ◽  
Alessandro Erba

The chemistry of f-electrons in lanthanide and actinide materials is yet to be fully rationalized. Quantum-mechanical simulations can provide useful complementary insight to that obtained from experiments. The quantum theory of atoms in molecules and crystals (QTAIMAC), through thorough topological analysis of the electron density (often complemented by that of its Laplacian) constitutes a general and robust theoretical framework to analyze chemical bonding features from a computed wave function. Here, we present the extension of the Topond module (previously limited to work in terms of s-, p- and d-type basis functions only) of the Crystal program to f- and g-type basis functions within the linear combination of atomic orbitals (LCAO) approach. This allows for an effective QTAIMAC analysis of chemical bonding of lanthanide and actinide materials. The new implemented algorithms are applied to the analysis of the spatial distribution of the electron density and its Laplacian of the cesium uranyl chloride, Cs2UO2Cl4, crystal. Discrepancies between the present theoretical description of chemical bonding and that obtained from a previously reconstructed electron density by experimental X-ray diffraction are illustrated and discussed.



Molecules ◽  
2020 ◽  
Vol 25 (19) ◽  
pp. 4409
Author(s):  
Carlo Gatti ◽  
Alessandro Dessì ◽  
Roberto Dallocchio ◽  
Victor Mamane ◽  
Sergio Cossu ◽  
...  

Positive electrostatic potential (V) values are often associated with σ- and π-holes, regions of lower electron density which can interact with electron-rich sites to form noncovalent interactions. Factors impacting σ- and π-holes may thus be monitored in terms of the shape and values of the resulting V. Further precious insights into such factors are obtained through a rigorous decomposition of the V values in atomic or atomic group contributions, a task here achieved by extending the Bader–Gatti source function (SF) for the electron density to V. In this article, this general methodology is applied to a series of 4,4′-bipyridine derivatives containing atoms from Groups VI (S, Se) and VII (Cl, Br), and the pentafluorophenyl group acting as a π-hole. As these molecules are characterized by a certain degree of conformational freedom due to the possibility of rotation around the two C–Ch bonds, from two to four conformational motifs could be identified for each structure through conformational search. On this basis, the impact of chemical and conformational features on σ- and π-hole regions could be systematically evaluated by computing the V values on electron density isosurfaces (VS) and by comparing and dissecting in atomic/atomic group contributions the VS maxima (VS,max) values calculated for different molecular patterns. The results of this study confirm that both chemical and conformational features may seriously impact σ- and π-hole regions and provide a clear analysis and a rationale of why and how this influence is realized. Hence, the proposed methodology might offer precious clues for designing changes in the σ- and π-hole regions, aimed at affecting their potential involvement in noncovalent interactions in a desired way.



Author(s):  
Chinnasamy Kalaiarasi ◽  
Mysore S Pavan ◽  
Poomani Kumaradhas

An experimental charge density distribution of 2-nitroimidazole was determined from high-resolution X-ray diffraction and the Hansen–Coppens multipole model. The 2-nitroimidazole compound was crystallized and a high-angle X-ray diffraction intensity data set has been collected at low temperature (110 K). The structure was solved and further, an aspherical multipole model refinement was performed up to octapole level; the results were used to determine the structure, bond topological and electrostatic properties of the molecule. In the crystal, the molecule exhibits a planar structure and forms weak and strong intermolecular hydrogen-bonding interactions with the neighbouring molecules. The Hirshfeld surface of the molecule was plotted, which explores different types of intermolecular interactions and their strength. The topological analysis of electron density at the bond critical points (b.c.p.) of the molecule was performed, from that the electron density ρbcp(r) and the Laplacian of electron density ∇2ρbcp(r) at the b.c.p.s of the molecule have been determined; these parameters show the charge concentration/depletion of the nitroimidazole bonds in the crystal. The electrostatic parameters like atomic charges and the dipole moment of the molecule were calculated. The electrostatic potential surface of the molecule has been plotted, and it displays a large electronegative region around the nitro group. All the experimental results were compared with the corresponding theoretical calculations performed usingCRYSTAL09.



2000 ◽  
Vol 56 (1) ◽  
pp. 155-165 ◽  
Author(s):  
Farid Benabicha ◽  
Virginie Pichon-Pesme ◽  
Christian Jelsch ◽  
Claude Lecomte ◽  
Ahmed Khmou

The experimental electron density distribution in glycyl-L-threonine dihydrate has been investigated using single-crystal X-ray diffraction data at 110 K to a resolution of (sin θ/λ) = 1.2 Å−1. Multipolar pseudo-atom refinement was carried out against 5417 observed data and the molecular electron density was analyzed using topological methods. The experimental electrostatic potential around the molecule is discussed in terms of molecular interactions. Crystal data: C6H12N2O4.2H2O, Mr = 212.2, orthorhombic, P212121, Z = 4, F(000) = 456 e, T = 110 K, a = 9.572 (3), b = 10.039 (3), c = 10.548 (2) Å, V = 1013.6 (4) Å3, Dx = 1.3 g cm−3, µ = 1.2 cm−1 for λMo = 0.7107 Å.



Sign in / Sign up

Export Citation Format

Share Document