Quantum Chemical Methods Research Articles

Exploiting Non-Abelian Point-Group Symmetry to Estimate the Exact Ground-State Correlation Energy of Benzene in a Polarized Split-Valence Triple-Zeta Basis Set.

Local electronic-structure methods in quantum chemistry operate on the ability to compress electron correlations more efficiently in a basis of spatially localized molecular orbitals than in a parent set of canonical orbitals. However, many typical choices of localized orbitals tend to be related by selected, near-exact symmetry operations whenever a molecule belongs to a point group, a feature which remains largely unexploited in most local correlation methods. The present Letter demonstrates how to leverage a recent unitary protocol for enforcing symmetry properties among localized orbitals to yield a high-accuracy estimate of the exact ground-state correlation energy of benzene (D6h) in correlation-consistent polarized basis sets of both double- and triple-ζ quality. Through an initial application to many-body expanded full configuration interaction (MBE-FCI) theory, we show how molecular point-group symmetry can lead to computational savings that are inversely proportional to the order of the point group in a manner generally applicable to the acceleration of modern local correlation methods.

Molecular Aspects of the Interactions between Selected Benzodiazepines and Common Adulterants/Diluents: Forensic Application of Theoretical Chemistry Methods.

Benzodiazepines are frequently encountered in crime scenes, often mixed with adulterants and diluents, complicating their analysis. This study investigates the interactions between two benzodiazepines, lorazepam (LOR) and alprazolam (ALP), with common adulterants/diluents (paracetamol, caffeine, glucose, and lactose) using infrared (IR) spectroscopy and quantum chemical methods. The crystallographic structures of LOR and ALP were optimized using several functionals (B3LYP, B3LYP-D3BJ, B3PW91, CAM-B3LYP, M05-2X, and M06-2X) combined with the 6-311++G(d,p) basis set. M05-2X was the most accurate when comparing experimental and theoretical bond lengths and angles. Vibrational and 13C NMR spectra were calculated to validate the functional's applicability. The differences between LOR's experimental and theoretical IR spectra were attributed to intramolecular interactions between LOR monomers, examined through density functional theory (DFT) optimization and quantum theory of atoms in molecules (QTAIM) analysis. Molecular dynamics simulations modeled benzodiazepine-adulterant/diluent systems, predicting the most stable structures, which were further analyzed using QTAIM. The strongest interactions and their effects on IR spectra were identified. Comparisons between experimental and theoretical spectra confirmed spectral changes due to interactions. This study demonstrates the potential of quantum chemical methods in analyzing complex mixtures, elucidating spectral changes, and assessing the structural stability of benzodiazepines in forensic samples.

Reaction-Free Energies for Complexation of Carbohydrates by Tweezer Diboronic Acids.

The accurate calculation of reaction-free energies (ΔrG°) for diboronic acids and carbohydrates is challenging due to reactant flexibility and strong solute-solvent interactions. In this study, these challenges are addressed with a semiautomatic workflow based on quantum chemistry methods to calculate conformational free energies, generate microsolvated solute structural ensembles, and compute ΔrG°. Workflow parameters were optimized for accuracy and precision while controlling computational costs. We assessed the accuracy by studying three reactions of diboronic acids with glucose and galactose, finding that the conformational entropy contributes significantly (by 3-5 kcal/mol at room temperature). Explicit solvent molecules improve the computed ΔrG° accuracy by about 4 kcal/mol compared to experimental data, though using 13 or more water molecules reduced precision and increased computational overhead. After fine-tuning, the workflow demonstrated remarkable accuracy, with an absolute error of about 2 kcal/mol compared to experimental ΔrG° and an average interquartile range of 2.4 kcal/mol. These results highlight the workflow's potential for designing and screening tweezer-like ligands with tailored selectivity for various carbohydrates.

Effect of adsorption of two green biopolymers on the corrosion of aluminum in 1.0 M NaCl solution: Physicochemical, spectroscopic, surface and quantum chemical exploration

Although aluminum is an extremely versatile material and can be manufactured and used for a wide range of industrial purposes, it has a limited resistance to corrosion. Hence, the efficacy of two linear biopolymers, keratan (Ker) and chitosan (Chi), was explored, for the first time, as ecologically safe metallic corrosion inhibitors (CIs) for Al in 1.0 M NaCl solution at different temperatures. For performing this study, several experimental and theoretical methods were applied. By combining these methods, it was confirmed that the investigated biopolymers were set to be effective inhibitors for Al corrosion in 1.0 M NaCl solution. Using a 500 mg/L dose of CIs, inhibition efficiencies (% IEs) of 92.7 % and 86.4 % were achieved for Ker and Chi, respectively, indicating their strong adsorption on the metal surface. Thermodynamic and kinetic investigations signified that the adsorption was founed to be physical and obeyed the Langmuir isotherm. The kinetics and mechanism of Al corrosion as well as its inhibition by the tested biopolymers were inspected. The results demonstrated that these biopolymers behaved as mixed-type and interface-type inhibitors. To investigate the inhibition mechanism, theoretical quantum chemical methods were also applied. The experimental investigations were validated by theoretical computations.

Vibrational and Rovibrational Spectroscopic Data for the Ground and First-Excited States of Phosgene (COCl2), Formic Acid (HCOOH), and Chloroformic Acid (ClCOOH).

While "black box" quantum chemical computations for the determination of rovibronic spectral data are not quite at hand, the present work utilizes the titular molecules to showcase how excited-state quantum chemical methods can be conjoined to quartic force field (QFF) anharmonic rovibrational treatments to provide novel and useful predictions for such data. This work employs hybrid QFFs with explicitly correlated coupled cluster theory along with the equation-of-motion formalism to generate harmonic force constants and time-dependent density functional theory (TD-DFT) to produce anharmonic force constants for the generation of electronically excited-state rovibrational spectral data, in effect, rovibronic spectral data. Specific spectroscopic results from this work show that the fundamental C═O stretch in phosgene as well as in cis- and trans-formic acid drop from the region of around 1800 cm-1 to close to 1100 cm-1 or less in the first excited states of each molecule. While such is expected for these n → π* excitations, this work provides quantitative predictions for these fundamental vibrational frequencies. The most notable theoretical result is that the TD-DFT-based QFFs can experience unexpected failures, and their inclusion in excited-state hybrid QFFs should require at least two functionals to be employed. The computation of DFT QFFs is relatively fast, and such a "doubling up" of the QFFs will not greatly increase the computational time.

Pyrolytic conversion of glucose into hydroxymethylfurfural and furfural: Benchmark quantum-chemical calculations.

Quantum chemical methods have been intensively applied to study the pyrolytic conversion of glucose into hydroxymethylfurfural (HMF) and furfural (FF). Herein, we collect the most relevant mechanistic proposals from the recent literature and organize them into a single reaction network. All the transition structures (TSs) and intermediates are characterized using highly accurate ab initio methods and the possible reaction pathways are assessed in terms of the Gibbs energies of the TSs and intermediates with respect to β-glucopyranose, selecting a 2D ideal-gas standard state at 773 K to represent the pyrolysis conditions. Several pathways can lead to the formation of both HMF and FF passing through rate-determining TSs that have ΔG‡ values of ~49-50 kcal/mol. Both water-assisted mechanisms and nonspecific environmental effects have a minor impact on the Gibbs energy profiles. We find that the HMF → FF + CH2O fragmentation has a small ΔrxnG value and an accessible ΔG‡ barrier. Our computational results, which are in consonance with the kinetic parameters derived from lumped models, the results of isotopic labeling experiments and the reported HMF/FF molecular ratios, could be useful for modeling studies including on nonequilibrium kinetic effects that may render more information about product yields and the relevance of the various pathways.

Structural and Electronic Evolution of Ethanolamine upon Microhydration: Insights from Hyperfine Resolved Rotational Spectroscopy

AbstractEthanolamine hydrates containing from one to seven water molecules were identified via rotational spectroscopy with the aid of accurate quantum chemical methods considering anharmonic vibrational corrections. Ethanolamine undergoes significant conformational changes upon hydration to form energetically favorable hydrogen bond networks. The final structures strongly resemble the pure (H2O)3–9 complexes reported before when replacing two water molecules by ethanolamine. The 14N nuclear quadrupole coupling constants of all the ethanolamine hydrates have been determined and show a remarkable correlation with the strength of hydrogen bonds involving the amino group. After addition of the seventh water molecule, both hydrogen atoms of the amino group actively contribute to hydrogen bond formation, reinforcing the network and introducing approximately 21–27 % ionicity towards the formation of protonated amine. These findings highlight the critical role of microhydration in altering the electronic environment of ethanolamine, enhancing our understanding of amine hydration dynamics.

Photoionization of aziridine: Nonadiabatic dynamics of the first six low-lying electronic states of the aziridine radical cation.

In this article, the theoretical photoionization spectroscopy of the aziridine (C2H5N) molecule is investigated. To start with, we have optimized the geometry of this molecule at the neutral electronic ground state at the density functional theory/augmented correlation-consistent polarized valence triple zeta level of theory using the G09 program. The electronic structure calculations were restricted to the first six low-lying electronic states in order to account for the experimental photoelectron spectrum of the C2H5N molecule. The first six low-lying electronic states (X̃2A', Ã2A', B̃2A″, C̃2A″, D̃2A', and Ẽ2A') of the potential energy surfaces (PESs) are calculated by both equation of motion-ionization potential-coupled cluster singles and doubles and multi-configuration quasi-degenerate perturbation theory abinitio quantum chemistry methods along the dimensionless normal displacement coordinates in which multiple conical intersections were established among the considered electronic states. A (6 × 6) model vibronic Hamiltonian is constructed on a diabatic electronic basis, using the symmetry selection rules and Taylor series expansion. The Cs symmetry point group of the aziridine molecule leads to electronic states symmetry of either A' or A″, and these states are close in energy, due to which the same symmetry electronic states avoid each other. To get a smooth diabatic PES, a fourfold diabatization scheme is used, which is implemented in the General Atomic and Molecular Electronic Structure Systems suite of programs. All the parameters used in the diabatic vibronic coupling model Hamiltonian are calculated in terms of the normal modes of vibrational coordinates. Finally, the vibronic model Hamiltonian constructed for the coupled six electronic states is used to solve both time-independent and time-dependent Schrödinger equations using the multi-configuration time-dependent Hartree program module to obtain the dynamical observables. The theoretical vibronic band structure is found to be in good accord with the available experimental results.

CO2 capture using choline chloride-based eutectic solvents. An experimental and theoretical investigation

Eutectic solvents (ESs) have attracted considerable attention as CO2 absorbents due to their tunable and unique properties. In this study, ESs based on choline chloride (ChCl) as hydrogen bond acceptor were combined with four different hydrogen bond donors (HBDs): urea, formamide, monoethanolamine and 1-aminopropan-2-ol, in different molar ratios. The impact of the primary amino group in HBD, with particular focus on the functional groups present in the vicinity, on the physicochemical and structural properties of the ESs was investigated. The CO2 capture in these ESs was measured. 13C NMR, FTIR and Raman analyses were performed to provide evidence of chemical and structural changes after CO2 capture. Quantum chemical methods were used to investigate the interaction mechanism between the ESs and CO2. The resulting ESs were characterised by measuring their basic physicochemical properties, including melting point, density and viscosity as a function of temperature. A three-parameter correlation equation was proposed to predict the solubility of CO2 in ChCl-based physical absorbents. It is shown that the use of the free volume parameter of ESs can help to identify alternative CO2 absorbents.

Temperature dependent stochastic dynamics mass spectrometric analysis of configurationally locked polyenes

This paper derived relation between mass spectrometric variable intensity of analyte peaks and three-dimensional molecular structures of ions toward temperature utilising innovative stochastic dynamics equations, which provide exact analyte quantification and 3D structural analysis and are attractive formulas for modelling of mass spectrometric phenomena. The same equations are opened to develop new relations accounting for temperature; thus, providing argumentation about temperature dependent variables. High-temperature measurement of vapours are important tool, amongst others and use the Arrhenius’s plot. However, there is a variation between experimental and theoretical linear correlation toward temperature; thus, highlighting soft ionisation methods showing significant standard deviations. It could be explained with temperature perturbation of charge values. Descriptive models, so far, provide semi-quantification. This study, first shows validity of new functional model. Straightforwardly, its empirical proof is obtained determining 3D structurally protomers and isotopomers of configurationally locked polyenes via ultra-high resolution electrospray ionization and atmospheric pressure chemical ionization mass spectrometry within the temperature range T=300-450°C. There are used high accuracy quantum chemical static methods, molecular dynamics and chemometrics. There are excellent method performances (|r|=0.99978). The tool could be regarded as an actually honest alternative of available theoretical models.

Structural and Electronic Evolution of Ethanolamine upon Microhydration: Insights from Hyperfine Resolved Rotational Spectroscopy.

Ethanolamine hydrates containing from one to seven water molecules were identified via rotational spectroscopy with the aid of accurate quantum chemical methods considering anharmonic vibrational corrections. Ethanolamine undergoes significant conformational changes upon hydration to form energetically favorable hydrogen bond networks. The final structures strongly resemble the pure (H2O)3-9 complexes reported before when replacing two water molecules by ethanolamine. The 14N nuclear quadrupole coupling constants of all the ethanolamine hydrates have been determined and show a remarkable correlation with the strength of hydrogen bonds involving the amino group. After addition of the seventh water molecule, both hydrogen atoms of the amino group actively contribute to hydrogen bond formation, reinforcing the network and introducing approximately 21-27 % ionicity towards the formation of protonated amine. These findings highlight the critical role of microhydration in altering the electronic environment of ethanolamine, enhancing our understanding of amine hydration dynamics.

Computer-Aided Drug Design Using the Fragment Molecular Orbital Method: Current Status and Future Applications for SBDD.

Owing to the increasing use of computers, computer-aided drug design (CADD) has become an essential component of drug discovery research. In structure-based drug design (SBDD), including inhibitor design and in silico screening of drug target molecules, concordance with wet experimental data is important to provide insights on unique perspectives derived from calculations. Fragment molecular orbital (FMO) method is a quantum chemical method that facilitates precise energy calculations. Fragmentation method makes it possible to apply the quantum chemical method to biological macromolecules for energy calculation based on the electron behavior. Furthermore, interaction energies calculated on a residue-by-residue basis via fragmentation aid in the analysis of interactions between the target and ligand molecule residues and molecular design. In this review, we outline the recent developments in SBDD and FMO methods and highlight the prospects of developing machine learning approaches for large computational data using the FMO method.

Molecular Docking of Chiral Drug Enantiomers With Different Bioactivities.

Chirality has an important role in the drug design because enantiomers may exhibit different bioactivity when interacting with macromolecules of a living organism. In our previous work, based on the analysis of a set of 100 chiral drugs, a relationship was established between the sign of chirality of enantiomers and their bioactivity. To understand the reasons for the observed patterns of chiral specificity of drug enantiomers, the interaction of 10 enantiomeric pairs of chiral drugs with the corresponding target proteins has been considered using molecular docking and further postprocessing by quantum chemistry methods. The data obtained confirm that the energetic aspect of the interaction between opposite enantiomers and target protein affects the enantiomer biological activity. In addition, the results show that molecular docking is able to distinguish between bioactive and inactive/less active enantiomers, although many docking programs are not accurate enough to distinguish a weak inhibitor from a strong one.

An algorithm for very high pressure molecular dynamics simulations.

We describe a method to run simulations of ground or excited state dynamics under extremely high pressures. The method is based on the introduction of a fictitious ideal gas that exerts the required pressure on a molecular sample and is therefore called XP-GAS (eXtreme Pressure by Gas Atoms in a Sphere). The algorithm is most suitable for approximately spherical clusters of molecules described by quantum chemistry methods, Molecular Mechanics or mixed QM/MM approaches. We compare the results obtained by the algorithm here presented and by the XP-PCM approach, based on a continuum description of the environment. As a test case, we study the conformational dynamics of 1,3-butadiene either as an isolated molecule ("naked" butadiene) or embedded in a cluster of argon atoms, under pressures up to 15 GPa. Overall, our results show that the XP-GAS QM/MM simulation method is in good agreement with the XP-PCM QM/Continuum model (Cammi model) in describing the effect of the pressure on static properties as the equilibrium geometry of butadiene in the ground state. Furthermore, the comparison of XP-GAS simulations with naked butadiene and butadiene in argon shows the importance, for XP-GAS and related methods, of a realistic representation of the medium in modelling pressure effects.

Van der Waals Radii of Free and Bonded Atoms from Hydrogen (Z = 1) to Oganesson (Z = 118).

Reliable numerical values of van der Waals (vdW) radii are required for constructing empirical force fields, vdW-inclusive density functional, and quantum-chemical methods, as well as for implicit solvent models. However, multiple definitions exist for vdW radii, involving either equilibrium or the closest contact distances between free or bonded atoms within molecules or crystals. For the paradigmatic case of the hydrogen atom, its reported vdW radius fluctuates between 2.15 and 3.70 Bohr depending on the definition, leading to a high uncertainty in calculations and different conceptual interpretations of noncovalent interactions. In this work, we systematically review different definitions and methodologies to establish the free and bonded vdW radii for hydrogen, based on equilibrium vdW distances in noncovalently bonded molecules, enveloping electron density cutoffs, noncovalent positron bonds in hydrogen anion dimer, vacuum virtual photon cloud caused by the hydrogen atom, and atomic dipole polarizability. By doing so, we show that the vdW radius of the free hydrogen atom is 3.16 ± 0.06 Bohr. By employing the most general and elegant definition of atomic vdW radius as a function of the atomic polarizability, we tabulate consistent values of vdW radii for all atoms in the periodic table up to Z = 118.

Thermodynamics of gaseous strontium and calcium cerates studied by Knudsen effusion mass spectrometry and estimation of relative electron ionization cross-section for CeO2(g).

Ceria-based systems are of great interest because of their unique properties. Such systems may be used as anode materials for SOFCs or in oxygen sensors. Theexploitation of these materials often requires high temperatures. In such conditions, the partial or complete evaporation of materials is possible. Therefore, knowledge of the values of partial pressures and thermodynamic properties is essential to predict and/or prevent possible consequences and accidents. Knudsen effusion mass spectrometry was used to determine partial pressures of vapor species over the SrO-CeO2 and CaO-CeO2 systems. Measurements of partial pressures were performed with a MS-1301 mass spectrometer. Vaporization was carried out using molybdenum or tungsten effusion cells. A theoretical study of gaseous strontium and calcium cerates was performed by several quantum chemical methods: density functional theory (DFT) M06, DFT PBE0, and MP2. The minimum value of relative electron ionization cross-section for CeO2 was estimated. In the temperature range of 2218-2249 K above the SrO-CeO2 system and of 2128-2208 K above the CaO-CeO2 system, gaseous M, MO, MO, CeO2, O, O2 and MCeO3 (M = Sr, Ca) were found. Energetically favorable structures of gaseous SrCeO3 and CaCeO3 were found, and vibrational frequencies were evaluated in the "rigid rotor-harmonic oscillator" approximation. On the basis of the equilibrium constants of gaseous reaction MO + CeO2 = MCeO3, the standard formation enthalpy of gaseous SrCeO3 (-913 ± 26 kJ/mol) and of gaseous CaCeO3 (-917 ± 26 kJ/mol) at 298 K were determined. The estimated value of relative electron ionization cross-section for CeO2 is in agreement with the general rule applicable for dioxides. The stability of SrCeO3 and CaCeO3 gaseous species was confirmed by KEMS. Gas-phase reactions involving gaseous SrO (CaO) and CeO2 with gaseous SrCeO3(CaCeO3) were studied. Enthalpy of formation reaction of gaseous SrCeO3(CaCeO3) from gaseous SrO (CaO) were evaluated theoretically, and the obtained value is in agreement with the experimental one.

Computationally Effective Approach for Studies of Mechanism and Thermodynamics of Heterogeneous Catalytic Processes on Metal Oxides

ABSTRACTTo understand the nature of heterogeneous catalytic processes and improve their efficiency, it is necessary to conduct both experimental and theoretical studies. At the same time, there is no unified approach to obtaining the necessary data using quantum chemistry methods. In this work, problems of the existing calculational approaches are analyzed. The obtained information is used to develop the original three‐layer embedded cluster model approach, which is shown to be the most effective. The general algorithm for obtaining such models for various oxides is formulated. The sufficient accuracy of the proposed models in predicting geometric and energy characteristics, vibrational frequencies, activation barriers, and thermodynamic characteristics is verified. The specifics of calculating the thermodynamic characteristics of heterogeneous processes using the proposed cluster models is studied in detail. The developed approach is an effective tool for studying the mechanism of heterogeneous catalytic processes both by itself and in combination with experiment.

Structure and Kinetics of Organic Acid Radicals Formed in Reaction with •OH Radicals.

The study investigated •OH-derived radicals from certain organic acids employed in nuclear fuel processing and separation using EPR spectroscopy and quantum chemistry methods. Hydroxyl radicals were generated through a Fenton-like reaction within the EPR resonator under both acidic and basic conditions, allowing for the detection of neutral and radical anions, respectively. The spectral assignment and analysis were conducted using a combination of literature data and quantum chemical calculations employing DFT theory with B3LYP or LPBE functionals and the L2a_3 basis set. The reaction of the •OH radical with lactic and glycolic acids yielded primary C-centered radicals through hydrogen abstraction from these acids. In contrast, the •OH radical exclusively generated secondary radicals from oxalic acid, whereas for citric acid, it resulted in both primary and secondary species induced by decarboxylation. The EPR spectrum of acetohydroxamic acid, upon reaction with the •OH radical, displayed a complex pattern featuring primary •N-type and N-O•-type radicals. The decay pathways of the generated radicals were primarily attributed to radical-radical reactions, with the extracted reaction rate constants generally falling within the typical range observed for such reactions. The EPR parameters calculated for potential radicals using B3LYP and LPBE functionals with L2a_3 basis set demonstrated good accuracy for neutral radicals, albeit requiring minor adjustments for radical anions.

Mechanism analysis and performance comparison of extractive distillation with different extractants for separating methyl tert-butyl ether/methanol mixture

Methyl tert-butyl ether (MTBE) is an important chemical raw material and fuel additive commonly used in the chemical industry and can form minimum boiling azeotrope with methanol (MeOH). In this article, extractive distillation processes with common-extractant (diethylene glycol, DG), counter-extractant (chlorobenzene, CB) and ionic liquid ([MIM][HSO4]) for the separation of MTBE and MeOH binary azeotrope are investigated simultaneously for the first time. Moreover, five quantum chemical calculation methods of σ-profiles, electrostatic potential, interaction energy, Hirshfeld surface and weak interaction are used simultaneously for the first time to reveal the extractive distillation mechanism of different extractants in detail. Finally, three extractive distillation processes (CED-IL, CED-DG, CED-CB) are established and optimized with total annual costs (TAC) and CO2 emissions as objectives. Their performances are comprehensively compared from the aspects of economy, environment, energy, exergy and inherent safety under the optimal conditions. The results demonstrate that, in comparison to the CED-DG process, the CED-IL process reduces TAC by 31.0%, total energy consumption (TEC) and gas emissions by 27.7%, process route index (PRI) by 36.4% and increase thermodynamic efficiency by 46.4%, while the CED-CB process increases TAC by 29.1%, TEC and gas emissions by 23.2%, PRI by 22.1% and reduce thermodynamic efficiency by 11.4%.

Effect of cytosine and adenine on graphene oxide for the sensing of dopamine: Electrochemical and theoretical studies

In this work, the electrochemical interaction between dopamine and cytosine/adenine on the graphene oxide (GO) is elucidated by quantum chemical simulation methods. Initially, GO, cytosine-functionalized GO (CGO), and adenine-functionalized GO (AGO) were prepared, and then utilized for the carbon paste electrode (CPE) modification. The as-prepared CGO-modified CPE (CGCPE) and AGO-modified CPE (AGCPE) were applied for the electrochemical sensing of dopamine (DA), and the obtained consequences illustrate the improved sensitivity of modified CPEs (MCPEs) compared to bare CPE (BCPE). Here, applied theoretical concepts based on analytical Fukui functions and frontier molecular orbitals (FMO) for cytosine and adenine functional groups in order to find active sites, and explained with evidence for the detection of DA by enhancing sensitivity. Moreover, the AGCPE has been shown higher electrochemical response compared to CGCPE for the detection of DA, and this result also agreed with theoretically from Global electron transfer results. The MCPEs were also utilized to examine the different electrochemical parameters and determine the DA concentration in the nanomolar range. The detection limit for CGCPE and AGCPE was found to be 9 nm and 7 nm, respectively.