1
|
Marzouk A, Papavasileiou KD, Peristeras LD, Bezemer L, van Bavel AP, Shenai PM, Economou IG. A systematic DFT study of structure and electronic properties of titanium dioxide. J Comput Chem 2024; 45:2153-2166. [PMID: 38785277 DOI: 10.1002/jcc.27376] [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/28/2024] [Revised: 04/02/2024] [Accepted: 04/08/2024] [Indexed: 05/25/2024]
Abstract
DFT functionals are of paramount importance for an accurate electronic and structural description of transition metal systems. In this work, a systematic analysis using some well-known and commonly used DFT functionals is performed. A comparison of the structural and energetic parameters calculated with the available experimental data is made in order to find the adequate functional for an accurate description of the TiO2 bulk and surface of both anatase and rutile structures. In the absence of experimental data on the surface energy, the theoretical predictions obtained using the high-accuracy HSE06 functional were used as a reference to compare against the surface energy values calculated with the other DFT functionals. A clear improvement in the electronic description of both anatase and rutile was observed by introducing the Hubbard U correction term to PBE, PW91, and OptPBE functionals. The OptPBE-U4 functional was found to offer a good compromise between accurately describing the structural and electronic properties of titania.
Collapse
Affiliation(s)
- Asma Marzouk
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| | - Konstantinos D Papavasileiou
- Molecular Thermodynamics and Modelling of Materials Laboratory, National Center for Scientific Research "Demokritos", Institute of Nanoscience and Nanotechnology, Athens, Greece
| | - Loukas D Peristeras
- Molecular Thermodynamics and Modelling of Materials Laboratory, National Center for Scientific Research "Demokritos", Institute of Nanoscience and Nanotechnology, Athens, Greece
| | - Leendert Bezemer
- GTL and XTL Research, Shell Global Solutions International BV, Amsterdam, The Netherlands
| | - Alexander P van Bavel
- Next Generation Breakthrough Research, Shell Global Solutions International BV, Amsterdam, The Netherlands
| | - Prathamesh M Shenai
- Computational Chemistry and Material Science, Shell India Markets Pvt. Ltd, Shell India Markets Pvt. Ltd, Banglore, India
| | - Ioannis G Economou
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| |
Collapse
|
2
|
Chufarov M, Vassilyeva YZ, Zhang X, Li S, Pak AY, Han W. Design of a long-lived Mo2C-MoO2@GC-N electrocatalyst by the ambient DC arc plasma for the hydrogen evolution reaction. iScience 2024; 27:110551. [PMID: 39228789 PMCID: PMC11369373 DOI: 10.1016/j.isci.2024.110551] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2024] [Revised: 05/02/2024] [Accepted: 07/16/2024] [Indexed: 09/05/2024] Open
Abstract
A crucial challenge in hydrogen production through electrolysis is developing inexpensive, earth-abundant, and highly efficient Pt-free electrocatalysts for the hydrogen evolution reaction (HER). Molybdenum carbide is ideal for this application because of its special electrical structure, low cost, and advantageous characteristics. Herein, the long-lived electrocatalysts for HER have been synthesized via the direct current (DC) arc discharge plasma method under ambient air conditions, and the relationship between the properties of materials and catalytic characteristics has been established. The samples differed in the ratio of molybdenum, graphite, and melamine. The sample with the highest proportion of melamine in the initial mixture has Mo2C-MoO2 heterointerfaces, which demonstrates the highest and most stable electrocatalytic activity with the overpotential of 148 mV at 10 mA·cm-2 and Tafel slope of 63 mV·dec-1 in alkaline electrolyte. Meanwhile, the electrodes demonstrated long-lived electrochemical durability for two weeks and investigated the features of forming a stable system for HER.
Collapse
Affiliation(s)
- Marian Chufarov
- College of Physics, the State Key Laboratory of Inorganic Synthesis and Preparative Chemistry, International Center of Future Science, Jilin University, Changchun 130012, China
| | - Yuliya Z. Vassilyeva
- College of Physics, the State Key Laboratory of Inorganic Synthesis and Preparative Chemistry, International Center of Future Science, Jilin University, Changchun 130012, China
- Laboratory of Advanced Materials for Energy Industry, Tomsk Polytechnic University, Tomsk 634050, Russian Federation
| | - Xinyu Zhang
- Department of Physics, Southern University of Science and Technology, Shenzhen 518055, China
| | - Shilin Li
- College of Physics, the State Key Laboratory of Inorganic Synthesis and Preparative Chemistry, International Center of Future Science, Jilin University, Changchun 130012, China
| | - Alexander Y. Pak
- Laboratory of Advanced Materials for Energy Industry, Tomsk Polytechnic University, Tomsk 634050, Russian Federation
| | - Wei Han
- College of Physics, the State Key Laboratory of Inorganic Synthesis and Preparative Chemistry, International Center of Future Science, Jilin University, Changchun 130012, China
| |
Collapse
|
3
|
Macke E, Timrov I, Marzari N, Ciacchi LC. Orbital-Resolved DFT +U for Molecules and Solids. J Chem Theory Comput 2024; 20:4824-4843. [PMID: 38820347 PMCID: PMC11171274 DOI: 10.1021/acs.jctc.3c01403] [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/22/2023] [Revised: 03/25/2024] [Accepted: 05/02/2024] [Indexed: 06/02/2024]
Abstract
We present an orbital-resolved extension of the Hubbard U correction to density-functional theory (DFT). Compared to the conventional shell-averaged approach, the prediction of energetic, electronic and structural properties is strongly improved, particularly for compounds characterized by both localized and hybridized states in the Hubbard manifold. The numerical values of all Hubbard parameters are readily obtained from linear-response calculations. The relevance of this more refined approach is showcased by its application to bulk solids pyrite (FeS2) and pyrolusite (β-MnO2), as well as to six Fe(II) molecular complexes. Our findings indicate that a careful definition of Hubbard manifolds is indispensable for extending the applicability of DFT+U beyond its current boundaries. The present orbital-resolved scheme aims to provide a computationally undemanding yet accurate tool for electronic structure calculations of charge-transfer insulators, transition-metal (TM) complexes and other compounds displaying significant orbital hybridization.
Collapse
Affiliation(s)
- Eric Macke
- Faculty
of Production Engineering, Bremen Center
for Computational Materials Science and MAPEX Center for Materials
and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, 28359 Bremen, Germany
| | - Iurii Timrov
- Theory
and Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique
Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Nicola Marzari
- Theory
and Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique
Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- University
of Bremen Excellence Chair, Bremen Center
for Computational Materials Science, Am Fallturm 1, 28359 Bremen, Germany
| | - Lucio Colombi Ciacchi
- Faculty
of Production Engineering, Bremen Center
for Computational Materials Science and MAPEX Center for Materials
and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, 28359 Bremen, Germany
| |
Collapse
|
4
|
Gallagher C, Siddiqui W, Arnold T, Cheng C, Su E, Zhao Q. Benchmarking a Molecular Flake Model on the Road to Programmable Graphene-Based Single-Atom Catalysts. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2024; 128:2876-2883. [PMID: 38414836 PMCID: PMC10895666 DOI: 10.1021/acs.jpcc.3c07681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Revised: 01/19/2024] [Accepted: 01/30/2024] [Indexed: 02/29/2024]
Abstract
Single-atom catalysts (SACs) of embedding an active metal in nitrogen-doped graphene are emergent catalytic materials in various applications. The rational design of efficient SACs necessitates an electronic and mechanistic understanding of those materials with reliable quantum mechanical simulations. Conventional computational methods of modeling SACs involve using an infinite slab model with periodic boundary condition, limiting to the selection of generalized gradient approximations as the exchange correlation (XC) functional within density functional theory (DFT). However, these DFT approximations suffer from electron self-interaction error and delocalization error, leading to errors in predicted charge-transfer energetics. An alternative strategy is using a molecular flake model, which carved out the important catalytic center by cleaving C-C bonds and employing a hydrogen capping scheme to saturate the innocent dangling bonds at the molecular boundary. By doing so, we can afford more accurate hybrid XC functionals, or even high-level correlated wavefunction theory, to study those materials. In this work, we compared the structural, electronic, and catalytic properties of SACs simulated using molecular flake models and periodic slab models with first-row transition metals as the active sites. Molecular flake models successfully reproduced structural properties, including both global distortion and local metal-coordination environment, as well as electronic properties, including spin magnetic moments and metal partial charges, for all transition metals studied. In addition, we calculated CO binding strength as a descriptor for electrochemical CO2 reduction reactivity and noted qualitatively similar trends between two models. Using the computationally efficient molecular flake models, we investigated the effect of tuning Hartree-Fock exchange in a global hybrid functional on the CO binding strength and observed system-dependent sensitivities. Overall, our calculations provide valuable insights into the development of accurate and efficient computational tools to simulate SACs.
Collapse
Affiliation(s)
- Colin Gallagher
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Wali Siddiqui
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Tyler Arnold
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Carmen Cheng
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Eric Su
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Qing Zhao
- Department of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| |
Collapse
|
5
|
Mohamed A, Shaban M, Kordy MGM, Al-Senani GM, Eissa MF, Hamdy H. Fabrication and characterization of NiCu/GO and NiCu/rGO nanocomposites for fuel cell application. RSC Adv 2024; 14:6776-6792. [PMID: 38405070 PMCID: PMC10884890 DOI: 10.1039/d3ra07822a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2023] [Accepted: 02/06/2024] [Indexed: 02/27/2024] Open
Abstract
This study investigated the electrochemical behavior of NiCu, NiCu/GO, and NiCu/rGO nanocomposites designed by combining a modified Hummers' method and hydrothermal technique. The prepared nanocomposites are tested as electrocatalysts in direct alcohol oxidation fuel cells (DAFCs) to identify the role of GO and rGO as catalyst supports for the enhancement of the NiCu composite performance. The production of the NiCu/GO and NiCu/rGO composites was demonstrated by FTIR spectroscopy, EDX, and SEM analyses. In DAFCs experiments, NiCu/rGO has better catalytic activity than pure NiCu and NiCu/GO composites, whereas the use of rGO and GO as supports enhances the performance of NiCu by 468.2% and 377.7% in methanol and 255.6% and 105.9% in ethanol, respectively. The higher performance was caused by the increased density of active dots and the combined electronic effects in the designed catalysts. The stabilities of the catalysts and charge carriers' dynamics are studied using chronoamperometry and electrochemical impedance spectroscopy.
Collapse
Affiliation(s)
- Aya Mohamed
- Nanophotonics and Applications Lab, Physics Department, Faculty of Science, Beni-Suef University 62514 Beni-Suef Egypt
| | - Mohamed Shaban
- Department of Physics, Faculty of Science, Islamic University of Madinah 42351 Madinah Saudi Arabia
| | - Mohamed G M Kordy
- Nanophotonics and Applications Lab, Physics Department, Faculty of Science, Beni-Suef University 62514 Beni-Suef Egypt
- Biochemistry Department, Faculty of Science, Beni-Suef University Beni-Suef 62521 Egypt
| | - Ghadah M Al-Senani
- Department of Chemistry, College of Science, Princess Nourah bint Abdulrahman University P.O. Box 84428 Riyadh 11671 Saudi Arabia
| | - M F Eissa
- Nanophotonics and Applications Lab, Physics Department, Faculty of Science, Beni-Suef University 62514 Beni-Suef Egypt
| | - Hany Hamdy
- Nanophotonics and Applications Lab, Physics Department, Faculty of Science, Beni-Suef University 62514 Beni-Suef Egypt
| |
Collapse
|
6
|
Yang Y, Li X, Liu G, Liu H, Shi Y, Ye C, Fang Z, Ye M, Shen J. Hierarchical Ohmic Contact Interface Engineering for Efficient Hydrazine-Assisted Hydrogen Evolution Reaction. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2024; 36:e2307979. [PMID: 37879754 DOI: 10.1002/adma.202307979] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Revised: 10/24/2023] [Indexed: 10/27/2023]
Abstract
Hydrazine oxidation reaction coupled with hydrogen evolution reaction (HER) is an effective strategy to achieve low energy water splitting for hydrogen production. In order to realize the application of hydrazine-assisted HER system, researchers have been focusing on the development of electrocatalysts with integrated dual active sites, while the performance under high current density is still unsatisfying. In this work, hierarchical Ohmic contact interface engineering is designed and used as a bridge between the NiMo and Ni2 P heterojunction toward industrial current density applications, with the charge transfer impedance greatly eliminated via such a pathway with low energy barrier. As a proof-of-concept, the importance of charge redistribution and energy barrier at the Ohmic contact interface is investigated by significantly reducing the voltage of overall hydrazine splitting (OHzS) at high current density. Intriguingly, the NiMo/Ni2 P hierarchical Ohmic contact heterojunction can drive current densities of 100 and 500 mA cm-2 with only 181 and 343 mV cell voltage in the OHzS electrolyzer with high electrocatalytic stability. The proposed hierarchical Ohmic contact interface engineering paves new avenue for hydrogen production with low energy consumption.
Collapse
Affiliation(s)
- Yifan Yang
- Institute of Special Materials and Technology, Fudan University, Shanghai, 200433, P. R. China
- Department of Materials Science, Fudan University, Shanghai, 200433, P. R. China
| | - Xuanyang Li
- Institute of Special Materials and Technology, Fudan University, Shanghai, 200433, P. R. China
- Department of Chemistry, Fudan University, Shanghai, 200433, P. R. China
| | - Guanglei Liu
- Institute of Special Materials and Technology, Fudan University, Shanghai, 200433, P. R. China
- Department of Materials Science, Fudan University, Shanghai, 200433, P. R. China
| | - Huixiang Liu
- Institute of Special Materials and Technology, Fudan University, Shanghai, 200433, P. R. China
- Department of Materials Science, Fudan University, Shanghai, 200433, P. R. China
| | - Yuehao Shi
- Institute of Special Materials and Technology, Fudan University, Shanghai, 200433, P. R. China
- Department of Materials Science, Fudan University, Shanghai, 200433, P. R. China
| | - Chuming Ye
- Institute of Special Materials and Technology, Fudan University, Shanghai, 200433, P. R. China
- Department of Materials Science, Fudan University, Shanghai, 200433, P. R. China
| | - Zhan Fang
- Institute of Special Materials and Technology, Fudan University, Shanghai, 200433, P. R. China
- Department of Materials Science, Fudan University, Shanghai, 200433, P. R. China
| | - Mingxin Ye
- Institute of Special Materials and Technology, Fudan University, Shanghai, 200433, P. R. China
| | - Jianfeng Shen
- Institute of Special Materials and Technology, Fudan University, Shanghai, 200433, P. R. China
| |
Collapse
|
7
|
Naveas N, Pulido R, Marini C, Gargiani P, Hernandez-Montelongo J, Brito I, Manso-Silván M. First-Principles Calculations of Magnetite (Fe 3O 4) above the Verwey Temperature by Using Self-Consistent DFT + U + V. J Chem Theory Comput 2023; 19:8610-8623. [PMID: 37974305 PMCID: PMC10720343 DOI: 10.1021/acs.jctc.3c00860] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2023] [Revised: 10/10/2023] [Accepted: 10/30/2023] [Indexed: 11/19/2023]
Abstract
In this report, we have used the DFT + U + V approach, an extension of the DFT + U approach that takes into account both on-site and intersite interactions, to simulate structural, magnetic, and electronic properties together with the Fe and O K-edge XAS spectra of Fe3O4 above the Verwey temperature (Tv). Moreover, we compared the simulated XAS spectra with experimental XAS data. We examined both orthogonalized and nonorthogonalized atomic orbital projectors and compared DFT + U + V to DFT, DFT + U, and HSE as a hybrid functional. It is noteworthy that, despite the widespread use of the same Hubbard U value for Feoct and Fetet at the DFT + U level in the literature, the HP code identified two distinct values for them using the Hubbard approaches (DFT + U and DFT + U + V). The resulting Hubbard U and V parameters are strongly dependent on the chosen orbital projectors. This study demonstrates how DFT + U + V can improve the structural, magnetic, and electronic properties of Fe3O4 compared to approximate DFT and DFT + U. In this context, DFT + U + V supports the half-metallic character of the bulk crystal Fe3O4 above Tv, since the Fermi level is found in the t2g band with a Feoct down-spin. Thus, the observations in the current study emphasize the significance of intersite interactions in the theoretical analysis of Fe3O4 above the Tv.
Collapse
Affiliation(s)
- Nelson Naveas
- Departamento
de Física Aplicada, Universidad Autónoma
de Madrid, 28049 Madrid, Spain
- Departamento
de Ingeniería Química y Procesos de Minerales, Universidad de Antofagasta, Avenida Angamos 601, 1270300 Antofagasta, Chile
- Instituto
Universitario de Ciencia de Materiales “Nicolás Cabrera”
(INC), Universidad Autónoma de Madrid, Campus de Cantoblanco, 28049 Madrid, Spain
| | - Ruth Pulido
- Instituto
Universitario de Ciencia de Materiales “Nicolás Cabrera”
(INC), Universidad Autónoma de Madrid, Campus de Cantoblanco, 28049 Madrid, Spain
- Departamento
de Química, Universidad de Antofagasta, Avenida Angamos 601, 1270300 Antofagasta, Chile
| | - Carlo Marini
- CELLS−ALBA
Synchrotron, 08290 Cerdanyola del Valles, Spain
| | | | | | - Ivan Brito
- Departamento
de Química, Universidad de Antofagasta, Avenida Angamos 601, 1270300 Antofagasta, Chile
| | - Miguel Manso-Silván
- Departamento
de Física Aplicada, Universidad Autónoma
de Madrid, 28049 Madrid, Spain
- Instituto
Universitario de Ciencia de Materiales “Nicolás Cabrera”
(INC), Universidad Autónoma de Madrid, Campus de Cantoblanco, 28049 Madrid, Spain
- Centro
de Microanálisis de Materiales, Universidad
Autónoma de Madrid, Campus de Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
8
|
Chen Z, Fan Q, Huang M, Cölfen H. The Structure, Preparation, Characterization, and Intercalation Mechanism of Layered Hydroxides Intercalated with Guest Anions. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2023; 19:e2300509. [PMID: 37271930 DOI: 10.1002/smll.202300509] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Revised: 05/06/2023] [Indexed: 06/06/2023]
Abstract
Since the intercalation of anions into layered hydroxides (LHs) has a great impact not only on their nucleation and growth but also on their structure, composition, and size, the intercalation chemistry of LHs has aroused the strong interest of researchers. However, the progress in the fundamental understanding of LHs intercalated with guest anions have not been paralleled by a concomitant development of the preparation and performance improvement of such materials. Considering the guidance of a timely in-depth review for scientists in this area, a systematic introduction about the development that is made on the above-mentioned issues is highly needed but yet missing so far. Herein, recent advances in understanding the chemical composition and structure of LHs intercalated with guest anions are systematically summarized. Meanwhile, typical and emerging bottom-up synthesis methods of LHs intercalated with anions are reviewed, and the potential impact of external reaction parameters on the intercalation of anions into LHs are discussed . Besides, different analytical characterization techniques employed in the examination of guest anion-intercalated LHs are deliberated upon. Finally, although progress is slow in exploring the intercalation mechanism, as many examples as possible are included in this review and inferred the possible intercalation mechanism.
Collapse
Affiliation(s)
- Zongkun Chen
- Physical Chemistry, University of Konstanz, Universitätsstraße 10, D-78457, Konstanz, Germany
| | - Qiqi Fan
- Physical Chemistry, University of Konstanz, Universitätsstraße 10, D-78457, Konstanz, Germany
| | - Minghua Huang
- School of Materials Science and Engineering, Ocean University of China, Qingdao, 266100, China
| | - Helmut Cölfen
- Physical Chemistry, University of Konstanz, Universitätsstraße 10, D-78457, Konstanz, Germany
| |
Collapse
|
9
|
Lindgren P, Kastlunger G, Peterson AA. Electrochemistry from the atomic scale, in the electronically grand-canonical ensemble. J Chem Phys 2022; 157:180902. [DOI: 10.1063/5.0123656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
The ability to simulate electrochemical reactions from first-principles has advanced significantly in recent years. Here, we discuss the atomistic interpretation of electrochemistry at three scales: from the electronic structure to elementary processes to constant-potential reactions. At each scale, we highlight the importance of the grand-canonical nature of the process and show that the grand-canonical energy is the natural thermodynamic state variable, which has the additional benefit of simplifying calculations. We show that atomic forces are the derivative of the grand-potential energy when the potential is fixed. We further examine the meaning of potential at the atomic scale and its link to the chemical potential and discuss the link between charge transfer and potential in several situations.
Collapse
Affiliation(s)
- Per Lindgren
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - Georg Kastlunger
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
- Department of Physics, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Andrew A. Peterson
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
10
|
Zhao Q, Martirez JMP, Carter EA. Electrochemical Hydrogenation of CO on Cu(100): Insights from Accurate Multiconfigurational Wavefunction Methods. J Phys Chem Lett 2022; 13:10282-10290. [PMID: 36305601 DOI: 10.1021/acs.jpclett.2c02444] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Copper (Cu) remains the most efficacious electrocatalyst for electrochemical CO2 reduction (CO2R). Its activity and selectivity are highly facet-dependent. We recently examined the commonly proposed rate-limiting CO hydrogenation step on Cu(111) via embedded correlated wavefunction (ECW) theory and demonstrated that only this higher-level theory yields predictions consistent with potential-dependent experimental kinetics. Here, to understand the differing activities of Cu(111) and Cu(100) in catalyzing CO2R, we explore CO hydrogenation on Cu(100) using ECW theory. We predict that the preferred pathway involves the reduction of adsorbed CO (*CO) to *COH via proton-coupled electron transfer (PCET) at working potentials, although *CHO also may form with a kinetically accessible but higher barrier. In contrast, our earlier work on Cu(111) concluded that *COH and *CHO formation via PCET are equally feasible. This work illustrates one possible origin of the facet dependence of CO2R mechanisms and products on Cu electrodes and sheds light on how the selectivity of CO2R electrocatalysts can be controlled by the surface morphology.
Collapse
Affiliation(s)
- Qing Zhao
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| | - John Mark P Martirez
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095-1592, United States
| | - Emily A Carter
- Department of Mechanical and Aerospace Engineering and the Andlinger Center for Energy and the Environment, Princeton University, Princeton, New Jersey 08544-5263, United States
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095-1592, United States
| |
Collapse
|
11
|
Nandy A, Adamji H, Kastner DW, Vennelakanti V, Nazemi A, Liu M, Kulik HJ. Using Computational Chemistry To Reveal Nature’s Blueprints for Single-Site Catalysis of C–H Activation. ACS Catal 2022. [DOI: 10.1021/acscatal.2c02096] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Husain Adamji
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - David W. Kastner
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Vyshnavi Vennelakanti
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Azadeh Nazemi
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Mingjie Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
12
|
Rana B, Coons MP, Herbert JM. Detection and Correction of Delocalization Errors for Electron and Hole Polarons Using Density-Corrected DFT. J Phys Chem Lett 2022; 13:5275-5284. [PMID: 35674719 DOI: 10.1021/acs.jpclett.2c01187] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Modeling polaron defects is an important aspect of computational materials science, but the description of unpaired spins in density functional theory (DFT) often suffers from delocalization error. To diagnose and correct the overdelocalization of spin defects, we report an implementation of density-corrected (DC-)DFT and its analytic energy gradient. In DC-DFT, an exchange-correlation functional is evaluated using a Hartree-Fock density, thus incorporating electron correlation while avoiding self-interaction error. Results for an electron polaron in models of titania and a hole polaron in Al-doped silica demonstrate that geometry optimization with semilocal functionals drives significant structural distortion, including the elongation of several bonds, such that subsequent single-point calculations with hybrid functionals fail to afford a localized defect even in cases where geometry optimization with the hybrid functional does localize the polaron. This has significant implications for traditional workflows in computational materials science, where semilocal functionals are often used for structure relaxation. DC-DFT calculations provide a mechanism to detect situations where delocalization error is likely to affect the results.
Collapse
Affiliation(s)
- Bhaskar Rana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Marc P Coons
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
13
|
Gonçalves JPZ, Seraglio J, Macuvele DLP, Padoin N, Soares C, Riella HG. Green synthesis of manganese based nanoparticles mediated by Eucalyptus robusta and Corymbia citriodora for agricultural applications. Colloids Surf A Physicochem Eng Asp 2022. [DOI: 10.1016/j.colsurfa.2021.128180] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
14
|
Kuklin MS, Eklund K, Linnera J, Ropponen A, Tolvanen N, Karttunen AJ. Structural Properties and Magnetic Ground States of 100 Binary d-Metal Oxides Studied by Hybrid Density Functional Methods. Molecules 2022; 27:molecules27030874. [PMID: 35164135 PMCID: PMC8838575 DOI: 10.3390/molecules27030874] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2021] [Revised: 01/23/2022] [Accepted: 01/24/2022] [Indexed: 02/04/2023] Open
Abstract
d-metal oxides play a crucial role in numerous technological applications and show a great variety of magnetic properties. We have systematically investigated the structural properties, magnetic ground states, and fundamental electronic properties of 100 binary d-metal oxides using hybrid density functional methods and localized basis sets composed of Gaussian-type functions. The calculated properties are compared with experimental information in all cases where experimental data are available. The used PBE0 hybrid density functional method describes the structural properties of the studied d-metal oxides well, except in the case of molecular oxides with weak intermolecular forces between the molecular units. Empirical D3 dispersion correction does not improve the structural description of the molecular oxides. We provide a database of optimized geometries and magnetic ground states to facilitate future studies on the more complex properties of the binary d-metal oxides.
Collapse
|
15
|
Bajaj A, Kulik HJ. Eliminating Delocalization Error to Improve Heterogeneous Catalysis Predictions with Molecular DFT + U. J Chem Theory Comput 2022; 18:1142-1155. [PMID: 35081711 DOI: 10.1021/acs.jctc.1c01178] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Approximate semilocal density functional theory (DFT) is known to underestimate surface formation energies yet paradoxically overbind adsorbates on catalytic transition-metal oxide surfaces due to delocalization error. The low-cost DFT + U approach only improves surface formation energies for early transition-metal oxides or adsorption energies for late transition-metal oxides. In this work, we demonstrate that this inefficacy arises due to the conventional usage of metal-centered atomic orbitals as projectors within DFT + U. We analyze electron density rearrangement during surface formation and O atom adsorption on rutile transition-metal oxides to highlight that a standard DFT + U correction fails to tune properties when the corresponding density rearrangement is highly delocalized across both metal and oxygen sites. To improve both surface properties simultaneously while retaining the simplicity of a single-site DFT + U correction, we systematically construct multi-atom-centered molecular-orbital-like projectors for DFT + U. We demonstrate this molecular DFT + U approach for tuning adsorption energies and surface formation energies of minimal two-dimensional models of representative early (i.e., TiO2) and late (i.e., PtO2) transition-metal oxides. Molecular DFT + U simultaneously corrects adsorption energies and surface formation energies of multilayer models of rutile TiO2(110) and PtO2(110) to resolve the paradoxical description of surface stability and surface reactivity of semilocal DFT.
Collapse
Affiliation(s)
- Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
16
|
Wu L, Lin J, Ren L, Li QN, Chi X, Luo L, Zhang Y, Zeng MH. Electronic Structures of Hydroxyled Low Index Surfaces of Rutile and Anatase-Type Titanium Dioxide. Phys Chem Chem Phys 2022; 24:15091-15102. [DOI: 10.1039/d1cp04729a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Different surface planes of various types of TiO2 crystals have diverse catalysis effects on the splitting of H2O and H2 and the electronic structures of the formed hydroxylated titanium dioxide...
Collapse
|
17
|
Zhao Q, Martirez JMP, Carter EA. Revisiting Understanding of Electrochemical CO 2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory. J Am Chem Soc 2021; 143:6152-6164. [PMID: 33851840 DOI: 10.1021/jacs.1c00880] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Copper (Cu) electrodes, as the most efficacious of CO2 reduction reaction (CO2RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO2 reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO2RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO2RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO2RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at -0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.
Collapse
Affiliation(s)
- Qing Zhao
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| | - John Mark P Martirez
- Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095-1592, United States
| | - Emily A Carter
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States.,Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, California 90095-1592, United States.,Office of the Chancellor, Box 951405, University of California, Los Angeles, Los Angeles, California 90095-1405, United States
| |
Collapse
|
18
|
The Role of Metal Ions in the Electron Transport through Azurin-Based Junctions. APPLIED SCIENCES-BASEL 2021. [DOI: 10.3390/app11093732] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We studied the coherent electron transport through metal–protein–metal junctions based on a blue copper azurin, in which the copper ion was replaced by three different metal ions (Co, Ni and Zn). Our results show that neither the protein structure nor the transmission at the Fermi level change significantly upon metal replacement. The discrepancy with previous experimental observations suggests that the transport mechanism taking place in these types of junctions is probably not fully coherent.
Collapse
|
19
|
Brown JJ, Page AJ. Reaction pathways in the solid state and the Hubbard U correction. J Chem Phys 2021; 154:124121. [PMID: 33810657 DOI: 10.1063/5.0045526] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
We investigate how the Hubbard U correction influences vacancy defect migration barriers in transition metal oxide semiconductors. We show that, depending on the occupation of the transition metal d orbitals, the Hubbard U correction can cause severe instabilities in the migration barrier energies predicted using generalized gradient approximation density functional theory (GGA DFT). For the d0 oxide SrTiO3, applying a Hubbard correction to the Ti4+ 3d orbitals below 4-5 eV yields a migration barrier of ∼0.4 eV. However, above this threshold, the barrier increases suddenly to ∼2 eV. This sudden increase in the transition state barrier arises from the Hubbard U correction changing the Ti4+ t2g/eg orbital occupation, and hence electron density localization, along the migration pathway. Similar results are observed in the d10 oxide ZnO; however, significantly larger Hubbard U corrections must be applied to the Zn2+ 3d orbitals for the same instability to be observed. These results highlight important limitations to the application of the Hubbard U correction when modeling reactive pathways in solid state materials using GGA DFT.
Collapse
Affiliation(s)
- Joshua J Brown
- School of Environmental and Life Sciences, The University of Newcastle, Callaghan, NSW 2308, Australia
| | - Alister J Page
- School of Environmental and Life Sciences, The University of Newcastle, Callaghan, NSW 2308, Australia
| |
Collapse
|
20
|
Morales-Vidal J, García-Muelas R, Ortuño MA. Defects as catalytic sites for the oxygen evolution reaction in Earth-abundant MOF-74 revealed by DFT. Catal Sci Technol 2021. [DOI: 10.1039/d0cy02163f] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The oxygen evolution reaction (OER) is the bottleneck of hydrogen production via water splitting and understanding electrocatalysts at atomic level becomes paramount to enhance the efficiency of this process.
Collapse
Affiliation(s)
- Jordi Morales-Vidal
- Institute of Chemical Research of Catalonia (ICIQ)
- Barcelona Institute of Science and Technology (BIST)
- 43007 Tarragona
- Spain
| | - Rodrigo García-Muelas
- Institute of Chemical Research of Catalonia (ICIQ)
- Barcelona Institute of Science and Technology (BIST)
- 43007 Tarragona
- Spain
| | - Manuel A. Ortuño
- Institute of Chemical Research of Catalonia (ICIQ)
- Barcelona Institute of Science and Technology (BIST)
- 43007 Tarragona
- Spain
| |
Collapse
|
21
|
Brown JJ, Page AJ. The Hubbard-U correction and optical properties of d 0 metal oxide photocatalysts. J Chem Phys 2020; 153:224116. [PMID: 33317276 DOI: 10.1063/5.0027080] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
We report a systematic investigation of individual and multisite Hubbard-U corrections for the electronic, structural, and optical properties of the metal titanate oxide d0 photocatalysts SrTiO3 and rutile/anatase TiO2. Accurate bandgaps for these materials can be reproduced with local density approximation and generalized gradient approximation exchange-correlation density functionals via a continuous series of empirically derived Ud and Up combinations, which are relatively insensitive to the choice of functional. On the other hand, lattice parameters are much more sensitive to the choice of Ud and Up, but in a systematic way that enables the Ud and Up corrections to be used to qualitatively gauge the extent of self-interaction error in the electron density. Modest Ud corrections (e.g., 4 eV-5 eV) yield the most reliable dielectric response functions for SrTiO3 and are comparable to the range of Ud values derived via linear response approaches. For r-TiO2 and a-TiO2, however, the Ud,p corrections that yield accurate bandgaps fail to accurately describe both the parallel and perpendicular components of the dielectric response function. Analysis of individual Ud and Up corrections on the optical properties of SrTiO3 suggests that the most consequential of the two individual corrections is Ud, as it predominately determines the accuracy of the dominant excitation from O-2p to the Ti-3d t2g/eg orbitals. Up, on the other hand, can be used to shift the entire optical response uniformly to higher frequencies. These results will assist high-throughput and machine learning approaches to screening photoactive materials based on d0 photocatalysts.
Collapse
Affiliation(s)
- Joshua J Brown
- School of Environmental and Life Sciences, The University of Newcastle, Callaghan 2308, NSW, Australia
| | - Alister J Page
- School of Environmental and Life Sciences, The University of Newcastle, Callaghan 2308, NSW, Australia
| |
Collapse
|
22
|
Zhao Q, Zhang X, Martirez JMP, Carter EA. Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H2 Desorption from and CH4 Dissociation on Cu(111). J Chem Theory Comput 2020; 16:7078-7088. [DOI: 10.1021/acs.jctc.0c00341] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Qing Zhao
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| | - Xing Zhang
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| | - John Mark P. Martirez
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| | - Emily A. Carter
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| |
Collapse
|
23
|
Choudhuri I, Truhlar DG. Calculating and Characterizing the Charge Distributions in Solids. J Chem Theory Comput 2020; 16:5884-5892. [PMID: 32544328 DOI: 10.1021/acs.jctc.0c00440] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Accurate estimation of the partial atomic charges on metal centers is useful for understanding electronic and catalytic properties of materials. However, different methods of calculating these charges may give quite different results; this issue has been more widely studied for molecules than for solids. Here we study the charges on the metal centers of a test set of 18 solids containing transition metals by using density functional theory with several density functionals (PBE, PBE+U, TPSS, revTPSS, HLE17, revM06-L, B3LYP, B3LYP*, and other exchange-modified B3LYP functionals) and four charge models (Bader, Hirshfeld, CM5, and DDEC6). The test set contains 12 systems with nonmagnetic metal centers (eight metal oxides (MO2), two metal sulfides (MS2), and two metal selenides (MSe2)) and six ferromagnetic transition metal complexes. Our study shows that, among the four types of charges, Bader charges are the highest and Hirshfeld charges are the lowest for all the systems, regardless of the functional being used. The CM5 charges are bigger than DDEC6 charges for MX2 with M = Ti or Mo and X = S or Se, but for the other 14 cases they are lower. We found that the most of the systems are sensitive to the Hubbard U parameters in PBE+U and to the percentage X of Hartree-Fock exchange in exchange-modified B3LYP; as we increase U or X, the charges on the metal atoms in MX2 increase steadily. Testing different density functionals shows charges calculated with higher Hubbard U parameters in PBE+U are comparable to B3LYP (with 20% Hartree-Fock exchange). Among four meta-GGA functionals studied, the charges with HLE17 have the closest agreement with B3LYP. The variation of charges with choice of charge model is greater than the variation with choice of density functional.
Collapse
Affiliation(s)
- Indrani Choudhuri
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
24
|
Rosen AS, Notestein JM, Snurr RQ. Comparing GGA, GGA+U, and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks. J Chem Phys 2020; 152:224101. [DOI: 10.1063/5.0010166] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrew S. Rosen
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, USA
| | - Justin M. Notestein
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, USA
| | - Randall Q. Snurr
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, USA
| |
Collapse
|
25
|
Nandy A, Chu DBK, Harper DR, Duan C, Arunachalam N, Cytter Y, Kulik HJ. Large-scale comparison of 3d and 4d transition metal complexes illuminates the reduced effect of exchange on second-row spin-state energetics. Phys Chem Chem Phys 2020; 22:19326-19341. [DOI: 10.1039/d0cp02977g] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The origin of distinct 3d vs. 4d transition metal complex sensitivity to exchange is explored over a large data set.
Collapse
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
- Department of Chemistry
| | - Daniel B. K. Chu
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| | - Daniel R. Harper
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
- Department of Chemistry
| | - Chenru Duan
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
- Department of Chemistry
| | - Naveen Arunachalam
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| | - Yael Cytter
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| | - Heather J. Kulik
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| |
Collapse
|
26
|
Liu F, Kulik HJ. Impact of Approximate DFT Density Delocalization Error on Potential Energy Surfaces in Transition Metal Chemistry. J Chem Theory Comput 2019; 16:264-277. [DOI: 10.1021/acs.jctc.9b00842] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
27
|
Laplaza R, Polo V, Contreras-García J. Localizing electron density errors in density functional theory. Phys Chem Chem Phys 2019; 21:20927-20938. [PMID: 31517339 DOI: 10.1039/c9cp02806d] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The accuracy of different density functional approximations is assessed through the use of quantum chemical topology on molecular electron densities. In particular, three simple yet ever-important systems are studied: N2, CO and ethane. Our results exemplify how real-space descriptors can help understand the sources of errors in density functional theory, avoiding unwanted error compensation present in simplified statistical metrics. Errors in "well-built" functionals are shown to be concentrated in chemically meaningful regions of space, and hence they are predictable. Conversely, strongly parametrized functionals show isotropic errors that cannot be traced back to chemically transferable units. Moreover, we will show that energetic corrections are mapped back into improvements in the density in chemically meaningful regions. These results point at the relevance of real-space perspectives when parametrizing or relating energy and density errors.
Collapse
Affiliation(s)
- Rubén Laplaza
- Sorbonne Université, CNRS, Laboratoire de Chimie Thèorique, LCT, F. 75005 Paris, France and Departamento de Química Física and Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza, E. 50009, Zaragoza, Spain
| | - Victor Polo
- Departamento de Química Física and Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza, E. 50009, Zaragoza, Spain
| | - Julia Contreras-García
- Sorbonne Université, CNRS, Laboratoire de Chimie Thèorique, LCT, F. 75005 Paris, France.
| |
Collapse
|
28
|
Lesage C, Devers E, Legens C, Fernandes G, Roudenko O, Briois V. High pressure cell for edge jumping X-ray absorption spectroscopy: Applications to industrial liquid sulfidation of hydrotreatment catalysts. Catal Today 2019. [DOI: 10.1016/j.cattod.2019.01.081] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
29
|
Zhao Q, Kulik HJ. Stable Surfaces That Bind Too Tightly: Can Range-Separated Hybrids or DFT+U Improve Paradoxical Descriptions of Surface Chemistry? J Phys Chem Lett 2019; 10:5090-5098. [PMID: 31411023 PMCID: PMC6748670 DOI: 10.1021/acs.jpclett.9b01650] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Accepted: 08/14/2019] [Indexed: 05/25/2023]
Abstract
Approximate, semilocal density functional theory (DFT) suffers from delocalization error that can lead to a paradoxical model of catalytic surfaces that both overbind adsorbates yet are also too stable. We investigate the effect of two widely applied approaches for delocalization error correction, (i) affordable DFT+U (i.e., semilocal DFT augmented with a Hubbard U) and (ii) hybrid functionals with an admixture of Hartree-Fock (HF) exchange, on surface and adsorbate energies across a range of rutile transition metal oxides widely studied for their promise as water-splitting catalysts. We observe strongly row- and period-dependent trends with DFT+U, which increases surface formation energies only in early transition metals (e.g., Ti and V) and decreases adsorbate energies only in later transition metals (e.g., Ir and Pt). Both global and local hybrids destabilize surfaces and reduce adsorbate binding across the periodic table, in agreement with higher-level reference calculations. Density analysis reveals why hybrid functionals correct both quantities, whereas DFT+U does not. We recommend local, range-separated hybrids for the accurate modeling of catalysis in transition metal oxides at only a modest increase in computational cost over semilocal DFT.
Collapse
Affiliation(s)
- Qing Zhao
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
- Department
of Mechanical Engineering, Massachusetts
Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
30
|
Gaggioli CA, Stoneburner SJ, Cramer CJ, Gagliardi L. Beyond Density Functional Theory: The Multiconfigurational Approach To Model Heterogeneous Catalysis. ACS Catal 2019. [DOI: 10.1021/acscatal.9b01775] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Carlo Alberto Gaggioli
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Samuel J. Stoneburner
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Christopher J. Cramer
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
31
|
Cirera J, Ruiz E. Computational Modeling of Transition Temperatures in Spin-Crossover Systems. COMMENT INORG CHEM 2019. [DOI: 10.1080/02603594.2019.1608967] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Jordi Cirera
- Departament de Química Inorgànica i Orgànica and Institut de Recerca de Química Teòrica i Computacional, Universitat de Barcelona, Barcelona, Spain
| | - Eliseo Ruiz
- Departament de Química Inorgànica i Orgànica and Institut de Recerca de Química Teòrica i Computacional, Universitat de Barcelona, Barcelona, Spain
| |
Collapse
|
32
|
Bajaj A, Liu F, Kulik HJ. Non-empirical, low-cost recovery of exact conditions with model-Hamiltonian inspired expressions in jmDFT. J Chem Phys 2019; 150:154115. [PMID: 31005112 DOI: 10.1063/1.5091563] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Density functional theory (DFT) is widely applied to both molecules and materials, but well known energetic delocalization and static correlation errors in practical exchange-correlation approximations limit quantitative accuracy. Common methods that correct energetic delocalization errors, such as the Hubbard U correction in DFT+U or Hartree-Fock exchange in global hybrids, do so at the cost of worsening static correlation errors. We recently introduced an alternate approach [Bajaj et al., J. Chem. Phys. 147, 191101 (2017)] known as judiciously modified DFT (jmDFT), wherein the deviation from exact behavior of semilocal functionals over both fractional spin and charge, i.e., the so-called flat plane, was used to motivate functional forms of second order analytic corrections. In this work, we introduce fully nonempirical expressions for all four coefficients in a DFT+U+J-inspired form of jmDFT, where all coefficients are obtained only from energies and eigenvalues of the integer-electron systems. We show good agreement for U and J coefficients obtained nonempirically as compared with the results of numerical fitting in a jmDFT U+J/J' correction. Incorporating the fully nonempirical jmDFT correction reduces and even eliminates the fractional spin error at the same time as eliminating the energetic delocalization error. We show that this approach extends beyond s-electron systems to higher angular momentum cases including p- and d-electrons. Finally, we diagnose some shortcomings of the current jmDFT approach that limit its ability to improve upon DFT results for cases such as weakly bound anions due to poor underlying semilocal functional behavior.
Collapse
Affiliation(s)
- Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
33
|
Liu F, Yang T, Yang J, Xu E, Bajaj A, Kulik HJ. Bridging the Homogeneous-Heterogeneous Divide: Modeling Spin for Reactivity in Single Atom Catalysis. Front Chem 2019; 7:219. [PMID: 31041303 PMCID: PMC6476907 DOI: 10.3389/fchem.2019.00219] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Accepted: 03/20/2019] [Indexed: 12/03/2022] Open
Abstract
Single atom catalysts (SACs) are emergent catalytic materials that have the promise of merging the scalability of heterogeneous catalysts with the high activity and atom economy of homogeneous catalysts. Computational, first-principles modeling can provide essential insight into SAC mechanism and active site configuration, where the sub-nm-scale environment can challenge even the highest-resolution experimental spectroscopic techniques. Nevertheless, the very properties that make SACs attractive in catalysis, such as localized d electrons of the isolated transition metal center, make them challenging to study with conventional computational modeling using density functional theory (DFT). For example, Fe/N-doped graphitic SACs have exhibited spin-state dependent reactivity that remains poorly understood. However, spin-state ordering in DFT is very sensitive to the nature of the functional approximation chosen. In this work, we develop accurate benchmarks from correlated wavefunction theory (WFT) for relevant octahedral complexes. We use those benchmarks to evaluate optimal DFT functional choice for predicting spin state ordering in small octahedral complexes as well as models of pyridinic and pyrrolic nitrogen environments expected in larger SACs. Using these guidelines, we determine Fe/N-doped graphene SAC model properties and reactivity as well as their sensitivities to DFT functional choice. Finally, we conclude with broad recommendations for computational modeling of open-shell transition metal single-atom catalysts.
Collapse
Affiliation(s)
- Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Tzuhsiung Yang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Jing Yang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Eve Xu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States.,Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, United States
| |
Collapse
|
34
|
Janet JP, Liu F, Nandy A, Duan C, Yang T, Lin S, Kulik HJ. Designing in the Face of Uncertainty: Exploiting Electronic Structure and Machine Learning Models for Discovery in Inorganic Chemistry. Inorg Chem 2019; 58:10592-10606. [PMID: 30834738 DOI: 10.1021/acs.inorgchem.9b00109] [Citation(s) in RCA: 59] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Recent transformative advances in computing power and algorithms have made computational chemistry central to the discovery and design of new molecules and materials. First-principles simulations are increasingly accurate and applicable to large systems with the speed needed for high-throughput computational screening. Despite these strides, the combinatorial challenges associated with the vastness of chemical space mean that more than just fast and accurate computational tools are needed for accelerated chemical discovery. In transition-metal chemistry and catalysis, unique challenges arise. The variable spin, oxidation state, and coordination environments favored by elements with well-localized d or f electrons provide great opportunity for tailoring properties in catalytic or functional (e.g., magnetic) materials but also add layers of uncertainty to any design strategy. We outline five key mandates for realizing computationally driven accelerated discovery in inorganic chemistry: (i) fully automated simulation of new compounds, (ii) knowledge of prediction sensitivity or accuracy, (iii) faster-than-fast property prediction methods, (iv) maps for rapid chemical space traversal, and (v) a means to reveal design rules on the kilocompound scale. Through case studies in open-shell transition-metal chemistry, we describe how advances in methodology and software in each of these areas bring about new chemical insights. We conclude with our outlook on the next steps in this process toward realizing fully autonomous discovery in inorganic chemistry using computational chemistry.
Collapse
Affiliation(s)
- Jon Paul Janet
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Fang Liu
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Aditya Nandy
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States.,Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Chenru Duan
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States.,Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Tzuhsiung Yang
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Sean Lin
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Heather J Kulik
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
35
|
Kröncke S, Herrmann C. Designing Long-Range Charge Delocalization from First-Principles. J Chem Theory Comput 2018; 15:165-177. [DOI: 10.1021/acs.jctc.8b00872] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Susanne Kröncke
- Department of Chemistry, University of Hamburg, Hamburg 20146, Germany
| | - Carmen Herrmann
- Department of Chemistry, University of Hamburg, Hamburg 20146, Germany
| |
Collapse
|
36
|
Kuklin MS, Karttunen AJ. Crystal Structure Prediction of Magnetic Transition-Metal Oxides by Using Evolutionary Algorithm and Hybrid DFT Methods. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2018; 122:24949-24957. [PMID: 30416641 PMCID: PMC6221369 DOI: 10.1021/acs.jpcc.8b08238] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/24/2018] [Revised: 10/10/2018] [Indexed: 06/09/2023]
Abstract
Although numerous crystal structures have been successfully predicted by using currently available computational techniques, prediction of strongly correlated systems such as transition-metal oxides remains a challenge. To overcome this problem, we have interfaced evolutionary algorithm-based USPEX method with the CRYSTAL code, enabling the use of Gaussian-type localized atomic basis sets and hybrid density functional (DFT) methods for the prediction of crystal structures. We report successful crystal structure predictions of several transition-metal oxides (NiO, CoO, α-Fe2O3, V2O3, and CuO) with correct atomic magnetic moments, spin configurations, and structures by using the USPEX method in combination with the CRYSTAL code and Perdew-Burke-Ernzerhof (PBE0) hybrid functional. Our benchmarking results demonstrate that USPEX + hybrid DFT is a suitable combination to reliably predict the magnetic structures of strongly correlated materials.
Collapse
Affiliation(s)
- Mikhail S. Kuklin
- Department of Chemistry and
Materials Science, Aalto University, P.O. Box 16100, FI-00076 Aalto, Finland
| | - Antti J. Karttunen
- Department of Chemistry and
Materials Science, Aalto University, P.O. Box 16100, FI-00076 Aalto, Finland
| |
Collapse
|
37
|
Valero R, Morales-García Á, Illas F. Theoretical Modeling of Electronic Excitations of Gas-Phase and Solvated TiO 2 Nanoclusters and Nanoparticles of Interest in Photocatalysis. J Chem Theory Comput 2018; 14:4391-4404. [PMID: 30011198 DOI: 10.1021/acs.jctc.8b00651] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The optical absorption spectra of (TiO2) n, nanoclusters ( n = 1-20) and nanoparticles ( n = 35, 84) have been calculated from the frequency-dependent dielectric function in the independent particle approximation under the framework of density functional theory. The PBE generalized gradient approach based functional, the so-called PBE+ U method and the PBE0 and PBEx hybrid functionals-containing 25% and 12.5% of nonlocal Fock exchange, respectively-have been used. The simulated spectra have been obtained in the gas phase and in water on previously PBE0 optimized atomic structures. The effect of the solvent has been accounted for by using an implicit water solvation model. For the smallest nanoclusters, the spectra show discrete peaks, whereas for the largest nanoclusters and for the nanoparticles they resemble a continuum absorption band. In the gas phase and for a given density functional, the onset of the absorption (optical gap, Ogap) remains relatively constant for all nanoparticle sizes although it increases with the percentage of nonlocal Fock exchange, as expected. For all tested functionals, the tendency of Ogap in water is very similar to that observed in the gas phase with an almost constant upshift. For comparison, the optical gap has also been calculated at the TD-DFT level with the PBEx functional in the gas phase and in water. Both approaches agree reasonably well although the TD-DFT gap values are lower than those derived from the dielectric-function. Overall, the position of the spectral maxima and the width of the spectra are relatively constant and independent of particle size which may have implications in the understanding of photocatalysis by TiO2.
Collapse
Affiliation(s)
- Rosendo Valero
- Departament de Ciència de Materials i Química Física & Institut de Química Teórica i Computacional (IQTCUB) , Universitat de Barcelona , c/Martí i Franquès 1 , 08028 , Barcelona , Spain
| | - Ángel Morales-García
- Departament de Ciència de Materials i Química Física & Institut de Química Teórica i Computacional (IQTCUB) , Universitat de Barcelona , c/Martí i Franquès 1 , 08028 , Barcelona , Spain
| | - Francesc Illas
- Departament de Ciència de Materials i Química Física & Institut de Química Teórica i Computacional (IQTCUB) , Universitat de Barcelona , c/Martí i Franquès 1 , 08028 , Barcelona , Spain
| |
Collapse
|
38
|
Romero-Muñiz C, Ortega M, Vilhena JG, Díez-Pérez I, Cuevas JC, Pérez R, Zotti LA. Ab initio electronic structure calculations of entire blue copper azurins. Phys Chem Chem Phys 2018; 20:30392-30402. [DOI: 10.1039/c8cp06862c] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
We present a theoretical study of the blue-copper azurin extracted from Pseudomonas aeruginosa and several of its single amino acid mutants.
Collapse
Affiliation(s)
- Carlos Romero-Muñiz
- Departamento de Física Teórica de la Materia Condensada
- Universidad Autónoma de Madrid
- E-28049 Madrid
- Spain
| | - María Ortega
- Departamento de Física Teórica de la Materia Condensada
- Universidad Autónoma de Madrid
- E-28049 Madrid
- Spain
| | - J. G. Vilhena
- Departamento de Física Teórica de la Materia Condensada
- Universidad Autónoma de Madrid
- E-28049 Madrid
- Spain
- Department of Physics
| | - I. Díez-Pérez
- Department of Materials Science and Physical Chemistry & Institute of Theoretical and Computational Chemistry (IQTCUB)
- University of Barcelona
- Barcelona 08028
- Spain
- Department of Chemistry
| | - Juan Carlos Cuevas
- Departamento de Física Teórica de la Materia Condensada
- Universidad Autónoma de Madrid
- E-28049 Madrid
- Spain
- Condensed Matter Physics Center (IFIMAC)
| | - Rubén Pérez
- Departamento de Física Teórica de la Materia Condensada
- Universidad Autónoma de Madrid
- E-28049 Madrid
- Spain
- Condensed Matter Physics Center (IFIMAC)
| | - Linda A. Zotti
- Departamento de Física Teórica de la Materia Condensada
- Universidad Autónoma de Madrid
- E-28049 Madrid
- Spain
- Condensed Matter Physics Center (IFIMAC)
| |
Collapse
|