A Database of Formation Enthalpies of Nitrogen Species by Compound Methods (CBS-QB3, CBS-APNO, G3, G4)

Modeling the thermochemistry of nitrogen-containing compounds via group additivity

First-principles based kinetic modeling is essential to gain insight into the governing chemistry of nitrogen-containing compounds over a wide range of technologically important processes, e.g. pyrolysis, oxidation and combustion. It also enables the development of predictive, fundamental models key to improving understanding of the influence of nitrogen-containing compounds present as impurities or process additives, considering safety, operability and quality of the product streams. A prerequisite for the generation of detailed fundamental kinetic models is the availability of accurate thermodynamic properties. To address the scarcity of thermodynamic properties for nitrogen-containing compounds, a consistent set of 91 group additive values and three non-nearest-neighbor interactions has been determined from a dataset of CBS-QB3 calculations for 300 species, including 104 radicals. This dataset contains a wide range of nitrogen-containing functionalities, i.e. imine, nitrile, nitro, nitroso, nitrite, nitrate and azo functional groups. The group additivity model enables the approximation of the standard enthalpy of formation and standard entropy at 298 K as well as the standard heat capacities over a large temperature range, i.e. 300-1500 K. For a test set of 27 nitrogen-containing compounds, the group additivity model succeeds in approximating the ab initio calculated values for the standard enthalpy of formation with a MAD of 2.3 kJ mol-1. The MAD for the standard entropy and heat capacity is lower than 4 and 2 J mol-1 K-1, respectively. For a test set of 11 nitrogen-containing compounds, the MAD between experimental and group additivity approximated values for the standard enthalpy of formation amounts to 2.8 kJ mol-1.

Isolated and solvated chiroptical behavior in conformationally flexible butanamines

The nonresonant optical activity of two highly flexible aliphatic amines, (2R)-3-methyl-2-butanamine (R-MBA) and (2R)-(3,3)-dimethyl-2-butanamine (R-DMBA), has been probed under isolated and solvated conditions to examine the roles of conformational isomerism and to explore the influence of extrinsic perturbations. The optical rotatory dispersion (ORD) measured in six solvents presented uniformly negative rotatory powers over the 320-590 nm region, with the long-wavelength magnitude of chiroptical response growing nearly monotonically as the dielectric constant of the surroundings diminished. The intrinsic specific optical rotation,α λ T (in deg dm-1 [g/mL]-1 ), extracted for ambient vapor-phase samples of R-MBA [-11.031(98) and -2.29 (11)] and R-DMBA [-9.434 (72) and -1.350 (48)] at 355 and 633 nm were best reproduced by counterintuitive solvents of high polarity (yet low polarizability) like acetonitrile and methanol. Attempts to interpret observed spectral signatures quantitatively relied on the linear-response frameworks of density-functional theory (B3LYP, cam-B3LYP, and dispersion-corrected analogs) and coupled-cluster theory (CCSD), with variants of the polarizable continuum model (PCM) deployed to account for the effects of implicit solvation. Building on the identification of several low-lying equilibrium geometries (nine for R-MBA and three for R-DMBA), ensemble-averaged ORD profiles were calculated at T = 300 K by means of the independent-conformer ansatz, which enabled response properties predicted for the optimized structure of each isomer to be combined through Boltzmann-weighted population fractions derived from corresponding relative internal-energy or free-energy values, the latter of which stemmed from composite CBS-APNO and G4 analyses. Although reasonable accord between theory and experiment was realized for the isolated (vapor-phase) species, the solution-phase results were less satisfactory and tended to degrade progressively as the solvent polarity increased. These trends were attributed to solvent-mediated changes in structural parameters and energy metrics for the transition states that separate and putatively isolate the equilibrium conformations supported by the ground electronic potential-energy surface, with the resulting displacement of barrier locations and/or decrease of barrier heights compromising the underlying premise of the independent-conformer ansatz.

C2H5NO Isomers: From Acetamide to 1,2-Oxazetidine and Beyond

This work documents the properties of a number of isomers of molecular formula C₂H₅NO from the most stable, acetamide, through 1,2-oxazetidine and including even higher energy species largely of a dipolar nature. Only two of the isomers have been detected in emissions from the interstellar medium (ISM); possible further candidates are identified, and the likelihood of their being detectable is considered. In general, hardly any of these compounds have been discussed in the existing chemical literature, so this work represents an important contribution extending the canon of chemical bonding which can contribute to machine learning, providing a more exacting test of AI applications. The presence in the ISM of acetamide, CH₃C(O)NH₂, is the subject of current debate with no clear and obvious paths to its formation; it is shown that a 1,3-[H]-transfer from (E,Z)-ethanimidic acid, CH₃C(OH)=NH, is feasible in spite of an energy barrier of 130 kJ mol^–1. It is speculated that imidic acid can itself be formed from abundant precursors, H₂O and CH₃C≡N, in an acid-induced, water addition, autocatalytic reaction on water–ice grains. H₃CC≡NH₃CC≡NH⁺ + H₂OH₃CC(O⁺H₂)=NHH₃CC(OH)=NH + H₃O⁺

Pyrolysis and Combustion Chemistry of Pyrrole, a Reference Component for Bio-oil Surrogates: Jet-Stirred Reactor Experiments and Kinetic Modeling

Fast-pyrolysis bio-oils (FPBOs) obtained from lignocellulosic biomass are gaining attention as sustainable fuels for various applications, including the transport sector and power production. A significant fraction of bio-oils is constituted by nitrogen-containing compounds (N fuels) that should be considered when developing surrogate models for FPBOs. Moreover, the content of N fuels in FPBOs is expected to strongly contribute to the production of nitrogen oxides (NO _x ) directly from fuel-bound nitrogen (fuel NO _x ), in addition to the thermal NO _x formation pathways typical of high-temperature combustion conditions. This work investigates the pyrolysis and combustion chemistry of pyrrole (C₄H₅N), a candidate reference fuel component for FPBO surrogate models. Speciation measurements in an atmospheric pressure jet-stirred reactor have been performed for both pyrolysis and oxidation conditions. Pyrolysis experiments have been performed for 1% pyrrole/helium mixtures over the temperature range T = 925-1200 K. Oxidation experiments were carried out for 1% pyrrole/oxygen/helium mixtures at three equivalence ratios (φ = 0.5, 1.0, and 2.0) over the temperature range T = 700-1200 K. These new data significantly extend the number of experimental targets for kinetic model validation available at present for pyrrole combustion. After a thorough revision of previous theoretical and kinetic modeling studies, a preliminary kinetic model is developed and validated by means of comparison to new experimental data and those previously reported in the literature. The rate of production and sensitivity analyses highlight important pathways deserving further investigations for a better understanding of pyrrole and, more in general, N fuel combustion chemistry. A critical discussion on experimental challenges to be faced when dealing with pyrrole is also reported, encouraging further experimental investigation with advanced diagnostics.

Density functional theory for the thermodynamic gas-phase investigation of butanol biofuel and its isomers mixed with gasoline and ethanol

Herein, we present the results of our study on the thermodynamic properties of the isomers of butanol (n-butanol, 2-butanol, i-butanol, and t-butanol) to evaluate their thermodynamic potential as a complementary biofuel and/or substitute for ethanol and gasoline. The Gaussian09W software was used to perform molecular geometry optimization calculations using density functional theory with the B3lyp hybrid function using the base set 6-311++g(d,p) and the compound methods G3, G4, and CBS-QB3. Calculations of the fundamental frequency of the molecules were performed to obtain the molecular vibration modes for the respective frequencies. These calculations provided thermodynamic parameters such as the entropy, enthalpy, and specific molar heat at constant pressure, all as a function of the temperature. The parameter values obtained by each method were compared to the experimental values available in the literature. The results showed good accuracy, especially those obtained at the B3lyp/6-311++g(d,p) level for n-butanol. The error between the theoretical and experimental values for the combustion enthalpy of n-butanol was less than 4% at 298.15 K; due to the good prediction of its thermodynamic properties, we used n-butanol as a model for the prediction of other thermodynamic properties. We started a molecular docking study of four ligands, namely, n-butanol, ethanol, propanol, heptane, isooctane, and methanol interacting with butanol isomers. The highest values of affinity energy found were for N-butanol. The possible formation of hydrogen bonds, associations by means of London forces, hydrogen, and alkyl interactions were analyzed. n-Butanol was added to ethanol-gasoline mixtures in the temperature range of 298.15 to 600 K and the results suggest that n-butanol has a higher calorific value than gasoline-ethanol mixtures in G30E, G40E, G50E, G60E, G70E, G80E, G90E, and E100 blends. As such, n-butanol releases greater amounts of heat during combustion and is thus a viable alternative to biofuels.

Energy-Geometry Dependency of Molecular Structures: A Multistep Machine Learning Approach

There is growing interest in estimating quantum observables while circumventing expensive computational overhead for facile in silico materials screening. Machine learning (ML) methods are implemented to perform such calculations in shorter times. Here, we introduce a multistep method based on machine learning algorithms to estimate total energy on the basis of spatial coordinates and charges for various chemical structures, including organic molecules, inorganic molecules, and ions. This method quickly calculates total energy with 0.76 au in root-mean-square error (RMSE) and 1.5% in mean absolute percent error (MAPE) when tested on a database of optimized and unoptimized structures. Using similar molecular representations, experimental thermochemical properties were estimated, with MAPE as low as 6% and RMSE of 8 cal/mol·K for heat capacity in a 10-fold cross-validation.

mHDFS-HoF: A generalized multilevel homodesmotic fragment-separation reaction based program for heat-of-formation calculation for acyclic hydrocarbons

Based on our modified classification of elemental species, a framework for automatic generation of multilevel Homodesmotic fragment-separation (mHDFS) reactions for chemical species was proposed. Combined the mHDFS framework with a database of heat of formation (HoF) and the calculated electronic structure data for the elemental mHD species, the mHDFS-HoF program was constructed in C/C++ language to calculate heat of formation for a species of interest on-the-fly. Using the electronic structure data calculated at CBS-QB3 level of theory for the elemental mHD species, applications and robustness of the code were discussed with several acyclic hydrocarbon systems including neutral and radical species. On-going work and extension to other systems were also discussed. The program and the supporting files can be freely downloaded at https://sites.google.com/view/mhdfs/. © 2019 Wiley Periodicals, Inc.

Ab initio dynamics of hydrogen abstraction from N2H4 by OH radicals: an RRKM-based master equation study

The detailed kinetic mechanism of the N₂H₄ + OH reaction is comprehensively reported for a wide condition range of conditions (i.e., 200–3000 K & 1–7600 Torr) using the CCSD(T)/CBS//M06-2X/6-311++G(3df,2p) level and the RRKM-based master equation rate model.

The Alexandria library, a quantum-chemical database of molecular properties for force field development

Data quality as well as library size are crucial issues for force field development. In order to predict molecular properties in a large chemical space, the foundation to build force fields on needs to encompass a large variety of chemical compounds. The tabulated molecular physicochemical properties also need to be accurate. Due to the limited transparency in data used for development of existing force fields it is hard to establish data quality and reusability is low. This paper presents the Alexandria library as an open and freely accessible database of optimized molecular geometries, frequencies, electrostatic moments up to the hexadecupole, electrostatic potential, polarizabilities, and thermochemistry, obtained from quantum chemistry calculations for 2704 compounds. Values are tabulated and where available compared to experimental data. This library can assist systematic development and training of empirical force fields for a broad range of molecules.

Systematic Search for Chemical Reactions in Gas Phase Contributing to Methanol Formation in Interstellar Space

A massive search for chemical routes leading to methanol formation in gas phase has been conducted using computational chemistry, at the CBS-QB3 level of theory. The calculations were performed at five different temperatures (100, 80, 50, 20, and 10 K) and at three pressures (0.1, 0.01, and 0.001 atm) for each temperature. The search was focused on identifying reactions with the necessary features to be viable in the interstellar medium (ISM). A searching strategy was applied to that purpose, which allowed to reduce an initial set of 678 possible reactions to a subset of 11 chemical routes that are recommended, for the first time, as potential candidates for contributing to methanol formation in the gas phase of the ISM. They are all barrier-less, and thus they are expected to take place at collision rates. Hopefully, including these reactions in the currently available models, for the gas-phase methanol formation in the ISM, would help improving the predicted fractional abundance of this molecule in dark clouds. Further investigations, especially those dealing with grain chemistry and electronic excited states, would be crucial to get a complete picture of the methanol formation in the ISM.

Thermodynamic DFT analysis of natural gas

Density functional theory was performed for thermodynamic predictions on natural gas, whose B3LYP/6-311++G(d,p), B3LYP/6-31+G(d), CBS-QB3, G3, and G4 methods were applied. Additionally, we carried out thermodynamic predictions using G3/G4 averaged. The calculations were performed for each major component of seven kinds of natural gas and to their respective air + natural gas mixtures at a thermal equilibrium between room temperature and the initial temperature of a combustion chamber during the injection stage. The following thermodynamic properties were obtained: internal energy, enthalpy, Gibbs free energy and entropy, which enabled us to investigate the thermal resistance of fuels. Also, we estimated an important parameter, namely, the specific heat ratio of each natural gas; this allowed us to compare the results with the empirical functions of these parameters, where the B3LYP/6-311++G(d,p) and G3/G4 methods showed better agreements. In addition, relevant information on the thermal and mechanic resistance of natural gases were investigated, as well as the standard thermodynamic properties for the combustion of natural gas. Thus, we show that density functional theory can be useful for predicting the thermodynamic properties of natural gas, enabling the production of more efficient compositions for the investigated fuels. Graphical abstract Investigation of the thermodynamic properties of natural gas through the canonical ensemble model and the density functional theory.

Efficient DLPNO–CCSD(T)-Based Estimation of Formation Enthalpies for C-, H-, O-, and N-Containing Closed-Shell Compounds Validated Against Critically Evaluated Experimental Data

An accurate and cost-efficient methodology for the estimation of the enthalpies of formation for closed-shell compounds composed of C, H, O, and N atoms is presented and validated against critically evaluated experimental data. The computational efficiency is achieved through the use of the resolution-of-identity (RI) and domain-based local pair-natural orbital coupled cluster (DLPNO-CCSD(T)) approximations, which results in a drastic reduction in both the computational cost and the number of necessary steps for a composite quantum chemical method. The expanded uncertainty for the proposed methodology evaluated using a data set of 45 thoroughly vetted experimental values for molecules containing up to 12 heavy atoms is about 3 kJ·mol-1, competitive with those of typical calorimetric measurements. For the compounds within the stated scope, the methodology is shown to be superior to a representative, more general, and widely used composite quantum chemical method, G4.

Modeling Nitrogen Species as Pollutants: Thermochemical Influences

To simulate emissions of nitrogen-containing compounds in practical combustion environments, it is necessary to have accurate values for their thermochemical parameters, as well as accurate kinetic values to describe the rates of their formation and decomposition. Significant disparity is observed in the literature for the former, and we therefore present herein high-accuracy ab initio gas-phase thermochemistry for 60 nitrogenous compounds, many of which are important in the formation and consumption chemistry of NOx species. Several quantum-chemical composite methods (CBS-APNO, G3, and G4) were utilized to derive enthalpies of formation via the atomization method. Entropies and heat capacities were calculated from traditional statistical thermodynamics, with oscillators treated as anharmonic based on ro-vibrational property analyses carried out at the B3LYP/cc-pVTZ level of theory. The use of quantum chemical methods, along with the treatments of anharmonicities and hindered rotors, ensures accurate enthalpy of formation, entropy, and heat capacity values across the temperature range 298.15-3000 K. The implications of these results for atmospheric and combustion modeling are discussed.