Benchmark Density Functional Theory Approach for the Calculation of Bond Dissociation Energies of the M–O2 Bond: A Key Step in Water Splitting Reactions
Ontology highlight
ABSTRACT:
SUBMITTER: Kosar N
PROVIDER: S-EPMC9348009 | biostudies-literature | 2022 Jun
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Transition metal silicides are promising materials for improved electronic devices, and this motivates achieving a better understanding of transition metal bonds to silicon. Here we model the ground and excited state bond dissociations of VSi, NbSi, and TaSi using a complete active space (CAS) wave function and a separated-pair (SP) wave function combined with two post-self-consistent field techniques: complete active space with perturbation theory at second order and multiconfiguration pair-density functional theory. The SP approximation is a multiconfiguration self-consistent field method with a selection of configurations based on generalized valence bond theory without the perfect pairing approximation. For both CAS and SP, the active-space composition corresponds to the nominal correlated-participating-orbital scheme. The ground state and low-lying excited states are explored to predict the state ordering for each molecule, and potential energy curves are calculated for the ground state to compare to experiment. The experimental bond dissociation energies of the three diatomic molecules are predicted with eight on-top pair-density functionals with a typical error of 0.2 eV for a CAS wave function and a typical error of 0.3 eV for the SP approximation. We also provide a survey of the accuracy achieved by the SP and extended separated-pair approximations for a broader set of 25 transition metal-ligand bond dissociation energies.
| S-EPMC8152470 | biostudies-literature
Project description:Thermodynamic properties, i.e., bond dissociation energies and enthalpy of formation, of chlorinated and brominated polycyclic aromatic hydrocarbons play a fundamental role in understanding their formation mechanisms and reactivity. Computational electronic structure calculations routinely used to predict thermodynamic properties of various species are limited for these compounds due to large computational cost to obtain accurate results by employing high-level wave function theory methods. In this work, a number of composite model chemistry methods (CBS-QB3, G3MP2, G3, and G4) are used to compute bond dissociation energies and enthalpies of formation of small to medium-size chlorinated and brominated polycyclic aromatic hydrocarbon compounds. The enthalpy of formation is derived via the atomization method and compared against the recommended values. Statistical analysis indicates that G4 is the best method. For comparison, three commonly used density functional theory (DFT) methods (M06-2X, ωB97X-D and B2PLYP-D3) with various basis sets including 6-311++G(d, p), cc-pVTZ, and cc-pVQZ in the prediction of bond dissociation energies and enthalpies of formation have been tested using the optimized geometries at the same M06-2X/6-311++G(d, p) level of theory. It is found that ωB97X-D/6-311++G(d, p) shows the best performance in computing the bond dissociation energies, while ωB97X-D/cc-pVTZ exhibits the best prediction in enthalpy of formation of the studied reaction systems. The structural effect on the bond dissociation energies and enthalpy of formation of chlorinated and brominated polycyclic aromatic hydrocarbons are then systematically analyzed. Based on comparisons of the various methods, reliable DFT methods are recommended for future theoretical studies on large chlorinated and brominated polycyclic aromatic hydrocarbons considering both accuracy and computational cost. This work, to the authors' knowledge, is the first to systematically benchmark theoretical methods for the accurate prediction of thermodynamic properties for chlorinated and brominated polycyclic aromatic hydrocarbons.
| S-EPMC9040899 | biostudies-literature
Project description:Knowledge of reliable X-H bond dissociation energies (X = C, N, O, S) for amino acids in proteins is key for studying the radical chemistry of proteins. X-H bond dissociation energies of model dipeptides were computed using the isodesmic reaction method at the BMK/6-31+G(2df,p) and G4(MP2)-6X levels of theory. The density functional theory values agree well with the composite-level calculations. By this high level of theory, combined with a careful choice of reference compounds and peptide model systems, our work provides a highly valuable data set of bond dissociation energies with unprecedented accuracy and comprehensiveness. It will likely prove useful to predict protein biochemistry involving radicals, e.g., by machine learning.
| S-EPMC9713614 | biostudies-literature
Project description:Selectivity of the palladium-catalyzed cross-coupling reactions of heterocycles bearing multiple identical halogens is mainly determined by the relative ease of oxidative addition. This is related to both the energy to distort the carbon halogen bond to the transition-state geometry (related to the CX bond-dissociation energy) and to the interaction between the heterocycle pi* (LUMO) and PdL(2) HOMO (J. Am. Chem. Soc. 2007, 129, 12664). The computed bond dissociation energies of a larger series of halo-heterocycles have been explored with B3LYP and higher accuracy G3B3 calculations. Quantitative trends in bond dissociation energies have been identified for five- and six-membered chloro and bromo substituted heterocycles with N, O, and S heteroatoms.
| S-EPMC2730358 | biostudies-literature
Project description:A DFT-SOFM-RBFNN method is proposed to improve the accuracy of DFT calculations on Y-NO (Y = C, N, O, S) homolysis bond dissociation energies (BDE) by combining density functional theory (DFT) and artificial intelligence/machine learning methods, which consist of self-organizing feature mapping neural networks (SOFMNN) and radial basis function neural networks (RBFNN). A descriptor refinement step including SOFMNN clustering analysis and correlation analysis is implemented. The SOFMNN clustering analysis is applied to classify descriptors, and the representative descriptors in the groups are selected as neural network inputs according to their closeness to the experimental values through correlation analysis. Redundant descriptors and intuitively biased choices of descriptors can be avoided by this newly introduced step. Using RBFNN calculation with the selected descriptors, chemical accuracy (?1 kcal·mol(-1)) is achieved for all 92 calculated organic Y-NO homolysis BDE calculated by DFT-B3LYP, and the mean absolute deviations (MADs) of the B3LYP/6-31G(d) and B3LYP/STO-3G methods are reduced from 4.45 and 10.53 kcal·mol(-1) to 0.15 and 0.18 kcal·mol(-1), respectively. The improved results for the minimal basis set STO-3G reach the same accuracy as those of 6-31G(d), and thus B3LYP calculation with the minimal basis set is recommended to be used for minimizing the computational cost and to expand the applications to large molecular systems. Further extrapolation tests are performed with six molecules (two containing Si-NO bonds and two containing fluorine), and the accuracy of the tests was within 1 kcal·mol(-1). This study shows that DFT-SOFM-RBFNN is an efficient and highly accurate method for Y-NO homolysis BDE. The method may be used as a tool to design new NO carrier molecules.
| S-EPMC3430220 | biostudies-literature
Project description:Adaptive multilevel splitting (AMS) is a rare event sampling method that requires minimal parameter tuning and allows unbiased sampling of transition pathways of a given rare event. Previous simulation studies have verified the efficiency and accuracy of AMS in the calculation of transition times for simple systems in both Monte Carlo and molecular dynamics (MD) simulations. Now, AMS is applied for the first time to an MD simulation of protein-ligand dissociation, representing a leap in complexity from the previous test cases. Of interest is the dissociation rate, which is typically too low to be accessible to conventional MD. The present study joins other recent efforts to develop advanced sampling techniques in MD to calculate dissociation rates, which are gaining importance in the pharmaceutical field as indicators of drug efficacy. The system investigated here, benzamidine bound to trypsin, is an example common to many of these efforts. The AMS estimate of the dissociation rate was found to be (2.6 ± 2.4) × 10(2) s(-1), which compares well with the experimental value.
| S-EPMC5724379 | biostudies-literature
Project description:In this study, the performance of 17 different density functional theory functionals was compared for the calculation of the bond dissociation energy (BDE) values of X-H (X=C, N, O, S) bonds of aromatic compounds. The effect of the size of the basis set (expansions of 6-31(G)) was also assessed for the initial geometry and zero-point energy calculations, followed by the single-point BDE calculations with different model chemistries with the 6-311 + (3df,2p) basis set. It was found that the size of the basis set for geometry optimization has a much smaller effect on the accuracy of BDE than the choice of functional for the following single-point calculations. The M06-2X, M05-2X and M08-HX functionals yielded highly accurate BDE values compared to experimental data (with the average mean unsigned error MUE = 1.2-1.5 kcal mol-1), performing better than any of the other functionals. The results suggest that geometry optimization may be performed with B3LYP functional and a small basis set, whereas the M06-2X, M05-2X and M08-HX density functionals with a suitably large basis set offer the best method for calculating BDEs of ArX-H (X=C, N, O, S) bonds.
| S-EPMC9174704 | biostudies-literature
Project description:A broad collection of technologies, including e.g. drug metabolism, biofuel combustion, photochemical decontamination of water, and interfacial passivation in energy production/storage systems rely on chemical processes that involve bond-breaking molecular reactions. In this context, a fundamental thermodynamic property of interest is the bond dissociation energy (BDE) which measures the strength of a chemical bond. Fast and accurate prediction of BDEs for arbitrary molecules would lay the groundwork for data-driven projections of complex reaction cascades and hence a deeper understanding of these critical chemical processes and, ultimately, how to reverse design them. In this paper, we propose a chemically inspired graph neural network machine learning model, BonDNet, for the rapid and accurate prediction of BDEs. BonDNet maps the difference between the molecular representations of the reactants and products to the reaction BDE. Because of the use of this difference representation and the introduction of global features, including molecular charge, it is the first machine learning model capable of predicting both homolytic and heterolytic BDEs for molecules of any charge. To test the model, we have constructed a dataset of both homolytic and heterolytic BDEs for neutral and charged (-1 and +1) molecules. BonDNet achieves a mean absolute error (MAE) of 0.022 eV for unseen test data, significantly below chemical accuracy (0.043 eV). Besides the ability to handle complex bond dissociation reactions that no previous model could consider, BonDNet distinguishes itself even in only predicting homolytic BDEs for neutral molecules; it achieves an MAE of 0.020 eV on the PubChem BDE dataset, a 20% improvement over the previous best performing model. We gain additional insight into the model's predictions by analyzing the patterns in the features representing the molecules and the bond dissociation reactions, which are qualitatively consistent with chemical rules and intuition. BonDNet is just one application of our general approach to representing and learning chemical reactivity, and it could be easily extended to the prediction of other reaction properties in the future.
| S-EPMC8179073 | biostudies-literature
Project description:Energetics of reactive intermediates underlies their reactivity. The availability of these data provides a rational basis for understanding and predicting a chemical reaction. We reported here a comprehensive computational study on the energetics of enamine intermediates that are fundamental in carbonyl chemistry. Accurate density functional theory (DFT) calculations were performed to determine the bond energies of enamines and their derived radical intermediates. These efforts led to the compilation of a database of enamine energetics including a thermodynamic index such as free-energy stability, bond dissociation energy (BDE), and acid dissociation constant (pK a) as well as a kinetic index such as nucleophilicity and electrophilicity. These data were validated by relating to experimentally determined parameters and their relevance and utility were discussed in the context of modern enamine catalysis. It was found that pK a values of enamine radical cations correlated well with redox potentials of their parent enamines, the former could be used to rationalize the proton-transfer behavior of enamine radical cations. An analysis of the BDE of enamine radical cations indicated that these species underwent facile β-C-H hydrogen transfer, in line with the known oxidative enamine catalysis. The enamine energetics offers the possibility of a systematic evaluation of the reactivities of enamines and related radicals, which would provide useful guidance in exploring new enamine transformations.
| S-EPMC8867489 | biostudies-literature
Project description:The combination of ion/ion chemistry with commercially available ion mobility/mass spectrometry systems has allowed rich structural information to be obtained for gaseous protein ions. Recently, the simple modification of such an instrument with an electrospray reagent source has allowed three-dimensional gas-phase interrogation of protein structures through covalent and noncovalent interactions coupled with collision cross section measurements. However, the energetics of these processes have not yet been studied quantitatively. In this work, previously developed Monte Carlo simulations of ion temperatures inside traveling wave ion guides are used to characterize the energetics of the transition state of activated ubiquitin cation/sulfo-benzoyl-HOAt reagent anion long-lived complexes formed via ion/ion reactions. The ΔH‡ and ΔS‡ of major processes observed from collisional activation of long-lived gas-phase ion/ion complexes, namely collision induced unfolding (CIU), covalent bond formation, or neutral loss of the anionic reagent via intramolecular proton transfer, were determined. Covalent bond formation via ion/ion complexes was found to be significantly lower energy compared to unfolding and bond cleavage. The ΔG‡ values of activation of all three processes lie between 55 and 75 kJ/mol, easily accessible with moderate collisional activation. Bond formation is favored over reagent loss at lower activation energies, whereas reagent loss becomes competitive at higher collision energies. Though the ΔG‡ values between CIU of a precursor ion and covalent bond formation of its ion/ion product complex are comparable, our data suggest covalent bond formation does not require extensive isomerization.
| S-EPMC9248411 | biostudies-literature