Search

The S66x8 noncovalent interactions benchmark has been re-evaluated at the "sterling silver" level, using explicitly correlated MP2-F12 near the complete basis set limit, CCSD(F12*)/aug-cc-pVTZ-F12, and a (T) correction from conventional CCSD(T)/sano-V{D,T}Z+ calculations. The revised reference value disagrees by 0.1 kcal/mol RMS with the original Hobza benchmark and its revision by Brauer et al, but by only 0.04 kcal/mol variety from the "bronze" level data in Kesharwani et al., Aust. J. Chem. 71, 238-248 (2018). We then used these to assess the performance of localized-orbital coupled cluster approaches with and without counterpoise corrections, such as PNO-LCCSD(T) as implemented in MOLPRO, DLPNO-CCSD (T1) as implemented in ORCA, and LNO-CCSD(T) as implemented in MRCC, for their respective "Normal", "Tight", and "very Tight" settings. We also considered composite approaches combining different basis sets and cutoffs. Furthermore, in order to isolate basis set convergence from domain truncation error, for the aug-cc-pVTZ basis set we compared PNO, DLPNO, and LNO approaches with canonical CCSD(T). We conclude that LNO-CCSD(T) with veryTight criteria performs very well for "raw" (CP-uncorrected), but struggles to reproduce counterpoise-corrected numbers even for veryVeryTight criteria: this means that accurate results can be obtained using either extrapolation from basis sets large enough to quench basis set superposition error (BSSE) such as aug-cc-pV{Q,5}Z, or using a composite scheme such as Tight{T,Q}+1.11[vvTight(T) - Tight(T)]. In contrast, PNO-LCCSD(T) works best with counterpoise, while performance with and without counterpoise is comparable for DLPNO-CCSD(T1). Among more economical methods, the highest accuracies are seen for dRPA75-D3BJ, {\omega}B97M-V, {\omega}B97M(2), revDSD-PBEP86-D4, and DFT(SAPT) with a TDEXX or ATDEXX kernel., Comment: Final published version with CC license

The computational modeling of fullerenes plays a fundamental role in designing low-dimension carbon nanostructures. Nevertheless, the relative energies of fullerenes larger than C

The structure-activity relationships (SAR) of 8-phenyl-1,3-dipropylxanthine derivatives in binding to recombinant human A

The relative free energies of the isomers formed upon N-chlorination of each nitrogen atom within the DNA nucleobases (adenine, guanine, and thymine) have been obtained using the high-level G4(MP2) composite ab initio method (the free energies of the N-chlorinated isomers of cytosine have been reported at the same level of theory previously). Having identified the lowest energy N-chlorinated derivatives for each nucleobase, we have computed the free energies associated with chlorine transfer from N-chlorinated nucleobases to other unsubstituted bases. Our results provide quantitative support pertaining to the results of previous experimental studies, which demonstrated that rapid chlorine transfer occurs from N-chlorothymidine to cytidine or adenosine. The results of our calculations in the gas-phase reveal that chlorine transfer from N-chlorothymine to either cytosine, adenine, or guanine proceed via exergonic processes with ΔG o values of −50.3 (cytosine), −28.0 (guanine), and −6.7 (adenine) kJ mol–1. Additionally, we consider the effect of aqueous solvation by augmenting our gas-phase G4(MP2) energies with solvation corrections obtained using the conductor-like polarizable continuum model. In aqueous solution, we obtain the following G4(MP2) free energies associated with chlorine transfer from N-chlorothymine to the three other nucleobases: −58.4 (cytosine), −26.4 (adenine), and −18.7 (guanine) kJ mol–1. Therefore, our calculations, whether in the gas phase or in aqueous solution, clearly indicate that chlorine transfer from any of the N-chlorinated nucleobases to cytosine provides a thermodynamic sink for the active chlorine. This thermodynamic preference for chlorine transfer to cytidine may be particularly deleterious since previous experimental studies have shown that nitrogen-centered radical formation (via N–Cl bond homolysis) is more easily achieved in N-chlorinated cytidine than in other N-chlorinated nucleosides.

We have examined the use of systematic bond-separation reactions and purposely constructed chemistry-preserving isodesmic reactions for the thermochemical calculation of aromatic hydrocarbon species. The bond-separation approach yields somewhat disappointing accuracy even when the reaction energies are obtained with generally robust composite and double-hybrid (DH) density functional theory (DFT) methods. In contrast, for the purposely constructed reactions, we find a dramatic improvement in the accuracy for energies calculated with all methods examined. Notably, for medium-sized aromatic hydrocarbons, we find that an effective approach for formulating a well-balanced reaction is to split the target species into two halves with an aromatic overlapping region. Overall, the G4(MP2)-XK, MPW2PLYP, MN15, PBE, and DC-DFTB3 methods are reasonable within their respective classes of methods for the calculation of bond-separation as well as chemistry-preserving isodesmic reactions. We have further computed per-carbon atomization energy (AE) for a series of D6h benzene-type molecules, and thus obtained a formula for extrapolation to the graphene limit [AEn = 711.5 × (1 − 1/n0.640) kJ mol−1, where n = number of carbons]. It suggests that nano-graphene with a length larger than 10 nm would resemble properties of bulk graphene, and conversely, downsizing a nano-graphene beyond this point may lead to considerably altered properties from the bulk.

Perylene bisimide (PBI) cyclophanes have been shown to form ideal hosts for large aromatic molecules, offering rigid cavities and ideal monomer spacing for π–π stacking interactions between the host and guest. Here, we show computationally that these characteristics also make PBI cyclophanes ideal receptors for the planar transition structures of a number of representative stereoinversion processes. The observed catalytic effect can hereby be explained using two fundamental concepts: transition state stabilization due to shape-complementarity between cyclophane and the transition structures and, to a smaller extent, ground state destabilization due to non-complementarity between cyclophane and equilibrium structures. Energy decomposition analyses of the noncovalent interactions between the host and guest structures reveal that dispersion and electrostatic interactions make the largest contributions to stabilizations and thus suggest that catalysis in this system is mainly driven by π–π stacking interactions.

The 21-homododecahedryl cation is a unique system in terms of its complete fluxionality based on two different rearrangements. In this work, we report the quantum tunneling effects that drive the reactions at temperatures where the semi-classical kinetics are impossible. We postulate that the tunnel effect in this system can serve to create a refrigerator that may operate at arbitrarily low temperatures.

Ethynylcyclopropenylidene (

The isomerisation energies and thermal stabilities of the entire set of C40 fullerene isomers are investigated using molecular dynamics simulations with the recently developed machine-learning-based Gaussian Approximation Potential (GAP-20) force field. Our simulations predict high statistical correlation between the relative isomerisation energy distribution of all C40 isomers and their thermal fragmentation temperatures. The most energetically and thermally stable isomer is the symmetric C40-D2 isomer and its fragmentation temperature is 3875 ± 25 K. In contrast, the least stable isomer is the D5d nanorod-shaped isomer which is 146.8 kcal mol−1 higher in energy and decomposes at 2975 ± 25 K, i.e., around 900 K below the most stable isomer. In addition, we show that most isomers lie in the range between 20 and 60 kcal mol−1 above the most stable isomer, and nearly 60% of isomers decompose in the temperature range of 3425–3675 K. These computational results provide valuable insights into the synthesis of C40 fullerenes at high-temperatures.

In youth ice hockey, girls are reported to suffer more concussions than boys, peaking around 13-14 years old, which may be related to differences in the level of brain trauma experienced by the players. The purpose of this research was to describe the differences in brain trauma characteristics, specifically the magnitude and frequency of head impacts between Peewee boys and girls from playing ice hockey. Thirty games of Peewee boys and Peewee girl's ice hockey were recorded to document the head impact events. These events were reconstructed using physical and computational techniques to estimate the strain to the brain tissue. The results found that Peewee boys experienced more head impacts than girls, specifically from the shoulder, ice, boards, and fist/punches (p 0.05). The boys also experienced more medium strain category impacts (p 0.05). These results identify that Peewee boys and girls engage in ice hockey differently, which affects the risk of brain trauma likely to be encountered while during game play, suggesting that the increased rate of concussion for girls may not be related to impact magnitudes within the sport but increased reporting of symptoms of concussion or gender differences in brain tissue response to an impact.

The need for sustainable and large-scale energy supply has led to significant development of renewable energy and energy storage technologies. Divalent metal ion (Mg, Ca, and Zn) batteries are promising energy storage technologies for the sustainable energy future, but the need for suitable electrode materials have limited their commercial development. This paper investigates, at the atomic scale, the adsorption and migration of Mg, Ca, and Zn on pristine and defective graphene surfaces, to bring insight into the metal storage and mobility in graphene and carbon-based anodes for divalent metal ion batteries. Such atomistic studies can help address the challenges facing the development of novel divalent metal battery technologies, and to understand the storage differences between divalent and monovalent metal-ion batteries. The adsorption of Ca on the graphene-based system is shown to be more energetically favorable than the adsorption of both Mg and Zn, with Ca showing adsorption behavior similar to the monovalent ions (Li, Na, and K). This was further investigated in terms of metal migration on the graphene surface, with much higher migration energy barriers for Ca than for Mg and Zn on the graphene systems, leading to the trapping of Ca at defect sites to a larger extent.

Gas-sensing properties of nitrogenated holey graphene (C2N), graphdiyne (GDY) and their van der Waals heterostructure (C2N…GDY) have been studied towards particular volatile organic compounds (VOCs) by means of spin-polarized, dispersion-corrected DFT calculations. We find that VOCs such as acetone, ethanol, propanal, and toluene interact weakly with the GDY monolayer; however, the bindings are significantly enhanced with the C2N monolayer and the hybrid C2N…GDY heterostructure in AB stacking. Electron localization function (ELF) analysis shows that all VOCs are van der Waals bound (physical binding) to the 2D materials, which result in significant changes of the charge density of C2N and GDY monolayers and the C2N…GDY heterostructure. These changes alter the electronic properties of C2N and GDY, and the C2N…GDY heterostructure, upon VOC adsorption, which are investigated by density-of-states plots. We further apply thermodynamic analysis to study the sensing characteristics of VOCs under varied conditions of pressure and temperature. Our findings clearly indicate that the C2N…GDY heterostructure is a promising material for sensing of certain VOCs.

Inverse vulcanization provides dynamic and responsive materials made from elemental sulfur and unsaturated cross-linkers. These polymers have been used in a variety of applications such as energy storage, infrared optics, repairable materials, environmental remediation, and precision fertilizers. In spite of these advances, there is a need for methods to recycle and reprocess these polymers. In this study, polymers prepared by inverse vulcanization are shown to undergo reactive compression molding. In this process, the reactive interfaces of sulfur polymers are brought into contact by mechanical compression. Upon heating these molds at relatively low temperatures (≈100 °C), chemical bonding occurs at the polymer interfaces by S-S metathesis. This method of processing is distinct from previous studies on inverse vulcanization because the polymers examined in this study do not form a liquid phase when heated. Neither compression nor heating alone was sufficient to mold these polymers into new architectures, so this is a new concept in the manipulation of sulfur polymers. Additionally, high-level ab initio calculations revealed that the weakest S-S bond in organic polysulfides decreases linearly in strength from a sulfur rank of 2 to 4, but then remains constant at about 100 kJ mol-1 for higher sulfur rank. This is critical information in engineering these polymers for S-S metathesis. Guided by this insight, polymer repair, recycling, and repurposing into new composites was demonstrated.

The purpose of this research was to compare the frequency and magnitude of head impact events between Pee Wee and Bantam ice hockey players. Videos of Pee Wee and Bantam boys' ice hockey were analysed to determine the frequency and type of head impact events. The head impact events were then reconstructed in the laboratory using physical and finite element models to determine the magnitude of strain in the brain tissues. The results showed that Pee Wee boys experienced more head impacts from elbows and boards, while Bantam players had more head impacts to the glass. Pee Wee and Bantam players experienced similar frequency and magnitudes of very low, low, and medium and above (med+) levels of strain to the brain. This research suggests to ice hockey leagues and coaches that to reduce the incidence of these levels of brain trauma, consideration must be given to either reducing the level of contact along the boards or the removal of body checking. In addition, companies who innovate in ice hockey should develop protective devices and equipment strategies that aim to reduce the risk of head injury from shoulder and glass impacts for Bantam players.

Polymers made by inverse vulcanization can be assembled, repaired, and recycled at room temperature through nucleophile-catalyzed S–S metathesis., Inverse vulcanization is a copolymerization of elemental sulfur and alkenes that provides unique materials with high sulfur content (typically ≥50% sulfur by mass). These polymers contain a dynamic and reactive polysulfide network that creates many opportunities for processing, assembly, and repair that are not possible with traditional plastics, rubbers and thermosets. In this study, we demonstrate that two surfaces of these sulfur polymers can be chemically joined at room temperature through a phosphine or amine-catalyzed exchange of the S–S bonds in the polymer. When the nucleophile is pyridine or triethylamine, we show that S–S metathesis only occurs at room temperature for a sulfur rank > 2—an important discovery for the design of polymers made by inverse vulcanization. This mechanistic understanding of the S–S metathesis was further supported with small molecule crossover experiments in addition to computational studies. Applications of this chemistry in latent adhesives, additive manufacturing, polymer repair, and recycling are also presented.

Two-dimensional nanoporous graphene (NPG) with uniformly distributed nanopores has been synthesized recently and shown remarkable electronic, mechanical, thermal, and optical properties with potent...

Eleven isomers of SiC4H2 lying within 50 kcal mol-1 have been theoretically investigated using density functional theory and high-level coupled-cluster methods. Among them, four isomers, 1-ethynyl-3-silacycloprop-1(2)-en-3-ylidene (1), diethynylsilylidene (2), 1-sila-1,2,3,4-pentatetraenylidene (4), and 1,3-butadiynylsilylidene (5), have already been identified in the laboratory. The current investigation reports three low-lying ( 1 eV) silylidenes [2-sila-(didehydrovinylidene)cyclopropene (8), an isomer with a planar tetracoordinate carbon (ptC) atom (10), and 1-ethynyl-1-silapropadienylidene (11)], which remain elusive in the laboratory to date. Isomer 9 also contains a ptC atom, which turned out to be a transition state at all levels. Though all isomers are polar (μ ≠ 0), rotational spectrum is available only for 4. Using matrix isolation, three isomers (1, 2, and 5) have been trapped in the laboratory at 10 K. Considering the astrochemical relevance of silicon-carbide clusters in the interstellar medium, the current theoretical data demand new molecular spectroscopic studies on SiC4H2. Surprisingly, unlike the isovalent C5H2 isomers, where the bent carbenes are yet to be identified in the laboratory, the bent silylidenes (2 and 5) have been trapped in the case of SiC4H2. In both the cases, molecules with transannular C-C and/or Si-C bonds remain elusive, though they lie in the low-lying region. Using suitable precursors, whether these peculiar geometries (especially 3 and 6) would be identified or not in the laboratory needs to be addressed by molecular spectroscopists. The present investigation documents structural and spectroscopic information of SiC4H2 isomers, which may compliment future molecular spectroscopic observations including radioastronomical searches.

The racemization of axially chiral biaryls is a fundamental step toward transforming kinetic resolutions into dynamic kinetic resolutions (DKRs). The high enantiomerization barriers of many biaryl compounds of synthetic relevance, however, may render DKR strategies challenging. Here, we computationally explore the potential of a