Catalysis enabled synthesis, structures, and reactivities of fluorinated S8-corona[n]arenes (n = 8-12)

Previously inaccessible large S₈-corona[n]arene macrocycles (n = 8-12) with alternating aryl and 1,4-C₆F₄ subunits are easily prepared on up to gram scales, without the need for chromatography (up to 45% yield, 10 different examples) through new high acceleration S_NAr substitution protocols (catalytic NR₄F in pyridine, R = H, Me, Bu). Macrocycle size and functionality are tunable by precursor and catalyst selection. Equivalent simple NR₄F catalysis allows facile late-stage S_NAr difunctionalisation of the ring C₆F₄ units with thiols (8 derivatives, typically 95+% yields) providing two-step access to highly functionalised fluoromacrocycle libraries. Macrocycle host binding supports fluoroaryl catalytic activation through contact ion pair binding of NR₄F and solvent inclusion. In the solid-state, solvent inclusion also intimately controls macrocycle conformation and fluorine-fluorine interactions leading to spontaneous self-assembly into infinite columns with honeycomb-like lattices.

Static Electron Correlation in Anharmonic Molecular Vibrations: A Hybrid TAO-DFT Study

Hybrid thermally-assisted-occupation density functional theory is used to examine the effects of static electron correlation on the prediction of a benchmark set of experimentally observed molecular vibrational frequencies. The B3LYP and B97-1 thermally-assisted-occupation measure of static electron correlation is important for describing the vibrations of many of the molecules that make up several popular test sets of experimental data. Shifts are seen for known multireference systems and for many molecules containing atoms from the second row of the periodic table of elements. Several molecules only show significant shifts in select vibrational modes, and significant improvements are seen for the prediction of hydrogen stretching frequencies throughout the test set.

Wavelength dependent photoextrusion and tandem photo-extrusion reactions of ninhydrin bis-acetals for the synthesis of 8-ring lactones, benzocyclobutenes and orthoanhydrides

Ninhydrin bis-acetals give access to 8-ring lactones, benzocyclo-butenes and spirocyclic orthoanhydrides through photoextrusion and tandem photoextrusion reactions. Syntheses of fimbricalyxlactone B, isoshihunine and numerous biologically-relevant heterocycles show the value of the methods, while TA-spectroscopy and TD-DFT studies provide mechanistic insights on their wavelength dependence.

Kernel Methods for Predicting Yields of Chemical Reactions

The use of machine learning methods for the prediction of reaction yield is an emerging area. We demonstrate the applicability of support vector regression (SVR) for predicting reaction yields, using combinatorial data. Molecular descriptors used in regression tasks related to chemical reactivity have often been based on time-consuming, computationally demanding quantum chemical calculations, usually density functional theory. Structure-based descriptors (molecular fingerprints and molecular graphs) are quicker and easier to calculate and are applicable to any molecule. In this study, SVR models built on structure-based descriptors were compared to models built on quantum chemical descriptors. The models were evaluated along the dimension of each reaction component in a set of Buchwald-Hartwig amination reactions. The structure-based SVR models outperformed the quantum chemical SVR models, along the dimension of each reaction component. The applicability of the models was assessed with respect to similarity to training. Prospective predictions of unseen Buchwald-Hartwig reactions are presented for synthetic assessment, to validate the generalizability of the models, with particular interest along the aryl halide dimension.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Reduced Two-Electron Interactions in Anharmonic Molecular Vibrational Calculations Involving Localized Normal Coordinates

Spatially localized vibrational normal mode coordinates are shown to reduce the importance of calculating the full set of two-electron terms in the molecular electronic Schrödinger equation. Electron correlation and dispersion interactions become less significant in (E,E)-1,3,5,7-octatetraene vibrational self-consistent field calculations when displacing remote atoms along multiple coordinates. Electron correlation interactions between spatially remote modes are also found to be less important compared to their corresponding uncorrelated interaction terms. Attenuation of the Coulomb operator indicates that the two-electron terms between remote electrons become less important for accurately describing the strongly contributing mode-coupling terms between sets of localized vibrational modes.

Influence of structure and solubility of chain transfer agents on the RAFT control of dispersion polymerisation in scCO2

Reversible addition–fragmentation chain transfer (RAFT) dispersion polymerisation of methyl methacrylate (MMA) is performed in supercritical carbon dioxide (scCO₂) with 2-(dodecylthiocarbonothioylthio)-2-methylpropionic acid (DDMAT) present as chain transfer agent (CTA) and surprisingly shows good control over PMMA molecular weight. Kinetic studies of the polymerisation in scCO₂ also confirm these data. By contrast, only poor control of MMA polymerisation is obtained in toluene solution, as would be expected for this CTA which is better suited for acrylates. In this regard, we select a range of CTAs and use them to determine the parameters that must be considered for good control in dispersion polymerisation in scCO₂. A thorough investigation of the nucleation stage during the dispersion polymerisation reveals an unexpected “in situ two-stage” mechanism that strongly determines how the CTA works. Finally, using a novel computational solvation model, we identify a correlation between polymerisation control and degree of solubility of the CTAs. All of this ultimately gives rise to a simple, elegant and counterintuitive guideline to select the best CTA for RAFT dispersion polymerisation in scCO₂.

RAFT dispersion polymerisation of methyl methacrylate is performed in scCO₂ with 2-(dodecylthiocarbonothioylthio)-2-methylpropionic acid (DDMAT) present as chain transfer agent (CTA) and surprisingly shows good control over PMMA molecular weight.

Calculating with Permanent Marker: How Blockchains Record Immutable Mistakes in Computational Chemistry

Computational science experiments within an open blockchain environment have recently been demonstrated and can improve transparency, reproducibility, and censorship resistance in theoretical scientific work. However, the append-only nature of these records also means that historical calculation errors cannot be effectively removed or changed. This process preserves otherwise unavailable data on the scientific process of error correction and is shown here for simulations of carbon monoxide.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Möbius and Hückel Cyclacenes with Dewar and Ladenburg Defects

Cyclacene nanobelts have not been synthesized in over 60 years and remain one of the last unsynthesized building blocks of carbon nanotubes. Recent work has predicted that Hückel-cyclacenes containing Dewar benzenoid ring isomers are the most stable isomeric forms for several of the smaller sizes of cyclacene belts. Here, we give a more complete picture of the isomers that are possible within these nanobelt systems by simulating embedded Ladenburg (prismane) benzenoid rings in Hückel-[n]cyclacenes (n = 5-14) and embedded Dewar benzenoid rings in twisted Möbius-[n]cyclacenes (n = 9-14). The Möbius-[9]cyclacene isomer containing one Dewar benzenoid defect and the Hückel-[5]cyclacene isomer containing two maximally spaced Ladenburg benzenoid defects are found to be more stable than their conventional Kekulé benzenoid ring counterparts. The isomers that contain Dewar and Ladenburg benzenoid rings have larger electronic singlet-triplet energy gaps and lower polyradical character when compared with the conventional isomers.

Dewar Benzenoids Discovered In Carbon Nanobelts

The synthesis of cyclacene nanobelts remains an elusive goal dating back over 60 years. These molecules represent the last unsynthesized building block of carbon nanotubes and may be useful both as seed molecules for the preparation of structurally well-defined carbon nanotubes and for understanding the behavior and formation of zigzag nanotubes more broadly. Here we report the discovery that isomers containing two Dewar benzenoid rings are the preferred form for several sizes of cyclacene. The predicted lower polyradical character and higher singlet-triplet stability that these isomers possess compared with their pure benzenoid counterparts suggest that they may be more stable synthetic targets than the structures that have previously been identified. Our findings should facilitate the exploration of new routes to cyclacene synthesis through Dewar benzene chemistry.

Computational chemistry experiments performed directly on a blockchain virtual computer

Blockchain technology has had a substantial impact across multiple disciplines, creating new methods for storing and processing data with improved transparency, immutability, and reproducibility. These developments come at a time when the reproducibility of many scientific findings has been called into question, including computational studies. Here we present a computational chemistry simulation run directly on a blockchain virtual machine, using a harmonic potential to model the vibration of carbon monoxide. The results demonstrate for the first time that computational science calculations are feasible entirely within a blockchain environment and that they can be used to increase transparency and accessibility across the computational sciences.

Influence of molecular design on radical spin multiplicity: characterisation of BODIPY dyad and triad radical anions

A strategy to create organic molecules with high degrees of radical spin multiplicity is reported in which molecular design is correlated with the behaviour of radical anions in a series of BODIPY dyads. Upon reduction of each BODIPY moiety radical anions are formed which are shown to have different spin multiplicities by electron paramagnetic resonance (EPR) spectroscopy and distinct profiles in their cyclic voltammograms and UV-visible spectra. The relationship between structure and multiplicity is demonstrated showing that the balance between singlet, biradical or triplet states in the dyads depends on relative orientation and connectivity of the BODIPY groups. The strategy is applied to the synthesis of a BODIPY triad which adopts an unusual quartet state upon reduction to its radical trianion.

Dimers of acetic acid in helium nanodroplets

Two metastable dimers are created inside superfluid helium and studied using infrared spectroscopy to provide insight into condensed phase structures.

The effect of coordination of alkanes, Xe and CO2 (η1-OCO) on changes in spin state and reactivity in organometallic chemistry: a combined experimental and theoretical study of the photochemistry of CpMn(CO)3

A combined experimental and theoretical study is presented of several ligand addition reactions of the triplet fragment 3CpMn(CO)2 formed upon photolysis of CpMn(CO)3. Experimental data are provided for reactions in n-heptane and perfluoromethylcyclohexane (PFMCH), as well as in PFMCH doped with C2H6, Xe and CO2. In PFMCH we find that the conversion of 3CpMn(CO)2 to 1CpMn(CO)2(PFMCH) is much slower (τ = 18 (±3) ns) than the corresponding reactions in conventional alkanes (τ = 111 (±10) ps). We measure the effect of the coordination ability by doping PFMCH with alkane, Xe and CO2; these doped ligands form the corresponding singlet adducts with significantly variable formation rates. The reactivity as measured by the addition timescale follows the order 1CpMn(CO)2(C5H10) (τ = 270 (±10) ps) > 1CpMn(CO)2Xe (τ = 3.9 (±0.4) ns) ∼ 1CpMn(CO)2(CO2) (τ = 4.7 (±0.5) ns) > 1CpMn(CO)2(C7F14) (τ = 18 (±3) ns). Electronic structure theory calculations of the singlet and triplet potential energy surfaces and of their intersections, together with non-adiabatic statistical rate theory, reproduce the observed rates semi-quantitatively. It is shown that triplet adducts of the ligand and 3CpMn(CO)2 play a role in the kinetics, and account for the variable timescales observed experimentally.

Reduced Basis Set Dependence in Anharmonic Frequency Calculations Involving Localized Coordinates

Localized normal coordinates are known to be effective in speeding up anharmonic frequency calculations by reducing the complexity of the nuclear Hamiltonian and wave function. Displacing atoms in localized coordinates can also cause relatively small changes in the electronic structure, which can be exploited for further computational efficiency improvements during ab initio electronic structure calculations of the potential energy surface by reducing the electronic basis set dependence. Three different schemes for reducing the basis set dependence have been investigated in this work. These include combining localized coordinate schemes with general mixed basis sets, distance based force-field reductions, and using coordinate specific basis sets. The importance of accurately describing electronic interactions is found to diminish both for multicoordinate terms involving the displacement of remote atoms and when describing the interactions between more remote atoms within specific coordinates.

The impact of sulfur functionalisation on nitrogen-based ionic liquid cations

XPS is used to investigate the impact of sulfur containing substituents on the electronic structure of a series of N-based cations, all with a common anion, [NTf₂]⁻. The experimental data is complex and cannot be easily deconstructed, DFT provides critical insight into bonding and electronic structure for each system studied.

Can aliphatic anchoring groups be utilised with dyes for p-type dye sensitized solar cells?

A series of novel laterally anchoring tetrahydroquinoline derivatives have been synthesized and investigated for their use in NiO-based p-type dye-sensitized solar cells. The kinetics of charge injection and recombination at the NiO-dye interface for these dyes have been thoroughly investigated using picosecond transient absorption and time-resolved infrared measurements. It was revealed that despite the anchoring unit being electronically decoupled from the dye structure, charge injection occurred on a sub picosecond timescale. However, rapid recombination was also observed due to the close proximity of the electron acceptor on the dyes to the NiO surface, ultimately limiting the performance of the p-DSCs.

Intermediate vibrational coordinate localization with harmonic coupling constraints

Optimized normal coordinates can significantly improve the speed and accuracy of vibrational frequency calculations. However, over-localization can occur when using unconstrained spatial localization techniques. The unintuitive mixtures of stretching and bending coordinates that result can make interpreting spectra more difficult and also cause artificial increases in mode-coupling during anharmonic calculations. Combining spatial localization with a constraint on the coupling between modes can be used to generate coordinates with properties in-between the normal and fully localized schemes. These modes preserve the diagonal nature of the mass-weighted Hessian matrix to within a specified tolerance and are found to prevent contamination between the stretching and bending vibrations of the molecules studied without a priori classification of the different types of vibration present. Relaxing the constraint can also be used to identify which normal modes form specific groups of localized modes. The new coordinates are found to center on more spatially delocalized functional groups than their fully localized counterparts and can be used to tune the degree of vibrational correlation energy during anharmonic calculations.

Examining the impact of harmonic correlation on vibrational frequencies calculated in localized coordinates

Carefully choosing a set of optimized coordinates for performing vibrational frequency calculations can significantly reduce the anharmonic correlation energy from the self-consistent field treatment of molecular vibrations. However, moving away from normal coordinates also introduces an additional source of correlation energy arising from mode-coupling at the harmonic level. The impact of this new component of the vibrational energy is examined for a range of molecules, and a method is proposed for correcting the resulting self-consistent field frequencies by adding the full coupling energy from connected pairs of harmonic and pseudoharmonic modes, termed vibrational self-consistent field (harmonic correlation). This approach is found to lift the vibrational degeneracies arising from coordinate optimization and provides better agreement with experimental and benchmark frequencies than uncorrected vibrational self-consistent field theory without relying on traditional correlated methods.

Calculating excited state properties using Kohn-Sham density functional theory

The accuracy of excited states calculated with Kohn-Sham density functional theory using the maximum overlap method has been assessed for the calculation of adiabatic excitation energies, excited state structures, and excited state harmonic and anharmonic vibrational frequencies for open-shell singlet excited states. The computed Kohn-Sham adiabatic excitation energies are improved significantly by post self-consistent field spin-purification, but remain too low compared with experiment with a larger error than time-dependent density functional theory. Excited state structures and vibrational frequencies are also improved by spin-purification. The structures show a comparable accuracy to time-dependent density functional theory, while the harmonic vibrational frequencies are found to be more accurate for the majority of vibrational modes. The computed harmonic vibrational frequencies are also further improved by perturbative anharmonic corrections, suggesting a good description of the potential energy surface. Overall, excited state Kohn-Sham density functional theory is shown to provide an efficient method for the calculation of excited state structures and vibrational frequencies in open-shell singlet systems and provides a promising technique that can be applied to study large systems.