1
|
Ventura E, Rodrigues GP, Leitão EFV, do Monte SA. Theoretical Study of an Authentic Hydrocarbon Ion Pair. ACS OMEGA 2024; 9:34981-34989. [PMID: 39157103 PMCID: PMC11325399 DOI: 10.1021/acsomega.4c04914] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/24/2024] [Revised: 07/17/2024] [Accepted: 07/22/2024] [Indexed: 08/20/2024]
Abstract
For many years researchers believed that hydrocarbons only contain covalent bonds. However, since 1985 Okamoto et al. demonstrated the formation of hydrocarbon salts in several systems, demolishing the structural principle that hydrocarbons only contain covalent bonds. Despite the great importance of this outcome to the study of chemical bonds, quantum chemical calculations on these systems are essentially nonexistent. The stability of the hydrocarbon ions along with the steric hindrance associated with the formation of the covalent bond contribute to their occurrence either in solution (dissociated) or in the solid state. These facts along with the common formation of ion pairs in solvents of low polarity motivated us to search for hydrocarbon ion pairs in the gas phase. Its energetics has also been studied in four nonprotic solvents, through a continuum solvation model (CPCM). DFT and CASSCF calculations indicate a metastable and highly polar ion pair between the tricyclopropylcyclopropenylium cation and a simplified Kuhn's anion. The barrier to the covalent structure varies from ∼4.8 to 14.4 kcal/mol, while the energy difference between the ion pair and the covalent form varies from ∼4.3 to 25.4 kcal/mol. The obtained theoretical results along with previous experimental results suggest the following strategy to obtain kinetically and thermodynamically stable hydrocarbon ion pairs: choose very stable hydrocarbon ions and systematically increase the steric hindrance between them.
Collapse
Affiliation(s)
- Elizete Ventura
- Departamento
de Química, CCEN, Universidade Federal
da Paraíba, 58059-900 João Pessoa, Brazil
| | | | - Ezequiel Fragoso Vieira Leitão
- Unidade
Acadêmica de Ciências Exatas e da Natureza, Universidade Federal de Campina Grande, Cajazeiras, PB 58900-000, Brazil
| | - Silmar Andrade do Monte
- Departamento
de Química, CCEN, Universidade Federal
da Paraíba, 58059-900 João Pessoa, Brazil
| |
Collapse
|
2
|
Chen Y, Liu F, Wei S, Xia Y, Li X, Liu S, Zhang X, Tang S, Shen D, Dong W, Yang S. First-Principles Investigation of Phosphorus-Doped Graphitic Carbon Nitride as Anchoring Material for the Lithium-Sulfur Battery. Molecules 2024; 29:2746. [PMID: 38930812 PMCID: PMC11206158 DOI: 10.3390/molecules29122746] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2024] [Revised: 06/05/2024] [Accepted: 06/06/2024] [Indexed: 06/28/2024] Open
Abstract
The utilization of lithium-sulfur battery is hindered by various challenges, including the "shuttle effect", limited sulfur utilization, and the sluggish conversion kinetics of lithium polysulfides (LiPSs). In the present work, a theoretical design for the viability of graphitic carbon nitride (g-C3N4) and phosphorus-doping graphitic carbon nitride substrates (P-g-C3N4) as promising host materials in a Li-S battery was conducted utilizing first-principles calculations. The PDOS shows that when the P atom is introduced, the 2p of the N atom is affected by the 2p orbital of the P atom, which increases the energy band of phosphorus-doping substrates. The energy bands of PC and Pi are 0.12 eV and 0.20 eV, respectively. When the lithium polysulfides are adsorbed on four substrates, the overall adsorption energy of PC is 48-77% higher than that of graphitic carbon nitride, in which the charge transfer of long-chain lithium polysulfides increase by more than 1.5-fold. It is found that there are powerful Li-N bonds between lithium polysulfides and P-g-C3N4 substrates. Compared with the graphitic carbon nitride monolayer, the anchoring effect of the LiPSs@P-g-C3N4 substrate is enhanced, which is beneficial for inhibiting the shuttle of high-order lithium polysulfides. Furthermore, the catalytic performance of the P-g-C3N4 substrate is assessed in terms of the S8 reduction pathway and the decomposition of Li2S; the decomposition energy barrier of the P-g-C3N4 substrate decrease by 10% to 18%. The calculated results show that P-g-C3N4 can promote the reduction of S8 molecules and Li-S bond cleavage within Li2S, thus improving the utilization of sulfur-active substances and the ability of rapid reaction kinetics. Therefore, the P-g-C3N4 substrates are a promising high-performance lithium-sulfur battery anchoring material.
Collapse
Affiliation(s)
- Yuehui Chen
- School of Mining, Liaoning Technical University, Fuxin 123000, China; (Y.C.); (F.L.); (S.W.); (Y.X.)
- College of Science, Liaoning Technical University, Fuxin 123000, China;
| | - Fengxia Liu
- School of Mining, Liaoning Technical University, Fuxin 123000, China; (Y.C.); (F.L.); (S.W.); (Y.X.)
| | - Shuang Wei
- School of Mining, Liaoning Technical University, Fuxin 123000, China; (Y.C.); (F.L.); (S.W.); (Y.X.)
| | - Yingkai Xia
- School of Mining, Liaoning Technical University, Fuxin 123000, China; (Y.C.); (F.L.); (S.W.); (Y.X.)
| | - Xiaodong Li
- College of Material Science and Engineering, Liaoning Technical University, Fuxin 123000, China; (X.L.); (X.Z.); (S.T.); (D.S.); (W.D.)
| | - Shengnan Liu
- College of Science, Liaoning Technical University, Fuxin 123000, China;
| | - Xu Zhang
- College of Material Science and Engineering, Liaoning Technical University, Fuxin 123000, China; (X.L.); (X.Z.); (S.T.); (D.S.); (W.D.)
| | - Shuwei Tang
- College of Material Science and Engineering, Liaoning Technical University, Fuxin 123000, China; (X.L.); (X.Z.); (S.T.); (D.S.); (W.D.)
| | - Ding Shen
- College of Material Science and Engineering, Liaoning Technical University, Fuxin 123000, China; (X.L.); (X.Z.); (S.T.); (D.S.); (W.D.)
| | - Wei Dong
- College of Material Science and Engineering, Liaoning Technical University, Fuxin 123000, China; (X.L.); (X.Z.); (S.T.); (D.S.); (W.D.)
| | - Shaobin Yang
- School of Mining, Liaoning Technical University, Fuxin 123000, China; (Y.C.); (F.L.); (S.W.); (Y.X.)
- College of Material Science and Engineering, Liaoning Technical University, Fuxin 123000, China; (X.L.); (X.Z.); (S.T.); (D.S.); (W.D.)
| |
Collapse
|
3
|
Britto NJ, Sen A, Rajaraman G. Unravelling the Effect of Acid-Driven Electron Transfer in High-Valent Fe IV =O/Mn IV =O Species and Its Implications for Reactivity. Chem Asian J 2023; 18:e202300773. [PMID: 37855305 DOI: 10.1002/asia.202300773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2023] [Revised: 10/12/2023] [Accepted: 10/16/2023] [Indexed: 10/20/2023]
Abstract
The electron transfer (ET) step is one of the crucial processes in biochemical redox reactions that occur in nature and has been established as a key step in dictating the reactivity of high-valent metal-oxo species. Although metalloenzymes possessing metal-oxo units at their active site are typically associated with outer-sphere electron transfer (OSET) processes, biomimetic models, in contrast, have been found to manifest either an inner-sphere electron transfer (ISET) or OSET mechanism. This distinction is clearly illustrated through the behaviour of [(N4Py)MnIV (O)]2+ (1) and [(N4Py)FeIV (O)]2+ (2) complexes, where complex 1 showcases an OSET mechanism, while complex 2 exhibits an ISET mechanism, especially evident in their reactions involving C-H bond activation and oxygen atom transfer reactions in the presence of a Lewis/Bronsted acid. However, the precise reason for this puzzling difference remains elusive. This work unveils the origin of the perplexing inner-sphere vs outer-sphere electron transfer process (ISET vs OSET) in [(N4Py)MnIV (O)]2+ (1) and [(N4Py)FeIV (O)]2+ (2) species in the presence of Bronsted acid. The calculations indicate that when the substrate (toluene) approaches both 1 and 2 that is hydrogen bonded with two HOTf molecules (denoted as 1-HOTf and 2-HOTf, respectively), proton transfer from one of the HOTf molecules to the metal-oxo unit is triggered and a simultaneous electron transfer occurs from toluene to the metal centre. Interestingly, the preference for OSET by 1-HOTf is found to originate from the choice of MnIV =O centre to abstract spin-down (β) electron from toluene to its δ(dxy ) orbital. On the other hand, in 2-HOTf, a spin state inversion from triplet to quintet state takes place during the proton (from HOTf) coupled electron transfer (from toluene) preferring a spin-up (α) electron abstraction to its σ* (dz 2 ) orbital mediated by HOTf giving rise to ISET. In addition, 2-HOTf was calculated to possess a larger reorganisation energy, which facilitates the ISET process via the acid. The absence of spin-inversion and smaller reorganisation energy switch the mechanism to OSET for 1-HOTf. Therefore, for the first time, the significance of spin-state and spin-inversion in the electron transfer process has been identified and demonstrated within the realm of high-valent metal-oxo chemistry. This discovery holds implications for the potential involvement of high-valent Mn-oxo species in performing similar transformative processes within Photosystem II.
Collapse
Affiliation(s)
| | - Asmita Sen
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, 400 076, India
| | - Gopalan Rajaraman
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, 400 076, India
| |
Collapse
|
4
|
Lutz JJ, Jensen DS, Hubbard JA. Deposition products predicted from conceptual DFT: The hydrolysis reactions of MoF6, WF6, and UF6. J Chem Phys 2023; 159:184305. [PMID: 37962449 DOI: 10.1063/5.0176552] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Accepted: 10/23/2023] [Indexed: 11/15/2023] Open
Abstract
Metal hexafluorides hydrolyze at ambient temperature to deposit compounds having fluorine-to-oxygen ratios that depend upon the identity of the metal. Uranium-hexafluoride hydrolysis, for example, deposits uranyl fluoride (UO2F2), whereas molybdenum hexafluoride (MoF6) and tungsten hexafluoride deposit trioxides. Here, we pursue general strategies enabling the prediction of depositing compounds resulting from multi-step gas-phase reactions. To compare among the three metal-hexafluoride hydrolyses, we first investigate the mechanism of MoF6 hydrolysis using hybrid density functional theory (DFT). Intermediates are then validated by performing anharmonic vibrational simulations and comparing with infrared spectra [McNamara et al., Phys. Chem. Chem. Phys. 25, 2990 (2023)]. Conceptual DFT, which is leveraged here to quantitatively evaluate site-specific electrophilicity and nucleophilicity metrics, is found to reliably predict qualitative deposition propensities for each intermediate. In addition to the nucleophilic potential of the oxygen ligands, several other contributing characteristics are discussed, including amphoterism, polyvalency, fluxionality, steric hindrance, dipolar strength, and solubility. To investigate the structure and composition of pre-nucleation clusters, an automated workflow is presented for the simulation of particle growth. The workflow entails a conformer search at the density functional tight-binding level, structural refinement at the hybrid DFT level, and computation of a composite free-energy profile. Such profiles can be used to estimate particle nucleation kinetics. Droplet formation is also considered, which helps to rationalize the different UO2F2 particle morphologies observed under varying levels of humidity. Development of predictive methods for simulating physical and chemical deposition processes is important for the advancement of material manufacturing involving coatings and thin films.
Collapse
Affiliation(s)
- Jesse J Lutz
- Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Daniel S Jensen
- Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Joshua A Hubbard
- Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| |
Collapse
|
5
|
Santos CV, Monteiro SA, Soares ASC, Souto ICA, Moura RT. Decoding Chemical Bonds: Assessment of the Basis Set Effect on Overlap Electron Density Descriptors and Topological Properties in Comparison to QTAIM. J Phys Chem A 2023; 127:7997-8014. [PMID: 37703453 DOI: 10.1021/acs.jpca.3c04504] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2023]
Abstract
Quantum chemical bonding descriptors based on the total and overlap density can provide valuable information about chemical interactions in different systems. However, these descriptors can be sensitive to the basis set used. To address this, different numerical treatments of electron density have been proposed to reduce the basis set dependency. In this work, we introduce overlap properties (OPs) obtained through numerical treatment of the electron density and present the topology of overlap density (TOP) for the first time. We compare the basis set dependency of numerical OP and TOP descriptors with their quantum theory of atoms in molecules (QTAIM) counterparts, considering the total electron density. Three single (C-C, C-O, and C-F) bonds in ethane, methanol, and fluoromethane and two double (C═C and C═O) bonds in ethene and formaldehyde were analyzed. Diatomic molecules Li-X with X = F, Cl, and Br were also analyzed. Eight parameters, including QTAIM descriptors and OP/TOP descriptors, are used to assess the basis dependency at the ωB97X-D level of theory using 28 basis sets from three classes: Pople, Ahlrichs, and Dunning. The study revealed that the topological overlap electron density properties exhibit comparatively lesser dependence on the basis set compared to their total electron density counterparts. Remarkably, these properties retain their chemical significance even with reduced basis set dependency. Similarly, numerical OP descriptors show less basis set dependency than their QTAIM counterparts. The excess of polarization functions increases charge concentration in the interatomic region and influences both QTAIM and OP descriptors. The basis sets Def2TZVP, 6-31++G(d,p), 6-311++G(d,p), cc-pVDZ, cc-pVTZ, and cc-pVQZ demonstrate reduced variability for the tested bond classes in this study, with particular emphasis on the triple-ζ quality Ahlrichs' basis set. We recommend against using basis sets with numerous polarization functions, such as augmented Dunning's and Ahlrichs' quadruple-ζ.
Collapse
Affiliation(s)
- Carlos V Santos
- Department of Chemistry, Federal University of Paraiba, Joao Pessoa, Paraiba 58051-970, Brazil
| | - Shirlene A Monteiro
- Department of Chemistry, State University of Paraiba, Campina Grande, Paraiba 58051-970, Brazil
| | - Amanda S C Soares
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraiba, Areia, Paraiba 58397-000, Brazil
| | - Isabeli C A Souto
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraiba, Areia, Paraiba 58397-000, Brazil
| | - Renaldo T Moura
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraiba, Areia, Paraiba 58397-000, Brazil
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| |
Collapse
|
6
|
Dakkouri M. A Theoretical Investigation of Novel Sila- and Germa-Spirocyclic Imines and Their Relevance for Electron-Transporting Materials and Drug Discovery. Molecules 2023; 28:6298. [PMID: 37687127 PMCID: PMC10489060 DOI: 10.3390/molecules28176298] [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: 07/23/2023] [Revised: 08/21/2023] [Accepted: 08/22/2023] [Indexed: 09/10/2023] Open
Abstract
A new class of spirocyclic imines (SCIs) has been theoretically investigated by applying a variety of quantum chemical methods and basis sets. The uniqueness of these compounds is depicted by various peculiarities, e.g., the incidence of planar six-membered rings each with two imine groups (two π bonds) and the incorporation of the isosteres carbon, silicon, or germanium spiro centers. Additional peculiarities of these novel SCIs are mirrored by their three-dimensionality, the simultaneous occurrence of nucleophilic and electrophilic centers, and the cross-hyperconjugative (spiro-conjugation) interactions, which provoke charge mobility along the spirocyclic scaffold. Substitution of SCIs with strong electron-withdrawing substituents, like the cyano group or fluorine, enhances their docking capability and impacts their reactivity and charge mobility. To gain thorough knowledge about the molecular properties of these SCIs, their structures have been optimized and various quantum chemical concepts and models were applied, e.g., full NBO analysis and the frontier molecular orbitals (FMOs) theory (HOMO-LUMO energy gap) and the chemical reactivity descriptors derived from them. For the assessment of the charge density distribution along the SCI framework, additional complementary quantum chemical methods were used, e.g., molecular electrostatic potential (MESP) and Bader's QTAIM. Additionally, using the aromaticity index NICS (nuclear independent chemical shift) and other criteria, it could be shown that the investigated cross-hyperconjugated sila and germa SCIs are spiro-aromatics of the Heilbronner Craig-type Möbius aromaticity.
Collapse
Affiliation(s)
- Marwan Dakkouri
- Department of Electrochemistry, University of Ulm, D-89069 Ulm, Germany
| |
Collapse
|
7
|
Kirchhoff B, Jung C, Gaissmaier D, Braunwarth L, Fantauzzi D, Jacob T. In silico characterization of nanoparticles. Phys Chem Chem Phys 2023; 25:13228-13243. [PMID: 37161752 DOI: 10.1039/d3cp01073b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Nanoparticles (NPs) make for intriguing heterogeneous catalysts due to their large active surface area and excellent and often size-dependent catalytic properties that emerge from a multitude of chemically different surface reaction sites. NP catalysts are, in principle, also highly tunable: even small changes to the NP size or surface facet composition, doping with heteroatoms, or changes of the supporting material can significantly alter their physicochemical properties. Because synthesis of size- and shape-controlled NP catalysts is challenging, the ability to computationally predict the most favorable NP structures for a catalytic reaction of interest is an in-demand skill that can help accelerate and streamline the material optimization process. Fundamentally, simulations of NP model systems present unique challenges to computational scientists. Not only must considerable methodological hurdles be overcome in performing calculations with hundreds to thousands of atoms while retaining appropriate accuracy to be able to probe the desired properties. Also, the data generated by simulations of NPs are typically more complex than data from simulations of, for example, single crystal surface models, and therefore often require different data analysis strategies. To this end, the present work aims to review analytical methods and data analysis strategies that have proven useful in extracting thermodynamic trends from NP simulations.
Collapse
Affiliation(s)
- Björn Kirchhoff
- Institute of Electrochemistry, Ulm University, Albert-Einstein-Allee 47, 89081 Ulm, Germany.
| | - Christoph Jung
- Institute of Electrochemistry, Ulm University, Albert-Einstein-Allee 47, 89081 Ulm, Germany.
- Helmholtz-Institute Ulm (HIU) Electrochemical Energy Storage, Helmholtz-Straße 16, 89081 Ulm, Germany
- Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Daniel Gaissmaier
- Institute of Electrochemistry, Ulm University, Albert-Einstein-Allee 47, 89081 Ulm, Germany.
- Helmholtz-Institute Ulm (HIU) Electrochemical Energy Storage, Helmholtz-Straße 16, 89081 Ulm, Germany
- Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Laura Braunwarth
- Institute of Electrochemistry, Ulm University, Albert-Einstein-Allee 47, 89081 Ulm, Germany.
| | - Donato Fantauzzi
- Institute of Electrochemistry, Ulm University, Albert-Einstein-Allee 47, 89081 Ulm, Germany.
| | - Timo Jacob
- Institute of Electrochemistry, Ulm University, Albert-Einstein-Allee 47, 89081 Ulm, Germany.
- Helmholtz-Institute Ulm (HIU) Electrochemical Energy Storage, Helmholtz-Straße 16, 89081 Ulm, Germany
- Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany
| |
Collapse
|
8
|
Rivic F, Lehr A, Schäfer R. Scaling of the permanent electric dipole moment in isolated silicon clusters with near-spherical shape. Phys Chem Chem Phys 2023; 25:13376-13382. [PMID: 37145168 DOI: 10.1039/d3cp00926b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
In silicon clusters a structural transition from prolate to near-spherical structures takes place at a size of about 25-30 atoms. While some of the prolate clusters are very polar, there has been no experimental evidence of the presence of dipole moments in larger silicon clusters with near-spherical shape. By means of electric molecular beam deflection experiments at cryogenic temperatures, it was possible to prove for the first time that SiN clusters with more than N = 30 atoms are also polar. Interestingly, the dipole moment per atom for clusters in the range between 30 and 80 or 90 atoms is almost constant and amounts to 0.02 D. This unusual behaviour manifests in effective polarizabilities increasing linearly with cluster size. The dipolar contribution to the polarizability means that SiN clusters with N = 80 atoms can be polarized more than twice as well as a correspondingly small sphere with the dielectric properties of bulk α-Si. This finding is analysed with quantum chemical calculations concerning the geometric structure as well as the charge distribution and is related to the dielectric behaviour of polar semiconductor nanocrystals.
Collapse
Affiliation(s)
- Filip Rivic
- Technical University of Darmstadt, Eduard-Zintl-Institute, Alarich-Weiss-Straße 8, 64287 Darmstadt, Germany.
| | - Andreas Lehr
- Technical University of Darmstadt, Eduard-Zintl-Institute, Alarich-Weiss-Straße 8, 64287 Darmstadt, Germany.
| | - Rolf Schäfer
- Technical University of Darmstadt, Eduard-Zintl-Institute, Alarich-Weiss-Straße 8, 64287 Darmstadt, Germany.
| |
Collapse
|
9
|
Wang Y, Liu M, Zhou J, Chen S, Li J, Liu J, Sun Y, Shi Z. DFT Study on the Oxidation Mechanism of Common Cyclic Carbonates in the Presence of BF 4- Anions. J Phys Chem A 2023; 127:3958-3965. [PMID: 37115673 DOI: 10.1021/acs.jpca.2c07827] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
Oxidative decomposition reactions of common cyclic carbonates in the presence of BF4- anions were investigated using density functional theory. A polarized continuum model was utilized to model solvent effects in the oxidation of ethylene carbonate (EC) and propylene carbonate (PC) clusters. We have found that the presence of BF4- significantly reduces EC and PC oxidation stability, from 7.11 to 6.17 and from 7.10 to 6.06 V (vs Li+/Li), respectively. The sequence of EC and PC oxidative decomposition paths and the oxidative products were affected by the BF4- anion. The decomposition products of the oxidized EC-BF4- contained CO2, vinyl alcohol, and acetaldehyde, while the decomposition products of the oxidized PC-BF4- contained CO2, acetone, and propanal, in agreement with the previous experimental studies. The oxidative decomposition reactions for PC-BF4- are compared with those for the isolated PC, PC2, PC-ClO4-, and PC-PF6-.
Collapse
Affiliation(s)
- Yating Wang
- Institute of Batteries, School of Materials and Energy, Guangdong University of Technology, Guangzhou 510006, China
| | - Mingzhu Liu
- Engineering Research Center of MTEES (Ministry of Education), Research Center of BMET (Guangdong Province), Engineering Laboratory of OFMHEB (Guangdong Province), Key Laboratory of ETESPG (GHEI), and Innovative Platform for ITBMD (Guangzhou Municipality), School of Chemistry, South China Normal University, Guangzhou 510006, China
| | - Jiasheng Zhou
- Institute of Batteries, School of Materials and Energy, Guangdong University of Technology, Guangzhou 510006, China
| | - Shaoru Chen
- School of Chemical Engineering and Light Industry, Guangdong University of Technology, Guangzhou 510006, Guangdong, China
| | - Jianhui Li
- Engineering Research Center of MTEES (Ministry of Education), Research Center of BMET (Guangdong Province), Engineering Laboratory of OFMHEB (Guangdong Province), Key Laboratory of ETESPG (GHEI), and Innovative Platform for ITBMD (Guangzhou Municipality), School of Chemistry, South China Normal University, Guangzhou 510006, China
| | - Jun Liu
- Institute of Batteries, School of Materials and Energy, Guangdong University of Technology, Guangzhou 510006, China
| | - Yang Sun
- School of Materials, Sun Yat-sen University, Shenzhen 518107, China
| | - Zhicong Shi
- Institute of Batteries, School of Materials and Energy, Guangdong University of Technology, Guangzhou 510006, China
| |
Collapse
|
10
|
Druchok M, Krasnov V, Krokhmalskii T, Cardoso E Bufalo T, Martins de Souza S, Rojas O, Derzhko O. Toward a quasiphase transition in the single-file chain of water molecules: Simple lattice model. J Chem Phys 2023; 158:104304. [PMID: 36922143 DOI: 10.1063/5.0133720] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023] Open
Abstract
Recently, Ma et al. [Phys. Rev. Lett. 118, 027402 (2017)] have suggested that water molecules encapsulated in (6,5) single-wall carbon nanotube experience a temperature-induced quasiphase transition around 150 K interpreted as changes in the water dipoles orientation. We discuss further this temperature-driven quasiphase transition performing quantum chemical calculations and molecular dynamics simulations and, most importantly, suggesting a simple lattice model to reproduce the properties of the one-dimensional confined finite arrays of water molecules. The lattice model takes into account not only the short-range and long-range interactions but also the rotations in a narrow tube, and both ingredients provide an explanation for a temperature-driven orientational ordering of the water molecules, which persists within a relatively wide temperature range.
Collapse
Affiliation(s)
- Maksym Druchok
- Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine, Svientsitskii Street 1, 79011 L'viv, Ukraine
| | - Volodymyr Krasnov
- Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine, Svientsitskii Street 1, 79011 L'viv, Ukraine
| | - Taras Krokhmalskii
- Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine, Svientsitskii Street 1, 79011 L'viv, Ukraine
| | | | - Sergio Martins de Souza
- Departamento de Fisica, Universidade Federal de Lavras, CP 3037, 37200-000 Lavras-MG, Brazil
| | - Onofre Rojas
- Departamento de Fisica, Universidade Federal de Lavras, CP 3037, 37200-000 Lavras-MG, Brazil
| | - Oleg Derzhko
- Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine, Svientsitskii Street 1, 79011 L'viv, Ukraine
| |
Collapse
|
11
|
Bhadoria P, Ramanathan V. Computational underpinnings for the dimerization of para-aminothiophenol to dimercaptoazobenzene on copper surface. Chem Phys 2023. [DOI: 10.1016/j.chemphys.2023.111910] [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]
|
12
|
Nanoscale Iron-Based Metal-Organic Frameworks: Incorporation of Functionalized Drugs and Degradation in Biological Media. Int J Mol Sci 2023; 24:ijms24043362. [PMID: 36834775 PMCID: PMC9965190 DOI: 10.3390/ijms24043362] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2023] [Revised: 02/03/2023] [Accepted: 02/06/2023] [Indexed: 02/11/2023] Open
Abstract
Metal-organic frameworks (MOFs) attract growing interest in biomedical applications. Among thousands of MOF structures, the mesoporous iron(III) carboxylate MIL-100(Fe) (MIL stands for the Materials of Lavoisier Institute) is among the most studied MOF nanocarrier, owing to its high porosity, biodegradability, and lack of toxicity. Nanosized MIL-100(Fe) particles (nanoMOFs) readily coordinate with drugs leading to unprecedented payloads and controlled release. Here, we show how the functional groups of the challenging anticancer drug prednisolone influence their interactions with the nanoMOFs and their release in various media. Molecular modeling enabled predicting the strength of interactions between prednisolone-bearing or not phosphate or sulfate moieties (PP and PS, respectively) and the oxo-trimer of MIL-100(Fe) as well as understanding the pore filling of MIL-100(Fe). Noticeably, PP showed the strongest interactions (drug loading up to 30 wt %, encapsulation efficiency > 98%) and slowed down the nanoMOFs' degradation in simulated body fluid. This drug was shown to bind to the iron Lewis acid sites and was not displaced by other ions in the suspension media. On the contrary, PS was entrapped with lower efficiencies and was easily displaced by phosphates in the release media. Noticeably, the nanoMOFs maintained their size and faceted structures after drug loading and even after degradation in blood or serum after losing almost the totality of the constitutive trimesate ligands. Scanning electron microscopy with high annular dark field (STEM-HAADF) in conjunction with X-Ray energy-dispersive spectrometry (XEDS) was a powerful tool enabling the unraveling of the main elements to gain insights on the MOF structural evolution after drug loading and/or upon degradation.
Collapse
|
13
|
North SC, Jorgensen KR, Pricetolstoy J, Wilson AK. Population analysis and the effects of Gaussian basis set quality and quantum mechanical approach: main group through heavy element species. Front Chem 2023; 11:1152500. [PMID: 37153525 PMCID: PMC10154537 DOI: 10.3389/fchem.2023.1152500] [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: 01/27/2023] [Accepted: 04/03/2023] [Indexed: 05/09/2023] Open
Abstract
Atomic charge and its distribution across molecules provide important insight into chemical behavior. Though there are many studies on various routes for the determination of atomic charge, there are few studies that examine the broader impact of basis set and quantum method used over many types of population analysis methods across the periodic table. Largely, such a study of population analysis has focused on main-group species. In this work, atomic charges were calculated using several population analysis methods including orbital-based methods (Mulliken, Löwdin, and Natural Population Analysis), volume-based methods (Atoms-in-Molecules (AIM) and Hirshfeld), and potential derived charges (CHELP, CHELPG, and Merz-Kollman). The impact of basis set and quantum mechanical method choices upon population analysis has been considered. The basis sets utilized include Pople (6-21G**, 6-31G**, 6-311G**) and Dunning (cc-pVnZ, aug-cc-pVnZ; n = D, T, Q, 5) basis sets for main group molecules. For the transition metal and heavy element species examined, relativistic forms of the correlation consistent basis sets were used. This is the first time the cc-pVnZ-DK3 and cc-pwCVnZ-DK3 basis sets have been examined with respect to their behavior across all levels of basis sets for atomic charges for an actinide. The quantum methods chosen include two density functional (PBE0 and B3LYP), Hartree-Fock, and second-order Møller-Plesset perturbation theory (MP2) approaches.
Collapse
Affiliation(s)
- Sasha C. North
- Department of Chemistry, Michigan State University, East Lansing, MI, United States
| | - Kameron R. Jorgensen
- Department of Biology and Chemistry, Texas A&M International University, Laredo, TX, United States
| | - Jason Pricetolstoy
- Department of Chemistry, Michigan State University, East Lansing, MI, United States
| | - Angela K. Wilson
- Department of Chemistry, Michigan State University, East Lansing, MI, United States
- *Correspondence: Angela K. Wilson,
| |
Collapse
|
14
|
Maciołek U, Mendyk E, Kamiński DM, Dranka M, Mazur L, Kuźniar A, Kalembkiewicz J, Kozioł AE. Structural and spectroscopic studies of Na+ – Quercetin-5′-sulfonic acid polymeric complexes obtained via solvothermal synthesis. Polyhedron 2022. [DOI: 10.1016/j.poly.2022.116083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/31/2022]
|
15
|
Hacıefendioǧlu T, Yildirim E. Design Principles for the Acceptor Units in Donor-Acceptor Conjugated Polymers. ACS OMEGA 2022; 7:38969-38978. [PMID: 36340112 PMCID: PMC9631416 DOI: 10.1021/acsomega.2c04713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Accepted: 10/06/2022] [Indexed: 06/16/2023]
Abstract
More than 50 different acceptor units from the experimental literature have been modeled, analyzed, and compared by using the computationally extracted data from the density functional theory (DFT) perspective for tetramer structures in the form of (D-B-A-B)4 (D, donor; A, acceptor; B, bridge) with fixed donor and bridge units. Comparison of dihedral angle between acceptor, donor, and bridge units, bond order, and hyperpolarizability reveals that these three structural properties have a dominant effect on the frontier electronic energy levels of the acceptor units. Systematic investigation of the structural properties has demonstrated the band gap energy dependency of the acceptor units on the planarity, conjugation, and the electron delocalization. Substitution effect, morphological alternation, and insertion of π-electron deficient atoms in A unit have also an important role to determine physical properties of the donor-acceptor conjugated polymers. This benchmark study will be beneficial for the band gap engineering and molecular design of the donor-acceptor copolymers using different acceptor units for the organic electronic applications.
Collapse
Affiliation(s)
| | - Erol Yildirim
- Department
of Chemistry, Middle East Technical University, 06800 Ankara, Turkey
- Department
of Polymer Science and Technology, Middle
East Technical University, 06800 Ankara, Turkey
- Department
of Micro- and Nanotechnology, Middle East
Technical University, 06800 Ankara, Turkey
| |
Collapse
|
16
|
An Overview of Basis Set Effects for Diatomic Boron Nitride Compounds (B2N(∓,0)): A Quantum Symmetry Breaking. QUANTUM REPORTS 2022. [DOI: 10.3390/quantum4030024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
The symmetry breaking (SB) of B2 not only exhibits an energy barrier for ionic or neutral forms dependent on various basis sets but it also exhibits a few SBs due to the asymmetry stretching and bending mode interactions. SB obeys the mechanical quantum theorem among discrete symmetries and their connection to the spin statistics in physical sciences. In this investigation, the unusual amount of energy barrier of SBs appeared upon the orbit–orbit coupling of BNB (both radical and ions) between transition states and the ground state. Our goal in this study is to understand the difference among the electromagnetic structures of the (B2N(∓,0)) variants due to effects of various basis sets and methods and also the quantum symmetry breaking phenomenon. In the D∞h point group of (B2N(∓,0)) variants, the unpaired electron is delocalized, while in the asymmetric C∞v point group, it is localized on either one of the B atoms. Structures with broken symmetry, C∞v, can be stable by interacting with the D∞h point group. In viewpoints of quantum chemistry, the second-order Jahn–Teller effect permits the unpaired electron to localize on boron atom, rather than being delocalized. In this study, we observed that the energy barrier of SB for BNB increases by post HF methods.
Collapse
|
17
|
Irregular structure of the hydrated Ag+ in aqueous solution and its Dynamics: An insight from perturbation theory hybrid forces molecular dynamics simulation. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119688] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
18
|
Liu S, Luan B. Benchmarking various types of partial atomic charges for classical all-atom simulations of metal-organic frameworks. NANOSCALE 2022; 14:9466-9473. [PMID: 35748335 DOI: 10.1039/d2nr00354f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The density derived electrostatic and chemical (DDEC) approach for calculating the charges of atoms in a metal-organic framework (MOF) is considered to be the most accurate (yet computationally costly) one among many charge-assignment methods. Here, we conducted a comparative study on five different types of atomic partial charges (namely CM5, Mulliken, Qeq, EQeq and PACMOF) prepared for a subset of MOFs with affordable computational costs and benchmarked them with respect to the DDEC charges, which is particularly relevant because currently most databases lack MOFs with pre-calculated DDEC charges. To find a suitable charge type alternative to the DDEC approach, we statistically ranked the five charge types based on two metrics, the relative standard deviation of charges and relative dipole moment difference, based on which we provide general guidance as well as suggestions for specific MOFs according to bond polarity analyses. Finally, we recommend a possible and more accurate parametrization scheme for future studies.
Collapse
Affiliation(s)
- Sizhe Liu
- IBM Thomas J. Watson Research Center, Yorktown Heights, New York 10598, USA.
| | - Binquan Luan
- IBM Thomas J. Watson Research Center, Yorktown Heights, New York 10598, USA.
| |
Collapse
|
19
|
Golchoubian D, Abbasi A. pH-Responsive Switching Properties of a Water-Soluble dinuclear Cu(II) Complex containing oximato bridges chelating ligand: Synthesis, crystal structure and computational study. Inorganica Chim Acta 2022. [DOI: 10.1016/j.ica.2022.121044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
20
|
Spiegel M. Current Trends in Computational Quantum Chemistry Studies on Antioxidant Radical Scavenging Activity. J Chem Inf Model 2022; 62:2639-2658. [PMID: 35436117 PMCID: PMC9198981 DOI: 10.1021/acs.jcim.2c00104] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
![]()
The antioxidative
nature of chemicals is now routinely studied
using computational quantum chemistry. Scientists are constantly proposing
new approaches to investigate those methods, and the subject is evolving
at a rapid pace. The goal of this review is to collect, consolidate,
and present current trends in a clear, methodical, and reference-rich
manner. This paper is divided into several sections, each of which
corresponds to a different stage of elaborations: preliminary concerns,
electronic structure analysis, and general reactivity (thermochemistry
and kinetics). The sections are further subdivided based on methodologies
used. Concluding remarks and future perspectives are presented based
on the remaining elements.
Collapse
Affiliation(s)
- Maciej Spiegel
- Department of Pharmacognosy and Herbal Medicines, Wroclaw Medical University, Borowska 211A, 50-556 Wroclaw, Poland
| |
Collapse
|
21
|
Wojtkowiak K, Jezierska A, Panek JJ. Revealing Intra- and Intermolecular Interactions Determining Physico-Chemical Features of Selected Quinolone Carboxylic Acid Derivatives. Molecules 2022; 27:2299. [PMID: 35408698 PMCID: PMC9000753 DOI: 10.3390/molecules27072299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2022] [Revised: 03/24/2022] [Accepted: 03/30/2022] [Indexed: 02/04/2023] Open
Abstract
The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car-Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion-the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.
Collapse
Affiliation(s)
| | - Aneta Jezierska
- Faculty of Chemistry, University of Wrocław, ul. F. Joliot-Curie 14, 50-383 Wrocław, Poland;
| | - Jarosław J. Panek
- Faculty of Chemistry, University of Wrocław, ul. F. Joliot-Curie 14, 50-383 Wrocław, Poland;
| |
Collapse
|
22
|
Badavath VN, Kumar A, Samanta PK, Maji S, Das A, Blum G, Jha A, Sen A. Determination of potential inhibitors based on isatin derivatives against SARS-CoV-2 main protease (m pro): a molecular docking, molecular dynamics and structure-activity relationship studies. J Biomol Struct Dyn 2022; 40:3110-3128. [PMID: 33200681 PMCID: PMC7682386 DOI: 10.1080/07391102.2020.1845800] [Citation(s) in RCA: 30] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Accepted: 10/29/2020] [Indexed: 12/27/2022]
Abstract
SARS-COV-2, the novel coronavirus and root of global pandemic COVID-19 caused a severe health threat throughout the world. Lack of specific treatments raised an effort to find potential inhibitors for the viral proteins. The recently invented crystal structure of SARS-CoV-2 main protease (Mpro) and its key role in viral replication; non-resemblance to any human protease makes it a perfect target for inhibitor research. This article reports a computer-aided drug design (CADD) approach for the screening of 118 compounds with 16 distinct heterocyclic moieties in comparison with 5 natural products and 7 repurposed drugs. Molecular docking analysis against Mpro protein were performed finding isatin linked with a oxidiazoles (A2 and A4) derivatives to have the best docking scores of -11.22 kcal/mol and -11.15 kcal/mol respectively. Structure-activity relationship studies showed a good comparison with a known active Mpro inhibitor and repurposed drug ebselen with an IC50 value of -0.67 μM. Molecular Dynamics (MD) simulations for 50 ns were performed for A2 and A4 supporting the stability of the two compounds within the binding pocket, largely at the S1, S2 and S4 domains with high binding energy suggesting their suitability as potential inhibitors of Mpro for SARS-CoV-2.
Collapse
Affiliation(s)
| | - Akhil Kumar
- Biotechnology Division, CSIR-Central Institute of Medicinal and Aromatic Plants, Lucknow, India
| | - Pralok K. Samanta
- School of Chemical and Bioprocess Engineering, University College Dublin, Dublin, Ireland
| | - Siddhartha Maji
- Department of Pharmaceutical Sciences and Technology, Birla Institute of Technology, Mesra, Ranchi, India
| | - Anik Das
- Department of Chemistry, Institute of Science, GITAM (Deemed to be University), Visakhapatnam, India
| | - Galia Blum
- Institute for Drug Research, The Hebrew University, Jerusalem, Israel
| | - Anjali Jha
- Department of Chemistry, Institute of Science, GITAM (Deemed to be University), Visakhapatnam, India
| | - Anik Sen
- Department of Chemistry, Institute of Science, GITAM (Deemed to be University), Visakhapatnam, India
| |
Collapse
|
23
|
Dwinandha D, Zhang B, Fujii M. Prediction of reaction mechanism for OH radical-mediated phenol oxidation using quantum chemical calculation. CHEMOSPHERE 2022; 291:132763. [PMID: 34740699 DOI: 10.1016/j.chemosphere.2021.132763] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Revised: 10/29/2021] [Accepted: 10/31/2021] [Indexed: 06/13/2023]
Abstract
Understanding the reaction mechanism of OH•-mediated oxidation of organic micropollutants (OMPs) contributes to the assessment and development of advanced oxidation processes (AOPs) for removal of OMPs in water environment. In this study, a theoretical approach using quantum chemical calculation (QCC) was employed to investigate the prediction accuracy of the reaction mechanism (i.e., reaction site and rate) for OH•-mediated oxidation of phenol, where the hydroquinone and catechol are generated as transformation products (TPs) via radical and electrophilic reactions. We compared three different levels of theory (Hartree-Fock, B3LYP, and M06-2X) with 6-311 + G (2d,2p)/SMD, and the reaction site and rate constants were predicted by the Fukui function and transition state theory, respectively. Overall, the prediction accuracy of the TPs formation mechanism was the highest in the calculations using M06-2X. For example, the initial OH• addition to phenol was predicted to occur with a probability of 77% for the ortho position and 23% for the para position, which was consistent with the experimental observation. By applying the transition state theory, the rate constants toward TPs formation pathway can be reasonably reproduced, suggesting that M06-2X has an effective function for polycyclic reactions. However, the observed discrepancies in rate constants are inferred from dispersion effects and the multi-reference property in the computational system or derived from mismatch of target reactions between theoretical calculations and experiments. Overall, this study provides an insight into QCC application for investigating the formation mechanism of TPs in AOPs for removal of OMPs in water environment.
Collapse
Affiliation(s)
- Dhimas Dwinandha
- Department of Civil and Environmental Engineering, Tokyo Institute of Technology, Meguro-ku, Tokyo, 152-8552, Japan
| | - Bei Zhang
- Department of Civil and Environmental Engineering, Tokyo Institute of Technology, Meguro-ku, Tokyo, 152-8552, Japan
| | - Manabu Fujii
- Department of Civil and Environmental Engineering, Tokyo Institute of Technology, Meguro-ku, Tokyo, 152-8552, Japan.
| |
Collapse
|
24
|
Salarieh S, Golchoubian H, Moayyedi G, Rezaee E. Copper(II) complexes containing mixed- chelate ligand, solvatochromism and density functional theory studies. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.131822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
25
|
Geoghegan BL, Liu Y, Peredkov S, Dechert S, Meyer F, DeBeer S, Cutsail GE. Combining Valence-to-Core X-ray Emission and Cu K-edge X-ray Absorption Spectroscopies to Experimentally Assess Oxidation State in Organometallic Cu(I)/(II)/(III) Complexes. J Am Chem Soc 2022; 144:2520-2534. [PMID: 35050605 PMCID: PMC8855422 DOI: 10.1021/jacs.1c09505] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
A series of organometallic
copper complexes in formal oxidation
states ranging from +1 to +3 have been characterized by a combination
of Cu K-edge X-ray absorption (XAS) and Cu Kβ valence-to-core
X-ray emission spectroscopies (VtC XES). Each formal oxidation state
exhibits distinctly different XAS and VtC XES transition energies
due to the differences in the Cu Zeff, concomitant with
changes in physical oxidation state from +1 to +2 to +3. Herein, we
demonstrate the sensitivity of XAS and VtC XES to the physical oxidation
states of a series of N-heterocyclic carbene (NHC) ligated organocopper
complexes. We then extend these methods to the study of the [Cu(CF3)4]− ion. Complemented by computational
methods, the observed spectral transitions are correlated with the
electronic structure of the complexes and the Cu Zeff.
These calculations demonstrate that a contraction of the Cu 1s orbitals
to deeper binding energy upon oxidation of the Cu center manifests
spectroscopically as a stepped increase in the energy of both XAS
and Kβ2,5 emission features with increasing formal
oxidation state within the [Cun+(NHC2)]n+ series. The newly synthesized Cu(III) cation
[CuIII(NHC4)]3+ exhibits spectroscopic
features and an electronic structure remarkably similar to [Cu(CF3)4]−, supporting a physical oxidation
state assignment of low-spin d8 Cu(III) for [Cu(CF3)4]−. Combining XAS and VtC XES
further demonstrates the necessity of combining multiple spectroscopies
when investigating the electronic structures of highly covalent copper
complexes, providing a template for future investigations into both
synthetic and biological metal centers.
Collapse
Affiliation(s)
- Blaise L. Geoghegan
- Max Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, 45470 Mülheim an der Ruhr, Germany
- Institute of Inorganic Chemistry, University of Duisburg-Essen, Universitätsstrasse 5-7, 45117 Essen, Germany
| | - Yang Liu
- Institute of Inorganic Chemistry, University of Göttingen, Tammannstrasse 4, 37077 Göttingen, Germany
| | - Sergey Peredkov
- Max Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Sebastian Dechert
- Institute of Inorganic Chemistry, University of Göttingen, Tammannstrasse 4, 37077 Göttingen, Germany
| | - Franc Meyer
- Institute of Inorganic Chemistry, University of Göttingen, Tammannstrasse 4, 37077 Göttingen, Germany
| | - Serena DeBeer
- Max Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, 45470 Mülheim an der Ruhr, Germany
| | - George E. Cutsail
- Max Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, 45470 Mülheim an der Ruhr, Germany
- Institute of Inorganic Chemistry, University of Duisburg-Essen, Universitätsstrasse 5-7, 45117 Essen, Germany
| |
Collapse
|
26
|
Zhou K, Zhang YY, Chen XT, Hu SX. 12-Membered Ring Carbides with Stabilization of Actinide Atoms. Inorg Chem 2022; 61:2119-2128. [PMID: 35041785 DOI: 10.1021/acs.inorgchem.1c03341] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Actinide (Th and U) carbides as the potential nuclear fuels in nuclear reactors require basic research in order to understand the thermodynamic stability and performance of these substances. Here we report the structural characterization and bonding analyses of [C12], ThC12, and UC12 clusters via a global-minimum search combined with relativistic quantum chemistry calculations to elucidate the stability and bonding nature of An-C bonds. We predict that these [C12], ThC12, and UC12 compounds have a planar structure with C6h, D12h, and D12h symmetry, respectively. [C12] has a hyperconjugation structure containing alternating single and double bonds. The significant stabilization when forming AnC12 predominantly comes from the electrostatic interaction between An4+ and [C12]4- and also from a certain degree of orbital interaction between the An 5f6d7s valence shell and [C12] π orbitals. The covalent character of the An-C bonds exhibits a direct in-plane σ-type overlap of the C 2p-derived MOs of [C12] and the An 5fϕ AO, thus leading to an unprecedented electronic configuration of d1f1 for U in UC12. Our results present an example of the novel properties that can be expected for actinide compounds and would provide the knowledge required to obtain novel structures of AnC12 in future experiments.
Collapse
Affiliation(s)
- Ke Zhou
- College of Chemistry and Environmental Science & Shaanxi Key Laboratory of Catalysis & Institute of Theoretical and Computational Chemistry, Shaanxi University of Technology, Hanzhong 723000, China
| | - Yang-Yang Zhang
- Department of Chemistry & Key Laboratory of Organic Optoelectronics, Tsinghua University, Beijing 100084, China
| | - Xiao-Tong Chen
- Institute of Nuclear and New Energy Technology & Collaborative Innovation Center of Advanced Nuclear Energy Technology, Tsinghua University, Beijing 100084, China
| | - Shu-Xian Hu
- Department of Physics, University of Science and Technology Beijing, Beijing 100083, China.,Beijing Computational Science Research Center, Beijing 100193, China
| |
Collapse
|
27
|
Hu SX, Zhang P, Zhang P. Electronic Structures and Properties of Bimetallic Plutonium Group 13 Carbonyl Compounds [XPuCO] (X = B, Al, and Ga). Inorg Chem 2021; 60:18794-18803. [PMID: 34841875 DOI: 10.1021/acs.inorgchem.1c02503] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Bonding features of heterobimetallic complexes containing f-block elements are fundamental content in actinide chemistry. In order to account for the structural periodicity of the X-Pu carbonyls and the formation of chemical bonds between bimetallic plutonium and group 13 carbonyl compounds, we report a comprehensively quantum-chemical study of the electronic structure and properties of XPuCO (X = B, Al, and Ga). With increasing atomic radii of the group 13 elements, the XPuCO structure alternates from cyclic [PuCBO] to linear [AlCPuO] and [GaCPuO]. The bonding analysis indicates that the donor-acceptor model is the best description for bonding interactions of metal and ligands with different donation patterns of CBO → Pu and XC → PuO (X = Al and Ga). The apparent XC ← PuO backdonation increases the C-Pu bond strength markedly and stabilizes the linear geometry of [AlCPuO] and [GaCPuO], while spin-orbit coupling is found to be significant in the stabilization of [PuCBO]. The ground electron configurations and natural orbital analysis indicate that cyclic [PuCBO] and linear [XCPuO] (X = Al and Ga) are considered as complexes of Pu(III) and Pu(V), respectively. The trend presents a valuable insight for the 5f/6d-np bonding interactions, especially for the fundamental understanding of transuranic elements.
Collapse
Affiliation(s)
- Shu-Xian Hu
- School of Mathematics and Physics, University of Science and Technology Beijing, Beijing 100083, China.,Beijing Computational Science Research Center, Beijing 100193, China
| | - Peng Zhang
- Beijing Computational Science Research Center, Beijing 100193, China
| | - Ping Zhang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
| |
Collapse
|
28
|
Zhang P, Zou W, Zhang P, Hu S. Electronic Structures and Properties of Actinide‐Bimetal Compounds An
2
O
2
(An=Th to Cf) and U
2
E
2
(E=N, F, S). Eur J Inorg Chem 2021. [DOI: 10.1002/ejic.202100444] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Peng Zhang
- School of Mathematics and Physics University of Science and Technology Beijing Beijing 100083 China
| | - Wen‐Li Zou
- Institute of Modern Physics Northwest University Xi'an 710127 China
| | - Ping Zhang
- Institute of Applied Physics and Computational Mathematics Haidian Beijing 100088 China
| | - Shu‐Xian Hu
- School of Mathematics and Physics University of Science and Technology Beijing Beijing 100083 China
- Beijing Computational Science Research Center Haidian Beijing 100193 China
| |
Collapse
|
29
|
Gittus OR, Bresme F. Thermophysical properties of water using reactive force fields. J Chem Phys 2021; 155:114501. [PMID: 34551553 DOI: 10.1063/5.0057868] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The widescale importance and rich phenomenology of water continue to motivate the development of computational models. ReaxFF force fields incorporate many characteristics desirable for modeling aqueous systems: molecular flexibility, polarization, and chemical reactivity (bond formation and breaking). However, their ability to model the general properties of water has not been evaluated in detail. We present comprehensive benchmarks of the thermophysical properties of water for two ReaxFF models, the water-2017 and CHON-2017_weak force fields. These include structural, electrostatic, vibrational, thermodynamic, coexistence, and transport properties at ambient conditions (300 K and 0.997 g cm-3) and along the standard pressure (1 bar) isobar. Overall, CHON-2017_weak predicts more accurate thermophysical properties than the water-2017 force field. Based on our results, we recommend potential avenues for improvement: the dipole moment to quadrupole moment ratio, the self-diffusion coefficient, especially for water-2017, and the gas phase vibrational frequencies with the aim to improve the vibrational properties of liquid water.
Collapse
Affiliation(s)
- Oliver R Gittus
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| | - Fernando Bresme
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| |
Collapse
|
30
|
Pang J, Jin W, Kuang X, Lu C. Two-Dimensional Fe 8N Nanosheets: Ferromagnets and Nitrogen Diffusion. J Phys Chem Lett 2021; 12:8453-8459. [PMID: 34448584 DOI: 10.1021/acs.jpclett.1c02242] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We perform a first-principles study and identify two intriguing ferromagnets, hollow-Fe8N (H-Fe8N) and bridge-Fe8N (B-Fe8N) monolayers, by extensive structural searches. Both H-Fe8N and B-Fe8N nanosheets are buckled triangular lattices with a similar motif, but they are distinguishable by the positions of N atoms. The magnetic and electronic properties show that H-Fe8N is a low-spin ferromagnet; in contrast, B-Fe8N is a high-spin ferromagnet, which originates from the 3d orbital splitting of the Fe atom due to the low/high symmetric crystal field. Surprisingly, two stable Fe8N monolayers can be transferred to each other by N atom diffusion from the bridge position to the hollow position with the migration energy barrier of 1.5 eV. The energy barrier is affected by introduced Fe defects and rare earth metal dopants. These findings introduce a new tactic to regulate the 2D Fe-nitride monolayers at the atomic scale.
Collapse
Affiliation(s)
- Jiafei Pang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Wenyuan Jin
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Xiaoyu Kuang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Cheng Lu
- School of Mathematics and Physics, China University of Geosciences (Wuhan), Wuhan 430074, China
| |
Collapse
|
31
|
Computational design of binder as the LC3-p62 protein-protein interaction. Bioorg Chem 2021; 115:105241. [PMID: 34426157 DOI: 10.1016/j.bioorg.2021.105241] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2021] [Revised: 07/05/2021] [Accepted: 08/02/2021] [Indexed: 12/22/2022]
Abstract
Cellular autophagy is an intracellular degradation pathway, which transports damaged, deformed, senescent or non-functional proteins and organelles to lysosome for digestion and degradation. Cellular autophagy is deeply evolutionarily conservedfromyeasttomammaliancells, and many homologous proteins of the autophahgy regulators are found in several species. This physiological process maintains the steady state of cells. Furtheremore, autophagy dysfunction is closely related to various diseases, such as neurodegenerative diseases, inflammation-related diseases, cardiovascular diseases, metabolic diseases, etc. The LC3 and p62 protein protein interaction (PPI) promotes the formation of autophagosomes and delivers polyubiquitinated "cargoes" to autophagic degradation. Therefore, LC3-p62 PPI plays an integral role in the formation of autophagosomes and effectively inhibits autophagy. However, there are still few studies on the LC3-p62 PPI inhibitors for its unclear molecular mechanism. Furthermore, most of these inhibitors are macromolecules with poorly active, and small molecules are particularly scarce. In this article, the computation method was used to identify the hot spot and design peptides as the binder of LC3-p62 PPI. Findings from this work provide a reference for the follow-up research of discovering small molecule inhibitors targeting LC3-p62 PPI.
Collapse
|
32
|
Computational design of β-fluorinated morphine derivatives for pH-specific binding. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138723] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
33
|
Le HM, Guagliardo M, Gorden AEV, Clark AE. Ensemble effects on allylic oxidation within explicit solvation environments. Dalton Trans 2021; 50:9259-9268. [PMID: 34128517 DOI: 10.1039/d1dt00785h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Umbrella-sampling density functional theory molecular dynamics (DFT-MD) has been employed to study the full catalytic cycle of the allylic oxidation of cyclohexene using a Cu(ii) 7-amino-6-((2-hydroxybenzylidene)amino)quinoxalin-2-ol complex in acetonitrile to create cyclohexenone and H2O as products. After the initial H-atom abstraction step, two different reaction pathways have been identified that are distinguished by the participation of alkyl hydroperoxide (referred to as the "open" cycle) versus the methanol side-product (referred to as the "closed" cycle) within the catalyst recovery process. Importantly, both pathways involve dehydrogenation and re-hydrogenation of the -NH2 group bound to the Cu-site - a feature that is revealed from the ensemble sampling of configurations of the reactive species that are stabilized within the explicit solvent environment of the simulation. Estimation of the energy span from the experimental turnover frequency yields an approximate value of 22.7 kcal mol-1 at 350 K. Whereas the closed cycle value is predicted to be 26.2 kcal mol-1, the open cycle value at 16.5 kcal mol-1. Both pathways are further consistent with the equilibrium between Cu(ii) and Cu(iii) that has previously been observed. In comparison to prior static DFT calculations, the ensemble of both solute and solvent configurations has helped to reveal a breadth of processes that underpin the full catalytic cycle yielding a more comprehensive understanding of the importance of radical reactions and catalysis recovery.
Collapse
Affiliation(s)
- Hung M Le
- Department of Chemistry, Washington State University, Pullman, Washington 99164, USA.
| | - Mariano Guagliardo
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79430, USA
| | - Anne E V Gorden
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79430, USA
| | - Aurora E Clark
- Department of Chemistry, Washington State University, Pullman, Washington 99164, USA. and Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| |
Collapse
|
34
|
Ullah Z, Kim K, Venkanna A, Kim HS, Kim MI, Kim MH. Plausible Pnicogen Bonding of epi-Cinchonidine as a Chiral Scaffold in Catalysis. Front Chem 2021; 9:669515. [PMID: 34295874 PMCID: PMC8290064 DOI: 10.3389/fchem.2021.669515] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2021] [Accepted: 06/10/2021] [Indexed: 12/11/2022] Open
Abstract
As a non-covalent interaction of a chiral scaffold in catalysis, pnicogen bonding of epi-cinchonidine (epi-CD), a cinchona alkaloid, was simulated to consider whether the interaction can have the potential controlling enantiotopic face like hydrogen bonding. Among five reactive functional groups in epi-CD, two stable complexes of the hydroxyl group (X-epi-CD1) at C17 and of the quinoline ring (X-epi-CD2) at N16 with pnictide family analytes [X = substituted phosphine (PX), i.e., F, Br, Cl, CF3, CN, HO, NO2, and CH3, and pnictide family analytes, i.e., PBr3, BiI3, SbI3, and AsI3] were predicted with intermolecular interaction energies, charge transfer (QMulliken and QNBO), and band gap energies of HOMO-LUMO (Eg) at the B3LYP/6-31G(d,p) level of density functional theory. It was found that the dominant site of pnicogen bonding in epi-CD is the quinoline ring (N16 atom) rather than the hydroxyl group (O36 atom). In addition, the UV-Vis spectra of the complex were calculated by time-dependent density functional theory (TD-DFT) at the B3LYP/6-31+G(d,p) level and compared with experimental measurements. Through these calculations, two intermolecular interactions (H-bond vs. pnicogen bond) of epi-CD were compared.
Collapse
Affiliation(s)
- Zakir Ullah
- Department of Pharmacy, College of Pharmacy, Gachon Institute of Pharmaceutical Science, Gachon University, Incheon, South Korea.,Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, South Korea
| | - Kang Kim
- Department of Pharmacy, College of Pharmacy, Gachon Institute of Pharmaceutical Science, Gachon University, Incheon, South Korea
| | - Arramshetti Venkanna
- Department of Pharmacy, College of Pharmacy, Gachon Institute of Pharmaceutical Science, Gachon University, Incheon, South Korea
| | - Hye Su Kim
- Department of BioNano Technology, Gachon University, Seongnam, South Korea
| | - Moon Il Kim
- Department of BioNano Technology, Gachon University, Seongnam, South Korea
| | - Mi-Hyun Kim
- Department of Pharmacy, College of Pharmacy, Gachon Institute of Pharmaceutical Science, Gachon University, Incheon, South Korea
| |
Collapse
|
35
|
Ploetz EA, Karunaweera S, Bentenitis N, Chen F, Dai S, Gee MB, Jiao Y, Kang M, Kariyawasam NL, Naleem N, Weerasinghe S, Smith PE. Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20. J Chem Theory Comput 2021; 17:2964-2990. [PMID: 33878263 DOI: 10.1021/acs.jctc.1c00075] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.
Collapse
Affiliation(s)
- Elizabeth A Ploetz
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Sadish Karunaweera
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nikolaos Bentenitis
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Feng Chen
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Shu Dai
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Moon B Gee
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Yuanfang Jiao
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Myungshim Kang
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nilusha L Kariyawasam
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nawavi Naleem
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | | | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| |
Collapse
|
36
|
Jabłoński M. Study of Beryllium, Magnesium, and Spodium Bonds to Carbenes and Carbodiphosphoranes. Molecules 2021; 26:2275. [PMID: 33920004 PMCID: PMC8071025 DOI: 10.3390/molecules26082275] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Revised: 04/09/2021] [Accepted: 04/11/2021] [Indexed: 11/17/2022] Open
Abstract
The aim of this article is to present results of theoretical study on the properties of C⋯M bonds, where C is either a carbene or carbodiphosphorane carbon atom and M is an acidic center of MX2 (M = Be, Mg, Zn). Due to the rarity of theoretical data regarding the C⋯Zn bond (i.e., the zinc bond), the main focus is placed on comparing the characteristics of this interaction with C⋯Be (beryllium bond) and C⋯Mg (magnesium bond). For this purpose, theoretical studies (ωB97X-D/6-311++G(2df,2p)) have been performed for a large group of dimers formed by MX2 (X = H, F, Cl, Br, Me) and either a carbene ((NH2)2C, imidazol-2-ylidene, imidazolidin-2-ylidene, tetrahydropyrymid-2-ylidene, cyclopropenylidene) or carbodiphosphorane ((PH3)2C, (NH3)2C) molecule. The investigated dimers are characterized by a very strong charge transfer effect from either the carbene or carbodiphosphorane molecule to the MX2 one. This may even be over six times as strong as in the water dimer. According to the QTAIM and NCI method, the zinc bond is not very different than the beryllium bond, with both featuring a significant covalent contribution. However, the zinc bond should be definitely stronger if delocalization index is considered.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University, 87-100 Toruń, Poland
| |
Collapse
|
37
|
Bella G, Santoro A, Nicolò F, Bruno G, Cordaro M. Do Secondary Electrostatic Interactions Influence Multiple Dihydrogen Bonds? AA-DD Array on an Amine-Borane Aza-Coronand: Theoretical Studies and Synthesis. Chemphyschem 2021; 22:593-605. [PMID: 33432782 DOI: 10.1002/cphc.202000906] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2020] [Revised: 01/09/2021] [Indexed: 11/06/2022]
Abstract
Hydrogen bond plays a key role in a wide range of inorganic, organic, as well as biological systems. The understanding on how the chemical environment can affect this kind of interaction is crucial to predict its binding strength and consequently the robustness and the dynamic properties of many supramolecular systems. In this paper a new donor-acceptor complex was synthesized and characterized by SCXRD, showing for the first time in an organic system an AA-DD pattern of a particular hydrogen interaction, called dihydrogen bond. Over 250 functionals were computationally evaluated to select the best method to reproduce the binding interaction geometry of this new pattern. Moreover, a new vector force model was used to split the contribution of primary and secondary electrostatic interactions (SEIs), in order to evaluate how the latter one can modify the binding strength of this unusual hydrogen-hydrogen interaction.
Collapse
Affiliation(s)
- Giovanni Bella
- Department of chemical, biological, pharmaceutical and environmental sciences, University of Messina, Viale F. Stagno d'Alcontres 31, 98166, Messina, Italy
| | - Antonio Santoro
- Department of chemical, biological, pharmaceutical and environmental sciences, University of Messina, Viale F. Stagno d'Alcontres 31, 98166, Messina, Italy
| | - Francesco Nicolò
- Department of chemical, biological, pharmaceutical and environmental sciences, University of Messina, Viale F. Stagno d'Alcontres 31, 98166, Messina, Italy
| | - Giuseppe Bruno
- Department of chemical, biological, pharmaceutical and environmental sciences, University of Messina, Viale F. Stagno d'Alcontres 31, 98166, Messina, Italy
| | - Massimiliano Cordaro
- Department of chemical, biological, pharmaceutical and environmental sciences, University of Messina, Viale F. Stagno d'Alcontres 31, 98166, Messina, Italy
| |
Collapse
|
38
|
Nicholas KM. DFT analysis of the structures, dynamics and reactivity of [(Propargyl)Co2(CO)5(L)]+ Z−. J Organomet Chem 2021. [DOI: 10.1016/j.jorganchem.2020.121667] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
39
|
Synthesis, Structure and In Vitro Anticancer Activity of Pd(II) Complex of Pyrazolyl-s-Triazine Ligand; A New Example of Metal-Mediated Hydrolysis of s-Triazine Pincer Ligand. CRYSTALS 2021. [DOI: 10.3390/cryst11020119] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The square planar complex [Pd(PT)Cl(H2O)]*H2O (HPT: 6-(3,5-dimethyl-1H-pyrazol-1-yl)-1,3,5-triazine-2,4(1H,3H)-dione) was obtained by the reaction of 2-methoxy-4,6-bis(3,5-dimethyl-1H-pyrazol-1-yl)-1,3,5-triazine (MBPT) pincer ligand with PdCl2 in a molar ratio (1:1) under thermal conditions and using acetone as a solvent. The reaction proceeded via C-N cleavage of one C-N moiety that connects the pyrazole and s-triazine combined with the hydrolysis of the O-CH3 group. The reaction of the chloride salt of its higher congener (PtCl2) gave [Pt(3,5-dimethyl-1H-pyrazole)2Cl2]. The crystal structure of [Pd(PT)Cl(H2O)]*H2O complex is stabilized by inter- and intra-molecular hydrogen bonding interactions. Hirshfeld analysis revealed that the H...H (34.6%), O...H (23.6%), and Cl...H (7.8%) interactions are the major contacts in the crystal. The charges at Pd, H2O, Cl and PT are changed to 0.4995, 0.2216, −0.4294 and −0.2917 instead of +2, 0, −1 and −1, respectively, using the MPW1PW91 method. [Pd(PT)Cl(H2O)]*H2O complex has almost equal activities against MDA-MB-231 and MCF-7 cell lines with IC50 of 38.3 µg/mL.
Collapse
|
40
|
Kenouche S, Belkadi A, Djebaili R, Melkemi N. High regioselectivity in the amination reaction of isoquinolinequinone derivatives using conceptual DFT and NCI analysis. J Mol Graph Model 2021; 104:107828. [PMID: 33444977 DOI: 10.1016/j.jmgm.2020.107828] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2020] [Revised: 12/17/2020] [Accepted: 12/21/2020] [Indexed: 11/27/2022]
Abstract
DFT-derived reactivity descriptors and Non-Covalent interaction (NCI) analysis were performed to rationalize the regioselectivity in the amination reaction of some isoquinolinequinone derivatives. Statistical analysis was performed to assess robustness of atomic charges to the basis set. Various electronic population schemes including Mulliken population analysis (MPA), electrostatic method (ChelpG), Hirshfeld population analysis (HPA) and Natural population analysis (NPA) have been considered. The results revealed that NPA was the most efficient for this purpose. NCI study using the reduced density gradient (RDG) was performed for revealing weak interactions. Domains defined by isosurfaces of RDG have been integrated to quantitatively study the strength of weak interactions and their stabilities have also been examined. Steric hindrance caused by coordination of ethanol with the neighboring carbonyl prevents the nucleophilic attack on C-6 and therefore leads to preferential C-7 substitution. The quantitative study of NCI clearly demonstrates that the hydrogen bond of carbonyl (2) is more stabilizing. Consequently, the high polarity of hydrogen bond on this carbonyl may explain the high electrophilicity of C-5 compared to C-8. Our work proved that the difference in local reactivity, as well as the steric hindrance are the key elements explaining the high regioselectivity exhibited by the amination reaction.
Collapse
Affiliation(s)
- S Kenouche
- Research Team: Modeling of Chemical Systems Using Quantum Calculations, LCA Laboratory, University Mohamed Khider of Biskra, 07000, Biskra, Algeria.
| | - A Belkadi
- Research Team: Computational and Pharmaceutical Chemistry, LMCE Laboratory, University Mohamed Khider of Biskra, 07000, Biskra, Algeria
| | - R Djebaili
- Research Team: Computational and Pharmaceutical Chemistry, LMCE Laboratory, University Mohamed Khider of Biskra, 07000, Biskra, Algeria
| | - N Melkemi
- Research Team: Computational and Pharmaceutical Chemistry, LMCE Laboratory, University Mohamed Khider of Biskra, 07000, Biskra, Algeria
| |
Collapse
|
41
|
Han B, Isborn CM, Shi L. Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models. J Chem Theory Comput 2021; 17:889-901. [PMID: 33405925 DOI: 10.1021/acs.jctc.0c01102] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree-Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.
Collapse
Affiliation(s)
- Bowen Han
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
42
|
Curtolo F, Arantes GM. Mechanisms for Flavin-Mediated Oxidation: Hydride or Hydrogen-Atom Transfer? J Chem Inf Model 2020; 60:6282-6287. [DOI: 10.1021/acs.jcim.0c00945] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Felipe Curtolo
- Department of Biochemistry, Instituto de Quı́mica, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-900 São Paulo, SP, Brazil
| | - Guilherme M. Arantes
- Department of Biochemistry, Instituto de Quı́mica, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-900 São Paulo, SP, Brazil
| |
Collapse
|
43
|
Meng Q, Huang Y, Deng D, Yang Y, Sha H, Zou X, Faller R, Yuan Y, Zhu G. Porous Aromatic Framework Nanosheets Anchored with Lewis Pairs for Efficient and Recyclable Heterogeneous Catalysis. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2020; 7:2000067. [PMID: 33240750 PMCID: PMC7675047 DOI: 10.1002/advs.202000067] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/06/2020] [Revised: 09/10/2020] [Indexed: 05/14/2023]
Abstract
Lewis pairs (LPs) with outstanding performance for nonmetal-mediated catalysis reactions have high fundamental interest and remarkable application prospects. However, their solubility characteristics lead to instability and deactivation upon recycling. Here, the layered porous aromatic framework (PAF-6), featuring two kinds of Lewis base sites (NPiperazine and NTriazine), is exfoliated into few-layer nanosheets to form the LP entity with the Lewis acid. After comparison with various porous networks and verification by density functional theory (DFT) calculations, the NTriazine atom in the specific spatial environment is determined to preferably coordinate with the electron-deficient boron compound in a sterically hindered pattern. LP-bare porous product displays high catalytic activity for the hydrogenation of both olefin and imine compounds, and demonstrates ≈100% activity after 10 successful cycles in hydrogenation reactions. Considering the natural advantage of porous organic frameworks to construct LP groups opens up novel prospects for preparing other nonmetallic heterogeneous catalysts for efficient and recyclable catalysis.
Collapse
Affiliation(s)
- Qinghao Meng
- Key Laboratory of Polyoxometalate Science of Ministry of EducationNortheast Normal UniversityRenmin AvenueChangchun130024China
| | - Yihan Huang
- Department of Materials Science and EngineeringUniversity of CaliforniaDavisDavisCA95616USA
| | - Dan Deng
- Key Laboratory of Polyoxometalate Science of Ministry of EducationNortheast Normal UniversityRenmin AvenueChangchun130024China
| | - Yajie Yang
- Key Laboratory of Polyoxometalate Science of Ministry of EducationNortheast Normal UniversityRenmin AvenueChangchun130024China
| | - Haoyan Sha
- Department of Chemical EngineeringUniversity of CaliforniaDavisDavisCA95616USA
| | - Xiaoqin Zou
- Key Laboratory of Polyoxometalate Science of Ministry of EducationNortheast Normal UniversityRenmin AvenueChangchun130024China
| | - Roland Faller
- Department of Chemical EngineeringUniversity of CaliforniaDavisDavisCA95616USA
| | - Ye Yuan
- Key Laboratory of Polyoxometalate Science of Ministry of EducationNortheast Normal UniversityRenmin AvenueChangchun130024China
| | - Guangshan Zhu
- Key Laboratory of Polyoxometalate Science of Ministry of EducationNortheast Normal UniversityRenmin AvenueChangchun130024China
| |
Collapse
|
44
|
Bai S, Kubelka J, Piri M. Relationship between molecular charge distribution and wettability reversal efficiency of cationic surfactants on calcite surfaces. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.114009] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
45
|
Monajjemi M. ANALYSIS OF LOCALIZED ORBITALS IN AZABORA DERIVATIVES OF [8] ANNULENE: IN THE VIEWPOINT OF AROMATICITY AND INDUCED RING CURRENTS. J STRUCT CHEM+ 2020. [DOI: 10.1134/s0022476620100078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
46
|
Zhang P, Liu H, Zou W, Zhang P, Hu SX. Relativistic Effects Stabilize the Planar Wheel-like Structure of Actinide-Doped Gold Clusters: An@Au 7 (An = Th to Cm). J Phys Chem A 2020; 124:8173-8183. [PMID: 32845148 DOI: 10.1021/acs.jpca.0c02148] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Despite the chemistry of actinide-ligand bonding is continuing and of burgeoning interest, investigations of the chemical bonding of bimetallic complexes involving transuranics remain relatively less, and there are rarely studies on the bonding features between actinide and coinage metals (CM). We present a systematic research on the series of An@Au7 (An = Th to Cm), UCM7 (CM = Cu, Ag, Au), and WAu7 clusters to investigate the unique geometries, electronic structures, and chemical bonding between An 5f6d orbitals and CM ns orbitals, and to find their periodicity across the actinides and within the group of transition metals. A unique planar wheel-like structure for An@Au7 clusters with the help of actinide metals encapsulation via spin-orbit coupling, resulting in An(III). Instead, the transition-metal (TM) element W retains its usual six-gold-coordination structure in WAu7, thus forcing the seventh Au out of plane. The An-CM interactions, depending on the ion radii, become stronger with the increase of the atomic number of the actinide metals, as well as the CM. These results show that the presence of actinides in clusters can lead to unique electronic and geometrical structures.
Collapse
Affiliation(s)
- Peng Zhang
- Beijing Computational Science Research Center, Beijing 100193, China
| | - Haitao Liu
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
| | - Wenli Zou
- Institute of Modern Physics, Northwest University and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an 710127, China
| | - Ping Zhang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
| | - Shu-Xian Hu
- Beijing Computational Science Research Center, Beijing 100193, China
| |
Collapse
|
47
|
Lu J, Jin Y, Dong J, Wang T, Zhang H, Song Z, Feng B. Preparation of Pentacyclo‐[8.2.1.1
4,7
0
2,9
0
3,8
] ‐tetradecane‐5,6,11,12‐tetracarboxylic Dianhydride Combined with Gaussian Calculation. ChemistrySelect 2020. [DOI: 10.1002/slct.202002702] [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]
Affiliation(s)
- Jianqiang Lu
- School of Chemical Engineering Qingdao University of Science & Technology Qingdao 266042 PR China
| | - Yan Jin
- School of Chemical Engineering Qingdao University of Science & Technology Qingdao 266042 PR China
| | - Jingjing Dong
- School of Chemical Engineering Qingdao University of Science & Technology Qingdao 266042 PR China
| | - Tielin Wang
- School of Chemical Engineering Qingdao University of Science & Technology Qingdao 266042 PR China
| | - Hexing Zhang
- School of Chemical Engineering Qingdao University of Science & Technology Qingdao 266042 PR China
| | - Zhanqian Song
- School of Chemical Engineering Qingdao University of Science & Technology Qingdao 266042 PR China
| | - Baicheng Feng
- School of Chemical Engineering Qingdao University of Science & Technology Qingdao 266042 PR China
| |
Collapse
|
48
|
Sen A, Vyas N, Pandey B, Rajaraman G. Deciphering the mechanism of oxygen atom transfer by non-heme Mn IV-oxo species: an ab initio and DFT exploration. Dalton Trans 2020; 49:10380-10393. [PMID: 32613212 DOI: 10.1039/d0dt01785j] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Oxygen atom transfer (OAT) reactions employing transition metal-oxo species have tremendous significance in homogeneous catalysis for industrial use. Understanding the structural and mechanistic aspects of OAT reactions using high-valent metal-oxo species is of great importance to fine-tune their reactivity. Herein we examine the reactivity of a non-heme high-valent oxo-manganese(iv) complex, [MnIVH3buea(O)]- towards a variety of substrates such as PPh2Me, PPhMe2, PCy3, PPh3, and PMe3 using density functional theory as well as ab initio CASSCF/NEVPT2 methods. We have initially explored the structure and bonding of [MnIVH3buea(O)]- and its congener [MnIVH3buea(S)]-. Our calculations affirm an S = 3/2 ground state of the catalyst with the S = 5/2 and S = 1/2 excited states predicted to be too high lying in energy to participate in the reaction mechanism. Our ab initio CASSCF/NEVPT2 calculations, however, reveal a strong multi-reference character for the ground S = 3/2 state with many low-lying quartets mixing significantly with the ground state. This opens up various reaction channels, and the admixed wave-function evolves during the reaction with the excited triplet dominating the ground state wave-function at the reactant complex. Our calculations predict the following pattern of reactivity, PCy3 < PMe3 < PPh3 < PPhMe2 < PPh2Me for the OAT reaction with the MnIV[double bond, length as m-dash]O species which correlates well with the experimental observations. Detailed electronic structure analysis of the transitions states reveal that these substrates react via an unusual low-energy δ-type pathway where a spin-up electron from the substrate is transferred to the δ*x2-y2 orbital of the MnIV[double bond, length as m-dash]O facilitated by its multi-reference character. The unusual reactivity observed here has implications in understanding the reactivity of [Mn4Ca] species in photosystem II.
Collapse
Affiliation(s)
- Asmita Sen
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India.
| | | | | | | |
Collapse
|
49
|
Choudhuri I, Truhlar DG. Calculating and Characterizing the Charge Distributions in Solids. J Chem Theory Comput 2020; 16:5884-5892. [PMID: 32544328 DOI: 10.1021/acs.jctc.0c00440] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Accurate estimation of the partial atomic charges on metal centers is useful for understanding electronic and catalytic properties of materials. However, different methods of calculating these charges may give quite different results; this issue has been more widely studied for molecules than for solids. Here we study the charges on the metal centers of a test set of 18 solids containing transition metals by using density functional theory with several density functionals (PBE, PBE+U, TPSS, revTPSS, HLE17, revM06-L, B3LYP, B3LYP*, and other exchange-modified B3LYP functionals) and four charge models (Bader, Hirshfeld, CM5, and DDEC6). The test set contains 12 systems with nonmagnetic metal centers (eight metal oxides (MO2), two metal sulfides (MS2), and two metal selenides (MSe2)) and six ferromagnetic transition metal complexes. Our study shows that, among the four types of charges, Bader charges are the highest and Hirshfeld charges are the lowest for all the systems, regardless of the functional being used. The CM5 charges are bigger than DDEC6 charges for MX2 with M = Ti or Mo and X = S or Se, but for the other 14 cases they are lower. We found that the most of the systems are sensitive to the Hubbard U parameters in PBE+U and to the percentage X of Hartree-Fock exchange in exchange-modified B3LYP; as we increase U or X, the charges on the metal atoms in MX2 increase steadily. Testing different density functionals shows charges calculated with higher Hubbard U parameters in PBE+U are comparable to B3LYP (with 20% Hartree-Fock exchange). Among four meta-GGA functionals studied, the charges with HLE17 have the closest agreement with B3LYP. The variation of charges with choice of charge model is greater than the variation with choice of density functional.
Collapse
Affiliation(s)
- Indrani Choudhuri
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
50
|
Gu B, Cunningham B, Muñoz Santiburcio D, Da Pieve F, Artacho E, Kohanoff J. Efficient ab initio calculation of electronic stopping in disordered systems via geometry pre-sampling: Application to liquid water. J Chem Phys 2020; 153:034113. [PMID: 32716198 DOI: 10.1063/5.0014276] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Knowledge of the electronic stopping curve for swift ions, Se(v), particularly around the Bragg peak, is important for understanding radiation damage. Experimentally, however, the determination of such a feature for light ions is very challenging, especially in disordered systems such as liquid water and biological tissue. Recent developments in real-time time-dependent density functional theory (rt-TDDFT) have enabled the calculation of Se(v) along nm-sized trajectories. However, it is still a challenge to obtain a meaningful statistically averaged Se(v) that can be compared to observations. In this work, taking advantage of the correlation between the local electronic structure probed by the projectile and the distance from the projectile to the atoms in the target, we devise a trajectory pre-sampling scheme to select, geometrically, a small set of short trajectories to accelerate the convergence of the averaged Se(v) computed via rt-TDDFT. For protons in liquid water, we first calculate the reference probability distribution function (PDF) for the distance from the proton to the closest oxygen atom, ϕR(rp→O), for a trajectory of a length similar to those sampled experimentally. Then, short trajectories are sequentially selected so that the accumulated PDF reproduces ϕR(rp→O) to increasingly high accuracy. Using these pre-sampled trajectories, we demonstrate that the averaged Se(vp) converges in the whole velocity range with less than eight trajectories, while other averaging methods using randomly and uniformly distributed trajectories require approximately ten times the computational effort. This allows us to compare the Se(vp) curve to experimental data and assess widely used empirical tables based on Bragg's rule.
Collapse
Affiliation(s)
- Bin Gu
- Department of Physics, Nanjing University of Information Science and Technology, Nanjing 210044, China
| | - Brian Cunningham
- Atomistic Simulation Centre, Queen's University Belfast, Belfast BT71NN, Northern Ireland, United Kingdom
| | | | - Fabiana Da Pieve
- Royal Belgian Institute for Space Aeronomy, Av Circulaire 3, 1180 Brussels, Belgium
| | - Emilio Artacho
- CIC Nanogune BRTA, Tolosa Hiribidea 76, 20018 San Sebastian, Spain
| | - Jorge Kohanoff
- Atomistic Simulation Centre, Queen's University Belfast, Belfast BT71NN, Northern Ireland, United Kingdom
| |
Collapse
|