At what chain length do unbranched alkanes prefer folded conformations?

The entropic penalty for associative reactions and their physical treatment during routine computations

A systematic study of the entropic penalty for associative reactions is presented. It is shown that computed solution-phase Gibbs free energies typically overestimate entropic contributions. This entropic penalty for associative reactions in solution, i.e., if the number of particles decreases along the reaction coordinate (sum of stoichiometric numbers ), originates from the insufficient treatment of entropic effects by implicit solvent models. We propose an additive correction scheme to Gibbs free energies that is suitable for routine applications by non-expert users. This correction is based on Garza's formalism for the solution-phase entropy [A. J. Garza, J. Chem. Theory Comput., 2019, 15, 3204.] that is physically sound and embedded into an efficient black-box type algorithm. To critically evaluate the entropic penalty and its proposed treatment, we compiled an experimental benchmark set of 31 Δ_rG and 22 in 15 different solvents. Using a representative best-practice computational protocol (at wave function theory (WFT) based DLPNO-CCSD(T) and density functional theory (DFT) based revDSD-PBEP86-D4 level with an implicit solvent model), we determined a sizeable entropic penalty ranging from 2-11 kcal mol^-1. Using the correction scheme presented herein, the entropic penalty is corrected to the chemical accuracy of ≤1 kcal mol^-1 (WFT and DFT). The same applies to at the WFT level. Barriers at the DFT level are overestimated by 2 kcal mol^-1 (classic) and underestimated by 2 kcal mol^-1 (corrected). This effect is attributed to the finding that barriers computed at the DFT level are systematically 2-3 kcal mol^-1 lower than barriers obtained with WFT.

Performance of Localized-Orbital Coupled-Cluster Approaches for the Conformational Energies of Longer n-Alkane Chains

We report an update and enhancement of the ACONFL (conformer energies of large alkanes [J. Phys. Chem. A2022,126, 3521-3535]) dataset. For the ACONF12 (n-dodecane) subset, we report basis set limit canonical coupled-cluster with singles, doubles, and perturbative triples [i.e., CCSD(T)] reference data obtained from the MP2-F12/cc-pV{T,Q}Z-F12 extrapolation, [CCSD(F12*)-MP2-F12]/aug-cc-pVTZ-F12, and a (T) correction from conventional CCSD(T)/aug-cc-pV{D,T}Z calculations. Then, we explored the performance of a variety of single and composite localized-orbital CCSD(T) approximations, ultimately finding an affordable localized natural orbital CCSD(T) [LNO-CCSD(T)]-based post-MP2 correction that agrees to 0.006 kcal/mol mean absolute deviation with the revised canonical reference data. In tandem with canonical MP2-F12 complete basis set extrapolation, this was then used to re-evaluate the ACONF16 and ACONF20 subsets for n-hexadecane and n-icosane, respectively. Combining those with the revised canonical reference data for the dodecane conformers (i.e., ACONF12 subset), a revised ACONFL set was obtained. It was then used to assess the performance of different localized-orbital coupled-cluster approaches, such as pair natural orbital localized CCSD(T) [PNO-LCCSD(T)] as implemented in MOLPRO, DLPNO-CCSD(T₀) and DLPNO-CCSD(T₁) as implemented in ORCA, and LNO-CCSD(T) as implemented in MRCC, at their respective "Normal", "Tight", "vTight", and "vvTight" accuracy settings. For a given accuracy threshold and basis set, DLPNO-CCSD(T₁) and DLPNO-CCSD(T₀) perform comparably. With "VeryTightPNO" cutoffs, explicitly correlated DLPNO-CCSD(T₁)-F12/VDZ-F12 is the best pick among all the DLPNO-based methods tested. To isolate basis set incompleteness from localized-orbital-related truncation errors (domain, LNOs), we have also compared the localized coupled-cluster approaches with canonical DF-CCSD(T)/aug-cc-pVTZ for the ACONF12 set. We found that gradually tightening the cutoffs improves the performance of LNO-CCSD(T), and using a composite scheme such as vTight + 0.50[vTight - Tight] improves things further. For DLPNO-CCSD(T₁), "TightPNO" and "VeryTightPNO" offer a statistically similar accuracy, which gets slightly better when T_CutPNO is extrapolated to the complete PNO space limit. Similar to Brauer et al.'s [Phys. Chem. Chem. Phys.2016,18 (31), 20905-20925] previous report for the S66x8 noncovalent interactions, the dispersion-corrected direct random phase approximation (dRPA)-based double hybrids perform remarkably well for the ACONFL set. While the revised reference data do not affect any conclusions on the less accurate methods, they may upend orderings for more accurate methods with error statistics on the same order as the difference between reference datasets.

Refinement of the Optimized Potentials for Liquid Simulations Force Field for Thermodynamics and Dynamics of Liquid Alkanes

Torsion and Lennard-Jones parameters of the optimized potentials for liquid simulations (OPLS) all-atom force field have been refined for describing thermodynamics and dynamics of a wide range of liquid alkanes. Monte Carlo statistical mechanics (MC) and molecular dynamics (MD) simulations were carried out. For thermodynamics properties, MC simulations with truncated electrostatic interactions performed very closely to MD simulations with a Verlet neighbor list and the particle mesh Ewald algorithm. The average errors in comparison with experimental data for computed properties were improved with the modified force field (OPLS/2020), especially for long-chain alkanes. For liquid densities, heats of vaporization, and free energies of hydration, the average errors are 0.01 g/cm³, 0.2 kcal/mol, and ca. 0.5 kcal/mol, respectively; significant gains were made for relative heats of vaporization of isomeric series. Results for self-diffusion coefficients also reproduce experimental data well for linear alkane liquids up to hexadecane. The new force field is suitable for use in improved modeling of myriad systems of importance in chemistry, biology, and materials science.

Conformational Energy Benchmark for Longer n-Alkane Chains

We present the first benchmark set focusing on the conformational energies of highly flexible, long n-alkane chains, termed ACONFL. Unbranched alkanes are ubiquitous building blocks in nature, so the goal is to be able to calculate their properties most accurately to improve the modeling of, e.g., complex (biological) systems. Very accurate DLPNO-CCSD(T1)/CBS reference values are provided, which allow for a statistical meaningful evaluation of even the best available density functional methods. The performance of established and modern (dispersion corrected) density functionals is comprehensively assessed. The recently introduced r²SCAN-V functional shows excellent performance, similar to efficient composite DFT methods like B97-3c and r²SCAN-3c, which provide an even better cost-accuracy ratio, while almost reaching the accuracy of much more computationally demanding hybrid or double hybrid functionals with large QZ AO basis sets. In addition, we investigated the performance of common wave function methods, where MP2/CBS surprisingly performs worse compared to the simple D4 dispersion corrected Hartree-Fock. Furthermore, we investigate the performance of several semiempirical and force field methods, which are commonly used for the generation of conformational ensembles in multilevel workflows or in large scale molecular dynamics studies. Outstanding performance is obtained by the recently introduced general force field, GFN-FF, while other commonly applied methods like the universal force field yield large errors. We recommend the ACONFL as a helpful benchmark set for parametrization of new semiempirical or force field methods and machine learning potentials as well as a meaningful validation set for newly developed DFT or dispersion methods.

Toward Chemical Accuracy in Predicting Enthalpies of Formation with General-Purpose Data-Driven Methods

Enthalpies of formation and reaction are important thermodynamic properties that have a crucial impact on the outcome of chemical transformations. Here we implement the calculation of enthalpies of formation with a general-purpose ANI-1ccx neural network atomistic potential. We demonstrate on a wide range of benchmark sets that both ANI-1ccx and our other general-purpose data-driven method AIQM1 approach the coveted chemical accuracy of 1 kcal/mol with the speed of semiempirical quantum mechanical methods (AIQM1) or faster (ANI-1ccx). It is remarkably achieved without specifically training the machine learning parts of ANI-1ccx or AIQM1 on formation enthalpies. Importantly, we show that these data-driven methods provide statistical means for uncertainty quantification of their predictions, which we use to detect and eliminate outliers and revise reference experimental data. Uncertainty quantification may also help in the systematic improvement of such data-driven methods.

Behavior of Long-Chain Hydrocarbons at High Pressures and Temperatures

Although long-chain aliphatic hydrocarbons are documented in meteorites, their origin is poorly understood. A key question is whether they are pristine or a byproduct of terrestrial alteration? To understand if these long-chain hydrocarbons are indigenous, it will be important to explore their thermodynamic and mechanical stability at conditions experienced by extraterrestrial objects during atmospheric entry and passage. Extreme pressures and temperatures experienced by meteorites are likely to alter the molecular organization of these long-chain hydrocarbons. These structural changes associated with extreme conditions are often documented via laboratory-based Raman spectroscopic measurements. So far, Raman spectroscopic measurements have investigated the effect of static compression on the aliphatic hydrocarbons. The effect of temperature on the structural changes remains poorly explored. To bridge this gap, in this study, we have explored the behavior of two aliphatic hydrocarbons at simultaneously high pressures and temperatures. We have used a resistively heated diamond anvil cell. On compression to moderate pressures, the appearance of new vibrational modes in the low-energy region confirms prior studies and is related to the bending of the linear chains. Upon heating to ∼220 °C, we note that the new low-energy mode undergoes softening. The mode softening is likely due to the combination of unbending of the alkane chain and mode anharmonicity.

Three-Component Self-Assembly Changes its Course: A Leap from Simple Polymers to 3D Networks of Spherical Host-Guest Assemblies

One‐pot self‐assembly reactions of the polyphosphorus complex [Cp*Fe(η⁵‐P₅)] (A), a coinage metal salt AgSbF₆, and flexible aliphatic dinitriles NC(CH₂)_xCN (x=1–10) yield 1D, 2D, and 3D coordination polymers. The seven‐membered backbone of the dinitrile was experimentally found as the borderline for the self‐assembly system furnishing products of different kinds. At x<7, various rather simple polymers are exclusively formed possessing either 0D or 1D Ag/A structural motifs connected by dinitrile spacers, while at x≥7, the self‐assembly switches to unprecedented extraordinary 3D networks of nano‐sized host–guest assemblies (SbF₆)@[(A)₉Ag₁₁]¹¹⁺ (x=7) or (A)@[(A)₁₂Ag₁₂]¹²⁺ (x=8–10) linked by dinitriles. The polycationic nodes represent the first superspheres based on A and silver and are host–guest able. All products are characterized by NMR spectroscopy, mass spectrometry, and single‐crystal X‐ray diffraction. The assemblies [(A)₁₂Ag₁₂]¹²⁺ were visualized by transmission electron microscopy.

Calculation of absolute molecular entropies and heat capacities made simple

We propose a fully-automated composite scheme for the accurate and numerically stable calculation of molecular entropies by efficiently combining density-functional theory (DFT), semi-empirical methods (SQM), and force-field (FF) approximations. The scheme is systematically expandable and can be integrated seamlessly with continuum-solvation models. Anharmonic effects are included through the modified rigid-rotor-harmonic-oscillator (msRRHO) approximation and the Gibbs-Shannon formula for extensive conformer ensembles (CEs), which are generated by a metadynamics search algorithm and are extrapolated to completeness. For the first time, variations of the ro-vibrational entropy over the CE are consistently accounted-for through a Boltzmann-population average. Extensive tests of the protocol with the two standard DFT approaches B97-3c and B3LYP-D3 reveal an unprecedented accuracy with mean deviations <1 cal mol-1 K-1 (about <1-2%) for the total gas phase molecular entropy of medium-sized molecules. Even for the hardship case of extremely flexible linear alkanes (C14H30-C16H34), errors are only about 3 cal mol-1 K-1. Comprehensive tests indicate a relatively strong variation of the conformational entropy on the underlying level of theory for typical drug molecules, inferring the complex potential energy surfaces as the main source of error. Furthermore, we show some application examples for the calculation of free energy differences in typical chemical reactions.

Low-Frequency (Gigahertz to Terahertz) Depolarized Raman Scattering Off n-Alkanes, Cycloalkanes, and Six-Membered Rings: A Physical Interpretation

Molecular liquids have long been known to undergo various distinct intermolecular motions, from fast librations and cage-rattling oscillations to slow orientational and translational diffusion. However, their resultant gigahertz to terahertz spectra are far from simple, appearing as broad shapeless bands that span many orders of magnitude of frequency, making meaningful interpretation troublesome. Ad hoc spectral line shape fitting has become a notoriously fine art in the field; a unified approach to handling such spectra is long overdue. Here we apply ultrafast optical Kerr-effect (OKE) spectroscopy to study the intermolecular dynamics of room-temperature n-alkanes, cycloalkanes, and six-carbon rings, as well as liquid methane and propane. This work provides stress tests and converges upon an experimentally robust model across simple molecular series and range of temperatures, providing a blueprint for the interpretation of the dynamics of van der Waals liquids. This will enable the interpretation of low-frequency spectra of more complex liquids.

Pyrolysis of Alkanes: A Computational Approach

In this study, we investigate the kinetic and thermodynamic aspects of thermal cracking reactions of short paraffin chains by density functional theory (DFT) methods. The thermal cracking reactions have been modeled for a series of shorter unbranched alkanes at 673 K by following a free-radical mechanism. Benchmark calculations have been carried out with different functionals (B3LYP, M06-2X, PBE0, BMK, B3PW91) and basis sets (6-31G(d,p), 6-311+G(d,p)) to determine the most suitable DFT method, and the results were compared to the available experimental data. Computations were also performed at the CBS-QB3 level to evaluate the accuracy of the DFT method. The thermodynamic and kinetic properties of the initiation, hydrogen atom transfer (HAT), and decomposition (β-scission) reactions are intensely discussed to better understand the trends in product distributions at high temperatures. Evans-Polanyi (EP) relations have been used to build a linear relationship between the enthalpy of reactions and their activation energies; this process may be useful for the determination of the kinetic parameters of longer paraffin chains as well. Finally, the preexponential factors of short-chain paraffin have been calculated and classified based on the identity of the radicalic products. The latter, together with the activation energies derived from the EP relations can be used safely for the prediction of the rate constants for long paraffin chains.

Semiempirical Quantum-Chemical Methods with Orthogonalization and Dispersion Corrections

We present two new semiempirical quantum-chemical methods with orthogonalization and dispersion corrections: ODM2 and ODM3 (ODM x). They employ the same electronic structure model as the OM2 and OM3 (OM x) methods, respectively. In addition, they include Grimme's dispersion correction D3 with Becke-Johnson damping and three-body corrections E ABC for Axilrod-Teller-Muto dispersion interactions as integral parts. Heats of formation are determined by adding explicitly computed zero-point vibrational energy and thermal corrections, in contrast to standard MNDO-type and OM x methods. We report ODM x parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine that are optimized with regard to a wide range of carefully chosen state-of-the-art reference data. Extensive benchmarks show that the ODM x methods generally perform better than the available MNDO-type and OM x methods for ground-state and excited-state properties, while they describe noncovalent interactions with similar accuracy as OM x methods with a posteriori dispersion corrections.

Gas-Phase Heat of Formation Values for Buckminsterfullerene (C60), C70 Fullerene (C70), Corannulene, Coronene, Sumanene, and Other Polycyclic Aromatic Hydrocarbons Calculated Using Density Functional Theory (M06 2X) Coupled with a Versatile Inexpensive Group-Equivalent Approach

A straightforward procedure using density functional theory (M06 2X) coupled with a group-equivalent approach is described that was used to calculate gas-phase heat of formation (Δ_f H°_g,298) values for buckminsterfullerene (C₆₀), C70 fullerene (C₇₀), corannulene, coronene, and sumanene. This procedure was also used to calculate exceptionally accurate Δ_f H°_g,298 values for a variety of single-ring aromatic and 2-7 ring polycyclic aromatic hydrocarbons (PAHs) as well as a large selection of other hydrocarbons and phenols. The approach described herein is internally consistent, and results for C₆₀, C₇₀, corannulene, coronene, and sumanene are in very close agreement with results reported by others who used higher-level computational theory. Statistical analysis of a test set containing benzene and 18 two to seven ring PAHs demonstrated that by using this approach a mean absolute deviation (MAD) and a root-mean-square deviation (RMSD) of 0.8 and 1.3 kJ/mol, respectively, were achieved for reference/experimental Δ_f H°_g,298 values versus calculated/predicted Δ_f H°_g,298 values. For statistical analysis of a larger test set containing 235 aromatic and aliphatic hydrocarbons and phenols, a MAD and a RMSD of 1.2 and 1.9 kJ/mol, respectively, were achieved for reference/experimental Δ_f H°_g,298 values versus calculated/predicted Δ_f H°_g,298 values.

Size extensivity of elastic properties of alkane fragments

Using MP2, CCSD, and B3LYP methods of computational chemistry, we show length dependence in the intrinsic elastic properties of short alkane fragments. For isolated alkane fragments of finite length in the gas phase and zero temperature, the intrinsic elasticity constants are found to vary with the number of carbon atoms and its parity. From extrapolation of the elasticity constants calculations to infinite chain length, and by comparing with in-situ elasticity constant of single poly(ethylene) molecule obtained with atomic force microscopy, we estimate the softening effect of environment on the extension response of the polymer.

Validation of Molecular Simulation: An Overview of Issues

Computer simulation of molecular systems enables structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic, or supra-molecular level. To interpret results of such simulations appropriately, the quality of the calculated properties must be evaluated. This depends on the way the simulations are performed and on the way they are validated by comparison to values Q^exp of experimentally observable quantities Q. One must consider 1) the accuracy of Q^exp , 2) the accuracy of the function Q(r^N ) used to calculate a Q-value based on a molecular configuration r^N of N particles, 3) the sensitivity of the function Q(r^N ) to the configuration r^N , 4) the relative time scales of the simulation and experiment, 5) the degree to which the calculated and experimental properties are equivalent, and 6) the degree to which the system simulated matches the experimental conditions. Experimental data is limited in scope and generally corresponds to averages over both time and space. A critical analysis of the various factors influencing the apparent degree of (dis)agreement between simulations and experiment is presented and illustrated using examples from the literature. What can be done to enhance the validation of molecular simulation is also discussed.

Benchmarking Quantum Chemical Methods: Are We Heading in the Right Direction?

Theoreticians and experimentalists should work together more closely to establish reliable rankings and benchmarks for quantum chemical methods. Comparison to carefully designed experimental benchmark data should be a priority. Guidelines to improve the situation for experiments and calculations are proposed.

Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Benchmarks for Ground-State Properties

The semiempirical orthogonalization-corrected OMx methods (OM1, OM2, and OM3) go beyond the standard MNDO model by including additional interactions in the electronic structure calculation. When augmented with empirical dispersion corrections, the resulting OMx-Dn approaches offer a fast and robust treatment of noncovalent interactions. Here we evaluate the performance of the OMx and OMx-Dn methods for a variety of ground-state properties using a large and diverse collection of benchmark sets from the literature, with a total of 13035 original and derived reference data. Extensive comparisons are made with the results from established semiempirical methods (MNDO, AM1, PM3, PM6, and PM7) that also use the NDDO (neglect of diatomic differential overlap) integral approximation. Statistical evaluations show that the OMx and OMx-Dn methods outperform the other methods for most of the benchmark sets.

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.

Stretching and folding of 2-nanometer hydrocarbon rods

Linear alkanes CnH2n+2 in vacuum isolation are finite models for an infinite polyethylene chain. Using spontaneous Raman scattering in supersonic jet expansions for n = 13-21 in different spectral ranges, we determine the minimal chain length nh for the cohesion-driven folding of the preferred extended all-trans conformation into a hairpin structure. We treat fully stretched all-trans alkanes as molecular "nanorods" and derive Young's modulus E for the stretching of an isolated single-strand polyethylene fibre by extrapolating the longitudinal acoustic mode to infinite chain length. Two key quality parameters for accurate intra- and intermolecular force fields of hydrocarbons (nh = 18 ± 1, E = 305 ± 5 GPa) are thus derived with high accuracy from experimental spectroscopy.

A concise review of computational studies of the carbon dioxide–epoxide copolymerization reactions

The production of polycarbonates from carbon dioxide and epoxides is an important route by which waste CO₂ is converted into useful products. This review surveys the use of computational chemistry toward understanding this reaction.