A new mixing of Hartree–Fock and local density‐functional theories

Rational design of molecularly imprinted electrochemical sensor based on Nb2C-MWCNTs heterostructures for highly sensitive and selective detection of Ochratoxin a

Developing an efficient method for screening Ochratoxin A (OTA) in agriculture products is vital to ensure food safety and human health. However, the complex food matrix seriously affects the sensitivity and accuracy. To address this issue, we designed a novel molecularly imprinted polymer (MIP) electrochemical sensor based on multiwalled carbon nanotube-modified niobium carbide (Nb2C-MWCNTs) with the aid of the density functional theory (DFT). In this design, a glassy carbon electrode (GCE) was first modified by Nb2C-MWCNTs heterostructure. Afterward, the MIP layer was prepared, with ortho-toluidine as a functional monomer selected via DFT and OTA acting as a template on the surface of Nb2C-MWCNTs/GCE using in-situ electropolymerization. Electrochemical tests and physical characterization revealed that Nb2C-MWCNTs improved the sensor's active surface area and electron transmission capacity. Nb2C-MWCNTs had a good synergistic effect on MIP, endowing the sensor with high sensitivity and specific recognition of OTA in complex food matrix systems. The MIP sensor showed a wide linear range from 0.04 to 10.0 μM with a limit of detection (LOD) of 3.6 nM. Moreover, it presented good repeatability and stability for its highly antifouling effect on OTA. In real sample analysis, the recoveries, ranging from 89.77% to 103.70%, agreed well with the results obtained by HPLC methods, suggesting the sensor has good accuracy and high potential in practical applications.

Vibrational circular dichroism of adenosine crystals

Adenosine is one of the building blocks of nucleic acids and other biologically important molecules. Spectroscopic methods have been among the most utilized techniques to study adenosine and its derivatives. However, most of them deal with adenosine in solution. Here, we present the first vibrational circular dichroism (VCD) spectroscopic study of adenosine crystals in solid state. Highly regular arrangement of adenosine molecules in a crystal resulted in a strongly enhanced supramolecular VCD signal originating from long-range coupling of vibrations. The data suggested that adenosine crystals, in contrast to guanosine ones, do not imbibe atmospheric water. Relatively large dimensions of the adenosine crystals resulted in scattering and substantial orientational artifacts affecting the spectra. Several strategies for tackling the artifacts have been proposed and tested. Atypical features in IR absorption spectra of crystalline adenosine (e.g., extremely low absorption in mid-IR spectral range) were observed and attributed to refractive properties of adenosine crystals.

Machine learning for accuracy in density functional approximations

Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. Here, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.

Theoretical Investigation on Water-Free, Water- and Self-Assisted H-Abstraction Reactions from Dimethylamine by Hydroxy Radicals

Accurate branching ratios of the H-abstraction reactions from dimethylamine (DMA) by OH radicals are important in understanding the atmospheric fate of DMA. In this work, the reaction kinetics of the water-free, water-assisted, and self-assisted H-abstraction reactions between DMA and OH radicals are accurately determined using the multipath canonical variational theory with the small-curvature tunneling correction, to explore the catalytic effects of the reactant (DMA) and product (water). To choose a suitable method that well describes the current reaction systems, various combinations with seven DFT methods and six basis sets are first evaluated, and the M08-HX/ma-TZVP method is identified as the most appropriate, with a mean unsigned deviation of 0.9 kcal mol-1 against the gold-standard CCSD(T)/CBS(T-Q) method. Based on the determined potential energy surfaces with the considerations of ground-state structures and specific-reaction parameters of zero-point energies, rate constants and branching ratios are calculated in a wide temperature range. The calculations show that the participation of water and DMA can lead to three-body complexes with a lower energy and influence the energy barriers, but neither of them shows the catalytic effect on the H-abstraction reactions in terms of kinetics. Additionally, the branching ratio analysis demonstrates that the product distribution is significantly altered in the presence of DMA and water.

Impact of the Acetyl Group on Cysteine: A Study of N-Acetyl-Cysteine through Rotational Spectroscopy

Herein, we report a spectroscopic study of N-acetyl-L-cysteine, an important antioxidant drug, using Fourier-transform microwave techniques and in isolated conditions. Two conformers are observed, where most stable structure adopts a cis disposition, and the second conformer has a lower abundance and adopts a trans disposition. The rotational constants and the barriers to methyl internal rotation are determined for each conformer, allowing a precise conformation identification. The results show that the cis form adopts an identical structure in the crystal, solution, and gas phases. Additionally, the structures are contrasted against those of cysteine.

Hey ho, where'd the proton go? Final deprotonation of O6 within the S3 state of photosystem II

The deprotonation of O6 within the S3 state marks the final deprotonation event before the formation of oxygen‑oxygen bond interactions and eventual production and release of dioxygen. Gaining a thorough understanding of this event, from the proton acceptors involved, to the exfiltration pathways available, is key in determining the nature of the resulting oxygen species, influencing the mechanism through which the first oxygen‑oxygen bond forms. Computational analysis, using BS-DFT methodologies, showed that proton abstraction by the local Glu189 residue provides consistent evidence against this being a viable mechanistic pathway due to the lack of a stable product structure. In contrast, abstraction via W3 shows an increasingly stable oxo-oxo product state between r[O5O6] = 2.1 Å & 1.9 Å. The resulting oxo-oxo state is stabilised through donation of β electron character from O6 to Mn1 and α electron character from O6 to O5. This donation from the O6 lone pair is shown to be a key factor in stabilising the oxo-oxo state, in addition to showing the initiation of first O5-O6 bond.

Benchmarking Quantum Mechanical Levels of Theory for Valence Parametrization in Force Fields

A wide range of density functional methods and basis sets are available to derive the electronic structure and properties of molecules. Quantum mechanical calculations are too computationally intensive for routine simulation of molecules in the condensed phase, prompting the development of computationally efficient force fields based on quantum mechanical data. Parametrizing general force fields, which cover a vast chemical space, necessitates the generation of sizable quantum mechanical data sets with optimized geometries and torsion scans. To achieve this efficiently, choosing a quantum mechanical method that balances computational cost and accuracy is crucial. In this study, we seek to assess the accuracy of quantum mechanical theory for specific properties such as conformer energies and torsion energetics. To comprehensively evaluate various methods, we focus on a representative set of 59 diverse small molecules, comparing approximately 25 combinations of functional and basis sets against the reference level coupled cluster calculations at the complete basis set limit.

Deciphering molecular structures: NMR spectroscopy and quantum mechanical insights of halogenated 4H-Chromenediones

Sesquiterpene lactones (SL) represent a class of secondary metabolites found in the Asteraceae family, notable for their unique structures. The SL α-santonin (1) and its derivatives are worthy of mention due to their diverse biological properties. Additionally, 4H-chromenes and 4H-chromones are appealing frameworks holding the capability to be used as structural motifs for new drugs. Furthermore, unambiguous structural elucidation is crucial for developing novel compounds for diverse applications. In this context, it is common to find in the literature molecules erroneously assigned. Therefore, the use of quantum mechanical calculations to simulate NMR chemical shifts has emerged as a valuable strategy. In this work, we conceived the synthesis of two halogenated 4H-chromenediones derived from photosantonic acid (2), a photoproduct arising from irradiation of α-santonin (1) in the ultraviolet region. The structure of the chlorinated and brominated products was determined by NMR analysis, with the aid of quantum mechanical calculations at the B3LYP/6-311 + G(2d,p)//M062x/6-31 + G(d,p) level of theory. All analyses were in agreement and led to the assignment of the brominated 4H-chromene-2,7-dione as (3S,3aS,5aR,9bS)-5a-(2-bromopropan-2-yl)-3-methyl-3,3a,5,5a,8,9b-hexahydro-4H-furo[2,3-f]chromene-2,7-dione (11b) and of the chlorinated 4H-chromene-2,7-dione as (3S,3aS,5aR,9bS)-5a-(2-chloropropan-2-yl)-3-methyl-3,3a,5,5a,8,9b-hexahydro-4H-furo[2,3-f]chromene-2,7-dione (12b). The diastereoselectivities of the reactions were explained based on products and intermediates formation energy calculated using B3LYP/6-31 + G(d,p) as the level of theory. Structures 11b and 12b were identified as the thermodynamic and kinetic products of the reaction among all candidates. Consequently, the strategy utilized in this study is robust and successfully illustrates the use of quantum mechanical calculations in the structural elucidation of new compounds with potential applications as novel drugs or products.

Detrimental Increase of Spin-Phonon Coupling in Molecular Qubits on Substrates

Molecular qubits are a promising platform for quantum information systems. Although single molecule and ensemble studies have assessed the performance of S = 1/2 molecules, it is understood that to function in devices, regular arrays of addressable qubits supported by a substrate are needed. The substrate imposes mechanical and electronic boundary conditions on the molecule; however, the impact of these effects on spin-lattice relaxation times is not well understood. Here we perform electronic structure calculations to assess the effects of a graphene (Cgr) substrate on the molecular qubit copper phthalocyanine (CuPc). We use a progressive Hessian approach to efficiently calculate and separate the substrate contributions. We also use a simple thermal model to predict the impact of these changes on the spin-phonon coupling from 0 to 200 K. Further analysis of the individual vibrational modes with and without Cgr shows that an overall increase in SPC between the vibrations modes of CuPc with the surface reduces the spin-lattice relaxation time T1. We explain these changes by examining how the substrate lifts symmetries of CuPc in the absorbed configuration. Our work shows that a surface can have a large unintentional impact on SPC and that ways to reduce this coupling need to be found to fully exploit arrays of molecular qubits in device architectures.

Synthesis of triphyrin(2.1.1)-triphyrin(2.1.1) homo dimers and bis-Re(I) triphyrin dimer complex

First examples of covalently linked triphyrin(2.1.1)-triphyrin(2.1.1) homo dimers connected via flexible alkoxy linkers at β-pyrrole positions were synthesized in 28-30% yields by reacting 2-bromo triphyrin(2.1.1) with three different alkane diols such as 1,2-ethanediol, 1,4-butanediol and 1,6-hexanediol in toluene under Pd2(dba)3/DPEphos catalyzed conditions by refluxing for 12-20 h. One of the free base triphyrin-triphyrin homo dimers was metallated by treating the homo dimer with Re(CO)5Cl in toluene under reflux and afforded bis-Re(I) homo dimer complex in 53% yield. Three free base homo dimers and one bis-Re(I) homo dimer complex were thoroughly characterized and studied using HR-MS, 1D and 2D NMR, absorption spectroscopy, cyclic voltammetry and DFT calculations. Our attempts to obtain the X-ray structure of homo dimers were not successful. However, DFT studies revealed that homo dimers prefer to exist in anti-conformation rather than syn-conformation and both triphyrin macrocycles in homo dimers were similar in terms of the deviation of pyrrole rings from the mean plane defined by the four meso carbons of the respective triphyrin unit. Absorption spectra of homo dimers showed similar absorption features such as in triphyrin(2.1.1) with slight changes in their absorption peak maxima. Electrochemical studies revealed that the homo dimers were electron rich, whereas the bis-Re(I) complex was electron deficient in nature.

A new parameterization of the DFT/CIS method with applications to core-level spectroscopy

Time-dependent density functional theory (TD-DFT) within a restricted excitation space is an efficient means to compute core-level excitation energies using only a small subset of the occupied orbitals. However, core-to-valence excitation energies are significantly underestimated when standard exchange-correlation functionals are used, which is partly traceable to systemic issues with TD-DFT's description of Rydberg and charge-transfer excited states. To mitigate this, we have implemented an empirically modified combination of configuration interaction with single substitutions (CIS) based on Kohn-Sham orbitals, which is known as "DFT/CIS." This semi-empirical approach is well-suited for simulating x-ray near-edge spectra, as it contains sufficient exact exchange to model charge-transfer excitations yet retains DFT's low-cost description of dynamical electron correlation. Empirical corrections to the matrix elements enable semi-quantitative simulation of near-edge x-ray spectra without the need for significant a posteriori shifts; this should be useful in complex molecules and materials with multiple overlapping x-ray edges. Parameter optimization for use with a specific range-separated hybrid functional makes this a black-box method intended for both core and valence spectroscopy. Results herein demonstrate that realistic K-edge absorption and emission spectra can be obtained for second- and third-row elements and 3d transition metals, with promising results for L-edge spectra as well. DFT/CIS calculations require absolute shifts that are considerably smaller than what is typical in TD-DFT.

Asymmetric photoenzymatic incorporation of fluorinated motifs into olefins

Enzymes capable of assimilating fluorinated feedstocks are scarce. This situation poses a challenge for the biosynthesis of fluorinated compounds used in pharmaceuticals, agrochemicals, and materials. We developed a photoenzymatic hydrofluoroalkylation that integrates fluorinated motifs into olefins. The photoinduced promiscuity of flavin-dependent ene-reductases enables the generation of carbon-centered radicals from iodinated fluoroalkanes, which are directed by the photoenzyme to engage enantioselectively with olefins. This approach facilitates stereocontrol through interaction between a singular fluorinated unit and the enzyme, securing high enantioselectivity at β, γ, or δ positions of fluorinated groups through enzymatic hydrogen atom transfer-a process that is notably challenging with conventional chemocatalysis. This work advances enzymatic strategies for integrating fluorinated chemical feedstocks and opens avenues for asymmetric synthesis of fluorinated compounds.

Rotational Spectroscopy as a Tool to Characterize Sweet Taste: The Study of Dulcin

According to old theories of sweetness, the perception of sweet substances is closely linked to the arrangement of atoms within them. To assess the validity of these theories, we conducted an analysis of the structure of the artificial sweetener dulcin for the first time, utilizing microwave spectroscopy and a laser ablation source. These techniques have enabled the identification of two conformers, which are stabilized by an intramolecular hydrogen bond between the amino group and the phenyl ring. The observed conformations were examined in light of the Shallenberger-Acree-Kier molecular theory of sweet taste, and they align with the hypothesized criteria. Furthermore, the study illustrates how conformational relaxation can alter the equilibrium conformational distribution, resulting in the absence of certain conformers in the conformational landscape.

Unraveling the Isotropic Hyperfine Coupling Constants of Nitroxide Radicals via Molecular Electrostatic Potential Analysis

Nitroxide radicals have wide and promising applications as organic magnetic materials. Modulating the isotropic hyperfine coupling constants (HFCCs) of these radicals through proper structural design is an effective strategy for their application as spin probes and spin labels. In the present work, density functional theory calculations were carried out to develop a robust descriptor based on the molecular electrostatic potential for nitrogen HFCCs of nitroxide radicals. Forty nitroxide radicals from five distinct classes, namely, derivatives of cyclic, acyclic, imino, nitronyl, and benzimidazole nitronyl nitroxides, were selected, and the molecular electrostatic potential (MESP) at the nitrogen atom (VN) of the NO moiety was calculated. The VN values efficiently capture the electronic changes associated with the steric and electronic nature of these systems. A significant correlation between VN values and the experimental HFCCs of nitrogen nuclei demonstrates the applicability of VN as a simple and efficient descriptor for monitoring HFCCs. Furthermore, a good correlation between VN and experimental nitrogen HFCCs for each class of nitroxide radicals indicates the use of VN in the evaluation of the magnetic nature of the nitroxide radicals. The findings in this work are expected to facilitate the design of novel nitroxide radicals with desirable magnetic properties based on MESP topology analysis.

Assessment of the applicability of DFT methods to [Cp*Rh]-catalyzed hydrogen evolution processes

The present computational study provides a benchmark of density functional theory (DFT) methods in describing hydrogen evolution processes catalyzed by [Cp*Rh]-containing organometallic complexes. A test set was composed of 26 elementary reactions featuring chemical transformations and bonding situations essential for the field, including the emerging concept of non-innocent Cp* behavior. Reference values were obtained from a highly accurate 3/4 complete basis set and 6/7 complete PNO space extrapolated DLPNO-CCSD(T) energies. The performance of lower-level extrapolation procedures was also assessed. We considered 84 density functionals (DF) (including 13 generalized gradient approximations (GGA), nine meta-GGAs, 33 hybrids, and 29 double-hybrids) and three composite methods (HF-3c, PBEh-3c, and r2SCAN-3c), combined with different types of dispersion corrections (D3(0), D3BJ, D4, and VV10). The most accurate approach is the PBE0-DH-D3BJ (MAD of 1.36 kcal mol-1) followed by TPSS0-D3BJ (MAD of 1.60 kcal mol-1). Low-cost r2SCAN-3c composite provides a less accurate but much faster alternative (MAD of 2.39 kcal mol-1). The widely used Minnesota-family M06-L, M06, and M06-2X DFs should be avoided (MADs of 3.70, 3.94, and 4.01 kcal mol-1, respectively).

Cyano-Polycyclic Aromatic Hydrocarbon Interstellar Candidates: Laboratory Identification, Equilibrium Structure and Astronomical Search of Cyanobiphenylene

The interplay between laboratory rotational spectroscopy and radio astronomical observations provides the most effective procedure for identifying molecules in the interstellar medium (ISM). Following the recent interstellar detections of several Polycyclic Aromatic Hydrocarbons (PAHs) and cyano derivatives in the dense molecular cloud TMC-1, it is reasonable to consider searching for other cyano-PAHs in this astronomical source. We present a rotational spectroscopy investigation of the two cyano derivatives of the PAH biphenylene, a plausible reaction product of interstellar benzyne. The rotational spectrum provided molecular parameters for the parent species and 14 monosubstituted isotopologues for each isomer. An accurate equilibrium structure was determined for both isomers using Watson's mass-dependence method (rm(2)), offering information on its uncommon ring union. Astronomical searches for the cyanobiphenylene isomers have been undertaken in TMC-1, using the QUIJOTE line survey. No lines of any isomer were found in this astronomical source, but the experimental data will serve to enable future searches for these species in the ISM.

Examination of Nafion membrane performance with acid-base functionalized graphene oxide

Functionalized graphene oxide (GO) as a good additive can improve the performance of proton exchange membrane (PEM) via the introduction of various functional groups. How to balance the proton conductivity and durability of membrane based on functionalized GO is a key issue. In this work, benzoic-acid-functionalized GO(BAF-GO) and 1,2,4-triazole-functionalized GO(TF-GO) are employed as doping candidates, and the co-doping effect on membrane performance is investigated by means of experiment and molecular dynamics simulation. Meanwhile, the quantum chemistry method is implemented to explore the interaction between TF-GO, membrane and BAF-GO. The results reveal that the composite membrane exhibits high durability and enhanced proton conductivity. When the doping mass ratio of BAF-GO to TF-GO is 3 : 1, the proton conductivity can be greatly improved, especially under low-humidity conditions. Excessive addition of basic groups does not enhance proton transport.

A Vacuum Ultraviolet Photoionization Mass Spectrometry and Density Functional Calculation Study of Formic Acid-Water Clusters

The interaction between formic acid (FA) and water (W) holds significant importance in various chemical processes. Our study combines vacuum-ultraviolet photoionization mass spectrometry with density functional calculations to investigate formic acid water clusters generated in supersonic molecular beams. The mass spectra obtained reveal the formation of protonated clusters as the major product. Enhanced intensities are observed in the mass spectra for a number of clusters holding the following composition, FA1W5H+, FA2W4H+, FA3W3H+, FA4W2H+, FA5W1H+ and FA6W2H+ compared to their neighbors with one less or one more water component. Our calculations shed light on these potentially stable structures, highlighting cyclic arrangements with molecules enclosed within the ring as the most stable structures, and demonstrate a decrease in the stability upon the addition of a water molecule. Comparing experimental appearance energies with calculated ionization energies suggests that fragmentation can occur from clusters of various sizes.

Influence of Electronic Correlations on Electron-Phonon Interactions of Molecular Systems with the GW and Coupled Cluster Methods

Electron-phonon interactions are of great importance to a variety of physical phenomena, and their accurate description is an important goal for first-principles calculations. Isolated examples of materials and molecular systems have emerged where electron-phonon coupling is enhanced over density functional theory (DFT) when using the Green's-function-based ab initio GW method, which provides a more accurate description of electronic correlations. It is, however, unclear how general this enhancement is and how employing high-end quantum chemistry methods, which further improve the description of electronic correlations, might further alter electron-phonon interactions over GW or DFT. Here, we address these questions by computing the renormalization of the highest occupied molecular orbital energies of Thiel's set of organic molecules by harmonic vibrations using DFT, GW, and equation-of-motion coupled-cluster calculations. We find that, depending on the amount of exact exchange included in the DFT starting point, GW can increase the magnitude of the electron-phonon coupling across Thiel's set of molecules by an average factor of 1.1-1.8 compared to the underlying DFT, while equation-of-motion coupled-cluster leads to an increase of 1.4-2. The electron-phonon coupling predicted with the ab initio GW method is generally in much closer agreement to coupled cluster values compared to DFT, establishing GW as a promising route for accurately computing electron-phonon phenomena in molecules and beyond at a much lower computational cost than higher-end quantum chemistry techniques.

Silylation-Driven Volatility Enhancement in Mononuclear Calcium, Magnesium, and Barium Complexes with Monomethyl or Silyl Ether and Bis(diketonate)

Magnesium, calcium, and barium heteroleptic complexes were synthesized by the substitution reaction of the bis(trimethylsilyl)amide of Mg(btsa)2·DME, Ca(btsa)2·DME, and Ba(btsa)2·2DME with an ethereal group and hfac ligands (btsa = bis(trimethylsilyl)amide, DME = dimethoxyethane). The compounds Mg(dts)(hfac)2 (1), Ca(dts)(hfac)2 (2), Mg(dmts)(hfac)2 (3), Ca(dmts)(hfac)2 (4), and Ba(dmts)(hfac)2 (5) were fabricated and analyzed using various techniques, including Fourier transform infrared spectroscopy, nuclear magnetic resonance spectroscopy, thermogravimetric analyses, and elemental analysis (dts = 2,2-dimethyl-3,6,9-trioxa-2-siladecane, dmts = 2,2-dimethyl-3,6,9,12-tetraoxa-2-silatridecane, hfac = hexafluoroacetylacetonate). The structures of complexes 2, 4, and 5 were confirmed using single-crystal X-ray crystallography; all complexes display monomeric structures. All compounds underwent trimethylsilylation of the coordinating ethereal alcohols (meeH and tmgeH) in the presence of HMDS as byproducts because of their increasing acidity originating from the electron-withdrawing hfac ligands. (meeH = 2-(2-methoxyethoxy)ethan-1-ol, tmgeH = tri(ethylene glycol) monoethyl ether, HMDS = hexamethyldisilazane).

Experimental and Computational Anticorrosion Behaviors of Pyrazole s-Triazine/anilino-morpholino Derivatives for Steel in Acidic Solutions

The corrosion inhibition of C-steel by two s-triazine/morpholino-anilino-pyrazole derivatives, namely, 4-(3,5-dimethyl-1H-pyrazol-1-yl)-6-morpholino-N-phenyl-1,3,5-triazin-2-amine (1) and N-(4-bromophenyl)-4-(3,5-dimethyl-1H-pyrazol-1-yl)-6-morpholino-1,3,5-triazin-2-amine (2) was investigated by impedimetric and potentiometric studies. It was found that (1) and (2) acted as cathodic-type corrosion inhibitors that retard the hydrogen evolution reaction. The percent corrosion inhibition, 98.5% for compound (2) (with bromo substituent) at 80 ppm, was slightly higher than 97.8% for (1) at 100 ppm. Thus, the replacement of a -H with -Br substituent increased the corrosion inhibition properties. Compound (2) exhibited Temkin isotherm adsorption, whereas compound (1) exhibited Langmuir adsorption. Scanning electron microscopy (SEM) analysis of the steel surface indicated that the inhibitors caused protection of the surface. The weight loss experiment also proved the decrease in the corrosion rate when inhibitors were added. The difference in inhibitory efficiency between compounds (1) and (2) was investigated by density functional theory (DFT) to study neutral and protonated species in gaseous and aqueous phases. The theoretical analysis demonstrated that compound (2) exhibited higher inhibitory activity on a metal surface compared to compound (1), aligning with the experimental results. The energy associated with the metal/adsorbate arrangement, represented by dE ads/dNi , was higher for (2) (-380.91 kcal mol-1) compared to (1) (-371.64 kcal mol-1). This indicated better adsorption of (2) over (1).

Proton Conducting Neuromorphic Materials and Devices

Neuromorphic computing and artificial intelligence hardware generally aims to emulate features found in biological neural circuit components and to enable the development of energy-efficient machines. In the biological brain, ionic currents and temporal concentration gradients control information flow and storage. It is therefore of interest to examine materials and devices for neuromorphic computing wherein ionic and electronic currents can propagate. Protons being mobile under an external electric field offers a compelling avenue for facilitating biological functionalities in artificial synapses and neurons. In this review, we first highlight the interesting biological analog of protons as neurotransmitters in various animals. We then discuss the experimental approaches and mechanisms of proton doping in various classes of inorganic and organic proton-conducting materials for the advancement of neuromorphic architectures. Since hydrogen is among the lightest of elements, characterization in a solid matrix requires advanced techniques. We review powerful synchrotron-based spectroscopic techniques for characterizing hydrogen doping in various materials as well as complementary scattering techniques to detect hydrogen. First-principles calculations are then discussed as they help provide an understanding of proton migration and electronic structure modification. Outstanding scientific challenges to further our understanding of proton doping and its use in emerging neuromorphic electronics are pointed out.

Fixed-node diffusion Monte Carlo shows promise for modeling reaction thermochemistry of hydrocarbon-based radicals

Reliable thermodynamic and kinetic properties of free radical polymerization reactions are essential for synthesizing both primary polymeric materials and specialty polymers. The computational generation of these data from quantum chemistry requires a time-efficient method capable of capturing the essential physics. One such method, fixed-node diffusion Monte Carlo (FN-DMC) (using single Slater-Jastrow trial wavefunctions), has demonstrated the capability to recover 90%-95% of missing dynamic correlation energy for typical systems. In this study, methyl radical addition to ethylene serves as a simple model to test FN-DMC's ability to calculate enthalpies of reaction and activation energies with different time steps, antisymmetric trial wavefunctions, basis set sizes, and effective core potentials. The FN-DMC computational protocol thus defined for methyl radical addition to ethylene is subsequently benchmarked against Weizmann-1 and experimental reaction enthalpies from Lin et al.'s test set of 21 radical addition and 28 hydrogen abstraction enthalpies. Our findings reveal that FN-DMC consistently generates reaction enthalpies with chemical accuracy, exhibiting mean absolute deviation of 3.5(7) and 1.4(8) kJ/mol from the Weizmann-1 reference for radical addition and hydrogen abstraction reactions, respectively. Given its favorable computational scaling and high degree of parallelizability, we, therefore, recommend more comprehensive testing of FN-DMC with effective core potentials to address more extensive and intricate polymerization reactions and reactions with other radicals.

Enhancing Bioactivity through the Transfer of the 2-(Hydroxymethoxy)Vinyl Moiety: Application in the Modification of Tyrosol and Hinokitiol

Utilizing Density Functional Theory (DFT) calculations at the B3LYP/QZVP level and incorporating the Conductor-like Polarizable Continuum Model (C-PCM) for solvation, the thermodynamic and chemical activity properties of 21-(hydroxymethoxy)henicosadecaenal, identified in cultured freshwater pearls from the mollusk Hyriopsis cumingii, have been elucidated. The study demonstrates that this compound releases formaldehyde, a potent antimicrobial agent, through dehydrogenation and deprotonation processes in both hydrophilic and lipophilic environments. Moreover, this polyenal exhibits strong anti-reductant properties, effectively scavenging free radicals. These critical properties classify the pearl-derived ingredient as a natural multi-functional compound, serving as a coloring, antiradical, and antimicrobial agent. The 2-(hydroxymethoxy)vinyl (HMV) moiety responsible for the formaldehyde release can be transferred to other compounds, thereby enhancing their biological activity. For instance, tyrosol (4-(2-hydroxyethyl)phenol) can be modified by substituting the less active 2-hydroxyethyl group with the active HMV one, and hinokitiol (4-isopropylotropolone) can be functionalized by attaching this moiety to the tropolone ring. A new type of meso-carrier, structurally modeled on pearls, with active substances loaded both in the layers and the mineral part, has been proposed.

Quantifying non-transferrin-bound iron (NTBI) in human plasma: incorporating BODIPY-pyridylhydrazone (BODIPY-PH) within a thin green film linked to a portable fluorescence-based device

Free iron in human serum or non-transferrin-bound iron (NTBI) can generate free radicals and lead to oxidative damage. Moreover, it is highly toxic to various tissues and a vital biomarker related to the iron-loading status of thalassemia and Alzheimer's patients. In NTBI in healthy individuals, NTBI levels are typically less than 1 µM; current NTBI analysis usually requires advanced instrumentation and many-step sample pretreatment. To address this issue, we employed our invented BODIPY derivative, BODIPY-PH, as a fluorescence probe and trapped it onto the microcentrifuge tube lid using tapioca starch. The fluorescence intensity of BODIPY-PH increased with increasing NTBI concentration (turn-on). The developed portable reaction chamber facilitates rapid analysis (∼5 min) using small sample volumes (10 μL sample in a total volume of 600 μL). Under optimum conditions, using the sample-developed portable fluorescence device and fluorescence spectrometer, we achieved impressive limits of detection (LOD) of 0.003 and 0.0015 μM, respectively. Furthermore, the developed sensors show relatively high selectivity toward Fe3+ over other metal ions and biomolecules (i.e., Fe2+, Cr3+, Cu2+, and glucose). The sensor performance in serum samples of thalassemia patients exhibited no significant difference compared to the labeled value (obtained from standard methods). Overall, the developed fluorescence sensor is suitable for determining NTBI and offers high sensitivity, high selectivity, and a short incubation time (5 min). Moreover, the method requires a limited number of reagents, is simple to use, and uses low-cost equipment to determine NTBI in human serum samples.

Identifying and embedding transferability in data-driven representations of chemical space

Transferability, especially in the context of model generalization, is a paradigm of all scientific disciplines. However, the rapid advancement of machine learned model development threatens this paradigm, as it can be difficult to understand how transferability is embedded (or missed) in complex models developed using large training data sets. Two related open problems are how to identify, without relying on human intuition, what makes training data transferable; and how to embed transferability into training data. To solve both problems for ab initio chemical modelling, an indispensable tool in everyday chemistry research, we introduce a transferability assessment tool (TAT) and demonstrate it on a controllable data-driven model for developing density functional approximations (DFAs). We reveal that human intuition in the curation of training data introduces chemical biases that can hamper the transferability of data-driven DFAs. We use our TAT to motivate three transferability principles; one of which introduces the key concept of transferable diversity. Finally, we propose data curation strategies for general-purpose machine learning models in chemistry that identify and embed the transferability principles.

D4-Symmetric Dirhodium Tetrakis(binaphthylphosphate) Catalysts for Enantioselective Functionalization of Unactivated C-H Bonds

Dirhodium tetrakis(2,2'-binaphthylphosphate) catalysts were successfully developed for asymmetric C-H functionalization with trichloroethyl aryldiazoacetates as the carbene precursors. The 2,2'-binaphthylphosphate (BNP) ligands were modified by introduction of aryl and/or chloro functionality at the 4,4',6,6' positions. As the BNP ligands are C2-symmetric, the resulting dirhodium tetrakis(2,2'-binaphthylphosphate) complexes were expected to be D4-symmetric, but X-ray crystallographic and computational studies revealed this is not always the case because of internal T-shaped CH-π and aryl-aryl interactions between the ligands. The optimum catalyst is Rh2(S-megaBNP)4, with 3,5-di(tert-butyl)phenyl substituents at the 4,4' positions and chloro substituents at the 6,6' positions. This catalyst adopts a D4-symmetric arrangement and is ideally suited for site-selective C-H functionalization at unactivated tertiary sites with high levels of enantioselectivity, outperforming the best dirhodium tetracarboxylate catalyst developed for this reaction. The standard reactions were conducted with a catalyst loading of 1 mol % but lower catalyst loadings can be used if desired, as illustrated in the C-H functionalization of cyclohexane in 91% ee with 0.0025 mol % catalyst loading (29,400 turnover numbers). These studies further illustrate the effectiveness of donor/acceptor carbenes in site-selective intermolecular C-H functionalization and expand the toolbox of catalysts available for catalyst-controlled C-H functionalization.

Quantum chemical calculations of nitrosamine activation and deactivation pathways for carcinogenicity risk assessment

N-nitrosamines and nitrosamine drug substance related impurities (NDSRIs) became a critical topic for the development and safety of small molecule medicines following the withdrawal of various pharmaceutical products from the market. To assess the mutagenic and carcinogenic potential of different N-nitrosamines lacking robust carcinogenicity data, several approaches are in use including the published carcinogenic potency categorization approach (CPCA), the Enhanced Ames Test (EAT), in vivo mutagenicity studies as well as read-across to analogue molecules with robust carcinogenicity data. We employ quantum chemical calculations as a pivotal tool providing insights into the likelihood of reactive ion formation and subsequent DNA alkylation for a selection of molecules including e.g., carcinogenic N-nitrosopiperazine (NPZ), N-nitrosopiperidine (NPIP), together with N-nitrosodimethylamine (NDMA) as well as non-carcinogenic N-nitrosomethyl-tert-butylamine (NTBA) and bis (butan-2-yl) (nitros)amine (BBNA). In addition, a series of nitroso-methylaminopyridines is compared side-by-side. We draw comparisons between calculated reaction profiles for structures representing motifs common to NDSRIs and those of confirmed carcinogenic and non-carcinogenic molecules with in vivo data from cancer bioassays. Furthermore, our approach enables insights into reactivity and relative stability of intermediate species that can be formed upon activation of several nitrosamines. Most notably, we reveal consistent differences between the free energy profiles of carcinogenic and non-carcinogenic molecules. For the former, the intermediate diazonium ions mostly react, kinetically controlled, to the more stable DNA adducts and less to the water adducts via transition-states of similar heights. Non-carcinogenic molecules yield stable carbocations as intermediates that, thermodynamically controlled, more likely form the statistically preferred water adducts. In conclusion, our data confirm that quantum chemical calculations can contribute to a weight of evidence approach for the risk assessment of nitrosamines.

Photoinduced Charge Centralization Quenches the Fluorescence of Conjugation-Fused Tetrazine Labels with Red-to-Near-Infrared Emissions

Tetrazine-derived fluorogenic labels are extensively studied for their potential in biological and medical imaging. Nonetheless, the fluorescence quenching mechanism in numerous precursors continues to be debated, particularly as the wavelengths extend into the red and near-infrared (NIR) regions. This challenge poses obstacles to systematically optimizing their fluorogenicity, i.e., achieving red-shifted wavelengths and improved fluorescence turn-on signals through click reactions. This paper highlights the significance of photoinduced charge centralization (PCC), a quenching mechanism observed in tetrazine-fused fluorogenic labels with integrated π-conjugations. PCC is primarily responsible for the quenching effects observed in such labels emitting in the red-to-NIR spectrum. Drawing from structure-property relationships, this study proposes two molecular design strategies that incorporate the PCC mechanism and constitutional isomerization to develop high-performance tetrazine-based labels. These strategies facilitate multiplex fluorescence imaging following click reactions, promising significant advancements in bio-orthogonal imaging techniques.

Theoretical Study of Dissolved Gas Molecules in Transformer Oil Adsorbed on Intrinsic and TM (Ta, V)-Doped MoTe2 Monolayer

In this paper, CH4, C2H2, H2, and CO adsorbed on intrinsic MoTe2 monolayer and transition metal atom (Ta, V)-doped MoTe2 monolayer have been investigated with density functional theory based on first-principles study. The adsorption energy, geometries, band structures, and density of states of four gases (CH4, C2H2, H2, and CO) adsorbed on the MoTe2 and doped MoTe2 surfaces were analyzed. The results shown that the gas adsorption performance of transition metal atom (Ta, V)-doped MoTe2 monolayers is more superior than that of intrinsic MoTe2, and the adsorption energy and charge transfer of the adsorbed gases on the TM-MoTe2 monolayer are significantly increased in comparison with both sides. Among them, Ta-MoTe2 has the largest Eads value in the adsorbed CO system with a very small adsorption distance, as well as a more suitable recovery time of CO at room temperature, so Ta-MoTe2 can be a candidate material for CO detection. New atoms were introduced during the doping process, which increased the carrier density and carrier mobility of the material, thus improving the charge transfer at the surface of the material. which provides a direction for the gas-sensitive properties of metal Ta-modified MoTe2 materials.

Dual-quartet phosphorescent emission in the open-shell M1Ag13 (M = Pt, Pd) nanoclusters

Dual emission (DE) in nanoclusters (NCs) is considerably significant in the research and application of ratiometric sensing, bioimaging, and novel optoelectronic devices. Exploring the DE mechanism in open-shell NCs with doublet or quartet emissions remains challenging because synthesizing open-shell NCs is difficult due to their inherent instability. Here, we synthesize two dual-emissive M1Ag13(PFBT)6(TPP)7 (M = Pt, Pd; PFBT = pentafluorobenzenethiol; TPP = triphenylphosphine) NCs with a 7-electron open-shell configuration to reveal the DE mechanism. Both NCs comprise a crown-like M1Ag11 kernel with Pt or Pd in the center surrounded by five PPh3 ligands and two Ag(SR)3(PPh3) motifs. The combined experimental and theoretical studies revealed the origin of DE in Pt1Ag13 and Pd1Ag13. Specifically, the high-energy visible emission and the low-energy near-infrared emission arise from two distinct quartet excited states: the core-shell charge transfer and core-based states, respectively. Moreover, PFBT ligands are found to play an important role in the existence of DE, as its low-lying π* levels result in energetically accessible core-shell transitions. This novel report on the dual-quartet phosphorescent emission in NCs with an open-shell electronic configuration advances insights into the origin of dual-emissive NCs and promotes their potential application in magnetoluminescence and novel optoelectronic devices.

Paramagnetic Properties of [AnIV(NO3)6]2- Complexes (An = U, Np, Pu) Probed by NMR Spectroscopy and Quantum Chemical Calculations

Actinide +IV complexes with six nitrates [AnIV(NO3)6]2- (An = Th, U, Np, and Pu) have been studied by 15N and 17O NMR spectroscopy in solution and first-principles calculations. Magnetic susceptibilities were evaluated experimentally using the Evans method and are in good agreement with the ab initio values. The evolution in the series of the crystal field parameters deduced from ab initio calculations is discussed. The NMR paramagnetic shifts are analyzed based on ab initio calculations. Because the cubic symmetry of the complex quenches the dipolar contribution, they are only of Fermi contact origin. They are evaluated from first-principles based on a complete active space/density functional theory (DFT) strategy, in good accordance with the experimental one. The ligand hyperfine coupling constants are deduced from paramagnetic shifts and calculated using unrestricted DFT. The latter are decomposed in terms of the contribution of molecular orbitals. It highlights two pathways for the delocalization of the spin density from the metallic open-shell 5f orbitals to the NMR active nuclei, either through the valence 5f hybridized with 6d to the valence 2p molecular orbitals of the ligands, or by spin polarization of the metallic 6p orbitals which interact with the 2s-based molecular orbitals of the ligands.

Discovery of a terpene synthase synthesizing a nearly non-flexible eunicellane reveals the basis of flexibility

Eunicellane diterpenoids, containing a typical 6,10-bicycle, are bioactive compounds widely present in marine corals, but rarely found in bacteria and plants. The intrinsic macrocycle exhibits innate structural flexibility resulting in dynamic conformational changes. However, the mechanisms controlling flexibility remain unknown. The discovery of a terpene synthase, MicA, that is responsible for the biosynthesis of a nearly non-flexible eunicellane skeleton, enable us to propose a feasible theory about the flexibility in eunicellane structures. Parallel studies of all eunicellane synthases in nature discovered to date, including 2Z-geranylgeranyl diphosphate incubations and density functional theory-based Boltzmann population computations, reveale that a trans-fused bicycle with a 2Z-configuration alkene restricts conformational flexibility resulting in a nearly non-flexible eunicellane skeleton. The catalytic route and the enzymatic mechanism of MicA are also elucidated by labeling experiments, density functional theory calculations, structural analysis of the artificial intelligence-based MicA model, and mutational studies.

SwissADME studies and Density Functional Theory (DFT) approaches of methyl substituted curcumin derivatives

Research suggests curcumin's safety and efficacy, prompting interest in its use for treating and preventing various human diseases. The current study aimed to predict drag ability of methyl substituted curcumin derivatives (BL1 to BL4) using SwissADME and Density Functional Theory (DFT) approaches. The curcumin derivatives investigated mostly adhere to Lipinski's rule of five, with molecular properties including MW, F. Csp3, nHBA, nHBD, and TPSA falling within acceptable limits. The compounds demonstrating high lipophilicity while poor water solubility. The pharmacokinetic evaluation revealed favorable gastrointestinal absorption and blood-brain barrier permeation while none were identified as substrates for P-glycoprotein, however, revealed inhibitory actions against various cytochrome P450 enzymes. Additionally, all derivatives exhibited a consistent bioavailability score of 0.55. Similarly, the DFT computations of the compounds of the curcumin derivatives were conducted at B3LYP/6-311 G** level to predict and then assess the key electronic characteristics underlying the bioactivity. Accordingly, the BL4 molecule (ΔEgap= 4.105 eV) would prefer to interact with the external molecular system more than the other molecules due to having the biggest energy gap. The ΔNmax (2.328 eV) and Δεback-donat. (-0.446 eV) scores implied that BL1 would have more charge transfer capability and the lowest stability via back donation among the compounds. In short, the derivative (BL1 to BL4) exhibited strong extrinsic therapeutic properties and therefore stand eligible for further in vitro and in vivo studies.

Adaptive Response to Solvation in Flexible Molecules: Oligo Hydrates of 4-Hydroxy-2-butanone

Structural changes induced by water play a pivotal role in chemistry and biology but remain challenging to predict, measure, and control at molecular level. Here we explore size-governed gas-phase water aggregation in the flexible molecule 4-hydroxy-2-butanone, modeling the conformational adaptability of flexible substrates to host water scaffolds and the preference for sequential droplet growth. The experiment was conducted using broadband rotational spectroscopy, rationalized with quantum chemical calculations. Two different isomers were observed experimentally from the di- to the pentahydrates (4-hydroxy-2-butanone-(H2O)n=2-5), including the 18O isotopologues for the di- and trihydrates. Interestingly, to accommodate water molecules effectively, the heavy atom skeleton of 4-hydroxy-2-butanone reshapes in every observed isomer and does not correspond to the stable conformer of the free monomer. All solvates initiate from the alcohol group (proton donor) but retain the carbonyl group as secondary binding point. The water scaffolds closely resemble those found in the pure water clusters, balancing between the capability of 4-hydroxy-2-butanone for steering the orientation and position of the water molecules and the ability of water to modulate the monomer's conformation. The present work thus provides an accurate molecular description on how torsionally flexible molecules dynamically adapt to water along progressing solvation.

The Difference between Plasmon Excitations in Chemically Heterogeneous Gold and Silver Atomic Clusters

Weak doping can broaden, shift, and quench plasmon peaks in nanoparticles, but the mechanistic intricacies of the diverse responses to doping remain unclear. In this study, we used the time-dependent density functional theory (TD-DFT) to compute the excitation properties of transition-metal Pd- or Pt-doped gold and silver atomic arrays and investigate the evolution characteristics and response mechanisms of their plasmon peaks. The results demonstrated that the Pd or Pt doping of the off-centered 10 × 2 atomic arrays broadened or shifted the plasmon peaks to varying degrees. In particular, for Pd-doped 10 × 2 Au atomic arrays, the broadened plasmon peak significantly blueshifted, whereas a slight red shift was observed for Pt-doped arrays. For the 10 × 2 Ag atomic arrays, Pd doping caused almost no shift in the plasmon peak, whereas Pt doping caused a substantial red shift in the broadened plasmon peak. The analysis revealed that the diversity in these doping responses was related to the energy positions of the d electrons in the gold and silver atomic clusters and the positions of the doping atomic orbitals in the energy bands. The introduction of doping atoms altered the symmetry and gap size of the occupied and unoccupied orbitals, so multiple modes of single-particle transitions were involved in the excitation. An electron transfer analysis indicated a close correlation between excitation energy and the electron transfer of doping atoms. Finally, the differences in the symmetrically centered 11 × 2 doped atomic array were discussed using electron transfer analysis to validate the reliability of this analytical method. These findings elucidate the microscopic mechanisms of the evolution of plasmon peaks in doped atomic clusters and provide new insights into the rational control and application of plasmons in low-dimensional nanostructures.

THF-Assisted CO2 Reduction Catalyzed by Electride Mg2EP: Insight from DFT Calculations

Hydroboration and hydrogenation reductions of CO2 catalyzed by a porphyrinoid-based dimagnesium(I) electride (Mg2EP) were investigated by density functional theory calculations. Herein, the presence of potentially excess electrons located at the Mg-Mg bond endows Mg2EP with the ability to activate small molecules such as CO2, HBpin, and H2, thus opening up the possibility for further CO2 conversion. The Mg2EP-catalyzed hydroboration of CO2 to HCOOBpin is predicted to have relatively higher activity in comparison to the hydrogenation reduction to formic acid (HCOOH). Interestingly, the common solvent molecule tetrahydrofuran as an auxiliary can coordinate with the Mg center to effectively weaken the bonding interaction between the dimagnesium center and the intermediate species from the CO2 conversion, thereby promoting the catalytic cycle for the CO2 hydroboration. The present results suggest that the electride Mg2EP is promising for the molecular catalyst in the CO2 transformation.

Investigating the stability of aromatic carboxylic acids in hydrated magnesium sulfate under UV irradiation to assist detection of organics on Mars

The Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) instrument onboard the Mars 2020 Perseverance rover detected so far some of the most intense fluorescence signals in association with sulfates analyzing abraded patches of rocks at Jezero crater, Mars. To assess the plausibility of an organic origin of these signals, it is key to understand if organics can survive exposure to ambient Martian UV after exposure by the Perseverance abrasion tool and prior to analysis by SHERLOC. In this work, we investigated the stability of organo-sulfate assemblages under Martian-like UV irradiation and we observed that the spectroscopic features of phthalic and mellitic acid embedded into hydrated magnesium sulfate do not change for UV exposures corresponding to at least 48 Martian sols and, thus, should still be detectable in fluorescence when the SHERLOC analysis takes place, thanks to the photoprotective properties of magnesium sulfate. In addition, different photoproduct bands diagnostic of the parent carboxylic acid molecules could be observed. The photoprotective behavior of hydrated magnesium sulfate corroborates the hypothesis that sulfates might have played a key role in the preservation of organics on Mars, and that the fluorescence signals detected by SHERLOC in association with sulfates could potentially arise from organic compounds.

In-depth exploration of catalytic sites on amorphous solid water: I. The astrosynthesis of aminomethanol

Chemical processes taking place on ice-grain mantles are pivotal to the complex chemistry of interstellar environments. In this study, we conducted a comprehensive analysis of the catalytic effects of an amorphous solid water (ASW) surface on the reaction between ammonia (NH3) and formaldehyde (H2CO) to form aminomethanol (NH2CH2OH) using density functional theory. We identified potential catalytic sites based on the binding energy distribution of NH3 and H2CO reactants, on a set-of-clusters surface model composed of 22 water molecules and found a total of 14 reaction paths. Our results indicate that the catalytic sites can be categorized into four groups, depending on the interactions of the carbonyl oxygen and the amino group with the ice surface in the reactant complex. A detailed analysis of the reaction mechanism using Intrinsic Reaction Coordinate and reaction force analysis, revealed three distinct chemical events for this reaction: formation of the C-N bond, breaking of the N-H bond, and formation of the O-H hydroxyl bond. Depending on the type of catalytic site, these events can occur within a single, concerted, albeit asynchronous, step, or can be isolated in a step-wise mechanism, with the lowest overall transition state energy observed at 1.3 kcal mol-1. A key requirement for the low-energy mechanism is the presence of a pair of dangling OH bonds on the surface, found at 5% of the potential catalytic sites on an ASW porous surface.

Visualization of Multiple-Resonance-Induced Frontier Molecular Orbitals in a Single Multiple-Resonance Thermally Activated Delayed Fluorescence Molecule

The spatial distribution and electronic properties of the frontier molecular orbitals (FMOs) in a thermally activated delayed fluorescence (TADF) molecule contribute significantly to the TADF properties, and thus, a detailed understanding and sophisticated control of the FMOs are fundamental to the design of TADF molecules. However, for multiple-resonance (MR)-TADF molecules that achieve spatial separation of FMOs by the MR effect, the distinctive distribution of these molecular orbitals poses significant challenges for conventional computational analysis and ensemble averaging methods to elucidate the FMOs' separation and the precise mechanism of luminescence. Therefore, the visualization and analysis of electronic states with the specific energy level of a single MR-TADF molecule will provide a deeper understanding of the TADF mechanism. Here, scanning tunneling microscopy/spectroscopy (STM/STS) was used to investigate the electronic states of the DABNA-1 molecule at the atomic scale. FMOs' visualization and local density of states analysis of the DABNA-1 molecule clearly show that MR-TADF molecules also have well-separated FMOs according to the internal heteroatom arrangement, providing insights that complement existing theoretical prediction methods.

Theoretical study of short-range exchange interaction based on semiconductor dielectric function model toward time-dependent dielectric density functional theory

This study explores various models of semiconductor dielectric functions, with a specific emphasis on the large wavenumber spectrum and the derivation of the screened exchange interaction. Particularly, we discuss the short-range effect of the screened exchange potential. Our investigation reveals that the short-range effect originating from the high wavenumber spectrum is contingent upon the dielectric constant of the targeted system. To incorporate dielectric-dependent behaviors concerning the short-range aspect into the dielectric density functional theory (DFT) framework, we utilize the local Slater term and the Yukawa-type term, adjusting the ratio between these terms based on the dielectric constant. Additionally, we demonstrate the efficacy of the time-dependent dielectric DFT method in accurately characterizing the electronic structure of excited states in dyes and functional molecules. Several theoretical approaches have incorporated parameters dependent on the system to elucidate short-range exchange interactions. Our theoretical analysis and discussions will be useful for those studies.

The Role of Lysozyme in the Formation of Bioinspired Silicon Dioxide

Several organisms are able to polycondensate tetraoxosilicic(IV) acid to form silicon(IV) dioxide using polycationic molecules. According to an earlier mechanistic proposal, these molecules undergo a phase separation and recent experimental evidence appears to confirm this model. At the same time, polycationic proteins like lysozyme can also promote polycondensation of silicon(IV) dioxide, and they do so under conditions that are not compatible with liquid-liquid phase separation. In this manuscript we investigate this conundrum by molecular simulations.

Combining Local Range Separation and Local Hybrids: A Step in the Quest for Obtaining Good Energies and Eigenvalues from One Functional

Some of the most successful exchange-correlation approximations in density functional theory are "hybrids", i.e., they rely on combining semilocal density functionals with exact nonlocal Fock exchange. In recent years, two classes of hybrid functionals have emerged as particularly promising: range-separated hybrids on the one hand, and local hybrids on the other hand. These functionals offer the hope to overcome a long-standing "observable dilemma", i.e., the fact that density functionals typically yield either a good description of binding energies, as obtained, e.g., in global and local hybrids, or physically interpretable eigenvalues, as obtained, e.g., in optimally tuned range-separated hybrids. Obtaining both of these characteristics from one and the same functional with the same set of parameters has been a long-standing challenge. We here discuss combining the concepts of local range separation and local hybrids as part of a constraint-guided quest for functionals that overcome the observable dilemma.

Conformational adaptation and large amplitude motions of 1-phenyl-2,2,2-trifluoroethanol with two water molecules: a rotational spectroscopic and ab initio investigation

The 1 : 2 adduct of 1-phenyl-2,2,2-trifluoroethanol (PhTFE), a chiral fluoroalcohol, with two water molecules (PhTFE⋯2H2O) was investigated via chirped pulse Fourier-transform microwave (CP-FTMW) spectroscopy and theoretical calculations. A systematic search of the PhTFE⋯2H2O conformational landscape identified 38 stable minima at the B3LYP-D3BJ/def2-TZVPPD level of theory, 27 of which are within an energy window of 10 kJ mol-1 after applying zero-point energy corrections. Rotational spectra of a single PhTFE⋯2H2O conformer along with eight deuterated and three oxygen-18 isotopologues were assigned. Interestingly, the observed PhTFE⋯2H2O conformer contains PhTFE II, the second most stable monomer conformer, and the most stable PhTFE I dihydrate is ca. 4 kJ mol-1 higher in energy. In contrast, PhTFE I⋯H2O was identified experimentally and theoretically as the most stable 1 : 1 conformer. Furthermore, the observed dihydrate structure experiences large amplitude motions connecting three theoretical minima which differ only in which water oxygen lone pairs are involved in the hydrogen-bonds, i.e., the free OH pointing directions. Additionally, the ortho and para-H2O tunnelling splittings were detected and attributed to the interchange water hydrogen atoms which interact with the aromatic part of PhTFE but not for the water interacting with PhTFE hydroxy group. Extensive theoretical modelling was carried out to gain insight into the associated large amplitude motions including tunnelling, supported by the experimental isotopic and tunnelling splitting data.

Electron-induced ligand loss from iron tetracarbonyl methyl acrylate

We probe the separation of ligands from iron tetracarbonyl methyl acrylate (Fe(CO)4(C4H6O2) or Fe(CO)4MA) induced by the interaction with free electrons. The motivation comes from the possible use of this molecule as a nanofabrication precursor and from the corresponding need to understand its elementary reactions fundamental to the electron-induced deposition. We utilize two complementary electron collision setups and support the interpretation of data by quantum chemical calculations. This way, both the dissociative ionization and dissociative electron attachment fragmentation channels are characterized. Considerable differences in the degree of precursor fragmentation in these two channels are observed. Interesting differences also appear when this precursor is compared to structurally similar iron pentacarbonyl. The present findings shed light on the recent electron-induced chemistry of Fe(CO)4MA on a surface under ultrahigh vacuum.

Predicting spin states of iron porphyrins with DFT methods including crystal packing effects and thermodynamic corrections

Accurate computational treatment of spin states for transition metal complexes, exemplified by iron porphyrins, lies at the heart of quantum bioinorganic chemistry, but at the same time represents a great challenge for approximate density functional theory (DFT) methods, which are predominantly used. Here, the accuracy of DFT methods for spin-state splittings in iron porphyrin is assessed by probing the ability to correctly predict the ground states for six FeIII or FeII complexes experimentally characterized in solid state. For each case, molecular and periodic DFT calculations are employed to quantify the effect of porphyrin side substituents and the crystal packing effect (CPE) on the spin-state splitting. It is proposed to partition the total CPE into additive components, the direct and structural one, the importance of which is shown to significantly vary from case to case. By knowing the substituent effect, the CPE, and the Gibbs free energy thermodynamic correction from calculations, one can employ the experimental ground-state information in order to derive a quantitative constraint on the electronic energy difference for a simplified (porphin) model of the experimentally characterized metalloporphyrin. The constraints derived in such a way-in the form of single or double inequalities-are used to assess the accuracy of dispersion-corrected DFT methods for 6 spin-state splittings of [FeIII(P)(2-MeIm)2]+, [FeIII(P)(2-MeIm)]+, [FeII(P)(THF)2] and [FeII(P)] models (where P is porphin, 2-MeIm is 2-methylimidazole, THF is tetrahydrofuran). These data constitute the new benchmark set of spin states for crystalline iron porphyrins (SSCIP6). The highest accuracy is obtained in the case of double-hybrid functionals (B2PLYP-D3, DSD-PBEB95-D3), whereas hybrid functionals, especially those with reduced admixture of the exact exchange (B3LYP*-D3, TPSSh-D3), are found to considerably overstabilize the intermediate spin state, leading to incorrect ground-state prediction in FeIII porphyrins. The present approach, which can be generalized to other transition metal complexes, is not only useful in method benchmarking, but also sheds light on the interpretations of experimental data for metalloporphyrins, which are important models to understand the electronic properties of heme proteins.

Impact of doping with organic dopants and mixed doping with alkali metals and organic dopants on the absorption, electronic, optoelectronic, thermodynamic and nonlinear optical properties of dibenzo[b,def]chrysene in gaseous media: DFT and TD-DFT studies

CONTEXT

In this study, we evaluate the geometrical, absorption, optoelectronic, electronic, nonlinear optical (NLO) and thermodynamic properties of dibenzo[b,def]chrysene molecule derivatives by means of DFT and TD-DFT simulations. In view of the aim of producing new high-performance materials for non-linear optics (NLO) by doping test, two types of doping were used. We obtained six derivatives by doping with organic dopants (Nitro, amide and ticyanoethenyl) and mixed alkali metal (potassium) and organic dopants. Doping with organic dopants produced molecules A, B and C, respectively when substituting one hydrogen with nitro (NO2), amide (CONH2) and tricyanoethenyl (C5N3) groups, while mixed doping involved considering A, B and C and then substituting two hydrogens with two potassiums to obtain compounds D, E and F respectively. The negative values of the various interaction energies calculated for all the doped molecules show that they are all stable, but also that molecules C and F are the most stable in the case of both dopings. The gap energies calculated at the B3LYP level of theory are all below 3 eV, which means that all the molecules obtained are semiconductors. Better still, compounds C and F, with gap energies of 1.852 eV and 1.204 eV, respectively, corresponding to decreases of 35.67% and 58.18% in gap energy compared with the pristine molecule, are more reactive than the other doped molecules. Mixed doping is therefore a highly effective way of narrowing the energy gap and boosting the semiconducting character and reactivity of organic materials. Optoelectronic properties have also been improved, with refractive index values higher than those of the reference material, glass. This shows that our compounds could be used under very high electric field conditions of the order of 4.164 × 109 V.m-1 for C and 7.410 × 109 V.m-1 for F the highest values at the B3LYP level of theory. The maximum first-order hyperpolarizability values for both types of doping are obtained at the CAM-B3LYP level of theory by C:β mol = 92.088 × 10-30esu and by F:β mol = 129.449 × 10-30esu, and second-order values are also given by these same compounds. These values are higher than the reference value, which is urea, making our compounds potential candidates for high-performance NLO applications. In dynamic mode and at a frequency of 1064 nm, at the CAM-B3LYP level of theory, the highest dynamic hyperpolarizability coefficients were obtained by C and F. Hyper-Rayleigh scattering β HRS , coefficients of the electro-optical Pockel effect (EOPE), EFISHG, third-order NLO-response degree four-wave mixingγ DFWM , quadratic nonlinear refractive index n2 were also calculated. The maximum values of n2 are obtained by C (6.13 × 10-20 m2/W) and F (6.60 × 10-20 m2/W), these values are 2.24 times higher than that of fused silica which is the reference for degenerate four-wave mixing so our molecules could also have applications in optoelectronics as wavelength converters, optical pulse modulators and optical switches.

METHODS

Using the DFT method, we were able to determine the optimized and stable electronic structures of doped dibenzo[b,def]chrysene derivatives in the gas phase. We limited ourselves to using the proven B3LYP and CAMB3LYP levels of theory for calculating electronic properties, and non-linear optics with the 6-311G +  + (d,p) basis set, which is a large basis set frequently used for these types of compound. Gaussian 09 software was used to run our calculations, and Gauss View 6.0.16 was used to visualize the output files. TD-DFT was also used to determine absorption properties at the B3LYP level of theory, using the same basis set.

Toward the Next Generation of Density Functionals: Escaping the Zero-Sum Game by Using the Exact-Exchange Energy Density

ConspectusKohn-Sham density functional theory (KS DFT) is arguably the most widely applied electronic-structure method with tens of thousands of publications each year in a wide variety of fields. Its importance and usefulness can thus hardly be overstated. The central quantity that determines the accuracy of KS DFT calculations is the exchange-correlation functional. Its exact form is unknown, or better "unknowable", and therefore the derivation of ever more accurate yet efficiently applicable approximate functionals is the "holy grail" in the field. In this context, the simultaneous minimization of so-called delocalization errors and static correlation errors is the greatest challenge that needs to be overcome as we move toward more accurate yet computationally efficient methods. In many cases, an improvement on one of these two aspects (also often termed fractional-charge and fractional-spin errors, respectively) generates a deterioration in the other one. Here we report on recent notable progress in escaping this so-called "zero-sum-game" by constructing new functionals based on the exact-exchange energy density. In particular, local hybrid and range-separated local hybrid functionals are discussed that incorporate additional terms that deal with static correlation as well as with delocalization errors. Taking hints from other coordinate-space models of nondynamical and strong electron correlations (the B13 and KP16/B13 models), position-dependent functions that cover these aspects in real space have been devised and incorporated into the local-mixing functions determining the position-dependence of exact-exchange admixture of local hybrids as well as into the treatment of range separation in range-separated local hybrids. While initial functionals followed closely the B13 and KP16/B13 frameworks, meanwhile simpler real-space functions based on ratios of semilocal and exact-exchange energy densities have been found, providing a basis for relatively simple and numerically convenient functionals. Notably, the correction terms can either increase or decrease exact-exchange admixture locally in real space (and in interelectronic-distance space), leading even to regions with negative admixture in cases of particularly strong static correlations. Efficient implementations into a fast computer code (Turbomole) using seminumerical integration techniques make such local hybrid and range-separated local hybrid functionals promising new tools for complicated composite systems in many research areas, where simultaneously small delocalization errors and static correlation errors are crucial. First real-world application examples of the new functionals are provided, including stretched bonds, symmetry-breaking and hyperfine coupling in open-shell transition-metal complexes, as well as a reduction of static correlation errors in the computation of nuclear shieldings and magnetizabilities. The newest versions of range-separated local hybrids (e.g., ωLH23tdE) retain the excellent frontier-orbital energies and correct asymptotic exchange-correlation potential of the underlying ωLH22t functional while improving substantially on strong-correlation cases. The form of these functionals can be further linked to the performance of the recent impactful deep-neural-network "black-box" functional DM21, which itself may be viewed as a range-separated local hybrid.

Methyl α-D-galactopyranosyl-(1→3)-β-D-galactopyranoside and methyl β-D-galactopyranosyl-(1→3)-β-D-galactopyranoside: Glycosidic linkage conformation determined from MA'AT analysis

MA'AT analysis has been applied to two biologically-important O-glycosidic linkages in two disaccharides, α-D-Galp-(1→3)-β-D-GalpOMe (3) and β-D-Galp-(1→3)-β-D-GalpOMe (4). Using density functional theory (DFT) to obtain parameterized equations relating a group of trans-O-glycosidic NMR spin-couplings to either phi (ϕ') or psi (ψ'), and experimental 3JCOCH, 2JCOC, and 3JCOCC spin-couplings measured in aqueous solution in 13C-labeled isotopomers, probability distributions of ϕ' and ψ' in each linkage were determined and compared to those determined by aqueous 1-μs molecular dynamics (MD) simulation. Good agreement was found between the MA'AT and single-state MD conformational models of these linkages for the most part, with modest (approximately <15°) differences in the mean values of ϕ' and ψ', although the envelope of allowed angles (encoded in circular standard deviations or CSDs) is consistently larger for ϕ' determined from MA'AT analysis than from MD for both linkages. The MA'AT model of the α-Galp-(1→3)-β-Galp linkage agrees well with those determined previously using conventional NMR methods (3JCOCH values and/or 1H-1H NOEs), but some discrepancy was observed for the β-Galp-(1→3)-β-Galp linkage, which may arise from errors in the conventions used to describe the linkage torsion angles. Statistical analyses of X-ray crystal structures show ranges of ϕ' and ψ' for both linkages that include the mean angles determined from MA'AT analyses, although both angles adopt a wide range of values in the crystalline state, with ϕ' in β-Galp-(1→3)-β-Galp linkages showing greater-than-expected conformational variability.

In Silico Chemical Experiments in the Age of AI: From Quantum Chemistry to Machine Learning and Back

Computational chemistry is an indispensable tool for understanding molecules and predicting chemical properties. However, traditional computational methods face significant challenges due to the difficulty of solving the Schrödinger equations and the increasing computational cost with the size of the molecular system. In response, there has been a surge of interest in leveraging artificial intelligence (AI) and machine learning (ML) techniques to in silico experiments. Integrating AI and ML into computational chemistry increases the scalability and speed of the exploration of chemical space. However, challenges remain, particularly regarding the reproducibility and transferability of ML models. This review highlights the evolution of ML in learning from, complementing, or replacing traditional computational chemistry for energy and property predictions. Starting from models trained entirely on numerical data, a journey set forth toward the ideal model incorporating or learning the physical laws of quantum mechanics. This paper also reviews existing computational methods and ML models and their intertwining, outlines a roadmap for future research, and identifies areas for improvement and innovation. Ultimately, the goal is to develop AI architectures capable of predicting accurate and transferable solutions to the Schrödinger equation, thereby revolutionizing in silico experiments within chemistry and materials science.