1
|
Kartashov SV, Fedonin AP, Fayzullin RR. Electronic Force Density Fields: Insights into Partial Bonds, Transition States, and Chemical Structure Evolution. J Phys Chem A 2024; 128:7471-7488. [PMID: 39186697 DOI: 10.1021/acs.jpca.4c03699] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/28/2024]
Abstract
This paper presents the quantum-topological binding approach, in which the electrostatic and total static force density fields, Fes(r) and F ( r ) , together with the electron density gradient field ∇ρ(r), are simultaneously analyzed to elucidate the chemical structure of transition states and the nature of interatomic interactions for semibroken semiformed partial chemical bonds. The approach attributes the discrepancies between the force fields F ( r ) and Fes(r) to the nonclassical electron-electron interaction effects. The internuclear gap between the zero-flux boundaries of Fes(r) and ∇ρ(r) indicates the interatomic charge transfer phenomenon (ICT) that occurs upon the formation of a system from free atoms. Concomitantly, the mismatch of the zero-flux surfaces defined in - F ( r ) and ∇ρ(r) can be interpreted as a phenomenon of the electron-transfer-induced quantum chemical response (QCR), which originates from the electron exchange correlation. Our study permits the assertion of parallels between partial bonds and noncovalent interactions, as both typically exhibit incomplete QCRs, indicating the partial electron sharing of the transferred density. The changes in atomic and pseudoatomic charges are employed to describe the evolution of the chemical structure upon the substitution reaction. It is observed that the acquired difference in the actual atomic electronegativity causes polarization upon the heterolytic breaking of virtually nonpolar bonds. It is further proposed that the proximity of closely related stationary states along the reaction path on a potential energy hypersurface implies their similarity in the manifestation of the ICT and sympathetic QCR. Furthermore, the involvement of an electron pair in a partial bond facilitates its delocalization through the attraction by the static forces F ( r ) and Fes(r) to a neighboring nucleus and through the smearing by the Pauli kinetic force FP(r).
Collapse
Affiliation(s)
- Sergey V Kartashov
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| | - Anton P Fedonin
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| | - Robert R Fayzullin
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| |
Collapse
|
2
|
Kartashov SV, Saifina AF, Fayzullin RR. Toward the Chemical Structure of Diborane: Electronic Force Density Fields, Effective Electronegativity, and Internuclear Turning Surface Properties. J Phys Chem A 2024; 128:7284-7297. [PMID: 39163513 DOI: 10.1021/acs.jpca.4c03492] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/22/2024]
Abstract
The chemical structure of diborane was elucidated through the superposition of the vector fields of the electron density gradient ∇ρ(r), the electrostatic force Fes(r), and the kinetic force Fk(r), together with the analysis of the cumulative charges of the atoms and pseudoatoms delimited in the aforementioned fields. It was proposed that the Fk-pseudoatomic charge could be employed as a metric for quantifying the ionic component of a related atomic charge. The electron permeability across an internuclear turning surface─specifically, the zero-flux surface in Fk(r)─was characterized by probing it through mapping the total static potential φem(r). The conceptualization of post hoc electronegativity was presented for consideration. Our analysis revealed that the ordinary B-H and bent B-μ-H bonds in diborane demonstrate the polar covalent character with minor contributions of the ionic component. The former bond exhibits a greater electron permeability through the internuclear turning surface, indicating a stronger tendency for electron sharing between the corresponding nucleus-dominated regions. The electron density accumulation along the bent B-μ-H bond path diverges from the minimum action trajectories of the forces Fes(r) and Fk(r). This phenomenon can be associated with the structural strain within the angled, three-center two-electron bonding B-μ-H-B. The internuclear B···B paths were identified in Fes(r) and Fk(r), in contrast to ∇ρ(r). This fact, in conjunction with a pronounced electron permeability through the mutual turning surface between the two boron nuclei, implies a certain degree of electron exchange between the boron-dominated pseudoatomic regions. Furthermore, the anomalous [B-]H···μ-H intermolecular polar interactions were described between the strongly negatively charged hydrogen atoms in the monoclinic crystalline β-phase of diborane. In fact, the ionic contributions to the charges of the hydrogen atoms are shown to be relatively small.
Collapse
Affiliation(s)
- Sergey V Kartashov
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| | - Alina F Saifina
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| | - Robert R Fayzullin
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan 420088, Russian Federation
| |
Collapse
|
3
|
Zhao G, Chung YG. PACMAN: A Robust Partial Atomic Charge Predicter for Nanoporous Materials Based on Crystal Graph Convolution Networks. J Chem Theory Comput 2024; 20:5368-5380. [PMID: 38822793 DOI: 10.1021/acs.jctc.4c00434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/03/2024]
Abstract
We report a fast and easy method (PACMAN) to assign partial atomic charges on metal-organic framework (MOF) and covalent-organic framework (COF) crystal structures based on graph convolution networks (GCNs) trained on >1.8 million high-fidelity partial atomic charge data obtained from the Quantum Metal-Organic Framework (QMOF) database. The developed model shows outstanding performance, achieving a mean absolute error (MAE) of 0.0055 e (test set performance) while maintaining consistency with DDEC6, Bader, and CM5 charges across diverse chemistry and topologies of MOFs and COFs. We find that the new method accurately assigns partial atomic charges for ion-containing nanoporous materials, which has not been possible in previous machine learning (ML) models. Grand canonical Monte Carlo (GCMC) simulation results for CO2 and N2 uptakes and the Widom particle insertion calculation for Henry's law constant of water results based on PACMAN and the original DDEC6 charges show excellent agreements compared to other ML models reported in the literature. The runtime analysis of the new method demonstrates that the partial atomic charges of MOF and COF structures with up to 500 atoms can be obtained in less than 10 s. An easy-to-use web interface has been developed to facilitate the adoption of the developed model.
Collapse
Affiliation(s)
- Guobin Zhao
- School of Chemical Engineering, Pusan National University, Busan 46241, South Korea
| | - Yongchul G Chung
- School of Chemical Engineering, Pusan National University, Busan 46241, South Korea
| |
Collapse
|
4
|
Bodo F, Erba A, Kraka E, Moura RT. Chemical bonding in Uranium-based materials: A local vibrational mode case study of Cs 2 UO 2 Cl 4 and UCl 4 crystals. J Comput Chem 2024; 45:1130-1142. [PMID: 38279637 DOI: 10.1002/jcc.27311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 12/10/2023] [Accepted: 01/08/2024] [Indexed: 01/28/2024]
Abstract
The Local Vibrational Mode Analysis, initially applied to diverse molecular systems, was extended to periodic systems in 2019. This work introduces an enhanced version of the LModeA software, specifically designed for the comprehensive analysis of two and three-dimensional periodic structures. Notably, a novel interface with the Crystal package was established, enabling a seamless transition from molecules to periodic systems using a unified methodology. Two distinct sets of uranium-based systems were investigated: (i) the evolution of the Uranyl ion (UO 2 2 + ) traced from its molecular configurations to the solid state, exemplified by Cs 2 UO 2 Cl 4 and (ii) Uranium tetrachloride (UCl 4 ) in both its molecular and crystalline forms. The primary focus was on exploring the impact of crystal packing on key properties, including IR and Raman spectra, structural parameters, and an in-depth assessment of bond strength utilizing local mode perspectives. This work not only demonstrates the adaptability and versatility of LModeA for periodic systems but also highlights its potential for gaining insights into complex materials and aiding in the design of new materials through fine-tuning.
Collapse
Affiliation(s)
- Filippo Bodo
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
- Dipartimento di Chimica, Università di Torino, Torino, Italy
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, Torino, Italy
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
| | - Renaldo T Moura
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, Texas, USA
- Department of Chemistry and Physics, Center of Agrarian Sciences, Federal University of Paraiba, Areia, Brazil
| |
Collapse
|
5
|
Richter WE. Symmetry-Constrained Properties Behave Differently for 2D or 3D Structures under the Same Point Group. J Phys Chem A 2024; 128:4308-4314. [PMID: 38758113 PMCID: PMC11145645 DOI: 10.1021/acs.jpca.4c02167] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2024] [Revised: 05/06/2024] [Accepted: 05/08/2024] [Indexed: 05/18/2024]
Abstract
In chemistry and physics, two molecules belonging to the same point group are expected to behave similarly regarding various properties, following their character tables. Here, we show that the derivative of the dipole moment with respect to the normal coordinate of vibration might show different symmetry constraints if the molecule is planar, even if these molecules belong to the same point group. Examples of pairs of molecules featuring these conditions are presented. These findings open a new path toward a much deeper understanding of how 2D materials behave so differently compared to 3D materials featuring the very same atoms and arrangements (graphene and graphite, for example); chemists and physicists dealing with 2D materials could benefit from looking more deeply into pure mathematical relations that might be governing 2D systems in a different way when compared to 3D systems. The aid from mathematicians is welcomed.
Collapse
Affiliation(s)
- Wagner Eduardo Richter
- Department of Chemistry, Federal
University of Technology − Paraná, Ponta Grossa, Paraná 84.017-220, Brazil
| |
Collapse
|
6
|
Barrera Y, Anderson JSM. Does the radical GPRI strongly depend on the population scheme? A comparative study to predict radical attack on unsaturated molecules with the radical general-purpose reactivity indicator. J Comput Chem 2024; 45:1152-1159. [PMID: 38299704 DOI: 10.1002/jcc.27314] [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/21/2023] [Revised: 07/15/2023] [Accepted: 01/09/2024] [Indexed: 02/02/2024]
Abstract
The reactivity of 22 unsaturated molecules undergoing attack by a methyl radical (⋅CH3) have been elucidated using the condensed radical general-purpose reactivity indicator (condensed radical GPRI) appropriate for relatively nucleophilic or electrophilic molecules. Using the appropriate radical GPRI equation for electrophilic attack or nucleophilic radical attack, seven different population schemes were used to assign the most reactive atoms in each of the 22 molecules. The results show that the condensed radical GPRI is sensitive to the population scheme chosen, but less sensitive than the radical Fukui function. Therefore, the reliability of these methods depends on the population scheme. Our investigation indicates that the condensed radical GPRI is most accurate in predicting the dominant products of the methyl radical addition reactions on a variety of unsaturated molecules when the Hirshfeld, Merz-Singh-Kollman, or Voronoi deformation density population schemes are used. Furthermore, for all populations schemes in the majority of instances where the radical Fukui function failed the radical GPRI was able to identify the most reactive atom under certain reactivity conditions.
Collapse
Affiliation(s)
- Yoshio Barrera
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Mexico City, Mexico
| | - James S M Anderson
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Mexico City, Mexico
| |
Collapse
|
7
|
Mehta N, Martin JML. On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment. J Comput Chem 2024; 45:1017-1032. [PMID: 38216516 DOI: 10.1002/jcc.27294] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2023] [Revised: 12/05/2023] [Accepted: 12/10/2023] [Indexed: 01/14/2024]
Abstract
Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges q on 1-particle basis set and (for WFT methods) n -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence.
Collapse
Affiliation(s)
- Nisha Mehta
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot, Israel
| |
Collapse
|
8
|
Bakó I, Pusztai L, Pothoczki S. Outstanding Properties of the Hydration Shell around β-d-Glucose: A Computational Study. ACS OMEGA 2024; 9:20331-20337. [PMID: 38737074 PMCID: PMC11080014 DOI: 10.1021/acsomega.4c00798] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/24/2024] [Revised: 04/10/2024] [Accepted: 04/12/2024] [Indexed: 05/14/2024]
Abstract
Ab initio molecular dynamics (AIMD) simulations have been performed on aqueous solutions of four simple sugars, α-d-glucose, β-d-glucose, α-d-mannose, and α-d-galactose. Hydrogen-bonding (HB) properties, such as the number of donor- and acceptor-type HB-s, and the lengths and strengths of hydrogen bonds between sugar and water molecules, have been determined. Related electronic properties, such as the dipole moments of water molecules and partial charges of the sugar O atoms, have also been calculated. The hydrophilic and hydrophobic shells were characterized by means of spatial distribution functions. β-d-Glucose was found to form the highest number of hydrophilic and the smallest number of hydrophobic connections to neighboring water molecules. The average sugar-water H-bond length was the shortest for β-d-glucose, which suggests that these are the strongest such H-bonds. Furthermore, β-d-glucose appears to stand out in terms of the symmetry properties of both its hydrophilic and hydrophobic hydration shells. In summary, in all aspects considered here, there seems to be a correlation between the distinct characteristics of β-d-glucose reported here and its outstanding solubility in water. Admittedly, our findings represent only some of the important factors that influence the solubility.
Collapse
Affiliation(s)
- Imre Bakó
- HUN-REN
Research Centre for Natural Sciences, Magyar tudósok körútja 2., H-1117 Budapest, Hungary
| | - László Pusztai
- HUN-REN
Wigner Research Centre for Physics, Konkoly-Thege M. út 29-33., H-1121 Budapest, Hungary
- International
Research Organization for Advanced Science and Technology (IROAST), Kumamoto University, 2-39-1 Kurokami, Chuo-ku, Kumamoto 860-8555, Japan
| | - Szilvia Pothoczki
- HUN-REN
Wigner Research Centre for Physics, Konkoly-Thege M. út 29-33., H-1121 Budapest, Hungary
| |
Collapse
|
9
|
Gallegos M, Isamura BK, Popelier PLA, Martín Pendás Á. An Unsupervised Machine Learning Approach for the Automatic Construction of Local Chemical Descriptors. J Chem Inf Model 2024; 64:3059-3079. [PMID: 38498942 PMCID: PMC11040729 DOI: 10.1021/acs.jcim.3c01906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2023] [Revised: 03/06/2024] [Accepted: 03/07/2024] [Indexed: 03/20/2024]
Abstract
Condensing the many physical variables defining a chemical system into a fixed-size array poses a significant challenge in the development of chemical Machine Learning (ML). Atom Centered Symmetry Functions (ACSFs) offer an intuitive featurization approach by means of a tedious and labor-intensive selection of tunable parameters. In this work, we implement an unsupervised ML strategy relying on a Gaussian Mixture Model (GMM) to automatically optimize the ACSF parameters. GMMs effortlessly decompose the vastness of the chemical and conformational spaces into well-defined radial and angular clusters, which are then used to build tailor-made ACSFs. The unsupervised exploration of the space has demonstrated general applicability across a diverse range of systems, spanning from various unimolecular landscapes to heterogeneous databases. The impact of the sampling technique and temperature on space exploration is also addressed, highlighting the particularly advantageous role of high-temperature Molecular Dynamics (MD) simulations. The reliability of the resulting features is assessed through the estimation of the atomic charges of a prototypical capped amino acid and a heterogeneous collection of CHON molecules. The automatically constructed ACSFs serve as high-quality descriptors, consistently yielding typical prediction errors below 0.010 electrons bound for the reported atomic charges. Altering the spatial distribution of the functions with respect to the cluster highlights the critical role of symmetry rupture in achieving significantly improved features. More specifically, using two separate functions to describe the lower and upper tails of the cluster results in the best performing models with errors as low as 0.006 electrons. Finally, the effectiveness of finely tuned features was checked across different architectures, unveiling the superior performance of Gaussian Process (GP) models over Feed Forward Neural Networks (FFNNs), particularly in low-data regimes, with nearly a 2-fold increase in prediction quality. Altogether, this approach paves the way toward an easier construction of local chemical descriptors, while providing valuable insights into how radial and angular spaces should be mapped. Finally, this work opens the possibility of encoding many-body information beyond angular terms into upcoming ML features.
Collapse
Affiliation(s)
- Miguel Gallegos
- Department
of Analytical and Physical Chemistry, University
of Oviedo, Oviedo E-33006, Spain
| | | | - Paul L. A. Popelier
- Department
of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, U.K.
| | - Ángel Martín Pendás
- Department
of Analytical and Physical Chemistry, University
of Oviedo, Oviedo E-33006, Spain
| |
Collapse
|
10
|
Xu C, Bao K, Wang S, Wu G, Ma S, Li L, Chu PK, Liu C. First-principles study of multifunctional Mn 2B 3 materials with high hardness and ferromagnetism. Phys Chem Chem Phys 2024; 26:12009-12015. [PMID: 38576324 DOI: 10.1039/d3cp04414a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/06/2024]
Abstract
Transition metal boride TM2B3 is widely studied in the field of physics and materials science. However, Mn2B3 has not been found in Mn-B systems so far. Mn2B3 undergoes phase transitions from Cmcm (0-28 GPa) to C2/m (28-80 GPa) and finally to C2/c (80-200 GPa) under pressure. Among these stable phases, Cmcm- and C2/m-Mn2B3s comprise six-membered boron rings and C2/c-Mn2B3 has wavy boron chains. They all have good mechanical properties and can become potential multifunctional materials. The strong B-B covalent bonding is mainly responsible for the structural stability and hardness. Comparison of the hardness of the five TM2B3s with different bonding strengths of TM-B and B-B bonds reveals a nonlinear change in the hardness. According to the Stoner model, these structures possess ferromagnetism, and the corresponding magnetic moments are almost the same as those of GGA and GGA + U (U = 3.9 eV, J = 1 eV).
Collapse
Affiliation(s)
- Chunhong Xu
- School of Physics and Electronic Engineering, Northeast Petroleum University, Daqing, 163318, China.
| | - Kuo Bao
- State Key Laboratory of Superhard Materials, College of Physics, Jilin University, Changchun, 130012, People's Republic of China.
| | - Sheng Wang
- School of Physics and Electronic Engineering, Northeast Petroleum University, Daqing, 163318, China.
| | - Gang Wu
- School of Physics and Electronic Engineering, Northeast Petroleum University, Daqing, 163318, China.
| | - Shuailing Ma
- Institute of High Pressure Physics, School of Physical Scientific and Technology, Ningbo University, Ningbo, 315211, China.
| | - Liangliang Li
- School of Physics and Electronic Engineering, Northeast Petroleum University, Daqing, 163318, China.
| | - Paul K Chu
- Department of Physics, Department of Materials Science and Engineering, and Department of Biomedical Engineering, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong, China
| | - Chao Liu
- School of Physics and Electronic Engineering, Northeast Petroleum University, Daqing, 163318, China.
| |
Collapse
|
11
|
Martins Filho PEC, Haiduke RLA. A Charge-Charge Flux-Dipole Flux Analysis of Simple Molecular Systems with Halogen Bonds. J Phys Chem A 2024; 128:2058-2071. [PMID: 38457530 DOI: 10.1021/acs.jpca.3c08229] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/10/2024]
Abstract
The presence of halogen bonds (R-X···B; R = substituent group, X = halogen, and B = Lewis base) provides quite amazing molecular systems for electronic structure investigations, presenting unique characteristics of fundamental relevance to supramolecular chemistry among other areas. Here, we use a double-hybrid approach from Density Functional Theory and triple-ζ basis sets augmented with diffuse functions (B2PLYP/def2-TZVPD) to deal with a large group of simple molecular systems containing halogen bonds (XBs), focusing on geometrical structures, binding energies, harmonic vibrational frequencies, and fundamental infrared intensities. Next, the electron densities and their variations on vibrations are carefully studied with the Quantum Theory of Atoms in Molecules (QTAIM) formalism and the charge-charge flux-dipole flux (CCFDF) model. We notice that the R-X stretching mode usually shows vibrational frequency decrements and infrared intensifications during the XB formation. Such features were also observed in hydrogen bonds, although the explanation for the band strengthening is different. Surprisingly, the most important contribution to these intensity increments due to complexation is now the interaction term between the charge flux and dipole flux (CF × DF). Thus, the use of atomic dipoles is mandatory to fully understand this phenomenon. In fact, the huge charge flux contributions to changes in dipole moment derivatives of R-X stretchings on halogen bonding are no longer accompanied by opposite variations of similar magnitudes in polarizations described by atomic dipole fluxes, which provided nearly unaltered values during the XB formation. Thus, the electronic charge flux direction change that takes place in complexes (from B to R) now reinforces dipole moment derivative terms from such atomic polarizations (mainly from the X atom). This intermolecular charge flux seems to be responsible for the unusual features noticed in the R-X stretching mode with the CCDDF/QTAIM model.
Collapse
Affiliation(s)
| | - Roberto Luiz Andrade Haiduke
- Departamento de Química e Física Molecular, Instituto de Química de São Carlos, Universidade de São Paulo, CP 780, São Carlos, SP 13560-970, Brazil
| |
Collapse
|
12
|
Zhang E, Hirao H. Synergistic Charge Transfer Effect in Ferrous Heme-CO Bonding within Cytochrome P450. Molecules 2024; 29:873. [PMID: 38398625 PMCID: PMC10891957 DOI: 10.3390/molecules29040873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2024] [Revised: 02/05/2024] [Accepted: 02/13/2024] [Indexed: 02/25/2024] Open
Abstract
We conducted ab initio valence bond (VB) calculations employing the valence bond self-consistent field (VBSCF) and breathing orbital valence bond (BOVB) methods to investigate the nature of the coordination bonding between ferrous heme and carbon monoxide (CO) within cytochrome P450. These calculations revealed the significant influence exerted by both proximal and equatorial ligands on the π-backdonation effect from the heme to the CO. Moreover, our VB calculations unveiled a phenomenon of synergistic charge transfer (sCT). In the case of ferrous heme-CO bonding, the significant stabilization in this sCT arises from cooperative resonance between the VB structures associated with σ donation and π backdonation. Unlike many other ligands, CO possesses the unique ability to establish two mutually perpendicular π-backdonation orbital interaction pairs, leading to an intensified stabilization attributed to σ-π resonance. Furthermore, while of a smaller energy magnitude, sCT due to one π-π pair is also present, contributing to the differential stabilization of ferrous heme-CO bonding.
Collapse
Affiliation(s)
| | - Hajime Hirao
- Warshel Institute for Computational Biology, School of Medicine, The Chinese University of Hong Kong, Shenzhen 518172, China;
| |
Collapse
|
13
|
Arismendi-Arrieta DJ, Sen A, Eriksson A, Broqvist P, Kullgren J, Hermansson K. H2O2(s) and H2O2·2H2O(s) crystals compared with ices: DFT functional assessment and D3 analysis. J Chem Phys 2023; 159:194701. [PMID: 37966002 DOI: 10.1063/5.0145203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2023] [Accepted: 07/27/2023] [Indexed: 11/16/2023] Open
Abstract
The H2O and H2O2 molecules resemble each other in a multitude of ways as has been noted in the literature. Here, we present density functional theory (DFT) calculations for the H2O2(s) and H2O2·2H2O(s) crystals and make selected comparisons with ice polymorphs. The performance of a number of dispersion-corrected density functionals-both self-consistent and a posteriori ones-are assessed, and we give special attention to the D3 correction and its effects. The D3 correction to the lattice energies is large: for H2O2(s) the D3 correction constitutes about 25% of the lattice energy using PBE, much more for RPBE, much less for SCAN, and it primarily arises from non-H-bonded interactions out to about 5 Å.The large D3 corrections to the lattice energies are likely a consequence of several effects: correction for missing dispersion interaction, the ability of D3 to capture and correct various other kinds of limitations built into the underlying DFT functionals, and finally some degree of cell-contraction-induced polarization enhancement. We find that the overall best-performing functionals of the twelve examined are optPBEvdW and RPBE-D3. Comparisons with DFT assessments for ices in the literature show that where the same methods have been used, the assessments largely agree.
Collapse
Affiliation(s)
| | - Anik Sen
- Department of Chemistry-Ångström, Uppsala University, P.O. Box 530, S-75121 Uppsala, Sweden
| | - Anders Eriksson
- Department of Chemistry-Ångström, Uppsala University, P.O. Box 530, S-75121 Uppsala, Sweden
| | - Peter Broqvist
- Department of Chemistry-Ångström, Uppsala University, P.O. Box 530, S-75121 Uppsala, Sweden
| | - Jolla Kullgren
- Department of Chemistry-Ångström, Uppsala University, P.O. Box 530, S-75121 Uppsala, Sweden
| | - Kersti Hermansson
- Department of Chemistry-Ångström, Uppsala University, P.O. Box 530, S-75121 Uppsala, Sweden
| |
Collapse
|
14
|
Benbella A, Jabraoui H, Matrane I, Mazroui M. Exploring adsorption behavior of sulfur and nitrogen compounds on transition metal-doped Cu(100) surfaces: insights from DFT and MD simulations. Phys Chem Chem Phys 2023; 25:27553-27565. [PMID: 37807806 DOI: 10.1039/d3cp04379g] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/10/2023]
Abstract
We conducted an extensive investigation using density functional theory (DFT) calculations and ReaxFF molecular dynamics (MD) simulations to elucidate the mechanisms of desulfurization and denitrogenation on Cu(100) surfaces. This study encompassed both pristine surfaces and those modified with Pt or Rh transition metals. Our primary objective was to gain a deep understanding of the adsorption behavior of thiophene (C4H4S) and pyridine (C5H5N) molecules on stepped Cu(100) surfaces, which serve as models for sulfur and nitrogen compounds. We systematically explored the interplay among water, adsorption efficiency, and surface regeneration capabilities. Using DFT, we thoroughly examined various aspects, including interaction energies, charge transfers, changes in electron density, and alterations in work function upon molecule adsorption. Notably, we observed a decrease in the interaction energy of thiophene, whereas that of pyridine increased when adsorbed on Pt/Rh-doped surfaces compared to pristine ones. Thiophene adsorption reduced the work function, potentially enhancing detectability, without causing inhibitory effects on any surface. Stepped Cu(100) surfaces demonstrated a strong affinity for thiophene, exhibiting an energy difference of approximately 86 kJ mol-1. However, this trend reversed on doped surfaces, where pyridine displayed stronger adsorption than thiophene, resulting in energy differences of around 123 kJ mol-1 and 62 kJ mol-1 on Pt-Cu and Rh-Cu surfaces, respectively. Moreover, our investigation highlighted the regeneration capacity of these surfaces, indicating that all surfaces can be considered promising candidates for desulfurization, while only Cu and Pt-Cu surfaces were found to be suitable for denitrogenation. Furthermore, results from MD simulations in combination with potential of mean force (PMF) simulations at 300 K, aligned with DFT calculations, confirmed the adsorption configurations of pyridine and thiophene. This analysis demonstrated the competitive advantage of thiophene over pyridine in adsorption and highlighted the inhibitory effect of water on pyridine adsorption on the Cu(100) surface.
Collapse
Affiliation(s)
- Achraf Benbella
- Laboratoire de Physique de la Matiere Condensee, Faculte des Sciences Ben Msik, Hassan II University of Casablanca, Casablanca B.P. 7955, Morocco.
| | | | - Imane Matrane
- Laboratoire de Physique de la Matiere Condensee, Faculte des Sciences Ben Msik, Hassan II University of Casablanca, Casablanca B.P. 7955, Morocco.
| | - M'hammed Mazroui
- Laboratoire de Physique de la Matiere Condensee, Faculte des Sciences Ben Msik, Hassan II University of Casablanca, Casablanca B.P. 7955, Morocco.
| |
Collapse
|
15
|
Hirao H, Zhang E. Bidirectional Charge Transfer at the Heme Iron in Reversible and Quasi-irreversible Cytochrome P450 Inhibition. Inorg Chem 2023; 62:16599-16608. [PMID: 37737847 DOI: 10.1021/acs.inorgchem.3c02541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/23/2023]
Abstract
The coordination bonding between inhibitor ligands and heme iron plays a critical role in disrupting the essential catalytic functions of cytochrome P450 enzymes (P450s). Despite its intrinsic importance and consequential implications for human health, our current understanding of coordination bonding in P450 inhibition remains limited. To address this knowledge gap, we conducted a systematic theoretical analysis of the complexes between a ferric or a ferrous heme model and representative inhibitor ligands. Specifically, we evaluated the charge-transfer (CT) effect within these complexes by employing a series of theoretical methods based on density functional theory (DFT). Through a comprehensive analysis, we unveiled the relative significance of ligand-to-heme forward CT in the ferric and ferrous complexes of reversible inhibitors. In contrast, backward CT dominates over forward CT in the ferrous heme complexes of quasi-irreversible inhibitors. Further analysis using the compact frontier orbital method underscores the elevated electron-accepting abilities of quasi-irreversible inhibitors for π backdonation, which greatly amplifies their binding affinity for the ferrous heme. This study sheds light on the intricate mechanisms underlying P450 inhibition and provides valuable insights for future inhibitor design and development.
Collapse
Affiliation(s)
- Hajime Hirao
- Warshel Institute for Computational Biology, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Shenzhen, Guangdong 518172, P. R. China
| | - Enhua Zhang
- Warshel Institute for Computational Biology, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Shenzhen, Guangdong 518172, P. R. China
| |
Collapse
|
16
|
Zhu Q, Wu Y, Zhao S, Cieplak P, Duan Y, Luo R. Streamlining and Optimizing Strategies of Electrostatic Parameterization. J Chem Theory Comput 2023; 19:6353-6365. [PMID: 37676646 PMCID: PMC10530599 DOI: 10.1021/acs.jctc.3c00659] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
Accurate characterization of electrostatic interactions is crucial in molecular simulation. Various methods and programs have been developed to obtain electrostatic parameters for additive or polarizable models to replicate electrostatic properties obtained from experimental measurements or theoretical calculations. Electrostatic potentials (ESPs), a set of physically well-defined observables from quantum mechanical (QM) calculations, are well suited for optimization efforts due to the ease of collecting a large amount of conformation-dependent data. However, a reliable set of QM ESP computed at an appropriate level of theory and atomic basis set is necessary. In addition, despite the recent development of the PyRESP program for electrostatic parameterizations of induced dipole-polarizable models, the time-consuming and error-prone input file preparation process has limited the widespread use of these protocols. This work aims to comprehensively evaluate the quality of QM ESPs derived by eight methods, including wave function methods such as Hartree-Fock (HF), second-order Møller-Plesset (MP2), and coupled cluster-singles and doubles (CCSD), as well as five hybrid density functional theory (DFT) methods, used in conjunction with 13 different basis sets. The highest theory levels CCSD/aug-cc-pV5Z (a5z) and MP2/aug-cc-pV5Z (a5z) were selected as benchmark data over two homemade data sets. The results show that the hybrid DFT method, ωB97X-D, combined with the aug-cc-pVTZ (a3z) basis set, performs well in reproducing ESPs while taking both accuracy and efficiency into consideration. Moreover, a flexible and user-friendly program called PyRESP_GEN was developed to streamline input file preparation. The restraining strengths, along with strategies for polarizable Gaussian multipole (pGM) model parameterizations, were also optimized. These findings and the program presented in this work facilitate the development and application of induced dipole-polarizable models, such as pGM models, for molecular simulations of both chemical and biological significance.
Collapse
Affiliation(s)
- Qiang Zhu
- Department of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, Irvine, California 92697, United States
| | - Yongxian Wu
- Department of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, Irvine, California 92697, United States
| | - Shiji Zhao
- Nurix Therapeutics, Inc., 1700 Owens St, Suite 205, San Francisco, California 94158, United States
| | - Piotr Cieplak
- SBP Medical Discovery Institute, La Jolla, California 92037, United States
| | - Yong Duan
- UC Davis Genome Center and Department of Biomedical Engineering, University of California, Davis, Davis, California 95616, United States
| | - Ray Luo
- Department of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, Irvine, California 92697, United States
| |
Collapse
|
17
|
Brinck T, Sahoo SK. Anomalous π-backbonding in complexes between B(SiR 3) 3 and N 2: catalytic activation and breaking of scaling relations. Phys Chem Chem Phys 2023; 25:21006-21019. [PMID: 37519222 DOI: 10.1039/d3cp00248a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/01/2023]
Abstract
Chemical transformations of molecular nitrogen (N2), including the nitrogen reduction reaction (NRR), are difficult to catalyze because of the weak Lewis basicity of N2. In this study, it is shown that Lewis acids of the types B(SiR3)3 and B(GeR3)3 bind N2 and CO with anomalously short and strong B-N or B-C bonds. B(SiH3)3·N2 has a B-N bond length of 1.48 Å and a complexation enthalpy of -15.9 kcal mol-1 at the M06-2X/jun-cc-pVTZ level. The selective binding enhancement of N2 and CO is due to π-backbonding from Lewis acid to Lewis base, as demonstrated by orbital analysis and density difference plots. The π-backbonding is found to be a consequence of constructive orbital interactions between the diffuse and highly polarizable B-Si and B-Ge bond regions and the π and π* orbitals of N2. This interaction is strengthened by electron donating substituents on Si or Ge. The π-backbonding interaction is predicted to activate N2 for chemical transformation and reduction, as it decreases the electron density and increases the length of the N-N bond. The binding of N2 and CO by the B(SiR3)3 and B(GeR3)3 types of Lewis acids also has a strong σ-bonding contribution. The relatively high σ-bond strength is connected to the highly positive surface electrostatic potential [VS(r)] above the B atom in the tetragonal binding conformation, but the σ-bonding also has a significant coordinate covalent (dative) contribution. Electron withdrawing substituents increase the potential and the σ-bond strength, but favor the binding of regular Lewis acids, such as NH3 and F-, more strongly than binding of N2 and CO. Molecules of the types B(SiR3)3 and B(GeR3)3 are chemically labile and difficult to synthesize. Heterogenous catalysts with the wanted B(Si-)3 or B(Ge-)3 bonding motif may be prepared by boron doping of nanostructured silicon or germanium compounds. B-doped and hydrogenated silicene is found to have promising properties as catalyst for the electrochemical NRR.
Collapse
Affiliation(s)
- Tore Brinck
- Department of Chemistry, CBH, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden.
| | - Suman Kalyan Sahoo
- Department of Chemistry, CBH, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden.
| |
Collapse
|
18
|
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
|
19
|
Marques HM. The inorganic chemistry of the cobalt corrinoids - an update. J Inorg Biochem 2023; 242:112154. [PMID: 36871417 DOI: 10.1016/j.jinorgbio.2023.112154] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2023] [Revised: 01/23/2023] [Accepted: 01/26/2023] [Indexed: 02/05/2023]
Abstract
The inorganic chemistry of the cobalt corrinoids, derivatives of vitamin B12, is reviewed, with particular emphasis on equilibrium constants for, and kinetics of, their axial ligand substitution reactions. The role the corrin ligand plays in controlling and modifying the properties of the metal ion is emphasised. Other aspects of the chemistry of these compounds, including their structure, corrinoid complexes with metals other than cobalt, the redox chemistry of the cobalt corrinoids and their chemical redox reactions, and their photochemistry are discussed. Their role as catalysts in non-biological reactions and aspects of their organometallic chemistry are briefly mentioned. Particular mention is made of the role that computational methods - and especially DFT calculations - have played in developing our understanding of the inorganic chemistry of these compounds. A brief overview of the biological chemistry of the B12-dependent enzymes is also given for the reader's convenience.
Collapse
Affiliation(s)
- Helder M Marques
- Molecular Sciences Institute, School of Chemistry, University of the Witwatersrand, Johannesburg 2050, South Africa.
| |
Collapse
|
20
|
Richter WE. On the quest for a molecular vibration whose absolute intensity is described solely by fluctuations of atomic dipoles: Can we find it? SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 291:122321. [PMID: 36621029 DOI: 10.1016/j.saa.2023.122321] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/22/2022] [Revised: 12/30/2022] [Accepted: 01/02/2023] [Indexed: 06/17/2023]
Abstract
The conditions for a molecular vibration to be active in the infrared spectrum while its absolute intensity is described solely by fluctuations of atomic dipoles are presented. A quick recipe of features of such a system is presented, to guide the seek for it. If such a system can be found, it will then be proved that population analyses based solely on a point-charge approximation are unrealistic as they cannot properly describe this intensity, a real, measurable, unambiguous property.
Collapse
Affiliation(s)
- Wagner Eduardo Richter
- Department of Chemistry, Federal University of Technology - Paraná, 84.017-220, Ponta Grossa, PR, Brazil.
| |
Collapse
|
21
|
Liu S, Xia S, Yue D, Sun H, Hirao H. The Bonding Nature of Fe–CO Complexes in Heme Proteins. Inorg Chem 2022; 61:17494-17504. [DOI: 10.1021/acs.inorgchem.2c02387] [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]
Affiliation(s)
- Shuyang Liu
- Warshel Institute for Computational Biology, School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Guangdong518172, P. R. China
| | - Songyan Xia
- Warshel Institute for Computational Biology, School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Guangdong518172, P. R. China
| | - Dongxiao Yue
- Warshel Institute for Computational Biology, School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Guangdong518172, P. R. China
| | - Haoran Sun
- Warshel Institute for Computational Biology, School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Guangdong518172, P. R. China
| | - Hajime Hirao
- Warshel Institute for Computational Biology, School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen, Guangdong518172, P. R. China
| |
Collapse
|
22
|
Theoretical insights into the nature of the bonding between carbon monoxide and iron(II) phthalocyanine: How do QTAIM descriptors change as a function of the Fe–CO distance? Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.139901] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
23
|
Richter WE, Duarte LJ. Infrared intensities of [Formula: see text]: a true challenge for DFT methods. J Mol Model 2022; 28:301. [PMID: 36066643 DOI: 10.1007/s00894-022-05275-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 08/11/2022] [Indexed: 10/14/2022]
Abstract
Absolute infrared intensities of [Formula: see text] were evaluated with a great variety of DFT and ab initio methods and basis sets. It is shown that the intensities calculated by different levels of theory may not agree with each other even in the qualitative (weak/strong) sense. Geometrical parameters, as well as net atomic charges evaluated from multiple partition schemes, did not vary as much as the intensities and thus cannot explain the tremendous differences found for the latter. As there are no experimental estimates for the intensities to guide the theoretical evaluation, deciding on the best level of theory is reduced to comparisons between the different DFT methods and QCISD or CCSD, believed to be the best theoretical estimates in the set. The differences found among the various DFT methods suggest the development of new methods, instead of converging to a focal point, is rather diverging.
Collapse
Affiliation(s)
- Wagner E Richter
- Department of Chemistry, Federal University of Technology - Paraná [UTFPR], 84.017-220, Ponta Grossa, PR, Brazil.
| | - Leonardo J Duarte
- Institute of Chemistry, State University of Campinas [UNICAMP], 13.083-970, Campinas, SP, Brazil
| |
Collapse
|
24
|
Richter WE, Duarte LJ, Bruns RE. Unavoidable failure of point charge descriptions of electronic density changes for out-of-plane distortions. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2022; 271:120891. [PMID: 35085995 DOI: 10.1016/j.saa.2022.120891] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Revised: 01/06/2022] [Accepted: 01/10/2022] [Indexed: 06/14/2023]
Abstract
Population analyses based on point charge approximations accurately estimating the equilibrium dipole moment will systematically fail when predicting infrared intensities of out-of-plane vibrations of planar molecules, whereas models based on both charges and dipoles will always succeed. It is not a matter of how the model is devised but rather how many degrees of freedom are available for the calculation. Population analyses based on point charges are very limited in terms of the amount of meaningful chemical information they provide, whereas models employing both atomic charges and atomic dipoles should be preferred for molecular distortions. A good model should be able to correctly describe not only static, equilibrium structures but also distorted geometries in order to correctly assess information from vibrating molecules. The limitations of point charge models also hold for distortions much larger than those encountered vibrationally.
Collapse
Affiliation(s)
- Wagner E Richter
- Department of Chemical Engineering, Federal University of Technology - Paraná, Postal Code 84.017-220, Ponta Grossa, Paraná, Brazil.
| | - Leonardo J Duarte
- Institute of Chemistry, State University of Campinas, P.O. Box 6154, Postal Code 13.083-970, Campinas, São Paulo, Brazil
| | - Roy E Bruns
- Institute of Chemistry, State University of Campinas, P.O. Box 6154, Postal Code 13.083-970, Campinas, São Paulo, Brazil
| |
Collapse
|
25
|
Auride ion interaction with borane: A theoretical study of AuBH3−. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2021.139288] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
26
|
Gallegos M, Guevara-Vela JM, Pendás ÁM. NNAIMQ: A neural network model for predicting QTAIM charges. J Chem Phys 2022; 156:014112. [PMID: 34998318 DOI: 10.1063/5.0076896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Atomic charges provide crucial information about the electronic structure of a molecular system. Among the different definitions of these descriptors, the one proposed by the Quantum Theory of Atoms in Molecules (QTAIM) is particularly attractive given its invariance against orbital transformations although the computational cost associated with their calculation limits its applicability. Given that Machine Learning (ML) techniques have been shown to accelerate orders of magnitude the computation of a number of quantum mechanical observables, in this work, we take advantage of ML knowledge to develop an intuitive and fast neural network model (NNAIMQ) for the computation of QTAIM charges for C, H, O, and N atoms with high accuracy. Our model has been trained and tested using data from quantum chemical calculations in more than 45 000 molecular environments of the near-equilibrium CHON chemical space. The reliability and performance of NNAIMQ have been analyzed in a variety of scenarios, from equilibrium geometries to molecular dynamics simulations. Altogether, NNAIMQ yields remarkably small prediction errors, well below the 0.03 electron limit in the general case, while accelerating the calculation of QTAIM charges by several orders of magnitude.
Collapse
Affiliation(s)
- Miguel Gallegos
- Depto. Química Física y Analítica, Universidad de Oviedo, 33006 Oviedo, Spain
| | - José Manuel Guevara-Vela
- Institute of Chemistry, National Autonomous University of Mexico, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán, Mexico City C.P. 04510, Mexico
| | - Ángel Martín Pendás
- Depto. Química Física y Analítica, Universidad de Oviedo, 33006 Oviedo, Spain
| |
Collapse
|
27
|
Zheng H, Li Y, Zhang X, Han J, Li S, Wu G, Liu Q, Liu X, Liao W. Interaction-determined extraction capacity between rare earth ions and extractants: taking lanthanum and lutetium as models through theoretical calculations. Inorg Chem Front 2022. [DOI: 10.1039/d2qi01261h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Extractant plays an important role in the separation and purification of rare earth elements (REEs), whereas, extraction performance is the most effective tactic to evaluate whether an extractant is complete...
Collapse
|
28
|
xue D, Chen Z, Liu JY, Wu D, Li ZR, Li Y. High electron affinity triggered by lithium coordination: quasi-chalcogen properties of Li2Sn8Be. Phys Chem Chem Phys 2022; 24:10611-10621. [DOI: 10.1039/d2cp00967f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This work puts forward an unusual but rational strategy to design superatoms mimicking the properties of group VIA elements. A stable dianion with closo-configuration, namely Li2Sn8Be2−, has been obtained by...
Collapse
|
29
|
Khalili F, Vafaee M, Shokri B. Attosecond charge migration following oxygen K-shell ionization in DNA bases and base pairs. Phys Chem Chem Phys 2021; 23:23005-23013. [PMID: 34611693 DOI: 10.1039/d1cp02920g] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Core ionization of DNA begins a cascade of events which could lead to cellular inactivation or death. The created core-hole following an impulse inner-shell ionization of molecules naturally decays in the auger timescale. We simulated charge migration (CM) phenomena following an impulsive core ionization of individual DNA bases at the oxygen K-edge which occurs before Auger decay of the oxygen. Our approach is based on real-time time dependent density functional theory (RT-TDDFT). It is shown that the pronounced hole fluctuation observed around bonds of the initial core-hole results in various valence orbital migrations. Also, the same photo-core-ionized dynamics is studied for the related base pairs. We investigate the role of base pairing and H-bonding interactions in the attosecond CM dynamics. In particular, the creation of a core-hole in the oxygen involved in H-bonding leads to an enhancement of charge migration relative to the respective single bases. Importantly, the hole oscillation of the adenine-thymine base pair upon creation of a core-hole at the oxygen, which does not contribute to the donor-acceptor interactions (not H-bonded), decreases compared to the single thymine base. Understanding the detailed dynamics of the localized core-hole initiating CM process would open the way for chemically controlling DNA damage/repair in the future.
Collapse
Affiliation(s)
- Fatemeh Khalili
- Department of Physics, Shahid Beheshti University, Velenjak, Tehran 19839, Iran.
| | - Mohsen Vafaee
- Department of Chemistry, Tarbiat Modares University, P. O. Box 14115-175, Tehran, Iran.
| | - Babak Shokri
- Department of Physics, Shahid Beheshti University, Velenjak, Tehran 19839, Iran. .,Laser-Plasma Research Institute, Shahid Beheshti University, Velenjak, Tehran 19839, Iran
| |
Collapse
|
30
|
Abstract
Generalized atomic polar tensor (GAPT) has turned into a very popular charge model since it was proposed three decades ago. During this period, several works aiming to compare different partition schemes have included it among their tested models. Nonetheless, GAPT exhibits a set of unique features that prevent it from being directly comparable to "standard" partition schemes. We take this opportunity to explore some of these features, mainly related to the need of evaluating multiple geometries and the dynamic character of GAPT, and show how to obtain the static and dynamic parts of GAPT from any static charge model in the literature. We also present a conceptual evaluation of charge models that aims to explain, at least partially, why GAPT and quantum theory of atoms in molecules (QTAIM) charges are strongly correlated with one another, even though they seem to be constructed under very different frameworks. Similar to GAPT, infrared charges (also derived from atomic polar tensors of planar molecules) are also shown to provide an improved interpretation if they are described as a combination of static charges and changing atomic dipoles rather than just experimental static atomic charges.
Collapse
Affiliation(s)
- Wagner
E. Richter
- Department
of Chemical Engineering, Federal University
of Technology—Paraná, Ponta Grossa 84017-220, PR, Brazil
| | - Leonardo J. Duarte
- Institute
of Chemistry, State University of Campinas, Campinas 13083-861, SP, Brazil
| | - Roy E. Bruns
- Institute
of Chemistry, State University of Campinas, Campinas 13083-861, SP, Brazil
| |
Collapse
|
31
|
Sauza-de la Vega A, Rocha-Rinza T, Guevara-Vela JM. Cooperativity and Anticooperativity in Ion-Water Interactions: Implications for the Aqueous Solvation of Ions. Chemphyschem 2021; 22:1269-1285. [PMID: 33635563 DOI: 10.1002/cphc.202000981] [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: 12/02/2020] [Revised: 02/24/2021] [Indexed: 01/03/2023]
Abstract
Non-additive effects in hydrogen bonds (HB) take place as a consequence of electronic charge transfers. Therefore, it is natural to expect cooperativity and anticooperativity in ion-water interactions. Nevertheless, investigations on this matter are scarce. This paper addresses the interactions of (i) the cations Li+ , Na+ , K+ , Be2+ , Mg2+ , and Ca2+ together with (ii) the anions F- , Cl- , Br- , NO3 - and SO4 2- with water clusters (H2 O)n , n=1-8, and the effects of these ions on the HBs within the complete molecular adducts. We used quantum chemical topology tools, specifically the quantum theory of atoms in molecules and the interacting quantum atoms energy partition to investigate non-additive effects among the interactions studied herein. Our results show a decrease on the interaction energy between ions and the first neighbouring water molecules with an increment of the coordination number. We also found strong cooperative effects in the interplay between HBs and ion-dipole interactions within the studied systems. Such cooperativity affects considerably the interactions among ions with their first and second solvation shells in aqueous environments. Overall, we believe this article provides valuable information about how ion-dipole contacts interact with each other and how they relate to other interactions, such as HBs, in the framework of non-additive effects in aqueous media.
Collapse
Affiliation(s)
- Arturo Sauza-de la Vega
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., 04510, CDMX, México
| | - Tomás Rocha-Rinza
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., 04510, CDMX, México
| | - José Manuel Guevara-Vela
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., 04510, CDMX, México
| |
Collapse
|
32
|
Explicit solvation effects on low-index Fe surfaces and small particles as adsorbents of Arsenic species: a DFT study. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02767-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
|
33
|
Janeček M, Kührová P, Mlýnský V, Otyepka M, Šponer J, Banáš P. W-RESP: Well-Restrained Electrostatic Potential-Derived Charges. Revisiting the Charge Derivation Model. J Chem Theory Comput 2021; 17:3495-3509. [PMID: 33999623 DOI: 10.1021/acs.jctc.0c00976] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Representation of electrostatic interactions by a Coulombic pairwise potential between atom-centered partial charges is a fundamental and crucial part of empirical force fields used in classical molecular dynamics simulations. The broad success of the AMBER force-field family originates mainly from the restrained electrostatic potential (RESP) charge model, which derives partial charges to reproduce the electrostatic field around the molecules. However, the description of the electrostatic potential around molecules by standard RESP may be biased for some types of molecules. In this study, we modified the RESP charge derivation model to improve its description of the electrostatic potential around molecules and thus electrostatic interactions in the force field. In particular, we reoptimized the atomic radii for definition of the grid points around the molecule, redesigned the restraining scheme, and included extra point (EP) charges. The RESP fitting was significantly improved for aromatic heterocyclic molecules. Thus, the suggested W-RESP(-EP) charge derivation model shows some potential for improving the performance of the nucleic acid force fields, for which the poor description of nonbonded interactions, such as the underestimated stability of base pairing, is well-established. We also report some preliminary simulation tests (around 1 ms of simulation data) on A-RNA duplexes, tetranucleotides, and tetraloops. The simulations reveal no adverse effects, while the description of base-pairing interactions might be improved. The new charges can thus be used in future attempts to improve the nucleic acid simulation force fields, in combination with reparametrization of the other terms.
Collapse
Affiliation(s)
- Michal Janeček
- Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic
| | - Vojtěch Mlýnský
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University, Křížkovského 511/8, 779 00 Olomouc, Czech Republic.,Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
34
|
Huang J, Cai X, Yin H, Li Y, Lin W, Huang S, Zhang Y. A New Candidate in Polyanionic Compounds for a Potassium-Ion Battery Cathode: KTiOPO 4. J Phys Chem Lett 2021; 12:2721-2726. [PMID: 33705136 DOI: 10.1021/acs.jpclett.1c00286] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
First-principles computations were performed to investigate the performance of KTiOPO4 (KTP) as a cathode material for potassium-ion batteries (PIBs), including the stability and electronic properties of depotassiated structures and mechanisms of K deintercalation and diffusion. As depotassiation proceeds, oxygen hole polarons are produced, and there are not peroxides or superoxides formed after deep depotassiation. The anionic oxygen redox in KTP provides a voltage vs K/K+ over 4 V by the PBE+U method and over 5 V with the more reliable HSE06 hybrid functional. When all K in KTP is removed, the calculated volume compression is only 1.528%. The AIMD simulations at 300 K for TiOPO4 verify its thermal stability. The PBE+U calculations predict a low ion diffusion barrier of 0.29 eV in bulk KTP, indicating a good charge-discharge rate for KTP as a cathode for PIBs. All of the calculated results indicate that KTP can be a promising cathode material for PIBs.
Collapse
Affiliation(s)
- Jiajia Huang
- College of Chemistry, Fuzhou University, Fuzhou, Fujian 350108, P. R. China
| | - Xu Cai
- College of Chemistry, Fuzhou University, Fuzhou, Fujian 350108, P. R. China
| | - Huimin Yin
- College of Chemistry, Fuzhou University, Fuzhou, Fujian 350108, P. R. China
| | - Yi Li
- College of Chemistry, Fuzhou University, Fuzhou, Fujian 350108, P. R. China
| | - Wei Lin
- College of Chemistry, Fuzhou University, Fuzhou, Fujian 350108, P. R. China
| | - Shuping Huang
- College of Chemistry, Fuzhou University, Fuzhou, Fujian 350108, P. R. China
- Fujian Provincial Key Laboratory of Electrochemical Energy Storage Materials, Fuzhou University, Fuzhou, Fujian 350108, P. R. China
| | - Yongfan Zhang
- College of Chemistry, Fuzhou University, Fuzhou, Fujian 350108, P. R. China
| |
Collapse
|
35
|
Kancharlapalli S, Gopalan A, Haranczyk M, Snurr RQ. Fast and Accurate Machine Learning Strategy for Calculating Partial Atomic Charges in Metal-Organic Frameworks. J Chem Theory Comput 2021; 17:3052-3064. [PMID: 33739834 DOI: 10.1021/acs.jctc.0c01229] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Computational high-throughput screening using molecular simulations is a powerful tool for identifying top-performing metal-organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges are often required to model the electrostatic interactions between the MOF and the adsorbate, especially when the adsorption involves molecules with dipole or quadrupole moments such as water and CO2. Although ab initio methods can be used to calculate accurate partial atomic charges, these methods are impractical for screening large material databases because of the high computational cost. We developed a random forest machine learning model to predict the partial atomic charges in MOFs using a small yet meaningful set of features that represent both the elemental properties and the local environment of each atom. The model was trained and tested on a collection of about 320 000 density-derived electrostatic and chemical (DDEC) atomic charges calculated on a subset of the Computation-Ready Experimental Metal-Organic Framework (CoRE MOF-2019) database and separately on charge model 5 (CM5) charges. The model predicts accurate atomic charges for MOFs at a fraction of the computational cost of periodic density functional theory (DFT) and is found to be transferable to other porous molecular crystals and zeolites. A strong correlation is observed between the partial atomic charge and the average electronegativity difference between the central atom and its bonded neighbors.
Collapse
Affiliation(s)
- Srinivasu Kancharlapalli
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, United States.,Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai-400085, India
| | - Arun Gopalan
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, United States
| | - Maciej Haranczyk
- IMDEA Materials Institute, C/Eric Kandel 2, 28906 Getafe, Madrid, Spain
| | - Randall Q Snurr
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, United States
| |
Collapse
|
36
|
Xue D, Wu D, Chen Z, Li Y, Sun W, Liu J, Li Z. On Close Parallels between the Zintl-Based Superatom Ge9Be and Chalcogen Elements. Inorg Chem 2021; 60:3196-3206. [DOI: 10.1021/acs.inorgchem.0c03531] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Duomei Xue
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Di Wu
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Zeren Chen
- School of Mechanical and Aerospace Engineering, Jilin University, Changchun 130025, P. R. China
| | - Ying Li
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China
| | - WeiMing Sun
- Department of Basic Chemistry, The School of Pharmacy, Fujian Medical University, Fuzhou 350108, P. R. China
| | - Jingyao Liu
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Zhiru Li
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China
| |
Collapse
|
37
|
Lebedev AV. Pyridine and Methylpyridines: Calculations of the Structure, Proton Affinity, Gas-Phase Basicity, and Mobility of Protonated Molecules and Proton-Bound Dimers. JOURNAL OF ANALYTICAL CHEMISTRY 2020. [DOI: 10.1134/s1061934820130079] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
38
|
Pei Z, Yang J, Deng J, Mao Y, Wu Q, Yang Z, Wang B, Aikens CM, Liang W, Shao Y. Analysis and visualization of energy densities. II. Insights from linear-response time-dependent density functional theory calculations. Phys Chem Chem Phys 2020; 22:26852-26864. [PMID: 33216085 PMCID: PMC8258743 DOI: 10.1039/d0cp04207b] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Inspired by the analysis of Kohn-Sham energy densities by Nakai and coworkers, we extended the energy density analysis to linear-response time-dependent density functional theory (LR-TDDFT) calculations. Using ethylene-tetrafluoroethylene and oxyluciferin-water complexes as examples, distinctive distribution patterns were demonstrated for the excitation energy densities of local excitations (within a molecular fragment) and charge-transfer excitations (between molecular fragments). It also provided a simple way to compute the effective energy of both hot carriers (particle and hole) from charge-transfer excitations via an integration of the excitation energy density over the donor and acceptor grid points.
Collapse
Affiliation(s)
- Zheng Pei
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, P. R. China.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
39
|
Yu D, Wu D, Liu JY, Li Y, Sun WM. Unveiling the potential of superalkali cation Li 3+ for capturing nitrogen. Phys Chem Chem Phys 2020; 22:26536-26543. [PMID: 33188670 DOI: 10.1039/d0cp03769a] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The potential of the superalkali cation Li3+ for capturing N2 and its behavior in gaseous nitrogen have been theoretically studied at the MP2/6-311+G(d) level. The evolution of structures and stability of the Li3+(N2)n (n = 1-7) complexes shows that the N2 molecules tend to bind to different vertices of the Li3+ core, and that Li3+ might have the capacity to capture up to twelve nitrogen molecules in the first coordination shell. Based on natural population and molecular orbital analyses, Li3+ keeps its superatom identity in the lowest-lying Li3+(N2)n (n = 1-4) complexes. The change in the Gibbs free energies of possible fragmentation channels also indicates the thermodynamic stability of Li3+ in the (N2)n clusters when n ≤ 4. Different from the case of Li3+(H2O)n, where the electrostatic interaction is dominant, the electrostatic and polarization components are found to make nearly equal contributions to Li3+(N2)n complex formation. In addition, it can be concluded that the superalkali cation Li3+ surpasses heavy alkali metal cations in capturing N2 molecules, since it has a larger binding energy with N2 than Na+ and K+ ions.
Collapse
Affiliation(s)
- Dan Yu
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China.
| | | | | | | | | |
Collapse
|
40
|
Yu D, Wu D, Liu JY, Li SY, Li Y. On single-electron magnesium bonding formation and the effect of methyl substitution. RSC Adv 2020; 10:34413-34420. [PMID: 35514394 PMCID: PMC9056782 DOI: 10.1039/d0ra06591a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Accepted: 09/08/2020] [Indexed: 12/14/2022] Open
Abstract
The complexes formed between MgX2 (X = F, H) molecules and alkyl radicals Y [Y = CH3, CH2CH3, CH(CH3)2, and C(CH3)3] have been characterized by using quantum chemical methods. The binding distance in all cases is less than the sum of vdW radii of Mg and C, indicating the formation of a non-covalent interaction, namely single-electron magnesium bond. Energy decomposition analysis reveals that electrostatic and polarization contributions are the major components responsible for the stability of the studied complexes. According to interaction energy, atoms in molecules, and independent gradient model analyses, methyl substitution on electron donor Y imposes a positive effect on its complexation with MgX2. When compared with other nonbonded interactions, the single-electron magnesium bond is found to have strength comparable to those of the single-electron beryllium bond and π-magnesium bond. The complexes formed between MgX2 (X = F, H) molecules and alkyl radicals Y [Y = CH3, CH2CH3, CH(CH3)2, and C(CH3)3] have been characterized by using quantum chemical methods.![]()
Collapse
Affiliation(s)
- Dan Yu
- Institute of Theoretical Chemistry, Laboratory of Theoretical and Computational Chemistry, College of Chemistry, Jilin University Changchun 130023 P. R. China
| | - Di Wu
- Institute of Theoretical Chemistry, Laboratory of Theoretical and Computational Chemistry, College of Chemistry, Jilin University Changchun 130023 P. R. China
| | - Jing-Yao Liu
- Institute of Theoretical Chemistry, Laboratory of Theoretical and Computational Chemistry, College of Chemistry, Jilin University Changchun 130023 P. R. China
| | - Si-Yi Li
- Department of Transdisciplinary Science and Engineering, Tokyo Institute of Technology 2-12-1 Ookayama, Meguro Tokyo 152-8551 Japan
| | - Ying Li
- Institute of Theoretical Chemistry, Laboratory of Theoretical and Computational Chemistry, College of Chemistry, Jilin University Changchun 130023 P. R. China
| |
Collapse
|
41
|
Młodzikowska-Pieńko K, Trzaskowski B. Rate-Limiting Steps in the Intramolecular C-H Activation of Ruthenium N-Heterocyclic Carbene Complexes. J Phys Chem A 2020; 124:3609-3617. [PMID: 32293892 PMCID: PMC7590980 DOI: 10.1021/acs.jpca.0c01354] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ruthenium (II) complexes with N-heterocyclic carbenes (NHC) are commonly used as efficient catalysts in hydrogenation of olefins with simultaneous intramolecular C-H activation. Using the DFT approach, we have investigated the entire hydrogenation reaction pathway for four new potential catalysts and ethylene, a model substrate. Our calculations imply that the dissociation of phosphine is the rate-limiting step of hydrogenation, contrary to recent computational results. We also found that catalysts bearing NHCs with aliphatic and aromatic side groups are energetically favorable over other aliphatic cyclohexyl-substituted NHC. To examine how electronic properties of various catalysts influence the energetic barrier in the crucial steps of the reaction, we applied the Noncovalent Interaction analysis, which allowed us to reveal crucial interactions which stabilize/destabilize important intermediates and transition states in the hydrogenation reaction.
Collapse
Affiliation(s)
- Katarzyna Młodzikowska-Pieńko
- Centre of New Technologies, University of Warsaw, Banacha 2C Street, 02-097 Warsaw, Poland.,Faculty of Chemistry, University of Warsaw, Pasteura 1 Street, 02-093 Warsaw, Poland
| | - Bartosz Trzaskowski
- Centre of New Technologies, University of Warsaw, Banacha 2C Street, 02-097 Warsaw, Poland
| |
Collapse
|
42
|
K. Macedo G, Haiduke RLA. A Quantum Theory Atoms in Molecules Study about the Inductive Effect of Substituents in Methane Derivatives. ACS OMEGA 2020; 5:9041-9045. [PMID: 32337469 PMCID: PMC7178782 DOI: 10.1021/acsomega.0c01081] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/10/2020] [Accepted: 03/30/2020] [Indexed: 06/11/2023]
Abstract
The substituent effect on the covalent character of C-H bonds in methane derivatives is evaluated by means of local descriptors based on the topology of the electron density. Halogens, -OH, -SH, =O, =S, -NO2, -NH2, and -OCH3 increase the covalent character of the remaining C-H bonds, while alkaline metals (-Li and -Na) result in the opposite trend. This study proposes that the inductive effect is due to polarization changes driven by substituent charges.
Collapse
|
43
|
Cho M, Sylvetsky N, Eshafi S, Santra G, Efremenko I, Martin JML. The Atomic Partial Charges Arboretum: Trying to See the Forest for the Trees. Chemphyschem 2020; 21:688-696. [PMID: 32052532 PMCID: PMC7317385 DOI: 10.1002/cphc.202000040] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Revised: 02/12/2020] [Indexed: 11/20/2022]
Abstract
Atomic partial charges are among the most commonly used interpretive tools in quantum chemistry. Dozens of different 'population analyses' are in use, which are best seen as proxies (indirect gauges) rather than measurements of a 'general ionicity'. For the GMTKN55 benchmark of nearly 2,500 main-group molecules, which span a broad swathe of chemical space, some two dozen different charge distributions were evaluated at the PBE0 level near the 1-particle basis set limit. The correlation matrix between the different charge distributions exhibits a block structure; blocking is, broadly speaking, by charge distribution class. A principal component analysis on the entire dataset suggests that nearly all variation can be accounted for by just two 'principal components of ionicity': one has all the distributions going in sync, while the second corresponds mainly to Bader QTAIM vs. all others. A weaker third component corresponds to electrostatic charge models in opposition to the orbital-based ones. The single charge distributions that have the greatest statistical similarity to the first principal component are iterated Hirshfeld (Hirshfeld-I) and a minimal-basis projected modification of Bickelhaupt charges. If three individual variables, rather than three principal components, are to be identified that contain most of the information in the whole dataset, one representative for each of the three classes of Corminboeuf et al. is needed: one based on partitioning of the density (such as QTAIM), a second based on orbital partitioning (such as NPA), and a third based on the molecular electrostatic potential (such as HLY or CHELPG).
Collapse
Affiliation(s)
- Minsik Cho
- Department of Organic ChemistryWeizmann Institute of Science76100RehovotIsrael
- Present Address: Department of ChemistryBrown UniversityProvidenceRhode Island02912USA
| | - Nitai Sylvetsky
- Department of Organic ChemistryWeizmann Institute of Science76100RehovotIsrael
| | - Sarah Eshafi
- Department of Organic ChemistryWeizmann Institute of Science76100RehovotIsrael
- Present Address: Integrated Science ProgramMcMaster UniversityHamiltonOntarioL8S 4 M1Canada
| | - Golokesh Santra
- Department of Organic ChemistryWeizmann Institute of Science76100RehovotIsrael
| | - Irena Efremenko
- Department of Organic ChemistryWeizmann Institute of Science76100RehovotIsrael
| | - Jan M. L. Martin
- Department of Organic ChemistryWeizmann Institute of Science76100RehovotIsrael
| |
Collapse
|
44
|
Ramírez-Palma LG, García-Jacas CR, García-Ramos JC, Almada-Monter R, Galindo-Murillo R, Cortés-Guzmán F. Pharmacophoric sites of anticancer metal complexes located using quantum topological atomic descriptors. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2019.127480] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
|
45
|
Zheng HF, Xu J, Ding YH. A sixteen-valence-electron carbon-group 13 family with global penta-atomic planar tetracoordinate carbon: an ionic strategy. Phys Chem Chem Phys 2020; 22:3975-3982. [PMID: 32022042 DOI: 10.1039/c9cp06577f] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The design of planar tetracoordinate carbon (ptC) has always been a challenge due to its unique bonding mode that necessitates the perfect balance between the carbon center and surrounding ligands both electronically and mechanically. A unique type of 18-valence-electron (18ve) template, i.e., CAl42-, has been found to be very effective in designing various novel 18ve-species upon skeletal substitution. In this work, we showed that though ptC is not the global structure for the parent 16ve-CAl4, suitable skeletal substitution can allow for a series of global minimum ptC species. Theoretical calculations at the level of CCSD(T)/def2-QZVP//B3LYP/def2-QZVP for 35 carbon-group 13 systems with 16-ve, i.e., CXaYbZcKd (X, Y, Z, K = Al/Ga/In/Tl; 0 ≤ a, b, c, d ≤ 4, a + b + c + d = 4), showed that 9 systems (CAl3Tl, CGa3Tl, CGa2Tl2, CAl2GaTl, CAl2InTl, CGa2InTl, CAlGa2Tl, CGa2InTl and CAlGaInTl) possess global minimum ptC and 2 systems (CAl3In and CAl2Tl2) have quasi-GM ptC. Except for CAl3Tl and CAl3In, all the ptCs were predicted for the first time. All these stable ptC structures have the same skeleton and can be described as the same ionic sub-structure, i.e., [A-]B+. This study not only enriches 16ve-ptC, but also directly demonstrates that utilizing an ionic strategy, non-ptC CAl4 also can be used as a template to extend the ptC family.
Collapse
Affiliation(s)
- Hai-Feng Zheng
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, P. R. China.
| | - Jing Xu
- Department of Optical Engineering Zhejiang A&F University, Lin'an, Zhejiang 311300, P. R. China.
| | - Yi-Hong Ding
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, P. R. China. and Key Laboratory of Carbon Materials of Zhejiang Province, College of Chemistry and Materials Engineering, Wenzhou University, Wenzhou 325035, People's Republic of China.
| |
Collapse
|
46
|
Zheng HF, Xu J, Ding YH. Mono-silicon isoelectronic replacement in CAl 4 : van't hoff/le bel carbon or not? J Comput Chem 2020; 41:119-128. [PMID: 31663141 DOI: 10.1002/jcc.26079] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2019] [Revised: 09/09/2019] [Accepted: 09/09/2019] [Indexed: 11/09/2022]
Abstract
In cluster studies, the isoelectronic replacement strategy has been successfully used to introduce new elements into a known structure while maintaining the desired topology. The well-known penta-atomic 18 valence electron (ve) species C Al 4 2 - and its Al- /Si or Al/Si+ isoelectronically replaced clusters CAl3 Si- , CAl2 Si2 , C AlSi 3 - , and C Si 4 2 + , all possess the same anti-van't Hoff/Le Bel skeletons, that is, nontraditional planar tetracoordinate carbon (ptC) structure. In this article, however, we found that such isoelectronic replacement between Si and Al does not work for the 16ve-CAl4 with the traditional van't Hoff/Le Bel tetrahedral carbon (thC) and its isoelectronic derivatives CAl3 X (X = Ga/In/Tl). At the level of CCSD(T)/def2-QZVP//B3LYP/def2-QZVP, none of the global minima of the 16ve mono-Si-containing clusters CAl2 SiX+ (X = Al/Ga/In/Tl) maintains thC as the parent CAl4 does. Instead, X = Al/Ga globally favors an unusual ptC structure that has one long C─X distance yet with significant bond index value, and X = In/Tl prefers the planar tricoordinate carbon. The frustrated formation of thC in these clusters is ascribed to the CSi bonding that prefers a planar fashion. Inclusion of chloride ion would further stabilize the ptC of CAl2 SiAl+ and CAl2 SiGa+ . The unexpectedly disclosed CAl2 SiAl+ and CAl2 SiGa+ represent the first type of 16ve-cationic ptCs with multiple bonds. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Hai-Feng Zheng
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun, 130023, People's Republic of China
| | - Jing Xu
- Department of Optical Engineering, Zhejiang A&F University, Lin'an, Zhejiang, 311300, People's Republic of China
| | - Yi-Hong Ding
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun, 130023, People's Republic of China
| |
Collapse
|
47
|
Ramírez-Palma DI, García-Jacas CR, Carpio-Martínez P, Cortés-Guzmán F. Predicting reactive sites with quantum chemical topology: carbonyl additions in multicomponent reactions. Phys Chem Chem Phys 2020; 22:9283-9289. [DOI: 10.1039/d0cp00300j] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The reactivity of an atom within a molecule depends mostly on the way the electron density polarizes reflected in the quadrupole moment of the reactive atom.
Collapse
Affiliation(s)
| | - Cesar R. García-Jacas
- Cátedras CONACYT – Departamento de Ciencias de la Computación
- Centro de Investigación Científica y de Educación Superior de Ensenada (CICESE)
- Ensenada
- Mexico
| | | | | |
Collapse
|
48
|
Huang PW, Wang CZ, Wu QY, Lan JH, Chai ZF, Shi WQ. Quantum chemical studies of selective back-extraction of Am(III) from Eu(III) and Cm(III) with two hydrophilic 1,10-phenanthroline-2,9-bis-triazolyl ligands. RADIOCHIM ACTA 2019. [DOI: 10.1515/ract-2019-3197] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Abstract
We theoretically investigated the selective back-extraction towards Am(III) over Eu(III) and Cm(III) with two water-soluble 2,9-bis-triazolyl-1,10-phenanthroline derivatives BTrzPhen1 (with two ethanol side chains) and BTrz-Phen2 (with two 1,2-butanediol side chains) by density functional theory (DFT). The molecular geometries and formation reactions of the metal-ligand complexes were modeled by using M(BTrzPhen)(NO3)3 and [M(BTrzPhen)2(NO3)]2+. Am(III) selectivity over Eu(III) and Cm(III) with BTrzPhen2 was successfully reproduced by back-extraction reaction free energy analysis. Moreover, bonding properties between the metal cations and coordinated ligands (model complexes) were studied in terms of Mayer bond order and quantum theory of atoms in molecule (QTAIM). The difference in covalency between An–N and Eu–N bonds were found to be the key factors for Am(III)/Eu(III) separation, while the Am(III) selectivity over Cm(III) of BTrzPhen2 might be attributed to the competition of donor atoms for cation binding preference toward Am(III) and Cm(III).
Collapse
Affiliation(s)
- Pin-Wen Huang
- Zhejiang University of Water Resources and Electric Power , Hangzhou, Zhejiang 310018 , China
| | - Cong-Zhi Wang
- Laboratory of Nuclear Energy Chemistry , Institute of High Energy Physics, Chinese Academy of Sciences , Beijing 100049 , China
| | - Qun-Yan Wu
- Laboratory of Nuclear Energy Chemistry , Institute of High Energy Physics, Chinese Academy of Sciences , Beijing 100049 , China
| | - Jian-Hui Lan
- Laboratory of Nuclear Energy Chemistry , Institute of High Energy Physics, Chinese Academy of Sciences , Beijing 100049 , China
| | - Zhi-Fang Chai
- Laboratory of Nuclear Energy Chemistry , Institute of High Energy Physics, Chinese Academy of Sciences , Beijing 100049 , China
- Engineering Laboratory of Advanced Energy Materials, Ningbo Institute of Industrial Technology , Chinese Academy of Sciences , Ningbo, Zhejiang 315201 , China
| | - Wei-Qun Shi
- Laboratory of Nuclear Energy Chemistry , Institute of High Energy Physics, Chinese Academy of Sciences , Beijing 100049 , China
| |
Collapse
|
49
|
Bakó I, Daru J, Pothoczki S, Pusztai L, Hermansson K. Effects of H-bond asymmetry on the electronic properties of liquid water – An AIMD analysis. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111579] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
50
|
Chang XT, Li Y, Liu JY, Ma HD, Wu D. Noble gas insertion compounds of hydrogenated and lithiated hyperhalogens. Phys Chem Chem Phys 2019; 21:20156-20165. [PMID: 31483426 DOI: 10.1039/c9cp01284b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Based on density functional theory (DFT) calculations, hydrogenated hyperhalogen HM(BO2)2, lithiated hyperhalogen LiM(BO2)2 (M = Cu, Ag, Au), and their compounds with xenon were studied. Different insertion sites of Xe resulted in various isomers. According to the natural population analysis, the Xe atom donated 0.12-0.77 electrons to HM(BO2)2 and 0.14-0.41 electrons to LiM(BO2)2 when they combined, leading to metastable charge-transfer compounds in most cases. The nature of bonding between xenon and HM(BO2)2/LiM(BO2)2 was found to be related to its location. Covalent bonds were formed when Xe bound with hydrogen atoms, as indicated by the large Wiberg bond indices of the Xe-H bonds. The same was true for most Xe-M bondings. When an Xe-O connection was formed, it was either an ionic or van der Waals force in nature depending on the specific structural feature of the isomer. A parallel study on hyperhalogen-supported Ar and Kr compounds indicated that they were not very stable and were less likely to exist at room temperature, which was in accordance with the high inertness of both Ar and Kr atoms.
Collapse
Affiliation(s)
- Xiao-Ting Chang
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, P. R. China.
| | | | | | | | | |
Collapse
|