Harvesting Chemical Understanding with Machine Learning and Quantum Computers

It is tenable to argue that nobody can predict the future with certainty, yet one can learn from the past and make informed projections for the years ahead. In this Perspective, we overview the status of how theory and computation can be exploited to obtain chemical understanding from wave function theory and density functional theory, and then outlook the likely impact of machine learning (ML) and quantum computers (QC) to appreciate traditional chemical concepts in decades to come. It is maintained that the development and maturation of ML and QC methods in theoretical and computational chemistry represent two paradigm shifts about how the Schrödinger equation can be solved. New chemical understanding can be harnessed in these two new paradigms by making respective use of ML features and QC qubits. Before that happens, however, we still have hurdles to face and obstacles to overcome in both ML and QC arenas. Possible pathways to tackle these challenges are proposed. We anticipate that hierarchical modeling, in contrast to multiscale modeling, will emerge and thrive, becoming the workhorse of in silico simulations in the next few decades.

Modulating Aluminum Solvation with Ionic Liquids for Improved Aqueous-Based Aluminum-Ion Batteries

Aqueous-based Al-ion batteries are attractive alternatives to Li-ion batteries due to their safety, high volumetric energy density, abundance, and recyclability. Although aluminum-ion batteries are attractive, there are major challenges to overcome, which include understanding the nature of the passive layer of aluminum oxide on the aluminum anode, the narrow electrochemical window of aqueous electrolytes, and lack of suitable cathodes. Here, we report using experiments in conjunction with DFT simulations to clarify the role of ionic liquids (ILs) in altering the Al solvation dynamics, which in turn affects the aluminum electrochemistry and aqueous-based battery performance significantly. DFT calculations showed that the addition of 1-ethyl-3-methylimidazolium trifluoromethylsulfonate (EMIMTfO) changes the aluminum solvation structure in the aqueous (Al(TfO)3) electrolyte to lower coordinated solvation shells, thereby influencing and improving Al deposition/stripping on the Zn/Al alloy anode. Furthermore, the addition of an IL reduces the strain in manganese oxide during intercalation/deintercalation, thereby improving the Zn/Al-MnOx battery performance. By optimizing the electrolyte composition, a battery potential of >1.7 V was achieved for the Zn/Al-MnOx system.

Nature of Catalytic Behavior of Cobalt Oxides for CO2 Hydrogenation

Cobalt oxide (CoO _x ) catalysts are widely applied in CO₂ hydrogenation but suffer from structural evolution during the reaction. This paper describes the complicated structure-performance relationship under reaction conditions. An iterative approach was employed to simulate the reduction process with the help of neural network potential-accelerated molecular dynamics. Based on the reduced models of catalysts, a combined theoretical and experimental study has discovered that CoO(111) provides active sites to break C-O bonds for CH₄ production. The analysis of the reaction mechanism indicated that the C-O bond scission of *CH₂O species plays a key role in producing CH₄. The nature of dissociating C-O bonds is attributed to the stabilization of *O atoms after C-O bond cleavage and the weakening of C-O bond strength by surface-transferred electrons. This work may offer a paradigm to explore the origin of performance over metal oxides in heterogeneous catalysis.

Active Learning Exploration of Transition-Metal Complexes to Discover Method-Insensitive and Synthetically Accessible Chromophores

Transition-metal chromophores with earth-abundant transition metals are an important design target for their applications in lighting and nontoxic bioimaging, but their design is challenged by the scarcity of complexes that simultaneously have well-defined ground states and optimal target absorption energies in the visible region. Machine learning (ML) accelerated discovery could overcome such challenges by enabling the screening of a larger space but is limited by the fidelity of the data used in ML model training, which is typically from a single approximate density functional. To address this limitation, we search for consensus in predictions among 23 density functional approximations across multiple rungs of "Jacob's ladder". To accelerate the discovery of complexes with absorption energies in the visible region while minimizing the effect of low-lying excited states, we use two-dimensional (2D)efficient global optimization to sample candidate low-spin chromophores from multimillion complex spaces. Despite the scarcity (i.e., ∼0.01%) of potential chromophores in this large chemical space, we identify candidates with high likelihood (i.e., >10%) of computational validation as the ML models improve during active learning, representing a 1000-fold acceleration in discovery. Absorption spectra of promising chromophores from time-dependent density functional theory verify that 2/3 of candidates have the desired excited-state properties. The observation that constituent ligands from our leads have demonstrated interesting optical properties in the literature exemplifies the effectiveness of our construction of a realistic design space and active learning approach.

High-Frequency Sheet Conductance of Nanolayered WS2 Crystals for Two-Dimensional Nanodevices

Time-resolved terahertz (THz) spectroscopy is a powerful technique for the determination of charge transport properties in photoexcited semiconductors. However, the relatively long wavelengths of THz radiation and the diffraction limit imposed by optical imaging systems reduce the applicability of THz spectroscopy to large samples with dimensions in the millimeter to centimeter range. Exploiting THz near-field spectroscopy, we present the first time-resolved THz measurements on a single exfoliated 2D nanolayered crystal of a transition metal dichalcogenide (WS₂). The high spatial resolution of THz near-field spectroscopy enables mapping of the sheet conductance for an increasing number of atomic layers. The single-crystalline structure of the nanolayered crystal allows for the direct observation of low-energy phonon modes, which are present in all thicknesses, coupling with free carriers. Density functional theory calculations show that the phonon mode corresponds to the breathing mode between atomic layers in the weakly bonded van der Waals layers, which can be strongly influenced by substrate-induced strain. The non-invasive and high-resolution mapping technique of carrier dynamics in nanolayered crystals by time-resolved THz time domain spectroscopy enables possibilities for the investigation of the relation between phonons and charge transport in nanoscale semiconductors for applications in two-dimensional nanodevices.

Elucidation of Cu-Zn Surface Alloying on Cu(997) by Machine-Learning Molecular Dynamics

The Cu-Zn surface alloy has been extensively involved in the investigation of the true active site of Cu/ZnO/Al₂O₃, the industrial catalyst for methanol synthesis which remains under controversy. The challenge lies in capturing the interplay between the surface and reaction under operating conditions, which can be overcome given that the explicit dynamics of the system is known. To provide a better understanding of the dynamic of Cu-Zn surface at the atomic level, the structure and the formation process of the Cu-Zn surface alloy on Cu(997) were investigated by machine-learning molecular dynamics (MD). Gaussian process regression aided with on-the-fly learning was employed to build the force field used in the MD. The simulation reveals atomistic details of the alloying process, that is, the incorporation of deposited Zn adatoms to the Cu substrate. The surface alloying is found to start at upper and lower terraces near the step edge, which emphasize the role of steps and kinks in the alloying. The incorporation of Zn at the middle terrace was found at the later stage of the simulation. The rationalization of alloying behavior was performed based on statistics and barriers of various elementary events that occur during the simulation. It was observed that the alloying scheme at the upper terrace is dominated by the confinement of Zn step adatoms by other adatoms, highlighting the importance of step fluctuations in the alloying process. On the other hand, the alloying scheme at the lower terrace is dominated by direct exchange between the Zn step adatom and the Cu atom underneath. The alloying at the middle terrace is dominated by the wave deposition mechanism and deep confinement of Zn adatoms. The short propagation of alloyed Zn in the middle terrace was observed to proceed by means of indirect exchange instead of local exchange as proposed in the previous scanning tunneling microscopy (STM) observation. The comparison of migration rate and activation energies to the result of STM observation is also made. We have found that at a certain distance from the surface, the STM tip significantly affects the elementary events such as vacancy formation and direct exchange.

Defect Engineering in Metal-Organic Framework Nanocrystals: Implications for Mechanical Properties and Performance

The growth process of metal-organic framework (MOF) nanocrystals defines their properties and functions. However, defects may be prevalent during the crystallization of even seemingly perfect MOFs, such as zeolitic imidazolate framework-8 (ZIF-8), and yet direct probing of such structural defects has been challenging because of the lack of nanoscale techniques to locally examine individual nanocrystals. Here, we directly study local defects, such as missing linkers or metal vacancies, in ZIF-8 nano- and microcrystals with near-field IR nanospectroscopy combined with density functional theory calculations. We track the chemical changes during crystallization and show that structural defects like zinc cations that are bound to molecules of the reactant gradually disappear with ripening of the crystals, while dangling and missing linker defects prevail. The resulting defect-terminating groups or open-metal sites produce mechanical anisotropy and reduce the Young's modulus, as measured via tip force microscopy with nanoscale resolution and supported by theoretical modeling. However, these structural defects also open the door for defect engineering to tune the performance of ZIF-8 by offering additional adsorption sites for targeted catalytic reactions, chemical sensing, or gas capture.

Di- and Tetrameric Molybdenum Sulfide Clusters Activate and Stabilize Dihydrogen as Hydrides

NaY zeolite-encapsulated dimeric (Mo₂S₄) and tetrameric (Mo₄S₄) molybdenum sulfide clusters stabilize hydrogen as hydride binding to Mo atoms. Density functional theory (DFT) calculations and adsorption measurements suggest that stabilization of hydrogen as sulfhydryl (SH) groups, as typical for layered MoS₂, is thermodynamically disfavored. Competitive adsorption of H₂ and ethene on Mo was probed by quantifying adsorbed CO on partly hydrogen and/or ethene covered samples with IR spectroscopy. During hydrogenation, experiment and theory suggest that Mo is covered predominately with ethene and sparsely with hydride. DFT calculations further predict that, under reaction conditions, each Mo _x S _y cluster can activate only one H₂, suggesting that the entire cluster (irrespective of its nuclearity) acts as one active site for hydrogenation. The nearly identical turnover frequencies (24.7 ± 3.3 mol_ethane·h^-1·mol_cluster ^-1), apparent activation energies (31-32 kJ·mol^-1), and reaction orders (∼0.5 in ethene and ∼1.0 in H₂) show that the active sites in both clusters are catalytically indistinguishable.

Understanding the Photocatalytic Activity of La5Ti2AgS5O7 and La5Ti2CuS5O7 for Green Hydrogen Production: Computational Insights

Green production of hydrogen is possible with photocatalytic water splitting, where hydrogen is produced while water is reduced by using energy derived from light. In this study, density functional theory (DFT) is employed to gain insights into the photocatalytic performance of La₅Ti₂AgS₅O₇ and La₅Ti₂CuS₅O₇-two emerging candidate materials for water splitting. The electronic structure of both bulk materials was calculated by using hybrid DFT, which indicated the band gaps and charge carrier effective masses are suitable for photocatalytic water splitting. Notably, the unique one-dimensional octahedral TiO _x S_6-x and tetragonal MS₄ channels formed provide a structural separation for photoexcited charge carriers which should inhibit charge recombination. Band alignments of surfaces that appear on the Wulff constructions of 12 nonpolar symmetric surface slabs were calculated by using hybrid DFT for each of the materials. All surfaces of La₅Ti₂AgS₅O₇ have band edge positions suitable for hydrogen evolution; however, the small overpotentials on the largest facets likely decrease the photocatalytic activity. In La₅Ti₂CuS₅O₇, 72% of the surface area can support oxygen evolution thermodynamically and kinetically. Based on their similar electronic structures, La₅Ti₂AgS₅O₇ and La₅Ti₂CuS₅O₇ could be effectively employed in Z-scheme photocatalytic water splitting.

Superatomic Au25(SC2H5)18 Nanocluster under Pressure

The past decade has witnessed significant advances in the synthesis and structure determination of atomically precise metal nanoclusters. However, little is known about the condensed matter properties of these nanosized metal nanoclusters packed in a crystal lattice under high pressure. Here using density functional theory calculations, we simulate the crystal of a representative superatomic gold cluster, Au₂₅(SR)₁₈ ⁰ (R = C₂H₅), under various pressures. At ambient conditions, Au₂₅(SC₂H₅)₁₈ ⁰ clusters are packed in a crystal via dispersion interactions; being a 7e superatom, each cluster carries a magnetic moment of 1 μ_B or one unpaired electron. Upon increasing compression (from 10 to 110 GPa), we observe the formation of intercluster Au-Au, Au-S, and S-S covalent bonds between staple motifs, thereby linking the clusters into a network. The pressure-induced structural change is accompanied by the vanishment of the magnetic moment and the semiconductor-to-metal transition. Our work shows that subjecting crystals of atomically precise metal nanoclusters to high pressures could lead to new crystalline states and physical properties.

Quantifying Effects of Active Site Proximity on Rates of Methanol Dehydration to Dimethyl Ether over Chabazite Zeolites through Microkinetic Modeling

Control of the spatial proximity of Brønsted acid sites within the zeolite framework can result in materials with properties that are distinct from materials synthesized through conventional crystallization methods or available from commercial sources. Recent experimental evidence has shown that turnover rates of different acid-catalyzed reactions increase with the fraction of proximal sites in chabazite (CHA) zeolites. The catalytic conversion of oxygenates is an important research area, and the dehydration of methanol to dimethyl ether (DME) is a well-studied reaction as part of methanol-to-olefin chemistry catalyzed by solid acids. Published experimental data have shown that DME formation rates (per acid site) increase systematically with the fraction of proximal acid sites in the six-membered ring of CHA. Here, we probe the effect of acid site proximity in CHA on methanol dehydration rates using electronic structure calculations and microkinetic modeling to identify the primary causes of this chemistry and their relationship to the local structure of the catalyst at the nanoscale. We report a density functional theory-parametrized microkinetic model of methanol dehydration to DME, catalyzed by acidic CHA zeolite with direct comparison to experimental data. Effects of proximal acid sites on reaction rates were captured quantitatively for a range of operating conditions and catalyst compositions, with a focus on total paired acid site concentration and reactant clustering to form higher nuclearity complexes. Next-nearest neighbor paired acid sites were identified as promoting the formation of methanol trimer clusters rather than the inhibiting tetramer or pentamer clusters, resulting in large increases in the rate for DME production due to the lower energy barriers present in the concerted methanol trimer reaction pathway. The model framework developed in this study can be extended to other zeolite materials and reaction chemistries toward the goal of rational design and development of next-generation catalytic materials and chemical processes.

Computational Screening of First-Row Transition-Metal Based Alloy Catalysts-Ligand Induced N2 Reduction Reaction Selectivity

Large-scale ammonia production through sustainable strategies from naturally abundant N₂ under ambient conditions represents a major challenge from a future perspective. Ammonia is one of the promising carbon-free alternative energy carriers. The high energy required for N≡N bond dissociation during the Haber-Bosch process demands extreme reaction conditions. This problem could be circumvented by tuning Fe catalyst composition with the help of an induced ligand effect on the surface. In this work, we utilized density functional theory calculations on the Fe(110) surface alloyed with first-row transition-metal (TM) series (Fe-TM) to understand the catalytic activity that facilitates the electrochemical nitrogen reduction reaction (NRR). We also calculated the selectivity against the competitive hydrogen evolution reaction (HER) under electrochemical conditions. The calculated results are compared with those from earlier reports on the periodic Fe(110) and Fe(111) surfaces, and also on the (110) surface of the Fe₈₅ nanocluster. Surface alloying with late TMs (Co, Ni, Cu) shows an improved NRR activity, whereas the low exchange current density observed for Fe-Co indicates less HER activity among them. Considering various governing factors, Fe-based alloys with Co (Fe-Co) showed enhanced overall performance compared to the periodic surface as well as other pure iron-based structures previously reported. Therefore, the iron-alloy based structured catalysts may also provide more opportunities in the future for enhancing NRR performance via electrochemical reduction pathways.

Elucidating the Mechanism of Ambient-Temperature Aldol Condensation of Acetaldehyde on Ceria

Using in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and density functional theory (DFT) calculations, we conclusively demonstrate that acetaldehyde (AcH) undergoes aldol condensation when flown over ceria octahedral nanoparticles, and the reaction is desorption-limited at ambient temperature. trans-Crotonaldehyde (CrH) is the predominant product whose coverage builds up on the catalyst with time on stream. The proposed mechanism on CeO2(111) proceeds via AcH enolization (i.e., α C-H bond scission), C-C coupling, and further enolization and dehydroxylation of the aldol adduct, 3-hydroxybutanal, to yield trans-CrH. The mechanism with its DFT-calculated parameters is consistent with reactivity at ambient temperature and with the kinetic behavior of the aldol condensation of AcH reported on other oxides. The slightly less stable cis-CrH can be produced by the same mechanism depending on how the enolate and AcH are positioned with respect to each other in C-C coupling. All vibrational modes in DRIFTS are identified with AcH or trans-CrH, except for a feature at 1620 cm-1 that is more intense relative to the other bands on the partially reduced ceria sample than on the oxidized sample. It is identified to be the C=C stretch mode of CH3CHOHCHCHO adsorbed on an oxygen vacancy. It constitutes a deep energy minimum, rendering oxygen vacancies an inactive site for CrH formation under given conditions.

Improved Pd/CeO2 Catalysts for Low-Temperature NO Reduction: Activation of CeO2 Lattice Oxygen by Fe Doping

Developing better three-way catalysts with improved low-temperature performance is essential for cold start emission control. Density functional theory in combination with microkinetics simulations is used to predict reactivity of CO/NO/H2 mixtures on a small Pd cluster on CeO2(111). At low temperatures, N2O formation occurs via a N2O2 dimer over metallic Pd3. Part of the N2O intermediate product re-oxidizes Pd, limiting NO conversion and requiring rich conditions to obtain high N2 selectivity. High N2 selectivity at elevated temperatures is due to N2O decomposition on oxygen vacancies. Doping CeO2 by Fe is predicted to lead to more oxygen vacancies and a higher N2 selectivity, which is validated by the lower onset of N2 formation for a Pd catalyst supported on Fe-doped CeO2 prepared by flame spray pyrolysis. Activating ceria surface oxygen by transition metal doping is a promising strategy to improve the performance of three-way catalysts.

Structural, Electronic, and Vibrational Properties of a Two-Dimensional Graphdiyne-like Carbon Nanonetwork Synthesized on Au(111): Implications for the Engineering of sp-sp2 Carbon Nanostructures

Graphdiyne, atomically thin two-dimensional (2D) carbon nanostructure based on sp-sp2 hybridization is an appealing system potentially showing outstanding mechanical and optoelectronic properties. Surface-catalyzed coupling of halogenated sp-carbon-based molecular precursors represents a promising bottom-up strategy to fabricate extended 2D carbon systems with engineered structure on metallic substrates. Here, we investigate the atomic-scale structure and electronic and vibrational properties of an extended graphdiyne-like sp-sp2 carbon nanonetwork grown on Au(111) by means of the on-surface synthesis. The formation of such a 2D nanonetwork at its different stages as a function of the annealing temperature after the deposition is monitored by scanning tunneling microscopy (STM), Raman spectroscopy, and combined with density functional theory (DFT) calculations. High-resolution STM imaging and the high sensitivity of Raman spectroscopy to the bond nature provide a unique strategy to unravel the atomic-scale properties of sp-sp2 carbon nanostructures. We show that hybridization between the 2D carbon nanonetwork and the underlying substrate states strongly affects its electronic and vibrational properties, modifying substantially the density of states and the Raman spectrum compared to the free standing system. This opens the way to the modulation of the electronic properties with significant prospects in future applications as active nanomaterials for catalysis, photoconversion, and carbon-based nanoelectronics.

Recent Developments in Enantioselective Transition Metal Catalysis Featuring Attractive Noncovalent Interactions between Ligand and Substrate

Enantioselective transition metal catalysis is an area very much at the forefront of contemporary synthetic research. The development of processes that enable the efficient synthesis of enantiopure compounds is of unquestionable importance to chemists working within the many diverse fields of the central science. Traditional approaches to solving this challenge have typically relied on leveraging repulsive steric interactions between chiral ligands and substrates in order to raise the energy of one of the diastereomeric transition states over the other. By contrast, this Review examines an alternative tactic in which a set of attractive noncovalent interactions operating between transition metal ligands and substrates are used to control enantioselectivity. Examples where this creative approach has been successfully applied to render fundamental synthetic processes enantioselective are presented and discussed. In many of the cases examined, the ligand scaffold has been carefully designed to accommodate these attractive interactions, while in others, the importance of the critical interactions was only elucidated in subsequent computational and mechanistic studies. Through an exploration and discussion of recent reports encompassing a wide range of reaction classes, we hope to inspire synthetic chemists to continue to develop asymmetric transformations based on this powerful concept.

Pyrazole Derivatives of Medically Relevant Phenolic Acids: Insight into Antioxidative and Anti-LOX Activity

BACKGROUND

From the point of view of medicinal chemistry, compounds containing phenolic and pyrazolic moiety are significant since they are often constituents of bioactive compounds.

OBJECTIVE

The aims of this study were to synthesize pyrazole derivatives of medically relevant phenolic acids, confirm their structure, and evaluate their antioxidative and anti-LOX activities.

METHODS

Phenolic pyrazole derivatives were obtained, starting from esters of medically relevant phenolic acids. The structures of all obtained compounds were determined by NMR and IR spectroscopy, and UV-Vis spectrophotometry. In addition, the single-crystal X-ray diffraction was used. Pyrazole derivatives were tested for their in vitro antioxidative (DPPH assay), and lipoxygenase (LOX) inhibitory activities. Radical quenching mechanism was estimated using DFT and thermodynamic approach, while molecular docking was used to estimate the binding mode within the enzyme.

RESULTS

Pyrazole derivatives were obtained in high yields. The crystal structure of a new compound 3e was determined. Pyrazole derivative with catechol moiety 3d exhibited excellent radical scavenging activity, while compound 3b exhibited the best anti-LOX activity. Molecular docking study revealed that there is no direct interaction of any ligand with the active site of LOX-Ib, but pyrazoles 3a-e behave as inhibitors blocking the approach of linoleic acid to the active site.

CONCLUSION

In this research, protocatechuic and vanillic acid pyrazole derivatives have been obtained for the first time. In vitro antioxidative assay suggests that pyrazole derivate of protocatechuic acid is a powerful radical scavenger, while anti-LOX assay indicates a pyrazole derivative with 4-hydroxyphenyl moiety.

Structure and Properties of Nanosilicates with Olivine (Mg2SiO4) N and Pyroxene (MgSiO3) N Compositions

Magnesium-rich silicates are ubiquitous both terrestrially and astronomically, where they are often present as small particles. Nanosized Mg-rich silicate particles are likely to be particularly important for understanding the formation, processing, and properties of cosmic dust grains. Although astronomical observations and laboratory studies have revealed much about such silicate dust, our knowledge of this hugely important class of nanosolids largely rests on top-down comparisons with the properties of bulk silicates. Herein, we provide a foundational bottom-up study of the structure and properties of Mg-rich nanosilicates based on carefully procured atomistic models. Specifically, we employ state-of-the-art global optimization methods to search for the most stable structures of silicate nanoclusters with olivine (Mg₂SiO₄) _N and pyroxene (MgSiO₃) _N compositions with N = 1-10. To ensure the reliability of our searches, we develop a new interatomic potential that has been especially tuned for nanosilicates. Subsequently, we refine these searches and calculate a range of physicochemical properties of the most stable nanoclusters using accurate density functional theory based electronic structure calculations. We report a detailed analysis of structural and energy properties, charge distributions, and infrared vibrational spectra, where in all cases we compare our finding for nanosilicates with those of the corresponding bulk silicate crystals. For most properties considered, we find large differences with respect to the bulk limit, underlining the limitations of a top-down approach for describing these species. Overall, our work provides a new platform for an accurate and detailed understanding of nanoscale silicates.

Atomistic origins of high-performance in hybrid halide perovskite solar cells

The performance of organometallic perovskite solar cells has rapidly surpassed that of both conventional dye-sensitized and organic photovoltaics. High-power conversion efficiency can be realized in both mesoporous and thin-film device architectures. We address the origin of this success in the context of the materials chemistry and physics of the bulk perovskite as described by electronic structure calculations. In addition to the basic optoelectronic properties essential for an efficient photovoltaic device (spectrally suitable band gap, high optical absorption, low carrier effective masses), the materials are structurally and compositionally flexible. As we show, hybrid perovskites exhibit spontaneous electric polarization; we also suggest ways in which this can be tuned through judicious choice of the organic cation. The presence of ferroelectric domains will result in internal junctions that may aid separation of photoexcited electron and hole pairs, and reduction of recombination through segregation of charge carriers. The combination of high dielectric constant and low effective mass promotes both Wannier-Mott exciton separation and effective ionization of donor and acceptor defects. The photoferroic effect could be exploited in nanostructured films to generate a higher open circuit voltage and may contribute to the current-voltage hysteresis observed in perovskite solar cells.