1
|
Ligório RF, Rodrigues JL, Krawczuk A, Dos Santos LHR. A building-block database of distributed polarizabilities and dipole moments to estimate optical properties of biomacromolecules in isolation or in an explicitly solvated medium. J Comput Chem 2023; 44:745-754. [PMID: 36433655 DOI: 10.1002/jcc.27037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Revised: 09/10/2022] [Accepted: 10/24/2022] [Indexed: 11/27/2022]
Abstract
Since atomic or functional-group properties in the bulk are generally not available from experimental methods, computational approaches based on partitioning schemes have emerged as a rapid yet accurate pathway to estimate the materials behavior from chemically meaningful building blocks. Among several applications, a comprehensive and systematically built database of atomic or group polarizabilities and related opto-electronic quantities would be very useful not only to envisage linear or non-linear optical properties of biomacromolecules but also to improve the accuracy of classical force fields devoted to simulate biochemical processes. In this work, we propose the first entries of such database that contains distributed polarizabilities and dipole moments extracted from fragments of peptides. Twenty three prototypical conformers of the dipeptides alanine-alanine and glycine-glycine were used to extract functional groups such as CH2 , CHCH3 , NH2 , COOH, CONH, thus allowing construction of a diversity of chemically relevant environments. To evaluate the accuracy of our database, reconstructed properties of larger peptides containing up to six residues of alanine and glycine were tested against density functional theory calculations at the M06-HF/aug-cc-pVDZ level of theory. The procedure is particularly accurate for the diagonal components of the polarizability tensor with errors up to 15%. In order to include solvent effects explicitly, the peptides were also surrounded by a box of water molecules whose distribution was optimized using the CHARMM force field. Solvent effects introduced by a classical dipole-dipole interaction model were compared to those obtained from polarizable-continuum model calculations.
Collapse
Affiliation(s)
- Raphael F Ligório
- Institut für Anorganische Chemie, Universität Göttingen, Göttingen, Germany
| | - José L Rodrigues
- Departamento de Química, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil.,Departamento de Química, Instituto Federal de Educação, Ciência e Tecnologia do Maranhão - Campus Grajaú, Grajaú, MA, Brazil
| | - Anna Krawczuk
- Institut für Anorganische Chemie, Universität Göttingen, Göttingen, Germany
| | | |
Collapse
|
2
|
Time-dependent DFT calculations of the dipole moment and polarizability for excited states. ADVANCES IN QUANTUM CHEMISTRY 2021. [DOI: 10.1016/bs.aiq.2021.05.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
3
|
Ligorio RF, Krawczuk A, Dos Santos LHR. Crystal Field Effects on Atomic and Functional-Group Distributed Polarizabilities of Molecular Materials. J Phys Chem A 2020; 124:10008-10018. [DOI: 10.1021/acs.jpca.0c09293] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Raphael F. Ligorio
- Departamento de Química, Universidade Federal de Minas Gerais, Av. Pres. Antônio Carlos 6627, Belo Horizonte 31270-901, Minas Gerais, Brazil
| | - Anna Krawczuk
- Institut für Anorganische Chemie, Unisersität Göttingen, Tammannstrasse 4, D-37077 Göttingen, Germany
| | - Leonardo H. R. Dos Santos
- Departamento de Química, Universidade Federal de Minas Gerais, Av. Pres. Antônio Carlos 6627, Belo Horizonte 31270-901, Minas Gerais, Brazil
| |
Collapse
|
4
|
|
5
|
Visentin G, Buchachenko AA. Polarizabilities, dispersion coefficients, and retardation functions at the complete basis set CCSD limit: From Be to Ba plus Yb. J Chem Phys 2019; 151:214302. [PMID: 31822071 DOI: 10.1063/1.5129583] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Static and dynamic polarizabilities of alkaline earth metal atoms Be-Ba and of the Yb atom, as well as dispersion coefficients and retardation functions for their long-range interactions, are used as a benchmark for the restricted coupled cluster method with singles and doubles (CCSD) and noniterative triples added [CCSD(T)] and related polarization propagator CCSD(3) methods at the complete basis set limit. The latter is attained through the sequence of the augmented correlation-consistent polarized weighted core valence n-zeta basis sets with the exact 2-component approximation for the scalar relativistic effects and with the small-core effective core potentials (for Ca, Sr, and Ba). At the converged level of core correlation treatment, the finite-field CCSD(T) method reproduces the best available data for the static dipole and quadrupole polarizabilities better than 1% and 4%, respectively. Systematic cancelation of the triple contribution in the CCSD(3) calculations of the dynamic polarizabilities of alkaline earth metal atoms makes their dispersion coefficients accurate within 3%. The retardation functions are computed and used for the analysis of the long-range interactions in the homonuclear dimers. Implications to accurate ab initio calculations of the global interaction potentials are discussed.
Collapse
Affiliation(s)
- Giorgio Visentin
- Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Moscow 121205, Russia
| | - Alexei A Buchachenko
- Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Moscow 121205, Russia
| |
Collapse
|
6
|
Ivanova B, Spiteller M. Noncentrosymmetric organic crystals of barbiturates as potential nonlinear optical phores: experimental and theoretical analyses. CHEMICAL PAPERS 2019. [DOI: 10.1007/s11696-019-00835-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
7
|
Milne AW, Jorge M. Polarization Corrections and the Hydration Free Energy of Water. J Chem Theory Comput 2018; 15:1065-1078. [DOI: 10.1021/acs.jctc.8b01115] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- Andrew W. Milne
- Department of Chemical and Process Engineering, University of Strathclyde, Glasgow G1 1XQ, United Kingdom
| | - Miguel Jorge
- Department of Chemical and Process Engineering, University of Strathclyde, Glasgow G1 1XQ, United Kingdom
| |
Collapse
|
8
|
Lao KU, Jia J, Maitra R, DiStasio RA. On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions. J Chem Phys 2018; 149:204303. [PMID: 30501247 DOI: 10.1063/1.5051458] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we investigate how geometric changes influence the static dipole polarizability (α) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-ζ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this ab initio surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order C 6 dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H2O. Using first-order perturbation theory, we compute the zero-point vibrational correction to α at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic (α iso) and nearly 50% to the anisotropic (α aniso) polarizability values. In doing so, we find that α iso = 9.8307 bohr3, which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr3 provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.
Collapse
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Junteng Jia
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Rahul Maitra
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
9
|
Uhlig F, Zeman J, Smiatek J, Holm C. First-Principles Parametrization of Polarizable Coarse-Grained Force Fields for Ionic Liquids. J Chem Theory Comput 2018; 14:1471-1486. [DOI: 10.1021/acs.jctc.7b00903] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Frank Uhlig
- Institute for Computational Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Johannes Zeman
- Institute for Computational Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Jens Smiatek
- Institute for Computational Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| | - Christian Holm
- Institute for Computational Physics, University of Stuttgart, D-70569 Stuttgart, Germany
| |
Collapse
|
10
|
Elking DM. Torque and atomic forces for Cartesian tensor atomic multipoles with an application to crystal unit cell optimization. J Comput Chem 2016; 37:2067-80. [PMID: 27349179 DOI: 10.1002/jcc.24427] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2016] [Accepted: 05/17/2016] [Indexed: 01/31/2023]
Abstract
New equations for torque and atomic force are derived for use in flexible molecule force fields with atomic multipoles. The expressions are based on Cartesian tensors with arbitrary multipole rank. The standard method for rotating Cartesian tensor multipoles and calculating torque is to first represent the tensor with n indexes and 3(n) redundant components. In this work, new expressions for directly rotating the unique (n + 1)(n + 2)/2 Cartesian tensor multipole components Θpqr are given by introducing Cartesian tensor rotation matrix elements X(R). A polynomial expression and a recursion relation for X(R) are derived. For comparison, the analogous rotation matrix for spherical tensor multipoles are the Wigner functions D(R). The expressions for X(R) are used to derive simple equations for torque and atomic force. The torque and atomic force equations are applied to the geometry optimization of small molecule crystal unit cells. In addition, a discussion of computational efficiency as a function of increasing multipole rank is given for Cartesian tensors. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Dennis M Elking
- Openeye Scientific Software, Santa Fe, New Mexico, 87508.,Elking Scientific Modeling, Ballwin, Missouri, 63102
| |
Collapse
|
11
|
Elking DM, Fusti-Molnar L, Nichols A. Crystal structure prediction of rigid molecules. ACTA CRYSTALLOGRAPHICA SECTION B-STRUCTURAL SCIENCE CRYSTAL ENGINEERING AND MATERIALS 2016; 72:488-501. [DOI: 10.1107/s2052520616010118] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2016] [Accepted: 06/21/2016] [Indexed: 11/11/2022]
Abstract
A non-polarizable force field based on atomic multipoles fit to reproduce experimental crystal properties andab initiogas-phase dimers is described. The Ewald method is used to calculate both long-range electrostatic and 1/r6dispersion energies of crystals. The dispersion energy of a crystal calculated by a cutoff method is shown to converge slowly to the exact Ewald result. A method for constraining space-group symmetry during unit-cell optimization is derived. Results for locally optimizing 4427 unit cells including volume, cell parameters, unit-cell r.m.s.d. and CPU timings are given for both flexible and rigid molecule optimization. An algorithm for randomly generating rigid molecule crystals is described. Using the correct experimentally determined space group, the average and maximum number of random crystals needed to find the correct experimental structure is given for 2440 rigid single component crystals. The force field energy rank of the correct experimental structure is presented for the same set of 2440 rigid single component crystals assuming the correct space group. A complete crystal prediction is performed for two rigid molecules by searching over the 32 most probable space groups.
Collapse
|
12
|
Loboda O, Ingrosso F, Ruiz-López MF, Reis H, Millot C. Dipole and quadrupole polarizabilities of the water molecule as a function of geometry. J Comput Chem 2016; 37:2125-32. [DOI: 10.1002/jcc.24431] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2016] [Revised: 05/24/2016] [Accepted: 05/27/2016] [Indexed: 11/11/2022]
Affiliation(s)
- Oleksandr Loboda
- SRSMC UMR 7565, Faculté des Sciences et Technologies, Boulevard des Aiguillettes, BP 70239; Université de Lorraine, CNRS; Vandoeuvre-lès-Nancy F-54506 France
| | - Francesca Ingrosso
- SRSMC UMR 7565, Faculté des Sciences et Technologies, Boulevard des Aiguillettes, BP 70239; Université de Lorraine, CNRS; Vandoeuvre-lès-Nancy F-54506 France
| | - Manuel F. Ruiz-López
- SRSMC UMR 7565, Faculté des Sciences et Technologies, Boulevard des Aiguillettes, BP 70239; Université de Lorraine, CNRS; Vandoeuvre-lès-Nancy F-54506 France
| | - Heribert Reis
- Institute of Biology, Medicinal Chemistry and Biotechnology, National Hellenic Research Foundation; 48 Vas. Constantinou Avenue Athens 11635 Greece
| | - Claude Millot
- SRSMC UMR 7565, Faculté des Sciences et Technologies, Boulevard des Aiguillettes, BP 70239; Université de Lorraine, CNRS; Vandoeuvre-lès-Nancy F-54506 France
| |
Collapse
|
13
|
Tafipolsky M, Ansorg K. Toward a Physically Motivated Force Field: Hydrogen Bond Directionality from a Symmetry-Adapted Perturbation Theory Perspective. J Chem Theory Comput 2016; 12:1267-79. [DOI: 10.1021/acs.jctc.5b01057] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Maxim Tafipolsky
- Institut
für Physikalische
und Theoretische Chemie, Universität Würzburg, Campus
Hubland Nord, Emil-Fischer-Strasse 42, D-97074 Würzburg, Germany
| | - Kay Ansorg
- Institut
für Physikalische
und Theoretische Chemie, Universität Würzburg, Campus
Hubland Nord, Emil-Fischer-Strasse 42, D-97074 Würzburg, Germany
| |
Collapse
|
14
|
Tao J, Rappe AM. Communication: Accurate higher-order van der Waals coefficients between molecules from a model dynamic multipole polarizability. J Chem Phys 2016; 144:031102. [DOI: 10.1063/1.4940397] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jianmin Tao
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Andrew M. Rappe
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, USA
| |
Collapse
|
15
|
Hatz R, Hänninen V, Halonen L. Pair-potential approach to accurate dispersion energies between group 12 (Zn, Cd, Hg) clusters. J Phys Chem A 2014; 118:12274-9. [PMID: 25427295 DOI: 10.1021/jp510622u] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Dispersion interactions between group 12 (Zn, Cd, Hg) metal clusters are studied at the CCSD(T) level with triple-ζ basis sets. We use atomic orientation dependent C6 dispersion coefficients derived from simple model systems to calculate the intermolecular dispersion energy between larger metal clusters. By using an atomic pair-potential model, we are able to obtain highly accurate results comparable to the CCSD(T) level of theory. For all the systems studied, our method results in a robust description of the dispersion energy at long distances, and it is able to outperform the examined dispersion corrected DFT calculations.
Collapse
Affiliation(s)
- Richard Hatz
- Laboratory of Physical Chemistry, Department of Chemistry, University of Helsinki , P.O. Box 55 (A.I. Virtasen aukio 1), FIN-00014 Helsinki, Finland
| | | | | |
Collapse
|
16
|
Hatz R, Hänninen V, Halonen L. Dispersion Interactions in Small Zinc, Cadmium, and Mercury Clusters. J Phys Chem A 2014; 118:5734-40. [DOI: 10.1021/jp505023g] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Richard Hatz
- Laboratory
of Physical Chemistry, Department of Chemistry, University of Helsinki, P.O. Box 55
(A.I. Virtasen aukio 1), FIN-00014 Helsinki, Finland
| | - Vesa Hänninen
- Laboratory
of Physical Chemistry, Department of Chemistry, University of Helsinki, P.O. Box 55
(A.I. Virtasen aukio 1), FIN-00014 Helsinki, Finland
| | - Lauri Halonen
- Laboratory
of Physical Chemistry, Department of Chemistry, University of Helsinki, P.O. Box 55
(A.I. Virtasen aukio 1), FIN-00014 Helsinki, Finland
| |
Collapse
|
17
|
Limiting assumptions in molecular modeling: electrostatics. J Comput Aided Mol Des 2013; 27:107-14. [PMID: 23354627 DOI: 10.1007/s10822-013-9634-x] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2012] [Accepted: 01/15/2013] [Indexed: 10/27/2022]
Abstract
Molecular mechanics attempts to represent intermolecular interactions in terms of classical physics. Initial efforts assumed a point charge located at the atom center and coulombic interactions. It is been recognized over multiple decades that simply representing electrostatics with a charge on each atom failed to reproduce the electrostatic potential surrounding a molecule as estimated by quantum mechanics. Molecular orbitals are not spherically symmetrical, an implicit assumption of monopole electrostatics. This perspective reviews recent evidence that requires use of multipole electrostatics and polarizability in molecular modeling.
Collapse
|
18
|
May F, Baumeier B, Lennartz C, Andrienko D. Can lattice models predict the density of states of amorphous organic semiconductors? PHYSICAL REVIEW LETTERS 2012; 109:136401. [PMID: 23030109 DOI: 10.1103/physrevlett.109.136401] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/03/2012] [Indexed: 06/01/2023]
Abstract
We extend existing lattice models of small-molecule amorphous semiconductors by accounting for changes in molecular polarizability upon charging or excitation. A compact expression of this contribution to the density of states is provided. Although the lattice model and the description based on a microscopic morphology both qualitatively predict an additional broadening, shift, and an exponential tail (traps) of the density of states, a quantitative agreement between the two cannot be achieved.
Collapse
Affiliation(s)
- Falk May
- Max Planck Institute for Polymer Research, Mainz, Germany
| | | | | | | |
Collapse
|
19
|
Elking DM, Perera L, Pedersen LG. HPAM: Hirshfeld Partitioned Atomic Multipoles. COMPUTER PHYSICS COMMUNICATIONS 2012; 183:390-397. [PMID: 22140274 PMCID: PMC3225920 DOI: 10.1016/j.cpc.2011.10.003] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
An implementation of the Hirshfeld (HD) and Hirshfeld-Iterated (HD-I) atomic charge density partitioning schemes is described. Atomic charges and atomic multipoles are calculated from the HD and HD-I atomic charge densities for arbitrary atomic multipole rank l(max) on molecules of arbitrary shape and size. The HD and HD-I atomic charges/multipoles are tested by comparing molecular multipole moments and the electrostatic potential (ESP) surrounding a molecule with their reference ab initio values. In general, the HD-I atomic charges/multipoles are found to better reproduce ab initio electrostatic properties over HD atomic charges/multipoles. A systematic increase in precision for reproducing ab initio electrostatic properties is demonstrated by increasing the atomic multipole rank from l(max) = 0 (atomic charges) to l(max) = 4 (atomic hexadecapoles). Both HD and HD-I atomic multipoles up to rank l(max) are shown to exactly reproduce ab initio molecular multipole moments of rank L for L ≤ l(max). In addition, molecular dipole moments calculated by HD, HD-I, and ChelpG atomic charges only (l(max) = 0) are compared with reference ab initio values. Significant errors in reproducing ab initio molecular dipole moments are found if only HD or HD-I atomic charges used.
Collapse
Affiliation(s)
- Dennis M. Elking
- University of North Carolina, Department of Chemistry, Chapel Hill, NC 27599, USA
| | - Lalith Perera
- Laboratory of Structural Biology, National Institute of Environmental Health Sciences, Research Triangle Park, North Carolina 27709, USA
| | - Lee G. Pedersen
- University of North Carolina, Department of Chemistry, Chapel Hill, NC 27599, USA
- Laboratory of Structural Biology, National Institute of Environmental Health Sciences, Research Triangle Park, North Carolina 27709, USA
| |
Collapse
|