Pople Basis Sets Research Articles

Proton permeation and selective separation of hydrogen isotopes through fullerene nanocages (X12Y12): A DFT insights

To resolve the problem of mobile electrolytes for the permeation of hydrogen isotopes in fuel cells, we have studied herein heteronuclear fullerene nanocages (Al12N12, Al12P12, B12N12, and B12P12) for the separation and permeation of hydrogen isotopes. The separation and permeation of fullerene nanocages for hydrogen isotopes are studied at ωB97XD functional of DFT along with 6–31g (d,p) pople basis set. Zero-point energy (ZPE) is calculated for Protium (H+) and its heavier isotopes including deuterium (D+) and tritium (T+). The Arrhenius equation is employed to calculate the selectivity of protium to deuterium (H+/D+) and deuterium to tritium (H+/T+) isotopes. The selectivity is calculated by using the zero-point energy difference (ZPE) of protium (H+) to its heavier isotopes. The proposed fullerene nanocages estimated better selectivity for proton isotopes. The phosphide-based nanocages (Al12P12, B12P12) provide superior selectivity values of H+/D+ and H+/T+ compared to nitrogen-based nanocages (Al12N12, B12N12). This study provides insights into separation pathways for proton isotopes through nanocages for many industrial applications in the energy sector.

LibERI-A portable and performant multi-GPU accelerated library for electron repulsion integrals via OpenMP offloading and standard language parallelism.

A portable and performant graphics processing unit (GPU)-accelerated library for electron repulsion integral (ERI) evaluation, named LibERI, has been developed and implemented via directive-based (e.g., OpenMP and OpenACC) and standard language parallelism (e.g., Fortran DO CONCURRENT). Offloaded ERIs consist of integrals over low and high contraction s, p, and d functions using the rotated-axis and Rys quadrature methods. GPU codes are factorized based on previous developments [Pham et al., J. Chem. Theory Comput. 19(8), 2213-2221 (2023)] with two layers of integral screening and quartet presorting. In this work, the density screening is moved to the GPU to enhance the computational efficacy for large molecular systems. The L-shells in the Pople basis set are also separated into pure S and P shells to increase the ERI homogeneity and reduce atomic operations and the memory footprint. LibERI is compatible with any quantum chemistry drivers supporting the MolSSI Driver Interface. Benchmark calculations of LibERI interfaced with the GAMESS software package were carried out on various GPU architectures and molecular systems. The results show that the LibERI performance is comparable to other state-of-the-art GPU-accelerated codes (e.g., TeraChem and GMSHPC) and, in some cases, outperforms conventionally developed ERI CUDA kernels (e.g., QUICK) while fully maintaining portability.

Exploring the molecular solvatochromism, stability, reactivity, and non-linear optical response of resveratrol.

This work analyzes the isomerization effects and solvent contributions to the stability, electronic excitations, reactivity, and non-linear optical properties (NLO) of resveratrol molecules within the formalism of the Density Functional Theory. The findings suggest that resveratrol solvatochromism is significantly influenced by solvent polarization. The electronic and free energies (E and G) indicate that trans is the most stable conformer. The system is classified as a strong nucleophile. However, the analysis of the Fukui functions and the Mulliken charges indicate that cis-trans isomerization jointly affects the reactive indices of the carbon and hydrogen atoms. The results also suggest that solvent is relevant to solvatochromism and the NLO response. Both cis and trans conformers present strong excitations that undergo a visible hypsochromic change when the polarity of the solvent increases. Once the absorption spectra are connected to the first hyperpolarization ( ) by the Oudar and Chemla relation, the hypsochromism of resveratrol is the reason for the drop in the generation of the second harmonic when the ambient polarity decreases. The CAM-B3LYP DFT results suggest that resveratrol is interesting for NLO applications. Depending on the choice of solvent, values 50 times those observed for urea ( esu), which is a standard NLO material. The optimized geometries of cis and trans isomers of resveratrol in vacuum were obtained using Density Functional Theory (DFT) with the hybrid exchange-correlation function (CAM-B3LYP) and Pople basis set functions, specifically 6-311++G(d,p). The solvent effect on the geometries of both isomers was included using the polarizable continuum model (PCM) with the same level of QM calculation. Vibrational analysis was conducted to confirm that all optimized geometries correspond to the minimum energy. Various electronic properties, including dipole moments, molecular orbitals, transition energy, dipole polarizabilities, and global reactivity parameters, were calculated using both continuum and discrete solvation models based on the sequential QM/MM methodology. All QM calculations were performed with the Gaussian 09 program and the MC simulations with the DICE program. All NLO analysis was carried out using the Multiwfn code.

Computational Study of the Kinetics and Mechanisms of Gas-Phase Decomposition of N-Diacetamides Using Density Functional Theory.

In this research work, we examined the decomposition mechanisms of N-substituted diacetamides. We focused on the substituent effect on the nitrogen lone-pair electron delocalization, with electron-withdrawing and electron donor groups. DFT functionals used the following: B1LYP, B3PW91, CAMB3LYP, LC-BLYP, and X3LYP. Dispersion corrections (d3bj) with Becke-Johnson damping were applied when necessary to improve non-covalent interactions in the transition state. Pople basis sets with higher angular moments and def2-TZVP basis sets were also applied and were crucial for obtaining consistent thermodynamic parameters. The proposed mechanism involves a six-membered transition state with the extraction of an α hydrogen. Several conformers of N-diacetamides were used to account for the decrease in entropy in the transition state in the rate-determining state. All calculations, including natural bond orbital (NBO) analyses, were performed using the Gaussian16 computational package and its GaussView 6.0 visualizer, along with VMD and GNUPLOT software. The isosurfaces and IBSIs were calculated using MultiWFN and IGMPlot, respectively.

Assessment of the usefulness of β-cyclodextrin in the removal of progesterone from the environment

The paper presents the results and interpretation of theoretical calculations for the progesterone–β-cyclodextrin (P@β-CD; G4MP2) and progesterone–β-cyclodextrin–β-cyclodextrin (P@β-CD–β-CD; G2) systems. The geometry of the progesterone molecule was optimized on basis of the DFT theory using the B3LYP, PBE1PBE and M06-2X functionals, for selected Pople basis sets [6-31G, 6‐31++G, 6‐31++G(d,p), 6-311G, 6‐311++G, 6‐311++G(d,p)] and the Dunning basis set (aug-cc-pVDZ). Presented results of theoretical calculations provide insight into the mechanism of formation of supramolecular systems of progesterone and β-cyclodextrin, allowing us to understand the impact of differences in the polarizability of specific fragments of the “guest” molecule, through the analysis of the Mulliken population distribution, on the tendency for equilibrium inclusion by “host” molecules characterized by selective affinity towards hydrophilic and hydrophobic molecular systems. The comparison of model structures of “guest–host” systems with 1:1 and 1:2 stoichiometry allows us to assess the contribution of a given type of non-covalent interactions (hydrogen bonds, van der Waals interactions, London dispersion forces) in the formation of supramolecular complexes. Due to the complexity of the real research object, including interactions between large molecules in solution, in this work it was decided to choose an approximation that reflects only the key effects. The results presented in this paper constitute a starting point for both theoretical research on analogous supramolecular systems based on β-cyclodextrin derivatives, as well as instrumental studies (NMR, FT–IR, HPLC) of preparations of real inclusion complexes. The approach adopted now can contribute to a better understanding of the phenomenon used in the development of water purification technologies.

Ab initio investigation on the mechanism of SO2 activation by P/B intermolecular frustrated Lewis pairs.

In silico study investigates the activation of sulfur dioxide by newly designed frustrated Lewis pairs, i.e., [P(tBu)3…B(C2NBSHF2)3], where the Lewis acid part is a super Lewis acid. The activation process involves the making of P-S and B-O bonds, leading to the formation of an FLP-SO2 adduct. The calculated results demonstrate that the activation of SO2 by the FLP is almost barrierless and exothermic. Exploration of the impact of the solvent environment on the feasibility and energetics of the reaction has been investigated. The exothermicity is increasing in nonpolar solvents. This study focuses on understanding the electronic activity of SO2 activation by FLP with the help of the Minnesota 06 functional, M06-2X (global hybrid functional with 54% HF exchange) along with Pople's basis set, 6-311G (d, p). Principal interacting orbital and extended transition state-natural orbitals for chemical valence studies, giving impactful insight into the favorable orbital interaction and electron transfer in this reaction. Furthermore, useful CDFT descriptors such as reaction force constant and reaction electronic flux profiles along the intrinsic reaction coordinate give insights into the synchronicity and total electronic activity of the reaction.

Ab Initio and Kinetic Modeling of β-d-Xylopyranose under Fast Pyrolysis Conditions.

Lignocellulosic biomass is an abundant renewable resource that can be upgraded to chemical and fuel products through a range of thermal conversion processes. Fast pyrolysis is a promising technology that uses high temperatures and fast heating rates to convert lignocellulose into bio-oils in high yields in the absence of oxygen. Hemicellulose is one of the three major components of lignocellulosic biomass and is a highly branched heteropolymer structure made of pentose, hexose sugars, and sugar acids. In this study, β-d-xylopyranose is proposed as a model structural motif for the essential chemical structure of hemicellulose. The gas-phase pyrolytic reactivity of β-d-xylopyranose is thoroughly investigated using computational strategies rooted in quantum chemistry. In particular, its thermal degradation potential energy surfaces are computed employing Minnesota global hybrid functional M06-2X in conjunction with the 6-311++G(d,p) Pople basis set. Electronic energies are further refined by performing DLPNO-CCSD(T)-F12 single-point calculations on top of M06-2X geometries using the cc-pVTZ-F12 basis set. Conformational analysis for minima and transition states is performed with state-of-the-art semiempirical quantum chemical methods coupled with metadynamics simulations. Key thermodynamic quantities (free energies, barrier heights, enthalpies of formation, and heat capacities) are computed. Rate coefficients for the initial steps of thermal decomposition are computed by means of reaction rate theory. For the first time, a detailed elementary reaction kinetic model for β-d-xylopyranose is developed by utilizing the thermodynamic and kinetic information acquired from the aforementioned calculations. This model specifically targets the initial stages of β-d-xylopyranose pyrolysis in the high-pressure limit, aiming to gain a deeper understanding of its reaction kinetics. This approach establishes a systematic strategy for exploring reactive pathways, evaluating competing parallel reactions, and selectively accepting or discarding pathways based on the analysis. The findings suggest that acyclic d-xylose plays a significant role as an intermediary in the production of key pyrolytic compounds during the pyrolysis of xylose. These compounds include furfural, anhydro-d-xylopyranose, glycolaldehyde, and dihydrofuran-3(2H)-one.

Potential energy surfaces of the [formula omitted] clusters in solution phase

In the present work, we have studied the effects of the implicit solvent on the structural and energetic parameters of Cu2+(MeOH)n=1-8 clusters. Initially, we explored the potential energy surfaces (PESs) obtained in the solution phase using the functional DFT M06-2X method, in conjunction with the 6–31++G** Pople basis set, and compared them with isomers of the same complexes optimized in the gaseous phase. The implicit solvent is described here by the polarized continuum model of integral equation formalism (IEF-PCM). Subsequently, we reported the temperature dependence and the binding electronic energies per methanol molecule. Our investigations reveal that the 6- and 5-coordinate isomers occupy the ground states on the PESs of the Cu2+(MeOH)n=1-8 clusters. The stability rules of the Cu2+ ion established in the gas phase are best adhered to in the implicit solvent. Regarding the temperature dependence of the solution phase conformer distribution, symmetric and compact structures, as well as those with a high coordination number (CN), dominate the population of the Cu2+(MeOH)n=1-8 clusters at all temperatures. Notably, the populations of clusters with sizes n=7-8 are exclusively dominated by hexa-coordinated structures at all temperatures, while penta- and tetra-coordinated isomers are disadvantaged. For the determination of the binding electronic energy, using a function fit, we computed these energies to be -1589kJ mol−1 in the implicit solvent and -7150kJ mol−1 in the gaseous phase. This significant difference is attributable to the interactions between the methanol molecules and the implicit solvent, which are present in the solvent phase but absent in a gaseous medium. This difference is further supported by the analysis of the Wiberg bond index of the metal–ligand interactions of the clusters in the two media.

Quantum-Chemical Simulation of 13C NMR Chemical Shifts of Fullerene C60 Exo-Derivatives

The 13C NMR chemical shifts of fullerene C60 exo-derivatives were calculated using quantum chemical hybrid functionals combined with Pople, Dunning correlation-consistent, and def2-TZVP split valence basis sets taking into account the solvent effect (polarizable continuum model). A relationship between theoretical and experimental 13C NMR chemical shifts (CSs) is assessed quantitatively to select a functional/basis set combination. It is found that the CAM-B3LYP/6-31G and M06L/6-31G combinations have the best convergence with experimental data in modeling the 13С NMR CSs of sp3 fullerene carbon atoms in С60 derivatives, whereas X3LYP/6-31G and CAM-B3LYP/6-31G(d) in modeling the 13С NMR CSs of their sp2 fullerene carbon atoms.

Multiple Intramolecular Hydrogen Bonding in Large Biomolecules: DFT Calculations and Deuterium Isotope Effects on 13C Chemical Shifts as a Tool in Structural Studies

Large biomolecules often have multiple intramolecular hydrogen bonds. In the cases where these interact, it requires special tools to disentangle the patterns. Such a tool could be deuterium isotope effects on chemical shifts. The use of theoretical calculations is an indispensable tool in such studies. The present paper illustrates how DFT calculations of chemical shifts and deuterium isotope effects on chemical shifts in combination with measurements of these effects can establish the complex intramolecular hydrogen bond patterns of rifampicin as an example) The structures were calculated using DFT theoretical calculations, performed with the Gaussian 16 software. The geometries were optimized using the B3LYP functional and the Pople basis set 6-31G(d) and the solvent (DMSO) was taken into account in the PCM approach. Besides the 6-31G(d) basis set, the 6-31 G(d,p) and the 6-3111G(d,p) basis sets were also tested. The nuclear shieldings were calculated using the GIAO approach. Deuteriation was simulated by shortening the X-H bond lengths by 0.01 Å.

Theoretical studies of interactions in cyprodinil–α-cyclodextrin and cyprodinil–β-cyclodextrin systems

The paper presents the results of theoretical calculations in terms of the G4MP2 composite method for cyprodinil–α-cyclodextrin (C 0@α-CD) and cyprodinil–β-cyclodextrin (C 0@β-CD) systems. Studies also covered analogous systems consisting of the anion (C −) and the cation (C +) of cyprodinil. The geometries of the cyprodinil molecule and ions were optimized on the basis of the DFT theory, using hybrid (B3LYP, PBE0), pure (B97-D) and “meta” (M06-2X) GGA functionals for selected Pople basis sets [6-311++G(d,p), 6-311++G(2d,p), 6-311++G(2d,2p)] and Dunning basis set (aug-cc-pVDZ). The research results suggest that the affinity of “guest” molecules for “hosts” is relatively low. Theoretical studies of the “guest-host” systems allow to predict the properties of the designed preparations.

Performance of Functionals and Basis Sets in Calculating Redox Potentials of Nitrile Alkenes and Aromatic Molecules using Density Functional Theory

AbstractTo calculate redox potentials of small organic compounds in acetonitrile, the performance of different functionals and basis sets has been evaluated using density functional theory with the polarized continuum model. Several Pople and Ahlrichs basis sets and functionals at distinct levels of the Jacob's Ladder have been studied: PW91, PBE, M06‐L, B3LYP, PBE0, M06‐2X, CAM‐B3LYP, ωB97X‐D3, RI‐B2PLYP. It is shown that redox studies should not be done considering oxidation and reduction potentials jointly but analyzing them separately. Functional M06‐2X has a more consistent and uniform response both in reductions and oxidations, having an adequate balance between accuracy and computational effort. Linear dependence between the basis set size and the accuracy of the calculated redox potential has not been found.

A Theoretical Study of Hydrogen Abstraction Reactions in Guanosine and Uridine.

All practically possible hydrogen abstraction reactions for guanosine and uridine have been investigated through quantum chemical calculations of energy barriers and rate constants. This was done at the level of density functional theory (DFT) with the ωB97X-D functional and the 6-311++G(2df,2pd) Pople basis set. Transition state theory with the Eckart tunneling correction was used to calculate the rate constants. The results show that the reaction involving the hydrogen labelled C4' in the ribofuranose part has the largest rate constant for guanosine with the value 5.097×1010 L mol-1s-1 and the largest for uridine with the value 1.62×1010 L mol-1s-1. Based on the results for these two nucleosides, there is a noticeable similarity between the rate constants in the ribofuranose part of the molecule, even though they are bound to two entirely different nucleobases.

Population analysis and the effects of Gaussian basis set quality and quantum mechanical approach: main group through heavy element species.

Atomic charge and its distribution across molecules provide important insight into chemical behavior. Though there are many studies on various routes for the determination of atomic charge, there are few studies that examine the broader impact of basis set and quantum method used over many types of population analysis methods across the periodic table. Largely, such a study of population analysis has focused on main-group species. In this work, atomic charges were calculated using several population analysis methods including orbital-based methods (Mulliken, Löwdin, and Natural Population Analysis), volume-based methods (Atoms-in-Molecules (AIM) and Hirshfeld), and potential derived charges (CHELP, CHELPG, and Merz-Kollman). The impact of basis set and quantum mechanical method choices upon population analysis has been considered. The basis sets utilized include Pople (6-21G**, 6-31G**, 6-311G**) and Dunning (cc-pVnZ, aug-cc-pVnZ; n = D, T, Q, 5) basis sets for main group molecules. For the transition metal and heavy element species examined, relativistic forms of the correlation consistent basis sets were used. This is the first time the cc-pVnZ-DK3 and cc-pwCVnZ-DK3 basis sets have been examined with respect to their behavior across all levels of basis sets for atomic charges for an actinide. The quantum methods chosen include two density functional (PBE0 and B3LYP), Hartree-Fock, and second-order Møller-Plesset perturbation theory (MP2) approaches.

Computational study of physicochemical, optical, and thermodynamic properties of 2,2-dimethylchromene derivatives.

A large number of heterocyclic compounds are used as drugs, mainly due to the duality of lipophilicity playing in hydrophobic interactions and solubility with at least one hydrogen bond acceptor. The study of electronic properties is then important to better understand not only these charge distribution effects but also some other physicochemical properties involved in bioactivity to directly assess the bioavailability of these compounds and a possible classification in related applications. Phytomolecules such as chromenes are very accessible molecules exhibiting a bioactivity. Our study is focused on the impact of a number of functional groups acting on some 2,2-dimethylchromene derivatives, namely their global reactivity from the frontier molecular orbitals and local reactivity from the Fukui functions, where the carbonyl group acting as an electron withdrawal group has the most relevant effect, the solubility from the partition coefficient Log P strongly depending on the charge distribution and electronegative sites, the optical effects from the delocalization in the vinyl group, as well as the evaluation of the entropy associated with the molecular flexibility also acting on the bioactivity. Despite the effects of the wave function or density methods on the order of magnitude of these properties, these compounds are consistent with the rules for a potential oral drug candidate. The calculations of the electronic properties were performed through two levels of theory: Hartree-Fock level as a wave function-based method as an ab initio reference including some physically consistent eigenvalues and density functional theory DFT as a correlation consistent method using different functionals: hybrid or with a long-range correction. The basis set used is a 6-311++G(d,p) Pople basis set including diffuse and polarization basis functions. The basis set is adapted to the size of the molecules and consequently to the degree of electronic localization. Gaussian 09 software was used for the computation.

Cation functional group effect on SO2 absorption in amino acid ionic liquids.

Introduction: The effect of the functional group of the cation on SO2 acidic gas absorption by some designed amino acid ionic liquids (AAILs) was studied. Methods: An isolated pair of glycinate anion and pristine imidazolium-based cation, as well as decorated cation functionalized by hydroxyl (OH), amine (NH2), carboxylic acid (COOH), methoxy (OCH3), and acetate (CH3COO) groups, were structurally optimized by density functional theory (DFT) using split-valence triple-zeta Pople basis set. Results and Discussion: The binding and Gibbs free energy (ΔGint) values of SO2 absorption show the AAIL functionalized by the COOH group is the most thermodynamically favorable green solvent and this functional group experiences the closest distance between anion and captured SO2 and vice versa in the case of cation … SO2 which may be the main reason for being the best absorbent; in addition, the highest net charge-transfer amount of SO2 is observed. Comparing the non-covalent interaction of the systems demonstrates that the strongest hydrogen bond between captured gas and anion, as well as π-hole, and van der Waals (vdW) interaction play critical roles in gas absorption; besides, the COOH functional group decreases the steric effect while the CH3COO functional group significantly increases steric effect after absorption that declines the hydrogen bond.

Benchmarking computational chemistry approaches on iminodiacetic acid

In the present study, theoretical harmonic vibrational frequencies (IR and Raman), carbon and proton NMR chemical shifts, geometric parameters, atomic charges (only for heteroatoms), reactivity indices (eLUMO, eHOMO, electronegativity, and hardness), and thermodynamical data (inner energy, enthalpy, Gibbs free energy, and entropy) of iminodiacetic acid (IDA) molecule have been investigated. We utilized ORCA software for B3LYP and HF (combined with Pople and Karlsruhe basis sets) calculations and MOPAC2016 software for semi-empirical calculations (AM1, PM3, and PM6). Theoretical vibrational frequencies and carbon and proton NMR chemical shifts have been compared with the corresponding experimental data. Although there was a strong correlation between the experimental and computational vibrational frequencies at low frequencies (<2200 cm−1), the computational predictions of vibrational frequencies were unsuccessful at high frequencies (>2200 cm−1). Distinctly, the studied computational approaches appeared to perform better in the prediction of carbon and proton NMR chemical shifts. Theoretical vibrational frequencies were also compared to each other to understand the impact of method choice (HF vs B3LYP D3 vs semi-empirical methods), dispersion correction (B3LYP D3 vs B3LYP), water solvation (SMD supplemented vs non-supplemented calculations), the family of basis set (Pople vs Karlsruhe basis sets), numbers of zeta (double vs triple zeta), polarization function (polarized vs nonpolarized basis sets), and diffusion function (diffusion supplemented vs non-supplemented basis sets). Moreover, geometric parameters, heteroatom charges, reactivity indices, and thermodynamical data produced by distinct computational approaches, as well, were compared to each other. Based on these comparisons, we detected critical factors (such as water solvation) acting on the computation of geometries, energies, and charges.

A computational study of the inclusion of β-cyclodextrin and nicotinic acid: DFT, DFT-D, NPA, NBO, QTAIM, and NCI-RDG studies.

Forming complexes with cyclodextrins can protect nicotinic acid (vitamin B3) from premature metabolism and enhance the solubility and stability of this drug. In this work, the formation of the inclusion complex of the neutral form of nicotinic acid and β-cyclodextrin was achieved. The complex is modeled using PM3, PM6-D4H3, and PM7, by considering two orientations of the guest: A and B, one is from wide to narrow rim, and the second is from narrow to the wide rim, respectively. The global minima positions were re-optimized using three density function methods: MN-15, B3LYP, and PW6B95-D3 with polarized Pople basis set 6-31G(d) in gas and aqueous phase. Orientation A showed the minimum complexation energy where the carboxylic functional group of nicotinic acid is located on the primary hydroxyl rim of β-cyclodextrin and the pyridine ring is totally embedded in the cavity. To further our study on the nature of complexation and the interactions of this host-guest system, different calculations were done. The reactivity indices showed that orientation A is harder than B and more electrophilic; the charge transfer occurred from the host to the guest and was confirmed by the natural population analysis (NPA). The natural bond orbitals (NBO) reveal the delocalization of orbitals between the host and the guest, quantum theory of atoms in molecules (QTAIM) analysis, and non-covalent interaction (NCI) analysis based on a reduced density gradient (RDG) give a detailed description of the nature of interactions between the host and the guest such as the hydrogen bonding and van der Waals interaction, and confirmed the stability of the complex given by the orientation A.

Quantum chemical hydrogenolysis strategy for elimination of heteroatoms in biomass homologous organic compounds based on oxolane and thiolane

Bio-oils obtained from biomass contain heteroatoms compounds, like oxolane and thiolane. It is quite difficult for industrialist to purify such refractory bio-oils. One of the efficient strategies for the elimination of heteroatoms is hydrogenolysis process, which results in the formation of H2O and H2S residues as by-products. In this work, quantum chemical studies have been used to analyse the reaction mechanism for the removal of hetero atoms (S and O) as H2O and H2S. We selected B3LYP functional of DFT with Pople's basis set 6-311G(d,p) for computing the hydrogenolysis steps without catalyst. LANL2DZ basis set, is used for studying hydrogenolysis steps involving AlCl3 and WS3H3+ as catalysts. All the reactions are analysed at the temperature of 600 K and pressure of 40 bars. Structural, thermodynamic, kinetic properties have been employed to study this process. The analysis of variations parameters during the hydrogenolysis process reveals that these two organic biomass compounds undergo sequential ring opening at C-X (X = O, S) bonds. Butanol and Butanethiol are obtained as a result of first hydrogenolysis process, and these compounds are converted to butane during second catalytic process while eliminating heteroatoms.

Benchmark Density Functional Theory Approach for the Calculation of Bond Dissociation Energies of the M-O2 Bond: A Key Step in Water Splitting Reactions.

A very fascinating aspect in quantum chemical research is to determine the accurate and cost-effective methods for the calculation of electronic and structural properties through a benchmark study. The current study focuses on the performance evaluation of density functional theory methods for the accurate measurement of bond dissociation energies (BDEs) of chemically important M–O2 bonds in water splitting reactions. The BDE measurement has got noteworthy attention due to its importance in all areas of chemistry. For BDE measurements of M–O2 bonds in five metal complexes with oxygen molecules, 14 density functionals (DFs) are chosen from seven classes of DFs with two series of mixed basis sets. A combination of pseudopotential and Pople basis sets [LANL2DZ & 6-31G(d) and SDD & 6-31+G(d)] are used as a series of mixed basis sets. The B3LYP-GD3BJ functional with LANL2DZ & 6-31G(d) gives outstanding results due to low deviations, error, and the best Pearson’s correlation (R) between the experimental and theoretical data. Our study suggested an efficient, low-cost, precise, and accurate B3LYP-GD3BJ/LANL2DZ & 6-31G(d) level of theory for BDE of the M–O2 bond, which may be useful for chemists working in the field of energy generation and utilization.