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.

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

Previously inaccessible large S8-corona[n]arene macrocycles (n = 8-12) with alternating aryl and 1,4-C6F4 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 SNAr substitution protocols (catalytic NR4F in pyridine, R = H, Me, Bu). Macrocycle size and functionality are tunable by precursor and catalyst selection. Equivalent simple NR4F catalysis allows facile late-stage SNAr difunctionalisation of the ring C6F4 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 NR4F 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.

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.

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.

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.

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.

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.

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.

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.

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 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 (scCO2) 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 scCO2 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 scCO2. 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 scCO2.

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.

Benchmarking DFT-D Dispersion Corrections for Anharmonic Vibrational Frequencies and Harmonic Scaling Factors.

Improvements in the form of the DFT-D empirical dispersion corrections to hybrid density functional theory are shown to have made corrections sufficiently accurate to improve the calculation of both anharmonic frequencies and scaled harmonic vibrational frequencies across a wide range of commonly tested molecules. The Becke-Johnson damping function is noted as being particularly versatile across the molecules tested, and the B3LYP-D3M(BJ) and B3LYP-D3(CSO) methods are found to be the most widely applicable. Dispersion corrections are shown to be important for accurately describing carbon-hydrogen bond stretching vibrations, and standard triple-dipole based three-body terms are found to cause large errors in these anharmonic frequencies. Preliminary results also indicate that there is a cancellation of error at this level of theory when using smaller finite difference step sizes to calculate anharmonic derivatives of the nuclear potential energy surface.

A scaled CIS(D) based method for the calculation of valence and core electron ionization energies.

The calculation of electron ionization energies is a key component for the simulation of photoelectron spectroscopy. CIS(D) is a perturbative doubles correction for the single excitation configuration interaction (CIS) method which provides a new approach for computing excitation energies. It is shown that by introducing a virtual orbital subspace that consists of a single "ghost" orbital, valence electron ionization energies can be computed using a scaled CIS(D) approach with an accuracy comparable with considerably more computationally intensive methods, such as ionization-potential equation of motion coupled cluster theory, and simulated spectra show a significant improvement relative to spectra based upon Koopmans' theorem. When the model is applied to the ionization energies for core orbitals, there is an increase in the error, particularly for the heavier nuclei considered (silicon to chlorine), although the relative energy of the ionization energies are predicted accurately. In addition to its inherent computational efficiency relative to other wavefunction based approaches, a significant advantage of this approach is that the ionization energies for all electrons can be obtained in a single calculation, in contrast to Δself-consistent field based methods.

Monitoring the Formation and Reactivity of Organometallic Alkane and Fluoroalkane Complexes with Silanes and Xe Using Time-Resolved X-ray Absorption Fine Structure Spectroscopy.

Complexes with weakly coordinating ligands are often formed in chemical reactions and can play key roles in determining the reactivity, particularly in catalytic reactions. Using time-resolved X-ray absorption fine structure (XAFS) spectroscopy in combination with time-resolved IR (TRIR) spectroscopy and tungsten hexacarbonyl, W(CO)6, we are able to structurally characterize the formation of an organometallic alkane complex, determine the W-C distances, and monitor the reactivity with silane to form an organometallic silane complex. Experiments in perfluorosolvents doped with xenon afford initially the corresponding solvated complex, which is sufficiently reactive in the presence of Xe that we can then observe the coordination of Xe to the metal center, providing a unique insight into the metal-xenon bonding. These results offer a step toward elucidating the structure, bonding, and chemical reactivity of transient species by X-ray absorption spectroscopy, which has sensitivity to small structural changes. The XAFS results indicate that the bond lengths of metal-alkane (W-H-C) bond in W(CO)5(heptane) as 3.07 (±0.06) Å, which is longer than the calculated W-C (2.86 Å) for binding of the primary C-H, but shorter than the calculated W-C (3.12 Å) for the secondary C-H. A statistical average of the calculated W-C alkane bond lengths is 3.02 Å, and comparison of this value indicates that the value derived from the XAFS measurements is averaged over coordination of all C-H bonds consistent with alkane chain walking. Photolysis of W(CO)6 in the presence of HSiBu3 allows the conversion of W(CO)5(heptane) to W(CO)5(HSiBu3) with an estimated W-Si distance of 3.20 (±0.03) Å. Time-resolved TRIR and XAFS experiments following photolysis of W(CO)6 in perfluoromethylcyclohexane (PFMCH) allows the characterization of W(CO)5(PFMCH) with a W-F distance of 2.65 (±0.06) Å, and doping PFMCH with Xe allows the characterization of W(CO)5Xe with a W-Xe bond length of 3.10 (±0.02) Å.

Excited-State Vibrational Frequencies: Restricted Virtual Space Time-Dependent Density Functional Theory.

Calculating accurate vibrational frequencies for molecules with electronically excited states has an important function in many areas of photochemistry. However, calculations are often limited to smaller molecules due to the rapid growth in the degrees of freedom that must be taken into account to accurately describe larger systems. The applicability of the restricted virtual space (RVS) approximation has been studied within adiabatic linear response time-dependent density functional theory when calculating excited-state nuclear vibrational frequencies. Using the S1 and T1 electronic states of CO, CN-, HOF, H2CS, and C2H4 as representative examples, it is found that vibrational frequency calculations are particularly sensitive to this approximation, with no more than 10-20% of orbitals recommended for safe removal without a priori knowledge when using the 6-311+G(d,p) and aug-cc-pVTZ basis sets. Higher-frequency vibrations such as those with a high degree of CH bond stretching character are found to be less sensitive to the RVS than the lower-frequency vibrations, and several of the triplet states are also found to be less sensitive to this approximation than their equivalent singlet states. Occupied core orbitals and high-energy virtual orbitals with core character can also be removed without introducing significant error.

Dimers of acetic acid in helium nanodroplets.

The structural arrangement of small carboxylic acid molecules in the liquid phase remains a controversial topic. Some studies indicate a dominance of the cyclic dimer that prevails in the gas phase, whilst other studies favor short fragments of the infinite catemer chains that are found in the crystalline phase. Furthermore, difficulties in preparing and probing size-selected catemer segments have resulted in a lack of benchmark data upon which theoretical models of the condensed phases can be built. To address these issues, we have combined infrared spectroscopy and quantum chemical calculations to explore regions of the intermolecular potential energy surface associated with the formation of metastable dimer isomers. The OH stretching region of the spectrum shows that aggregation of acetic acid molecules inside liquid helium nanodroplets yields two distinct metastable dimers, whilst negligible signal is observed from the cyclic dimer that typically overwhelms this spectral region. We deduce that the most abundant isomer in superfluid helium has one strong O-HO[double bond, length as m-dash]C and one weak C-HO[double bond, length as m-dash]C hydrogen bond. Since this bonding motif is common to the dimeric repeating unit of the catemer, it is of fundamental importance for understanding intermolecular interactions in the condensed phases of small carboxylic acids.

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

It has been demonstrated that bonding and interactions within ionic liquids (ILs) can be elegantly tuned by manipulation of structure and the introduction of functional groups. Here we use XPS to investigate the impact of sulfur containing substituents on the electronic structure of a series of N-based cations, all with a common anion, [NTf2]-. The experiments reveal complexity and perturbation of delocalised systems which cannot be easily interpreted by NMR and XPS alone, DFT provides critical insight into bonding and underpins the assignment of spectra and development of deconstruction models for each system studied.