Search

A thorough DFT investigation of over 800 small fullerenes reveals the ones with promising thermochemical properties for electronics applications.

In the present study, we have investigated the performance of RIJCOSX DLPNO-CCSD(T)-F12 methods for a wide range of systems. Calculations with a high-accuracy option ["DefGrid3 RIJCOSX DLPNO-CCSD(T

Density functional theory (DFT) is currently experiencing a golden age. The past two decades witnessed remarkable advances in the general applicability of density functionals in the top rungs of Jacob’s Ladder. Nevertheless, Jacob’s Ladder may have reached its highest rung in terms of dependencies on occupied (rung four) and unoccupied orbitals (rung five). Moreover, the fifth rung is associated with a computational cost far greater than the lower rungs. Another limitation is that each rung includes dozens of different functionals, and at present, there is no clear pathway for systematic improvements within each rung of the ladder. This highlight provides an overview of the exchange–correlation (XC) hole and how it could be used in developing new density functionals. We begin with a brief overview of the current status and challenges in developing better density functionals, followed by the intimate relationship between the XC functional and hole. We present a conceptually simple and computationally economical method for calculating the XC hole and how this method could offer new directions in developing better exchange–correlation functionals.

Directing groups play an important role in controlling the selectivity of C-H bond activation. Here we demonstrate that for iron, the nature of the directing group (

Metal–ligand cooperation is an important aspect in earth-abundant metal catalysis. Utilizing ligands as electron reservoirs to supplement the redox chemistry of the metal has resulted in many new exciting discoveries. Here, we demonstrate that iron bipyridine-diimine (BDI) complexes exhibit an extensive electron-transfer series that spans a total of five oxidation states, ranging from the trication [Fe(BDI)]3+ to the monoanion [Fe(BDI]−1. Structural characterization by X-ray crystallography revealed the multifaceted redox noninnocence of the BDI ligand, while spectroscopic (e.g., 57Fe Mössbauer and EPR spectroscopy) and computational studies were employed to elucidate the electronic structure of the isolated complexes, which are further discussed in this report., In this study, the multifaceted electrochemistry of a bis(imino)bipyridine iron(II) complex is described. Spectroscopic studies demonstrate that throughout the various reduction processes, the oxidation state on the metal center remains constant (Fe(II); S = 1), whereas the ligand accesses four distinct oxidation states: [BDI]0 ⇄ [BDI]3−. The observed redox noninnocence of the ligand was investigated by computational methods, which is further discussed in this report.

We evaluate the performance of ten common carbon force fields for the interaction energies in double and triple layered carbon onions. In particular, we consider the C20@C60, C20@C80, C20@C180, C80@C240, C60@C240 and C240@C540 double-layer carbon onions and C60@C240@C540 and C80@C240@C540 triple-layered carbon onions. We consider the following carbon force fields: Tersoff, REBO-II, AIREBO, AIREBO-M, screened versions of Tersoff and REBO-II, LCBOP-I, 2015 and 2020 versions of ReaxFF, and the machine-learning GAP force field. We show that the ReaxFF force fields give the best performance for the interaction energies of the cabon onions relative to density functional theory interaction energies obtained at the PBE0-D3/def2-TZVP level of theory. We proceed to use the ReaxFF-15 force field to explore the interaction energies in a giant ten-layered carbon onion with a C60 core and show that the interaction energy between the outer layer and the inner layers increases linearly with the number of layers in the carbon onion (with a squared correlation coefficient of R2 = 0.9996). This linear increase in the stabilization energy with each consecutive layer may have important thermodynamic consequences for describing the formation and growth of large carbon onions.

It is well established that the basis set convergence of the correlation consistent (cc-pV

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 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 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.

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

Stable equilibrium compounds containing a planar antiaromatic cyclooctatetraene (COT) ring are promising candidates for organic electronic devices such as organic semiconductor transistors. The planarization of COT by incorporation into rigid planar π-systems, as well as by oxidation or reduction has attracted considerable attention in recent years. Using dispersion-corrected density functional theory calculations, we explore an alternative approach of planarizing COT derivatives by adsorption onto graphene. We show that strong π-π stacking interactions between graphene and COT derivatives induce a planar structure with an antiaromatic central COT ring. In addition to being reversible, this strategy provides a novel approach for planarizing COT without the need for incorporation into a rigid structure, atomic substitution, oxidation, or reduction.

Density functional theory (DFT) has been extensively benchmarked for energetic properties; however, less attention has been given to equilibrium structures and the effect of using a certain DFT geometry on subsequent energetic properties. We evaluate the performance of 52 contemporary DFT methods for obtaining the structures of 122 species in the W4-11-GEOM database. This dataset includes a total of 246 unique bonds: 117 H─X, 65 X─Y, 49 X═Y, and 15 XY bonds (where X and Y are first- and second-row atoms) and 133 key bond angles: 96 X-Y-H, 22 X-Y-Z, and 15 H-X-H angles. The reference geometries are optimized at the CCSD(T)/jul-cc-pV(n+d)Z level of theory (n = 5, 6). The performance of DFT is evaluated in conjunction with the Def2-nZVPP (n = T, Q), cc-pV(T+d)Z, and jul-cc-pV(T+d)Z basis sets. The root-mean-square deviations (RMSDs) over the bond distances of the best performing functionals from each rung of Jacob's Ladder are 0.0086 (SOGGA11), 0.0088 (τ-HCTH), 0.0059 (B3LYP), 0.0054 (TPSSh), and 0.0032 (DSD-PBEP86) A. We evaluate the effect of the choice of the DFT geometry on subsequent molecular energies calculated with W1-F12 theory. Geometries obtained with GGA and MGGA methods result in large RMSDs in the subsequent W1-F12 energies; however, six hybrid GGA functionals (B3LYP, B3P86, mPW3PBE, B3PW91, mPW1LYP, and X3LYP) result in an excellent performance with RMSDs between 0.25 and 0.30 kJ mol-1 relative to the CCSD(T)/CBS reference geometries. The B2GP-PLYP and mPW2-PLYP DHDFT methods result in near-CCSD(T) accuracy with RMSDs of 0.11 and 0.10 kJ mol-1 , respectively.

Reactions involving Criegee intermediates (CIs, R1 R2 COO) are important in atmospheric ozonolysis models. In recent years, density functional theory (DFT) and CCSD(T)-based ab initio methods are increasingly being used for modeling reaction profiles involving CIs. We obtain highly accurate CCSDT(Q)/CBS reaction energies and barrier heights for ring-closing reactions involving atmospherically important CIs (R1 /R2 = H, Me, OH, OMe, F, CN, cyclopropene, ethylene, acetaldehyde, and acrolein). We use this benchmark data to evaluate the performance of DFT, double-hybrid DFT (DHDFT), and ab initio methods for the kinetics and thermodynamics of these reactions. We find that reaction energies are more challenging for approximate theoretical procedures than barrier heights. Overall, taking both reaction energies and barrier heights into account, only one of the 58 considered DFT methods (the meta-GGA MN12-L) attains near chemical accuracy, with root-mean-square deviations (RMSDs) of 3.5 (barrier heights) and 4.7 (reaction energies) kJ mol-1 . Therefore, MN12-L is recommended for investigations where CCSD(T)-based methods are not computationally feasible. For reaction barrier heights performance does not strictly follow Jacob's Ladder, for example, DHDFT methods do not perform better than conventional DFT methods. Of the ab initio methods, the cost-effective CCSD(T)/CBS(MP2) approach gives the best performance for both reaction energies and barrier heights, with RMSDs of 1.7 and 1.4 kJ mol-1 , respectively. All the considered Gaussian-n methods show good performance with RMSDs below the threshold of chemical accuracy for both reaction energies and barrier heights, where G4(MP2) shows the best overall performance with RMSDs of 2.9 and 1.5 kJ mol-1 , respectively. © 2019 Wiley Periodicals, Inc.

Inspired by the promise of hydrogen (H2) as a clean alternate to the existing energy sources, we have employed spin-polarized density functional theory calculations on a recently designed two-dimensional C 4 N monolayer as a promising H 2 storage material. By means of first principles DFT calculations, we have comprehensively studied the geometric and electronic properties of pristine, defected and metal-doped C 4 N nanosheets and further explored their H 2 storage properties. We found that light metal dopants such as Li, Na, K, Mg, and Ca bind strongly to defects on a C 4 N nanosheet with binding energies of 3–4 eV per dopant. These binding energies are sufficiently strong to surpass metal clustering. Thermal stability of the metal-doped C 4 N nanosheets has been further verified by means of ab initio molecular dynamics simulations. The bonding nature of the metal dopants with the C 4 N nanosheet has been studied through Bader analysis and Roby-Gould methods and the electronic properties were studied through density of states. We found that each dopant in the metal-doped C 4 N nanosheet can bind up to five H 2 molecules with adsorption energies ranging between 0.15 and 0.60 eV/ H 2 , which results in optimal H 2 storage capacities. Finally, we employed thermodynamic analysis to investigate the H 2 adsorption/desorption mechanism under practical operating conditions.

The quest for a suitable material with the potential of capturing toxic nitrogen-containing gases (NH3, NO, and NO2) has motivated us to explore the structural, electronic, and gas-sensing properti ...

Based on spin-polarized DFT calculations, we have studied the interaction mechanism of recently synthesized blue phosphorene (BlueP) monolayers towards selected key volatile organic compounds (VOCs) such as acetone, ethanol and propanal. Our binding energy analysis shows that pristine BlueP weakly binds the VOCs and that this binding does not appreciably change the electronic properties of the monolayer – a prerequisite for any sensing material. However, mono, di, and tri-vacancy defects and Si/S-substitutional doping significantly enhance the binding energies with VOCs. Density of state (DOS) calculations show that upon adsorption of VOCs, mono-vacancy and S-substituted BlueP monolayers undergo a major change in electronic structure, which make them potential candidates for VOCs sensing materials. By contrast, binding of VOCs to di- and tri-vacancy and Si-substitution sites does not alter the electronic structure of BlueP monolayers drastically, therefore, are not qualified for VOCs sensing applications.

The ability to accurately calculate reaction barrier heights is of central importance to many areas of chemistry. We report an extensive study examining the basis set convergence of post-CCSD(T) contributions (up to CCSDT(Q)) for a diverse set of 28 reaction barrier heights. In contrast to previous studies, we focus here on larger transition structures (TSs) involving 4-7 non-hydrogen atoms. The set of reaction barrier heights includes pericyclic, bipolar cycloaddition, cycloreversion, and multiple-proton transfer reactions. We find that in most cases post-CCSD(T) contributions converge rapidly toward the basis set limit, such that even double-ζ and truncated double-ζ basis sets provide useful estimates of the T-(T) and (Q) contributions, respectively. In addition, we find that due to the tendency of these small basis sets to systematically underestimate the T-(T) and (Q) components, scaling is an effective approach for improving performance. For example, scaling the T-(T)/cc-pVDZ contribution by 1.25 results in an RMSD of merely 0.4 kJ mol

The equilibrium geometries and spectroscopic properties of two key singlet carbenes, buta-1,3-diynylcarbene (6) and 2-methylenebicyclo[1.1.0]but-1(3)-en-4-ylidene (9), which have not been experimentally observed to date, are investigated using high-level coupled-cluster (CC) methods. The current theoretical study necessitates new experimental data on C5H2 isomers considering the relevance of these molecules to interstellar chemistry. Bent-pentadiynylidene (4) has been missing in the laboratory and the prime focus of our earlier theoretical work. The present theoretical study indicates that isomers 6 and 9 are also viable experimental targets. Apart from ethynylcyclopropenylidene (2), pentatetraenylidene (3), ethynylpropadienylidene (5), and 3-(didehydrovinylidene)cyclopropene (8), which are identified by Fourier transform microwave spectroscopy, the dipole moments of elusive 4, 6, and 9 are also nonzero (μ ≠ 0). The relative energies of these isomers, calculated at the CCSDT(Q)/CBS level of theory, with respect to linear triplet pentadiynylidene (1) reveal that they all lie within 25.1 kcal mol-1. Therefore, geometric, energetic, aromatic, and spectroscopic parameters are reported here, which may assist the efforts of molecular spectroscopists in the future. Anharmonic vibrational calculations on isomers 6 and 9 indicate that the former is loosely bound and would be challenging to be detected experimentally. Among the undetected carbenes, 9 may be considered as a potential target molecule considering its higher polarity and aromatic nature.

The observation of an unusual crystal habit in the common diuretic drug hydrochlorothiazide (HCT), and identification of its subtle conformational chirality, has stimulated a detailed investigation of its crystalline forms. Enantiomeric conformers of HCT resolve into an unusual structure of conjoined enantiomorphic twin crystals comprising enantiopure domains of opposite chirality. The purity of the domains and the chiral molecular conformation are confirmed by spatially revolved synchrotron micro-XRD experiments and neutron diffraction, respectively. Macroscopic inversion twin symmetry observed between the crystal wings suggests a pseudoracemic structure that is not a solid solution or a layered crystal structure, but an unusual structural variant of conglomerates and racemic twins. Computed interaction energies for molecular pairs in the racemic and enantiopure polymorphs of HCT, and the observation of large opposing unit-cell dipole moments for the enantiopure domains in these twin crystals, suggest a plausible crystal nucleation mechanism for this unusual crystal habit.

International audience; Nitrogen Monoxide (NO) is observed in the gas phase of molecular clouds. It may accrete on dust grains and there its hydrogenation should lead to hydroxylamine (NH 2 OH), the same way that CO is transformed in methanol (CH 3 OH) on the surface dust grains. NO hydrogenation has been said barrier-less, whereas CO hydrogenation proceed through quantum tunneling, and is thus slower. However, CH 3 OH is widely observed and is considered as a proxy of complex organic molecules while hydroxy-lamine remains undetected. We aim at studying, analyzing, and understanding the 1 chemical network of NO hydrogenation on cold surfaces. Experiments are carried out using a new Ultra-Hight Vacuum (UHV) setup named VENUS. NO molecules and H atoms are co-deposited on a golden mirror at different temperatures. Infrared spectroscopy as well as Temperature Programmed Desorption (TPD) are used to follow the NO reactivity, with both H and D, and in presence or absence of water substrate. Quantum calculations on water ice cluster models are computed separately. During the hydrogenation of NO, 10 reactions proceed concurrently. They are identified and constrained by changing physical conditions in experiments or in calculations. Among them, we demonstrate that the HNO+H addition reaction has a barrier which is probably crossed via quantum tunneling at 10 K. Moreover, abstraction reactions are occurring although they are limited by H and O bonding with their environment. Chemical desorption should occur especially in absence of water which enhances the total production of hydroxylamine. The chemical network of the hydrogenation of NO has been re-investigated. Each of the 10 reactions are sorted by efficiency. We exclude the possibility of a chemical loop between NO and HNO, especially in presence of water. Therefore hydroxylamine remains the main product of the hydrogenation of NO on grains and the question of its non-detection in ices or in the gas phase, specifically in shocked regions where ice mantles should be sputtered, is still open.

The enthalpies of formation and isomerization energies of P4Sn molecular cages are not experimentally (or theoretically) well known. We obtain accurate enthalpies of formation and isomerization energies for P4Sn cages (n = 3, 4, 5, 6, and 10) by means of explicitly correlated high-level thermochemical procedures approximating the CCSD(T) and CCSDT(Q) energies at the complete basis set (CBS) limit. The atomization reactions have very significant contribution from post-CCSD(T) correlation effects and, due to the presence of many second-row atoms, the CCSD and (T) correlation energies converge exceedingly slowly with the size of the one-particle basis set. As a result, these cage structures are challenging targets for thermochemical procedures approximating the CCSD(T) energy (e.g., W1-F12 and G4). Our best enthalpies of formation at 298 K (∆fH°298) are obtained from thermochemical cycles in which the P4Sn cages are broken down into P2S2 and S2 fragments for which highly accurate ∆fH°298 values are available from W4 theory. For the smaller P4S3 and P4S4 cages, the reaction energies are calculated at the CCSDT(Q)/CBS level and for the larger P4S5, P4S6, and P4S10 cages, they are obtained at the CCSD(T)/CBS level. Our best ∆fH°298 values are − 94.5 (P4S3), − 108.4 (α-P4S4), − 98.7 (β-P4S4), − 126.2 (α-P4S5), − 126.1 (β-P4S5), − 112.7 (γ-P4S5), − 144.7 (α-P4S6), − 153.9 (β-P4S6), − 134.4 (γ-P4S6), − 136.3 (δ-P4S6), − 118.7 (e-P4S6), and − 215.4 (P4S10) kJ mol−1. Interestingly, we find a linear correlation (R2 = 0.992) between the enthalpies of formation of the most stable isomers of each molecular formula and the number of atoms in the P4Sn cages. We use our best ∆fH°298 values to assess the performance of a number of lower-cost composite ab initio methods. For absolute enthalpies of formation, G4(MP2) and G3(MP2)B3 result in the best overall performance with root-mean-square deviations (RMSDs) of 10.6 and 12.9 kJ mol−1, respectively, whereas G3, G3B3, and CBS-QB3 result in the worst performance with RMSDs of 27.0–38.8 kJ mol−1. In contrast to absolute enthalpies of formation, all of the considered composite procedures give a good-to-excellent performance for the isomerization energies with RMSDs below the 5 kJ mol−1 mark.

Motivated by the promise of two-dimensional nanostructures in the field of gas sensing, we have employed van der Waals corrected density functional theory calculations to study the structural, electronic and gas sensing propensities of the recently designed Si2BN monolayer. Our rigorous simulations reveal that the representative members of nitrogen-containing gases (NCGs) such as NO, NO2 and NH3 binds extremely strongly on pristine Si2BN monolayer. However, a strong dissociative adsorption in case of NO and NO2 would poison the Si2BN and ultimately reversibility of the monolayer would be compromised. Exploring the sensing mechanism in more realistic pre-oxidized conditions, the binding characteristics of O2@Si2BN changed dramatically, resulting into much lower adsorption in associative manner for all NO, NO2 and NH3. A visible change in work function indicates the variation in conductivity of O2@Si2BN upon the exposure of incident gases. Sustainable values of binding energies would also ensure a quick recovery time that makes O2@Si2BN an efficient nano sensor for pollutants like NCGs.

The pursuit of a material capable of storing a high capacity of hydrogen (H2) efficiently has prompted us to study the structural, electronic and H2 storage properties of recently designed two-dimensional BN2 nanosheets. Our spin-polarized density functional theory based calculations have revealed that the pristine BN2 barely anchor H2 molecules, however, alkali metal (AM) doping enhances the binding energies drastically. Van der Waals corrected energetics analysis established a uniform distribution of AMs over the BN2 monolayers even at a high doping concentration of 12.50%, which ensure the reversibility of the systems. Bader charge analysis, Roby-Gould bond index method, and electron localization function isosurfaces conclude the transfer of charges from AMs to BN2, which has resulted into strong ionic bonds between the former and the latter. The presence of partial positive charges on each of the AMs would adsorb multiple H2 molecules with binding energies that are ideal for mobile H2 storage applications. Considerably high H2 storage capacities of 6.75%, 6.87% and 6.55% could be achieved with 3Li@BN2, 3Na@BN2 and 3K@BN2 systems, respectively that guarantees the promise of AMs decorated BN2 as a promising H2 storage material.

Nitroxide radicals have previously been successfully used as electrodes in all-organic radical batteries. However, one drawback of these batteries is significantly reduced redox potentials, in comparison to that of widely used lithium-ion batteries, making their energy-producing capacity rather small for use as a primary battery. In addition, strong propensity of nitroxide radicals to engage in side reactions with traditional electrolytes based on molecular solvents give rise to a series of undesirable and irreversible byproducts, thus significantly reducing the life of nitroxide batteries. Ionic liquids (ILs) have previously demonstrated their ability to reduce the reactivity of radicals through strong intermolecular interactions. In this study, we investigate the use of ILs as electrolytes with the view of increasing redox potentials of nitroxide radicals. A series of imidazolium, phosphonium, and pyrrolidinium-based ILs coupled with widely used anions were chosen to predict redox potentials of the 2,2,...

Hypochlorous acid (HOCl) is a highly reactive, toxic species generated by neutrophils via the action of myeloperoxidase in order to destroy invading pathogens. However, when HOCl is produced inappropriately, it can damage host tissue and proteins and plays a role in the initiation and progression of disease. Carnosine, a peptide of β-alanine and histidine, has been shown to react rapidly with HOCl yielding monochloramines and can undergo intramolecular transchlorination. The current study examines the kinetics and pH dependence of the reactions of carnosine and novel structural derivatives with HOCl and the occurrence of intra- and intermolecular transchlorination processes. We demonstrate that the transchlorination reactions of carnosine are pH dependent, with intramolecular transfer favored at higher pH. Carcinine, having a structure identical to carnosine though lacking the carboxylic acid group of the histidine residue, reacts with HOCl and forms monochloramines though intramolecular transfer reactions are not observed, and this is supported by computational modeling. Novel analogues with one (carnosine+1) and two (carnosine+2) methylene groups in the alkyl chain of the β-alanine react with HOCl to yield monochloramines that undergo transchlorinations to yield a mixture of mono- and dichloramines. The latter are stable over 24 h. The ability of carnosine and derivatives to react rapidly with HOCl to give long-lived, poorly reactive, species may prevent damage to proteins and other targets at sites of inflammation.

Polymers are an attractive anchoring platform for the synthesis of radioimmunoconjugates. They enable independent control over the amount of radioisotope loading and antibody attachment, which is pivotal in developing tailorable formulations for personalised medicine. Herein, we report the synthesis of p(HEMA-ran-GMA) for the conjugation of lutetium ions and rituximab as a functional platform for radioimmunotherapy. We demonstrate the suitability of this platform using non-Hodgkin's lymphoma cells.

Perylenequinones are a class of aromatic polyketides characterised by a highly conjugated pentacyclic core, which confers them with potent light-induced bioactivities and unique photophysical properties. Despite the biosynthetic gene clusters for the perylenequinones elsinochrome A (1), cercosporin (4) and hypocrellin A (6) being recently identified, key biosynthetic aspects remain elusive. Here, we first expressed the intact elc gene cluster encoding 1 from the wheat pathogen Parastagonospora nodorum heterologously in Aspergillus nidulans on a yeast-fungal artificial chromosome (YFAC). This led to the identification of a novel flavin-dependent monooxygenase, ElcH, responsible for oxidative enolate coupling of a perylenequinone intermediate to the hexacyclic dihydrobenzo(ghi)perylenequinone in 1. In the absence of ElcH, the perylenequione intermediate formed a hexacyclic cyclohepta(ghi)perylenequinone system via an intramolecular aldol reaction resulting in 6 and a novel hypocrellin 12 with opposite helicity to 1. Theoretical calculations supported that 6 and 12 resulted from atropisomerisation upon formation of the 7-membered ring. Using a bottom-up pathway reconstruction approach on a tripartite YFAC system developed in this study, we uncovered that both a berberine bridge enzyme-like oxidase ElcE and a laccase-like multicopper oxidase ElcG are involved in the double coupling of two naphthol intermediates to form the perylenequinone core. Gene swapping with the homologs from the biosynthetic pathway of 4 showed that cognate pairing of the two classes of oxidases is required for the formation of the perylenequinone core, suggesting the involvement of protein–protein interactions.

Herein we demonstrated that the shape and phase of colloidal ZnSe nanocrystals can be spontaneously tuned through tailoring the selenium precursor reactivity in a phosphine-free reaction system. Selenium species with diverse reaction activities, i.e. Se22− or Se2−, were produced by the addition of different volumes of reductant superhydride (lithium triethylborohydride). Theoretical calculation of ΔGr indicates that superhydride-reduced Se22− is less active than Se2− for the reaction with a Zn precursor. Nanoparticle growth using Se22− produces wurtzite ZnSe nanorods whereas the reaction of more reactive Se2− with the zinc precursor leads to the formation of spherical zinc blende ZnSe nanodots. This work not only provides a facile synthetic approach for the preparation of high quality ZnSe nanocrystals but also gives insights into the shape and phase control of other colloidal nanocrystals.