1
|
Mairesse F, Maschio L, Champagne B. First-principles evaluation of the second harmonic generation response of reference organic and inorganic crystals. J Chem Phys 2023; 158:064707. [PMID: 36792490 DOI: 10.1063/5.0137274] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Using the CRYSTAL17 package at the coupled-perturbed Kohn-Sham (CPKS) level, periodic boundary conditions first-principles calculations are enacted to predict the second harmonic generation second-order nonlinear optical (NLO) susceptibility, χ(2), values of six historical NLO crystals. This selection allowed the comparison between state-of-the-art calculations and experiment. Several computational aspects are tackled to define conditions where the results are converged with respect to the range of lattice summations, to the number of k-points in the first Brillouin zone, to the order of the multipole expansions for evaluating the long-range part of the electrostatic interactions, as well as to the atomic basis set size. A valence triple zeta basis set supplemented with polarization functions has been selected. Then, χ(2) calculations have been performed using a range of exchange-correlation functionals (XCFs). Results show the large impact of the amount of Hartree-Fock (HF) exchange on the amplitude but also on the sign on the χ(2) tensor components. To a given extent, these amplitude effects are consistent with results on molecules, but the sign reversal effects and the non-monotonic behavior of the χ(2) tensor components as a function of the amount of HF exchange are scarcely found for molecules. Then, using the recommended range-separated hybrid XCFs, the CPKS scheme leads to good agreement with experimental data for potassium dihydrogenophosphate, urea, and χZXX (2) of LiNbO3. The agreement is more questionable for χZZZ (2) of LiNbO3 whereas it remains poor for ammonium dihydrogenophosphate and 2-methyl-4-nitroaniline, with large underestimations by about a factor of 3, opening a path to further fine-tuning of the ranges of inclusion of HF exchange.
Collapse
Affiliation(s)
- François Mairesse
- Theoretical Chemistry Laboratory, Unit of Theoretical and Structural Physical Chemistry, NISM (Namur Institute of Structured Matter), University of Namur (UNamur), B-5000 Namur, Belgium
| | - Lorenzo Maschio
- Dipartimento di Chimica, Università degli Studi di Torino (UNITO), Torino, Italy
| | - Benoît Champagne
- Theoretical Chemistry Laboratory, Unit of Theoretical and Structural Physical Chemistry, NISM (Namur Institute of Structured Matter), University of Namur (UNamur), B-5000 Namur, Belgium
| |
Collapse
|
2
|
Weatherly J, Macchi P, Volkov A. On the calculation of the electrostatic potential, electric field and electric field gradient from the aspherical pseudoatom model. II. Evaluation of the properties in an infinite crystal. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2021; 77:399-419. [PMID: 34473095 DOI: 10.1107/s2053273321005532] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Accepted: 05/28/2021] [Indexed: 11/10/2022]
Abstract
The previously reported exact potential and multipole moment (EP/MM) method for fast and precise evaluation of the intermolecular electrostatic interaction energies in molecular crystals using the pseudoatom representation of the electron density [Nguyen, Macchi & Volkov (2020), Acta Cryst. A76, 630-651] has been extended to the calculation of the electrostatic potential (ESP), electric field (EF) and electric field gradient (EFG) in an infinite crystal. The presented approach combines an efficient Ewald-type summation (ES) of atomic multipoles up to the hexadecapolar level in direct and reciprocal spaces with corrections for (i) the net polarization of the sample (the `surface term') due to a net dipole moment of the crystallographic unit cell (if present) and (ii) the short-range electron-density penetration effects. The rederived and reported closed-form expressions for all terms in the ES algorithm have been augmented by the expressions for the surface term available in the literature [Stenhammar, Trulsson & Linse (2011), J. Chem. Phys. 134, 224104] and the exact potential expressions reported in a previous study [Volkov, King, Coppens & Farrugia (2006), Acta Cryst. A62, 400-408]. The resulting algorithm, coded using Fortran in the XDPROP module of the software package XD, was tested on several small molecular crystal systems (formamide, benzene, L-dopa, paracetamol, amino acids etc.) and compared with a series of EP/MM-based direct-space summations (DS) performed within a certain number of unit cells generated along both the positive and negative crystallographic directions. The EP/MM-based ES technique allows for a noticeably more precise determination of the EF and EFG and significantly better precision of the evaluated ESP when compared with the DS calculations, even when the latter include contributions from an array of symmetry-equivalent atoms generated within four additional unit cells along each crystallographic direction. In terms of computational performance, the ES/EP/MM method is significantly faster than the DS calculations performed within the extended unit-cell limits but trails the DS calculations within the reduced summation ranges. Nonetheless, the described EP/MM-based ES algorithm is superior to the direct-space summations as it does not require the user to monitor continuously the convergence of the evaluated properties as a function of the summation limits and offers a better precision-performance balance.
Collapse
Affiliation(s)
- Jessie Weatherly
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| | - Piero Macchi
- Department of Chemistry, Materials and Chemical Engineering, Politecnico di Milano, Via Mancinelli 7, Milano, 20131, Italy
| | - Anatoliy Volkov
- Department of Chemistry and Computational Science Program, Middle Tennessee State University, Murfreesboro, TN 37132, USA
| |
Collapse
|
3
|
Konovalova IS, Muzyka EN, Urzhuntseva VV, Shishkina SV. Role of intermolecular interactions in formation of mono- and diaminopyridine crystals: study from the energetic viewpoint. Struct Chem 2021. [DOI: 10.1007/s11224-020-01625-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
4
|
Medviediev V, Shishkina S, Ribalka AO, Zaręba JK, Drozd M, Daszkiewicz M. Revisiting 2-chloro-4-nitroaniline: analysis of intricate supramolecular ordering of a triclinic polymorph featuring a high Z value and strong second harmonic generation. CrystEngComm 2020. [DOI: 10.1039/d0ce00582g] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
A new polymorph of 2-chloro-4-nitroaniline is more stable than its 55 year-old antecedent.
Collapse
Affiliation(s)
- Volodymyr Medviediev
- Institute of Low Temperature and Structure Research
- Polish Academy of Sciences
- 50-422 Wrocław
- Poland
| | - Svitlana Shishkina
- SSI “Institute for Single Crystals”
- National Academy of Science of Ukraine
- Kharkiv 61001
- Ukraine
- V. N. Karazin Kharkiv National University
| | - A. O. Ribalka
- V. N. Karazin Kharkiv National University
- Kharkiv 61077
- Ukraine
| | - Jan K. Zaręba
- Advanced Materials Engineering and Modelling Group
- Faculty of Chemistry
- Wroclaw University of Science and Technology
- 50-370 Wrocław
- Poland
| | - Marek Drozd
- Institute of Low Temperature and Structure Research
- Polish Academy of Sciences
- 50-422 Wrocław
- Poland
| | - Marek Daszkiewicz
- Institute of Low Temperature and Structure Research
- Polish Academy of Sciences
- 50-422 Wrocław
- Poland
| |
Collapse
|
5
|
Shi MW, Thomas SP, Hathwar VR, Edwards AJ, Piltz RO, Jayatilaka D, Koutsantonis GA, Overgaard J, Nishibori E, Iversen BB, Spackman MA. Measurement of Electric Fields Experienced by Urea Guest Molecules in the 18-Crown-6/Urea (1:5) Host-Guest Complex: An Experimental Reference Point for Electric-Field-Assisted Catalysis. J Am Chem Soc 2019; 141:3965-3976. [PMID: 30761898 DOI: 10.1021/jacs.8b12927] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
High-resolution synchrotron and neutron single-crystal diffraction data of 18-crown-6/(pentakis)urea measured at 30 K are combined, with the aim of better appreciating the electrostatics associated with intermolecular interactions in condensed matter. With two 18-crown-6 molecules and five different urea molecules in the crystal, this represents the most ambitious combined X-ray/synchrotron and neutron experimental charge density analysis to date on a cocrystal or host-guest system incorporating such a large number of unique molecules. The dipole moments of the five urea guest molecules in the crystal are enhanced considerably compared to values determined for isolated molecules, and 2D maps of the electrostatic potential and electric field show clearly how the urea molecules are oriented with dipole moments aligned along the electric field exerted by their molecular neighbors. Experimental electric fields in the range of 10-19 GV m-1, obtained for the five different urea environments, corroborate independent measurements of electric fields in the active sites of enzymes and provide an important experimental reference point for recent discussions focused on electric-field-assisted catalysis.
Collapse
Affiliation(s)
- Ming W Shi
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| | - Sajesh P Thomas
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia.,Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark
| | - Venkatesha R Hathwar
- Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark.,Division of Physics, Faculty of Pure and Applied Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8571 , Japan
| | - Alison J Edwards
- Australian Nuclear Science and Technology Organization , Australian Centre for Neutron Scattering , New Illawarra Road , Lucas Heights , New South Wales 2234 , Australia
| | - Ross O Piltz
- Australian Nuclear Science and Technology Organization , Australian Centre for Neutron Scattering , New Illawarra Road , Lucas Heights , New South Wales 2234 , Australia
| | - Dylan Jayatilaka
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| | - George A Koutsantonis
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| | - Jacob Overgaard
- Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark
| | - Eiji Nishibori
- Division of Physics, Faculty of Pure and Applied Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8571 , Japan
| | - Bo B Iversen
- Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark
| | - Mark A Spackman
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| |
Collapse
|
6
|
Raikwar MM, Rhyman L, Ramasami P, Sekar N. Theoretical Investigation of Difluoroboron Complex of Curcuminoid Derivatives with and without Phenyl Substituent (at Meso Position): Linear and Non-Linear Optical Study. ChemistrySelect 2018. [DOI: 10.1002/slct.201802231] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Manish M. Raikwar
- Dyestuff Technology Department; Institute of Chemical Technology Nathalal Parekh Marg; Matunga, Mumbai 400019 India
| | - Lydia Rhyman
- Department of Chemistry; Computational Chemistry Group; Faculty of Science University of Mauritius; Réduit 80837 Mauritius
- Department of Applied Chemistry; University of Johannesburg; Doornfontein Campus; Johannesburg 2028 South Africa
| | - Ponnadurai Ramasami
- Department of Chemistry; Computational Chemistry Group; Faculty of Science University of Mauritius; Réduit 80837 Mauritius
- Department of Applied Chemistry; University of Johannesburg; Doornfontein Campus; Johannesburg 2028 South Africa
| | - Nagaiyan Sekar
- Dyestuff Technology Department; Institute of Chemical Technology Nathalal Parekh Marg; Matunga, Mumbai 400019 India
| |
Collapse
|
7
|
Hübschle CB, Ruhmlieb C, Burkhardt A, van Smaalen S, Dittrich B. On avoiding negative electron density in Gram-Charlier refinements of anharmonic motion: the example of glutathione. Z KRIST-CRYST MATER 2018. [DOI: 10.1515/zkri-2018-2060] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Abstract
The structure of glutathione, γ-l-Glutamyl-l-cysteinyl-glycine (C10H17N3O6S), was studied by multi-temperature single-crystal X-ray diffraction. Residual density maps from conventional independent atom model refinement gave indication of anharmonic motion in the molecule. This was further investigated by invariom refinement with anisotropic displacement parameters for all atoms, which described asphericity due to chemical bonding and lone pairs; afterwards only the residual-density signal of anharmonic motion remained. Treating anharmonicity with third-order Gram-Charlier displacement parameters led to regions with unphysical negative electron density. In contrast, a maximum entropy method (MEM) determination of the electron density successfully takes the features into account. Respective difference electron density plots (MEM minus prior and [Invariom+GC] minus invariom) agree well with each other. Challenges in treating and understanding the phenomenon are discussed. A procedure is proposed how unphysical negative electron density can be avoided. It is closely related to the free lunch algorithm.
Collapse
Affiliation(s)
- Christian B. Hübschle
- Laboratory of Crystallography , University of Bayreuth, BGI-Gebäude , 95440 Bayreuth , Germany
| | - Charlotte Ruhmlieb
- Institut für Physikalische Chemie der Universität Hamburg , Martin-Luther-King-Platz 6 , 20146 Hamburg , Germany
| | - Anja Burkhardt
- Deutsches Elektronensynchrotron DESY, Hamburg , Notkestrasse 85 , 22607 Hamburg , Germany
| | - Sander van Smaalen
- Laboratory of Crystallography , University of Bayreuth, BGI-Gebäude , 95440 Bayreuth , Germany
| | - Birger Dittrich
- Heinrich-Heine Universität Düsseldorf , Institut für Anorganische Chemie und Strukturchemie , Universitätsstraße 1 , 40225 Düsseldorf , Germany
| |
Collapse
|
8
|
Demircioğlu Z, Özdemir FA, Dayan O, Şerbetçi Z, Özdemir N. Synthesis, X-ray diffraction method, spectroscopic characterization (FT-IR, 1H and 13C NMR), antimicrobial activity, Hirshfeld surface analysis and DFT computations of novel sulfonamide derivatives. J Mol Struct 2018. [DOI: 10.1016/j.molstruc.2018.02.063] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
9
|
Lu SI. Application of discrete solvent reaction field model with self-consistent atomic charges and atomic polarizabilities to calculate the χ(1) and χ(2) of organic molecular crystals. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2017.10.036] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
10
|
Woińska M, Jayatilaka D, Dittrich B, Flaig R, Luger P, Woźniak K, Dominiak PM, Grabowsky S. Validation of X-ray Wavefunction Refinement. Chemphyschem 2017; 18:3334-3351. [PMID: 29168318 DOI: 10.1002/cphc.201700810] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2017] [Revised: 09/22/2017] [Indexed: 11/10/2022]
Abstract
In this work, the quality of the electron density in crystals reconstructed by the multipolar model (MM) and by X-ray wavefunction refinement (XWR) is tested on a set of high-resolution X-ray diffraction data sets of four amino acids and six tripeptides. It results in the first thorough validation of XWR. Agreement statistics, figures of merit, residual- and deformation-density maps, as well as atomic displacement parameters are used to measure the quality of the reconstruction relative to the measured structure factors. Topological analysis of the reconstructed density is carried out to obtain atomic and bond-topological properties, which are subsequently compared to the values derived from benchmarking periodic DFT geometry optimizations. XWR is simultaneously in better agreement than the MM with both benchmarking theory and the measured diffraction pattern. In particular, the obvious problems with the description of polar bonds in the MM are significantly reduced by using XWR. Similarly, modeling of electron density in the vicinity of hydrogen atoms with XWR is visibly improved.
Collapse
Affiliation(s)
- Magdalena Woińska
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Zwirki i Wigury 101, 02-089, Warsaw, Poland
| | - Dylan Jayatilaka
- School of Molecular Sciences, The University of Western Australia, 35 Stirling Highway, Perth, WA, 6009, Australia
| | - Birger Dittrich
- Heinrich-Heine-Universität Düsseldorf, Anorganische Chemie und Strukturchemie, Gebäude 26.42.01.21, Universitätsstraße 1, 40225, Düsseldorf, Germany
| | - Ralf Flaig
- Diamond Light Source, Harwell Science and Innovation Campus, Didcot, Oxon, OX11 0DE, UK
| | - Peter Luger
- Freie Universität Berlin, Institut für Chemie und Biochemie, Fabeckstraße 36a, 14195, Berlin, Germany
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Zwirki i Wigury 101, 02-089, Warsaw, Poland
| | - Paulina M Dominiak
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Zwirki i Wigury 101, 02-089, Warsaw, Poland
| | - Simon Grabowsky
- Institut für Anorganische Chemie und Kristallographie, Fachbereich 2-Biologie/Chemie, Universität Bremen, Leobener Straße NW2, 28359, Bremen, Germany
| |
Collapse
|
11
|
Daszkiewicz M, Mielcarek A. Comparison of hydrogen bonds and diverse weak interactions of the nitro group in 2-methyl-4-nitroanilinium nitrate, bisulfate and two hexafluoridosilicates: elementary graph-set approach. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2016; 72:916-926. [PMID: 27910842 DOI: 10.1107/s2052520616016139] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2016] [Accepted: 10/11/2016] [Indexed: 06/06/2023]
Abstract
Crystal structures of (H2m4na)NO3 (1), (H2m4na)HSO4 (2), (H2m4na)2SiF6 (3) and (H2m4na)2SiF6·2H2O (4), where 2m4na = 2-methyl-4-nitroaniline, are presented. Two layers of interactions occur in the structures, N-H...O/F hydrogen bonds and interactions with the nitro group. Although diverse, hydrogen-bonding patterns are compared with each other by means of interrelations among elementary graph-set descriptors and descriptors of hydrogen-bonding patterns. Using mathematical relations, the gradual expansion of the ring patterns was shown in the crystal structures. Parallel and perpendicular arranged nitro groups form weak π(N)NO2...π(O)NO2 and π(N)NO2...ONO2 interactions, respectively. The πNO2...πring interaction has an impact to the stabilization of parallel oriented nitro groups. Generally, weak interactions constructed by the nitro group occur in the studied crystals as follows: π(N)NO2...π(O)NO2, πring...πring, C-H...O (1); π(N)NO2...π(O)NO2, π(N)NO2...ONO2 (2); π(N)NO2...π(O)NO2, π(N)NO2...ONO2 (3); C-H...O (4).
Collapse
Affiliation(s)
- Marek Daszkiewicz
- Institute of Low Temperature and Structure Research, Polish Academy of Sciences, PO Box 1410, Wroclaw 50-950, Poland
| | - Agnieszka Mielcarek
- Department of Inorganic Chemistry, Gdańsk University of Technology, Narutowicza str. 11/12, Gdańsk 80-233, Poland
| |
Collapse
|
12
|
Kalaiarasi C, Pavan MS, Kumaradhas P. Topological characterization of electron density, electrostatic potential and intermolecular interactions of 2-nitroimidazole: an experimental and theoretical study. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2016; 72:775-786. [PMID: 27698320 DOI: 10.1107/s2052520616010581] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2016] [Accepted: 06/29/2016] [Indexed: 06/06/2023]
Abstract
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 using CRYSTAL09.
Collapse
Affiliation(s)
- Chinnasamy Kalaiarasi
- Laboratory of Biocrystallography and Computational Molecular Biology, Department of Physics, Periyar University, Salem 636 011, India
| | - Mysore S Pavan
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560 012, India
| | - Poomani Kumaradhas
- Laboratory of Biocrystallography and Computational Molecular Biology, Department of Physics, Periyar University, Salem 636 011, India
| |
Collapse
|
13
|
Electrostatic Properties of N-Heterocyclic Carbenes Obtained by Experimental Charge-Density Analysis of Two Selenium Adducts. Eur J Inorg Chem 2016. [DOI: 10.1002/ejic.201600292] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
14
|
Santos OL, Fonseca TL, Sabino JR, Georg HC, Castro MA. Polarization effects on the electric properties of urea and thiourea molecules in solid phase. J Chem Phys 2016; 143:234503. [PMID: 26696062 DOI: 10.1063/1.4937481] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present theoretical results for the dipole moment, linear polarizability, and first hyperpolarizability of the urea and thiourea molecules in solid phase. The in-crystal electric properties were determined by applying a supermolecule approach in combination with an iterative electrostatic scheme, in which the surrounding molecules are represented by point charges. It is found for both urea and thiourea molecules that the influence of the polarization effects is mild for the linear polarizability, but it is marked for the dipole moment and first hyperpolarizability. The replacement of oxygen atoms by sulfur atoms increases, in general, the electric responses. Our second-order Møller-Plesset perturbation theory based iterative scheme predicts for the in-crystal dipole moment of urea and thiourea the values of 7.54 and 9.19 D which are, respectively, increased by 61% and 58%, in comparison with the corresponding isolated values. The result for urea is in agreement with the available experimental result of 6.56 D. In addition, we present an estimate of macroscopic quantities considering explicit unit cells of urea and thiourea crystals including environment polarization effects. These supermolecule calculations take into account partially the exchange and dispersion effects. The results illustrate the role played by the electrostatic interactions on the static second-order nonlinear susceptibility of the urea crystal.
Collapse
Affiliation(s)
- O L Santos
- Instituto de Física, Universidade Federal de Goiás, Campus Samambaia, 74.690-900 Goiânia, GO, Brazil
| | - T L Fonseca
- Instituto de Física, Universidade Federal de Goiás, Campus Samambaia, 74.690-900 Goiânia, GO, Brazil
| | - J R Sabino
- Instituto de Física, Universidade Federal de Goiás, Campus Samambaia, 74.690-900 Goiânia, GO, Brazil
| | - H C Georg
- Instituto de Física, Universidade Federal de Goiás, Campus Samambaia, 74.690-900 Goiânia, GO, Brazil
| | - M A Castro
- Instituto de Física, Universidade Federal de Goiás, Campus Samambaia, 74.690-900 Goiânia, GO, Brazil
| |
Collapse
|
15
|
|
16
|
Sovago I, Gutmann MJ, Senn HM, Thomas LH, Wilson CC, Farrugia LJ. Electron density, disorder and polymorphism: high-resolution diffraction studies of the highly polymorphic neuralgic drug carbamazepine. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2016; 72:39-50. [PMID: 26830795 DOI: 10.1107/s2052520615019538] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2015] [Accepted: 10/15/2015] [Indexed: 06/05/2023]
Abstract
Analysis of neutron and high-resolution X-ray diffraction data on form (III) of carbamazepine at 100 K using the atoms in molecules (AIM) topological approach afforded excellent agreement between the experimental results and theoretical densities from the optimized gas-phase structure and from multipole modelling of static theoretical structure factors. The charge density analysis provides experimental confirmation of the partially localized π-bonding suggested by the conventional structural formula, but the evidence for any significant C-N π bonding is not strong. Hirshfeld atom refinement (HAR) gives H atom positional and anisotropic displacement parameters that agree very well with the neutron parameters. X-ray and neutron diffraction data on the dihydrate of carbemazepine strongly indicate a disordered orthorhombic crystal structure in the space group Cmca, rather than a monoclinic crystal structure in space group P2(1)/c. This disorder in the dihydrate structure has implications for both experimental and theoretical studies of polymorphism.
Collapse
Affiliation(s)
- Ioana Sovago
- WESTChem School of Chemistry, University of Glasgow, Glasgow G12 8QQ, Scotland
| | - Matthias J Gutmann
- ISIS Facility, STFC Rutherford Appleton Laboratory, Harwell Science and Innovation Campus, Chilton, Didcot, Oxfordshire OX11 0QX, England
| | - Hans Martin Senn
- WESTChem School of Chemistry, University of Glasgow, Glasgow G12 8QQ, Scotland
| | - Lynne H Thomas
- Department of Chemistry, University of Bath, Claverton Down, Bath BA2 7AY, England
| | - Chick C Wilson
- Department of Chemistry, University of Bath, Claverton Down, Bath BA2 7AY, England
| | - Louis J Farrugia
- WESTChem School of Chemistry, University of Glasgow, Glasgow G12 8QQ, Scotland
| |
Collapse
|
17
|
Gryl M. Charge density and optical properties of multicomponent crystals containing active pharmaceutical ingredients or their analogues. ACTA CRYSTALLOGRAPHICA SECTION B-STRUCTURAL SCIENCE CRYSTAL ENGINEERING AND MATERIALS 2015. [DOI: 10.1107/s2052520615013505] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Active pharmaceutical ingredients (APIs), through their favourable donor/acceptor spatial distribution and synthon formation flexibility, are attractive building blocks in modern materials crystallography. The optical properties of a crystal strongly depend on two factors,i.e.the spatial distribution of molecules in the crystal structure and the electronic properties of molecular building blocks (dipole moments, polarizabilities, hyperpolarizabilities). Although the latter are easy to predict throughab initiocalculations, the former are not. Only a combination of experimental and theoretical charge density studies together with prediction and measurement of optical properties enable full analysis of the obtained functional material in terms of its usefulness in practical applications. This article presents design strategies of optical materials based on selected pharmaceutical molecules. Factors that contribute to molecular recognition in the four selected polar/chiral crystal phases (derived through charge density and Hirshfeld surfaces analysis) have been determined. Theoretically predicted optical properties of the molecular/ionic building blocks as well as bulk effects have been confirmed experimentally. This research is a first step in the design of novel optical materials based on push–pull molecules and APIs.
Collapse
|
18
|
Rawat P, Singh RN. Assessment of conformational, spectral, antimicrobial activity, chemical reactivity and NLO application of Pyrrole-2,5-dicarboxaldehyde bis(oxaloyldihydrazone). SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2015; 140:344-355. [PMID: 25615679 DOI: 10.1016/j.saa.2014.12.080] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2014] [Revised: 12/15/2014] [Accepted: 12/17/2014] [Indexed: 06/04/2023]
Abstract
An orange colored pyrrole dihydrazone: Pyrrole-2,5-dicarboxaldehyde bis(oxaloyldihydrazone) (PDBO) has been synthesized by reaction of oxalic acid dihydrazide with 2,5 diformyl-1H-pyrrole and has been characterized by spectroscopic analysis (1H, 13C NMR, UV-visible, FT-IR and DART Mass). The properties of the compound has been evaluated using B3LYP functional and 6-31G(d,p)/6-311+G(d,p) basis set. The symmetric (3319, 3320 cm(-1)) and asymmetric (3389, 3382 cm(-1)) stretching wave number confirm free NH2 groups in PDBO. NBO analysis shows, inter/intra molecular interactions within the molecule. Topological parameters have been analyzed by QTAIM theory and provide the existence of intramolecular hydrogen bonding (N-H⋯O). The local reactivity descriptors analyses determine the reactive sites within molecule. The calculated first hyperpolarizability value (β0=23.83×10(-30) esu) of pyrrole dihydrazone shows its suitability for non-linear optical (NLO) response. The preliminary bioassay suggested that the PDBO exhibits relatively good antibacterial and fungicidal activity against Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus, Streptococcus pyogenes, Candida albicans, Aspergillus niger. The local reactivity descriptors--Fukui functions (fk+, fk-), local softnesses (sk+, sk-) and electrophilicity indices (ωk+, ωk-) analyses have been used to determine the reactive sites within molecule.
Collapse
Affiliation(s)
- Poonam Rawat
- Department of Chemistry, University of Lucknow, Lucknow 226007, U.P., India
| | - R N Singh
- Department of Chemistry, University of Lucknow, Lucknow 226007, U.P., India.
| |
Collapse
|
19
|
Mermer A, Lis T, Starynowicz P. Charge density distribution and theoretical analysis of low and high energy phosphate esters. RSC Adv 2015. [DOI: 10.1039/c5ra12605c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
There is a strict relation between the energy of hydrolysis of phosphate esters and the extent of interactions between the p ester oxygen lone pair and the antibonding orbitals of the rest of the molecule. Its impact on experimental charge density distribution is analyzed.
Collapse
Affiliation(s)
- Adrian Mermer
- Wydział Chemii
- Uniwersytet Wrocławski
- 50-383 Wrocław
- Poland
| | - Tadeusz Lis
- Wydział Chemii
- Uniwersytet Wrocławski
- 50-383 Wrocław
- Poland
| | | |
Collapse
|
20
|
Seidler T, Champagne B. Which charge definition for describing the crystal polarizing field and the χ(1) and χ(2) of organic crystals? Phys Chem Chem Phys 2015; 17:19546-56. [DOI: 10.1039/c5cp03248b] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Crystal optical susceptibilities are probes to assess the performance of the charge definition employed to describe the crystal polarizing field.
Collapse
Affiliation(s)
- Tomasz Seidler
- Laboratoire de Chimie Théorique
- University of Namur
- Namur
- Belgium
- Faculty of Chemistry
| | | |
Collapse
|
21
|
Madsen AØ, Hoser AA. SHADE3server: a streamlined approach to estimate H-atom anisotropic displacement parameters using periodicab initiocalculations or experimental information. J Appl Crystallogr 2014. [DOI: 10.1107/s1600576714022973] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
A major update of theSHADEserver (http://shade.ki.ku.dk) is presented. In addition to all of the previous options for estimating H-atom anisotropic displacement parameters (ADPs) that were offered bySHADE2, the newest version offers two new methods. The first method combines the original translation–libration–screw analysis with input from periodicab initiocalculations. The second method allows the user to input experimental information from spectroscopic measurements or from neutron diffraction experiments on related structures and utilize this information to evaluate ADPs of H atoms. Tools are provided to set up theab initiocalculations and to derive the internal motion from the calculations. The new server was tested on a range of compounds where neutron diffraction data were available. In most cases, the results are significantly better than previous estimates, and for strong hydrogen bonds in proton sponges, theab initiocalculations become crucial.
Collapse
|
22
|
Starynowicz P, Lis T. Charge-density distribution in sodium bis(4-nitrophenyl)phosphate. ACTA CRYSTALLOGRAPHICA SECTION B-STRUCTURAL SCIENCE CRYSTAL ENGINEERING AND MATERIALS 2014; 70:723-31. [DOI: 10.1107/s2052520614010439] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/20/2014] [Accepted: 05/07/2014] [Indexed: 11/10/2022]
Abstract
The electron-density distribution in sodium bis(4-nitrophenyl)phosphate has been analyzed using the multipole refinement of X-ray diffraction data and of theoretical density-functional theory (DFT) calculations. The ester P—O bonds are particularly long and their topological parameters (density at the bond critical point, Laplacian) are lower than for other P—O bonds. Some disagreement between the experimental and theoretical charges of atoms constituting the nitro groups has been observed and the possible reasons are discussed. Weak polarization effects produced by sodium cations may be observed within the phosphate fragment; they are more manifest in the case of the nitro groups.
Collapse
|
23
|
Seidler T, Stadnicka K, Champagne B. Evaluation of the Linear and Second-Order NLO Properties of Molecular Crystals within the Local Field Theory: Electron Correlation Effects, Choice of XC Functional, ZPVA Contributions, and Impact of the Geometry in the Case of 2-Methyl-4-nitroaniline. J Chem Theory Comput 2014; 10:2114-24. [DOI: 10.1021/ct5001654] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tomasz Seidler
- Faculty
of Chemistry, Jagiellonian University, ul, Ingardena 3, 30-060, Kraków, Poland
- Laboratoire
de Chimie Théorique, Unité de Chimie Physique Théorique
et Structurale, University of Namur, rue de Bruxelles, 61, B-5000, Namur, Belgium
| | - Katarzyna Stadnicka
- Faculty
of Chemistry, Jagiellonian University, ul, Ingardena 3, 30-060, Kraków, Poland
| | - Benoît Champagne
- Laboratoire
de Chimie Théorique, Unité de Chimie Physique Théorique
et Structurale, University of Namur, rue de Bruxelles, 61, B-5000, Namur, Belgium
| |
Collapse
|
24
|
Sovago I, Gutmann M, Hill JG, Senn HM, Thomas LH, Wilson CC, Farrugia LJ. Experimental Electron Density and Neutron Diffraction Studies on the Polymorphs of Sulfathiazole. CRYSTAL GROWTH & DESIGN 2014; 14:1227-1239. [PMID: 24672285 PMCID: PMC3963452 DOI: 10.1021/cg401757z] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/22/2013] [Revised: 01/13/2014] [Indexed: 05/31/2023]
Abstract
High resolution X-ray diffraction data on forms I-IV of sulfathiazole and neutron diffraction data on forms II-IV have been collected at 100 K and analyzed using the Atoms in Molecules topological approach. The molecular thermal motion as judged by the anisotropic displacement parameters (adp's) is very similar in all four forms. The adp of the thiazole sulfur atom had the greatest amplitude perpendicular to the five-membered ring, and analysis of the temperature dependence of the adps indicates that this is due to genuine thermal motion rather than a concealed disorder. A minor disorder (∼1-2%) is evident for forms I and II, but a statistical analysis reveals no deleterious effect on the derived multipole populations. The topological analysis reveals an intramolecular S-O···S interaction, which is consistently present in all experimental topologies. Analysis of the gas-phase conformation of the molecule indicates two low-energy theoretical conformers, one of which possesses the same intramolecular S-O···S interaction observed in the experimental studies and the other an S-O···H-N intermolecular interaction. These two interactions appear responsible for "locking" the molecular conformation. The lattice energies of the various polymorphs computed from the experimental multipole populations are highly dependent on the exact refinement model. They are similar in magnitude to theoretically derived lattice energies, but the relatively high estimated errors mean that this method is insufficiently accurate to allow a definitive stability order for the sulfathiazole polymorphs at 0 K to be determined.
Collapse
Affiliation(s)
- Ioana Sovago
- WESTChem
School of Chemistry, University of Glasgow, Glasgow G12 8QQ, U.K.
| | - Matthias
J. Gutmann
- ISIS
Facility, STFC Rutherford Appleton Laboratory, Harwell Science and Innovation Campus, Chilton,
Didcot, Oxfordshire OX11 0QX, U.K.
| | - J. Grant Hill
- WESTChem
School of Chemistry, University of Glasgow, Glasgow G12 8QQ, U.K.
| | - Hans Martin Senn
- WESTChem
School of Chemistry, University of Glasgow, Glasgow G12 8QQ, U.K.
| | - Lynne H. Thomas
- Department
of Chemistry, University of Bath, Claverton Down, Bath BA2 7AY, U.K.
| | - Chick C. Wilson
- Department
of Chemistry, University of Bath, Claverton Down, Bath BA2 7AY, U.K.
| | - Louis J. Farrugia
- WESTChem
School of Chemistry, University of Glasgow, Glasgow G12 8QQ, U.K.
| |
Collapse
|
25
|
Seidler T, Stadnicka K, Champagne B. Investigation of the linear and second-order nonlinear optical properties of molecular crystals within the local field theory. J Chem Phys 2013; 139:114105. [DOI: 10.1063/1.4819769] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Pozzi CG, Fantoni AC, Goeta AE, de Matos Gomes E, McIntyre GJ, Punte G. On the use of crystal vibrational modes in the estimation of the anisotropic displacement parameters of hydrogen atoms in molecular crystals: para-Nitroaniline as a test case. Chem Phys 2013. [DOI: 10.1016/j.chemphys.2013.07.004] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
27
|
Herbst-Irmer R, Henn J, Holstein JJ, Hübschle CB, Dittrich B, Stern D, Kratzert D, Stalke D. Anharmonic motion in experimental charge density investigations. J Phys Chem A 2013; 117:633-41. [PMID: 23241030 DOI: 10.1021/jp309985e] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
In the charge density study of 9-diphenylthiophosphinoylanthracene the thermal motion of several atoms needed an anharmonic description via Gram-Charlier coefficients even for data collected at 15 K. As several data sets at different temperatures were measured, this anharmonic model could be proved to be superior to a disorder model. Refinements against theoretical data showed the resemblance of an anharmonic model and a disorder model with two positions very close to each other (~0.2 Å), whereas these two models could be clearly distinguished if the second position is 0.5 Å apart. The refined multipole parameters were distorted when the anharmonic motion was not properly refined. Therefore, this study reveals the importance of detecting and properly handling anharmonic motion. Unrefined anharmonic motion leads to typical shashlik-like residual density patterns. Therefore, careful analysis of the residual density and the derived probability density function after the refinement of the Gram-Charlier coefficients proved to be the most useful tools to indicate the presence of anharmonic motion.
Collapse
Affiliation(s)
- Regine Herbst-Irmer
- Institut für Anorganische Chemie, Universität Göttingen, Tammannstrasse 4, 37077 Göttingen, Germany.
| | | | | | | | | | | | | | | |
Collapse
|
28
|
Yang J, Waller MP. Revealing noncovalent interactions in quantum crystallography: Taurine revisited. J Comput Chem 2012; 34:466-70. [DOI: 10.1002/jcc.23155] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2012] [Revised: 08/10/2012] [Accepted: 08/30/2012] [Indexed: 11/07/2022]
|
29
|
Dittrich B, Jayatilaka D. Reliable Measurements of Dipole Moments from Single-Crystal Diffraction Data and Assessment of an In-Crystal Enhancement. ELECTRON DENSITY AND CHEMICAL BONDING II 2012. [DOI: 10.1007/430_2012_78] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
30
|
Szostak MM, Chojnacki H, Piela K, Okwieka-Lupa U, Bidzińska E, Dyrek K. Helical Superstructure and Charged Polarons Contributions to Optical Nonlinearity of 2-Methyl-4-nitroaniline Crystals Studied by Resonance Raman, Electron Paramagnetic Resonance, Circular Dichroism Spectroscopies, and Quantum Chemical Calculations. J Phys Chem A 2011; 115:7448-55. [DOI: 10.1021/jp204168s] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- M. Magdalena Szostak
- Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wybrzeże Wyspiańskieg 27, PL-50-370 Wrocław, Poland
| | - Henryk Chojnacki
- Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wybrzeże Wyspiańskieg 27, PL-50-370 Wrocław, Poland
| | - Katarzyna Piela
- Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wybrzeże Wyspiańskieg 27, PL-50-370 Wrocław, Poland
| | - Urszula Okwieka-Lupa
- Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wybrzeże Wyspiańskieg 27, PL-50-370 Wrocław, Poland
| | - Ewa Bidzińska
- Faculty of Chemistry, Jagiellonian University, ul. Ingardena 3, PL-30-060 Kraków, Poland
| | - Krystyna Dyrek
- Faculty of Chemistry, Jagiellonian University, ul. Ingardena 3, PL-30-060 Kraków, Poland
| |
Collapse
|
31
|
Farrugia LJ, Khalaji AD. Evidence for Side-Chain π-Delocalization in a Planar Substituted Benzene: An Experimental and Theoretical Charge Density Study on 2,5-Dimethoxybenzaldehyde Thiosemicarbazone. J Phys Chem A 2011; 115:12512-22. [DOI: 10.1021/jp2026169] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Louis J. Farrugia
- WestCHEM, School of Chemistry, University of Glasgow, Glasgow G12 8QQ, Scotland
| | | |
Collapse
|
32
|
Stephen AD, Thomas R, Srinivasan P, Narayayanasamy V, Kumaradhas P. Exploring the bond topological and electrostatic properties of benzimidazole molecule via experimental and theoretical charge density study. J Mol Struct 2011. [DOI: 10.1016/j.molstruc.2011.01.014] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
33
|
Farrugia LJ, Senn HM. Metal-metal and metal-ligand bonding at a QTAIM catastrophe: a combined experimental and theoretical charge density study on the alkylidyne cluster Fe3(μ-H)(μ-COMe)(CO)10. J Phys Chem A 2010; 114:13418-33. [PMID: 21182291 DOI: 10.1021/jp1098624] [Citation(s) in RCA: 79] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The charge density in the tri-iron methoxymethylidyne cluster Fe(3)(μ-H)(μ-COMe)(CO)(10) (1) has been studied experimentally at 100 K and by DFT calculations on the isolated molecule using the Quantum Theory of Atoms in Molecules (QTAIM). The COMe ligand acts as a nearly symmetric bridge toward two of the Fe atoms (Fe-C = 1.8554(4), 1.8608(4) Å) but with a much longer interaction to the third Fe atom, Fe-C = 2.6762(4) Å. Complex 1 provides a classic example where topological QTAIM catastrophes render an exact structure description ambiguous. While all experimental and theoretical studies agree in finding no direct metal-metal interaction for the doubly bridged Fe-Fe vector, the chemical bonding between the Fe(CO)(4) unit and the Fe(2)(μ-H)(μ-COMe)(CO)(6) moiety in terms of conventional QTAIM descriptors is much less clear. Bond paths implying direct Fe-Fe interactions and a weak interaction between the COMe ligand and the Fe(CO)(4) center are observed, depending on the experimental or theoretical density model examined. Theoretical studies using the Electron Localizability Indicator (ELI-D) suggest the metal-metal bonding is more significant, while the delocalization indices imply that both Fe-Fe bonding and Fe···C(alkylidyne) bonding are equally important. The source functions at various interfragment reference points are similar and highly delocalized. The potential-energy surface (PES) for the migration of the alkylidyne group from a μ(2) to a semi-μ(3) coordination mode has been explored by DFT calculations on 1 and the model complexes M(3)(μ-H)(μ-CH)(CO)(10) (M = Fe, 2; Ru, 3; and Os, 4). These calculations confirm a semi-μ(3) bridging mode for the alkylidyne ligand as the minimum-energy geometry for compounds 2-4 and demonstrate that, for 1, both Fe-Fe and Fe···C(alkylidyne) interactions are important in the cluster bonding. The PES between μ(2) and semi-μ(3) alkylidyne coordination for 1 is extremely soft, and the interconversion between several topological isomers is predicted to occur with almost no energy cost. Analysis of the density ρ(r) and the Laplacian of the density ▽(2)ρ(r(b)) in the methoxymethylidyne ligand is consistent with a partial π-bond character of the C-O bond, associated with an sp(2) hybridization for these atoms.
Collapse
Affiliation(s)
- Louis J Farrugia
- WestCHEM, School of Chemistry, University of Glasgow, Glasgow G12 8QQ, Scotland
| | | |
Collapse
|
34
|
Meindl K, Herbst-Irmer R, Henn J. On the effect of neglecting anharmonic nuclear motion in charge density studies. Acta Crystallogr A 2010; 66:362-71. [PMID: 20404442 DOI: 10.1107/s0108767310006343] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2009] [Accepted: 02/17/2010] [Indexed: 11/10/2022] Open
Abstract
The effect of neglecting anharmonic nuclear motion when it is definitely present is studied. To ensure the presence of anharmonic nuclear motion a model was used that was previously refined against experimental data including anharmonic nuclear motion, and these calculated structure factors were used as observed data for a multipole refinement. It was then studied how the neglect of anharmonic nuclear motion and noise in the data affects the usual crystallographic quality measure R, the density parameters and the residual density distribution. It is demonstrated that the neglect of anharmonic nuclear motion leads to a characteristic imprint onto the residual density distribution in terms of residual density peaks and holes, in terms of the whole residual density distribution and in terms of the number, location and strength of valence shell charge concentrations (VSCCs). These VSCCs differ from that of the input model in a way which heavily influences and misleads the chemical interpretation of the charge density. This imprint vanishes after taking anharmonic nuclear motion into account. Also the input model VSCCs are restored. The importance of modeling anharmonic nuclear motion is furthermore shown by the characteristic imprint on the residual density distribution, even in the case of a numerically almost unaffected R value.
Collapse
Affiliation(s)
- Kathrin Meindl
- Georg-August Universität Göttingen, Tammannstrasse 4, 37077 Göttingen, Germany
| | | | | |
Collapse
|
35
|
Okwieka U, Hołderna-Natkaniec K, Misiaszek T, Medycki W, Baran J, Magdalena Szostak M. Dynamical disorder in 2-methyl-4-nitroaniline and its deuterated analogue crystals studied by Fourier transform infrared and nuclear magnetic resonance. J Chem Phys 2009; 131:144505. [DOI: 10.1063/1.3243855] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
36
|
Netzel J, van Smaalen S. Topological properties of hydrogen bonds and covalent bonds from charge densities obtained by the maximum entropy method (MEM). ACTA CRYSTALLOGRAPHICA. SECTION B, STRUCTURAL SCIENCE 2009; 65:624-38. [PMID: 19767685 PMCID: PMC2749645 DOI: 10.1107/s0108768109026767] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/05/2009] [Accepted: 07/08/2009] [Indexed: 11/10/2022]
Abstract
Charge densities have been determined by the Maximum Entropy Method (MEM) from the high-resolution, low-temperature (T approximately 20 K) X-ray diffraction data of six different crystals of amino acids and peptides. A comparison of dynamic deformation densities of the MEM with static and dynamic deformation densities of multipole models shows that the MEM may lead to a better description of the electron density in hydrogen bonds in cases where the multipole model has been restricted to isotropic displacement parameters and low-order multipoles (l(max) = 1) for the H atoms. Topological properties at bond critical points (BCPs) are found to depend systematically on the bond length, but with different functions for covalent C-C, C-N and C-O bonds, and for hydrogen bonds together with covalent C-H and N-H bonds. Similar dependencies are known for AIM properties derived from static multipole densities. The ratio of potential and kinetic energy densities |V(BCP)|/G(BCP) is successfully used for a classification of hydrogen bonds according to their distance d(H...O) between the H atom and the acceptor atom. The classification based on MEM densities coincides with the usual classification of hydrogen bonds as strong, intermediate and weak [Jeffrey (1997). An Introduction to Hydrogen Bonding. Oxford University Press]. MEM and procrystal densities lead to similar values of the densities at the BCPs of hydrogen bonds, but differences are shown to prevail, such that it is found that only the true charge density, represented by MEM densities, the multipole model or some other method can lead to the correct characterization of chemical bonding. Our results do not confirm suggestions in the literature that the promolecule density might be sufficient for a characterization of hydrogen bonds.
Collapse
Affiliation(s)
- Jeanette Netzel
- Laboratory of Crystallography, University of Bayreuth, D-95440 Bayreuth, Germany
| | - Sander van Smaalen
- Laboratory of Crystallography, University of Bayreuth, D-95440 Bayreuth, Germany
| |
Collapse
|
37
|
Fantoni AC, Pozzi CG, Punte G. On the Electron Density Topology and Electrostatic Properties of Nitroanilines. A Theoretical Investigation on m-Nitroaniline and 2-Methyl-5-nitroaniline Crystals. J Phys Chem A 2009; 113:9527-32. [DOI: 10.1021/jp904318q] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Adolfo C. Fantoni
- IFLP (CCT-La Plata), Departamento de Física, Facultad de Ciencias Exactas, UNLP CC 67, 1900 La Plata, Argentina
| | - C. Gustavo Pozzi
- IFLP (CCT-La Plata), Departamento de Física, Facultad de Ciencias Exactas, UNLP CC 67, 1900 La Plata, Argentina
| | - Graciela Punte
- IFLP (CCT-La Plata), Departamento de Física, Facultad de Ciencias Exactas, UNLP CC 67, 1900 La Plata, Argentina
| |
Collapse
|
38
|
|
39
|
Farrugia LJ, Evans C, Lentz D, Roemer M. The QTAIM Approach to Chemical Bonding Between Transition Metals and Carbocyclic Rings: A Combined Experimental and Theoretical Study of (η5-C5H5)Mn(CO)3, (η6-C6H6)Cr(CO)3, and (E)-{(η5-C5H4)CF═CF(η5-C5H4)}(η5-C5H5)2Fe2. J Am Chem Soc 2008; 131:1251-68. [DOI: 10.1021/ja808303j] [Citation(s) in RCA: 142] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Louis J. Farrugia
- WestCHEM, Department of Chemistry, University of Glasgow, Glasgow G12 8QQ, Scotland, U.K., and Institut für Chemie - Anorganische und Analytische Chemie, Fachbereich Biologie, Chemie, Pharmazie, Freie Universität, D-14195 Berlin, Germany
| | - Cameron Evans
- WestCHEM, Department of Chemistry, University of Glasgow, Glasgow G12 8QQ, Scotland, U.K., and Institut für Chemie - Anorganische und Analytische Chemie, Fachbereich Biologie, Chemie, Pharmazie, Freie Universität, D-14195 Berlin, Germany
| | - Dieter Lentz
- WestCHEM, Department of Chemistry, University of Glasgow, Glasgow G12 8QQ, Scotland, U.K., and Institut für Chemie - Anorganische und Analytische Chemie, Fachbereich Biologie, Chemie, Pharmazie, Freie Universität, D-14195 Berlin, Germany
| | - Max Roemer
- WestCHEM, Department of Chemistry, University of Glasgow, Glasgow G12 8QQ, Scotland, U.K., and Institut für Chemie - Anorganische und Analytische Chemie, Fachbereich Biologie, Chemie, Pharmazie, Freie Universität, D-14195 Berlin, Germany
| |
Collapse
|
40
|
Farrugia LJ, Middlemiss DS, Sillanpää R, Seppälä P. A Combined Experimental and Theoretical Charge Density Study of the Chemical Bonding and Magnetism in 3-Amino-propanolato Cu(II) Complexes Containing Weakly Coordinated Anions. J Phys Chem A 2008; 112:9050-67. [DOI: 10.1021/jp804865j] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Louis J. Farrugia
- WestCHEM, Department of Chemistry, University of Glasgow, Glasgow G12 8QQ Scotland U.K., and Department of Chemistry, P.O. Box 35, FIN-40014, University of Jyväskylä, Jyväskylä, Finland
| | - Derek S. Middlemiss
- WestCHEM, Department of Chemistry, University of Glasgow, Glasgow G12 8QQ Scotland U.K., and Department of Chemistry, P.O. Box 35, FIN-40014, University of Jyväskylä, Jyväskylä, Finland
| | - Reijo Sillanpää
- WestCHEM, Department of Chemistry, University of Glasgow, Glasgow G12 8QQ Scotland U.K., and Department of Chemistry, P.O. Box 35, FIN-40014, University of Jyväskylä, Jyväskylä, Finland
| | - Petri Seppälä
- WestCHEM, Department of Chemistry, University of Glasgow, Glasgow G12 8QQ Scotland U.K., and Department of Chemistry, P.O. Box 35, FIN-40014, University of Jyväskylä, Jyväskylä, Finland
| |
Collapse
|
41
|
Qiu YQ, Liu XD, Sun SL, Fan M, Su ZM, Wang RS. DFT study on second-order nonlinear optical properties of the derivatives of 7-vertex cobalt–carborane metallocenyl. ACTA ACUST UNITED AC 2008. [DOI: 10.1016/j.theochem.2008.05.021] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
42
|
Destro R, Soave R, Barzaghi M. Physicochemical properties of zwitterionic L- and DL-alanine crystals from their experimental and theoretical charge densities. J Phys Chem B 2008; 112:5163-74. [PMID: 18373373 DOI: 10.1021/jp710496q] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The total experimental electron density distributions rho(r) of zwitterionic L- and DL-alanine crystals, as derived from extensive sets of X-ray diffracted intensities collected at 23 and 19 K, are compared to gain an insight into the different physical properties of the two related chiral compounds in the solid state and to explore the extent of the rho(r) transferability. Relevant parameters that characterize the two crystal forms are obtained, showing differences and similarities in terms of (i) geometric descriptors, (ii) topological indexes, (iii) molecular electrostatic potential Phi(r) distributions, (iv) atomic volumes and charges, (v) molecular electric moments, and (vi) electrostatic interaction energies. To assess the relative stability of the racemate with respect to the pure enantiomer, the crystal lattice energies, as obtained through DFT fully periodic calculations, are also discussed and compared with the experimental sublimation enthalpies after correction for the proton-transfer energies. In-crystal group charges, evaluated with the quantum theory of atoms in molecules, are found to be transferable between the racemic and the pure enantiomer, at variance with group volumes. Similarly, molecular first and third moments are not strictly transferable and indicate that for the zwitterionic alanine molecule the molecular charge distribution in the DL-crystal is more polarized in the c direction by about 10%. By contrast, quantitative agreement is observed for second and fourth moments. Significant differences arise from (1) the crystal packing of the dipole vectors, which are aligned in an antiparallel fashion in the L-crystal, to be compared with a parallel alignment in the racemate, due the polar space group Pna21 of the latter, (2) the strongly attractive electrostatic energy of a homochiral pair in the L-crystal, which is opposed to the corresponding heterochiral pair in the DL-crystal form. The difference between these Ees values amounts to 135-150 kJ mol(-1). Despite this, the two crystal forms are predicted as equally thermodynamically favored by the theoretical P-B3LYP estimates of the crystal lattice energies. Finally, the necessity of an upgrading of the dispersion and exchange-repulsion terms currently adopted within the experimental charge density approach to intermolecular interactions is recognized and discussed.
Collapse
Affiliation(s)
- Riccardo Destro
- Dipartimento di Chimica Fisica ed Elettrochimica, Università di Milano, Via Golgi 19, 20133 Milano, Italy
| | | | | |
Collapse
|
43
|
Blomberg E, Verrall R, Claesson PM. Interactions between adsorbed layers of cationic gemini surfactants. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2008; 24:1133-1140. [PMID: 18052227 DOI: 10.1021/la702940p] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
The forces acting between glass and between mica surfaces in the presence of two cationic gemini surfactants, 1,4 diDDAB (1,4-butyl-bis(dimethyldodecylammonium bromide)) and 1,12 diDDAB (1,12-dodecyl-bis(dimethyldodecylammonium bromide)), have been investigated below the critical micelle concentration (cmc) of the surfactants using two different surface force techniques. In both cases, it was found that a recharging of the surfaces occurred at a surfactant concentration of about 0.1 x cmc, and at all surfactant concentrations investigated repulsive double-layer forces dominated the interaction at large separations. At smaller separations, attractive forces, or regions of separation with (close to) constant force, were observed. This was interpreted as being due to desorption and rearrangement in the adsorbed layer induced by the proximity of a second surface. Analysis of the decay length of the repulsive double-layer force showed that the majority of the gemini surfactants were fully dissociated. However, the degree of ion pair formation, between a gemini surfactant and a bromide counterion, increased with increasing surfactant concentration and was larger for the gemini surfactant with a shorter spacer length.
Collapse
Affiliation(s)
- Eva Blomberg
- Department of Chemistry, Surface Chemistry, Royal Institute of Technology, Drottning Kristinas väg 51, Stockholm, Sweden.
| | | | | |
Collapse
|
44
|
Tummala NR, Striolo A. Role of Counterion Condensation in the Self-Assembly of SDS Surfactants at the Water−Graphite Interface. J Phys Chem B 2008; 112:1987-2000. [DOI: 10.1021/jp077678m] [Citation(s) in RCA: 88] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Naga Rajesh Tummala
- The University of Oklahoma School of Chemical Biological and Materials Engineering, Norman, Oklahoma 73019
| | - Alberto Striolo
- The University of Oklahoma School of Chemical Biological and Materials Engineering, Norman, Oklahoma 73019
| |
Collapse
|
45
|
Soave R, Barzaghi M, Destro R. Progress in the understanding of drug-receptor interactions, part 2: experimental and theoretical electrostatic moments and interaction energies of an angiotensin II receptor antagonist (C30H30N6(O)3S). Chemistry 2007; 13:6942-56. [PMID: 17539033 DOI: 10.1002/chem.200601516] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
A combined experimental and theoretical charge density study of an angiotensin II receptor antagonist (1) is presented focusing on electrostatic properties such as atomic charges, molecular electric moments up to the fourth rank and energies of the intermolecular interactions, to gain an insight into the physical nature of the drug-receptor interaction. Electrostatic properties were derived from both the experimental electron density (multipole refinement of X-ray data collected at T=17 K) and the ab initio wavefunction (single molecule and fully periodic calculations at the DFT level). The relevance of SO and SN intramolecular interactions on the activity of 1 is highlighted by using both the crystal and gas-phase geometries and their electrostatic nature is documented by means of QTAIM atomic charges. The derived electrostatic properties are consistent with a nearly spherical electron density distribution, characterised by an intermingling of electropositive and -negative zones rather than by a unique electrophilic region opposed to a nucleophilic area. This makes the first molecular moment scarcely significant and ill-determined, whereas the second moment is large, significant and highly reliable. A comparison between experimental and theoretical components of the third electric moment shows a few discrepancies, whereas the agreement for the fourth electric moment is excellent. The most favourable intermolecular bond is show to be an NHN hydrogen bond with an energy of about 50 kJ mol(-1). Key pharmacophoric features responsible for attractive electrostatic interactions include CHX hydrogen bonds. It is shown that methyl and methylene groups, known to be essential for the biological activity of the drug, provide a significant energetic contribution to the total binding energy. Dispersive interactions are important at the thiophene and at both the phenyl fragments. The experimental estimates of the electrostatic contribution to the intermolecular interaction energies of six molecular pairs, obtained by a new model proposed by Spackman, predict the correct relative electrostatic energies with no exceptions.
Collapse
Affiliation(s)
- Raffaella Soave
- CNR-ISTM, Istituto di Scienze e Tecnologie Molecolari, Via Golgi 19, 20133 Milano, Italy.
| | | | | |
Collapse
|
46
|
Spackman MA, Munshi P, Dittrich B. Dipole Moment Enhancement in Molecular Crystals from X-ray Diffraction Data. Chemphyschem 2007; 8:2051-63. [PMID: 17676648 DOI: 10.1002/cphc.200700339] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Although reliable determination of the molecular dipole moment from experimental charge density analyses on molecular crystals is a challenging undertaking, these values are becoming increasingly common experimental results. We collate all known experimental determinations and use this database to identify broad trends in the dipole moment enhancements implied by these measurements as well as outliers for which enhancements are pronounced. Compelling evidence emerges that molecular dipole moments from X-ray diffraction data can provide a wealth of information on the change in the molecular charge distribution that results from crystal formation. Most importantly, these experiments are unrivalled in their potential to provide this information in such detail and deserve to be exploited to a much greater extent. The considerable number of experimental determinations now available has enabled us to pinpoint those studies that merit further attention, either because they point unequivocally to a considerable enhancement in the crystal (of 50 % or more), or because the experimental determinations suggest enhancements of 100 % or more--much larger than independent theoretical estimates. In both cases further detailed experimental and theoretical studies are indicated.
Collapse
Affiliation(s)
- Mark A Spackman
- School of Biomedical, Biomolecular & Chemical Sciences, University of Western Australia, Crawley WA 6009, Australia.
| | | | | |
Collapse
|
47
|
Spackman MA, Munshi P, Jayatilaka D. The use of dipole lattice sums to estimate electric fields and dipole moment enhancement in molecular crystals. Chem Phys Lett 2007. [DOI: 10.1016/j.cplett.2007.06.058] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
48
|
Whitten AE, Jayatilaka D, Spackman MA. Effective molecular polarizabilities and crystal refractive indices estimated from x-ray diffraction data. J Chem Phys 2007; 125:174505. [PMID: 17100452 DOI: 10.1063/1.2364897] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Although it was proposed some time ago that (hyper)polarizabilities might be estimated from the results of x-ray charge density refinements, early results were unconvincing. In this work we show that the one particle density obtained from the usual multipole refinement model does not contain sufficient information to determine these response properties and instead pursue the "constrained wave function" approach of fitting to x-ray structure factors. Simplified sum-over-states expressions are derived for determining the dipole polarizability from these wave functions, and these clearly show that the earlier work ignored important two-electron expectation values for the dipole polarizability, and two- and three-electron terms for Beta, etc. Correction factors for the simplified sum-over-states polarizability tensors from the constrained wave function are obtained by calibration against coupled Hartree-Fock ab initio results to yield in-crystal effective polarizability tensors. Results obtained for benzene, urea, and 2-methyl-4-nitroaniline demonstrate that the effective molecular polarizabilities clearly include the effects of intermolecular interactions and electron correlation, especially for urea where the effects on the polarizability are known to be quite large. We also carefully consider the way in which the linear bulk susceptibility, chi((1)), and refractive indices are determined from the x-ray fitted polarizabilities, employing three models based on a rigorous treatment of the local field. Incorrect results are obtained for the sort of molecules that are of interest in nonlinear optical applications if the molecules are approximated by single point dipoles. In contrast, the use of Lorentz-factor tensors averaged over several sites yields excellent results, with refractive indices obtained using this model in remarkably good agreement with optical measurements extrapolated to zero frequency.
Collapse
Affiliation(s)
- Andrew E Whitten
- Chemistry, University of New England, Armidale, New South Wales 2351, Australia
| | | | | |
Collapse
|