In Search of the Best Low-Cost Methods for Efficient Screening of Conformers.

Locating the lowest energy conformer is crucial for the accurate computation of equilibrium properties of molecular systems. This paper examines the performance of efficient low-cost methods in terms of the alignment and relative energies of their energy minima against the benchmark revDSD-PBEP86-D4/def2-TZVPP//MP2/cc-pVTZ potential energy surface. The low-cost methods considered include GFN-FF, GFN2-xTB, DFTB3, HF-3c, B97-3c, PBEh-3c, and r2SCAN-3c composite methods against a diverse test set of 20 compounds including alkanes, perfluoroalkyl molecules, peptides, open-shell radicals, and Zn(II) complexes of varying sizes. The "3c" composite methods are generally more accurate, but are at least 2-3 orders of magnitude more expensive than tight-binding methods which have energy minima that align well with the benchmark potential energy surface. The findings of this paper were further exploited to introduce a simple strategy involving Grimme's CENSO energy-sorting algorithm that resulted in up to an order of magnitude reduction in computational time for locating the lowest energy conformer on the revDSD-PBEP86-D4/def2-TZVPP//MP2/cc-pVTZ surface.

High-Level Quantum Chemical Prediction of C-F Bond Dissociation Energies of Perfluoroalkyl Substances.

In this study, 550 C-F bond dissociation energies (BDEs) of a variety of per- and polyfluoroalkyl substances (PFASs) obtained from high-level DLPNO-CCSD(T)/CBS calculations were used to assess the accuracy of contemporary density functional theory (DFT) and semiempirical methods. DLPNO-CCSD(T)/CBS gas phase C-F BDEs fall between 404.9-550.7 kJ mol-1 and M06-2X and ωB97M-V in conjunction with the aug-cc-pVTZ basis set predicted BDEs closest to the benchmark level with a mean absolute deviation (MAD) of 7.3 and 8.3 kJ mol-1, respectively. It was observed that DFT prediction errors increase with the degree of fluorination and system size. As such, previous model chemistry recommendations based on smaller nonfluorinated systems may not be carried over to modeling the energetics of PFASs and related systems.

Molecular Geometries and Vibrational Contributions to Reaction Thermochemistry Are Surprisingly Insensitive to the Choice of Basis Sets.

Calculation of molecular geometries and harmonic vibrational frequencies are pre-requisites for thermochemistry calculations. Contrary to conventional wisdom, this paper demonstrates that quantum chemical predictions of the thermochemistry of many gas and solution phase chemical reactions appear to be very insensitive to the choice of basis sets. For a large test set of 80 diverse organic and transition-metal-containing reactions, variations in reaction free energy based on geometries and frequencies calculated using a variety of double and triple-zeta basis sets from the Pople, Jensen, Ahlrichs, and Dunning families are typically less than 4 kJ mol-1, especially when the quasiharmonic oscillator correction is applied to mitigate the effects of low-frequency modes. Our analysis indicates that for many organic molecules and their transition states, high-level revDSD-PBEP86-D4 and DLPNO-CCSD(T)/(aug-)cc-pVTZ single-point energies usually vary by less than 2 kJ mol-1 on density functional theory geometries optimized using basis sets ranging from 6-31+G(d) to aug-pcseg-2 and aug-cc-pVTZ. In cases where these single-point energies vary significantly, indicating sensitivity of molecular geometries to the choice of basis set, there is often substantial cancellation of errors when the reaction energy or barrier is calculated. The study concludes that the choice of basis set for molecular geometry and frequencies, particularly those considered in this study, is not critical for the accuracy of thermochemistry calculations in the gas or solution phase.

Incorporation of Cation Affects the Redox Reactivity of Fe-NNN Complexes on C-H Oxidation.

Cations such as Lewis acids have been shown to enhance the catalytic activity of high-valent Fe-oxygen intermediates. Herein, we present a pyridine diamine ethylene glycol macrocycle, which can form Zn(II)- or Fe(III)-complex with the NNN site, while allowing redox-inactive cations to bind to the ethylene glycol moiety. The addition of alkali, alkali earth, and lanthanum ions resulted in positive shifts to the Fe(III/II) redox potential. Calculation of dissociation constants showed the tightest binding with a Ba2+ ion. Density functional theory calculations were used to elucidate the effects of redox inactive cations toward the electronic structures of Fe complexes. Although the Fe-NNN complexes, both in the absence and presence of cations, can catalyze C-H oxidation of 9,10-dihydroanthracene, to give anthracene [hydrogen atom transfer (HAT) product], anthrone, and anthraquinone [oxygen atom transfer (OAT) products], highest overall activity and OAT/HAT product ratios were obtained in the presence of dications, that is, Ba2+ and Mg2+, respectively.

Z-matrix template-based substitution approach for enumeration of 3D molecular structures.

The exhaustive enumeration of 3D chemical structures based on Z-matrix templates has recently been used in the quantum chemical investigation of constitutional isomers, diastereomers and rotamers. This simple yet powerful initial structure generation approach can apply beyond the investigation of compounds of identical formula by quantum chemical methods. This paper provides a comprehensive description of the overall concept followed by a practical tutorial to the approach.•The four steps required for Z-matrix template-based substitution are template construction, generation of tuples for substitution sites, removal of duplicate tuples and substitution on the template.•The generated tuples can be used to create chemical identifiers to query compound properties from chemical databases.•All of these steps are demonstrated in this paper by common model compounds and are very straightforward for an undergraduate audience to reproduce. A comparison of the approach in this paper and other options is also discussed.

Halogen bond-induced electrophilic aromatic halogenations.

In recent years, there has been increasing interest in utilising halogen bonds in organic synthesis, especially in aromatic halogenation reactions. N-Halosuccinimides and 1,3-dihalo-5,5-dimethylhydantoins are popular sources of halonium ions due to their ease of handling and low toxicities. Traditionally, these N-haloimides are activated by electrophiles, namely Brønsted and Lewis acids. The recent discovery of possible activation by nucleophilic Lewis base catalysts led to a paradigm shift in aromatic halogenation. Active functional motifs in Lewis base catalysts such as CS, R-S-R1, Ar-S-S-Ar, SO, Ar-NH2, and R2NH+Cl- form halogen bonds with the positively charged σ-hole of the halogen atoms: an essential interaction to produce halonium ions. This review highlights the evolution of the two modes of activation. Evidence of halogen bond formation from mechanistic studies of nucleophilic activation is also discussed herein.

Geometric and energetic data from ab initio calculations of haloethene, haloimine, halomethylenephosphine, haloiminophosphine, halodiazene, halodiphosphene and halocyclopropane.

This article presents theoretical data on geometric and energetic features of halogenated compounds of ethene (C[bond, double bond]C), imine (C[bond, double bond]N), methylenephosphine (C[bond, double bond]P), iminophosphine (N[bond, double bond]P), diazene (N[bond, double bond]N), diphosphene (P[bond, double bond]P) and cyclopropane (Δ). The data were obtained from ab initio geometric optimization and frequency calculations at HF, B3LYP, MP2 and CCSD levels of theory on 6-311++G(d,p) basis set. Input structures were generated by shell scripts and run by Q-Chem quantum chemical package. The output files were processed to extract geometric and energetic information by Wolfram Mathematica.

An in silico investigation of menthol metabolism.

Prevalence of mentholated products for consumption has brought great importance to studies on menthol's metabolic pathways to ensure safety, design more potent derivatives, and identify therapeutic benefits. Proposed pathways of (-)-menthol metabolism based on metabolites found experimentally in previous works by Yamaguchi, Caldwell & Farmer, Madyastha & Srivatsan and Hiki et al. were not in agreement. This in silico approach is based on the three in vivo studies and aims to resolve the discrepancies. Reactions in the pathways are conjugation with glucuronic acid/sulfate, oxidation to alcohol, aldehyde & carboxylic acid, and formation of a four-membered/five-membered ring. Gas-phase structures, standard Gibbs energies and SMD solvation energies at B3LYP/6-311++G(d,p) level were obtained for 102 compounds in the pathways. This study provides a more complete picture of menthol metabolism by combining information from three experimental studies and filling missing links in previously published pathways.