51
|
Nakata M, Maeda T. PubChemQC B3LYP/6-31G*//PM6 Data Set: The Electronic Structures of 86 Million Molecules Using B3LYP/6-31G* Calculations. J Chem Inf Model 2023; 63:5734-5754. [PMID: 37677147 DOI: 10.1021/acs.jcim.3c00899] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/09/2023]
Abstract
The presented "PubChemQC B3LYP/6-31G*//PM6" data set is composed of the electronic properties of 85,938,443 molecules, encompassing a broad spectrum of molecules from essential compounds to biomolecules with a molecular weight up to 1000. These molecules account for 94.0% of the original PubChem Compound catalog as of August 29, 2016. The electronic properties, including orbitals, orbital energies, total energies, dipole moments, and other pertinent properties, were computed by using the B3LYP/6-31G* and PM6 methods. The data set, available in three formats, namely, GAMESS quantum chemistry program files, selected JSON output files, and a PostgreSQL database, provides researchers with the ability to query molecular properties. It is further subdivided into five subdata sets for more specific data. The first two subsets encompass molecules with carbon, hydrogen, oxygen, and nitrogen with molecular weights under 300 and 500, respectively. The third and fourth subsets incorporate molecules with carbon, hydrogen, nitrogen, oxygen, phosphorus, sulfur, fluorine, and chlorine, with molecular weights under 300 and 500, respectively. The fifth subset comprises molecules with carbon, hydrogen, nitrogen, oxygen, phosphorus, sulfur, fluorine, chlorine, sodium, potassium, magnesium, and calcium, with a molecular weight of under 500. The coefficients of determination for the highest occupied molecular orbital-lowest unoccupied molecular orbital energy gap range from 0.892 (for CHON500) to 0.803 (for the whole data set). These comprehensive results pave the way for applications in drug discovery and materials science, among others. The data sets can be accessed under the Creative Commons Attribution 4.0 International license at the following web address: https://nakatamaho.riken.jp/pubchemqc.riken.jp/b3lyp_pm6_datasets.html.
Collapse
Affiliation(s)
- Maho Nakata
- RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Toshiyuki Maeda
- Software Technology and Artificial Intelligence Research Laboratory, Chiba Institute of Technology, 2-17-1 Tsudanuma, Narashino, Chiba 275-0016, Japan
| |
Collapse
|
52
|
Lehtola S, Marques MAL. Reproducibility of density functional approximations: How new functionals should be reported. J Chem Phys 2023; 159:114116. [PMID: 37725491 DOI: 10.1063/5.0167763] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2023] [Accepted: 08/28/2023] [Indexed: 09/21/2023] Open
Abstract
Density functional theory is the workhorse of chemistry and materials science, and novel density functional approximations are published every year. To become available in program packages, the novel density functional approximations (DFAs) need to be (re)implemented. However, according to our experience as developers of Libxc [Lehtola et al., SoftwareX 7, 1 (2018)], a constant problem in this task is verification due to the lack of reliable reference data. As we discuss in this work, this lack has led to several non-equivalent implementations of functionals such as Becke-Perdew 1986, Perdew-Wang 1991, Perdew-Burke-Ernzerhof, and Becke's three-parameter hybrid functional with Lee-Yang-Parr correlation across various program packages, yielding different total energies. Through careful verification, we have also found many issues with incorrect functional forms in recent DFAs. The goal of this work is to ensure the reproducibility of DFAs. DFAs must be verifiable in order to prevent the reappearance of the above-mentioned errors and incompatibilities. A common framework for verification and testing is, therefore, needed. We suggest several ways in which reference energies can be produced with free and open source software, either with non-self-consistent calculations with tabulated atomic densities or via self-consistent calculations with various program packages. The employed numerical parameters-especially the quadrature grid-need to be converged to guarantee a ≲0.1 μEh precision in the total energy, which is nowadays routinely achievable in fully numerical calculations. Moreover, as such sub-μEh level agreement can only be achieved when fully equivalent implementations of the DFA are used, the source code of the reference implementation should also be made available in any publication describing a new DFA.
Collapse
Affiliation(s)
- Susi Lehtola
- Molecular Sciences Software Institute, Blacksburg, Virginia 24061, USA
- Department of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
| | - Miguel A L Marques
- Research Center Future Energy Materials and Systems of the University Alliance Ruhr, Faculty of Mechanical Engineering, Ruhr University Bochum, Universitätsstraße 150, D-44801 Bochum, Germany
| |
Collapse
|
53
|
Duan X, Li Y, Zhao C, Shen Y, Guo Q, Huang Z, Shan D, Gao Y, Zhang K, Shi J, Liu J, Chen Y, Yuan CG. Efficient immobilization and detoxification of gaseous elemental mercury by nanoflower/rod WSe 2/halloysite composite: Performance and mechanisms. JOURNAL OF HAZARDOUS MATERIALS 2023; 458:131898. [PMID: 37354718 DOI: 10.1016/j.jhazmat.2023.131898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/04/2023] [Revised: 05/29/2023] [Accepted: 06/18/2023] [Indexed: 06/26/2023]
Abstract
Gaseous mercury pollution control technologies with low stability and high releasing risks always face with great challenges. Herein, we developed one halloysite nanotubes (HNTs)-supported tungsten diselenide (WSe2) composite (WSe2/HNTs) by one-pot solvothermal approach, curing Hg0 from complicated flue gas (CFG) and reducing second environment risks. WSe2 as a monolayer with nano-flower structure and HNTs with rod shapes in the as-prepared sorbent exhibited outstanding synergy efficiency, resulting in exceptional performance for Hg0 removal with high capture capacity of 30.6 mg·g-1 and rate of 9.09 μg·g-1·min-1, which benefited from the high affinity of selenium and mercury (1 ×1045) and the adequate exposure of Se-terminated. The adsorbent showed beneficial tolerance to high amount of NOx and SOx. An online lab-built thermal decomposition system (TPD-AFS) was employed to explore Hg species on the used-sorbent, finding that the adsorbed-mercury species were principally mercury selenide (HgSe). Density functional theory calculations indicated that the hollow-sites were the major adsorption sites and exhibited excellent selectivity for Hg0, as well as HgSe generation needed to overcome the 0.32 eV energy barrier. The adsorbed mercury displayed high environmental stability after the leaching toxicity test, which significantly decreased its secondary environmental risks. With these advantages, WSe2/HNTs possess enormous potential to achieve the effective and permanent immobilization of gaseous mercury from CFG in the future.
Collapse
Affiliation(s)
- Xuelei Duan
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China
| | - Yuan Li
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China
| | - Changxian Zhao
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China
| | - Yiwen Shen
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China
| | - Qi Guo
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China
| | - Zhihao Huang
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China
| | - Dexu Shan
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China
| | - Yue Gao
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China
| | - Kegang Zhang
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China
| | - Jianbo Shi
- State Key Laboratory of Environmental Chemistry and Ecotoxicology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, P. O. Box 2871, Beijing 100085, China
| | - Jingfu Liu
- State Key Laboratory of Environmental Chemistry and Ecotoxicology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, P. O. Box 2871, Beijing 100085, China
| | - Yongsheng Chen
- Department of Civil & Environmental Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Chun-Gang Yuan
- Hebei Key Lab of Power Plant Flue Gas Multi-Pollutants Control, Department of Environmental Science & Engineering, North China Electric Power University, Baoding 071000, China; MOE Key Laboratory of Resources and Environmental Systems Optimization, College of Environmental Science and Engineering, North China Electric Power University, Beijing 102206, China.
| |
Collapse
|
54
|
Sitha S. Better performance of Hartree-Fock over DFT: a quantum mechanical investigation on pyridinium benzimidazolate types of zwitterions in the light of localization/delocalization issues. J Mol Model 2023; 29:313. [PMID: 37704866 PMCID: PMC10499969 DOI: 10.1007/s00894-023-05706-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 08/22/2023] [Indexed: 09/15/2023]
Abstract
CONTEXT With the advent of fast computing facilities, combined with rapid emerges of many new and intricate quantum mechanical functionals, computations with pure Hartree-Fock (HF) theory are now-a-days regarded as trivial or obsolete, or even considered as not reliable by many researchers. Consequently, current trends in computational chemistry show extensive use of post-HF theories for smaller molecular systems and various DFT methods for organic and inorganic chemistry related problems (larger molecules/systems). In this contribution, I have tried to show that sometimes, HF might be more suitable over DFT methodologies in addressing structure-property correlations. Molecules studied here were previously synthesized by Boyd in 1966 and important experimental data were produced by Alcalde and co-workers in 1987. Comparison of computed and experimental results clearly shows that HF method was more effective in reproducing the experimental data compared to especially the DFT methodologies. Reliability of HF method was further assured from the very similar results shown by the CCSD, CASSCF, CISD and QCISD methods. Current study also indicates that the localization issue associated with HF proved to be advantageous over delocalization issue of DFT based methodologies, in correctly describing the structure-property correlation for zwitterion systems. METHODS All computations were performed with Gaussian 09. A wide-range of quantum mechanical methodologies, HF, B3LYP, CAM-B3LYP, BMK, B3PW91, TPSSh, LC-ωPBE, M06-2X, M06-HF, ωB97xD, MP2, CASSCF, CCSD, QCISD, CISD and semi-empirical methods like, Huckel, CNDO, AM1, PM3MM and PM6, were used for investigations.
Collapse
Affiliation(s)
- Sanyasi Sitha
- Department of Chemical Sciences, University of Johannesburg, Auckland Park, PO Box 524, Johannesburg, 2006, South Africa.
| |
Collapse
|
55
|
França TCC, Valle da Silva JA, Dos Santos MC, Cavalcante SFDA, Kuca K. Applications of the Near Attack Conformation (NAC) approach in the search for Acetylcholinesterase reactivators. Chem Biol Interact 2023; 382:110619. [PMID: 37406983 DOI: 10.1016/j.cbi.2023.110619] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2023] [Revised: 06/25/2023] [Accepted: 06/30/2023] [Indexed: 07/07/2023]
Abstract
The Near Attack Conformation (NAC) approach states that the efficiency of an enzyme-catalyzed reaction depends on the prior attainment of optimal conditions for substrate atom organization and positioning for bond formation. These conditions are prerequisites for the transition state (TS) in which the involved atoms are within the van der Waals range of contact and positioned at an angle similar to that achieved after bond formation. The successful application of this approach to investigate the reactivation mechanism of acetylcholinesterase inhibited by nerve agents has contributed to a better understanding of this mechanism and demonstrated consistent corroboration with experimental data. In this article, we summarize the accomplishments achieved thus far and outline future perspectives.
Collapse
Affiliation(s)
- Tanos Celmar Costa França
- Laboratory of Molecular Modeling Applied to Chemical and Biological Defense, Military Institute of Engineering, Rio de Janeiro, Brazil; Department of Chemistry, Faculty of Science, University of Hradec Kralove, Rokitanskeho 62, 50003, Hradec Kralove, Czech Republic.
| | - Jorge Alberto Valle da Silva
- Laboratory of Molecular Modeling Applied to Chemical and Biological Defense, Military Institute of Engineering, Rio de Janeiro, Brazil
| | - Marcelo Carneiro Dos Santos
- Laboratory of Molecular Modeling Applied to Chemical and Biological Defense, Military Institute of Engineering, Rio de Janeiro, Brazil
| | - Samir Frontino de Almeida Cavalcante
- Department of Chemistry, Faculty of Science, University of Hradec Kralove, Rokitanskeho 62, 50003, Hradec Kralove, Czech Republic; Institute of Chemical, Biological, Radiological and Nuclear Defense (IDQBRN), Brazilian Army Technological Center (CTEx), Rio de Janeiro, RJ, Brazil
| | - Kamil Kuca
- Department of Chemistry, Faculty of Science, University of Hradec Kralove, Rokitanskeho 62, 50003, Hradec Kralove, Czech Republic
| |
Collapse
|
56
|
Dufils T, Knijff L, Shao Y, Zhang C. PiNNwall: Heterogeneous Electrode Models from Integrating Machine Learning and Atomistic Simulation. J Chem Theory Comput 2023; 19:5199-5209. [PMID: 37477645 PMCID: PMC10413855 DOI: 10.1021/acs.jctc.3c00359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2023] [Indexed: 07/22/2023]
Abstract
Electrochemical energy storage always involves the capacitive process. The prevailing electrode model used in the molecular simulation of polarizable electrode-electrolyte systems is the Siepmann-Sprik model developed for perfect metal electrodes. This model has been recently extended to study the metallicity in the electrode by including the Thomas-Fermi screening length. Nevertheless, a further extension to heterogeneous electrode models requires introducing chemical specificity, which does not have any analytical recipes. Here, we address this challenge by integrating the atomistic machine learning code (PiNN) for generating the base charge and response kernel and the classical molecular dynamics code (MetalWalls) dedicated to the modeling of electrochemical systems, and this leads to the development of the PiNNwall interface. Apart from the cases of chemically doped graphene and graphene oxide electrodes as shown in this study, the PiNNwall interface also allows us to probe polarized oxide surfaces in which both the proton charge and the electronic charge can coexist. Therefore, this work opens the door for modeling heterogeneous and complex electrode materials often found in energy storage systems.
Collapse
Affiliation(s)
- Thomas Dufils
- Department of Chemistry-Ångström
Laboratory, Uppsala University, Lägerhyddsvägen 1, P. O. Box 538, 75121 Uppsala, Sweden
| | - Lisanne Knijff
- Department of Chemistry-Ångström
Laboratory, Uppsala University, Lägerhyddsvägen 1, P. O. Box 538, 75121 Uppsala, Sweden
| | - Yunqi Shao
- Department of Chemistry-Ångström
Laboratory, Uppsala University, Lägerhyddsvägen 1, P. O. Box 538, 75121 Uppsala, Sweden
| | - Chao Zhang
- Department of Chemistry-Ångström
Laboratory, Uppsala University, Lägerhyddsvägen 1, P. O. Box 538, 75121 Uppsala, Sweden
| |
Collapse
|
57
|
Ko HY, Calegari Andrade MF, Sparrow ZM, Zhang JA, DiStasio RA. High-Throughput Condensed-Phase Hybrid Density Functional Theory for Large-Scale Finite-Gap Systems: The SeA Approach. J Chem Theory Comput 2023. [PMID: 37385014 DOI: 10.1021/acs.jctc.2c00827] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/01/2023]
Abstract
High-throughput electronic structure calculations (often performed using density functional theory (DFT)) play a central role in screening existing and novel materials, sampling potential energy surfaces, and generating data for machine learning applications. By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal DFT and furnish a more accurate description of the underlying electronic structure, albeit at a computational cost that often prohibits such high-throughput applications. To address this challenge, we have constructed a robust, accurate, and computationally efficient framework for high-throughput condensed-phase hybrid DFT and implemented this approach in the PWSCF module of Quantum ESPRESSO (QE). The resulting SeA approach (SeA = SCDM + exx + ACE) combines and seamlessly integrates: (i) the selected columns of the density matrix method (SCDM, a robust noniterative orbital localization scheme that sidesteps system-dependent optimization protocols), (ii) a recently extended version of exx (a black-box linear-scaling EXX algorithm that exploits sparsity between localized orbitals in real space when evaluating the action of the standard/full-rank V^xx operator), and (iii) adaptively compressed exchange (ACE, a low-rank V^xx approximation). In doing so, SeA harnesses three levels of computational savings: pair selection and domain truncation from SCDM + exx (which only considers spatially overlapping orbitals on orbital-pair-specific and system-size-independent domains) and low-rank V^xx approximation from ACE (which reduces the number of calls to SCDM + exx during the self-consistent field (SCF) procedure). Across a diverse set of 200 nonequilibrium (H2O)64 configurations (with densities spanning 0.4-1.7 g/cm3), SeA provides a 1-2 order-of-magnitude speedup in the overall time-to-solution, i.e., ≈8-26× compared to the convolution-based PWSCF(ACE) implementation in QE and ≈78-247× compared to the conventional PWSCF(Full) approach, and yields energies, ionic forces, and other properties with high fidelity. As a proof-of-principle high-throughput application, we trained a deep neural network (DNN) potential for ambient liquid water at the hybrid DFT level using SeA via an actively learned data set with ≈8,700 (H2O)64 configurations. Using an out-of-sample set of (H2O)512 configurations (at nonambient conditions), we confirmed the accuracy of this SeA-trained potential and showcased the capabilities of SeA by computing the ground-truth ionic forces in this challenging system containing >1,500 atoms.
Collapse
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Marcos F Calegari Andrade
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
- Quantum Simulations Group, Materials Science Division, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Zachary M Sparrow
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Ju-An Zhang
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
58
|
Pogrebetsky J, Siklitskaya A, Kubas A. MP2-Based Correction Scheme to Approach the Limit of a Complete Pair Natural Orbitals Space in DLPNO-CCSD(T) Calculations. J Chem Theory Comput 2023. [PMID: 37338422 DOI: 10.1021/acs.jctc.3c00444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
The domain-based local pair natural orbital (PNO) coupled-cluster DLPNO-CCSD(T) method has been proven to provide accurate single-point energies at a fraction of the cost of canonical CCSD(T) calculations. However, the desired "chemical accuracy" can only be obtained with a large PNO space and extended basis set. We present a simple yet accurate and efficient correction scheme based on a perturbative approach. Here, in addition to DLPNO-CCSD(T) energy, one calculates DLPNO-MP2 correlation energy with the same settings as in the preceding coupled-cluster calculation. In the next step, the canonical MP2 correlation energy is obtained in the same orbital basis. This can be efficiently performed for essentially all molecule sizes accessible with the DLPNO-CCSD(T) method. By taking the difference between the canonical MP2 and DLPNO-MP2 energies, we obtain a correction term that can be added to the DLPNO-CCSD(T) correlation energy. This way, one can obtain the total correlation energy close to the limit of the complete PNO space (cPNO). The presented approach allows us to significantly increase the accuracy of the DLPNO-CCSD(T) method for both closed- and open-shell systems. The latter are known to be especially challenging for locally correlated methods. Unlike the previously developed PNO extrapolation procedure by Altun, Neese, and Bistoni ( J. Chem. Theory Comput. 2020, 16, 6142-6149), this strategy allows us to get the DLPNO-CCSD(T) correlation energy at the cPNO limit in a cost-efficient way, resulting in a minimal overall increase in calculation time as compared to the uncorrected method.
Collapse
Affiliation(s)
- James Pogrebetsky
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, Warszawa 01-224, Poland
| | - Alexandra Siklitskaya
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, Warszawa 01-224, Poland
| | - Adam Kubas
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, Warszawa 01-224, Poland
| |
Collapse
|
59
|
Ren M, Zhu X, Luo Q, Li X, Yang J. Proposal of spin crossover as a reversible switch of catalytic activity for the oxygen evolution reaction in two dimensional metal-organic frameworks. NANOSCALE 2023. [PMID: 37314098 DOI: 10.1039/d3nr01708g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The development of oxygen evolution reaction (OER) catalysts with high activity and controllability is crucial for clean energy conversion and storage but remains a challenge. Here, based on first-principles calculations, we propose to utilize spin crossover (SCO) in two-dimensional (2D) metal-organic frameworks (MOFs) to achieve reversible control of OER catalytic activity. The theoretical design of a 2D square lattice MOF with Co as nodes and tetrakis-substituted cyanimino squaric acid (TCSA) as ligands, which transforms between the high spin (HS) and the low spin (LS) state by applying an external strain (∼2%), confirms our proposal. In particular, the HS-LS spin state transition of Co(TCSA) considerably regulates the adsorption strength of the key intermediate HO* in the OER process, resulting in a significant reduction of the overpotential from 0.62 V in the HS state to 0.32 V in the LS state, thus realizing a reversible switch for the activity of the OER. Moreover, the high activity of the LS state is confirmed by microkinetic and constant potential method simulations.
Collapse
Affiliation(s)
- Min Ren
- Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiangyu Zhu
- Institutes of Physical Science and Information, Department of Chemistry, Anhui University, Hefei, Anhui 230601, China
| | - Qiquan Luo
- Institutes of Physical Science and Information, Department of Chemistry, Anhui University, Hefei, Anhui 230601, China
| | - Xingxing Li
- Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei 230088, China
| | - Jinglong Yang
- Key Laboratory of Precision and Intelligent Chemistry, University of Science and Technology of China, Hefei, Anhui 230026, China
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei 230088, China
| |
Collapse
|
60
|
Sumithra M, Sundaraganesan N, Rajesh R, Ilangovan V, Irfan A, Muthu S. Electron density, charge transfer, solvent effect and molecular spectroscopic studies on 2,2-Dimethyl-N-pyridin-4-yl-propionamide – A potential antioxidant. COMPUT THEOR CHEM 2023. [DOI: 10.1016/j.comptc.2023.114103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
|
61
|
Moore CC, Staroverov VN, Konermann L. Using Density Functional Theory for Testing the Robustness of Mobile-Proton Molecular Dynamics Simulations on Electrosprayed Ions: Structural Implications for Gaseous Proteins. J Phys Chem B 2023; 127:4061-4071. [PMID: 37116098 DOI: 10.1021/acs.jpcb.3c01581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/30/2023]
Abstract
Current experiments only provide low-resolution information on gaseous protein ions generated by electrospray ionization (ESI). Molecular dynamics (MD) simulations can yield complementary insights. Unfortunately, conventional MD does not capture the mobile nature of protons in gaseous proteins. Mobile-proton MD (MPMD) overcomes this limitation. Earlier MPMD data at 300 K indicated that protein ions generated by "native" ESI retain solution-like structures with a hydrophobic core and zwitterionic exterior [Bakhtiari, M.; Konermann, L. J. Phys. Chem. B 2019, 123, 1784-1796]. MPMD redistributes protons using electrostatic and proton affinity calculations. The robustness of this approach has never been scrutinized. Here, we close this gap by benchmarking MPMD against density functional theory (DFT) at the B3LYP/6-31G* level, which is well suited for predicting proton affinities. The computational cost of DFT necessitated the use of small peptides. The MPMD energetic ranking of proton configurations was found to be consistent with DFT single-point energies, implying that MPMD can reliably identify favorable protonation sites. Peptide MPMD runs converged to DFT-optimized structures only when applying 300-500 K temperature cycling, which was necessary to prevent trapping in local minima. Temperature cycling MPMD was then applied to gaseous protein ions. Native ubiquitin converted to slightly expanded structures with a zwitterionic core and a nonpolar exterior. Our data suggest that such inside-out protein structures are intrinsically preferred in the gas phase, and that they form in ESI experiments after moderate collisional excitation. This is in contrast to native ESI (with minimal collisional excitation, simulated by MPMD at 300 K), where kinetic trapping promotes the survival of solution-like structures. In summary, this work validates the MPMD approach for simulations on gaseous peptides and proteins.
Collapse
Affiliation(s)
- Conrad C Moore
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Viktor N Staroverov
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Lars Konermann
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
62
|
Ichibha T, Saritas K, Krogel JT, Luo Y, Kent PRC, Reboredo FA. Existence of La-site antisite defects in [Formula: see text] ([Formula: see text], Fe, and Co) predicted with many-body diffusion quantum Monte Carlo. Sci Rep 2023; 13:6703. [PMID: 37185382 PMCID: PMC10130183 DOI: 10.1038/s41598-023-33578-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2022] [Accepted: 04/14/2023] [Indexed: 05/17/2023] Open
Abstract
The properties of [Formula: see text] (M: 3d transition metal) perovskite crystals are significantly dependent on point defects, whether introduced accidentally or intentionally. The most studied defects in La-based perovskites are the oxygen vacancies and doping impurities on the La and M sites. Here, we identify that intrinsic antisite defects, the replacement of La by the transition metal, M, can be formed under M-rich and O-poor growth conditions, based on results of an accurate many-body ab initio approach. Our fixed-node diffusion Monte Carlo (FNDMC) calculations of [Formula: see text] ([Formula: see text], Fe, and Co) find that such antisite defects can have low formation energies and are magnetized. Complementary density functional theory (DFT)-based calculations show that Mn antisite defects in [Formula: see text] may cause the p-type electronic conductivity. These features could affect spintronics, redox catalysis, and other broad applications. Our bulk validation studies establish that FNDMC reproduces the antiferromagnetic state of [Formula: see text], whereas DFT with PBE (Perdew-Burke-Ernzerhof), SCAN (strongly constrained and appropriately normed), and the LDA+U (local density approximation with Coulomb U) functionals all favor ferromagnetic states, at variance with experiment.
Collapse
Affiliation(s)
- Tom Ichibha
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
- School of Information Science, Japan Advanced Institute of Science and Technology, Asahidai 1-1, Nomi, Ishikawa 923-1292 Japan
| | - Kayahan Saritas
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Jaron T. Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Ye Luo
- Computational Sciences Division, Argonne National Laboratory, Argonne, IL 60439 USA
| | - Paul R. C. Kent
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Fernando A. Reboredo
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| |
Collapse
|
63
|
Marques EA, De Gendt S, Pourtois G, van Setten MJ. Benchmarking First-Principles Reaction Equilibrium Composition Prediction. Molecules 2023; 28:molecules28093649. [PMID: 37175062 PMCID: PMC10179931 DOI: 10.3390/molecules28093649] [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: 03/31/2023] [Revised: 04/19/2023] [Accepted: 04/20/2023] [Indexed: 05/15/2023] Open
Abstract
The availability of thermochemical properties allows for the prediction of the equilibrium compositions of chemical reactions. The accurate prediction of these can be crucial for the design of new chemical synthesis routes. However, for new processes, these data are generally not completely available. A solution is the use of thermochemistry calculated from first-principles methods such as Density Functional Theory (DFT). Before this can be used reliably, it needs to be systematically benchmarked. Although various studies have examined the accuracy of DFT from an energetic point of view, few studies have considered its accuracy in predicting the temperature-dependent equilibrium composition. In this work, we collected 117 molecules for which experimental thermochemical data were available. From these, we constructed 2648 reactions. These experimentally constructed reactions were then benchmarked against DFT for 6 exchange-correlation functionals and 3 quality of basis sets. We show that, in reactions that do not show temperature dependence in the equilibrium composition below 1000 K, over 90% are predicted correctly. Temperature-dependent equilibrium compositions typically demonstrate correct qualitative behavior. Lastly, we show that the errors are equally caused by errors in the vibrational spectrum and the DFT electronic ground state energy.
Collapse
Affiliation(s)
- Esteban A Marques
- Department of Chemistry, KU Leuven (University of Leuven), Celestijnenlaan 200 F, 3001 Heverlee, Belgium
- IMEC, Kapeldreef 75, 3001 Leuven, Belgium
| | - Stefan De Gendt
- Department of Chemistry, KU Leuven (University of Leuven), Celestijnenlaan 200 F, 3001 Heverlee, Belgium
- IMEC, Kapeldreef 75, 3001 Leuven, Belgium
| | | | - Michiel J van Setten
- IMEC, Kapeldreef 75, 3001 Leuven, Belgium
- ETSF European Theoretical Spectroscopy Facility, Institut de Physique, Université de Liège, Allée du 6 août 17, 4000 Liège, Belgium
| |
Collapse
|
64
|
Ma Y, Zeng C, Zeng P, Hu Y, Li F, Zheng Z, Qin M, Lu X, Liu M. How Do Surface Polar Molecules Contribute to High Open-Circuit Voltage in Perovskite Solar Cells? ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2023:e2205072. [PMID: 37078797 DOI: 10.1002/advs.202205072] [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: 10/04/2022] [Revised: 12/20/2022] [Indexed: 05/03/2023]
Abstract
To date, the improvement of open-circuit voltage (VOC ) offers a breakthrough for the performance of perovskite solar cells (PSCs) toward their theoretical limit. Surface modification through organic ammonium halide salts (e.g., phenethylammonium ions PEA+ and phenmethylammonium ions PMA+ ) is one of the most straightforward strategies to suppress defect density, thereby leading to improved VOC . However, the mechanism underlying the high voltage remains unclear. Here, polar molecular PMA+ is applied at the interface between perovskite and hole transporting layer and a remarkably high VOC of 1.175 V is obtained which corresponds to an increase of over 100 mV in comparison to the control device. It is revealed that the equivalent passivation effect of surface dipole effectively improves the splitting of the hole quasi-Fermi level. Ultimately the combined effect of defect suppression and surface dipole equivalent passivation effect leads to an overall increase in significantly enhanced VOC . The resulted PSCs device reaches an efficiency of up to 24.10%. Contributions are identified here by the surface polar molecules to the high VOC in PSCs. A fundamental mechanism is suggested by use of polar molecules which enables further high voltage, leading ways to highly efficient perovskite-based solar cells.
Collapse
Affiliation(s)
- Yinyi Ma
- School of Materials and Energy, University of Electronic Science and Technology of China, Chengdu, 611731, P. R. China
| | - Chengsong Zeng
- School of Materials and Energy, University of Electronic Science and Technology of China, Chengdu, 611731, P. R. China
| | - Peng Zeng
- School of Materials and Energy, University of Electronic Science and Technology of China, Chengdu, 611731, P. R. China
| | - Yuchao Hu
- School of Materials and Energy, University of Electronic Science and Technology of China, Chengdu, 611731, P. R. China
| | - Faming Li
- School of Materials and Energy, University of Electronic Science and Technology of China, Chengdu, 611731, P. R. China
| | - Zhonghao Zheng
- School of Materials and Energy, University of Electronic Science and Technology of China, Chengdu, 611731, P. R. China
| | - Minchao Qin
- Department of Physics, The Chinese University of Hong Kong, Shatin, Hong Kong SAR, 999077, China
| | - Xinhui Lu
- Department of Physics, The Chinese University of Hong Kong, Shatin, Hong Kong SAR, 999077, China
| | - Mingzhen Liu
- School of Materials and Energy, University of Electronic Science and Technology of China, Chengdu, 611731, P. R. China
- State Key Laboratory Electronic Thin Film and Integrated Devices, University of Science and Technology of China, Chengdu, 611731, P. R. China
| |
Collapse
|
65
|
Nordholm S. Analysis of Bonding by Quantum Chemistry─Resolving Delocalization Stabilization in a Mechanistic Basis and New Hückel Model. J Phys Chem A 2023; 127:3449-3471. [PMID: 37040546 PMCID: PMC10123672 DOI: 10.1021/acs.jpca.2c08497] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/13/2023]
Abstract
No general and unique understanding of the mechanism of covalent bonding in physical terms is provided by current computational methods or by a consensus among experts. Bonding is studied by energy decomposition analysis but may also be related to the interatomic motion of valence electrons within the molecule. This dynamical view of the mechanism of bonding is not widely appreciated. The aim here is to make it accessible by translation into a corresponding form of quantum chemical energy analysis. The interatomic electron motion is directly related to the delocalization taking place when atomic basis functions are combined into molecular orbitals. A "tribasis method" is introduced, allowing an atomic basis set to form subsets of (1) strictly localized atomic functions and (2) interatomic bridge functions which allow delocalization. Calculations can then identify ground states without (no bridge functions) and with delocalization. The scheme is based on exact quantum mechanics but demonstrated by a minimal basis treatment of H2+ and H2 in Hartree-Fock and valence bond approximations which show that the bond energy is a sum of repulsive localization and more strongly attractive delocalization energies. The tribasis method is used to reconstruct the Hückel theory of π-electron delocalization in planar hydrocarbon molecules to account for the "overlap problem". In its empirically fitted form, the new theory can accurately resolve both π → π* transition energy and aromatic stabilization energy. The picture of covalent bonding emerging from both hydrogenic and Hückel calculations is that there is a presence of a Pauli repulsion of localization which is overcome by a roughly twice as strong delocalization stabilization to form the bond.
Collapse
Affiliation(s)
- Sture Nordholm
- Department of Chemistry and Molecular Biology, The University of Gothenburg, 412 96 Gothenburg, Sweden
| |
Collapse
|
66
|
Cytter Y, Nandy A, Duan C, Kulik HJ. Insights into the deviation from piecewise linearity in transition metal complexes from supervised machine learning models. Phys Chem Chem Phys 2023; 25:8103-8116. [PMID: 36876903 DOI: 10.1039/d3cp00258f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2023]
Abstract
Virtual high-throughput screening (VHTS) and machine learning (ML) with density functional theory (DFT) suffer from inaccuracies from the underlying density functional approximation (DFA). Many of these inaccuracies can be traced to the lack of derivative discontinuity that leads to a curvature in the energy with electron addition or removal. Over a dataset of nearly one thousand transition metal complexes typical of VHTS applications, we computed and analyzed the average curvature (i.e., deviation from piecewise linearity) for 23 density functional approximations spanning multiple rungs of "Jacob's ladder". While we observe the expected dependence of the curvatures on Hartree-Fock exchange, we note limited correlation of curvature values between different rungs of "Jacob's ladder". We train ML models (i.e., artificial neural networks or ANNs) to predict the curvature and the associated frontier orbital energies for each of these 23 functionals and then interpret differences in curvature among the different DFAs through analysis of the ML models. Notably, we observe spin to play a much more important role in determining the curvature of range-separated and double hybrids in comparison to semi-local functionals, explaining why curvature values are weakly correlated between these and other families of functionals. Over a space of 187.2k hypothetical compounds, we use our ANNs to pinpoint DFAs for which representative transition metal complexes have near-zero curvature with low uncertainty, demonstrating an approach to accelerate screening of complexes with targeted optical gaps.
Collapse
Affiliation(s)
- Yael Cytter
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| |
Collapse
|
67
|
Parambath S, Narayanan S J J, Parameswaran P. Five-membered N-heterocyclic beryllium(I) compounds: fluctuating electronic structures with ambiphilic reactivity. Dalton Trans 2023; 52:3378-3385. [PMID: 36810658 DOI: 10.1039/d2dt03263e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The structure, bonding, and reactivity of the five-membered N-heterocyclic beryllium compounds (NHBe), BeN2C2H4 (1) and BeN2(CH3)2C2H2 (2) were studied at the M06/def2-TZVPP//BP86/def2-TZVPP level of theory. The molecular orbital analysis indicates that NHBe is an aromatic 6π-electron system with an unoccupied σ-type spn-hybrid orbital on Be. Energy decomposition analysis combined with natural orbitals for chemical valence has been carried out with Be and L (L = N2C2H4 (1), N2(CH3)2C2H2 (2)) in their different electronic states as fragments at the BP86/TZ2P level of theory. The results indicate that the best bonding representation can be considered as an interaction between Be+ having the 2s02px12py02pz0 electronic configuration and L-. Accordingly, L- forms two donor-acceptor σ-bonds and one electron sharing π-bond with Be+. Compounds 1 and 2 show high proton and hydride affinity at beryllium, indicating its ambiphilic reactivity. The protonated structure results from adding a proton on the lone pair of electrons in the doubly excited state. On the other hand, the hydride adduct is formed by donating electrons from the hydride to an unoccupied σ-type spn-hybrid orbital on Be. These compounds show very high exothermic reaction energy for adduct formation with two electron donor ligands such as cAAC, CO, NHC, and PMe3.
Collapse
Affiliation(s)
- Sneha Parambath
- Department of Chemistry, National Institute of Technology, Calicut, Kerala, India.
| | - Jishnu Narayanan S J
- Department of Chemistry, National Institute of Technology, Calicut, Kerala, India.
| | | |
Collapse
|
68
|
DFT and TD-DFT study of hydrogen bonded complexes of aspartic acid and n water (n = 1 and 2). J Mol Model 2023; 29:94. [PMID: 36905452 DOI: 10.1007/s00894-023-05500-z] [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: 01/18/2023] [Accepted: 03/02/2023] [Indexed: 03/12/2023]
Abstract
CONTEXT Hydrogen bonds (HB) influence the conformational preferences of biomolecules and their optical and electronic properties. The directional interaction of molecules of water can be a prototype to understand the effects of HBs on biomolecules. Among the neurotransmitters (NT), L-aspartic acid (ASP) stands out due to its importance in health and as a precursor of several biomolecules. As it presents different functional groups and readily forms inter- and intramolecular HBs, ASP can be considered a prototype for understanding the behavior of NTs when interacting by HB with other substances. Although several theoretical studies have been performed in the past on isolated ASP and its formed complexes with water, both in gas and liquid phases, using DFT and TD-DFT formalisms, these works did not perform large basis set calculations or study electronic transitions of ASP-water complexes. We investigated the HB interactions in complexes of ASP and water molecules. The results show that the interactions between the carboxylic groups of ASP with water molecules, forming cyclic structures with two HBs, lead to more stable and less polar complexes than other conformers formed between water and the NH2 group. It was observed that there is a relationship between the deviation in the UV-Vis absorption band of the ASP and the interactions of water with the HOMO and LUMO orbitals with the stabilization/destabilization of the S1 state to the S0 of the complexes. However, in some cases, such as 1:1 complex ASP-W2, this analysis may be inaccurate due to small changes in ΔE. METHODS We studied the landscapes of the ground state surface of different conformers of isolated L-ASP and the L-ASP-(H2O)n complexes (n = 1 and 2) using the DFT formalism, with the B3LYP functional, and six different basis sets: 6-31 + + G(d,p), 6-311 + + G(d,p), D95 + + (d,p), D95V + + (d,p), cc-pVDZ, and, cc-pVTZ basis sets. The cc-pVTZ basis set provides the minimum energy of all conformers, and therefore, we performed the analysis with this basis set. We evaluated the stabilization of the ASP and complexes using the minimum ground state energy, corrected by the zero point energy and the interaction energy between the ASP and the water molecules. We also calculated the vertical electronic transitions S1 ← S0, and their properties using the TD-DFT formalism at B3LYP/cc-pVTZ level with the optimized geometries for S0 state with the same basis set. For the analysis of the vertical transitions of isolated ASP and the ASP-(H2O)n complexes, we calculated the electrostatic energy in the S0 and S1 states. We performed the calculations with the Gaussian 09 software package. We used the VMD software package to visualize the geometries and shapes of the molecule and complexes.
Collapse
|
69
|
Roy PO, Cuierrier E, Ernzerhof M. Generating Exchange-Correlation Functionals with a Simplified, Self-Consistent Correlation Factor Model. J Phys Chem A 2023; 127:2026-2033. [PMID: 36802604 DOI: 10.1021/acs.jpca.2c08397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2023]
Abstract
We focus on the spherically averaged exchange-correlation hole ρXC(r, u) of density functional theory, which describes the reduction in the electron density at a distance u due to the reference electron localized at position r. The correlation factor (CF) approach, where the model exchange hole ρXmodel(r, u) is multiplied by a CF (fC(r, u)) to yield an approximation to the exchange-correlation hole ρXC(r, u) = fC(r, u) ρXmodel(r, u), has proven to be a powerful tool for the development of new approximations. One of the remaining challenges within the CF approach is the self-consistent implementation of the resulting functionals. To address this issue, here we propose a simplification of the previously developed CFs such that self-consistent implementations become feasible. As an illustration of the simplified CF model, we develop a new meta-GGA functional, and using only a minimum of empiricism, we provide an easy derivation of an approximation that is of an accuracy similar to more involved meta-GGA functionals.
Collapse
Affiliation(s)
- Pierre-Olivier Roy
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| | - Etienne Cuierrier
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| | - Matthias Ernzerhof
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| |
Collapse
|
70
|
Huang P, Lu J, Jin L, Liu E, Li L. A DFT/TDDFT Investigation on Fluorescence and Electronic Properties of Chromone Derivatives. J Fluoresc 2023; 33:453-458. [PMID: 36441339 DOI: 10.1007/s10895-022-03095-x] [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: 10/13/2022] [Accepted: 11/18/2022] [Indexed: 11/29/2022]
Abstract
The development of quick and precise detection technologies for active compounds in vivo is critical for disease prevention, diagnosis and pathological investigation. The fluorescence signal of the fluorophore usually defines the probe's sensitivity to the chemical being examined. Many natural compounds containing flavone and isoflavone scaffolds exhibit a certain amount fluorescence, albeit with poor fluorescence quantum yields. Therefore, we used density functional theory (DFT) and time-dependent density functional theory (TDDFT) calculations to investigate the fluorescence characteristics of chromium-derived fluorophores in more depth. Different substituents are introduced at different positions of the chromone. As weak electron donor groups, alkyl and aromatic groups were discovered to have varying quantum yields on the fluorophore scaffold, and longer alkyl chains are favorable to enhance fluorescence quantum yield. In comparison to the amino group, substituted amino group can avoid group rotation, and the introduction of cyclic amines such as pyrrolidine and heterocyclic amines can improve optical characteristics. The electron-donating methoxy group at position 6 helps to increase the fluorescence quantum yield.
Collapse
Affiliation(s)
- Pei Huang
- Shaanxi Key Laboratory of Catalysis, College of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, 723001, China.
| | - Jiufu Lu
- Shaanxi Key Laboratory of Catalysis, College of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, 723001, China
| | - Lingxia Jin
- Shaanxi Key Laboratory of Catalysis, College of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, 723001, China
| | - Ernu Liu
- Shaanxi Key Laboratory of Catalysis, College of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, 723001, China
| | - Li Li
- Shaanxi Key Laboratory of Catalysis, College of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, 723001, China
| |
Collapse
|
71
|
An exotic 3-center/4-electron carbon–carbon pi long-bond: Is it tangible? Theor Chem Acc 2023. [DOI: 10.1007/s00214-023-02968-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
|
72
|
A density functional theory for ecology across scales. Nat Commun 2023; 14:1089. [PMID: 36841818 PMCID: PMC9968302 DOI: 10.1038/s41467-023-36628-4] [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: 06/16/2021] [Accepted: 02/09/2023] [Indexed: 02/27/2023] Open
Abstract
Ecology lacks a holistic approach that can model phenomena across temporal and spatial scales, largely because of the challenges in modelling systems with a large number of interacting constituents. This hampers our understanding of complex ecosystems and the impact that human interventions (e.g., deforestation, wildlife harvesting and climate change) have on them. Here we use density functional theory, a computational method for many-body problems in physics, to develop a computational framework for ecosystem modelling. Our methods accurately fit experimental and synthetic data of interacting multi-species communities across spatial scales and can project to unseen data. As the key concept we establish and validate a cost function that encodes the trade-offs between the various ecosystem components. We show how this single general modelling framework delivers predictions on par with established, but specialised, approaches for systems from predatory microbes to territorial flies to tropical tree communities. Our density functional framework thus provides a promising avenue for advancing our understanding of ecological systems.
Collapse
|
73
|
Li H, Jiao Y, Davey K, Qiao SZ. Data-Driven Machine Learning for Understanding Surface Structures of Heterogeneous Catalysts. Angew Chem Int Ed Engl 2023; 62:e202216383. [PMID: 36509704 DOI: 10.1002/anie.202216383] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Revised: 12/11/2022] [Accepted: 12/12/2022] [Indexed: 12/15/2022]
Abstract
The design of heterogeneous catalysts is necessarily surface-focused, generally achieved via optimization of adsorption energy and microkinetic modelling. A prerequisite is to ensure the adsorption energy is physically meaningful is the stable existence of the conceived active-site structure on the surface. The development of improved understanding of the catalyst surface, however, is challenging practically because of the complex nature of dynamic surface formation and evolution under in-situ reactions. We propose therefore data-driven machine-learning (ML) approaches as a solution. In this Minireview we summarize recent progress in using machine-learning to search and predict (meta)stable structures, assist operando simulation under reaction conditions and micro-environments, and critically analyze experimental characterization data. We conclude that ML will become the new norm to lower costs associated with discovery and design of optimal heterogeneous catalysts.
Collapse
Affiliation(s)
- Haobo Li
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, SA 5005, Australia
| | - Yan Jiao
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, SA 5005, Australia
| | - Kenneth Davey
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, SA 5005, Australia
| | - Shi-Zhang Qiao
- School of Chemical Engineering and Advanced Materials, The University of Adelaide, Adelaide, SA 5005, Australia
| |
Collapse
|
74
|
Bakht MA, Pooventhiran T, Thomas R, Kamal M, Din IU, Rehman NU, Ali I, Ajmal N, Ahsan MJ. Synthesis and Biological Evaluation of Octahydroquinazolinones as Phospholipase A2, and Protease Inhibitors: Experimental and Theoretical Exploration. Molecules 2023; 28:molecules28041944. [PMID: 36838935 PMCID: PMC9963251 DOI: 10.3390/molecules28041944] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 02/07/2023] [Accepted: 02/13/2023] [Indexed: 02/22/2023] Open
Abstract
Phospholipase A2 (PLA2) promotes inflammation via lipid mediators and releases arachidonic acid (AA), and these enzymes have been found to be elevated in a variety of diseases, including rheumatoid arthritis, sepsis, and atherosclerosis. The mobilization of AA by PLA2 and subsequent synthesis of prostaglandins are regarded as critical events in inflammation. Inflammatory processes may be treated with drugs that inhibit PLA2, thereby blocking the COX and LOX pathways in the AA cascade. To address this issue, we report herein an efficient method for the synthesis of a series of octahydroquinazolinone compounds (4a-h) in the presence of the catalyst Pd-HPW/SiO2 and their phospholipase A2, as well as protease inhibitory activities. Among eight compounds, two of them exhibited overwhelming results against PLA2 and protease. By using FT-IR, Raman, NMR, and mass spectroscopy, two novel compounds were thoroughly studied. After carefully examining the SAR of the investigated compounds against these enzymes, it was found that compounds (4a, 4b) containing both electron-donating and electron-withdrawing groups on the phenyl ring exhibited higher activity than compounds with only one of these groups. DFT studies were employed to study the electronic nature and reactivity properties of the molecules by optimizing at the BLYP/cc-pVDZ. Natural bond orbitals helped to study the various electron delocalizations in the molecules, and the frontier molecular orbitals helped with the reactivity and stability parameters. The nature and extent of the expressed biological activity of the molecule were studied using molecular docking with human non-pancreatic secretory phospholipase A2 (hnps-PLA2) (PDB ID: 1DB4) and protease K (PDB ID: 2PWB). The drug-ability of the molecule has been tested using ADMET, and pharmacodynamics data have been extracted. Both the compounds qualify for ADME properties and follow Lipinski's rule of five.
Collapse
Affiliation(s)
- Md. Afroz Bakht
- Department of Chemistry, College of Science and Humanity Studies, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
- Correspondence: (M.A.B.); (M.J.A.)
| | - Thangaiyan Pooventhiran
- Department of Chemistry, St Berchmans College (Autonomous), Mahatma Gandhi University, Changanassery 686101, Kerala, India
- Department of Mechanical Engineering, University Centre for Research & Development, Chandigarh University, Gharuan, Mohali 140413, Punjab, India
| | - Renjith Thomas
- Department of Chemistry, St Berchmans College (Autonomous), Mahatma Gandhi University, Changanassery 686101, Kerala, India
| | - Mehnaz Kamal
- Department of Pharmaceutical Chemistry, College of Pharmacy, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
| | - Israf Ud Din
- Department of Chemistry, College of Science and Humanity Studies, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
| | - Najeeb Ur Rehman
- Department of Pharmacology & Toxicology, College of Pharmacy, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
| | - Imtiaz Ali
- Preparatory College, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
| | - Noushin Ajmal
- Department of Basic Science and Humanities, Pratap University, Jaipur 303104, Rajasthan, India
| | - Mohamed Jawed Ahsan
- Department of Pharmaceutical Chemistry, Maharishi Arvind College of Pharmacy, Ambabari Circle, Jaipur 302039, Rajasthan, India
- Correspondence: (M.A.B.); (M.J.A.)
| |
Collapse
|
75
|
Lahiri N, Song D, Zhang X, Huang X, Stoerzinger KA, Carvalho OQ, Adiga PP, Blum M, Rosso KM. Interplay between Facets and Defects during the Dissociative and Molecular Adsorption of Water on Metal Oxide Surfaces. J Am Chem Soc 2023; 145:2930-2940. [PMID: 36696237 DOI: 10.1021/jacs.2c11291] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Surface terminations and defects play a central role in determining how water interacts with metal oxides, thereby setting important properties of the interface that govern reactivity such as the type and distribution of hydroxyl groups. However, the interconnections between facets and defects remain poorly understood. This limits the usefulness of conventional notions such as that hydroxylation is controlled by metal cation exposure at the surface. Here, using hematite (α-Fe2O3) as a model system, we show how oxygen vacancies overwhelm surface cation-dependent hydroxylation behavior. Synchrotron-based ambient-pressure X-ray photoelectron spectroscopy was used to monitor the adsorption of molecular water and its dissociation to form hydroxyl groups in situ on (001), (012), or (104) facet-engineered hematite nanoparticles. Supported by density functional theory calculations of the respective surface energies and oxygen vacancy formation energies, the findings show how oxygen vacancies are more prone to form on higher energy facets and induce surface hydroxylation at extremely low relative humidity values of 5 × 10-5%. When these vacancies are eliminated, the extent of surface hydroxylation across the facets is as expected from the areal density of exposed iron cations at the surface. These findings help answer fundamental questions about the nature of reducible metal oxide-water interfaces in natural and technological settings and lay the groundwork for rational design of improved oxide-based catalysts.
Collapse
Affiliation(s)
- Nabajit Lahiri
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington99352, United States
| | - Duo Song
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington99352, United States
| | - Xin Zhang
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington99352, United States
| | - Xiaopeng Huang
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington99352, United States
| | - Kelsey A Stoerzinger
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington99352, United States.,Department of Chemical, Biological and Environmental Engineering, Oregon State University, Corvallis, Oregon97331, United States
| | - O Quinn Carvalho
- Department of Chemical, Biological and Environmental Engineering, Oregon State University, Corvallis, Oregon97331, United States
| | - Prajwal P Adiga
- Department of Chemical, Biological and Environmental Engineering, Oregon State University, Corvallis, Oregon97331, United States
| | - Monika Blum
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Kevin M Rosso
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington99352, United States
| |
Collapse
|
76
|
Sun X, Hansen T, Poater J, Hamlin TA, Bickelhaupt FM. Rational design of iron catalysts for C-X bond activation. J Comput Chem 2023; 44:495-505. [PMID: 35137432 PMCID: PMC10078697 DOI: 10.1002/jcc.26818] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2021] [Revised: 01/16/2022] [Accepted: 01/20/2022] [Indexed: 01/10/2023]
Abstract
We have quantum chemically studied the iron-mediated CX bond activation (X = H, Cl, CH3 ) by d8 -FeL4 complexes using relativistic density functional theory at ZORA-OPBE/TZ2P. We find that by either modulating the electronic effects of a generic iron-catalyst by a set of ligands, that is, CO, BF, PH3 , BN(CH3 )2 , or by manipulating structural effects through the introduction of bidentate ligands, that is, PH2 (CH2 )n PH2 with n = 6-1, one can significantly decrease the reaction barrier for the CX bond activation. The combination of both tuning handles causes a decrease of the CH activation barrier from 10.4 to 4.6 kcal mol-1 . Our activation strain and Kohn-Sham molecular orbital analyses reveal that the electronic tuning works via optimizing the catalyst-substrate interaction by introducing a strong second backdonation interaction (i.e., "ligand-assisted" interaction), while the mechanism for structural tuning is mainly caused by the reduction of the required activation strain because of the pre-distortion of the catalyst. In all, we present design principles for iron-based catalysts that mimic the favorable behavior of their well-known palladium analogs in the bond-activation step of cross-coupling reactions.
Collapse
Affiliation(s)
- Xiaobo Sun
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Departament de Química Inorgànica i Orgànica & IQTCUB, Universitat de Barcelona, Barcelona, Spain
| | - Thomas Hansen
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Departament de Química Inorgànica i Orgànica & IQTCUB, Universitat de Barcelona, Barcelona, Spain.,Leiden Institute of Chemistry, Leiden University, Leiden, The Netherlands
| | - Jordi Poater
- Departament de Química Inorgànica i Orgànica & IQTCUB, Universitat de Barcelona, Barcelona, Spain.,ICREA, Barcelona, Spain
| | - Trevor A Hamlin
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Friedrich Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Institute for Molecules and Materials (IMM), Radboud University, Nijmegen, The Netherlands
| |
Collapse
|
77
|
Monte MBDM, Pimentel DA, Albuquerque MDDFD, Neumann R, Silva LA, Correia JC, Uliana A. Synergism of mixed cationic collectors in the flotation of quartz unveiled by AFM, solution chemistry and quantum chemical calculations. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
|
78
|
First-principles calculations of hematite (α-Fe 2O 3) by self-consistent DFT+U+V. iScience 2023; 26:106033. [PMID: 36824287 PMCID: PMC9941207 DOI: 10.1016/j.isci.2023.106033] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Revised: 09/22/2022] [Accepted: 01/18/2023] [Indexed: 01/25/2023] Open
Abstract
Owing to the confined Fe-3d orbitals and self-interaction error of exchange-correlation functionals, approximate DFT fails to describe iron oxides electronic structure and magnetic properties accurately. Hybrid DFT or DFT + U can solve these problems, but the former is expensive, and the latter only considers on-site interactions. Here, we used DFT + U + V, a DFT + U extension including inter-site interactions, to simulate the structural, magnetic, and electronic properties, along with Fe and O K-edge XAS spectra of α-Fe2O3. Two types of atomic orbital projectors were studied, orthogonalized and non-orthogonalized. DFT + U + V improves the description of the structural, magnetic, and electronic properties of α-Fe2O3 compared to approximate DFT. The accuracy of the correction depends on the orbital projector used. DFT + U + V with orthogonalized projectors achieves the best experimental agreement at a fraction of hybrid DFT cost. This work emphasizes the importance of inter-site interactions and the type of atomic orbital projectors used in the theoretical research of α-Fe2O3.
Collapse
|
79
|
Molecular engineering on D-π-A organic dyes with flavone-based different acceptors for highly efficient dye-sensitized solar cells using experimental and computational study. J Mol Model 2023; 29:45. [PMID: 36656395 DOI: 10.1007/s00894-023-05445-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Accepted: 01/04/2023] [Indexed: 01/20/2023]
Abstract
CONTEXT The improvement of new organic flavone-based donor-spacer-acceptor (D-π-A) type dye molecules of the 3-(4-hydroxypiperidin-2-yloxy)-7-hydroxy-2-(3,4-dihydroxyphenyl)-4H-chromen-4-one (D1), 7-hydroxy-2-(3,4-dihydroxyphenyl)-3-(piperidin-4-yloxy)-4H-chromen-4-one (D2), and 3-((2-aminopyridin-4-yloxy)methoxy)-7-hydroxy-2-(3,4-dihydroxyphenyl)-4H-chromen-4-one (D3) were successfully designed and synthesized for dye-sensitized solar cells (DSSCs). METHODS Here, we discuss the synthesis of flavone compounds as well as their photophysical and electrochemical characterization. Using the Gaussian 09w software, the electronic structures and apsorption spectra have been calculated at the B3LYP, B3PW91, CAM-B3LYP, MPW1PW91, PBEPBE, and ωB97XD theory with the 6-311G(d,p) basis sets. RESULTS The computed values of the D2 molecule ground state optimized HOMOs-LUMOs energy is well positioned for advantageous charge transfer (CT) into the semiconducting material (TiO2) as well as the electron injection process. With a high power conversion efficiency (PCE) of 3.46% (VOC = 0.718 V, JSC = 7.07 mA cm-2, and FF = 0.68), the D2 compound also demonstrated good photovoltaic (PV) properties. CONCLUSION These findings unequivocally demonstrate that altering the D-π-A metal-free organic material electron-withdrawing capacity is a useful strategy for enhancing the optical and electrical characteristics of the organic PV system.
Collapse
|
80
|
Fürst S, Haasler M, Grotjahn R, Kaupp M. Full Implementation, Optimization, and Evaluation of a Range-Separated Local Hybrid Functional with Wide Accuracy for Ground and Excited States. J Chem Theory Comput 2023; 19:488-502. [PMID: 36625881 DOI: 10.1021/acs.jctc.2c00782] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
We report the first full and efficient implementation of range-separated local hybrid functionals (RSLHs) into the TURBOMOLE program package. This enables the computation of ground-state energies and nuclear gradients as well as excitation energies. Regarding the computational effort, RSLHs scale like regular local hybrid functionals (LHs) with system or basis set size and increase timings by a factor of 2-3 in total. An advanced RSLH, ωLH22t, has been optimized for atomization energies and reaction barriers. It is an extension of the recent LH20t local hybrid and is based on short-range PBE and long-range HF exchange-energy densities, a pig2 calibration function to deal with the gauge ambiguity of exchange-energy densities, and reoptimized B95c correlation. ωLH22t has been evaluated for a wide range of ground-state and excited-state quantities. It further improves upon the already successful LH20t functional for the GMTKN55 main-group energetics test suite, and it outperforms any global hybrid while performing close to the top rung-4 functional, ωB97M-V, for these evaluations when augmented by D4 dispersion corrections. ωLH22t performs excellently for transition-metal reactivity and provides good balance between delocalization errors and left-right correlation for mixed-valence systems, with a somewhat larger bias toward localized states compared to LH20t. It approaches the accuracy of the best local hybrids to date for core, valence singlet and triplet, and Rydberg excitation energies while improving strikingly on intra- and intermolecular charge-transfer excitations, comparable to the most successful range-separated hybrids available.
Collapse
Affiliation(s)
- Susanne Fürst
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Matthias Haasler
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Robin Grotjahn
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
81
|
Alkhodier RA, Mishra SK, Doerksen RJ, Colby DA. Comparison of Conformational Analyses of Naturally Occurring Flavonoid- O-Glycosides with Unnatural Flavonoid-CF 2-Glycosides Using Molecular Modeling. J Chem Inf Model 2023; 63:375-386. [PMID: 36512328 PMCID: PMC9904208 DOI: 10.1021/acs.jcim.2c01147] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Many glycosylated natural products display biological activity and are deglycosylated by the metabolic processes of the body. Although unnatural CF2-glycosides have been proposed as nonhydrolyzable analogues, CF2-derivatives of natural products are exceedingly challenging to synthesize and few examples exist. These difluorinated molecules may have unique conformational behavior as a consequence of changing the glycosidic linkage. In this study, we performed conformational searches using MacroModel followed by molecular dynamics simulations to investigate the conformational behavior of the glycosidic bonds in flavonoid-O-glycosides and in corresponding CF2-glycosylated derivatives. Compared to their O-glycosylated analogues, flavonoid-3-CF2-glycosides and flavonoid-5-CF2-glycosides showed conformational bias, whereas flavonoid-7-CF2-glycosides showed more flexibility. Flavonoid-5-CF2-glycosides were the least flexible compared to all others. Our results show that the site of the glycosylation and the substitution pattern on the flavonoid determine the conformational properties of these molecules. These two factors influence the steric destabilization and/or stereoelectronic stabilization which govern the conformational behavior of the flavonoid glycosides. Moreover, a docking study of quercitrin and its CF2-analogue into murine ribosomal kinase RSK2 demonstrated the potential for flavonoid-CF2-glycosides to retain a similar binding pose as the parent O-glycoside. These findings will assist in designing stable flavonoid-CF2-glycosides for carbohydrate research.
Collapse
Affiliation(s)
- Reem A Alkhodier
- Department of BioMolecular Sciences, University of Mississippi, University, MS 38677
| | - Sushil K. Mishra
- Department of BioMolecular Sciences, University of Mississippi, University, MS 38677
| | - Robert J. Doerksen
- Department of BioMolecular Sciences, University of Mississippi, University, MS 38677
| | - David A. Colby
- Department of BioMolecular Sciences, University of Mississippi, University, MS 38677
| |
Collapse
|
82
|
Brocławik E, Fidrysiak M, Hendzel M, Spałek J. Interparticle correlations and chemical bonding from physical side: Covalency vs atomicity and ionicity. ADVANCES IN QUANTUM CHEMISTRY 2023. [DOI: 10.1016/bs.aiq.2023.02.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/18/2023]
|
83
|
Liu N, Li R, Zhu J, Liu Q, Chen R, Yu J, Li Y, Zhang H, Wang J. Z-scheme heterojunction ZnS/WO3 composite: Photocatalytic reduction of uranium and band gap regulation mechanism. J Colloid Interface Sci 2023; 630:727-737. [DOI: 10.1016/j.jcis.2022.10.151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Revised: 10/22/2022] [Accepted: 10/29/2022] [Indexed: 11/08/2022]
|
84
|
Hu H, Tao Y, Wang D, Li C, Jiang Q, Shi Y, Wang J, Qin J, Zhou S, Kong Y. Rational modification of hydroxy-functionalized covalent organic frameworks for enhanced photocatalytic hydrogen peroxide evolution. J Colloid Interface Sci 2023; 629:750-762. [PMID: 36193619 DOI: 10.1016/j.jcis.2022.09.111] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 09/12/2022] [Accepted: 09/20/2022] [Indexed: 11/16/2022]
Abstract
Covalent organic frameworks (COFs), a class of flexibly tunable crystalline materials, have fascinating potential in photocatalytic hydrogen peroxide (H2O2) evolution under visible light irradiation. However, achieving efficient catalytic activity by tuning the composition of COFs and the linkages of building blocks is still a challenge. Herein, four imine-linked COFs with different numbers of hydroxy-functionalized are constructed to unveil the latent structure-activity relationship between the reversibility of bonding in supramolecular chemistry and the photocatalytic H2O2 performance. As the optimized material, TAPT-HTA-COF (1H-COF) containing single hydroxy group in aldehyde node exhibits a highest ordered structure and conjugation degree along and across the plane in the extended frameworks originating from the flexibly reversible iminol-to-ketoenamine tautomerism than others, which broadens the visible light absorption and accelerates the dissociation of photogenerated carriers in 1H-COF. These merits ensure that 1H-COF has the highest H2O2 yield (44.5 μmol L-1) and O2 two-electron reduction pathway among the four COFs under visible light irradiation (λ > 420 nm, 10 vol% isopropanol aqueous solution). At the same time, the long-range ordered framework of 1H-COF is well preserved during the photocatalytic H2O2 evolution process assisted by the proton-induced tautomerization. This work facilitates the design and development of COF-based photocatalysts in the evolution of H2O2.
Collapse
Affiliation(s)
- Hao Hu
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China
| | - Yinglong Tao
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China
| | - Di Wang
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China
| | - Changlai Li
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China
| | - Qichao Jiang
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China
| | - Yuexin Shi
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China
| | - Jian Wang
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China
| | - Jinping Qin
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China
| | - Shijian Zhou
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China.
| | - Yan Kong
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical ·Engineering, Nanjing Tech University, Nanjing 210009, China.
| |
Collapse
|
85
|
Qi Y, Liu Y, Zhang B, Wang M, Cao L, Song L, Jin N, Zhang H. Comparative antibacterial analysis of the anthraquinone compounds based on the AIM theory, molecular docking, and dynamics simulation analysis. J Mol Model 2022; 29:16. [PMID: 36547745 DOI: 10.1007/s00894-022-05406-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2022] [Accepted: 12/04/2022] [Indexed: 12/24/2022]
Abstract
BACKGROUND Hydroxyanthraquinones and anthraquinone glucoside derivatives are always considered as the active antibacterial components. METHODS Comparison of structure characteristics and antibacterial effect of these compounds was performed by applying quantum chemical calculations, atoms in molecules theory, molecular docking, and dynamics simulation procedure. Density functional theory calculation with B3LYP using 6-31G (d, p) basis set has been used to determine ground state molecular geometries. RESULTS The molecular geometric stability, electrostatic potential, frontier orbital energies, and topological properties were analyzed at the active site. Once glucose ring is introduced into the hydroxyanthraquinone rings, almost all of the positive molecular potentials are distributed among the hydroxyl hydrogen atoms of the glucose rings. In addition, low electron density ρ (r) and positive Laplacian value of the O-H bond of the anthraquinone glucoside are the evidences of the highly polarized and covalently decreased bonding interactions. The anthraquinone glucoside compounds have generally higher intermolecular binding energies than the corresponding aglycones due to the strong interaction between the glucose rings and the surrounding amino acids. Molecular dynamics simulations further explored the stability and dynamic behavior of the anthraquinone compound and protein complexes through RMSD, RMSF, SASA, and Rg. CONCLUSION The type of carboxyl, hydroxyl, and hydroxymethyl groups on phenyl ring and the substituent glucose rings is important to the interactions with the topoisomerase type II enzyme DNA gyrase B.
Collapse
Affiliation(s)
- Yanjiao Qi
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, 730000, Lanzhou, People's Republic of China.,Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in Universities of Gansu Province, 730000, Lanzhou, People's Republic of China
| | - Yue Liu
- Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in Universities of Gansu Province, 730000, Lanzhou, People's Republic of China
| | - Bo Zhang
- Key Laboratory for Utility of Environment-Friendly Composite Materials and Biomass in Universities of Gansu Province, 730000, Lanzhou, People's Republic of China
| | - Mingyang Wang
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, 730000, Lanzhou, People's Republic of China
| | - Long Cao
- Department of Chemical Engineering, Northwest Minzu University, 730124, Lanzhou, People's Republic of China
| | - Li Song
- Gansu Hualing Dairy Co. LTD., Lanzhou, People's Republic of China
| | - Nengzhi Jin
- Gansu Province Computing Center, 730000, Lanzhou, People's Republic of China
| | - Hong Zhang
- Key Laboratory of Environment-Friendly Composite Materials of the State Ethnic Affairs Commission, 730000, Lanzhou, People's Republic of China.
| |
Collapse
|
86
|
Hesselmann A, Werner HJ, Knowles PJ. Thermochemical evaluation of adaptive and fixed density functional theory quadrature schemes. J Chem Phys 2022; 157:234106. [PMID: 36550055 DOI: 10.1063/5.0119622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
A systematic study is made of the accuracy and efficiency of a number of existing quadrature schemes for molecular Kohn-Sham Density-Functional Theory (DFT) using 408 molecules and 254 chemical reactions. Included are the fixed SG-x (x = 0-3) grids of Gill et al., Dasgupta, and Herbert, the 3-zone grids of Treutler and Ahlrichs, a fixed five-zone grid implemented in Molpro, and a new adaptive grid scheme. While all methods provide a systematic reduction of errors upon extension of the grid sizes, significant differences are observed in the accuracies for similar grid sizes with various approaches. For the tests in this work, the SG-x fixed grids are less suitable to achieve high accuracies in the DFT integration, while our new adaptive grid performed best among the schemes studied in this work. The extra computational time to generate the adaptive grid scales linearly with molecular size and is negligible compared with the time needed for the self-consistent field iterations for large molecules. A comparison of the grid accuracies using various density functionals shows that meta-GGA functionals need larger integration grids than GGA functionals to reach the same degree of accuracy, confirming previous investigations of the numerical stability of meta-GGA functionals. On the other hand, the grid integration errors are almost independent of the basis set, and the basis set errors are mostly much larger than the errors caused by the numerical integrations, even when using the smallest grids tested in this work.
Collapse
Affiliation(s)
- Andreas Hesselmann
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
87
|
Fedorov DG. Parametrized quantum-mechanical approaches combined with the fragment molecular orbital method. J Chem Phys 2022; 157:231001. [PMID: 36550057 DOI: 10.1063/5.0131256] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022] Open
Abstract
Fast parameterized methods such as density-functional tight-binding (DFTB) facilitate realistic calculations of large molecular systems, which can be accelerated by the fragment molecular orbital (FMO) method. Fragmentation facilitates interaction analyses between functional parts of molecular systems. In addition to DFTB, other parameterized methods combined with FMO are also described. Applications of FMO methods to biochemical and inorganic systems are reviewed.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
88
|
Lu Y, Gao J. Fundamental Variable and Density Representation in Multistate DFT for Excited States. J Chem Theory Comput 2022; 18:7403-7411. [PMID: 36346908 DOI: 10.1021/acs.jctc.2c00859] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Complementary to the theorems of Hohenberg and Kohn for the ground state, Theophilou's subspace theory establishes a one-to-one relationship between the total eigenstate energy and density ρV(r) of the subspace spanned by the lowest N eigenstates. However, the individual eigenstate energies are not directly available from such a subspace density functional theory. Lu and Gao (J. Phys. Chem. Lett. 2022, 13, 7762) recently proved that the Hamiltonian projected on to this subspace is a matrix functional H[D] of the multistate matrix density D(r) and that variational optimization of the trace of the Hamiltonian matrix functional yields exactly the individual eigenstates and densities. This study shows that the matrix density D(r) is the necessary fundamental variable in order to determine the exact energies and densities of the individual eigenstates. Furthermore, two ways of representing the matrix density are introduced, making use of nonorthogonal and orthogonal orbitals. In both representations, a multistate active space of auxiliary states can be constructed to exactly represent D(r) with which an explicit formulation of the Hamiltonian matrix functional H[D] is presented. Importantly, the use of a common set of orthonormal orbitals makes it possible to carry out multistate self-consistent-field optimization of the auxiliary states with singly and doubly excited configurations (MS-SDSCF).
Collapse
Affiliation(s)
- Yangyi Lu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen518055, China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen518055, China.,Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota55455, United States
| |
Collapse
|
89
|
Teale AM, Helgaker T, Savin A, Adamo C, Aradi B, Arbuznikov AV, Ayers PW, Baerends EJ, Barone V, Calaminici P, Cancès E, Carter EA, Chattaraj PK, Chermette H, Ciofini I, Crawford TD, De Proft F, Dobson JF, Draxl C, Frauenheim T, Fromager E, Fuentealba P, Gagliardi L, Galli G, Gao J, Geerlings P, Gidopoulos N, Gill PMW, Gori-Giorgi P, Görling A, Gould T, Grimme S, Gritsenko O, Jensen HJA, Johnson ER, Jones RO, Kaupp M, Köster AM, Kronik L, Krylov AI, Kvaal S, Laestadius A, Levy M, Lewin M, Liu S, Loos PF, Maitra NT, Neese F, Perdew JP, Pernal K, Pernot P, Piecuch P, Rebolini E, Reining L, Romaniello P, Ruzsinszky A, Salahub DR, Scheffler M, Schwerdtfeger P, Staroverov VN, Sun J, Tellgren E, Tozer DJ, Trickey SB, Ullrich CA, Vela A, Vignale G, Wesolowski TA, Xu X, Yang W. DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science. Phys Chem Chem Phys 2022; 24:28700-28781. [PMID: 36269074 PMCID: PMC9728646 DOI: 10.1039/d2cp02827a] [Citation(s) in RCA: 62] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 08/09/2022] [Indexed: 12/13/2022]
Abstract
In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.
Collapse
Affiliation(s)
- Andrew M. Teale
- School of Chemistry, University of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andreas Savin
- Laboratoire de Chimie Théorique, CNRS and Sorbonne University, 4 Place Jussieu, CEDEX 05, 75252 Paris, France.
| | - Carlo Adamo
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany.
| | - Alexei V. Arbuznikov
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7Straße des 17. Juni 13510623Berlin
| | | | - Evert Jan Baerends
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56125 Pisa, Italy.
| | - Patrizia Calaminici
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Eric Cancès
- CERMICS, Ecole des Ponts and Inria Paris, 6 Avenue Blaise Pascal, 77455 Marne-la-Vallée, France.
| | - Emily A. Carter
- Department of Mechanical and Aerospace Engineering and the Andlinger Center for Energy and the Environment, Princeton UniversityPrincetonNJ 08544-5263USA
| | | | - Henry Chermette
- Institut Sciences Analytiques, Université Claude Bernard Lyon1, CNRS UMR 5280, 69622 Villeurbanne, France.
| | - Ilaria Ciofini
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - T. Daniel Crawford
- Department of Chemistry, Virginia TechBlacksburgVA 24061USA,Molecular Sciences Software InstituteBlacksburgVA 24060USA
| | - Frank De Proft
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | | | - Claudia Draxl
- Institut für Physik and IRIS Adlershof, Humboldt-Universität zu Berlin, 12489 Berlin, Germany. .,Fritz-Haber-Institut der Max-Planck-Gesellschaft, 14195 Berlin, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany. .,Beijing Computational Science Research Center (CSRC), 100193 Beijing, China.,Shenzhen JL Computational Science and Applied Research Institute, 518110 Shenzhen, China
| | - Emmanuel Fromager
- Laboratoire de Chimie Quantique, Institut de Chimie, CNRS/Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France.
| | - Patricio Fuentealba
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile.
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute, and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, USA.
| | - Giulia Galli
- Pritzker School of Molecular Engineering and Department of Chemistry, The University of Chicago, Chicago, IL, USA.
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China. .,Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA
| | - Paul Geerlings
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | - Nikitas Gidopoulos
- Department of Physics, Durham University, South Road, Durham DH1 3LE, UK.
| | - Peter M. W. Gill
- School of Chemistry, University of SydneyCamperdown NSW 2006Australia
| | - Paola Gori-Giorgi
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Andreas Görling
- Chair of Theoretical Chemistry, University of Erlangen-Nuremberg, Egerlandstrasse 3, 91058 Erlangen, Germany.
| | - Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University, Gold Coast, Qld 4222, Australia.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany.
| | - Oleg Gritsenko
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Hans Jørgen Aagaard Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark.
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie UniversityHalifaxNova ScotiaB3H 4R2Canada
| | - Robert O. Jones
- Peter Grünberg Institut PGI-1, Forschungszentrum Jülich52425 JülichGermany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, 10623, Berlin.
| | - Andreas M. Köster
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav)CDMX07360Mexico
| | - Leeor Kronik
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth, 76100, Israel.
| | - Anna I. Krylov
- Department of Chemistry, University of Southern CaliforniaLos AngelesCalifornia 90089USA
| | - Simen Kvaal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andre Laestadius
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Mel Levy
- Department of Chemistry, Tulane University, New Orleans, Louisiana, 70118, USA.
| | - Mathieu Lewin
- CNRS & CEREMADE, Université Paris-Dauphine, PSL Research University, Place de Lattre de Tassigny, 75016 Paris, France.
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC 27599-3420, USA. .,Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599-3290, USA
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France.
| | - Neepa T. Maitra
- Department of Physics, Rutgers University at Newark101 Warren StreetNewarkNJ 07102USA
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany.
| | - John P. Perdew
- Departments of Physics and Chemistry, Temple UniversityPhiladelphiaPA 19122USA
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.
| | - Pascal Pernot
- Institut de Chimie Physique, UMR8000, CNRS and Université Paris-Saclay, Bât. 349, Campus d'Orsay, 91405 Orsay, France.
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA. .,Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA
| | - Elisa Rebolini
- Institut Laue Langevin, 71 avenue des Martyrs, 38000 Grenoble, France.
| | - Lucia Reining
- Laboratoire des Solides Irradiés, CNRS, CEA/DRF/IRAMIS, École Polytechnique, Institut Polytechnique de Paris, F-91120 Palaiseau, France. .,European Theoretical Spectroscopy Facility
| | - Pina Romaniello
- Laboratoire de Physique Théorique (UMR 5152), Université de Toulouse, CNRS, UPS, France.
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA.
| | - Dennis R. Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS – Centre for Molecular Simulation, IQST – Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary2500 University Drive NWCalgaryAlbertaT2N 1N4Canada
| | - Matthias Scheffler
- The NOMAD Laboratory at the FHI of the Max-Planck-Gesellschaft and IRIS-Adlershof of the Humboldt-Universität zu Berlin, Faradayweg 4-6, D-14195, Germany.
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632 Auckland, New Zealand.
| | - Viktor N. Staroverov
- Department of Chemistry, The University of Western OntarioLondonOntario N6A 5B7Canada
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA.
| | - Erik Tellgren
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - David J. Tozer
- Department of Chemistry, Durham UniversitySouth RoadDurhamDH1 3LEUK
| | - Samuel B. Trickey
- Quantum Theory Project, Deptartment of Physics, University of FloridaGainesvilleFL 32611USA
| | - Carsten A. Ullrich
- Department of Physics and Astronomy, University of MissouriColumbiaMO 65211USA
| | - Alberto Vela
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Giovanni Vignale
- Department of Physics, University of Missouri, Columbia, MO 65203, USA.
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, Université de Genève30 Quai Ernest-Ansermet1211 GenèveSwitzerland
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovation Materials, Collaborative Innovation Centre of Chemistry for Energy Materials, MOE Laboratory for Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, China.
| | - Weitao Yang
- Department of Chemistry and Physics, Duke University, Durham, NC 27516, USA.
| |
Collapse
|
90
|
Cuierrier E, Roy PO, Wang R, Ernzerhof M. The fourth-order expansion of the exchange hole and neural networks to construct exchange–correlation functionals. J Chem Phys 2022; 157:171103. [DOI: 10.1063/5.0122761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
The curvature Q σ of spherically averaged exchange (X) holes ρX, σ(r, u) is one of the crucial variables for the construction of approximations to the exchange–correlation energy of Kohn–Sham theory, the most prominent example being the Becke–Roussel model [A. D. Becke and M. R. Roussel, Phys. Rev. A 39, 3761 (1989)]. Here, we consider the next higher nonzero derivative of the spherically averaged X hole, the fourth-order term T σ. This variable contains information about the nonlocality of the X hole and we employ it to approximate hybrid functionals, eliminating the sometimes demanding calculation of the exact X energy. The new functional is constructed using machine learning; having identified a physical correlation between T σ and the nonlocality of the X hole, we employ a neural network to express this relation. While we only modify the X functional of the Perdew–Burke–Ernzerhof functional [Perdew et al., Phys. Rev. Lett. 77, 3865 (1996)], a significant improvement over this method is achieved.
Collapse
Affiliation(s)
- Etienne Cuierrier
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| | - Pierre-Olivier Roy
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| | - Rodrigo Wang
- Good Chemistry Company, Vancouver, British Columbia V6E 4B1, Canada
| | - Matthias Ernzerhof
- Département de Chimie, Université de Montréal, C.P. 6128 Succursale A, Montréal, Québec H3C 3J7, Canada
| |
Collapse
|
91
|
Synthesis and Spectroscopic Characterization of Selected Phenothiazines and Phenazines Rationalized Based on DFT Calculation. Molecules 2022; 27:molecules27217519. [PMID: 36364378 PMCID: PMC9653876 DOI: 10.3390/molecules27217519] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Revised: 10/26/2022] [Accepted: 10/28/2022] [Indexed: 11/06/2022] Open
Abstract
Two unique structures were isolated from the phosphorylation reaction of 10H-phenothiazine. The 5,5-dimethyl-2-(10H-phenothiazin-10-yl)-1,3,2-dioxaphosphinane 2-oxide (2a) illustrates the product of N-phosphorylation of phenothiazine. Moreover, a potential product of 2a instability, a thiophosphoric acid 2b, was successfully isolated and structurally characterized. Molecule 2a, similarly to sulfoxide derivative 3, possesses interesting phosphorescence properties due to the presence of d-pπ bonds. The X-ray, NMR, and DFT computational studies indicate that compound 2a exhibits an anomeric effect. Additionally, the syntheses of selected symmetrical and unsymmetrical pyridine-embedded phenazines were elaborated. To compare the influence of phosphorus and sulfur atoms on the structural characteristics of 10H-phenothiazine derivatives, the high-quality crystals of (4a,12a-dihydro-12H-benzo[5,6][1,4]thiazino[2,3-b]quinoxalin-12-yl)(phenyl)methanone (1) and selected phenazines 5,12-diisopropyl-3,10-dimethyldipyrido[3,2-a:3′,2′-h]phenazine (5) and 5-isopropyl-N,N,3-trimethylpyrido[3,2-a]phenazin-10-amine (6a) were obtained. The structures of molecules 1, 2a, 2-mercapto-5,5-dimethyl-1,3,2-dioxaphosphinane 2-oxide (2b), 3,7-dinitro-10H-phenothiazine 5-oxide (3), 5 and 6a were determined by single-crystal X-ray diffraction measurements.
Collapse
|
92
|
Chen S, Zhang S. A structural optimization algorithm with stochastic forces and stresses. NATURE COMPUTATIONAL SCIENCE 2022; 2:736-744. [PMID: 38177372 DOI: 10.1038/s43588-022-00350-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/07/2022] [Accepted: 10/07/2022] [Indexed: 01/06/2024]
Abstract
We propose an algorithm for optimizations in which the gradients contain stochastic noise. This arises, for example, in structural optimizations when computations of forces and stresses rely on methods involving Monte Carlo sampling, such as quantum Monte Carlo or neural network states, or are performed on quantum devices that have intrinsic noise. Our proposed algorithm is based on the combination of two ingredients: an update rule derived from the steepest-descent method, and a staged scheduling of the targeted statistical error and step size, with position averaging. We compare it with commonly applied algorithms, including some of the latest machine learning optimization methods, and show that the algorithm consistently performs efficiently and robustly under realistic conditions. Applying this algorithm, we achieve full-degree optimizations in solids using ab initio many-body computations, by auxiliary-field quantum Monte Carlo with plane waves and pseudopotentials. A potential metastable structure in Si is discovered using density-functional calculations with synthetic noisy forces.
Collapse
Affiliation(s)
- Siyuan Chen
- Department of Physics, College of William & Mary, Williamsburg, VA, USA.
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, New York, NY, USA
| |
Collapse
|
93
|
Rojas-Cuervo A, Rey-González R. Electronic band gap on graphene induced by interaction with hydrogen cyanide. An DFT analysis. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111744] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
94
|
Lu Y, Zhao R, Zhang J, Liu M, Gao J. Minimal Active Space: NOSCF and NOSI in Multistate Density Functional Theory. J Chem Theory Comput 2022; 18:6407-6420. [DOI: 10.1021/acs.jctc.2c00908] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yangyi Lu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen518055, China
| | - Ruoqi Zhao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen518055, China
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province130023, China
| | - Jun Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen518055, China
| | - Meiyi Liu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen518055, China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen518055, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota55455, United States
| |
Collapse
|
95
|
Intramolecular Interactions in Derivatives of Uracil Tautomers. Molecules 2022; 27:molecules27217240. [PMID: 36364066 PMCID: PMC9656941 DOI: 10.3390/molecules27217240] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Revised: 10/17/2022] [Accepted: 10/21/2022] [Indexed: 11/17/2022] Open
Abstract
The influence of solvents on intramolecular interactions in 5- or 6-substituted nitro and amino derivatives of six tautomeric forms of uracil was investigated. For this purpose, the density functional theory (B97-D3/aug-cc-pVDZ) calculations were performed in ten environments (1 > ε > 109) using the polarizable continuum model (PCM) of solvation. The substituents were characterized by electronic (charge of the substituent active region, cSAR) and geometric parameters. Intramolecular interactions between non-covalently bonded atoms were investigated using the theory of atoms in molecules (AIM) and the non-covalent interaction index (NCI) method, which allowed discussion of possible interactions between the substituents and N/NH endocyclic as well as =O/−OH exocyclic groups. The nitro group was more electron-withdrawing in the 5 than in the 6 position, while the opposite effect was observed in the case of electron donation of the amino group. These properties of both groups were enhanced in polar solvents; the enhancement depended on the ortho interactions. Substitution or solvation did not change tautomeric preferences of uracil significantly. However, the formation of a strong NO∙∙∙HO intramolecular hydrogen bond in the 5-NO2 derivative stabilized the dienol tautomer from +17.9 (unsubstituted) to +5.4 kcal/mol (substituted, energy relative to the most stable diketo tautomer).
Collapse
|
96
|
Bursch M, Mewes J, Hansen A, Grimme S. Best-Practice DFT Protocols for Basic Molecular Computational Chemistry. Angew Chem Int Ed Engl 2022; 61:e202205735. [PMID: 36103607 PMCID: PMC9826355 DOI: 10.1002/anie.202205735] [Citation(s) in RCA: 139] [Impact Index Per Article: 69.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2022] [Indexed: 01/11/2023]
Abstract
Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.
Collapse
Affiliation(s)
- Markus Bursch
- Max-Planck-Institut für KohlenforschungKaiser-Wilhelm-Platz 145470Mülheim an der RuhrGermany
| | - Jan‐Michael Mewes
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| | - Andreas Hansen
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| | - Stefan Grimme
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| |
Collapse
|
97
|
Satalkar V, Benassi E, Mao Y, Pan X, Ran C, Chen X, Shao Y. Computational Investigation of Substituent Effects on the Fluorescence Wavelengths of Oxyluciferin Analogs. J Photochem Photobiol A Chem 2022; 431:114018. [PMID: 36407037 PMCID: PMC9673899 DOI: 10.1016/j.jphotochem.2022.114018] [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] [Indexed: 11/21/2022]
Abstract
Oxyluciferin, which is the light emitter for firefly bioluminescence, has been subjected to extensive chemical modifications to tune its emission wavelength and quantum yield. However, the exact mechanisms for various electron-donating and withdrawing groups to perturb the photophysical properties of oxyluciferin analogs are still not fully understood. To elucidate the substituent effects on the fluorescence wavelength of oxyluciferin analogs, we applied the absolutely localized molecular orbitals (ALMO)-based frontier orbital analysis to assess various types of interactions (i.e. permanent electrostatics/exchange repulsion, polarization, occupied-occupied orbital mixing, virtual-virtual orbital mixing, and charge-transfer) between the oxyluciferin and substituent orbitals. We suggested two distinct mechanisms that can lead to red-shifted oxyluciferin emission wavelength, a design objective that can help increase the tissue penetration of bioluminescence emission. Within the first mechanism, an electron-donating group (such as an amino or dimethylamino group) can contribute its highest occupied molecular orbital (HOMO) to an out-of-phase combination with oxyluciferin's HOMO, thus raising the HOMO energy of the substituted analog and narrowing its HOMO-LUMO gap. Alternatively, an electron-withdrawing group (such as a nitro or cyano group) can participate in an in-phase virtual-virtual orbital mixing of fragment LUMOs, thus lowering the LUMO energy of the substituted analog. Such an ALMO-based frontier orbital analysis is expected to lead to intuitive principles for designing analogs of not only the oxyluciferin molecule, but also many other functional dyes.
Collapse
Affiliation(s)
- Vardhan Satalkar
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, OK 73019, USA
| | - Enrico Benassi
- Novosibirsk State University, Novosibirsk, 630090, Russia
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, CA 94305, USA
| | - Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, OK 73019, USA
| | - Chongzhao Ran
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital/Harvard Medical School, Charlestown, Boston, MA 02129, USA
| | - Xiaoyuan Chen
- Yong Loo Lin School of Medicine and Faculty of Engineering, National University of Singapore, 117597, Singapore
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, OK 73019, USA
| |
Collapse
|
98
|
Advances in the development of novel green liquids: thymol/water, thymol/urea and thymol/phenylacetic acid as innovative hydrophobic natural deep eutectic solvents. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
99
|
Fedorov IA, Korabelnikov DV. AB INITIO STUDY OF THE COMPRESSIBILITY AND ELECTRONIC PROPERTIES OF CRYSTALLINE PURINE. J STRUCT CHEM+ 2022. [DOI: 10.1134/s0022476622100134] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
100
|
Li SZ, Wei YX, Huang Y, Dong WK. Counteranion-driven self-assembly of di- and tetra-nuclear Zn(II) single‐armed salamo‐type complexes. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.133473] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|