1
|
Hernández-Segura LI, Olvera-Rubalcava FA, Flores-Moreno R, Calaminici P, Köster AM. Exchange-correlation kernel for perturbation dependent auxiliary functions in auxiliary density perturbation theory. J Mol Model 2024; 30:302. [PMID: 39115689 PMCID: PMC11310252 DOI: 10.1007/s00894-024-06091-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2024] [Accepted: 07/23/2024] [Indexed: 08/11/2024]
Abstract
CONTEXT Analytic exchange-correlation kernel formulations are of the outermost importance for density functional theory (DFT) perturbation calculations. In this paper, the working equation for the exchange-correlation kernel of the generalized gradient approximation (GGA) for perturbation dependent auxiliary functions is derived and discussed in the framework of auxiliary density functional theory (ADFT). The presented new formulation is extended to the unrestricted approach, too. A comprehensive discussion of the implementation of the GGA ADFT kernel, using either the native exchange-correlation functional implementations in deMon2k or the ones from the LibXC library, is given. Calculations with analytic exchange-correlation kernels are compared to their finite difference counterparts. The obtained results are in quantitative agreement. Nevertheless, analytic GGA ADFT kernel implementations show substantial improvement in the computational performance. Similar results are reported for analytic second derivatives of effective core potential (ECP) and model core potential (MCP) matrix elements when compared to their finite difference counterparts in molecular frequency analyses. METHOD All calculations are performed in the framework of ADFT as implemented in deMon2k. In the ADFT analytic frequency calculations, auxiliary density perturbation theory was used. The underlying two-center exchange-correlation kernel matrix elements are calculated by numerical integration either with analytic or finite difference kernel expressions. Validation calculations are performed with the VWN and PBE functionals employing DFT-optimized DZVP basis sets in conjunction with automatically generated GEN-A2 auxiliary density function sets. In the (Pt3Cu)n cluster benchmark calculations, the RPBE functional was used. For Pt atoms, the quasi-relativistic LANL2DZ effective core potential with the corresponding valence basis set was employed, whereas for Cu atoms, the all-electron DFT-optimized TZVP basis was applied. The auxiliary density was expanded by the automatically generated GEN-A2* auxiliary function set. We run all benchmark calculations in parallel on 24 cores.
Collapse
Affiliation(s)
- Luis I Hernández-Segura
- Chemistry Department, CINVESTAV, Av. Instituto Politecnico Nacional 2508, Col. San Pedro Zacatenco, Del. Gustavo A. Madero, Mexico City, C.P. 07360, Mexico.
| | - Flor A Olvera-Rubalcava
- Chemistry Department, CINVESTAV, Av. Instituto Politecnico Nacional 2508, Col. San Pedro Zacatenco, Del. Gustavo A. Madero, Mexico City, C.P. 07360, Mexico
| | - Roberto Flores-Moreno
- Departamento de Química, Universidad de Guadalajara, Blvd. Gral. Marcelino García Barragán 1421, Guadalajara, Jal., C.P. 44430, Mexico
| | - Patrizia Calaminici
- Chemistry Department, CINVESTAV, Av. Instituto Politecnico Nacional 2508, Col. San Pedro Zacatenco, Del. Gustavo A. Madero, Mexico City, C.P. 07360, Mexico
| | - Andreas M Köster
- Chemistry Department, CINVESTAV, Av. Instituto Politecnico Nacional 2508, Col. San Pedro Zacatenco, Del. Gustavo A. Madero, Mexico City, C.P. 07360, Mexico.
| |
Collapse
|
2
|
Luna-Valenzuela A, Pedroza-Montero JN, Köster AM, Calaminici P, Gálvez-González LE, Posada-Amarillas A. Pd 8 Cluster: Too Small to Melt? A BOMD Study. J Phys Chem A 2024; 128:572-580. [PMID: 38207112 DOI: 10.1021/acs.jpca.3c06173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
The question of whether a solid-liquid phase transition occurs in small clusters poses a fundamental challenge. In this study, we attempt to elucidate this phenomenon through a thorough examination of the thermal behavior and structural stability of Pd8 clusters employing ab initio simulations. Initially, a systematic global search is carried out to identify the various isomers of the Pd8 cluster. This is accomplished by employing an ab initio basin-hopping algorithm and using the PBE/SDD scheme integrated in the Gaussian code. The resulting isomers are further refined through reoptimization using the deMon2k package. To ensure the structural firmness of the lowest-energy isomer, we calculated normal modes. The structural stability as a function of temperature is analyzed through the Born-Oppenheimer molecular dynamics (BOMD) approach. Multiple BOMD trajectories at distinct simulated temperatures are examined with data clustering analysis to determine cluster isomers. This analysis establishes a connection between the potential energy landscape and the simulated temperature. To address the question of cluster melting, canonical parallel-tempering BOMD runs are performed and analyzed with the multiple-histogram method. A broad maximum in the heat capacity curve indicates a melting transition between 500 and 600 K. To further examine this transition, the mean-squared displacement and the pair-distance distribution function are calculated. The results of these calculations confirm the existence of a solid-liquid phase transition, as indicated by the heat capacity curve.
Collapse
Affiliation(s)
- Analila Luna-Valenzuela
- Departamento de Ciencias de la Salud, Unidad Regional Los Mochis, Universidad Autónoma de Occidente, Blvd. Macario Gaxiola S/N, 81217 Los Mochis, Sinaloa, México
| | - Jesús N Pedroza-Montero
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional, 2508 Ciudad de Mexico, México
| | - Andreas M Köster
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional, 2508 Ciudad de Mexico, México
| | - Patrizia Calaminici
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional, 2508 Ciudad de Mexico, México
| | - Luis E Gálvez-González
- Departamento de Investigación en Física, Universidad de Sonora, Blvd. Luis Encinas & Rosales, 83000 Hermosillo, Sonora, México
| | - Alvaro Posada-Amarillas
- Departamento de Investigación en Física, Universidad de Sonora, Blvd. Luis Encinas & Rosales, 83000 Hermosillo, Sonora, México
| |
Collapse
|
3
|
Diop M, El-Hayek M, Attard J, Muhieddine A, Veremeienko V, Soorkia S, Carbonnière P, de la Lande A, Soep B, Shafizadeh N. Chlorophyll and pheophytin protonated and deprotonated ions: Observation and theory. J Chem Phys 2023; 159:194308. [PMID: 37987519 DOI: 10.1063/5.0174351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Accepted: 10/30/2023] [Indexed: 11/22/2023] Open
Abstract
Pheophytin a and chlorophyll a have been investigated by electrospray mass spectrometry in the positive and negative modes, in view of the importance of the knowledge of their properties in photosynthesis. Pheophytin and chlorophyll are both observed intensely in the protonated mode, and their main fragmentation route is the loss of their phytyl chain. Pheophytin is observed intact in the negative mode, while under collisions, it is primarily cleaved beyond the phytyl chain and loses the attaching propionate group. Chlorophyll is not detected in normal conditions in the negative mode, but addition of methanol solvent molecule is detected. Fragmentation of this adduct primarily forms a product (-30 amu) that dissociates into dephytyllated deprotonated chlorophyll. Semi-empirical molecular dynamics calculations show that the phytyl chain is unfolded from the chlorin cycle in pheophytin a and folded in chlorophyll a. Density functional theory calculations have been conducted to locate the charges on protonated and deprotonated pheophytin a and chlorophyll a and have found the major location sites that are notably more stable in energy by more than 0.5 eV than the others. The deprotonation site is found identical for pheophytin a and the chlorophyll a-methanol adduct. This is in line with experiment and calculation locating the addition of methanol on a double bond of deprotonated chlorophyll a.
Collapse
Affiliation(s)
- M Diop
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d' Orsay, 91405 Orsay, France
| | - M El-Hayek
- Université Paris Saclay, CNRS, Institut de Chimie Physique, UMR 8000, 91405 Orsay, France
| | - J Attard
- Université Paris Saclay, CNRS, Institut de Chimie Physique, UMR 8000, 91405 Orsay, France
- Institut des sciences analytiques et de physico-chimie pour l'environnement et les materiaux (IPREM), Pau, France
| | - A Muhieddine
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d' Orsay, 91405 Orsay, France
| | - V Veremeienko
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), 91198 Gif-sur-Yvette, France
| | - S Soorkia
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d' Orsay, 91405 Orsay, France
| | - Ph Carbonnière
- Institut des sciences analytiques et de physico-chimie pour l'environnement et les materiaux (IPREM), Pau, France
| | - A de la Lande
- Université Paris Saclay, CNRS, Institut de Chimie Physique, UMR 8000, 91405 Orsay, France
| | - B Soep
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d' Orsay, 91405 Orsay, France
| | - N Shafizadeh
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d' Orsay, 91405 Orsay, France
| |
Collapse
|
4
|
Mejia-Rodriguez D, Aprà E, Autschbach J, Bauman NP, Bylaska EJ, Govind N, Hammond JR, Kowalski K, Kunitsa A, Panyala A, Peng B, Rehr JJ, Song H, Tretiak S, Valiev M, Vila FD. NWChem: Recent and Ongoing Developments. J Chem Theory Comput 2023; 19:7077-7096. [PMID: 37458314 DOI: 10.1021/acs.jctc.3c00421] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2023]
Abstract
This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.
Collapse
Affiliation(s)
- Daniel Mejia-Rodriguez
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Edoardo Aprà
- Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Nicholas P Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Eric J Bylaska
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jeff R Hammond
- Accelerated Computing, NVIDIA Helsinki Oy, Porkkalankatu 1, 00180 Helsinki, Finland
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Alexander Kunitsa
- Zapata Computing, Inc., 100 Federal Street, Boston, Massachusetts 02110, United States
| | - Ajay Panyala
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - John J Rehr
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Marat Valiev
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Fernando D Vila
- Department of Physics, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
5
|
Liang J, Feng X, Liu X, Head-Gordon M. Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments. J Chem Phys 2023; 158:204109. [PMID: 37218699 PMCID: PMC10208678 DOI: 10.1063/5.0152838] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2023] [Accepted: 04/27/2023] [Indexed: 05/24/2023] Open
Abstract
VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. This study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).
Collapse
Affiliation(s)
- Jiashu Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, USA
| | | | - Xiao Liu
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, USA
| | | |
Collapse
|
6
|
Galindo-Uribe CD, Calaminici P, Solorza-Feria O. First-principle investigation of structures and energy properties of $$({\textbf {Pt}}_{3}{} {\textbf {Cu}})_{{\varvec{n}}}$$, n = 10–11 nanoclusters. Theor Chem Acc 2023. [DOI: 10.1007/s00214-023-02963-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
|
7
|
Teale AM, Helgaker T, Savin A, Adamo C, Aradi B, Arbuznikov AV, Ayers PW, Baerends EJ, Barone V, Calaminici P, Cancès E, Carter EA, Chattaraj PK, Chermette H, Ciofini I, Crawford TD, De Proft F, Dobson JF, Draxl C, Frauenheim T, Fromager E, Fuentealba P, Gagliardi L, Galli G, Gao J, Geerlings P, Gidopoulos N, Gill PMW, Gori-Giorgi P, Görling A, Gould T, Grimme S, Gritsenko O, Jensen HJA, Johnson ER, Jones RO, Kaupp M, Köster AM, Kronik L, Krylov AI, Kvaal S, Laestadius A, Levy M, Lewin M, Liu S, Loos PF, Maitra NT, Neese F, Perdew JP, Pernal K, Pernot P, Piecuch P, Rebolini E, Reining L, Romaniello P, Ruzsinszky A, Salahub DR, Scheffler M, Schwerdtfeger P, Staroverov VN, Sun J, Tellgren E, Tozer DJ, Trickey SB, Ullrich CA, Vela A, Vignale G, Wesolowski TA, Xu X, Yang W. DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science. Phys Chem Chem Phys 2022; 24:28700-28781. [PMID: 36269074 PMCID: PMC9728646 DOI: 10.1039/d2cp02827a] [Citation(s) in RCA: 62] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 08/09/2022] [Indexed: 12/13/2022]
Abstract
In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.
Collapse
Affiliation(s)
- Andrew M. Teale
- School of Chemistry, University of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andreas Savin
- Laboratoire de Chimie Théorique, CNRS and Sorbonne University, 4 Place Jussieu, CEDEX 05, 75252 Paris, France.
| | - Carlo Adamo
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany.
| | - Alexei V. Arbuznikov
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7Straße des 17. Juni 13510623Berlin
| | | | - Evert Jan Baerends
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56125 Pisa, Italy.
| | - Patrizia Calaminici
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Eric Cancès
- CERMICS, Ecole des Ponts and Inria Paris, 6 Avenue Blaise Pascal, 77455 Marne-la-Vallée, France.
| | - Emily A. Carter
- Department of Mechanical and Aerospace Engineering and the Andlinger Center for Energy and the Environment, Princeton UniversityPrincetonNJ 08544-5263USA
| | | | - Henry Chermette
- Institut Sciences Analytiques, Université Claude Bernard Lyon1, CNRS UMR 5280, 69622 Villeurbanne, France.
| | - Ilaria Ciofini
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - T. Daniel Crawford
- Department of Chemistry, Virginia TechBlacksburgVA 24061USA,Molecular Sciences Software InstituteBlacksburgVA 24060USA
| | - Frank De Proft
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | | | - Claudia Draxl
- Institut für Physik and IRIS Adlershof, Humboldt-Universität zu Berlin, 12489 Berlin, Germany. .,Fritz-Haber-Institut der Max-Planck-Gesellschaft, 14195 Berlin, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany. .,Beijing Computational Science Research Center (CSRC), 100193 Beijing, China.,Shenzhen JL Computational Science and Applied Research Institute, 518110 Shenzhen, China
| | - Emmanuel Fromager
- Laboratoire de Chimie Quantique, Institut de Chimie, CNRS/Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France.
| | - Patricio Fuentealba
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile.
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute, and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, USA.
| | - Giulia Galli
- Pritzker School of Molecular Engineering and Department of Chemistry, The University of Chicago, Chicago, IL, USA.
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China. .,Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA
| | - Paul Geerlings
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | - Nikitas Gidopoulos
- Department of Physics, Durham University, South Road, Durham DH1 3LE, UK.
| | - Peter M. W. Gill
- School of Chemistry, University of SydneyCamperdown NSW 2006Australia
| | - Paola Gori-Giorgi
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Andreas Görling
- Chair of Theoretical Chemistry, University of Erlangen-Nuremberg, Egerlandstrasse 3, 91058 Erlangen, Germany.
| | - Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University, Gold Coast, Qld 4222, Australia.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany.
| | - Oleg Gritsenko
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Hans Jørgen Aagaard Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark.
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie UniversityHalifaxNova ScotiaB3H 4R2Canada
| | - Robert O. Jones
- Peter Grünberg Institut PGI-1, Forschungszentrum Jülich52425 JülichGermany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, 10623, Berlin.
| | - Andreas M. Köster
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav)CDMX07360Mexico
| | - Leeor Kronik
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth, 76100, Israel.
| | - Anna I. Krylov
- Department of Chemistry, University of Southern CaliforniaLos AngelesCalifornia 90089USA
| | - Simen Kvaal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andre Laestadius
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Mel Levy
- Department of Chemistry, Tulane University, New Orleans, Louisiana, 70118, USA.
| | - Mathieu Lewin
- CNRS & CEREMADE, Université Paris-Dauphine, PSL Research University, Place de Lattre de Tassigny, 75016 Paris, France.
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC 27599-3420, USA. .,Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599-3290, USA
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France.
| | - Neepa T. Maitra
- Department of Physics, Rutgers University at Newark101 Warren StreetNewarkNJ 07102USA
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany.
| | - John P. Perdew
- Departments of Physics and Chemistry, Temple UniversityPhiladelphiaPA 19122USA
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.
| | - Pascal Pernot
- Institut de Chimie Physique, UMR8000, CNRS and Université Paris-Saclay, Bât. 349, Campus d'Orsay, 91405 Orsay, France.
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA. .,Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA
| | - Elisa Rebolini
- Institut Laue Langevin, 71 avenue des Martyrs, 38000 Grenoble, France.
| | - Lucia Reining
- Laboratoire des Solides Irradiés, CNRS, CEA/DRF/IRAMIS, École Polytechnique, Institut Polytechnique de Paris, F-91120 Palaiseau, France. .,European Theoretical Spectroscopy Facility
| | - Pina Romaniello
- Laboratoire de Physique Théorique (UMR 5152), Université de Toulouse, CNRS, UPS, France.
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA.
| | - Dennis R. Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS – Centre for Molecular Simulation, IQST – Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary2500 University Drive NWCalgaryAlbertaT2N 1N4Canada
| | - Matthias Scheffler
- The NOMAD Laboratory at the FHI of the Max-Planck-Gesellschaft and IRIS-Adlershof of the Humboldt-Universität zu Berlin, Faradayweg 4-6, D-14195, Germany.
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632 Auckland, New Zealand.
| | - Viktor N. Staroverov
- Department of Chemistry, The University of Western OntarioLondonOntario N6A 5B7Canada
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA.
| | - Erik Tellgren
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - David J. Tozer
- Department of Chemistry, Durham UniversitySouth RoadDurhamDH1 3LEUK
| | - Samuel B. Trickey
- Quantum Theory Project, Deptartment of Physics, University of FloridaGainesvilleFL 32611USA
| | - Carsten A. Ullrich
- Department of Physics and Astronomy, University of MissouriColumbiaMO 65211USA
| | - Alberto Vela
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Giovanni Vignale
- Department of Physics, University of Missouri, Columbia, MO 65203, USA.
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, Université de Genève30 Quai Ernest-Ansermet1211 GenèveSwitzerland
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovation Materials, Collaborative Innovation Centre of Chemistry for Energy Materials, MOE Laboratory for Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, China.
| | - Weitao Yang
- Department of Chemistry and Physics, Duke University, Durham, NC 27516, USA.
| |
Collapse
|
8
|
Pedroza-Montero JN, Calaminici P, Köster AM. First-principle polarizabilities of nanosystems from auxiliary density perturbation theory with MINRES. Theor Chem Acc 2022. [DOI: 10.1007/s00214-021-02864-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
9
|
Tandiana R, Sicard-Roselli C, Van-Oanh NT, Steinmann S, Clavaguéra C. In-depth theoretical understanding of the chemical interaction of aromatic compounds with a gold nanoparticle. Phys Chem Chem Phys 2022; 24:25327-25336. [DOI: 10.1039/d2cp02654f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
The orientations of aromatic molecules at the surface of gold nanoparticles are probed and characterized by a combination of several topological analyses, energy decomposition analyses, and infrared spectroscopy.
Collapse
Affiliation(s)
- Rika Tandiana
- Institut de Chimie Physique, Université Paris-Saclay – CNRS, UMR 8000, 91405 Orsay, France
| | - Cécile Sicard-Roselli
- Institut de Chimie Physique, Université Paris-Saclay – CNRS, UMR 8000, 91405 Orsay, France
| | - Nguyen-Thi Van-Oanh
- Institut de Chimie Physique, Université Paris-Saclay – CNRS, UMR 8000, 91405 Orsay, France
| | | | - Carine Clavaguéra
- Institut de Chimie Physique, Université Paris-Saclay – CNRS, UMR 8000, 91405 Orsay, France
| |
Collapse
|
10
|
Flores-Moreno R, Cortes-Llamas SA, Pineda-Urbina K, Medel VM, Jayaprakash GK. Analytic Alchemical Derivatives for the Analysis of Differential Acidity Assisted by the h Function. J Phys Chem A 2021; 125:10463-10474. [PMID: 34812636 DOI: 10.1021/acs.jpca.1c07364] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Analytical calculation of alchemical derivatives based on auxiliary density perturbation theory is described, coded, and validated. For the case where the nucleus is a hydrogen atom and the nuclear charge is changed from 1 to 0, it turns out that a good estimate of the proton binding energies can be obtained very efficiently. First-order results correspond exactly to the molecular electrostatic potential evaluated at the hydrogen nucleus location (removing self-repulsion), in agreement with previously reported extensive studies. Therefore, the second-order results reported here are refinements in accuracy that finally allow a quantitative exploration of differential acidity. Furthermore, the recently reported h function is produced in its analytical form as a byproduct and local descriptor associated with the proton binding energy values found with this approach. In an example application, proton binding energies are computed for a family of imidazolium derivatives to demonstrate the capabilities and the stability of the method with respect to changes in basis set or exchange-correlation functional.
Collapse
Affiliation(s)
- Roberto Flores-Moreno
- Departamento de Química, Universidad de Guadalajara, Blvd. Marcelino García Barragán 1421, Col. Olímpica, 44430 Guadalajara, Jalisco, México
| | - Sara A Cortes-Llamas
- Departamento de Química, Universidad de Guadalajara, Blvd. Marcelino García Barragán 1421, Col. Olímpica, 44430 Guadalajara, Jalisco, México
| | - Kayim Pineda-Urbina
- Facultad de Ciencias Químicas, Universidad de Colima, Carretera Colima-Coquimatlan Km. 9, 28400 Coquimatlan, Colima, México
| | - Victor M Medel
- Departamento de Química, Universidad de Guadalajara, Blvd. Marcelino García Barragán 1421, Col. Olímpica, 44430 Guadalajara, Jalisco, México
| | | |
Collapse
|
11
|
Geudtner G. Parallelization of deMon2k: an overview. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02786-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
12
|
Galindo-Uribe CD, Calaminici P, Cruz-Martínez H, Cruz-Olvera D, Solorza-Feria O. First-principle study of the structures, growth pattern, and properties of (Pt 3Cu) n, n = 1-9, clusters. J Chem Phys 2021; 154:154302. [PMID: 33887945 DOI: 10.1063/5.0045203] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this work, a first-principles systematic study of (Pt3Cu)n, n = 1-9, clusters was performed employing the linear combination of Gaussian-type orbital auxiliary density functional theory approach. The growth of the clusters has been achieved by increasing the previous cluster by one Pt3Cu unit at a time. To explore in detail the potential energy surface of these clusters, initial structures were obtained from Born-Oppenheimer molecular dynamics trajectories generated at different temperatures and spin multiplicities. For each cluster size, several dozens of structures were optimized without any constraints. The most stable structures were characterized by frequency analysis calculations. This study demonstrates that the obtained most stable structures prefer low spin multiplicities. To gain insight into the growing pattern of these systems, average bond lengths were calculated for the lowest stable structures. This work reveals that the Cu atoms prefer to be together and to localize inside the cluster structures. Moreover, these systems tend to form octahedra moieties in the size range of n going from 4 to 9 Pt3Cu units. Magnetic moment per atom and spin density plots were obtained for the neutral, cationic, and anionic ground state structures. Dissociation energies, ionization potential, and electron affinity were calculated, too. The dissociation energy and the electron affinity increase as the number of Pt3Cu units grows, whereas the ionization potential decreases.
Collapse
Affiliation(s)
- Carlos Daniel Galindo-Uribe
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional 2508, San Pedro Zacatenco, Gustavo A. Madero, CP 07360 Mexico City, Mexico
| | - Patrizia Calaminici
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional 2508, San Pedro Zacatenco, Gustavo A. Madero, CP 07360 Mexico City, Mexico
| | - Heriberto Cruz-Martínez
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional 2508, San Pedro Zacatenco, Gustavo A. Madero, CP 07360 Mexico City, Mexico
| | - Domingo Cruz-Olvera
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional 2508, San Pedro Zacatenco, Gustavo A. Madero, CP 07360 Mexico City, Mexico
| | - Omar Solorza-Feria
- Departamento de Química, CINVESTAV, Av. Instituto Politécnico Nacional 2508, San Pedro Zacatenco, Gustavo A. Madero, CP 07360 Mexico City, Mexico
| |
Collapse
|
13
|
Becker M, Sierka M. Density functional theory for molecular and periodic systems using density fitting and continuous fast multipole method: Stress tensor. J Comput Chem 2019; 40:2563-2570. [PMID: 31322769 DOI: 10.1002/jcc.26033] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2019] [Revised: 06/24/2019] [Accepted: 06/27/2019] [Indexed: 11/06/2022]
Abstract
A full implementation of the analytical stress tensor for periodic systems is reported in the TURBOMOLE program package within the framework of Kohn-Sham density functional theory using Gaussian-type orbitals as basis functions. It is the extension of the implementation of analytical energy gradients (Lazarski et al., Journal of Computational Chemistry 2016, 37, 2518-2526) to the stress tensor for the purpose of optimization of lattice vectors. Its key component is the efficient calculation of the Coulomb contribution by combining density fitting approximation and continuous fast multipole method. For the exchange-correlation (XC) part the hierarchical numerical integration scheme (Burow and Sierka, Journal of Chemical Theory and Computation 2011, 7, 3097-3104) is extended to XC weight derivatives and stress tensor. The computational efficiency and favorable scaling behavior of the stress tensor implementation are demonstrated for various model systems. The overall computational effort for energy gradient and stress tensor for the largest systems investigated is shown to be at most two and a half times the computational effort for the Kohn-Sham matrix formation. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Martin Becker
- Otto-Schott-Institut für Materialforschung, Friedrich-Schiller-Universität Jena, Löbdergraben 32, Jena, D-07743, Germany
| | - Marek Sierka
- Otto-Schott-Institut für Materialforschung, Friedrich-Schiller-Universität Jena, Löbdergraben 32, Jena, D-07743, Germany
| |
Collapse
|
14
|
Pérez-Figueroa SE, Calaminici P, Köster AM. Hybrid ADFT Study of the C 104 and C 106 IPR Isomers. J Phys Chem A 2019; 123:4565-4574. [PMID: 31021089 DOI: 10.1021/acs.jpca.9b00665] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
This work presents a hybrid auxiliary density functional theory (ADFT) study of the neutral and hexaanionic C104 and C106 fullerenes with the aim to determine their ground state structures. To this end, all C104 and C106 fullerene structures that obey the isolated pentagon rule (IPR) were optimized with the Perdew-Burke-Ernzerhof generalized gradient approximation followed by a single-point energy calculation with the PBE0 hybrid functional. Our studies show that this composite approach yields relative energies of giant fullerenes that are accurate to around 1 kcal/mol. As a result, the ground states of C104, C1046-, and C1066- can be assigned to the isomers 234:Cs, 821:D2, and 891:Cs, respectively. On the other hand, the energetically lowest lying IPR isomers of C106, 331:Cs, 1194:C2, 534:C1 are separated by less than 1 kcal/mol which makes an unequivocal ground state assignment by hybrid DFT methods impossible. To guide future experiments, we also report the simulated IR and Raman spectra of the most stable neutral and hexaanionic C104 and C106 fullerenes.
Collapse
|
15
|
Delgado-Venegas RI, Calaminici P, Köster AM. Mixed second and third energy derivatives from auxiliary density perturbation theory. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1549339] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
| | - Patrizia Calaminici
- Departamento de Química, CINVESTAV. Avenida Instituto Politécnico Nacional 2508, México D.F., México
| | - Andreas M. Köster
- Departamento de Química, CINVESTAV. Avenida Instituto Politécnico Nacional 2508, México D.F., México
| |
Collapse
|
16
|
Chou CP, Witek H, Irle S. When finite becomes infinite: convergence properties of vibrational spectra of oligomer chains. J Mol Model 2018; 24:288. [PMID: 30242484 DOI: 10.1007/s00894-018-3824-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Accepted: 09/11/2018] [Indexed: 11/29/2022]
Abstract
We present a computational study of convergence properties of vibrational IR and Raman spectra for a series of increasingly long units of polyethylene, cis- and trans-polyacetylenes, and polyynes. Convergent behavior to the spectra of infinitely long polymers was observed in all cases when chains reached lengths of approximately 60 carbon atoms, both with respect to the positions of the bands and to their intensities. The vibrational spectra of longer chains are practically indistinguishable. The convergence rate depends on the degree of the π conjugation in a studied system: Vibrational spectra for oligoethylenes converge noticeably faster than the spectra for the conjugated systems. The slowest convergence is observed for skeletal motions of the oligomer chains, which may require more than a hundred carbon atoms in the chain to show deviations smaller than 1 cm-1 to the corresponding solid-state calculations. The results suggest that the boundary between the properties of finite and infinite molecular systems fades away for a surprisingly small number of atoms.
Collapse
Affiliation(s)
- Chien-Pin Chou
- Institute of Molecular Science and Department of Applied Chemistry, National Chiao Tung University, Hsinchu, 30010, Taiwan.,Waseda Research Institute for Science and Engineering (WISE), Waseda University, Tokyo, 169-8555, Japan
| | - Henryk Witek
- Institute of Molecular Science and Department of Applied Chemistry, National Chiao Tung University, Hsinchu, 30010, Taiwan. .,Center for Emergent Functional Matter Science, National Chiao Tung University, Hsinchu, 30010, Taiwan.
| | - Stephan Irle
- Institute for Advanced Research and Department of Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, 464-8602, Japan. .,Computational Sciences and Engineering Division & Chemical Sciences Division, ORNL, Oak Ridge, TN, USA.
| |
Collapse
|
17
|
Gómez-Pérez JR, Delesma FA, Calaminici P, Köster AM. Accuracy of auxiliary density functional theory hybrid calculations for activation and reaction enthalpies of pericyclic reactions. J Mol Model 2018; 24:223. [PMID: 30078124 DOI: 10.1007/s00894-018-3759-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2018] [Accepted: 07/24/2018] [Indexed: 11/28/2022]
Abstract
Auxiliary density functional theory (ADFT) hybrid calculations are based on the variational fitting of the Coulomb and Fock potential and, therefore, are free of four-center electron repulsion integrals. So far, ADFT hybrid calculations have been validated successfully for standard enthalpies of formation. In this work the accuracy of ADFT hybrid calculations for the description of pericyclic reactions was quantitatively validated at the B3LYP/6-31G*/GEN-A2* level of theory. Our comparison with conventional Kohn-Sham density functional theory (DFT) results shows that the DFT and ADFT activation and reaction enthalpies are practically indistinguishable. A systematic study of various functionals (PBE, B3LYP, PBE0, CAMB3LYP, CAMPBE0 and HSE06) and basis sets (6-31G*, DZVP-GGA and aug-cc-pVXZ; X = D, T and Q) revealed that the ADFT HSE06/aug-cc-pVTZ/GEN-A2* level of theory yields best balanced accuracy for the activation and reaction enthalpies of the studied pericyclic reactions. With the successfully validate ADFT composite approach consisting of PBE/DZVP-GGA/GEN-A2* structure and transition state optimizations and single-point HSE06/aug-cc-pVTZ/GEN-A2* energy calculations, an accurate, reliable and efficient computational approach for the study of pericyclic reactions in systems at the nanometer scale is proposed.
Collapse
Affiliation(s)
- José R Gómez-Pérez
- Departamento de Química, CINVESTAV, Centro de Investigación y de Estudios Avanzados, Av. Instituto Politécnico Nacional 2508, A.P. 14-740, 07000, México D.F, Mexico
| | - Francisco A Delesma
- Programa de Doctorado en Nanociencias y Nanotecnología, CINVESTAV, Instituto Politécnico Nacional 2508, A.P. 14-740, 07000, Ciudad de México, Mexico
| | - Patrizia Calaminici
- Departamento de Química, CINVESTAV, Centro de Investigación y de Estudios Avanzados, Av. Instituto Politécnico Nacional 2508, A.P. 14-740, 07000, México D.F, Mexico. .,Programa de Doctorado en Nanociencias y Nanotecnología, CINVESTAV, Instituto Politécnico Nacional 2508, A.P. 14-740, 07000, Ciudad de México, Mexico.
| | - Andreas M Köster
- Departamento de Química, CINVESTAV, Centro de Investigación y de Estudios Avanzados, Av. Instituto Politécnico Nacional 2508, A.P. 14-740, 07000, México D.F, Mexico. .,Programa de Doctorado en Nanociencias y Nanotecnología, CINVESTAV, Instituto Politécnico Nacional 2508, A.P. 14-740, 07000, Ciudad de México, Mexico.
| |
Collapse
|