1
|
Han Y, Xu H, Li Q, Du A, Yan X. DFT-assisted low-dimensional carbon-based electrocatalysts design and mechanism study: a review. Front Chem 2023; 11:1286257. [PMID: 37920412 PMCID: PMC10619919 DOI: 10.3389/fchem.2023.1286257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Accepted: 09/27/2023] [Indexed: 11/04/2023] Open
Abstract
Low-dimensional carbon-based (LDC) materials have attracted extensive research attention in electrocatalysis because of their unique advantages such as structural diversity, low cost, and chemical tolerance. They have been widely used in a broad range of electrochemical reactions to relieve environmental pollution and energy crisis. Typical examples include hydrogen evolution reaction (HER), oxygen evolution reaction (OER), oxygen reduction reaction (ORR), carbon dioxide reduction reaction (CO2RR), and nitrogen reduction reaction (NRR). Traditional "trial and error" strategies greatly slowed down the rational design of electrocatalysts for these important applications. Recent studies show that the combination of density functional theory (DFT) calculations and experimental research is capable of accurately predicting the structures of electrocatalysts, thus revealing the catalytic mechanisms. Herein, current well-recognized collaboration methods of theory and practice are reviewed. The commonly used calculation methods and the basic functionals are briefly summarized. Special attention is paid to descriptors that are widely accepted as a bridge linking the structure and activity and the breakthroughs for high-volume accurate prediction of electrocatalysts. Importantly, correlated multiple descriptors are used to systematically describe the complicated interfacial electrocatalytic processes of LDC catalysts. Furthermore, machine learning and high-throughput simulations are crucial in assisting the discovery of new multiple descriptors and reaction mechanisms. This review will guide the further development of LDC electrocatalysts for extended applications from the aspect of DFT computations.
Collapse
Affiliation(s)
- Yun Han
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan Campus, Brisbane, QLD, Australia
- School of Engineering and Built Environment, Griffith University, Nathan Campus, Brisbane, QLD, Australia
| | - Hongzhe Xu
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan Campus, Brisbane, QLD, Australia
- School of Engineering and Built Environment, Griffith University, Nathan Campus, Brisbane, QLD, Australia
| | - Qin Li
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan Campus, Brisbane, QLD, Australia
- School of Engineering and Built Environment, Griffith University, Nathan Campus, Brisbane, QLD, Australia
| | - Aijun Du
- School of Chemistry and Physics and Centre for Materials Science, Queensland University of Technology, Gardens Point Campus, Brisbane, QLD, Australia
| | - Xuecheng Yan
- Queensland Micro- and Nanotechnology Centre, Griffith University, Nathan Campus, Brisbane, QLD, Australia
| |
Collapse
|
2
|
Brémond E, Rodríguez-Mayorga M, Pérez-Jiménez AJ, Adamo C, Sancho-García JC. Assessment of the nonempirical r2SCAN-QIDH double-hybrid density functional against large and diverse datasets. J Chem Phys 2023; 159:141101. [PMID: 37811824 DOI: 10.1063/5.0157259] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2023] [Accepted: 09/21/2023] [Indexed: 10/10/2023] Open
Abstract
We update the Quadratic Integrand Double-Hybrid (QIDH) model [J. Chem. Phys. 141, 031101 (2014)] by incorporating the nonempirical restored-regularized Strongly Constrained and Appropriately Normed (r2SCAN) meta-generalized gradient approximation exchange-correlation functional, thus devising a robust density functional approximation free of any empirical parameter and incorporating all the constraints so far known for the exchange-correlation kernel. We assessed the new r2SCAN-QIDH expression on the GMTKN55 database and further extend its application to various types of non-covalent interactions (e.g., S66 × 8, O24 × 5). The assessment done shows that the model becomes very competitive in accuracy with respect to parent exchange-correlation functionals of any type, but without relying on any fitted parameter or numerical training.
Collapse
Affiliation(s)
- E Brémond
- Université de Paris Cité, ITODYS, CNRS, F-75006 Paris, France
| | - M Rodríguez-Mayorga
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - A J Pérez-Jiménez
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - C Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), FRE 2027, F-75005 Paris, France
- Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| | - J C Sancho-García
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| |
Collapse
|
3
|
Olatomiwa A, Adam T, Edet C, Adewale A, Chik A, Mohammed M, Gopinath SC, Hashim U. Recent advances in density functional theory approach for optoelectronics properties of graphene. Heliyon 2023; 9:e14279. [PMID: 36950613 PMCID: PMC10025043 DOI: 10.1016/j.heliyon.2023.e14279] [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: 12/08/2022] [Revised: 02/28/2023] [Accepted: 03/01/2023] [Indexed: 03/09/2023] Open
Abstract
Graphene has received tremendous attention among diverse 2D materials because of its remarkable properties. Its emergence over the last two decades gave a new and distinct dynamic to the study of materials, with several research projects focusing on exploiting its intrinsic properties for optoelectronic devices. This review provides a comprehensive overview of several published articles based on density functional theory and recently introduced machine learning approaches applied to study the electronic and optical properties of graphene. A comprehensive catalogue of the bond lengths, band gaps, and formation energies of various doped graphene systems that determine thermodynamic stability was reported in the literature. In these studies, the peculiarity of the obtained results reported is consequent on the nature and type of the dopants, the choice of the XC functionals, the basis set, and the wrong input parameters. The different density functional theory models, as well as the strengths and uncertainties of the ML potentials employed in the machine learning approach to enhance the prediction models for graphene, were elucidated. Lastly, the thermal properties, modelling of graphene heterostructures, the superconducting behaviour of graphene, and optimization of the DFT models are grey areas that future studies should explore in enhancing its unique potential. Therefore, the identified future trends and knowledge gaps have a prospect in both academia and industry to design future and reliable optoelectronic devices.
Collapse
Affiliation(s)
- A.L. Olatomiwa
- Institute of Nano Electronic Engineering, Universiti Malaysia Perlis, 01000, Kangar, Perlis, Malaysia
- Faculty of Electronic Engineering and Technology, Universiti Malaysia Perlis, 02600, Arau, Perlis, Malaysia
| | - Tijjani Adam
- Institute of Nano Electronic Engineering, Universiti Malaysia Perlis, 01000, Kangar, Perlis, Malaysia
- Faculty of Electronic Engineering and Technology, Universiti Malaysia Perlis, 02600, Arau, Perlis, Malaysia
- Micro System Technology, Centre of Excellence (CoE), Universiti Malaysia Perlis (UniMAP), Pauh Campus, 02600, Arau, Perlis, Malaysia
| | - C.O. Edet
- Faculty of Electronic Engineering and Technology, Universiti Malaysia Perlis, 02600, Arau, Perlis, Malaysia
- Institute of Engineering Mathematics, Universiti Malaysia Perlis, 02600, Arau, Perlis, Malaysia
- Department of Physics, Cross River University of Technology, Calabar, Nigeria
| | - A.A. Adewale
- Department of Pure and Applied Physics, Ladoke Akintola University of Technology, Ogbomoso, Nigeria
| | - Abdullah Chik
- Centre for Frontier Materials Research, Universiti Malaysia Perlis, 01000, Kangar, Perlis, Malaysia
- Faculty of Chemical Engineering and Technology, Universiti Malaysia Perlis (UniMAP), Taman Muhibbah, Jejawi, 02600, Arau, Perlis, Malaysia
| | - Mohammed Mohammed
- Faculty of Chemical Engineering and Technology, Universiti Malaysia Perlis (UniMAP), Taman Muhibbah, Jejawi, 02600, Arau, Perlis, Malaysia
- Center of Excellence Geopolymer & Green Technology (CEGeoGTech), Universiti Malaysia Perlis, 02600, Arau, Perlis, Malaysia
| | - Subash C.B. Gopinath
- Institute of Nano Electronic Engineering, Universiti Malaysia Perlis, 01000, Kangar, Perlis, Malaysia
- Micro System Technology, Centre of Excellence (CoE), Universiti Malaysia Perlis (UniMAP), Pauh Campus, 02600, Arau, Perlis, Malaysia
- Faculty of Chemical Engineering and Technology, Universiti Malaysia Perlis (UniMAP), Taman Muhibbah, Jejawi, 02600, Arau, Perlis, Malaysia
| | - U. Hashim
- Institute of Nano Electronic Engineering, Universiti Malaysia Perlis, 01000, Kangar, Perlis, Malaysia
| |
Collapse
|
4
|
Hersh WH, Chan TY. Improving the accuracy of 31P NMR chemical shift calculations by use of scaling methods. Beilstein J Org Chem 2023; 19:36-56. [PMID: 36726479 PMCID: PMC9843238 DOI: 10.3762/bjoc.19.4] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2022] [Accepted: 12/14/2022] [Indexed: 01/12/2023] Open
Abstract
Calculation of 31P NMR chemical shifts for a series of tri- and tetracoordinate phosphorus compounds using several basis sets and density functional theory (DFT) functionals gave a modest fit to experimental chemical shifts, but an excellent linear fit when plotted against the experimental values. The resultant scaling methods were then applied to a variety of "large" compounds previously selected by Latypov et al. and a set of stereoisomeric and unusual compounds selected here. No one method was best for all structural types. For compounds that contain P-P bonds and P-C multiple bonds, the Latypov et al. method using the PBE0 functional was best (mean absolute deviation/root mean square deviation (MAD/RMSD) = 6.9/8.5 ppm and 6.6/8.2 ppm, respectively), but for the full set of compounds gave higher deviations (MAD/RMSD = 8.2/12.3 ppm), and failed by over 60 ppm for a three-membered phosphorus heterocycle. Use of the M06-2X functional for both the structural optimization and NMR chemical shift calculation was best overall for the compounds without P-C multiple bonds (MAD/RMSD = 5.4/7.1 ppm), but failed by 30-49 ppm for compounds having any P-C multiple-bond character. Failures of these magnitudes have not been reported previously for these widely used functionals. These failures were then used to screen a variety of recommended functionals, leading to better overall methods for calculation of these chemical shifts: optimization with the M06-2X functional and NMR calculation with the PBE0 or ωB97x-D functionals gave values for MAD/RMSD = 6.9/8.5 ppm and 6.8/9.1 ppm, respectively, over an experimental chemical shift range of -181 to 356 ppm. Due to the unexplained failures observed, we recommend use of more than one method when looking at novel structures.
Collapse
Affiliation(s)
- William H Hersh
- Department of Chemistry and Biochemistry, Queens College, Queens, NY 11367-1597, USA,Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, USA
| | - Tsz-Yeung Chan
- Department of Chemistry and Biochemistry, Queens College, Queens, NY 11367-1597, USA,Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, USA
| |
Collapse
|
5
|
Kraka E, Quintano M, La Force HW, Antonio JJ, Freindorf M. The Local Vibrational Mode Theory and Its Place in the Vibrational Spectroscopy Arena. J Phys Chem A 2022; 126:8781-8798. [DOI: 10.1021/acs.jpca.2c05962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| | - Mateus Quintano
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| | - Hunter W. La Force
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| | - Juliana J. Antonio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| | - Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, Texas75275-0314, United States
| |
Collapse
|
6
|
Tan S, Wang Y, Zhang IY, Xu X. Assessment of advanced xDH@B3LYP methods in describing various potential energy curves driven by π- π, CH/ π, and SH/ π non-bonded interactions. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2206098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Accurate description of potential energy curves driven by non-bonded interactions remains a great challenge for pure density functional approximations (DFAs). It is because the R−6 decay behavior of dispersion cannot be intrinsically captured by the (semi)-local ingredients and the exact-exchange used in the popular hybrid DFAs. Overemphasizing the accuracy on the equilibrium region for the functional construction would likely deteriorate the overall performance on the other regions of potential energy surfaces. In consequence, the empirical dispersion correction becomes the standard component in DFAs to treat the non-bonded interactions. In this Letter, we demonstrate that without the use of empirical dispersion correction, doubly hybrid approximations, in particular two recently proposed revXYG3 and XYG7 functionals, hold the promise to have a balanced description of non-bonded interactions on the whole potential energy curves for several prototypes of π- π, CH/ π, and SH/ π interactions. The error of revXYG3 and XYG7 for non-bonded interactions is around 0.1 kcal/mol, and their potential energy curves almost coincide with the accurate CCSD(T)/CBS curves.
Collapse
Affiliation(s)
- Shiqian Tan
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Yizhen Wang
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Igor Ying Zhang
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, China
| |
Collapse
|
7
|
Śmiga S, Della Sala F, Gori-Giorgi P, Fabiano E. Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit. J Chem Theory Comput 2022; 18:5936-5947. [PMID: 36094908 PMCID: PMC9558377 DOI: 10.1021/acs.jctc.2c00352] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Adiabatic connection
models (ACMs), which interpolate between the
limits of weak and strong interaction, are powerful tools to build
accurate exchange–correlation functionals. If the exact weak-interaction
expansion from the second-order perturbation theory is included, a
self-consistent implementation of these functionals is challenging
and still absent in the literature. In this work, we fill this gap
by presenting a fully self-consistent-field (SCF) implementation of
some popular ACM functionals. While using second-order perturbation
theory at weak interactions, we have also introduced new generalized
gradient approximations (GGAs), beyond the usual point-charge-plus-continuum
model, for the first two leading terms at strong interactions, which
are crucial to ensure robustness and reliability. We then assess the
SCF–ACM functionals for molecular systems and for prototypical
strong-correlation problems. We find that they perform well for both
the total energy and the electronic density and that the impact of
SCF orbitals is directly connected to the accuracy of the ACM functional
form. For the H2 dissociation, the SCF–ACM functionals
yield significant improvements with respect to standard functionals
also thanks to the use of the new GGAs for the strong-coupling functionals.
Collapse
Affiliation(s)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudziądzka 5, 87-100 Toruń, Poland
| | - Fabio Della Sala
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Via Monteroni 73100, Italy.,Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano, Lecce 73010, Italy
| | - 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, 1081HV Amsterdam, The Netherlands
| | - Eduardo Fabiano
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Via Monteroni 73100, Italy.,Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano, Lecce 73010, Italy
| |
Collapse
|
8
|
Fu Z, Liu W, Huang C, Mei T. A Review of Performance Prediction Based on Machine Learning in Materials Science. NANOMATERIALS (BASEL, SWITZERLAND) 2022; 12:nano12172957. [PMID: 36079994 PMCID: PMC9457802 DOI: 10.3390/nano12172957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Revised: 08/07/2022] [Accepted: 08/24/2022] [Indexed: 05/11/2023]
Abstract
With increasing demand in many areas, materials are constantly evolving. However, they still have numerous practical constraints. The rational design and discovery of new materials can create a huge technological and social impact. However, such rational design and discovery require a holistic, multi-stage design process, including the design of the material composition, material structure, material properties as well as process design and engineering. Such a complex exploration using traditional scientific methods is not only blind but also a huge waste of time and resources. Machine learning (ML), which is used across data to find correlations in material properties and understand the chemical properties of materials, is being considered a new way to explore the materials field. This paper reviews some of the major recent advances and applications of ML in the field of properties prediction of materials and discusses the key challenges and opportunities in this cross-cutting area.
Collapse
Affiliation(s)
- Ziyang Fu
- School of Computer Science and Information Engineering, Hubei University, Wuhan 430062, China
- Hubei Software Engineering Technology Research Center, Wuhan 430062, China
- Hubei Engineering Research Center for Smart Government and Artificial Intelligence Application, Wuhan 430062, China
| | - Weiyi Liu
- School of Materials Science and Engineering, Hubei University, Wuhan 430062, China
| | - Chen Huang
- School of Computer Science and Information Engineering, Hubei University, Wuhan 430062, China
- Hubei Software Engineering Technology Research Center, Wuhan 430062, China
- Hubei Engineering Research Center for Smart Government and Artificial Intelligence Application, Wuhan 430062, China
- Correspondence: (C.H.); (T.M.)
| | - Tao Mei
- School of Materials Science and Engineering, Hubei University, Wuhan 430062, China
- Hubei Collaborative Innovation Center for Advanced Organic Chemical Materials, Wuhan 430062, China
- Key Laboratory for the Green Preparation and Application of Functional Materials, Wuhan 430062, China
- Correspondence: (C.H.); (T.M.)
| |
Collapse
|
9
|
Xie W, Xu J, Chen J, Wang H, Hu P. Achieving Theory-Experiment Parity for Activity and Selectivity in Heterogeneous Catalysis Using Microkinetic Modeling. Acc Chem Res 2022; 55:1237-1248. [PMID: 35442027 PMCID: PMC9069691 DOI: 10.1021/acs.accounts.2c00058] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Microkinetic modeling based on density functional
theory (DFT)
energies plays an essential role in heterogeneous catalysis because
it reveals the fundamental chemistry for catalytic reactions and bridges
the microscopic understanding from theoretical calculations and experimental
observations. Microkinetic modeling requires building a set of ordinary
differential equations (ODEs) based on the calculation results of
thermodynamic properties of adsorbates and kinetic parameters for
the reaction elementary steps. Solving a microkinetic model can extract
information on catalytic chemistry, including critical reaction intermediates,
reaction pathways, the surface species distribution, activity, and
selectivity, thus providing vital guidelines for altering catalysts. However, the quantitative reliability of traditional microkinetic
models is often insufficient to conclusively extrapolate the mechanistic
details of complex reaction systems. This can be attributed to several
factors, the most important of which is the limitation of obtaining
an accurate estimation of the energy inputs via traditional calculation
methods. These limitations include the difficulty of using static
DFT methods to calculate reaction energies of adsorption/desorption
processes, often rate-controlling or selectivity-determining steps,
and the inadequate consideration of surface coverage effects. In addition,
the robust microkinetic software is rare, which also complicates the
resolution of complex catalytic systems. In this Account, we
review our recent works toward refining the
predictions of microkinetic modeling in heterogeneous catalysis and
achieving theory–experiment parity for activity and selectivity.
First, we introduce CATKINAS, a microkinetic software developed in
our group, and show how it disentangles the problem that traditional
microkinetic software has and how it can now be applied to obtain
kinetic results for more sophisticated reaction systems. Second, we
describe a molecular dynamics method developed recently to obtain
the free-energy changes for the adsorption/desorption process to fill
in the missing energy inputs. Third, we show that a rigorous consideration
of surface coverage effects is pivotal for building more realistic
models and obtaining accurate kinetic results. Following a series
of studies on acetylene hydrogenation reactions on Pd catalysts, we
demonstrate how this new approach can provide an improved quantitative
understanding of the mechanism, active site, and intrinsic structural
sensitivity. Finally, we conclude with a brief outlook and the remaining
challenges in this field.
Collapse
Affiliation(s)
- Wenbo Xie
- School of Chemistry and Chemical Engineering, The Queen’s University of Belfast, Belfast BT9 5AG, U.K
| | - Jiayan Xu
- School of Chemistry and Chemical Engineering, The Queen’s University of Belfast, Belfast BT9 5AG, U.K
| | - Jianfu Chen
- Key Laboratory for Advanced Materials, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, East China University of Science and Technology, Shanghai 200237, P. R. China
| | - Haifeng Wang
- Key Laboratory for Advanced Materials, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, East China University of Science and Technology, Shanghai 200237, P. R. China
| | - P. Hu
- School of Chemistry and Chemical Engineering, The Queen’s University of Belfast, Belfast BT9 5AG, U.K
- Key Laboratory for Advanced Materials, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, East China University of Science and Technology, Shanghai 200237, P. R. China
| |
Collapse
|
10
|
Yeh SH, Yang W, Hsu CP. Reformulation of Thermally-Assisted-Occupation Density Functional Theory in the Kohn-Sham Framework. J Chem Phys 2022; 156:174108. [DOI: 10.1063/5.0087012] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We reformulate the thermally-assisted-occupation density functional theory (TAO-DFT) into the Kohn-Sham single-determinant framework and construct two new post-self-consistent field (post-SCF) static correlation correction schemes, named rTAO and rTAO-1. In contrast to the original TAO-DFT with the density in an ensemble form, in which each orbital density is weighted with a fractional occupation number, the ground-state density is given by a single-determinant wavefunction, a regular Kohn-Sham density, and total ground state energy is expressed in the normal Kohn-Sham form with a static correlation energy formulated in terms of the Kohn-Sham orbitals. In post-SCF calculations with rTAO functionals, an efficient energy scanning to quantitatively determine $\theta$ is also proposed. The rTAOs provide a promising method to simulate systems with strong static correlation as original TAO, but simpler and more efficient.We show that both rTAO and rTAO-1 is capable of reproducing most results from TAO-DFT without the additional functional Eθ used in TAO-DFT. Furthermore, our numerical results support that, without the functional Eθ, both rTAO and rTAO-1 can capture correct static correlation profiles in various systems.
Collapse
Affiliation(s)
- Shu-Hao Yeh
- Institute of Chemistry Academia Sinica, Taiwan
| | - Weitao Yang
- Department of Chemistry, Duke University, United States of America
| | - Chao-Ping Hsu
- Institute of Chemistry, Institute of Chemistry Academia Sinica, Taiwan
| |
Collapse
|
11
|
Prokopiou G, Hartstein M, Govind N, Kronik L. Optimal Tuning Perspective of Range-Separated Double Hybrid Functionals. J Chem Theory Comput 2022; 18:2331-2340. [PMID: 35369687 PMCID: PMC9009176 DOI: 10.1021/acs.jctc.2c00082] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2022] [Indexed: 11/29/2022]
Abstract
We study the optimal tuning of the free parameters in range-separated double hybrid functionals, based on enforcing the exact conditions of piecewise linearity and spin constancy. We find that introducing the range separation in both the exchange and the correlation terms allows for the minimization of both fractional charge and fractional spin errors for singlet atoms. The optimal set of parameters is system specific, underlining the importance of the tuning procedure. We test the performance of the resulting optimally tuned functionals for the dissociation curves of diatomic molecules. We find that they recover the correct dissociation curve for the one-electron system, H2+, and improve the dissociation curves of many-electron molecules such as H2 and Li2, but they also yield a nonphysical maximum and only converge to the correct dissociation limit at very large distances.
Collapse
Affiliation(s)
- Georgia Prokopiou
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Michal Hartstein
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Niranjan Govind
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Leeor Kronik
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| |
Collapse
|
12
|
Brémond É, Tognetti V, Chermette H, Sancho-García JC, Joubert L, Adamo C. Electronic Energy and Local Property Errors at QTAIM Critical Points while Climbing Perdew's Ladder of Density-Functional Approximations. J Chem Theory Comput 2021; 18:293-308. [PMID: 34958205 DOI: 10.1021/acs.jctc.1c00981] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We investigate the relationships between electron-density and electronic-energy errors produced by modern exchange-correlation density-functional approximations belonging to all of the rungs of Perdew's ladder. To this aim, a panel of relevant (semi)local properties evaluated at critical points of the electron-density field (as defined within the framework of Bader's atoms-in-molecules theory) are computed on a large selection of molecular systems involved in thermodynamic, kinetic, and noncovalent interaction chemical databases using density functionals developed in a nonempirical and minimally and highly parametrized fashion. The comparison of their density- and energy-based performance, also discussed in terms of density-driven errors, casts light on the strengths and weaknesses of the most recent and efficient density-functional approximations.
Collapse
Affiliation(s)
- Éric Brémond
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France
| | - Vincent Tognetti
- Normandy University, COBRA UMR 6014 and FR 3038, Université de Rouen INSA Rouen, CNRS, F-76821 Mont St Aignan, France
| | - Henry Chermette
- Université de Lyon, Institut des Sciences Analytiques, UMR 5280, CNRS, Université Lyon 1, 5 rue de la Doua, F-69100 Villeurbanne, France
| | | | - Laurent Joubert
- Normandy University, COBRA UMR 6014 and FR 3038, Université de Rouen INSA Rouen, CNRS, F-76821 Mont St Aignan, France
| | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), UMR 8060, F-75005 Paris, France.,Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| |
Collapse
|
13
|
Grofe A, Gao J, Li X. Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF. J Chem Phys 2021; 155:014103. [PMID: 34241404 DOI: 10.1063/5.0054227] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to a relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time-reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break the symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time-reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.
Collapse
Affiliation(s)
- Adam Grofe
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China; Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA; and Beijing University Shenzhen Graduate School, Shenzhen 518055, China
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
14
|
Wang Y, Li Y, Chen J, Zhang IY, Xu X. Doubly Hybrid Functionals Close to Chemical Accuracy for Both Finite and Extended Systems: Implementation and Test of XYG3 and XYGJ-OS. JACS AU 2021; 1:543-549. [PMID: 34467317 PMCID: PMC8395692 DOI: 10.1021/jacsau.1c00011] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
While being widely used to understand the chemical reactions in heterogeneous catalysis or other multidisciplinary systems, a great challenge that semilocal and hybrid density functional approximations (DFAs) are facing is to deliver a uniformly accurate description for both finite and extended systems. Herein, we perform reliable and well-converged periodic calculations of two doubly hybrid approximations (DHAs), XYG3 and XYGJ-OS, and demonstrate that the good accuracy of DHAs achieved for molecules is transferable to the semiconductors and insulators. Such an accuracy is not only for energetic properties but also for the first- and second-order response properties, which is general for different kinds of chemical environments, including simple cubic bulks, perovskite-type transition metal oxides like TiO2, and heterogeneous systems like CO adsorption on the NaCl(100) surface. The present finding has strengthened the predictive power of DFT, which not only will inspire the future development of the top-rung DFAs but also will boost their applications in multidisciplinary studies with high accuracy and efficiency.
Collapse
Affiliation(s)
- Yizhen Wang
- Collaborative
Innovation Center of Chemistry for Energy Materials, Shanghai, Key
Laboratory of Molecular Catalysis and Innovative Materials, MOE Key
Laboratory of Computational Physical Sciences, Shanghai Key Laboratory
of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, People’s Republic of China
| | - Yajing Li
- Collaborative
Innovation Center of Chemistry for Energy Materials, Shanghai, Key
Laboratory of Molecular Catalysis and Innovative Materials, MOE Key
Laboratory of Computational Physical Sciences, Shanghai Key Laboratory
of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, People’s Republic of China
| | - Jun Chen
- State
Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese
Academy of Sciences, Fuzhou 350002, People’s Republic
of China
| | - Igor Ying Zhang
- Collaborative
Innovation Center of Chemistry for Energy Materials, Shanghai, Key
Laboratory of Molecular Catalysis and Innovative Materials, MOE Key
Laboratory of Computational Physical Sciences, Shanghai Key Laboratory
of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, People’s Republic of China
| | - Xin Xu
- Collaborative
Innovation Center of Chemistry for Energy Materials, Shanghai, Key
Laboratory of Molecular Catalysis and Innovative Materials, MOE Key
Laboratory of Computational Physical Sciences, Shanghai Key Laboratory
of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, People’s Republic of China
| |
Collapse
|
15
|
Cao X, Tian P. "Dividing and Conquering" and "Caching" in Molecular Modeling. Int J Mol Sci 2021; 22:5053. [PMID: 34068835 PMCID: PMC8126232 DOI: 10.3390/ijms22095053] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2021] [Revised: 04/26/2021] [Accepted: 04/27/2021] [Indexed: 11/17/2022] Open
Abstract
Molecular modeling is widely utilized in subjects including but not limited to physics, chemistry, biology, materials science and engineering. Impressive progress has been made in development of theories, algorithms and software packages. To divide and conquer, and to cache intermediate results have been long standing principles in development of algorithms. Not surprisingly, most important methodological advancements in more than half century of molecular modeling are various implementations of these two fundamental principles. In the mainstream classical computational molecular science, tremendous efforts have been invested on two lines of algorithm development. The first is coarse graining, which is to represent multiple basic particles in higher resolution modeling as a single larger and softer particle in lower resolution counterpart, with resulting force fields of partial transferability at the expense of some information loss. The second is enhanced sampling, which realizes "dividing and conquering" and/or "caching" in configurational space with focus either on reaction coordinates and collective variables as in metadynamics and related algorithms, or on the transition matrix and state discretization as in Markov state models. For this line of algorithms, spatial resolution is maintained but results are not transferable. Deep learning has been utilized to realize more efficient and accurate ways of "dividing and conquering" and "caching" along these two lines of algorithmic research. We proposed and demonstrated the local free energy landscape approach, a new framework for classical computational molecular science. This framework is based on a third class of algorithm that facilitates molecular modeling through partially transferable in resolution "caching" of distributions for local clusters of molecular degrees of freedom. Differences, connections and potential interactions among these three algorithmic directions are discussed, with the hope to stimulate development of more elegant, efficient and reliable formulations and algorithms for "dividing and conquering" and "caching" in complex molecular systems.
Collapse
Affiliation(s)
- Xiaoyong Cao
- School of Life Sciences, Jilin University, Changchun 130012, China;
| | - Pu Tian
- School of Life Sciences, Jilin University, Changchun 130012, China;
- School of Artificial Intelligence, Jilin University, Changchun 130012, China
| |
Collapse
|
16
|
Zhang IY, Xu X. Exploring the Limits of the XYG3-Type Doubly Hybrid Approximations for the Main-Group Chemistry: The xDH@B3LYP Model. J Phys Chem Lett 2021; 12:2638-2644. [PMID: 33689361 DOI: 10.1021/acs.jpclett.1c00360] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Despite being the most accurate class of density functional approximations for the main-group chemistry, doubly hybrid approximations (DHAs) are generally considered to be incomplete in describing the medium- to long-range dispersive interactions. The existing DHAs are often supplemented with empirical long-range dispersion corrections. By using the extensive and chemically diverse GMTKN55 database, we explore the limits of the XYG3-type DHAs using the B3LYP reference orbitals, namely, xDH@B3LYP, with a gradually relaxed constraint on the mixing parameters of DHAs. Our results demonstrate that the xDH@B3LYP model can provide a balanced description of both covalent and noncovalent interactions with the accuracy and robustness comparable to or even better than the very expensive composite methods in wave function theory. Such an accuracy can be achieved without resorting to the use of any long-range correction scheme, shedding new light on the development of DHAs.
Collapse
Affiliation(s)
- Igor Ying Zhang
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, People's Republic of China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, People's Republic of China
| |
Collapse
|
17
|
Galushko AS, Prima DO, Burykina JV, Ananikov VP. Comparative study of aryl halides in Pd-mediated reactions: key factors beyond the oxidative addition step. Inorg Chem Front 2021. [DOI: 10.1039/d0qi01133a] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The comparative experimental study of Ar–X (X = Cl, Br, I) reactivity and analysis reported herein suggest that oxidative addition cannot be considered the sole reason of the observed low reactivity of aryl chlorides.
Collapse
Affiliation(s)
- Alexey S. Galushko
- Zelinsky Institute of Organic Chemistry
- Russian Academy of Sciences
- Moscow
- Russia
| | - Darya O. Prima
- Zelinsky Institute of Organic Chemistry
- Russian Academy of Sciences
- Moscow
- Russia
| | - Julia V. Burykina
- Zelinsky Institute of Organic Chemistry
- Russian Academy of Sciences
- Moscow
- Russia
| | | |
Collapse
|