1
|
Salta Z, Ventura ON, Rais N, Tasinato N, Barone V. A new chapter in the never ending story of cycloadditions: The puzzling case of SO 2 and acetylene. J Comput Chem 2024; 45:1587-1602. [PMID: 38517313 DOI: 10.1002/jcc.27350] [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: 01/30/2024] [Revised: 03/04/2024] [Accepted: 03/08/2024] [Indexed: 03/23/2024]
Abstract
A comprehensive study of the different classes of cycloaddition reactions ([3+2], [2+2], and [2+1]) of SO2 to acetylene and ethylene has been performed using density functional theory (DFT) and composite wavefunction methods. The [3+2] cycloaddition reaction, that was previously explored in the context of the cycloaddition of thioformaldehyde S-methylide (TSM) to ethylene and acetylene, proceeds in a concerted way to the formation of stable heterocycles. In this paper, we extend our study to the [2+2] and [2+1] cycloadditions of SO2 to acetylene, which would produce 1,1-oxathiete-2-oxide and thiirene-1,1-dioxide, respectively. One of the main conclusions is that cyclic 1,1-oxathiete-2-oxide can open through a relatively easy breaking of the SO single bond and rearrange toward sulfinyl acetaldehyde (SA). The SA molecule can easily undergo several internal rearrangements, which eventually lead to sulfenic acid and sulfoxide derivatives of ethenone, 1,2,3-dioxathiole, and CO plus sulfinylmethane. The most probable path, however, produces 2-thioxoacetic acid, whose derivatives (or those of the corresponding acetate) are usually obtained by Willgerodt-Kindler-type sulfuration of acetates. This product can in turn decompose, leading to the final products CO2 and H2CS. Comparison of this decomposition path with that of 2-amino-2-thioxoacetic acid shows that the process occurs through different H-transfer processes.
Collapse
Affiliation(s)
- Zoi Salta
- Scuola Normale Superiore, Pisa, Italy
| | - Oscar N Ventura
- Computational Chemistry and Biology Group, CCBG, DETEMA, Facultad de Química, Universidad de la República, Montevideo, Uruguay
| | - Nadjib Rais
- Scuola Normale Superiore, Pisa, Italy
- IUSS Scuola Universitaria Superiore, Pavia, Italy
| | | | | |
Collapse
|
2
|
Damour Y, Scemama A, Jacquemin D, Kossoski F, Loos PF. State-Specific Coupled-Cluster Methods for Excited States. J Chem Theory Comput 2024; 20:4129-4145. [PMID: 38749498 DOI: 10.1021/acs.jctc.4c00034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.
Collapse
Affiliation(s)
- Yann Damour
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| |
Collapse
|
3
|
Aldossary A, Campos-Gonzalez-Angulo JA, Pablo-García S, Leong SX, Rajaonson EM, Thiede L, Tom G, Wang A, Avagliano D, Aspuru-Guzik A. In Silico Chemical Experiments in the Age of AI: From Quantum Chemistry to Machine Learning and Back. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2024:e2402369. [PMID: 38794859 DOI: 10.1002/adma.202402369] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Revised: 04/28/2024] [Indexed: 05/26/2024]
Abstract
Computational chemistry is an indispensable tool for understanding molecules and predicting chemical properties. However, traditional computational methods face significant challenges due to the difficulty of solving the Schrödinger equations and the increasing computational cost with the size of the molecular system. In response, there has been a surge of interest in leveraging artificial intelligence (AI) and machine learning (ML) techniques to in silico experiments. Integrating AI and ML into computational chemistry increases the scalability and speed of the exploration of chemical space. However, challenges remain, particularly regarding the reproducibility and transferability of ML models. This review highlights the evolution of ML in learning from, complementing, or replacing traditional computational chemistry for energy and property predictions. Starting from models trained entirely on numerical data, a journey set forth toward the ideal model incorporating or learning the physical laws of quantum mechanics. This paper also reviews existing computational methods and ML models and their intertwining, outlines a roadmap for future research, and identifies areas for improvement and innovation. Ultimately, the goal is to develop AI architectures capable of predicting accurate and transferable solutions to the Schrödinger equation, thereby revolutionizing in silico experiments within chemistry and materials science.
Collapse
Affiliation(s)
- Abdulrahman Aldossary
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
| | | | - Sergio Pablo-García
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, ON, M5S 2E4, Canada
| | - Shi Xuan Leong
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
| | - Ella Miray Rajaonson
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- Vector Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto, ON, M5G 1M1, Canada
| | - Luca Thiede
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, ON, M5S 2E4, Canada
- Vector Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto, ON, M5G 1M1, Canada
| | - Gary Tom
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- Vector Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto, ON, M5G 1M1, Canada
| | - Andrew Wang
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
| | - Davide Avagliano
- Chimie ParisTech, PSL University, CNRS, Institute of Chemistry for Life and Health Sciences (iCLeHS UMR 8060), Paris, F-75005, France
| | - Alán Aspuru-Guzik
- Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, ON, M5S 2E4, Canada
- Vector Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto, ON, M5G 1M1, Canada
- Department of Materials Science & Engineering, University of Toronto, 184 College St., Toronto, ON, M5S 3E4, Canada
- Department of Chemical Engineering & Applied Chemistry, University of Toronto, 200 College St., Toronto, ON, M5S 3E5, Canada
- Lebovic Fellow, Canadian Institute for Advanced Research (CIFAR), 66118 University Ave., Toronto, M5G 1M1, Canada
- Acceleration Consortium, 80 St George St, Toronto, M5S 3H6, Canada
| |
Collapse
|
4
|
Daas K, Klute E, Seidl M, Gori-Giorgi P. Møller-Plesset Adiabatic Connection at Large Coupling Strengths for Open-Shell Systems. J Phys Chem A 2024; 128:4138-4149. [PMID: 38717868 PMCID: PMC11129316 DOI: 10.1021/acs.jpca.4c00788] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2024] [Revised: 03/29/2024] [Accepted: 04/05/2024] [Indexed: 05/24/2024]
Abstract
We study the adiabatic connection that has as weak-coupling expansion the Møller-Plesset perturbation series, generalizing to the open-shell case previous closed-shell results for the large-coupling limit. We first focus on the hydrogen atom with fractional spins, providing results along the adiabatic connection from small to large coupling strengths. We reveal an intriguing phase diagram and an equation for the large-coupling leading order that has closed-form solutions for specific choices of its relevant quantum numbers. We then show that the hydrogen atom results provide variational estimates for the large-coupling leading terms for the general many-electron open-shell case in terms of functionals of the Hartree-Fock α-spin and β-spin densities.
Collapse
Affiliation(s)
- Kimberly
J. Daas
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, Amsterdam 1081 HV, The Netherlands
| | - Eveline Klute
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, Amsterdam 1081 HV, The Netherlands
| | - Michael Seidl
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, Amsterdam 1081 HV, The Netherlands
| | - Paola Gori-Giorgi
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, Amsterdam 1081 HV, The Netherlands
- Microsoft
Research AI for Science, Evert van de Beekstraat 354, Schiphol 1118 CZ, The Netherlands
| |
Collapse
|
5
|
Mrovec M, Gill PMW. How delocalized are the polyacenes? J Comput Chem 2024; 45:701-709. [PMID: 38100265 DOI: 10.1002/jcc.27258] [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: 09/30/2023] [Revised: 10/31/2023] [Accepted: 11/01/2023] [Indexed: 12/17/2023]
Abstract
In an attempt to quantify electron delocalization in polyacenes with up to 50 carbon atoms, we have performed self-consistent field calculations in which the π electrons are constrained to occupy highly localized molecular orbitals (HILOs) centered on a maximum of two, six or ten adjacent carbon atoms. We have also performed similar calculations on simple polyacene analogs consisting only of hydrogen atoms and exhibiting electron delocalization in the σ framework. We find that the energetic cost of localizing the π electrons in the polyacenes is roughly 60, 5 or 0.1 kJ/mol per ring atom for the two-, six- and ten-atom HILOs, respectively, and the use of these localized models overestimates the predicted hydrogenation energies of the acenes by roughly 50%, 4% and 0.1%, respectively. We conclude that the chemistry of polyacenes can be modeled well using highly localized descriptions of the π electrons.
Collapse
Affiliation(s)
- Martin Mrovec
- School of Chemistry, University of Sydney, Camperdown, New South Wales, Australia
| | - Peter M W Gill
- School of Chemistry, University of Sydney, Camperdown, New South Wales, Australia
| |
Collapse
|
6
|
Mendolicchio M, Barone V. Unbiased Comparison between Theoretical and Experimental Molecular Structures and Properties: Toward an Accurate Reduced-Cost Evaluation of Vibrational Contributions. J Chem Theory Comput 2024; 20:2842-2857. [PMID: 38556752 DOI: 10.1021/acs.jctc.4c00023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/02/2024]
Abstract
The tremendous development of hardware and software is constantly increasing the role of quantum chemical (QC) computations in the assignment and interpretation of experimental results. However, an unbiased comparison between theory and experiment requires the proper account of vibrational averaging effects. In particular, high-resolution spectra in the gas phase are now available for molecules containing up to about 50 atoms, which are too large for a brute-force approach with the available QC methods of sufficient accuracy. In the present paper, we introduce hybrid approaches, which allow the accurate evaluation of vibrational averaging effects for molecules of this size beyond the harmonic approximation, with special attention being devoted to rotational constants. After the validation of new tools for relatively small molecules, the β-estradiol hormone and a prototypical molecular motor have been considered to witness the feasibility of accurate computations for large molecules.
Collapse
|
7
|
Zhao D, Zhao Y, Xu T, He X, Hu S, Ayers PW, Liu S. Chiral Jahn-Teller Distortion in Quasi-Planar Boron Clusters. Molecules 2024; 29:1624. [PMID: 38611903 PMCID: PMC11013085 DOI: 10.3390/molecules29071624] [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: 02/09/2024] [Revised: 04/03/2024] [Accepted: 04/03/2024] [Indexed: 04/14/2024] Open
Abstract
In this work, we have observed that some chiral boron clusters (B16-, B20-, B24-, and B28-) can simultaneously have helical molecular orbitals and helical spin densities; these seem to be the first compounds discovered to have this intriguing property. We show that chiral Jahn-Teller distortion of quasi-planar boron clusters drives the formation of the helical molecular spin densities in these clusters and show that elongation/enhancement in helical molecular orbitals can be achieved by simply adding more building blocks via a linker. Aromaticity of these boron clusters is discussed. Chiral boron clusters may find potential applications in spintronics, such as molecular magnets.
Collapse
Affiliation(s)
- Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, China
| | - Yilin Zhao
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
| | - Tianlv Xu
- College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha 410081, China
| | - Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, China
| | - Shankai Hu
- Institute of Biomedical Research, Yunnan University, Kunming 650500, China
| | - Paul W. Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
| | - 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
| |
Collapse
|
8
|
Uribe L, Di Grande S, Crisci L, Lazzari F, Mendolicchio M, Barone V. Accurate Structures and Rotational Constants of Steroid Hormones at DFT Cost: Androsterone, Testosterone, Estrone, β-Estradiol, and Estriol. J Phys Chem A 2024; 128:2629-2642. [PMID: 38530336 DOI: 10.1021/acs.jpca.4c00573] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/27/2024]
Abstract
A comprehensive analysis of the structural, conformational, and spectroscopic properties in the gas phase has been performed for five prototypical steroid hormones, namely, androsterone, testosterone, estrone, β-estradiol, and estriol. The revDSD-PBEP86 double-hybrid functional in conjunction with the D3BJ empirical dispersion and a suitable triple-ζ basis set provides accurate conformational energies and equilibrium molecular structures, with the latter being further improved by proper account of core-valence correlation. Average deviations within 0.1% between computed and experimental ground state rotational constants are reached when adding to those equilibrium values vibrational corrections obtained at the cost of standard harmonic frequencies thanks to the use of a new computational tool. Together with the intrinsic interest of the studied hormones, the accuracy of the results obtained at DFT cost for molecules containing about 50 atoms paves the way toward the accurate investigations of other flexible bricks of life.
Collapse
Affiliation(s)
- Lina Uribe
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Scuola Superiore Meridionale, Largo San Marcellino 10, 80138 Napoli, Italy
| | - Silvia Di Grande
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Scuola Superiore Meridionale, Largo San Marcellino 10, 80138 Napoli, Italy
| | - Luigi Crisci
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Federico Lazzari
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Mendolicchio
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | |
Collapse
|
9
|
Liu S. Harvesting Chemical Understanding with Machine Learning and Quantum Computers. ACS PHYSICAL CHEMISTRY AU 2024; 4:135-142. [PMID: 38560751 PMCID: PMC10979482 DOI: 10.1021/acsphyschemau.3c00067] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Revised: 12/29/2023] [Accepted: 01/02/2024] [Indexed: 04/04/2024]
Abstract
It is tenable to argue that nobody can predict the future with certainty, yet one can learn from the past and make informed projections for the years ahead. In this Perspective, we overview the status of how theory and computation can be exploited to obtain chemical understanding from wave function theory and density functional theory, and then outlook the likely impact of machine learning (ML) and quantum computers (QC) to appreciate traditional chemical concepts in decades to come. It is maintained that the development and maturation of ML and QC methods in theoretical and computational chemistry represent two paradigm shifts about how the Schrödinger equation can be solved. New chemical understanding can be harnessed in these two new paradigms by making respective use of ML features and QC qubits. Before that happens, however, we still have hurdles to face and obstacles to overcome in both ML and QC arenas. Possible pathways to tackle these challenges are proposed. We anticipate that hierarchical modeling, in contrast to multiscale modeling, will emerge and thrive, becoming the workhorse of in silico simulations in the next few decades.
Collapse
|
10
|
Patra A, Pipim GB, Krylov AI, Mallikarjun Sharada S. Performance of Density Functionals for Excited-State Properties of Isolated Chromophores and Exciplexes: Emission Spectra, Solvatochromic Shifts, and Charge-Transfer Character. J Chem Theory Comput 2024; 20:2520-2537. [PMID: 38488640 DOI: 10.1021/acs.jctc.4c00005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/27/2024]
Abstract
This study assesses the performance of various meta-generalized gradient approximation (meta-GGA), global hybrid, and range-separated hybrid (RSH) density functionals in capturing the excited-state properties of organic chromophores and their excited-state complexes (exciplexes). Motivated by their uses in solar energy harvesting and photoredox CO2 reduction, we use oligo-(p-phenylenes) and their excited-state complexes with triethylamine as model systems. We focus on the fluorescence properties of these systems, specifically emission energies. We also consider solvatochromic shifts and wave function characteristics. The latter is described by using reduced quantities such as natural transition orbitals (NTOs) and exciton descriptors. The functionals are benchmarked against the experimental fluorescence spectra and the equation-of-motion coupled-cluster method with single and double excitations. Both in isolated chromophores and in exciplexes, meta-GGA functionals drastically underestimate the emission energies and exhibit significant exciton delocalization and anticorrelation between electron and hole motion. The performance of global hybrid functionals is strongly dependent on the percentage of exact exchange. Our study identifies RSH GGAs as the best-performing functionals, with ωPBE demonstrating the best agreement with experimental results. RSH meta-GGAs often overestimate emission energies in exciplexes and yield larger hole NTOs. Their performance can be improved by optimally tuning the range-separation parameter.
Collapse
Affiliation(s)
- Abhilash Patra
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089, United States
| | - George Baffour Pipim
- Department of Chemistry, University of Southern California, Los Angeles ,California 90089, United States
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles ,California 90089, United States
| | - Shaama Mallikarjun Sharada
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089, United States
- Department of Chemistry, University of Southern California, Los Angeles ,California 90089, United States
| |
Collapse
|
11
|
Zhao D, Zhao Y, Xu E, Liu W, Ayers PW, Liu S, Chen D. Fragment-Based Deep Learning for Simultaneous Prediction of Polarizabilities and NMR Shieldings of Macromolecules and Their Aggregates. J Chem Theory Comput 2024; 20:2655-2665. [PMID: 38441881 DOI: 10.1021/acs.jctc.3c01415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2024]
Abstract
Simultaneous prediction of the molecular response properties, such as polarizability and the NMR shielding constant, at a low computational cost is an unresolved issue. We propose to combine a linear-scaling generalized energy-based fragmentation (GEBF) method and deep learning (DL) with both molecular and atomic information-theoretic approach (ITA) quantities as effective descriptors. In GEBF, the total molecular polarizability can be assembled as a linear combination of the corresponding quantities calculated from a set of small embedded subsystems in GEBF. In the new GEBF-DL(ITA) protocol, one can predict subsystem polarizabilities based on the corresponding molecular wave function (thus electron density and ITA quantities) and DL model rather than calculate them from the computationally intensive coupled-perturbed Hartree-Fock or Kohn-Sham equations and finally obtain the total molecular polarizability via a linear combination equation. As a proof-of-concept application, we predict the molecular polarizabilities of large proteins and protein aggregates. GEBF-DL(ITA) is shown to be as accurate enough as GEBF, with mean absolute percentage error <1%. For the largest protein aggregate (>4000 atoms), GEBF-DL(ITA) gains a speedup ratio of 3 compared with GEBF. It is anticipated that when more advanced electronic structure methods are used, this advantage will be more appealing. Moreover, one can also predict the NMR chemical shieldings of proteins with reasonably good accuracy. Overall, the cost-efficient GEBF-DL(ITA) protocol should be a robust theoretical tool for simultaneously predicting polarizabilities and NMR shieldings of large systems.
Collapse
Affiliation(s)
- Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming, Yunnan 650500, P. R. China
| | - Yilin Zhao
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton ONL8S4M1, Canada
| | - Enhua Xu
- Graduate School of System Informatics, Kobe University, Nada-ku, Kobe, Hyogo 657-8501, Japan
| | - Wenqi Liu
- Institute of Biomedical Research, Yunnan University, Kunming, Yunnan 650500, P. R. China
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton ONL8S4M1, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| | - Dahua Chen
- Institute of Biomedical Research, Yunnan University, Kunming, Yunnan 650500, P. R. China
| |
Collapse
|
12
|
Panchagnula K, Graf D, Albertani FEA, Thom AJW. Translational eigenstates of He@C60 from four-dimensional ab initio potential energy surfaces interpolated using Gaussian process regression. J Chem Phys 2024; 160:104303. [PMID: 38465682 DOI: 10.1063/5.0197903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2024] [Accepted: 02/22/2024] [Indexed: 03/12/2024] Open
Abstract
We investigate the endofullerene system 3He@C60 with a four-dimensional potential energy surface (PES) to include the three He translational degrees of freedom and C60 cage radius. We compare second order Møller-Plesset perturbation theory (MP2), spin component scaled-MP2, scaled opposite spin-MP2, random phase approximation (RPA)@Perdew, Burke, and Ernzerhof (PBE), and corrected Hartree-Fock-RPA to calibrate and gain confidence in the choice of electronic structure method. Due to the high cost of these calculations, the PES is interpolated using Gaussian Process Regression (GPR), owing to its effectiveness with sparse training data. The PES is split into a two-dimensional radial surface, to which corrections are applied to achieve an overall four-dimensional surface. The nuclear Hamiltonian is diagonalized to generate the in-cage translational/vibrational eigenstates. The degeneracy of the three-dimensional harmonic oscillator energies with principal quantum number n is lifted due to the anharmonicity in the radial potential. The (2l + 1)-fold degeneracy of the angular momentum states is also weakly lifted, due to the angular dependence in the potential. We calculate the fundamental frequency to range between 96 and 110 cm-1 depending on the electronic structure method used. Error bars of the eigenstate energies were calculated from the GPR and are on the order of ∼±1.5 cm-1. Wavefunctions are also compared by considering their overlap and Hellinger distance to the one-dimensional empirical potential. As with the energies, the two ab initio methods MP2 and RPA@PBE show the best agreement. While MP2 has better agreement than RPA@PBE, due to its higher computational efficiency and comparable performance, we recommend RPA as an alternative electronic structure method of choice to MP2 for these systems.
Collapse
Affiliation(s)
- K Panchagnula
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - D Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - F E A Albertani
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - A J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
13
|
Wodyński A, Lauw B, Reimann M, Kaupp M. Spin-Symmetry Breaking and Hyperfine Couplings in Transition-Metal Complexes Revisited Using Density Functionals Based on the Exact-Exchange Energy Density. J Chem Theory Comput 2024; 20:2033-2048. [PMID: 38411554 PMCID: PMC10938646 DOI: 10.1021/acs.jctc.3c01422] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Revised: 02/02/2024] [Accepted: 02/07/2024] [Indexed: 02/28/2024]
Abstract
A small set of mononuclear manganese complexes evaluated previously for their Mn hyperfine couplings (HFCs) has been analyzed using density functionals based on the exact-exchange energy density─in particular, the spin symmetry breaking (SSB) found previously when using hybrid functionals. Employing various strong-correlation corrected local hybrids (scLHs) and strong-correlation corrected range-separated local hybrids (scRSLHs) with or without additional corrections to their local mixing functions (LMFs) to mitigate delocalization errors (DE), the SSB and the associated dipolar HFCs of [Mn(CN)4]2-, MnO3, [Mn(CN)4N]-, and [Mn(CN)5NO]2- (the latter with cluster embedding) have been examined. Both strong-correlation (sc)-correction and DE-correction terms help to diminish SSB and correct the dipolar HFCs. The DE corrections are more effective, and the effects of the sc corrections depend on their damping factors. Interestingly, the DE-corrections reduce valence-shell spin polarization (VSSP) and thus SSB by locally enhancing exact-exchange (EXX) admixture near the metal center and thereby diminishing spin-density delocalization onto the ligand atoms. In contrast, sc corrections diminish EXX admixture locally, mostly on specific ligand atoms. This then reduces VSSP and SSB as well. The performance of scLHs and scRSLHs for the isotropic Mn HFCs has also been analyzed, with particular attention to core-shell spin-polarization contributions. Further sc-corrected functionals, such as the KP16/B13 construction and the DM21 deep-neural-network functional, have been examined.
Collapse
Affiliation(s)
- Artur Wodyński
- Technische Universität
Berlin, Institut für Chemie, Theoretische
Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, Berlin, D-10623, Germany
| | - Bryan Lauw
- Technische Universität
Berlin, Institut für Chemie, Theoretische
Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, Berlin, D-10623, Germany
| | - Marc Reimann
- Technische Universität
Berlin, Institut für Chemie, Theoretische
Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, Berlin, D-10623, Germany
| | - Martin Kaupp
- Technische Universität
Berlin, Institut für Chemie, Theoretische
Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, Berlin, D-10623, Germany
| |
Collapse
|
14
|
Mortensen JJ, Larsen AH, Kuisma M, Ivanov AV, Taghizadeh A, Peterson A, Haldar A, Dohn AO, Schäfer C, Jónsson EÖ, Hermes ED, Nilsson FA, Kastlunger G, Levi G, Jónsson H, Häkkinen H, Fojt J, Kangsabanik J, Sødequist J, Lehtomäki J, Heske J, Enkovaara J, Winther KT, Dulak M, Melander MM, Ovesen M, Louhivuori M, Walter M, Gjerding M, Lopez-Acevedo O, Erhart P, Warmbier R, Würdemann R, Kaappa S, Latini S, Boland TM, Bligaard T, Skovhus T, Susi T, Maxson T, Rossi T, Chen X, Schmerwitz YLA, Schiøtz J, Olsen T, Jacobsen KW, Thygesen KS. GPAW: An open Python package for electronic structure calculations. J Chem Phys 2024; 160:092503. [PMID: 38450733 DOI: 10.1063/5.0182685] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 01/15/2024] [Indexed: 03/08/2024] Open
Abstract
We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.
Collapse
Affiliation(s)
- Jens Jørgen Mortensen
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Ask Hjorth Larsen
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Mikael Kuisma
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Aleksei V Ivanov
- Riverlane Ltd., St Andrews House, 59 St Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Alireza Taghizadeh
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Andrew Peterson
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - Anubhab Haldar
- Department of Electrical and Computer Engineering, Boston University, Boston, Massachusetts 02215, USA
| | - Asmus Ougaard Dohn
- Department of Physics, Technical University of Denmark, 2800 Lyngby, Denmark and Science Institute and Faculty of Physical Sciences, VR-III, University of Iceland, Reykjavík 107, Iceland
| | - Christian Schäfer
- Department of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
| | - Elvar Örn Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Eric D Hermes
- Quantum-Si, 29 Business Park Drive, Branford, Connecticut 06405, USA
| | | | - Georg Kastlunger
- CatTheory, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Gianluca Levi
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Hannu Häkkinen
- Departments of Physics and Chemistry, Nanoscience Center, University of Jyväskylä, FI-40014 Jyväskylä, Finland
| | - Jakub Fojt
- Department of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
| | - Jiban Kangsabanik
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Joachim Sødequist
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Jouko Lehtomäki
- Department of Applied Physics, Aalto University, P.O. Box 11100, 00076 Aalto, Finland
| | - Julian Heske
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Jussi Enkovaara
- CSC-IT Center for Science Ltd., P.O. Box 405, FI-02101 Espoo, Finland
| | - Kirsten Trøstrup Winther
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Marcin Dulak
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Marko M Melander
- Department of Chemistry, Nanoscience Center, University of Jyväskylä, FI-40014 Jyväskylä, Finland
| | - Martin Ovesen
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Martti Louhivuori
- CSC-IT Center for Science Ltd., P.O. Box 405, FI-02101 Espoo, Finland
| | - Michael Walter
- FIT Freiburg Centre for Interactive Materials and Bioinspired Technologies, University of Freiburg, Georges-Köhler-Allee 105, 79110 Freiburg, Germany
| | - Morten Gjerding
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Olga Lopez-Acevedo
- Biophysics of Tropical Diseases, Max Planck Tandem Group, University of Antioquia UdeA, 050010 Medellin, Colombia
| | - Paul Erhart
- Department of Physics, Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
| | - Robert Warmbier
- School of Physics and Mandelstam Institute for Theoretical Physics, University of the Witwatersrand, 1 Jan Smuts Avenue, 2001 Johannesburg, South Africa
| | - Rolf Würdemann
- Freiburger Materialforschungszentrum, Universität Freiburg, Stefan-Meier-Straße 21, D-79104 Freiburg, Germany
| | - Sami Kaappa
- Computational Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| | - Simone Latini
- Nanomade, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Tara Maria Boland
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Thomas Bligaard
- Department of Energy Conversion and Storage, Technical University of Denmark, DK-2800 Lyngby, Denmark
| | - Thorbjørn Skovhus
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Toma Susi
- Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
| | - Tristan Maxson
- Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, Alabama 35487, USA
| | - Tuomas Rossi
- CSC-IT Center for Science Ltd., P.O. Box 405, FI-02101 Espoo, Finland
| | - Xi Chen
- School of Physical Science and Technology, Lanzhou University, Lanzhou, Gansu 730000, China
| | | | - Jakob Schiøtz
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Thomas Olsen
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | | | | |
Collapse
|
15
|
Zhang H, Liu S, You J, Liu C, Zheng S, Lu Z, Wang T, Zheng N, Shao B. Overcoming the barrier of orbital-free density functional theory for molecular systems using deep learning. NATURE COMPUTATIONAL SCIENCE 2024; 4:210-223. [PMID: 38467870 DOI: 10.1038/s43588-024-00605-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Accepted: 02/07/2024] [Indexed: 03/13/2024]
Abstract
Orbital-free density functional theory (OFDFT) is a quantum chemistry formulation that has a lower cost scaling than the prevailing Kohn-Sham DFT, which is increasingly desired for contemporary molecular research. However, its accuracy is limited by the kinetic energy density functional, which is notoriously hard to approximate for non-periodic molecular systems. Here we propose M-OFDFT, an OFDFT approach capable of solving molecular systems using a deep learning functional model. We build the essential non-locality into the model, which is made affordable by the concise density representation as expansion coefficients under an atomic basis. With techniques to address unconventional learning challenges therein, M-OFDFT achieves a comparable accuracy to Kohn-Sham DFT on a wide range of molecules untouched by OFDFT before. More attractively, M-OFDFT extrapolates well to molecules much larger than those seen in training, which unleashes the appealing scaling of OFDFT for studying large molecules including proteins, representing an advancement of the accuracy-efficiency trade-off frontier in quantum chemistry.
Collapse
Affiliation(s)
- He Zhang
- National Key Laboratory of Human-Machine Hybrid Augmented Intelligence, National Engineering Research Center for Visual Information and Applications, and Institute of Artificial Intelligence and Robotics, Xi'an Jiaotong University, Xi'an, China
- Microsoft Research AI4Science, Beijing, China
| | - Siyuan Liu
- Microsoft Research AI4Science, Beijing, China
| | | | - Chang Liu
- Microsoft Research AI4Science, Beijing, China.
| | | | - Ziheng Lu
- Microsoft Research AI4Science, Beijing, China
| | - Tong Wang
- Microsoft Research AI4Science, Beijing, China
| | - Nanning Zheng
- National Key Laboratory of Human-Machine Hybrid Augmented Intelligence, National Engineering Research Center for Visual Information and Applications, and Institute of Artificial Intelligence and Robotics, Xi'an Jiaotong University, Xi'an, China
| | - Bin Shao
- Microsoft Research AI4Science, Beijing, China.
| |
Collapse
|
16
|
Laestadius A, Csirik MA, Penz M, Tancogne-Dejean N, Ruggenthaler M, Rubio A, Helgaker T. Exchange-only virial relation from the adiabatic connection. J Chem Phys 2024; 160:084115. [PMID: 38421067 DOI: 10.1063/5.0184934] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Accepted: 01/25/2024] [Indexed: 03/02/2024] Open
Abstract
The exchange-only virial relation due to Levy and Perdew is revisited. Invoking the adiabatic connection, we introduce the exchange energy in terms of the right-derivative of the universal density functional w.r.t. the coupling strength λ at λ = 0. This agrees with the Levy-Perdew definition of the exchange energy as a high-density limit of the full exchange-correlation energy. By relying on v-representability for a fixed density at varying coupling strength, we prove an exchange-only virial relation without an explicit local-exchange potential. Instead, the relation is in terms of a limit (λ ↘ 0) involving the exchange-correlation potential vxcλ, which exists by assumption of v-representability. On the other hand, a local-exchange potential vx is not warranted to exist as such a limit.
Collapse
Affiliation(s)
- Andre Laestadius
- Department of Computer Science, Oslo Metropolitan University, 0130 Oslo, Norway
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, 0315 Oslo, Norway
| | - Mihály A Csirik
- Department of Computer Science, Oslo Metropolitan University, 0130 Oslo, Norway
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, 0315 Oslo, Norway
| | - Markus Penz
- Department of Computer Science, Oslo Metropolitan University, 0130 Oslo, Norway
- Max Planck Institute for the Structure and Dynamics of Matter and Center for Free-Electron Laser Science and Department of Physics, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Nicolas Tancogne-Dejean
- Max Planck Institute for the Structure and Dynamics of Matter and Center for Free-Electron Laser Science and Department of Physics, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Michael Ruggenthaler
- Max Planck Institute for the Structure and Dynamics of Matter and Center for Free-Electron Laser Science and Department of Physics, Luruper Chaussee 149, 22761 Hamburg, Germany
- The Hamburg Center for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Angel Rubio
- Max Planck Institute for the Structure and Dynamics of Matter and Center for Free-Electron Laser Science and Department of Physics, Luruper Chaussee 149, 22761 Hamburg, Germany
- The Hamburg Center for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, USA
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, 0315 Oslo, Norway
| |
Collapse
|
17
|
He X, Li M, Rong C, Zhao D, Liu W, Ayers PW, Liu S. Some Recent Advances in Density-Based Reactivity Theory. J Phys Chem A 2024; 128:1183-1196. [PMID: 38329898 DOI: 10.1021/acs.jpca.3c07997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Establishing a chemical reactivity theory in density functional theory (DFT) language has been our intense research interest in the past two decades, exemplified by the determination of steric effect and stereoselectivity, evaluation of electrophilicity and nucleophilicity, identification of strong and weak interactions, and formulation of cooperativity, frustration, and principle of chirality hierarchy. In this Featured Article, we first overview the four density-based frameworks in DFT to appreciate chemical understanding, including conceptual DFT, use of density associated quantities, information-theoretic approach, and orbital-free DFT, and then present a few recent advances of these frameworks as well as new applications from our studies. To that end, we will introduce the relationship among these frameworks, determining the entire spectrum of interactions with Pauli energy derivatives, performing topological analyses with information-theoretic quantities, and extending the density-based frameworks to excited states. Applications to examine physiochemical properties in external electric fields and to evaluate polarizability for proteins and crystals are discussed. A few possible directions for future development are followed, with the special emphasis on its merger with machine learning.
Collapse
Affiliation(s)
- Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Meng Li
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Chunying Rong
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton ONL8S, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| |
Collapse
|
18
|
Hölzer C, Gordiy I, Grimme S, Bursch M. Hybrid DFT Geometries and Properties for 17k Lanthanoid Complexes─The LnQM Data Set. J Chem Inf Model 2024; 64:825-836. [PMID: 38238264 DOI: 10.1021/acs.jcim.3c01832] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/13/2024]
Abstract
The unique properties of lanthanoids and their diverse applications make them an indispensable part of modern research and industry. While the field has garnered attention, there remains a gap in available molecule data sets that facilitate both classical quantum chemistry calculations and the burgeoning field of machine learning in data science applications. This research addresses the need for a comprehensive data set that allows for a comparative analysis of various lanthanoids. The herein presented, curated data set includes 17269 monolanthanoid complexes derived from 1205 distinct ligand motifs. Structures encompass all 15 lanthanoids in the +3 oxidation state and exhibit molecular charges ranging from -1 to +3, including structures with a high spin multiplicity up to 8. Starting from lanthanum complexes, samples were processed with a permutation of the central lanthanoid atom, resulting in highly comparable subsets, facilitating comparative studies in which the influence of the lanthanoid can be investigated independently of ligand effects. The data set provides a broad range of features such as PBE0-D4/def2-SVP optimized geometries and optimization trajectories, while also covering ωB97M-V/def2-SVPD energies, rotational constants, dipole moments, highest occupied molecular orbital-lowest-unoccupied molecular orbital (HOMO-LUMO) energies, and Mulliken, Löwdin, and Hirshfeld population analyses. Additionally, coordination numbers, polarizabilities, and partial charges from D4, electronegativity equilibration (EEQ), GFN2-xTB, and charge extended Hückel (CEH) calculations are included. The data set is openly accessible and may serve as a basis for further investigations into the properties of lanthanoids.
Collapse
Affiliation(s)
- Christian Hölzer
- Mulliken Center for Theoretical Chemistry, University of Bonn, 53115 Bonn, Germany
| | - Igor Gordiy
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, Zürich 8093, Switzerland
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, 53115 Bonn, Germany
| | - Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
19
|
Wang B, Geerlings P, Liu S, De Proft F. Extending the Scope of Conceptual Density Functional Theory with Second Order Analytical Methodologies. J Chem Theory Comput 2024. [PMID: 38310523 DOI: 10.1021/acs.jctc.3c01184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2024]
Abstract
In the context of the growing impact of conceptual density functional theory (DFT) as one of the most successful chemical reactivity theories, response functions up to second order have now been widely applied; in recent years, among others, particular attention has been focused on the linear response function and also extensions to higher order have been put forward. As the larger part of these studies have been carried using a finite difference approach to compute these concepts, we now embarked on (an extension of) an analytical approach to conceptual DFT. With the ultimate aim of providing a complete set of analytically computable second order properties, including the softness and hardness kernels, the hardness as the simplest second order response function is scrutinized again with numerical results highlighting the difference in nature between the analytical hardness (referred to as hardness condition) and the Parr-Pearson absolute chemical hardness. The hardness condition is investigated for its capability to gauge the (de)localization error of density functional approximations (DFAs). The analytical Fukui function, besides overcoming the difficulties in the finite difference approach in treating negatively charged systems, also showcases the errors of deviating from the straight-line behavior using fractional occupation number calculations. Subsequently, the softness kernel and its atom-condensed inverse, the hardness matrix, are accessed through the Berkowitz-Parr relation. Revisiting the softness kernel confirms and extends previous discussions on how Kohn's Nearsightedness of Electronic Matter principle can be retrieved and identified as the physicist's version of the chemist's "transferability of functional groups" concept. The accurate, analytical hardness matrix evaluation on the other hand provides further support for the basics of Nalewajski's charge sensitivity analysis. Based on Parr and Liu's functional expansion of the energy functional, a new energy decomposition is introduced with an order of magnitude analysis of the different terms for a series of simple molecules both at their equilibrium geometry and upon variation in bond length and dihedral angle. Finally, for the first time, the perturbation expansion of the energy functional is studied numerically up to second order now that all response functions and integration techniques are at hand. The perturbation expansion energies are in excellent agreement with those obtained directly from DFA calculations giving confidence in the convergence of the perturbation series and its use in judging the importance of the different terms in reactivity investigations.
Collapse
Affiliation(s)
- Bin Wang
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| | - Paul Geerlings
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| | - Frank De Proft
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
| |
Collapse
|
20
|
Ayers PW. Energy is not a convex function of particle number for r-k interparticle potentials with k > log34. J Chem Phys 2024; 160:044110. [PMID: 38275194 DOI: 10.1063/5.0179137] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Accepted: 12/11/2023] [Indexed: 01/27/2024] Open
Abstract
The energy of a many-particle system is not convex with respect to particle number for r-k interparticle repulsion potentials if k > log34 ≈ 1.262. With such potentials, some finite electronic systems have ionization potentials that are less than the electron affinity: they have negative band gap (chemical hardness). Although the energy may be a convex function of the number of electrons (for which k = 1), it suggests that finding an analytic proof of convexity will be very difficult. The bound on k is postulated to be tight. An apparent signature of non-convex behavior is that the Dyson orbital corresponding to the lowest-energy mode of electron attachment has a vanishingly small amplitude.
Collapse
Affiliation(s)
- Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, Ontario L8S 4M1, Canada
| |
Collapse
|
21
|
Seenithurai S, Chai JD. Electronic Properties of Graphene Nano-Parallelograms: A Thermally Assisted Occupation DFT Computational Study. Molecules 2024; 29:349. [PMID: 38257262 DOI: 10.3390/molecules29020349] [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: 12/18/2023] [Revised: 01/05/2024] [Accepted: 01/08/2024] [Indexed: 01/24/2024] Open
Abstract
In this computational study, we investigate the electronic properties of zigzag graphene nano-parallelograms (GNPs), which are parallelogram-shaped graphene nanoribbons of various widths and lengths, using thermally assisted occupation density functional theory (TAO-DFT). Our calculations revealed a monotonic decrease in the singlet-triplet energy gap as the GNP length increased. The GNPs possessed singlet ground states for all the cases examined. With the increase of GNP length, the vertical ionization potential and fundamental gap decreased monotonically, while the vertical electron affinity increased monotonically. Some of the GNPs studied were found to possess fundamental gaps in the range of 1-3 eV, lying in the ideal region relevant to solar energy applications. Besides, as the GNP length increased, the symmetrized von Neumann entropy increased monotonically, denoting an increase in the degree of the multi-reference character associated with the ground state GNPs. The occupation numbers and real-space representation of active orbitals indicated that there was a transition from the nonradical nature of the shorter GNPs to the increasing polyradical nature of the longer GNPs. In addition, the edge/corner localization of the active orbitals was found for the wider and longer GNPs.
Collapse
Affiliation(s)
- Sonai Seenithurai
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Jeng-Da Chai
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
- Center for Theoretical Physics and Center for Quantum Science and Engineering, National Taiwan University, Taipei 10617, Taiwan
- Physics Division, National Center for Theoretical Sciences, Taipei 10617, Taiwan
| |
Collapse
|
22
|
Huynh BC, Wibowo-Teale M, Wibowo-Teale AM. QSym 2: A Quantum Symbolic Symmetry Analysis Program for Electronic Structure. J Chem Theory Comput 2024; 20:114-133. [PMID: 38145888 PMCID: PMC10782455 DOI: 10.1021/acs.jctc.3c01118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Revised: 11/16/2023] [Accepted: 11/29/2023] [Indexed: 12/27/2023]
Abstract
Symmetry provides a powerful machinery to classify, interpret, and understand quantum-mechanical theories and results. However, most contemporary quantum chemistry packages lack the ability to handle degeneracy and symmetry breaking effects, especially in non-Abelian groups, and they are not able to characterize symmetry in the presence of external magnetic or electric fields. In this article, a program written in Rust entitled QSym2 that makes use of group and representation theories to provide symmetry analysis for a wide range of quantum-chemical calculations is introduced. With its ability to generate character tables symbolically on-the-fly and by making use of a generic symmetry-orbit-based representation analysis method formulated in this work, QSym2 is able to address all of these shortcomings. To illustrate these capabilities of QSym2, four sets of case studies are examined in detail in this article: (i) high-symmetry C84H64, C60, and B9- to demonstrate the analysis of degenerate molecular orbitals (MOs); (ii) octahedral Fe(CN)63- to demonstrate the analysis of symmetry-broken determinants and MOs; (iii) linear hydrogen fluoride in a magnetic field to demonstrate the analysis of magnetic symmetry; and (iv) equilateral H3+ to demonstrate the analysis of density symmetries.
Collapse
Affiliation(s)
- Bang C. Huynh
- School
of Chemistry, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Meilani Wibowo-Teale
- School
of Chemistry, University of Nottingham, Nottingham NG7 2RD, United Kingdom
| | - Andrew M. Wibowo-Teale
- School
of Chemistry, University of Nottingham, Nottingham NG7 2RD, United Kingdom
- Hylleraas
Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| |
Collapse
|
23
|
Robles-Navarro A, Jerabek P, Schwerdtfeger P. Tipping the Balance Between the bcc and fcc Phase Within the Alkali and Coinage Metal Groups. Angew Chem Int Ed Engl 2024; 63:e202313679. [PMID: 37877444 DOI: 10.1002/anie.202313679] [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: 09/14/2023] [Revised: 10/24/2023] [Accepted: 10/25/2023] [Indexed: 10/26/2023]
Abstract
Why the Group 1 elements crystallize in the body-centered cubic (bcc) structure, and the iso-electronic Group 11 elements in the face-centered cubic (fcc) structure, remains a mystery. Here we show that a delicate interplay between many-body effects, vibrational contributions and dispersion interactions obtained from relativistic density functional theory offers an answer to this long-standing controversy. It also sheds light on the Periodic Table of Crystal Structures. A smooth diffusionless transition through cuboidal lattices gives a detailed insight into the bcc→fcc phase transition for the Groups 1 and 11 elements.
Collapse
Affiliation(s)
- Andrés Robles-Navarro
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study (NZIAS), Massey University Albany, Private Bag 102904, Auckland, 0745, New Zealand
| | - Paul Jerabek
- Institute of Hydrogen Technology, Helmholtz-Zentrum Hereon, Max-Planck-Str. 1, 21502, Geesthacht, Germany
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study (NZIAS), Massey University Albany, Private Bag 102904, Auckland, 0745, New Zealand
| |
Collapse
|
24
|
Ponra A, Bakasa C, Etindele AJ, Casida ME. Diagrammatic multiplet sum method (MSM) density functional theory (DFT): Investigation of the transferability of integrals in "simple" DFT-based approaches to multideterminantal problems. J Chem Phys 2023; 159:244306. [PMID: 38149739 DOI: 10.1063/5.0173572] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Accepted: 11/30/2023] [Indexed: 12/28/2023] Open
Abstract
Kohn-Sham density functional theory (DFT) typically works well for describing dynamic correlation. Two other types of correlation, arising in the cases of degenerate (static) or quasidegenerate (nondynamic) zero-order states, represent a difficult problem for DFT. When symmetry is present, multiplet sum method (MSM) DFT [Ziegler et al., Theor. Chim. Acta 4, 877 (1977)] provides one of the earliest and simplest ways to include static correlation in DFT. MSM-DFT assumes that DFT provides a good description of single-determinant energies and uses symmetry and simple ansätze to include the effects of static correlation. This is equivalent to determining the off-diagonal matrix elements in a small configuration interaction (CI) eigenvalue problem. Our ultimate goal, however, is nondynamic correlation in cases where symmetry is inadequate for fixing the dynamic-correlation limitation of DFT. To this end, we have developed a diagrammatic approach to MSM-DFT, which does not, by itself, solve the nondynamic correlation problem in DFT but which facilitates comparison with wave function CI and so allows educated guesses of off-diagonal CI matrix elements even in the absence of symmetry. In every case, an additional exchange-only ansatz (EXAN) allows the MSM-DFT formulas to be transformed into wave function formulas. This EXAN also works for transforming time-dependent DFT into time-dependent Hartree-Fock. Although not enough to uniquely guess DFT formulas from wave function formulas, the diagrammatic approach and the EXAN provide important constraints on any guesses that might be used. We illustrate how diagrammatic MSM-DFT may be used to guess a nondynamic correlation correction for the dissociation of H2 and how diagrammatic MSM-DFT may be used to guess a nonsymmetry-based coupling element in the O2 multiplet problem, which is reasonably close to a previous symmetry-derived result.
Collapse
Affiliation(s)
- Abraham Ponra
- Department of Physics, Faculty of Science, University of Maroua, P.O. Box 814, Maroua, Cameroon
| | - Carolyne Bakasa
- Technical University of Kenya, P.O. Box 52428-00200, Haile Selassie Avenue, Nairobe, Kenya
| | - Anne Justine Etindele
- Higher Teachers Training College, University of Yaounde I, P.O. Box 47, Yaounde, Cameroon
| | - Mark E Casida
- Laboratoire de Spectrométrie, Interactions et Chimie théorique (SITh), Département de Chimie Moléculaire (DCM, UMR CNRS/UGA 5250), Institut de Chimie Moléculaire de Grenoble (ICMG, FR2607), Université Grenoble Alpes (UGA) 301 rue de la Chimie, BP 53, F-38041 Grenoble Cedex 9, France
| |
Collapse
|
25
|
Alonso M, Bettens T, Eeckhoudt J, Geerlings P, De Proft F. Wandering through quantum-mechanochemistry: from concepts to reactivity and switches. Phys Chem Chem Phys 2023; 26:21-35. [PMID: 38086672 DOI: 10.1039/d3cp04907h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
Mechanochemistry has experienced a renaissance in recent years witnessing, at the molecular level, a remarkable interplay between theory and experiment. Molecular mechanochemistry has welcomed a broad spectrum of quantum-chemical methods to evaluate the influence of an external mechanical force on molecular properties. In this contribution, an overview is given on recent work on quantum mechanochemistry in the Brussels Quantum Chemistry group (ALGC). The effect of an external force was scrutinized both in fundamental topics, like reactivity descriptors in Conceptual DFT, and in applied topics, such as designing molecular force probes and tuning the stereoselectivity of certain types of reactions. In the conceptual part, a brief overview of the techniques introducing mechanical forces into a quantum-mechanical description of a molecule is followed by an introduction to conceptual DFT. The evolution of the electronic chemical potential (or electronegativity), chemical hardness and electrophilicity are investigated when a chemical bond in a series of diatomics is put under mechanical stress. Its counterpart, the influence of mechanical stress on bond angles, is analyzed by varying the strain present in alkyne triple bonds by applying a bending force, taking the strain promoted alkyne-azide coupling cycloaddition as an example. The increase of reactivity of the alkyne upon bending is probed by Fukui functions and the local softness. In the applied part, a new molecular force probe is presented based on an intramolecular 6π-electrocyclization in constrained polyenes operating under thermal conditions. A cyclic process is conceived where ring opening and closure are triggered by applying or removing an external pulling force. The efficiency of mechanical activation strongly depends on the magnitude of the applied force and the distance between the pulling points. The idea of pulling point distances as a tool to identify new mechanochemical processes is then tested in [28]hexaphyrins with an intricate equilibrium between Möbius aromatic and Hückel antiaromatic topologies. A mechanical force is shown to trigger the interconversion between the two topologies, using the distance matrix as a guide to select appropriate pulling points. In a final application, the Felkin-Anh model for the addition of nucleophiles to chiral carbonyls under the presence of an external mechanical force is scrutinized. By applying a force for restricting the conformational freedom of the chiral ketone, otherwise inaccessible reaction pathways are promoted on the force-modified potential energy surfaces resulting in a diastereoselectivity different from the force-free reaction.
Collapse
Affiliation(s)
- Mercedes Alonso
- Eenheid Algemene Chemie (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium.
| | - Tom Bettens
- Eenheid Algemene Chemie (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium.
| | - Jochen Eeckhoudt
- Eenheid Algemene Chemie (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium.
| | - Paul Geerlings
- Eenheid Algemene Chemie (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium.
| | - Frank De Proft
- Eenheid Algemene Chemie (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium.
| |
Collapse
|
26
|
Daas KJ, Kooi DP, Benyahia T, Seidl M, Gori-Giorgi P. Large-Z atoms in the strong-interaction limit of DFT: Implications for gradient expansions and for the Lieb-Oxford bound. J Chem Phys 2023; 159:234114. [PMID: 38112505 DOI: 10.1063/5.0174592] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 11/26/2023] [Indexed: 12/21/2023] Open
Abstract
We numerically study the strong-interaction limit of the exchange-correlation functional for neutral atoms and Bohr atoms as the number of electrons increases. Using a compact representation, we analyze the second-order gradient expansion, comparing it with the one for exchange (weak interaction limit). The two gradient expansions, at strong and weak interaction, turn out to be very similar in magnitude but with opposite signs. We find that the point-charge plus continuum model is surprisingly accurate for the gradient expansion coefficient at strong coupling, while generalized gradient approximations, such as Perdew-Burke-Ernzerhof (PBE) and PBEsol, severely underestimate it. We then use our results to analyze the Lieb-Oxford bound from the point of view of slowly varying densities, clarifying some aspects on the bound at a fixed number of electrons.
Collapse
Affiliation(s)
- Kimberly J Daas
- 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
| | - Derk P Kooi
- 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
- Microsoft Research AI4Science, Amsterdam, The Netherlands
| | - Tarik Benyahia
- 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
| | - Michael Seidl
- 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
| | - 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
- Microsoft Research AI4Science, Amsterdam, The Netherlands
| |
Collapse
|
27
|
Burgess AC, Linscott E, O'Regan DD. The convexity condition of density-functional theory. J Chem Phys 2023; 159:211102. [PMID: 38038199 DOI: 10.1063/5.0174159] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2023] [Accepted: 11/08/2023] [Indexed: 12/02/2023] Open
Abstract
It has long been postulated that within density-functional theory (DFT), the total energy of a finite electronic system is convex with respect to electron count so that 2Ev[N0] ≤ Ev[N0 - 1] + Ev[N0 + 1]. Using the infinite-separation-limit technique, this Communication proves the convexity condition for any formulation of DFT that is (1) exact for all v-representable densities, (2) size-consistent, and (3) translationally invariant. An analogous result is also proven for one-body reduced density matrix functional theory. While there are known DFT formulations in which the ground state is not always accessible, indicating that convexity does not hold in such cases, this proof, nonetheless, confirms a stringent constraint on the exact exchange-correlation functional. We also provide sufficient conditions for convexity in approximate DFT, which could aid in the development of density-functional approximations. This result lifts a standing assumption in the proof of the piecewise linearity condition with respect to electron count, which has proven central to understanding the Kohn-Sham bandgap and the exchange-correlation derivative discontinuity of DFT.
Collapse
Affiliation(s)
- Andrew C Burgess
- School of Physics, Trinity College Dublin, The University of Dublin, Dublin, Ireland
| | - Edward Linscott
- Theory and Simulation of Materials (THEOS), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - David D O'Regan
- School of Physics, Trinity College Dublin, The University of Dublin, Dublin, Ireland
| |
Collapse
|
28
|
Penz M, Tellgren EI, Csirik MA, Ruggenthaler M, Laestadius A. The Structure of the Density-Potential Mapping. Part II: Including Magnetic Fields. ACS PHYSICAL CHEMISTRY AU 2023; 3:492-511. [PMID: 38034040 PMCID: PMC10683500 DOI: 10.1021/acsphyschemau.3c00006] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Revised: 07/13/2023] [Accepted: 07/13/2023] [Indexed: 12/02/2023]
Abstract
The Hohenberg-Kohn theorem of density-functional theory (DFT) is broadly considered the conceptual basis for a full characterization of an electronic system in its ground state by just one-body particle density. In this Part II of a series of two articles, we aim at clarifying the status of this theorem within different extensions of DFT including magnetic fields. We will in particular discuss current-density-functional theory (CDFT) and review the different formulations known in the literature, including the conventional paramagnetic CDFT and some nonstandard alternatives. For the former, it is known that the Hohenberg-Kohn theorem is no longer valid due to counterexamples. Nonetheless, paramagnetic CDFT has the mathematical framework closest to standard DFT and, just like in standard DFT, nondifferentiability of the density functional can be mitigated through Moreau-Yosida regularization. Interesting insights can be drawn from both Maxwell-Schrödinger DFT and quantum-electrodynamic DFT, which are also discussed here.
Collapse
Affiliation(s)
- Markus Penz
- Basic
Research Community for Physics, Innsbruck 6020, Austria
| | - Erik I. Tellgren
- Hylleraas
Centre for Quantum Molecular Sciences, University
of Oslo, Oslo 0315, Norway
| | - Mihály A. Csirik
- Hylleraas
Centre for Quantum Molecular Sciences, University
of Oslo, Oslo 0315, Norway
- Department
of Computer Science, Oslo Metropolitan University, Oslo 0130, Norway
| | - Michael Ruggenthaler
- Max
Planck Institute for the Structure and Dynamics of Matter, Hamburg 22761, Germany
| | - Andre Laestadius
- Hylleraas
Centre for Quantum Molecular Sciences, University
of Oslo, Oslo 0315, Norway
- Department
of Computer Science, Oslo Metropolitan University, Oslo 0130, Norway
| |
Collapse
|
29
|
Zhang W, He X, Li M, Zhang J, Zhao D, Liu S, Rong C. Simultaneous identification of strong and weak interactions with Pauli energy, Pauli potential, Pauli force, and Pauli charge. J Chem Phys 2023; 159:184104. [PMID: 37942871 DOI: 10.1063/5.0173666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Accepted: 10/12/2023] [Indexed: 11/10/2023] Open
Abstract
Strong and weak interatomic interactions in chemical and biological systems are ubiquitous, yet how to identify them on a unified theoretical foundation is still not well established. Recently, we proposed employing Pauli energy-based indexes, such as strong covalent interaction and bonding and noncovalent interaction indexes, in the framework of density functional theory for the purpose. In this work, we extend our previous theoretical work by directly employing Pauli energy, Pauli potential, Pauli force, and Pauli charge to simultaneously identify both strong covalent bonding and weak noncovalent interactions. Our results from this work elucidate that using their signature isosurfaces, we can identify different types of interactions, either strong or weak, including single, double, triple, and quadruple covalent bonds, ionic bond, metallic bond, hydrogen bonding, and van der Waals interaction. We also discovered strong linear correlations between Pauli energy derived quantities and different covalent bond orders. These qualitative and quantitative results from our present study solidify the viewpoint that a unified approach to simultaneously identify both strong and weak interactions is possible. In our view, this work signifies one step forward towards the goal of establishing a density-based theory of chemical reactivity in density functional theory.
Collapse
Affiliation(s)
- Wenbiao Zhang
- Key Laboratory of Light Energy Conversion Materials of Hunan Province College, Hunan Normal University, Changsha 410081, Hunan, People's Republic of China
| | - Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Meng Li
- Key Laboratory of Light Energy Conversion Materials of Hunan Province College, Hunan Normal University, Changsha 410081, Hunan, People's Republic of China
| | - Jingwen Zhang
- Key Laboratory of Light Energy Conversion Materials of Hunan Province College, Hunan Normal University, Changsha 410081, Hunan, People's Republic of China
| | - Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, People's Republic of China
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, USA
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, USA
| | - Chunying Rong
- Key Laboratory of Light Energy Conversion Materials of Hunan Province College, Hunan Normal University, Changsha 410081, Hunan, People's Republic of China
| |
Collapse
|
30
|
Hostaš J, Pérez-Becerra KO, Calaminici P, Barrios-Herrera L, Lourenço MP, Tchagang A, Salahub DR, Köster AM. How important is the amount of exact exchange for spin-state energy ordering in DFT? Case study of molybdenum carbide cluster, Mo4C2. J Chem Phys 2023; 159:184301. [PMID: 37947508 DOI: 10.1063/5.0169409] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Accepted: 10/23/2023] [Indexed: 11/12/2023] Open
Abstract
Since the form of the exact functional in density functional theory is unknown, we must rely on density functional approximations (DFAs). In the past, very promising results have been reported by combining semi-local DFAs with exact, i.e. Hartree-Fock, exchange. However, the spin-state energy ordering and the predictions of global minima structures are particularly sensitive to the choice of the hybrid functional and to the amount of exact exchange. This has been already qualitatively described for single conformations, reactions, and a limited number of conformations. Here, we have analyzed the mixing of exact exchange in exchange functionals for a set of several hundred isomers of the transition metal carbide, Mo4C2. The analysis of the calculated energies and charges using PBE0-type functional with varying amounts of exact exchange yields the following insights: (1) The sensitivity of spin-energy splitting is strongly correlated with the amount of exact exchange mixing. (2) Spin contamination is exacerbated when correlation is omitted from the exchange-correlation functional. (3) There is not one ideal value for the exact exchange mixing which can be used to parametrize or choose among the functionals. Calculated energies and electronic structures are influenced by exact exchange at a different magnitude within a given distribution; therefore, to extend the application range of hybrid functionals to the full periodic table the spin-energy splitting energies should be investigated.
Collapse
Affiliation(s)
- Jiří Hostaš
- Department of Chemistry, CMS - Centre for Molecular Simulation, IQST - Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Kevin O Pérez-Becerra
- Departamento de Química, Cinvestav, Avenida Instituto Politécnico Nacional 2508, A.P. 14-740, CDMX C.P. 07360, Mexico
| | - Patrizia Calaminici
- Departamento de Química, Cinvestav, Avenida Instituto Politécnico Nacional 2508, A.P. 14-740, CDMX C.P. 07360, Mexico
| | - Lizandra Barrios-Herrera
- Department of Chemistry, CMS - Centre for Molecular Simulation, IQST - Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Maicon Pierre Lourenço
- Departamento de Química e Física - Centro de Ciências Exatas, Naturais e da Saúde - CCENS - Universidade Federal do Espírito Santo, 29500-000 Alegre, Espírito Santo, Brazil
| | - Alain Tchagang
- Digital Technologies Research Centre, National Research Council of Canada, 1200 Montréal Road, Ottawa, Ontario K1A 0R6, Canada
| | - Dennis R Salahub
- Department of Chemistry, CMS - Centre for Molecular Simulation, IQST - Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Andreas M Köster
- Departamento de Química, Cinvestav, Avenida Instituto Politécnico Nacional 2508, A.P. 14-740, CDMX C.P. 07360, Mexico
| |
Collapse
|
31
|
Langer MF, Frank JT, Knoop F. Stress and heat flux via automatic differentiation. J Chem Phys 2023; 159:174105. [PMID: 37921248 DOI: 10.1063/5.0155760] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2023] [Accepted: 09/25/2023] [Indexed: 11/04/2023] Open
Abstract
Machine-learning potentials provide computationally efficient and accurate approximations of the Born-Oppenheimer potential energy surface. This potential determines many materials properties and simulation techniques usually require its gradients, in particular forces and stress for molecular dynamics, and heat flux for thermal transport properties. Recently developed potentials feature high body order and can include equivariant semi-local interactions through message-passing mechanisms. Due to their complex functional forms, they rely on automatic differentiation (AD), overcoming the need for manual implementations or finite-difference schemes to evaluate gradients. This study discusses how to use AD to efficiently obtain forces, stress, and heat flux for such potentials, and provides a model-independent implementation. The method is tested on the Lennard-Jones potential, and then applied to predict cohesive properties and thermal conductivity of tin selenide using an equivariant message-passing neural network potential.
Collapse
Affiliation(s)
- Marcel F Langer
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
- BIFOLD-Berlin Institute for the Foundations of Learning and Data, Berlin, Germany
- The NOMAD Laboratory at the Fritz Haber Institute of the Max Planck Society and Humboldt University, Berlin, Germany
| | - J Thorben Frank
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
- BIFOLD-Berlin Institute for the Foundations of Learning and Data, Berlin, Germany
| | - Florian Knoop
- Theoretical Physics Division, Department of Physics, Chemistry and Biology (IFM), Linköping University, SE-581 83 Linköping, Sweden
| |
Collapse
|
32
|
Tsai HY, Chai JD. Real-Time Extension of TAO-DFT. Molecules 2023; 28:7247. [PMID: 37959667 PMCID: PMC10647330 DOI: 10.3390/molecules28217247] [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: 09/27/2023] [Revised: 10/19/2023] [Accepted: 10/21/2023] [Indexed: 11/15/2023] Open
Abstract
Thermally assisted occupation density functional theory (TAO-DFT) has been an efficient electronic structure method for studying the ground-state properties of large electronic systems with multi-reference character over the past few years. To explore the time-dependent (TD) properties of electronic systems (e.g., subject to an intense laser pulse), in this work, we propose a real-time (RT) extension of TAO-DFT, denoted as RT-TAO-DFT. Moreover, we employ RT-TAO-DFT to study the high-order harmonic generation (HHG) spectra and related TD properties of molecular hydrogen H2 at the equilibrium and stretched geometries, aligned along the polarization of an intense linearly polarized laser pulse. The TD properties obtained with RT-TAO-DFT are compared with those obtained with the widely used time-dependent Kohn-Sham (TDKS) method. In addition, issues related to the possible spin-symmetry breaking effects in the TD properties are discussed.
Collapse
Affiliation(s)
- Hung-Yi Tsai
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan;
| | - Jeng-Da Chai
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan;
- Center for Theoretical Physics and Center for Quantum Science and Engineering, National Taiwan University, Taipei 10617, Taiwan
- Physics Division, National Center for Theoretical Sciences, Taipei 10617, Taiwan
| |
Collapse
|
33
|
Di Felice R, Mayes ML, Richard RM, Williams-Young DB, Chan GKL, de Jong WA, Govind N, Head-Gordon M, Hermes MR, Kowalski K, Li X, Lischka H, Mueller KT, Mutlu E, Niklasson AMN, Pederson MR, Peng B, Shepard R, Valeev EF, van Schilfgaarde M, Vlaisavljevich B, Windus TL, Xantheas SS, Zhang X, Zimmerman PM. A Perspective on Sustainable Computational Chemistry Software Development and Integration. J Chem Theory Comput 2023; 19:7056-7076. [PMID: 37769271 PMCID: PMC10601486 DOI: 10.1021/acs.jctc.3c00419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Indexed: 09/30/2023]
Abstract
The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.
Collapse
Affiliation(s)
- Rosa Di Felice
- Departments
of Physics and Astronomy and Quantitative and Computational Biology, University of Southern California, Los Angeles, California 90089, United States
- CNR-NANO
Modena, Modena 41125, Italy
| | - Maricris L. Mayes
- Department
of Chemistry and Biochemistry, University
of Massachusetts Dartmouth, North Dartmouth, Massachusetts 02747, United States
| | | | | | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Wibe A. de Jong
- Lawrence
Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Niranjan Govind
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Martin Head-Gordon
- Pitzer Center
for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Karol Kowalski
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409, United States
| | - Karl T. Mueller
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced
Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Anders M. N. Niklasson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Mark R. Pederson
- Department
of Physics, The University of Texas at El
Paso, El Paso, Texas 79968, United States
| | - Bo Peng
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, United States
| | - Edward F. Valeev
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | | | - Bess Vlaisavljevich
- Department
of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Theresa L. Windus
- Department
of Chemistry, Iowa State University and
Ames Laboratory, Ames, Iowa 50011, United States
| | - Sotiris S. Xantheas
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Advanced
Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Xing Zhang
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Paul M. Zimmerman
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
34
|
Bosko IP, Staroverov VN. Derivation and reinterpretation of the Fermi-Amaldi functional. J Chem Phys 2023; 159:131101. [PMID: 37800642 DOI: 10.1063/5.0166358] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2023] [Accepted: 09/13/2023] [Indexed: 10/07/2023] Open
Abstract
The Fermi-Amaldi correction to the electrostatic self-repulsion of the particle density is usually regarded as a semi-classical exchange functional that happens to be exact only for one- and closed-shell two-electron systems. We show that this functional can be derived quantum-mechanically and is exact for any number of fermions or bosons of arbitrary spin as long as the particles occupy the same spatial orbital. The Fermi-Amaldi functional is also size-consistent for such systems, provided that the factor N in its expression is understood as an orbital occupation number rather than the total number of particles. These properties of the Fermi-Amaldi functional are ultimately related to the fact that it is a special case of the self-exchange energy formula. Implications of our findings are discussed.
Collapse
Affiliation(s)
- Ivan P Bosko
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Viktor N Staroverov
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
35
|
Giarrusso S, Loos PF. Exact Excited-State Functionals of the Asymmetric Hubbard Dimer. J Phys Chem Lett 2023; 14:8780-8786. [PMID: 37739406 PMCID: PMC10561271 DOI: 10.1021/acs.jpclett.3c02052] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2023] [Accepted: 08/31/2023] [Indexed: 09/24/2023]
Abstract
The exact functionals associated with the (singlet) ground state and the two singlet excited states of the asymmetric Hubbard dimer at half-filling are calculated using both Levy's constrained search and Lieb's convex formulation. While the ground-state functional is, as is commonly known, a convex function with respect to the density, the functional associated with the doubly excited state is found to be concave. Also, because the density-potential mapping associated with the first excited state is noninvertible, its "functional" is a partial, multivalued function composed of one concave and one convex branch that correspond to two separate domains of the external potential. Remarkably, it is found that, although the one-to-one mapping between density and external potential may not apply (as in the case of the first excited state), each state-specific energy and corresponding universal functional are "functions" whose derivatives are each other's inverse, just as in the ground state formalism.
Collapse
Affiliation(s)
- Sara Giarrusso
- Laboratoire de Chimie et
Physique Quantiques (UMR 5626), Université
de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et
Physique Quantiques (UMR 5626), Université
de Toulouse, CNRS, UPS, 31062 Toulouse, France
| |
Collapse
|
36
|
Hermann J, Spencer J, Choo K, Mezzacapo A, Foulkes WMC, Pfau D, Carleo G, Noé F. Ab initio quantum chemistry with neural-network wavefunctions. Nat Rev Chem 2023; 7:692-709. [PMID: 37558761 DOI: 10.1038/s41570-023-00516-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/16/2023] [Indexed: 08/11/2023]
Abstract
Deep learning methods outperform human capabilities in pattern recognition and data processing problems and now have an increasingly important role in scientific discovery. A key application of machine learning in molecular science is to learn potential energy surfaces or force fields from ab initio solutions of the electronic Schrödinger equation using data sets obtained with density functional theory, coupled cluster or other quantum chemistry (QC) methods. In this Review, we discuss a complementary approach using machine learning to aid the direct solution of QC problems from first principles. Specifically, we focus on quantum Monte Carlo methods that use neural-network ansatzes to solve the electronic Schrödinger equation, in first and second quantization, computing ground and excited states and generalizing over multiple nuclear configurations. Although still at their infancy, these methods can already generate virtually exact solutions of the electronic Schrödinger equation for small systems and rival advanced conventional QC methods for systems with up to a few dozen electrons.
Collapse
Affiliation(s)
- Jan Hermann
- Microsoft Research AI4Science, Berlin, Germany
- FU Berlin, Department of Mathematics and Computer Science, Berlin, Germany
| | | | - Kenny Choo
- Department of Physics, University of Zurich, Zurich, Switzerland
- IBM Quantum, IBM Research Zurich, Ruschlikon, Switzerland
| | | | - W M C Foulkes
- Imperial College London, Department of Physics, London, UK
| | - David Pfau
- DeepMind, London, UK.
- Imperial College London, Department of Physics, London, UK.
| | | | - Frank Noé
- Microsoft Research AI4Science, Berlin, Germany.
- FU Berlin, Department of Mathematics and Computer Science, Berlin, Germany.
- FU Berlin, Department of Physics, Berlin, Germany.
- Department of Chemistry,Rice University, Houston, TX, USA.
| |
Collapse
|
37
|
Zhao D, Zhao Y, He X, Li Y, Ayers PW, Liu S. Accurate and Efficient Prediction of Post-Hartree-Fock Polarizabilities of Condensed-Phase Systems. J Chem Theory Comput 2023; 19:6461-6470. [PMID: 37676647 DOI: 10.1021/acs.jctc.3c00646] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
To accurately and efficiently predict the molecular response properties (such as polarizability) at post-Hartree-Fock levels for condensed-phase systems under periodic boundary conditions (PBC) is still an unaccomplished and ongoing task. We demonstrate that static isotropic polarizabilities can be cost-effectively predicted at post-Hartree-Fock levels by combining the linear-scaling generalized energy-based fragmentation (GEBF) and information-theoretic approach (ITA) quantities. In PBC-GEBF, the total molecular polarizability of an extended system is obtained as a linear combination of the corresponding quantities of a series of small embedded subsystems of several monomers. Here, we show that in the PBC-GEBF-ITA framework, one can obtain the molecular polarizabilities and establish linear relations to ITA quantities. Once these relations are established for smaller subsystems, one can predict the polarizabilities of larger subsystems directly from the molecular wavefunction (or electron density) via ITA quantities. Alternatively, one can determine the total molecular polarizability via a linear combination equation in PBC-GEBF. We have corroborated that this newly proposed PBC-GEBF-ITA protocol is much more efficient than the original PBC-GEBF approach but is not much less accurate and that this conclusion holds for both many-body perturbation theory and the coupled cluster calculations. Good efficiency and transferability of the PBC-GEBF-ITA protocol are demonstrated for periodic systems with several hundred atoms in a unit cell.
Collapse
Affiliation(s)
- Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, P. R. China
| | - Yilin Zhao
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
| | - Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, P. R. China
| | - Yunzhi Li
- School of Chemistry and Chemical Engineering, Linyi University, Linyi 276000, P. R. China
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| |
Collapse
|
38
|
Sitha S. Better performance of Hartree-Fock over DFT: a quantum mechanical investigation on pyridinium benzimidazolate types of zwitterions in the light of localization/delocalization issues. J Mol Model 2023; 29:313. [PMID: 37704866 PMCID: PMC10499969 DOI: 10.1007/s00894-023-05706-1] [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: 07/19/2023] [Accepted: 08/22/2023] [Indexed: 09/15/2023]
Abstract
CONTEXT With the advent of fast computing facilities, combined with rapid emerges of many new and intricate quantum mechanical functionals, computations with pure Hartree-Fock (HF) theory are now-a-days regarded as trivial or obsolete, or even considered as not reliable by many researchers. Consequently, current trends in computational chemistry show extensive use of post-HF theories for smaller molecular systems and various DFT methods for organic and inorganic chemistry related problems (larger molecules/systems). In this contribution, I have tried to show that sometimes, HF might be more suitable over DFT methodologies in addressing structure-property correlations. Molecules studied here were previously synthesized by Boyd in 1966 and important experimental data were produced by Alcalde and co-workers in 1987. Comparison of computed and experimental results clearly shows that HF method was more effective in reproducing the experimental data compared to especially the DFT methodologies. Reliability of HF method was further assured from the very similar results shown by the CCSD, CASSCF, CISD and QCISD methods. Current study also indicates that the localization issue associated with HF proved to be advantageous over delocalization issue of DFT based methodologies, in correctly describing the structure-property correlation for zwitterion systems. METHODS All computations were performed with Gaussian 09. A wide-range of quantum mechanical methodologies, HF, B3LYP, CAM-B3LYP, BMK, B3PW91, TPSSh, LC-ωPBE, M06-2X, M06-HF, ωB97xD, MP2, CASSCF, CCSD, QCISD, CISD and semi-empirical methods like, Huckel, CNDO, AM1, PM3MM and PM6, were used for investigations.
Collapse
Affiliation(s)
- Sanyasi Sitha
- Department of Chemical Sciences, University of Johannesburg, Auckland Park, PO Box 524, Johannesburg, 2006, South Africa.
| |
Collapse
|
39
|
Zhan J, Govoni M, Galli G. Nonempirical Range-Separated Hybrid Functional with Spatially Dependent Screened Exchange. J Chem Theory Comput 2023; 19:5851-5862. [PMID: 37591004 DOI: 10.1021/acs.jctc.3c00580] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/19/2023]
Abstract
Electronic structure calculations based on density functional theory (DFT) have successfully predicted numerous ground-state properties of a variety of molecules and materials. However, exchange and correlation functionals currently used in the literature, including semilocal and hybrid functionals, are often inaccurate to describe the electronic properties of heterogeneous solids, especially systems composed of building blocks with large dielectric mismatch. Here, we present a dielectric-dependent range-separated hybrid functional, screened-exchange range-separated hybrid (SE-RSH), for the investigation of heterogeneous materials. We define a spatially dependent fraction of exact exchange inspired by the static Coulomb-hole and screened-exchange (COHSEX) approximation used in many-body perturbation theory, and we show that the proposed functional accurately predicts the electronic structure of several nonmetallic interfaces, three- and two-dimensional, pristine, and defective solids and nanoparticles.
Collapse
Affiliation(s)
- Jiawei Zhan
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Govoni
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Department of Physics, Computer Science, and Mathematics, University of Modena and Reggio Emilia, Modena 41125, Italy
| | - Giulia Galli
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Materials Science Division and Center for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
40
|
Acke G, Van Hende D, De Vriendt X, Bultinck P. Extending Conceptual Density Functional Theory toward First-Order Reduced Density Matrices: An Open Subsystems Viewpoint on the Fukui Matrix. J Chem Theory Comput 2023; 19:5418-5426. [PMID: 37531218 DOI: 10.1021/acs.jctc.3c00080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2023]
Abstract
As a matrix extension of the Fukui function, a reactivity descriptor grounded within Conceptual Density Functional Theory, the Fukui matrix extends Frontier Molecular Orbital Theory to correlated regimes with its eigendecomposition in Fukui occupations and Fukui naturals. Despite successful applications, the questions remain as to whether replacing a quantity derived from a purely density-based framework by its matrix extension is theoretically well-founded and what chemical information is contained in the corresponding eigendecomposition. In this study, we show that the matrix extension of the Fukui function is only well-defined if one also generalizes the external potential to become nonlocal, leading to the introduction of Conceptual First-Order Reduced Density Matrix Functional Theory. By interpreting the Anderson impurity model from an interacting open subsystem perspective, we show how Fukui occupations and Fukui naturals reflect the influence of an increasing (static) correlation and which characteristic patterns we should expect within a molecular context. This study represents a step in generalizing Conceptual Density Functional Theory beyond its density-based perspective.
Collapse
Affiliation(s)
- Guillaume Acke
- Department of Chemistry - Ghent Quantum Chemistry Group, Ghent University, Krijgslaan 281 (S3), Ghent B-9000, Belgium
| | - Daria Van Hende
- Department of Chemistry - Ghent Quantum Chemistry Group, Ghent University, Krijgslaan 281 (S3), Ghent B-9000, Belgium
| | - Xeno De Vriendt
- Department of Chemistry - Ghent Quantum Chemistry Group, Ghent University, Krijgslaan 281 (S3), Ghent B-9000, Belgium
| | - Patrick Bultinck
- Department of Chemistry - Ghent Quantum Chemistry Group, Ghent University, Krijgslaan 281 (S3), Ghent B-9000, Belgium
| |
Collapse
|
41
|
He X, Lu T, Rong C, Liu S, Ayers PW, Liu W. Topological analysis of information-theoretic quantities in density functional theory. J Chem Phys 2023; 159:054112. [PMID: 37548307 DOI: 10.1063/5.0159941] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2023] [Accepted: 07/06/2023] [Indexed: 08/08/2023] Open
Abstract
We have witnessed considerable research interest in the recent literature about the development and applications of quantities from the information-theoretic approach (ITA) in density functional theory. These ITA quantities are explicit density functionals, whose local distributions in real space are continuous and well-behaved. In this work, we further develop ITA by systematically analyzing the topological behavior of its four representative quantities, Shannon entropy, two forms of Fisher information, and relative Shannon entropy (also called information gain or Kullback-Leibler divergence). Our results from their topological analyses for 103 molecular systems provide new insights into bonding interactions and physiochemical properties, such as electrophilicity, nucleophilicity, acidity, and aromaticity. We also compare our results with those from the electron density, electron localization function, localized orbital locator, and Laplacian functions. Our results offer a new methodological approach and practical tool for applications that are especially promising for elucidating chemical bonding and reactivity propensity.
Collapse
Affiliation(s)
- Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Tian Lu
- Beijing Kein Research Center for Natural Sciences, Beijing 100022, China
| | - Chunying Rong
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, USA
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, USA
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, Ontario L8S 4M1, Canada
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| |
Collapse
|
42
|
Penz M, Tellgren EI, Csirik MA, Ruggenthaler M, Laestadius A. The Structure of Density-Potential Mapping. Part I: Standard Density-Functional Theory. ACS PHYSICAL CHEMISTRY AU 2023; 3:334-347. [PMID: 37520314 PMCID: PMC10375883 DOI: 10.1021/acsphyschemau.2c00069] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Revised: 02/18/2023] [Accepted: 02/21/2023] [Indexed: 08/01/2023]
Abstract
The Hohenberg-Kohn theorem of density-functional theory (DFT) is broadly considered the conceptual basis for a full characterization of an electronic system in its ground state by just the one-body particle density. Part I of this review aims at clarifying the status of the Hohenberg-Kohn theorem within DFT and Part II at different extensions of the theory that include magnetic fields. We collect evidence that the Hohenberg-Kohn theorem does not so much form the basis of DFT, but is rather the consequence of a more comprehensive mathematical framework. Such results are especially useful when it comes to the construction of generalized DFTs.
Collapse
Affiliation(s)
- Markus Penz
- Basic
Research Community for Physics, Innsbruck 6020, Austria
| | - Erik I. Tellgren
- Hylleraas
Centre for Quantum Molecular Sciences, University
of Oslo, Oslo 0315, Norway
| | - Mihály A. Csirik
- Hylleraas
Centre for Quantum Molecular Sciences, University
of Oslo, Oslo 0315, Norway
- Department
of Computer Science, Oslo Metropolitan University, Oslo 0130, Norway
| | - Michael Ruggenthaler
- Max
Planck Institute for the Structure and Dynamics of Matter, Hamburg 22761, Germany
| | - Andre Laestadius
- Hylleraas
Centre for Quantum Molecular Sciences, University
of Oslo, Oslo 0315, Norway
- Department
of Computer Science, Oslo Metropolitan University, Oslo 0130, Norway
| |
Collapse
|
43
|
Ariai J, Gellrich U. The entropic penalty for associative reactions and their physical treatment during routine computations. Phys Chem Chem Phys 2023; 25:14005-14015. [PMID: 37161492 DOI: 10.1039/d3cp00970j] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
A systematic study of the entropic penalty for associative reactions is presented. It is shown that computed solution-phase Gibbs free energies typically overestimate entropic contributions. This entropic penalty for associative reactions in solution, i.e., if the number of particles decreases along the reaction coordinate (sum of stoichiometric numbers ), originates from the insufficient treatment of entropic effects by implicit solvent models. We propose an additive correction scheme to Gibbs free energies that is suitable for routine applications by non-expert users. This correction is based on Garza's formalism for the solution-phase entropy [A. J. Garza, J. Chem. Theory Comput., 2019, 15, 3204.] that is physically sound and embedded into an efficient black-box type algorithm. To critically evaluate the entropic penalty and its proposed treatment, we compiled an experimental benchmark set of 31 ΔrG and 22 in 15 different solvents. Using a representative best-practice computational protocol (at wave function theory (WFT) based DLPNO-CCSD(T) and density functional theory (DFT) based revDSD-PBEP86-D4 level with an implicit solvent model), we determined a sizeable entropic penalty ranging from 2-11 kcal mol-1. Using the correction scheme presented herein, the entropic penalty is corrected to the chemical accuracy of ≤1 kcal mol-1 (WFT and DFT). The same applies to at the WFT level. Barriers at the DFT level are overestimated by 2 kcal mol-1 (classic) and underestimated by 2 kcal mol-1 (corrected). This effect is attributed to the finding that barriers computed at the DFT level are systematically 2-3 kcal mol-1 lower than barriers obtained with WFT.
Collapse
Affiliation(s)
- Jama Ariai
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany.
| | - Urs Gellrich
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany.
| |
Collapse
|
44
|
Wang B, Rong C, Lei M, Liu S, De Proft F. Mechanistic Study and Conceptual Chemical Reactivity Analysis of Hydroboration of Carbon Dioxide Catalyzed by a Manganese(I)-PNP-Pincer Complex. Inorg Chem 2023; 62:7366-7375. [PMID: 37129648 DOI: 10.1021/acs.inorgchem.3c00553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Designing efficient and selective catalysts for carbon dioxide reduction is an intensive research area in the recent literature on homogeneous catalysis. In this work, we study the catalytic activity of a newly reported Mn(I)-PNP-pincer catalyst with an embedded aromatic ring. First, we systematically examine its capability to yield different products and highlight the importance of ligand aromaticity and steric effects on metal-ligand cooperativity. We then further conceptually probe its reactivity with descriptors from both conceptual density functional theory and an information-theoretic approach, thereby proposing a novel partitioning of the reaction coordinate into three relevant regions. Our results show that the reactivity in these different regions is governed by different properties such as steric effects, electrophilicity/nucleophilicity, or aromaticity. We anticipate that this methodology, with the analytical tools employed in this study, can be generalized and extended to other catalytic systems and find applications in designing better catalysts.
Collapse
Affiliation(s)
- Bin Wang
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha 410081, Hunan, P. R. China
| | - Chunying Rong
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha 410081, Hunan, P. R. China
| | - Ming Lei
- State Key Laboratory of Chemical Resource Engineering, Institute of Computational Chemistry, College of Chemistry, Beijing University of Chemical Technology, Beijing 100029, China
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| | - Frank De Proft
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium
| |
Collapse
|
45
|
Kim JH, Kim D, Yang W, Baik MH. Fractional Charge Density Functional Theory and Its Application to the Electro-inductive Effect. J Phys Chem Lett 2023; 14:3329-3334. [PMID: 36989527 DOI: 10.1021/acs.jpclett.3c00323] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
We employed the chemical potential equalization principle to demonstrate that fractional electrons are involved in the electro-inductive effect as well as the vibrational Stark effect. By the chemical potential model, we were able to deduce that the frontier molecular orbitals of immobilized molecules can provide valuable insight into these effects. To further understand and quantify these findings, we introduced fractional charge density functional theory (FC-DFT), a canonical ensemble approach for open systems. This method allows for the calculation of electronic energies, nuclear gradients, and the Hessian matrix of fractional electronic systems. To correct the spurious delocalization error commonly found in approximate density functionals for small systems, we imposed the Perdew-Parr-Levy-Balduz (PPLB) condition through linear interpolation of two adjacent integer points (LI-FC-DFT). Although this approach is relatively simple in terms of molecular modeling, the results obtained through LI-FC-DFT calculations predict the same trend seen in experimental reactivity and the frequency change of immobilized molecules.
Collapse
Affiliation(s)
- Jun-Hyeong Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
- Center for Catalytic Hydrocarbon Functionalizations, Institute for Basic Science (IBS), Daejeon 34141, Republic of Korea
| | - Dongju Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
- Center for Catalytic Hydrocarbon Functionalizations, Institute for Basic Science (IBS), Daejeon 34141, Republic of Korea
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Mu-Hyun Baik
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
- Center for Catalytic Hydrocarbon Functionalizations, Institute for Basic Science (IBS), Daejeon 34141, Republic of Korea
| |
Collapse
|
46
|
Kotov AA, Kozhedub YS, Glazov DA, Iliaš M, Pershina V, Shabaev VM. Relativistic Coupled-Cluster Calculations of Spectroscopic Properties of Copernicium and Flerovium Monoxides. Chemphyschem 2023; 24:e202200680. [PMID: 36383485 DOI: 10.1002/cphc.202200680] [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: 09/08/2022] [Revised: 11/11/2022] [Accepted: 11/13/2022] [Indexed: 11/17/2022]
Abstract
Calculations of spectroscopic properties of the CnO and FlO molecules are performed using ab initio all-electron 4c- and 2c-relativistic coupled-cluster approaches with single, double, and perturbative triple excitations. The corresponding calculation for HgO is also accomplished for comparison with the published data. The dependence of the results on the parameters of the basis set and approximations used is investigated in detail. The overall relative uncertainties of the recommended values on the level of 1-2 % are reached. The calculated spectroscopic constants are indicative of the following trend in the reactivity of the oxides HgO>FlO>CnO. This is confirmed by the trend in the adsorption energies, Eads , of these molecules on the surfaces of gold, quartz, and Teflon. The predicted rather low Eads values for the latter case should guarantee their delivery from the recoil chamber to the chemistry set up in gas-phase experiments.
Collapse
Affiliation(s)
- Artem A Kotov
- Department of Physics, St. Petersburg State University, 199034, St. Petersburg, Russia
| | - Yury S Kozhedub
- Department of Physics, St. Petersburg State University, 199034, St. Petersburg, Russia
| | - Dmitry A Glazov
- Department of Physics, St. Petersburg State University, 199034, St. Petersburg, Russia
| | - Miroslav Iliaš
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University, Tajovského 40, 97401, Banská Bystrica, Slovakia.,Helmholtz-Institut Mainz, Johannes Gutenberg-Universität, 55099, Mainz, Germany.,GSI Helmholtzzentrum für Schwerionenforschung GmbH, Planckstr. 1, 64291, Darmstadt, Germany
| | - Valeria Pershina
- GSI Helmholtzzentrum für Schwerionenforschung GmbH, Planckstr. 1, 64291, Darmstadt, Germany
| | - Vladimir M Shabaev
- Department of Physics, St. Petersburg State University, 199034, St. Petersburg, Russia
| |
Collapse
|
47
|
Zhao D, He X, Ayers PW, Liu S. Excited-State Polarizabilities: A Combined Density Functional Theory and Information-Theoretic Approach Study. Molecules 2023; 28:molecules28062576. [PMID: 36985548 PMCID: PMC10058485 DOI: 10.3390/molecules28062576] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2023] [Revised: 03/07/2023] [Accepted: 03/10/2023] [Indexed: 03/14/2023] Open
Abstract
Accurate and efficient determination of excited-state polarizabilities (α) is an open problem both experimentally and computationally. Following our previous work, (Phys. Chem. Chem. Phys. 2023, 25, 2131−2141), in which we employed simple ground-state (S0) density-related functions from the information-theoretic approach (ITA) to accurately and efficiently evaluate the macromolecular polarizabilities, in this work we aimed to predict the lowest excited-state (S1) polarizabilities. The philosophy is to use density-based functions to depict excited-state polarizabilities. As a proof-of-principle application, employing 2-(2′-hydroxyphenyl)benzimidazole (HBI), its substituents, and some other commonly used ESIPT (excited-state intramolecular proton transfer) fluorophores as model systems, we verified that either with S0 or S1 densities as an input, ITA quantities can be strongly correlated with the excited-state polarizabilities. When transition densities are considered, both S0 and S1 polarizabilities are in good relationships with some ITA quantities. The transferability of the linear regression model is further verified for a series of molecules with little or no similarity to those molecules in the training set. Furthermore, the excitation energies can be predicted based on multivariant linear regression equations of ITA quantities. This study also found that the nature of both the ground-state and excited-state polarizabilities of these species are due to the spatial delocalization of the electron density.
Collapse
Affiliation(s)
- Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, China
| | - Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao 266237, China
| | - Paul W. Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON L8S 4M1, Canada
- Correspondence: (P.W.A.); (S.L.)
| | - 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
- Correspondence: (P.W.A.); (S.L.)
| |
Collapse
|
48
|
Oozeer NB, Ponra A, Etindele AJ, Casida ME. A new freely-downloadable hands-on density-functional theory workbook using a freely-downloadable version of deMon2k. PURE APPL CHEM 2023. [DOI: 10.1515/pac-2022-1109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2023]
Abstract
Abstract
A hands-on workbook for density-functional theory (DFT) has been developed that can be used to provide practical teaching for students at the Masters or advanced undergraduate level that is free, can be used on a student’s own personal computer, and complements formal course work. The workbook is also very much intended to encourage students to explore program options, discover theory limitations, puzzle out what to do when the program does not work as expected, and to help students transition to thinking and using quantum chemistry programs as a researcher might do. After describing the structure of the workbook, we describe how the workbook has been used thus far as a teaching tool and as a useful step towards research-level problems.
Collapse
Affiliation(s)
- Nabila B. Oozeer
- Département de Chimie Moléculaire (DCM, UMR CNRS/UGA 5250), Laboratoire de Spectrométrie, Interactions et Chimie théorique (SITh), Institut de Chimie Moléculaire de Grenoble (ICMG, FR2607) , Université Grenoble Alpes (UGA) , 301 rue de la Chimie, BP 53 , F-38041 Grenoble Cedex 9 , France
| | | | - Anne Justine Etindele
- Higher Teachers Training College , University of Yaounde I , P.O. Box 47 , Yaounde , Cameroon
| | - Mark E. Casida
- Département de Chimie Moléculaire (DCM, UMR CNRS/UGA 5250), Laboratoire de Spectrométrie, Interactions et Chimie théorique (SITh), Institut de Chimie Moléculaire de Grenoble (ICMG, FR2607) , Université Grenoble Alpes (UGA) , 301 rue de la Chimie, BP 53 , F-38041 Grenoble Cedex 9 , France
| |
Collapse
|
49
|
Jha R, Kumar Chattaraj P. Effect of confinement on the structure, stability and aromaticity of Be32-. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2023.140390] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023]
|
50
|
Zhou K, Qian C, Liu Y. Quantifying the Structure of Water and Hydrated Monovalent Ions by Density Functional Theory-Based Molecular Dynamics. J Phys Chem B 2022; 126:10471-10480. [PMID: 36451081 DOI: 10.1021/acs.jpcb.2c05330] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
The accurate description of the structures of water and hydrated ions is important in electrochemical desalination, ion separation, and supercapacitors. In this work, we present an ab initio atomistic simulation-based study to explore the structure of water and hydrated monovalent ions (Li+, Na+, K+, Rb+, F-, and Cl-) at ambient conditions using generalized gradient approximation (GGA)-based methods with and without van der Waals correction (PBE, PBE + D3, and revPBE + D3) and recently developed strongly constrained and appropriately normed (SCAN) meta-GGA. We find that both revPBE + D3 and SCAN can well capture the structure of bulk water with +30 K artificial high temperature in contrast to overstructuring water using PBE and PBE + D3. However, being the same as PBE + D3, revPBE + D3 overestimates the structure of the hydration shell, especially for monovalent cations. Surprisingly, SCAN can well match the experimental results of hydrated monovalent ions. Detailed structure analyzes of entropy reveal that the hydration shell under the level of PBE + D3 and revPBE + D3 is more disordered and looser than SCAN. The successful prediction of the flexible SCAN functional could facilitate the exploration of complex ionic processes in the aqueous phase, the interactions of hydrated ions with surfaces, and solvation states in nanopores at an accurate, efficient, predictive, and ab initio level.
Collapse
Affiliation(s)
- Ke Zhou
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou215006, China.,Laboratory for Multiscale Mechanics and Medical Science, SV LAB, School of Aerospace, Xi'an Jiaotong University, Xi'an710049, China
| | - Chen Qian
- Department of Mechanical Engineering, Zhejiang University, Hangzhou310058, China
| | - Yilun Liu
- Laboratory for Multiscale Mechanics and Medical Science, SV LAB, School of Aerospace, Xi'an Jiaotong University, Xi'an710049, China
| |
Collapse
|