1
|
Domenichini G. Extending the definition of atomic basis sets to atoms with fractional nuclear charge. J Chem Phys 2024; 160:124107. [PMID: 38526100 DOI: 10.1063/5.0196383] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2024] [Accepted: 03/10/2024] [Indexed: 03/26/2024] Open
Abstract
Alchemical transformations showed that perturbation theory can be applied also to changes in the atomic nuclear charges of a molecule. The alchemical path that connects two different chemical species involves the conceptualization of a non-physical system in which an atom possess a non-integer nuclear charge. A correct quantum mechanical treatment of these systems is limited by the fact that finite size atomic basis sets do not define exponents and contraction coefficients for fractional charge atoms. This paper proposes a solution to this problem and shows that a smooth interpolation of the atomic orbital coefficients and exponents across the periodic table is a convenient way to produce accurate alchemical predictions, even using small size basis sets.
Collapse
Affiliation(s)
- Giorgio Domenichini
- Faculty of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| |
Collapse
|
2
|
Kaya S, Robles-Navarro A, Mejía E, Gómez T, Cardenas C. On the Prediction of Lattice Energy with the Fukui Potential: Some Supports on Hardness Maximization in Inorganic Solids. J Phys Chem A 2022; 126:4507-4516. [PMID: 35766899 PMCID: PMC9289887 DOI: 10.1021/acs.jpca.1c09898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Using perturbation theory within the framework of conceptual density functional theory, we derive a lower bound for the lattice energy of the ionic solids. The main element of the lower bound is the Fukui potential in the nuclei of the molecule corresponding to the unit formula of the solid. Thus, we propose a model to calculate the lattice energy in terms of the Fukui potential. Our method, which is extremely simple, performs well as other methods using the crystal structure information of alkali halide solids. The method proposed here correlates surprisingly well with the experimental data on the lattice energy of a diverse series of solids having even a non-negligible covalent characteristic. Finally, the validity of the maximum hardness principle (MHP) is assessed, showing that in this case, the MHP is limited.
Collapse
Affiliation(s)
- Savaş Kaya
- Health Services Vocational School, Department of Pharmacy, Sivas Cumhuriyet University, Sivas58140, Turkey
| | - Andrés Robles-Navarro
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Las Palmeras 3425, Santiago Casilla653, Chile
| | - Erica Mejía
- Facultad de Ingeniería-(Medellin-Colombia), Institución Universitaria Pascual Bravo, Medellín050025, Colombia
| | - Tatiana Gómez
- Theoretical and Computational Chemistry Center, Institute of Applied Chemical Sciences, Faculty of Engineering, Universidad Autonoma de Chile, Santiago9170124, Chile
| | - Carlos Cardenas
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Las Palmeras 3425, Santiago Casilla653, Chile.,Centro para el Desarrollo de la Nanociencia y la Nanotecnología (CEDENNA), Avda. Ecuador 3493, Santiago9170124, Chile
| |
Collapse
|
3
|
Miranda-Quintana RA, Heidar-Zadeh F, Fias S, Chapman AEA, Liu S, Morell C, Gómez T, Cárdenas C, Ayers PW. Molecular Interactions From the Density Functional Theory for Chemical Reactivity: The Interaction Energy Between Two-Reagents. Front Chem 2022; 10:906674. [PMID: 35769444 PMCID: PMC9234655 DOI: 10.3389/fchem.2022.906674] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 04/19/2022] [Indexed: 12/13/2022] Open
Abstract
Reactivity descriptors indicate where a reagent is most reactive and how it is most likely to react. However, a reaction will only occur when the reagent encounters a suitable reaction partner. Determining whether a pair of reagents is well-matched requires developing reactivity rules that depend on both reagents. This can be achieved using the expression for the minimum-interaction-energy obtained from the density functional reactivity theory. Different terms in this expression will be dominant in different circumstances; depending on which terms control the reactivity, different reactivity indicators will be preferred.
Collapse
Affiliation(s)
- Ramón Alain Miranda-Quintana
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, FL, United States
- *Correspondence: Ramón Alain Miranda-Quintana, ; Carlos Cárdenas, ; Paul W. Ayers, ; Tatiana Gómez,
| | | | - Stijn Fias
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON, Canada
| | - Allison E. A. Chapman
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC, United States
| | - Christophe Morell
- Université de Lyon, Université Claude Bernard Lyon 1, Institut des Sciences Analytiques-UMR CNRS 5280, Villeurbanne, France
| | - Tatiana Gómez
- Theoretical and Computational Chemistry Center, Institute of Applied Chemical Sciences, Faculty of Engineering, Universidad Autonoma de Chile, Santiago, Chile
- *Correspondence: Ramón Alain Miranda-Quintana, ; Carlos Cárdenas, ; Paul W. Ayers, ; Tatiana Gómez,
| | - Carlos Cárdenas
- Departamento de Fisica, Facultad de Ciencias, Universidad de Chile, Santiago, Chile
- Centro para el desarrollo de la Nanociencias y Nanotecnologia, CEDENNA, Santiago, Chile
- *Correspondence: Ramón Alain Miranda-Quintana, ; Carlos Cárdenas, ; Paul W. Ayers, ; Tatiana Gómez,
| | - Paul W. Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON, Canada
- *Correspondence: Ramón Alain Miranda-Quintana, ; Carlos Cárdenas, ; Paul W. Ayers, ; Tatiana Gómez,
| |
Collapse
|
4
|
Abstract
We propose to relax geometries throughout chemical compound space (CCS) using alchemical perturbation density functional theory (APDFT). APDFT refers to perturbation theory involving changes in nuclear charges within approximate solutions to Schr\"odinger's equation. We give an analytical formula to calculate the mixed second order energy derivatives with respect to both, nuclear charges and nuclear positions (named "alchemical force"), within the restricted Hartree-Fock case.We have implemented and studied the formula for its use in geometry relaxation of various reference and target molecules.We have also analysed the convergence of the alchemical force perturbation series, as well as basis set effects.Interpolating alchemically predicted energies, forces, and Hessian to a Morse potential yields more accurate geometries and equilibrium energies than when performing a standard Newton Raphson step. Our numerical predictions for small molecules including BF, CO, N2, CH$_4$, NH$_3$, H$_2$O, and HF yield mean absolute errors of of equilibrium energies and bond lengths smaller than 10 mHa and 0.01 Bohr for 4$^\text{th}$ order APDFT predictions, respectively.Our alchemical geometry relaxation still preserves the combinatorial efficiency of APDFT: Based on a single coupled perturbed Hartree Fock derivative for benzene we provide numerical predictions of equilibrium energies and relaxed structures of all the 17 iso-electronic charge-netural BN-doped mutants with averaged absolute deviations of $\sim$27 mHa and $\sim$0.12 Bohr, respectively.
Collapse
|
5
|
Protonation of Borylated Carboxonium Derivative [2,6-B10H8O2CCH3]−: Theoretical and Experimental Investigation. Int J Mol Sci 2022; 23:ijms23084190. [PMID: 35457007 PMCID: PMC9025682 DOI: 10.3390/ijms23084190] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2022] [Revised: 04/07/2022] [Accepted: 04/08/2022] [Indexed: 01/03/2023] Open
Abstract
The process of protonation of [2,6-B10H8O2CCH3]− was investigated both theoretically and experimentally. The most suitable conditions for protonation of the derivative [2,6-B10H8O2CCH3]− were found. The process of protonation was carried out in the presence of an excess of trifluoromethanesulfonic acid CF3SO3H at room temperature in dichloromethane solution. The structure of the resulting complex [2,6-B10H8O2CCH3*Hfac]0 was established using NMR data and the results of DFT calculations. An additional proton atom Hfac was found to be localized on one of the facets that was opposite the boron atom in a substituted position, and which bonded mainly with one apical boron atom. The main descriptors of the B-Hfac bond were established theoretically using QTAIM and NBO approaches. In addition, the mechanism of [2,6-B10H8O2CCH3]− protonation was investigated.
Collapse
|
6
|
Klyukin IN, Vlasova YS, Novikov AS, Zhdanov AP, Hagemann HR, Zhizhin KY, Kuznetsov NT. B-F bonding and reactivity analysis of mono- and perfluoro-substituted derivatives of closo-borate anions (6, 10, 12): A computational study. Polyhedron 2022. [DOI: 10.1016/j.poly.2021.115559] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|