Modeling the Conformational Preference of the Lignocellulose Interface and Its Interaction with Weak Acids

Reliable Quantum-Chemistry Heats of Formation for an Extensive Set of C-, H-, N-, O-, F-, S-, Cl-, Br-Containing Molecules in the NIST Chemistry Webbook

In the present study, we have computed the heat of formation (HOF) for over 500 C-, H-, N-, O-, F-, S-, Cl-, Br-containing molecules in the NIST Chemistry Webbook with a previously established methodology [from the highest- to lowest-level methods, W1X-2, CCSD(T)-F12b, DSD-PBEP86, and ωB97M-V, with the lower levels calibrated against higher levels for the atomic energies, see: J. Phys. Chem. A 2022, 126, 4981-4990]. We find a reasonable level of agreement between the computed and NIST values for the present set of species. However, the set of F-containing compounds shows considerably larger discrepancies, which can in part be attributed to dubious experimental values, as we have demonstrated in some cases. With our highest-level computed HOFs, we validated the lower-level methods used in our protocol. Specifically, CCSD(T)-F12b yields chemically accurate (±4.2 kJ mol-1) values for all types of molecules, while DSD-PBEP86 and ωB97M-V yield similar levels of accuracy for most systems, with key exceptions being molecules with numerous electron-withdrawing F and NO2 groups. Our results further support the use of the protocol for the computation of HOFs, particularly for systems with few reliable reference values.

State-of-the-art local correlation methods enable affordable gold standard quantum chemistry for up to hundreds of atoms

In this feature, we review the current capabilities of local electron correlation methods up to the coupled cluster model with single, double, and perturbative triple excitations [CCSD(T)], which is a gold standard in quantum chemistry. The main computational aspects of the local method types are assessed from the perspective of applications, but the focus is kept on how to achieve chemical accuracy (i.e., <1 kcal mol-1 uncertainty), as well as on the broad scope of chemical problems made accessible. The performance of state-of-the-art methods is also compared, including the most employed DLPNO and, in particular, our local natural orbital (LNO) CCSD(T) approach. The high accuracy and efficiency of the LNO method makes chemically accurate CCSD(T) computations accessible for molecules of hundreds of atoms with resources affordable to a broad computational community (days on a single CPU and 10-100 GB of memory). Recent developments in LNO-CCSD(T) enable systematic convergence and robust error estimates even for systems of complicated electronic structure or larger size (up to 1000 atoms). The predictive power of current local CCSD(T) methods, usually at about 1-2 order of magnitude higher cost than hybrid density functional theory (DFT), has become outstanding on the palette of computational chemistry applicable for molecules of practical interest. We also review more than 50 LNO-based and other advanced local-CCSD(T) applications for realistic, large systems across molecular interactions as well as main group, transition metal, bio-, and surface chemistry. The examples show that properly executed local-CCSD(T) can contribute to binding, reaction equilibrium, rate constants, etc. which are able to match measurements within the error estimates. These applications demonstrate that modern, open-access, and broadly affordable local methods, such as LNO-CCSD(T), already enable predictive computations and atomistic insight for complicated, real-life molecular processes in realistic environments.

Accurate Computation of Aqueous pKas of Biologically Relevant Organic Acids: Overcoming the Challenges Posed by Multiple Conformers, Tautomeric Equilibria, and Disparate Functional Groups with the Fully Black-Box pK-Yay Method

The use of static electronic structure calculations to compute solution-phase pKas offers a great advantage in that a macroscopic bulk property could be computed via microscopic computations involving very few molecules. There are various sources of errors in the quantum chemical calculations though. Overcoming these errors to accurately compute pKas of a plethora of acids is an active area of research in physical chemistry pursued by both computational as well as experimental chemists. We recently developed the pK-Yay method in our attempt to accurately compute aqueous pKas of strong and weak acids. The method is fully black-box, computationally inexpensive, and is very easy for even a nonexpert to use. However, the method was thus far tested on very few molecules (only 16 in all). Herein, in order to assess the future applicability of pK-Yay, we study the effect of multiple conformers, the presence of tautomers under equilibrium, and the impact of a wide variety of functional groups (derivatives of acetic acid with substituents at various positions, dicarboxylic acids, aromatic carboxylic acids, amines and amides, phenols and thiols, and fluorine bearing organic acids). Starting with more than 1000 conformers and tautomers, this study establishes that overall errors of ∼ 1.0 pKa units are routinely obtained for a majority of the molecules. Larger errors are noted in cases where multiple charges, intramolecular hydrogen bonding, and several ionizable functional groups are simultaneously present. An important conclusion to emerge from this work is that, the computed pKas are insensitive (difference <0.5) to whether we consider multiple conformers/tautomers or only choose the most stable conformer/tautomer. Further, pK-Yay captures the stereoelectronic effects arising due to differing axial vs equatorial pattern, and is useful to predict the dominant acid-base equilibrium in a system featuring several equilibria. Overall, pK-Yay may be employed in several chemical applications featuring organic molecules and biomonomers.

Bond dissociation energies of ethyl valerate and tripropionin

CONTEXT

Due to the expected decrease in the availability of conventional oils, numerous studies are currently underway to find complementary sources of energy. Among the explored avenue is that of biofuels. Ethyl valerate (ETV) and tripropionin (TPP) are two biofuels whose thermal decomposition has not received the attention it deserves. Herein, we have evaluated the bond dissociation enthalpies (BDHs) to predict how easy it is to break some bonds in these compounds, and subsequently contribute to revealing the initiation step in their combustion reactions. Our computations consistently predict C4-C5 and C1-C2 bonds in ETV and TPP as the weakest bonds, likely to break first and initiate the thermal decomposition of these two compounds, respectively. The conformational changes in ETV and TPP have only a small influence on the BDHs of 1 kcal/mol at M06-2X/6-311 + G(3df,2p). B3LYP and ωB97XD appear to be the most affordable methods for estimating BDHs at 6-31G(d,p) as they give good results for ETV (RMSD: 2.94 kcal/mol and 3.22 kcal/mol) and performed better than CBS-QB3 (RMSD: 3.64 kcal/mol). Using a larger basis set, the M06-2X (RMSD: 3.61 kcal/mol) and ωB97XD (RMSD: 3.51 kcal/mol) functionals are found to provide the most accurate predictions at 6-311 + G(3df,2p) as compared to G4MP2.

METHODS

BDHs of ETV and TPP are computed using density functional theory (DFT) and quantum chemistry composite methods at 6-31G(d,p) and 6-311 + G(3df,2p) levels. Because of its reliability and accuracy in thermochemical calculations, the G4MP2 theory is used as a reference to gauge the performance of DFT methods. All the calculations were carried out using the Gaussian 09 program.

Compilation of Ionic Clusters with the Rock Salt Structure: Accurate Benchmark Thermochemical Data, Assessment of Quantum Chemistry Methods, and the Convergence Behavior of Lattice Energies

In the present study, computational quantum chemistry is used to obtain lattice energies (LEs) for a range of ionic clusters with the NaCl structure. Specifically, the compounds include NaF, NaCl, MgO, MgS, KF, CaO, and CaS clusters, (MX)_n, with n = 1, 2, 4, 6, 8, 12, 16, 24, 32, 40, 50, 60, 75, 90, and 108. The highest-level W2 and W1X-2 methods are applied to the small clusters with n = 1 to 8 (the MX35 data set). The assessment with MX35 shows that, for the calculation of geometries and vibrational frequencies, the PBE0-D3(BJ) and PBE-D3(BJ) DFT methods are reasonable, but the calculation of atomization energies is more challenging. This is a result of different systematic deviations for clusters of different species. Thus, species-specific adjustments are applied for larger clusters, which are calculated with the DuT-D3 double-hybrid DFT method, the MN15 DFT method, and the PM7 semi-empirical method. They yield smoothly converging LEs to the bulk values. It is also found that, for the alkali-metal species, the LEs for a single molecule are ∼70% of the bulk values, while for the alkali-earth species, they are ∼80%. This has enabled a straightforward means to the first-principles estimation of LEs for similarly structured ionic compounds.

Advances in machine learning for high value-added applications of lignocellulosic biomass

Lignocellulose can be converted into biofuel or functional materials to achieve high value-added utilization. Biomass utilization process is complex and multi-dimensional. This paper focuses on the biomass conversion reaction conditions, the preparation of biomass-based functional materials, the combination of biomass conversion and traditional wet chemistry, molecular simulation and process simulation. This paper analyzes the mechanism, advantages and disadvantages of important machine learning (ML) methods. The application examples of ML in different aspects of high value utilization of lignocellulose are summarized in detail. The challenges and future prospects of ML in this field are analyzed.

Divide-and-Conquer Linear-Scaling Quantum Chemical Computations

Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.

High-Level Quantum Chemistry Reference Heats of Formation for a Large Set of C, H, N, and O Species in the NIST Chemistry Webbook and the Identification and Validation of Reliable Protocols for Their Rapid Computation

A recent study has examined the accuracy of NIST heats of formation for a set of C, H, and O-containing species with a proposed low-cost quantum chemistry approach. In the present study, we have used high-level methods such as W1X-2 to obtain these data more rigorously, which we have then used to assess the NIST and the previously computed values. We find that many of these NIST data that are as suggested to be unreliable by the previous study are indeed inconsistent with our high-level reference values. However, we also find substantial deviations for the previously computed values from our benchmark. Thus, we have assessed the performance of alternative low-cost methods. In our assessment, we have additionally examined C, H, N, and O-containing species for which heats of formation are available from the NIST database. We find the ωB97M-V/ma-def2-TZVP, DSD-PBEP86/ma-def2-TZVP, and CCSD(T)-F12b/aug'-cc-pVDZ methods to be adequate for obtaining heats of formation with the atomization approach, once their atomic energies are optimized with our benchmark. Notably, the low-cost ωB97M-V method yields values that agree to be within 10 kJ mol^-1 for more than 90% of the (∼1500) species. A higher 20 kJ mol^-1 threshold captures 98% of the data. The outlier species typically contain many electron-withdrawing (nitro) groups. In these cases, the use of isodesmic-type reactions rather than the atomization approach is more reliable. Our assessment has also identified significant outliers from the NIST database, for which experimental re-determination of the heats of formation would be desirable.