The ORCA quantum chemistry program package

Single Excitation Energies Obtained from the Ensemble "HOMO-LUMO Gap": Exact Results and Approximations

In calculations based on density functional theory, the "HOMO-LUMO gap" (difference between the highest occupied and lowest unoccupied molecular orbital energies) is often used as a low-cost, ad hoc approximation for the lowest excitation energy. Here we show that a simple correction based on rigorous ensemble density functional theory makes the HOMO-LUMO gap exact in principle and significantly more accurate in practice. The introduced perturbative ensemble density functional theory approach predicts different and useful values for singlet-singlet and singlet-triplet excitations, using semilocal and hybrid approximations. Excitation energies are similar in quality to time-dependent density functional theory, especially at high fractions of exact exchange. The approach therefore offers an easy-to-implement and low-cost route to robust prediction of molecular excitation energies.

The Apparently Unreactive Substrate Facilitates the Electron Transfer for Dioxygen Activation in Rieske Dioxygenases

Rieske dioxygenases belong to the non-heme iron family of oxygenases and catalyze important cis-dihydroxylation as well as O-/N-dealkylation and oxidative cyclization reactions for a wide range of substrates. The lack of substrate coordination at the non-heme ferrous iron center, however, makes it particularly challenging to delineate the role of the substrate for productive O 2 activation. Here, we studied the role of the substrate in the key elementary reaction leading to O 2 activation from a theoretical perspective by systematically considering (i) the 6-coordinate to 5-coordinate conversion of the non-heme FeII upon abstraction of a water ligand, (ii) binding of O 2 , and (iii) transfer of an electron from the Rieske cluster. We systematically evaluated the spin-state-dependent reaction energies and structural effects at the active site for all combinations of the three elementary processes in the presence and absence of substrate using naphthalene dioxygenase as a prototypical Rieske dioxygenase. We find that reaction energies for the generation of a coordination vacancy at the non-heme FeII center through thermoneutral H2 O reorientation and exothermic O 2 binding prior to Rieske cluster oxidation are largely insensitive to the presence of naphthalene and do not lead to formation of any of the known reactive Fe-oxygen species. By contrast, the role of the substrate becomes evident after Rieske cluster oxidation and exclusively for the 6-coordinate non-heme FeII sites in that the additional electron is found at the substrate instead of at the iron and oxygen atoms. Our results imply an allosteric control of the substrate on Rieske dioxygenase reactivity to happen prior to changes at the non-heme FeII in agreement with a strategy that avoids unproductive O 2 activation.

A Hierarchical Theoretical Study of the Hydrogen Abstraction Reactions of H2/C1-C4 Molecules by the Methyl Peroxy Radical and Implications for Kinetic Modeling

The hydrogen atom abstraction by the methyl peroxy radical (CH₃O₂) is an important reaction class in detailed chemical kinetic modeling of the autoignition properties of hydrocarbon fuels. Systematic theoretical studies are performed on this reaction class for H₂/C₁-C₄ fuels, which is critical in the development of a base model for large fuels. The molecules include hydrogen, alkanes, alkenes, and alkynes with a carbon number from 1 to 4. The B2PLYP-D3/cc-pVTZ level of theory is employed to optimize the geometries of all of the reactants, transition states, and products and also the treatments of hindered rotation for lower frequency modes. Accurate benchmark calculations for abstraction reactions of hydrogen, methane, and ethylene with CH₃O₂ are performed by using the coupled cluster method with explicit inclusion of single and double electron excitations and perturbative inclusion of triple electron excitations (CCSD(T)), the domain-based local pair-natural orbital coupled cluster method (DLPNO-CCSD(T)), and the explicitly correlated CCSD(T)-F12 method with large basis sets. Reaction rate constants are computed via conventional transition state theory with quantum tunneling corrections. The computed rate constants are compared with literature values and those employed in detailed chemical kinetic mechanisms. The calculated rate constants are implemented into the recently developed NUIGMECH1.1 base model for kinetic modeling of ignition properties.

Photoinduced bond oscillations in ironpentacarbonyl give delayed synchronous bursts of carbonmonoxide release

Early excited state dynamics in the photodissociation of transition metal carbonyls determines the chemical nature of short-lived catalytically active reaction intermediates. However, time-resolved experiments have not yet revealed mechanistic details in the sub-picosecond regime. Hence, in this study the photoexcitation of ironpentacarbonyl Fe(CO)₅ is simulated with semi-classical excited state molecular dynamics. We find that the bright metal-to-ligand charge-transfer (MLCT) transition induces synchronous Fe-C oscillations in the trigonal bipyramidal complex leading to periodically reoccurring release of predominantly axial CO. Metaphorically the photoactivated Fe(CO)₅ acts as a CO geyser, as a result of dynamics in the potential energy landscape of the axial Fe-C distances and non-adiabatic transitions between manifolds of bound MLCT and dissociative metal-centered (MC) excited states. The predominant release of axial CO ligands and delayed release of equatorial CO ligands are explained in a unified mechanism based on the σ^*(Fe-C) anti-bonding character of the receiving orbital in the dissociative MC states.

The photodissociation of transition metal carbonyls is involved in catalysis and synthetic processes. Here the authors, using semi-classical excited state molecular dynamics, observe details of the early stage dynamics in the photodissociation of Fe(CO)₅, including synchronous bursts of CO at periodic intervals of 90 femtoseconds.

Dissociative electron attachment to hexachlorobenzene

Gas phase molecules of hexachlorobenzene (C 6 Cl 6 ) were investigated by means of Dissociative Electron Attachment Spectroscopy (DEAS). Three channels of molecular negative ions decay have been recorded: Cl- and Cl 2 - as well as electron detachment (t a ~250 ms at 343 K). All of these decay channels exhibit temperature dependence. The electron affinity was estimated within the framework of the simple Arrhenius model (EA=1.6-1.9 eV), as well as quantum-chemically at various levels (EA a =0.9-1.0 eV). The reasons for the discrepancy between the obtained estimates of the EA value and the results of previous works (EA from 1 to 2 eV) are discussed.

Tunable Excimer Circularly Polarized Luminescence in Isohexide Derivatives from Renewable Resources

Organic compounds showing circularly polarized luminescence (CPL) are at the forefront of novel applications and technologies. Here we show the synthesis and chiroptical properties of pyrene and perylene derivatives of inexpensive chiral scaffolds: isomannide and isosorbide. Low‐intensity ECD spectra were obtained, suggesting the absence of chromophore interaction in the ground state, except in the case of isomannide bis‐perylenecarboxylate, whose ECD spectrum showed a positive exciton couplet. All isomannide derivatives, with the only exception of the one containing a pyrenecarboxylate and a perylenecarboxylate, exhibited excimer CPL spectra, whereas isosorbide derivatives did not show any CPL. Isomannide derivatives bearing two pyrenecarboxylate or two pyrenylacetate groups showed positive CPL emission with dissymmetry factors up to 10⁻², which depends on the conformational freedom of the appended units. The CPL sign, Stokes shift and order of magnitude of dissymmetry factor were reproduced by excited‐state calculations on a representative compound. Interestingly, the mixed derivative containing pyrenic units with different spacing from the isomannide scaffold showed an oppositely signed excimer band with respect to the homo‐substituted derivatives.

Effects of Salt Concentration on the Water and Ion Self-Diffusion Coefficients of a Model Aqueous Sodium-Ion Battery Electrolyte

The aqueous sodium-ion battery is a promising alternative to the well-known lithium-ion battery owing to the large abundance of sodium ion resources. Although it is safer than the lithium-ion battery, the voltage window of the sodium-ion battery is narrower than that of the lithium-ion battery, thus limiting its practical implementation. Therefore, a highly concentrated electrolyte is required to address this issue. In the present work, the effect of the salt concentration on the transport properties of water molecules is investigated via theoretical analyses at the quantum mechanical level. A molecular dynamics simulation at the quantum mechanical level revealed that as the salt concentration increases, the ion-water interactions became stronger, leading to a lower diffusivity and a lower electronic band gap. These imply that the superconcentrated aqueous-based electrolytes have high potentials for the sodium-ion battery applications.

Vibration-Assisted Intersystem Crossing in the Ultrafast Excited-State Relaxation Dynamics of Halocoumarins

Due to its numerous applications, triplet formation and resulting phosphorescence remain a frontier area of research for over eight decades. Facile intersystem crossing (ISC) is the primary requirement for triplet formation and observation of phosphorescence. The incorporation of a heavy atom in molecules is one of the common approaches employed to facilitate ISC. A detailed study of the excited state dynamics that governs ISC is necessary to understand the mechanism of heavy atom effect (HAE). Incorporation of iodine at the 3 position of coumarin-1 reduces fluorescence quantum yield (ϕ_f) drastically as expected, whereas bromine substitution at the same position increased the ϕ_f. Such a contrasting effect of the two heavy atoms suggests that there are other features yet to be discovered to fully understand the HAE. Detailed steady state and femtosecond transient absorption studies along with theoretical calculations suggest that the C3-X (X = Br, I) bond vibration plays an important role in the ISC process. The study reveals that while in the case of the iodo-derivative there is no energy barrier in the singlet triplet crossing path, there is a barrier in the case of the bromo-derivative, which slows the ISC process. Such an unexpected phenomenon is not limited to halocoumarins as this rationalizes the photobehavior of 1-bromo-/iodo-substituted naphthalenes as well.

Comparison of rhodamine 6G, rhodamine B and rhodamine 101 spirolactam based fluorescent probes: A case of pH detection

Ring-opening reaction of rhodamine spirolactam has been widely applied to construct fluorescent probes. The fluorescence properties of the probe were finely tuned for specific purpose through changing the rhodamine fluorophore. However, the influence on response range and kinetic parameters of the probe during the change has been seldom discussed. Herein, we took pH detection as an example and constructed spirolactam based probes (RLH A-C) with Rhodamine 6G, Rhodamine B and Rhodamine 101. The pK_a values and observed rate constant k_obs of RLH A-C were determined and found to negatively correlated with the calculated Gibbs free energy differences ΔG_C-O and ΔG_TS respectively. The potential applications of RLH A-C in imaging acidic microenvironment were also investigated in cells. We expect the comparison of rhodamine fluorophores will facilitate the quantitative optimization of rhodamine spirolactam based fluorescent probes.

The PM6-FGC Method: Improved Corrections for Amines and Amides

Recently, we reported a new approach to develop pairwise analytical corrections to improve the description of noncovalent interactions, by approximate methods of electronic structures, such as semiempirical quantum mechanical (SQM) methods. In particular, and as a proof of concept, we used the PM6 Hamiltonian and we named the method PM6-FGC, where the FGC acronym, corresponding to Functional Group Corrections, emphasizes the idea that the corrections work for specific functional groups rather than for individual atom pairs. The analytical corrections were derived from fits to B3LYP-D3/def2-TZVP (reference). PM6 interaction energy differences, evaluated for a reduced set of small bimolecular complexes, were chosen as representatives of saturated hydrocarbons, carboxylic, amine and, tentatively, amide functional groups. For the validation, the method was applied to several complexes of well-known databases, as well as to complexes of diglycine and dialanine, assuming the transferability of amine group corrections to amide groups. The PM6-FGC method showed great potential but revealed significant inaccuracies for the description of some interactions involving the –NH₂ group in amines and amides, caused by the inadequate selection of the model compound used to represent these functional groups (an NH₃ molecule). In this work, methylamine and acetamide are used as representatives of amine and amide groups, respectively. This new selection leads to significant improvements in the calculation of noncovalent interactions in the validation set.

Improved First-Principles Model of Differential Mobility Using Higher Order Two-Temperature Theory

Differential mobility spectrometry is a separation technique that may be applied to a variety of analytes ranging from small molecule drugs to peptides and proteins. Although rudimentary theoretical models of differential mobility exist, these models are often only applied to small molecules and atomic ions without considering the effects of dynamic microsolvation. Here, we advance our theoretical description of differential ion mobility in pure N₂ and microsolvating environments by incorporating higher order corrections to two-temperature theory (2TT) and a pseudoequilibrium approach to describe ion-neutral interactions. When comparing theoretical predictions to experimentally measured dispersion plots of over 300 different compounds, we find that higher order corrections to 2TT reduce errors by roughly a factor of 2 when compared to first order. Model predictions are accurate especially for pure N₂ environments (mean absolute error of 4 V at SV = 4000 V). For strongly clustering environments, accurate thermochemical corrections for ion-solvent clustering are likely required to reliably predict differential ion mobility behavior. Within our model, general trends associated with clustering strength, solvent vapor concentration, and background gas temperature are well reproduced, and fine structure visible in some dispersion plots is captured. These results provide insight into the dynamic ion-solvent clustering process that underpins the phenomenon of differential ion mobility.

A comparative study on the stability of the furfural molecule on the low index Ni, Pd and Pt surfaces

We present a comparative density functional theory investigation of the furfural (Ff) molecule on the low index Ni, Pd and Pt surfaces to understand its geometrical and electronic properties to gain mechanistic insights into the experimentally measured catalytic reactivities of these metal catalysts. We show that the number of metal d-states, which hybridize with the nearest C and O p-orbitals of the Ff molecule, can be used to explain the stability of the Ff molecule on these surfaces. We find that the hybridization between atoms with higher electronegativity and the metal d-states plays a crucial role in determining the stability of these systems. Furthermore, we also find electron transfer from metal to the Ff molecule on the Ni and Pd surfaces, with a reverse process occurring on the Pt surface.

CalcUS: An Open-Source Quantum Chemistry Web Platform

Computational chemistry is an increasingly active field due to the improvement of computing resources and theoretical tools. However, its use remains usually limited to technically inclined users due to the technical challenges of preparing, launching, and analyzing calculations. In this context, we have developed CalcUS, an open-source platform to streamline quantum chemistry studies. Its objective is to democratize access to computational chemistry by providing a user-friendly web interface to simplify running and analyzing quantum mechanical calculations. It is freely available, expandable, and customizable. It promotes connectivity to multiple software packages and algorithms, thus providing state-of-the-art techniques to all practitioners. We propose CalcUS as a standalone tool and infrastructure to support other open-source packages.

Wearable Light Sensors Based on Unique Features of a Natural Biochrome

Overexposure to complete solar radiation (combined ultraviolet, visible, and infrared) is correlated with several harmful biological consequences including hyperpigmentation, skin cancer, eye damage, and immune suppression. With limited effective therapeutic options available for these conditions, significant efforts have been directed toward promoting preventative habits. Recently, wearable solar radiometers have emerged as practical tools for managing personal exposure to sunlight. However, designing simple and inexpensive sensors that can measure energy across multiple spectral regions without incorporating electronic components remains challenging, largely due to inherent spectral limitations of photoresponsive indicators. In this work, we report the design, fabrication, and characterization of wearable radiation sensors that leverage an unexpected feature of a natural biochrome, xanthommatin-its innate sensitivity to both ultraviolet and visible through near-infrared radiation. We found that xanthommatin-based sensors undergo a visible shift from yellow to red in the presence of complete sunlight. This color change is driven by intrinsic photoreduction of the molecule, which we investigated using computational modeling and supplemented by radiation-driven formation of complementary reducing agents. These sensors are responsive to dermatologically relevant doses of erythemally weighted radiation, as well as cumulative doses of high-energy ultraviolet radiation used for germicidal sterilization. We incorporated these miniature sensors into pressure-activated microfluidic systems to illustrate on-demand activation of a wearable and mountable form factor. When taken together, our findings encompass an important advancement toward accessible, quantitative measurements of UVC and complete solar radiation for a variety of use cases.

How Correlations and Spin-Orbit Coupling Work within Extended Orbitals of Transition-Metal Tetrahedra of 4d/5d Lacunar Spinels

Spin-orbit quartet ground states are associated with rich phenomenology, ranging from multipolar phases in f¹ rare-earth borides to magnetism emerging through covalency and vibronic couplings in d¹ transition-metal compounds. The latter effect has been studied since the 1960s on t_2g¹ octahedral ML₆ units in both molecular complexes and extended solid-state lattices. Here we analyze the J_eff = 3/2 quartet ground state of larger cubane-like M₄L₄ entities in lacunar spinels, composed of transition-metal (M) tetrahedra caged by chalcogenide ligands (L). These represent a unique platform where spin-orbit coupling acts on molecular-like, delocalized t₂ orbitals. Using quantum chemical methods, we pin down the interplay of spin-orbit couplings in such a setting and many-body physics related to other molecular-like single-electron levels, both below and above the reference t₂¹. We provide a different interpretation of resonant inelastic X-ray scattering data on GaTa₄Se₈ and, by comparing magnetic susceptibility data with calculated g factors, valuable insights into the important role of vibronic couplings.

Trigonal Bipyramidal Rhodium(I) Methyl and Phenyl Complexes: Precursors of Oxidative Methyl and Phenyl Radical Generation

The new complexes [Rh(Me)(trop3P)] (2) and [Rh(Ph)(trop3P)] (3) (trop = 5H-dibenzo[a,d]cyclohepten-5-yl) were synthesised by addition of organolithium reagents (MeLi and PhLi) to the parent pentacoordinated chloride complex [RhCl(trop3P)]. The compounds have a trigonal bipyramidal structure with olefin-only ligands in the equatorial position and the methyl or phenyl substituent in the axial position. Oxidation of complexes 2 and 3 leads to the liberation of methyl and phenyl radicals, which were indirectly detected by reaction with common spin trapping reagents.

Mapping the Finite-Temperature Behavior of Conformations to Their Potential Energy Barriers: Case Studies on Si6B and Si5B Clusters

Dynamical simulations of molecules and materials have been the route to understand the rearrangement of atoms within them at different temperatures. Born-Oppenheimer molecular dynamical simulations have further helped to comprehend the reaction dynamics at various finite temperatures. We take a case study of Si₆B and Si₅B clusters and demonstrate that their finite-temperature behavior is rather mapped to the potential energy surface. The study further brings forth the fact that an accurate description of the dynamics is rather coupled with the accuracy of the method in defining the potential energy surface. A more precise potential energy surface generated through the coupled cluster method is finally used to identify the most accurate description of the potential energy surface and the interconnected finite-temperature behavior.

Ionization Energies and Redox Potentials of Hydrated Transition Metal Ions: Evaluation of Domain-Based Local Pair Natural Orbital Coupled Cluster Approaches

Hydrated transition metal ions are prototypical systems that can be used to model properties of transition metals in complex chemical environments. These seemingly simple systems present challenges for computational chemistry and are thus crucial in evaluations of quantum chemical methods for spin-state and redox energetics. In this work, we explore the applicability of the domain-based pair natural orbital implementation of coupled cluster (DLPNO-CC) theory to the calculation of ionization energies and redox potentials for hydrated ions of all first transition row (3d) metals in the 2+/3+ oxidation states, in connection with various solvation approaches. In terms of model definition, we investigate the construction of a minimally explicitly hydrated quantum cluster with a first and second hydration layer. We report on the convergence with respect to the coupled cluster expansion and the PNO space, as well as on the role of perturbative triple excitations. A recent implementation of the conductor-like polarizable continuum model (CPCM) for the DLPNO-CC approach is employed to determine self-consistent redox potentials at the coupled cluster level. Our results establish conditions for the convergence of DLPNO-CCSD(T) energetics and stress the absolute necessity to explicitly consider the second solvation sphere even when CPCM is used. The achievable accuracy for redox potentials of a practical DLPNO-based approach is, on average, 0.13 V. Furthermore, multilayer approaches that combine a higher-level DLPNO-CCSD(T) description of the first solvation sphere with a lower-level description of the second solvation layer are investigated. The present work establishes optimal and transferable methodological choices for employing DLPNO-based coupled cluster theory, the associated CPCM implementation, and cost-efficient multilayer derivatives of the approach for open-shell transition metal systems in complex environments.

Mechanisms of Electron-Induced Chemistry in Molecular Ices

Electron-induced chemistry is relevant to many processes that occur when ionizing radiation interacts with matter. This includes radiation damage, curing of polymers, and nanofabrication processes but also the formation of complex molecules in molecular ices grown on dust particles in space. High-energy radiation liberates from such materials an abundance of secondary electrons of which most have energies below 20 eV. These electrons efficiently trigger reactions when they attach to molecules or induce electronic excitation and further ionization. This review focuses on the present state of insight regarding the mechanisms of reactions induced by electrons with energies between 0 and 20 eV that lead to formation of larger products in binary ice layers consisting of small molecules (H2O, CO, CH3OH, NH3, CH4, C2H4, CH3CN, C2H6) or some derivatives thereof (C2H5NH2 and (C2H5)2NH, CH2=CHCH3). It summarizes our approach to identify products and quantify their amounts based on thermal desorption spectrometry (TDS) and electron-stimulated desorption (ESD) experiments performed in ultrahigh vacuum (UHV). The overview of the results demonstrates that, although the initial electron-molecule interaction is a non-thermal process, product formation from the resulting reactive species is often governed by subsequent reactions that follow well-known thermal and radical-driven mechanisms of organic chemistry.

Computational studies of the magneto-structural correlations in a manganese dimer with Jahn-Teller distortions

Inspired by the crystal structure of a Mn^III dinuclear complex we obtained featuring both Jahn-Teller (JT) elongation and compression distortions, we have modelled a series of complex cations based on the disordered crystal formulation; [Mn₂(L1)₂(μ₂-OH)₂)⁴⁺ (1), [Mn₂(L1)(L2)(μ₂-OH)₂)⁴⁺ (2), [Mn₂(L2)(L1)(μ₂-OH)₂)⁴⁺ (3), and [Mn₂(L2)₂(μ₂-OH)₂)⁴⁺ (4) (where L1 = (1E,1'E)-5-tert-butyl-3-(((4-(((5-tert-butyl-2-hydroxy-3-((E)-(hydroxyimino)methyl)benzyl)(methyl)amino)methyl)benzyl)(methyl)amino)methyl)-2-hydroxybenzaldehyde and L2 = 3,3'-(1,4-phenylenebis(methylene))bis(methylazanediyl)bis(methylene)bis(5-tert-butyl-2-hydroxybenzaldehyde)) with different geometries to investigate the effects of the distortions on the magnetic coupling parameter. All computationally modelled dimers had a ferromagnetic interaction between the Mn^III centres, with greater magnetic coupling calculated for complexes with both JT elongation and compression present. The ferromagnetic contribution to the J coupling was ascribed to the orthogonality of the singly occupied magnetic orbitals along with the cross-interaction between the unfilled Mn1(d_x²-y²) and singly occupied Mn2(d_x²-y²) orbitals. Constrained calculations showed that reducing the extent of the compression at Mn2 results in a concomitant increase in the dihedral angle between the JT axes, thereby reducing the relative magnitude of the magnetic coupling between Mn^III centres.

Photostability of 2,6-diaminopurine and its 2'-deoxyriboside investigated by femtosecond transient absorption spectroscopy

Ultraviolet radiation (UVR) from the sun is essential for the prebiotic syntheses of nucleotides, but it can also induce photolesions such as the cyclobutane pyrimidine dimers (CPDs) to RNA or DNA oligonucleotide in prebiotic Earth. 2,6-Diaminopurine (26DAP) has been proposed to repair CPDs in high yield under prebiotic conditions and be a key component in enhancing the photostability of higher-order prebiotic DNA structures. However, its electronic relaxation pathways have not been studied, which is necessary to know whether 26DAP could have survived the intense UV fluxes of the prebiotic Earth. We investigate the electronic relaxation mechanism of both 26DAP and its 2'-deoxyribonucleoside (26DAP-d) in aqueous solution using steady-state and femtosecond transient absorption measurements that are complemented with electronic-structure calculations. The results demonstrate that both purine derivatives are significantly photostable to UVR. It is shown that upon excitation at 287 nm, the lowest energy ¹ππ* state is initially populated. The population then branches following two relaxation coordinates in the ¹ππ* potential energy surface, which are identified as the C2- and C6-relaxation coordinates. The population following the C6-coordinate internally converts to the ground state nonradiatively through a nearly barrierless conical intersection within 0.7 ps in 26DAP or within 1.1 ps in 26DAP-d. The population that follows the C2-relaxation coordinate decays back to the ground state by a combination of nonradiative internal conversion via a conical intersection and fluorescence emission from the ¹ππ* minimum in 43 ps and 1.8 ns for the N9 and N7 tautomers of 26DAP, respectively, or in 70 ps for 26DAP-d. Fluorescence quantum yields of 0.037 and 0.008 are determined for 26DAP and 26DAP-d, respectively. Collectively, it is demonstrated that most of the excited state population in 26DAP and 26DAP-d decays back to the ground state via both nonradiative and radiative relaxation pathways. This result lends support to the idea that 26DAP could have accumulated in large enough quantities during the prebiotic era to participate in the formation of prebiotic RNA or DNA oligomers and act as a key component in the protection of the prebiotic genetic alphabet.

Theoretical Calculations and Experiments on the Thermal Properties of Fluorinated Graphene and Its Effects on the Thermal Decomposition of Nitrate Esters

Fluorinated graphene contains F atoms with high levels of chemical activity, and the application of fluorinated graphene in energetic materials may greatly contribute to the progress of combustion reactions. However, there is a lack of research on the thermal properties of fluorinated graphene and its application on nitrate esters. In this paper, theoretical calculations and experiments were used to study the thermal properties of fluorinated graphene and its application on nitrate esters. The anaerobicity and poor thermal stability of fluorinated graphene were proved by ab initio molecular dynamics (AIMD) calculations and TG-DSC experiments. The ester weakening effect of fluorinated graphene on nitroglycerin was determined via wavefunction analysis, with the greater the fluorination degree, the stronger the ester weakening effect. The existence of fluorinated graphene can significantly increase the heat dissipation of the composites, which was concluded by TG-DSC experiments and TG-DSC-MS-FTIR. The research in this article provides an important reference for the application of fluorinated graphene in energetic materials.

On the Photostability of Cyanuric Acid and Its Candidature as a Prebiotic Nucleobase

Cyanuric acid is a triazine derivative that has been identified from reactions performed under prebiotic conditions and has been proposed as a prospective precursor of ancestral RNA. For cyanuric acid to have played a key role during the prebiotic era, it would have needed to survive the harsh electromagnetic radiation conditions reaching the Earth’s surface during prebiotic times (≥200 nm). Therefore, the photostability of cyanuric acid would have been crucial for its accumulation during the prebiotic era. To evaluate the putative photostability of cyanuric acid in water, in this contribution, we employed density functional theory (DFT) and its time-dependent variant (TD-DFT) including implicit and explicit solvent effects. The calculations predict that cyanuric acid has an absorption maximum at ca. 160 nm (7.73 eV), with the lowest-energy absorption band extending to ca. 200 nm in an aqueous solution and exhibiting negligible absorption at longer wavelengths. Excitation of cyanuric acid at 160 nm or longer wavelengths leads to the population of S_5,6 singlet states, which have ππ* character and large oscillator strengths (0.8). The population reaching the S_5,6 states is expected to internally convert to the S_1,2 states in an ultrafast time scale. The S_1,2 states, which have nπ* character, are predicted to access a conical intersection with the ground state in a nearly barrierless fashion (ca. ≤ 0.13 eV), thus efficiently returning the population to the ground state. Furthermore, based on calculated spin–orbit coupling elements of ca. 6 to 8 cm⁻¹, the calculations predict that intersystem crossing to the triplet manifold should play a minor role in the electronic relaxation of cyanuric acid. We have also calculated the vertical ionization energy of cyanuric acid at 8.2 eV, which predicts that direct one-photon ionization of cyanuric acid should occur at ca. 150 nm. Collectively, the quantum-chemical calculations predict that cyanuric acid would have been highly photostable under the solar radiation conditions reaching the Earth’s surface during the prebiotic era in an aqueous solution. Of relevance to the chemical origin of life and RNA-first theories, these observations lend support to the idea that cyanuric acid could have accumulated in large quantities during the prebiotic era and thus strengthens its candidature as a relevant prebiotic nucleobase.

The origin of enhanced [Formula: see text] production from photoionized CO2 clusters

CO2-rich planetary atmospheres are continuously exposed to ionising radiation driving major photochemical processes. In the Martian atmosphere, CO2 clusters are predicted to exist at high altitudes motivating a deeper understanding of their photochemistry. In this joint experimental-theoretical study, we investigate the photoreactions of CO2 clusters (≤2 nm) induced by soft X-ray ionisation. We observe dramatically enhanced production of [Formula: see text] from photoionized CO2 clusters compared to the case of the isolated molecule and identify two relevant reactions. Using quantum chemistry calculations and multi-coincidence mass spectrometry, we pinpoint the origin of this enhancement: A size-dependent structural transition of the clusters from a covalently bonded arrangement to a weakly bonded polyhedral geometry that activates an exothermic reaction producing [Formula: see text]. Our results unambiguously demonstrate that the photochemistry of small clusters/particles will likely have a strong influence on the ion balance in atmospheres.

Gas-phase thermochemistry of polycyclic aromatic hydrocarbons: an approach integrating the quantum chemistry composite scheme and reaction generator

We introduce a protocol aimed at predicting the accurate gas-phase enthalpies of formation of polycyclic aromatic hydrocarbons (PAHs). Automatic generation of a dataset of equilibrated chemical reactions preserving the number of carbon atoms in each hybridization state on each side of equations is at the core of our scheme. The performed tests suggest the recommended enthalpy of formation to be derived via a two-step scheme. First, we consider the reactions with a minimal sum of the total number of particles involved, N, and the absolute difference between the total number of products and reactants, |ΔN|. Second, among these reactions, we identify the one with the smallest absolute reaction enthalpy change, . This approach has been applied to predict the gas-phase enthalpies of formation of 113 PAHs via the Feller-Peterson-Dixon approach. Our calculated values provide the mean absolute deviations of 1.7, 1.9, 4.2, 8.1, and 18.5 kJ mol^-1 with respect to the literature group-based error corrected (GBEC) G3MP2B3, ATOMIC (HC), group equivalent M06-2X, GBEC B3LYP, and G4MP2 values. Our predicted values give the mean signed and mean absolute errors of -7.5 and 12.9 kJ mol^-1 with respect to the experimental enthalpies of formation. The combination of our predicted and the experimental values provide the solid-state enthalpies of formation, , which are not available for a few species. Approaching these values as well as , producing large discrepancies from the experimental side, would be indispensable for testing and further tuning of computational chemistry approaches.

Adsorption and catalytic degradation of preservative parabens by graphene-family nanomaterials

Parabens pose increasing threats to human health due to endocrine disruption activity. Adsorption and degradation of parabens by three types of graphene-family nanomaterials (GFNs) were therefore investigated. For a given paraben, the maximum adsorption capacities (Q⁰) followed the order of reduced graphene oxide (RGO) > multilayered graphene (MG) > graphene oxide (GO); for a given GFN, Q⁰ followed the order of butylparaben (BuP) > propylparaben (PrP) > ethylparaben (EtP) > methylparaben (MeP), dominated by hydrophobic interaction. MeP removal by all the three GFNs was highly enhanced (0.55-4.37 times) with the assistance of H₂O₂ due to additional catalytic degradation process, and MG showed the highest removal enhancement. ∙OH was confirmed as the dominant radicals responsible for parabens degradation. For MG and RGO, the metal impurities (Fe, Cu, Mn, and Co) initiated Fenton-like reaction with H₂O₂ to generate ∙OH. GO contained oxygen-centered free radicals, which were responsible for ∙OH formation via transferring electron to H₂O₂. Four degradation byproducts of MeP were identified, including oxalic, propanedioic, fumaric, and 2,5-dihydroxybenzoic acids. Combined with density function theory calculations, the degradation sites and pathways were identified and confirmed. These findings provide useful information on mechanistic understanding towards the adsorption and degradation of parabens by GFNs.

Mechanistic Investigation of H2 O2 -dependent Chemiluminescence from Tetrabromo-1,4-Benzoquinone

As a H₂ O₂ -dependent bioluminescent substrate, tetrabromo-1,4-benzoquinone (TBBQ) was first isolated from acorn worm. The mechanism of chemiluminescence (CL) corresponding to the bioluminescence (BL) of acorn worm is largely unknown, let alone the mechanism of BL. In this article, we firstly studied the chemical and physical processes, and mechanism of H₂ O₂ -dependent CL from TBBQ by theoretical and experimental methods. The research results indicate: the CL process is initiated by a nucleophilic substitution reaction, which leads to the formation of an anionic dioxetane through five consecutive reactions; the anionic dioxetane decomposes to the first singlet excited state (S₁ ) via a conical interaction of the potential energy surfaces (PESs) between the ground (S₀ ) and S₁ state; the anionic S₁ -state changes to its neutral form by a proton transfer from the solvent and this neutral product is assigned as the actual luminophore. Moreover, the experimental detection of CL, ^. OH and the identifications of 2,3-dibromo maleic acid and 2-bromo malonic acid as the major final products provide direct evidence of the theoretically suggested mechanism. Finally, this study proves that the activity of the H₂ O₂ -dependent CL from TBBQ is significantly lower than the one from tetrachloro-1,4-benzoquinone (TCBQ), which is caused by the weaker electron withdrawing effect and the stronger heavy atomic effect of bromine.

Trigonally Distorted Hexacoordinate Co(II) Single-Ion Magnets

By simple reactions involving various cobalt(II) carboxylates (acetate and in situ prepared pivalate and 4-hydroxybenzoate salts) and neocuproine (neo), we were able to prepare three different carboxylate complexes with the general formula [Co(neo)(RCOO)₂] (R = –CH₃ for 1, (CH₃)₃C– for 2, and 4OH-C₄H₆– for 3). The [Co(neo)(RCOO)₂] molecules in the crystal structures of 1–3 adopt a rather distorted coordination environment, with the largest trigonal distortion observed for 1, whereas 2 and 3 are similarly distorted from ideal octahedral geometry. The combined theoretical and experimental investigations of magnetic properties revealed that the spin Hamiltonian formalism was not a valid approach and the L-S Hamiltonian had to be used to reveal very large magnetic anisotropies for 1–3. The measurements of AC susceptibility showed that all three compounds exhibited slow-relaxation of magnetization in a weak external static magnetic field, and thus can be classified as field-induced single-ion magnets. It is noteworthy that 1 also exhibits a weak AC signal in a zero-external magnetic field.

The Adsorption of Small Molecules on the Copper Paddle-Wheel: Influence of the Multi-Reference Ground State

We report a theoretical study of the adsorption of a set of small molecules (C2H2, CO, CO2, O2, H2O, CH3OH, C2H5OH) on the metal centers of the "copper paddle-wheel"-a key structural motif of many MOFs. A systematic comparison between DFT of different rungs, single-reference post-HF methods (MP2, SOS-MP2, MP3, DLPNO-CCSD(T)), and multi-reference approaches (CASSCF, DCD-CAS(2), NEVPT2) is performed in order to find a methodology that correctly describes the complicated electronic structure of paddle-wheel structure together with a reasonable description of non-covalent interactions. Apart from comparison with literature data (experimental values wherever possible), benchmark calculations with DLPNO-MR-CCSD were also performed. Despite tested methods show qualitative agreement in the majority of cases, we showed and discussed reasons for quantitative differences as well as more fundamental problems of specific cases.

From the Physicochemical Characteristic of Novel Hesperetin Hydrazone to Its In Vitro Antimicrobial Aspects

Microorganisms are able to give rise to biofilm formation on food matrixes and along food industry infrastructures or medical equipment. This growth may be reduced by the application of molecules preventing bacterial adhesion on these surfaces. A new Schiff base ligand, derivative of hesperetin, HABH (2-amino-N'-(2,3-dihydro-5,7-dihydroxy-2-(3-hydroxy-4-methoxyphenyl)chromen-4-ylidene)benzohydrazide), and its copper complex, CuHABH [CuLH₂(OAc)], were designed, synthesized and analyzed in terms of their structure and physicochemical properties, and tested as antibacterial agents. Their structures both in a solid state and in solution were established using several methods: FT-IR, ¹H NMR, ¹³C NMR, UV-Vis, FAB MS, EPR, ESI-MS and potentiometry. Coordination binding of the copper(II) complex dominating at the physiological pH region in the solution was found to be the same as that detected in the solid state. Furthermore, the interaction between the HABH and CuHABH with calf-thymus DNA (CT-DNA) were investigated. These interactions were tracked by UV-Vis, CD (circular dichroism) and spectrofluorimetry. The results indicate a stronger interaction of the CuHABH with the CT-DNA than the HABH. It can be assumed that the nature of the interactions is of the intercalating type, but in the high concentration range, the complex can bind to the DNA externally to phosphate residues or to a minor/major groove. The prepared compounds possess antibacterial and antibiofilm activities against Gram-positive and Gram-negative bacteria. Their antagonistic activity depends on the factor-strain test system. The glass was selected as a model surface for the experiments on antibiofilm activity. The adhesion of bacterial cells to the glass surface in the presence of the compounds was traced by luminometry and the best antiadhesive action against both bacterial strains was detected for the CuHABH complex. This molecule may play a crucial role in disrupting exopolymers (DNA/proteins) in biofilm formation and can be used to prevent bacterial adhesion especially on glass equipment.

Increased Photostability of the Integral mRNA Vaccine Component N1 -Methylpseudouridine Compared to Uridine

N1 -Methylation of pseudouridine (m1 ψ) replaces uridine (Urd) in several therapeutics, including the Moderna and BioNTech-Pfizer COVID-19 vaccines. Importantly, however, it is currently unknown if exposure to electromagnetic radiation can affect the chemical integrity and intrinsic stability of m1 ψ. In this study, the photochemistry of m1 ψ is compared to that of uridine by using photoirradiation at 267 nm, steady-state spectroscopy, and quantum-chemical calculations. Furthermore, femtosecond transient absorption measurements are collected to delineate the electronic relaxation mechanisms for both nucleosides under physiologically relevant conditions. It is shown that m1 ψ exhibits a 12-fold longer 1 ππ* decay lifetime than uridine and a 5-fold higher fluorescence yield. Notably, however, the experimental results also demonstrate that most of the excited state population in both molecules decays back to the ground state in an ultrafast time scale and that m1 ψ is 6.7-fold more photostable than Urd following irradiation at 267 nm.

Multi-fidelity prediction of molecular optical peaks with deep learning

Optical properties are central to molecular design for many applications, including solar cells and biomedical imaging. A variety of ab initio and statistical methods have been developed for their prediction, each with a trade-off between accuracy, generality, and cost. Existing theoretical methods such as time-dependent density functional theory (TD-DFT) are generalizable across chemical space because of their robust physics-based foundations but still exhibit random and systematic errors with respect to experiment despite their high computational cost. Statistical methods can achieve high accuracy at a lower cost, but data sparsity and unoptimized molecule and solvent representations often limit their ability to generalize. Here, we utilize directed message passing neural networks (D-MPNNs) to represent both dye molecules and solvents for predictions of molecular absorption peaks in solution. Additionally, we demonstrate a multi-fidelity approach based on an auxiliary model trained on over 28 000 TD-DFT calculations that further improves accuracy and generalizability, as shown through rigorous splitting strategies. Combining several openly-available experimental datasets, we benchmark these methods against a state-of-the-art regression tree algorithm and compare the D-MPNN solvent representation to several alternatives. Finally, we explore the interpretability of the learned representations using dimensionality reduction and evaluate the use of ensemble variance as an estimator of the epistemic uncertainty in our predictions of molecular peak absorption in solution. The prediction methods proposed herein can be integrated with active learning, generative modeling, and experimental workflows to enable the more rapid design of molecules with targeted optical properties.

Electrolyte Additives for Improving the High-Temperature Storage Performance of Li-Ion Battery NCM523∥Graphite with Overcharge Protection

The overcharge safety performance of lithium-ion batteries has been the major bottleneck in the widespread deployment of this promising technology. Pushing the limitations further may jeopardize cell safety when it is performed at high-temperature storage. On the basis of the lacking systematic research on overcharge protection electrolyte additives with high-temperature storage capacity, we explore the promotion effect of overcharge additives on electrolyte decomposition at 60 °C. Specifically, the addition of tris(trimethylsily) phosphite (TMSP) and lithium difluoro(oxalato)borate (LiDFOB) in the electrolyte can not only form the robust cathode electrolyte interface/solid electrolyte interphase (CEI/SEI) but also improve the thermal stability of the electrolyte. Therefore, we promote the electrolyte system to realize the 18,650 LIB storage at 60 °C for 50 days by optimizing the formula in the electrolyte containing biphenyl (BP) and cyclohexylbenzene (CHB) overcharge protection additives, and the capacity retention rate can reach more than 90% with overcharge safety. Further, the optimized electrolyte system has also been implemented to commercial 18,650 LIBs and demonstrates the widening of the route to the widespread application of the electrolyte under extreme conditions.