Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex chemical systems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Editorial: Particulate Matter Emissions From Conventional and Reformulated Fuel Combustion: Advances in Experiments and Simulations

Particulate matter (PM or soot) emissions generated by combustion of conventional and reformulated fuels adversely impact the environment and human health (Bond et al., 2013; Landrigan et al., 2018; World Health Organization (WHO), 2016), generating public awareness and motivating efforts towards the mitigation of their harmful effects. Soot formation is one of the most complex phenomena in combustion, involving interactions between combustion chemistry, fluid mechanics, mass/heat transport, and particle dynamics, spanning different spatial and temporal scales. A fundamental understanding of soot formation process is thus necessary to achieve a strong reduction of PM emissions and design cleaner and more efficient combustion systems. This has motivated a longstanding and ongoing research activity aimed at improving our understanding of the physical and chemical processes involved in soot formation, well reviewed in these recent works (Wang and Chung, 2019; Michelsen et al., 2020; Martin et al., 2022). Despite the wide interest, the transition from gas-phase molecules to incipient soot particles is still elusive and the successive particle growth and oxidation processes are far from being fully understood, especially in conditions relevant to real-world applications. Recent advances in combustion PM emission diagnostic and computational capabilities helped in improving the predictability of fundamental chemical and aerosol models for practical applications, thus tackling some of the above-mentioned challenges. The aim of this research topic is to display the ongoing research efforts in addressing the existing gaps on particulate formation from various fuel sources and in conditions typical of practical combustion applications (e.g. flames, engines, pool fires), through both experimental and numerical approaches.

42 ENGINEERING↗

A non-orthogonal variational quantum eigensolver

Variational algorithms for strongly correlated chemical and materials systems are one of the most promising applications of near-term quantum computers. We present an extension to the variational quantum eigensolver that approximates the ground state of a system by solving a generalized eigenvalue problem in a subspace spanned by a collection of parametrized quantum states. This allows for the systematic improvement of a logical wavefunction ansatz without a significant increase in circuit complexity. To minimize the circuit complexity of this approach, we propose a strategy for efficiently measuring the Hamiltonian and overlap matrix elements between states parametrized by circuits that commute with the total particle number operator. This strategy doubles the size of the state preparation circuits but not their depth, while adding a small number of additional two-qubit gates relative to standard variational quantum eigensolver. We also propose a classical Monte Carlo scheme to estimate the uncertainty in the ground state energy caused by a finite number of measurements of the matrix elements. We explain how this Monte Carlo procedure can be extended to adaptively schedule the required measurements, reducing the number of circuit executions necessary for a given accuracy. We apply these ideas to two model strongly correlated systems, a square configuration of H4 and the π-system of hexatriene (C6H8).

97 MATHEMATICS AND COMPUTING↗

Utilisation of flavonoids from Lycium barbarum L. leaves cross‐linked with whey protein for the stabilisation and delivery of β ‐carotene emulsions

Summary Although β ‐carotene is an essential nutrient for the maintenance of human health, it is insoluble in water and extremely unstable, which limits its application. Here, we prepared a stable emulsification system of whey protein (WP)‐ Lycium barbarum L. leaves flavonoids (LBLF) complex as a carrier to load β ‐carotene. The rheological properties, antioxidant properties and stability of emulsion were studied, and the effects of emulsion on the biological retention of β ‐carotene and stability were investigated. The results indicate that, compared with WP, the WP‐LBLF emulsion has smaller particles, a smaller creaming index, and significantly enhanced chemical stability under light and heat treatment conditions. The thiobarbituric acid value of the WP‐LBLF emulsion (0.031 ± 0.016 mg/kg) was dramatically lower than that of the WP emulsion (0.045 ± 0.012 mg/kg) ( P < 0.05). In the simulated gastrointestinal digestion experiment, the retention rate of β ‐carotene was higher in WP‐LBLF emulsion. Additionally, the WP‐LBLF emulsion loaded with β ‐carotene has good salt tolerance and heat stability. In our study, WP‐LBLF was used as an emulsifier to load β ‐carotene emulsion to improve its bioavailability and stability, which provides a new preparation and application method for a complex functional emulsifier.

Fan, Yanli↗

Accelerating the identification of novel secondary metabolites in bioenergy plant root exudates using MicroED

Small molecule metabolites drive inter- and intraspecies communication and dependencies in diverse biological systems, yet a large proportion of these important chemical compounds remain uncharacterized in plants and microbes. Approximately 90% of the metabolites in root exudate profiles are unknown compounds, despite the importance of root exudate composition in plant-microbe interactions. We need advanced analytical capabilities that will support rapid discovery and structural elucidation of metabolites from biological samples that may be limited in quantity and high in complexity. To fill this gap, this project aimed to develop an integrated workflow involving metabolite extraction, separation, and crystallization from plant root exudates followed by characterization using nuclear magnetic resonance (NMR) spectroscopy, mass spectrometry, and microcrystal electron diffraction (MicroED). Using crude root exudates from sorghum, this project successfully developed higher throughput exudate fractionation strategies to obtain pure compounds for crystallization and identified crystals in multiple fractions that diffracted. Additional efforts to increase the throughput of high-quality crystal generation for MicroED, such as crystallization screening and crystallization chaperone exploration, will be needed to further advance root exudate metabolite identification. The overall optimized sample preparation process can then be integrated with the existing data collection and data analysis pipelines for MicroED at PNNL to facilitate more rapid natural product discovery.

59 BASIC BIOLOGICAL SCIENCES↗

Advancing Our Understanding of Pyranopterin-Dithiolene Contributions to Moco Enzyme Catalysis

The pyranopterin dithiolene ligand is remarkable in terms of its geometric and electronic structure and is uniquely found in mononuclear molybdenum and tungsten enzymes. The pyranopterin dithiolene is found coordinated to the metal ion, deeply buried within the protein, and non-covalently attached to the protein via an extensive hydrogen bonding network that is enzyme-specific. However, the function of pyranopterin dithiolene in enzymatic catalysis has been difficult to determine. This focused account aims to provide an overview of what has been learned from the study of pyranopterin dithiolene model complexes of molybdenum and how these results relate to the enzyme systems. This work begins with a summary of what is known about the pyranopterin dithiolene ligand in the enzymes. We then introduce the development of inorganic small molecule complexes that model aspects of a coordinated pyranopterin dithiolene and discuss the results of detailed physical studies of the models by electronic absorption, resonance Raman, X-ray absorption and NMR spectroscopies, cyclic voltammetry, X-ray crystallography, and chemical reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient fully-coherent quantum signal processing algorithms for real-time dynamics simulation

Simulating the unitary dynamics of a quantum system is a fundamental problem of quantum mechanics, in which quantum computers are believed to have significant advantage over their classical counterparts. One prominent such instance is the simulation of electronic dynamics, which plays an essential role in chemical reactions, non-equilibrium dynamics, and material design. These systems are time-dependent, which requires that the corresponding simulation algorithm can be successfully concatenated with itself over different time intervals to reproduce the overall coherent quantum dynamics of the system. In this paper, we quantify such simulation algorithms by the property of being fully-coherent: the algorithm succeeds with arbitrarily high success probability 1 − δ while only requiring a single copy of the initial state. Here we subsequently develop fully-coherent simulation algorithms based on quantum signal processing (QSP), including a novel algorithm that circumvents the use of amplitude amplification while also achieving a query complexity additive in time t, ln(1/δ), and ln(1/ϵ) for error tolerance ϵ: $Θ‖\mathscr{H}‖|t|+ln(1/ϵ)+ln(1/δ)$. Furthermore, we numerically analyze these algorithms by applying them to the simulation of the spin dynamics of the Heisenberg model and the correlated electronic dynamics of an H2 molecule. Since any electronic Hamiltonian can be mapped to a spin Hamiltonian, our algorithm can efficiently simulate time-dependent ab initio electronic dynamics in the circuit model of quantum computation. Accordingly, it is also our hope that the present work serves as a bridge between QSP-based quantum algorithms and chemical dynamics, stimulating a cross-fertilization between these exciting fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Catalytic Resonance Theory: The Catalytic Mechanics of Programmable Ratchets

Catalytic reaction networks of multiple elementary steps operating under dynamic conditions via a programmed input oscillation are difficult to interpret and optimize due to reaction system complexity. To understand these dynamic systems, individual elementary catalytic reactions oscillating between catalyst states were evaluated to identify their three fundamental characteristics that define their ability to promote reactions away from equilibrium. First, elementary catalytic reactions exhibit directionality to promote reactions forward or backward from equilibrium as determined by a ratchet directionality metric comprised of the input oscillation duty cycle and the reaction rate constants. Second, catalytic ratchets are defined by the catalyst state of strong or weak binding that permit reactants to proceed through the transition state. Third, elementary catalytic ratchets exhibit a cutoff frequency which defines the transition in applied frequency for which the catalytic ratchet functions to promote chemistry away from equilibrium. All three ratchet characteristics are calculated from chemical reaction parameters including rate constants derived from linear scaling parameters, reaction conditions, and catalyst state. The characteristics of the reaction network’s constituent elementary catalytic reactions provided an interpretation of complex reaction networks and a method of predicting the behavior of dynamic surface chemistry on oscillating catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stakeholder-Engaged Structured Decision Making for the Los Alamos Legacy Cleanup Mission - 20501

The Los Alamos National Laboratory (LANL) environmental legacy cleanup program requires decisions to be made for environmental remediation, decommissioning and disposal or management of radioactive waste. This legacy cleanup program was established to address groundwater contamination, material disposal areas (MDAs) that have been used to dispose of radioactive and other waste material, and 'aggregate areas' that might produce radioactive or other chemical waste as a result of remediation activities. The LANL site is regulated for environmental concerns under the Resource Conservation and Recovery Act (RCRA). However, some parts of LANL, such as material disposal area G (MDA G), have disposed of radioactive waste under DOE Order 435.1, and are subject to other regulations. For example, decommissioning the remote-handled TRU material stored in 33 shafts at MDA G falls under EPA's 40 CFR 191. Collectively, the regulations are all aimed in the same direction of finding the best solution, either through constructs such as 'as low as reasonably achievable' (ALARA), considering balancing factors as opposed to only cost and human health risk, and, under EPA regulations such as RCRA and NEPA, evaluating impact from all chemicals and both human health and ecological endpoints. Despite the basic goals and objectives of the regulations or their guidance, the main challenge is in their implementation. Arguably perhaps, but really in principle, all of these (and similar) regulations are asking for a decision analysis to be performed. Implementation challenges encountered have included lack of understanding of decision analysis in the industry, lack of effective stakeholder engagement in the decision analysis process, and lack of appreciation of the need to separate value judgments from science, the latter leading to developing conservative, or protective, science-based models. Conservative models lead to poor solutions, lack of effective stakeholder engagement leads to long drawn out protracted approaches to finding a solution (which also might never be found with this approach), and lack of understanding of decision analysis and Bayesian statistics causes poor models to be built, which creates unfortunate situations of 'garbage in, garbage out' becoming the basis for decision making. Stakeholder-engaged structured decision making (SDM) is an approach to solving problems that relies on the theory of decision science to involve stakeholders in the decision-making process. This approach incorporates stakeholder values using a scientifically rigorous methodology that separates value judgments from science in a way that helps avoid the pitfalls of biased, protective, or conservative modeling. This approach has its foundation in Keeney's 1992 treatise on value-focused thinking [1]. Keeney advocated a paradigm shift in decision making based on the idea that the standard way of thinking about decisions is backwards. The standard approach of focusing first on identifying alternatives rather than on articulating values results in a reactive approach with the emphasis on mechanics and fixed choices instead of the core values that have meaning to stakeholders. This paradigm shift effectively engages all stakeholders in the decision-making process while using a values focused thinking approach that can lead to the identification of decision opportunities and the creation of better alternatives. The intent is to be proactive and generate solutions that are related directly to values and objectives as identified by stakeholders. There are, perhaps, two overarching reasons why SDM can be used to benefit LANL's environmental legacy cleanup. Some of LANL's remaining waste management and environmental management problems are challenging and complex (for example, the Cr and RDX plumes, and MDAs) and while the traditional approach has, arguably, worked well for relatively simple risk-based problems, it cannot, or should not, be applied to more complex problems if the most effective and efficient solutions are desired. The second reason is cost. This has perhaps become more critical since publication of the Government Accountability Office (GAO) reports that DoE's environmental liability is considered a high-risk concern for the nation [2]. The focus of SDM is on structuring solutions to decision risk problems by first addressing stakeholder and decision maker values and subsequently developing decision objectives and ways to measure those objectives, preference weighting across objectives, identifying decision alternatives that best achieve those values, and characterizing uncertainty in predictions of the measures. Because a complete decision model is created using SDM, it can be evaluated to find the main elements of the model that drive, or predict, the best solution. This approach creates complete decision models that are transparent, traceable, reproducible and technically defensible. The science behind SDM, or decision analysis, is well founded, yet it is not unusual to see ad hoc approaches to decision making implemented under various environmental regulations that are pertinent to the LANL site, including NEPA, RCRA and DOE Order 435.1. Such ad hoc approaches are often not transparent or traceable, and lack reproducibility and technical defensibility. The LANL legacy cleanup program has embarked on using SDM to address the complex problems that remain. Stakeholder meetings have been held, and a prototype version of the stakeholder value system has been developed. Further meetings are expected in the future to address specific project needs. This is a long-term endeavor considering the complex environmental problems faced by DOE EM in Los Alamos (EM-LA), and careful planning, consideration of stakeholder value systems, and engagement with stakeholders throughout the SDM process is expected to lead to a successful endpoint. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Crystallization pathways and interfacial drivers for the formation of hierarchical architectures

The development of structural hierarchy on various length scales during crystallization process is ubiquitous in biological systems and is also observed in synthetic nanomaterials. The driving forces for the formations of complex architectures range from local interfacial interactions, that modify interfacial speciation, local supersaturation, and nucleation barriers, to macroscopic interparticle forces. Although it is enticing to interpret the formation of hierarchical architectures as the assembly of independently nucleated building blocks, often crystallization pathways follow monomer-by-monomer addition with structural complexity arising from interfacial chemical coupling and strongly correlated fluctuation dynamics in the electric double layers. Here, the mechanism of the development of structural hierarchy through heterogeneous nucleation, coupled interfacial nucleation and assembly, and oriented attachment of independently nucleated particles is discussed. In this article, the emphasis is made on the discussion of the underlying interfacial forces and chemical coupling that drives crystallization pathways towards the formation of structural hierarchy.

36 MATERIALS SCIENCE↗

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗

H 2 Evolution at an Electrochemical “Underpotential” with an Iridium-Based Molecular Photoelectrocatalyst

Molecular platforms capable of both photochemical H 2 evolution and electrochemical catalyst regeneration hold promise for the efficient generation of solar fuels. Systems studied to date require a significant electrochemical overpotential (η), so no photon energy is being chemically stored. A detailed study of the parameters that govern the electrochemical overpotential of molecular photoelectrocatalysts is presented. By understanding how tuning the solvent, catalyst structure, and acid pK a affect the kinetics and thermodynamics of the potential- and light-driven H 2 evolution reaction, conditions were identified wherein the solvento complex [Cp*Ir(bpy)(NCCH 3 )][PF 6 ] 2 (Cp* is pentamethylcyclopentadienyl; bpy is 2,2'-bipyridine) mediates photoelectrocatalytic H 2 formation at an electrochemical “underpotential” (η < 0) driven by visible light. Under 460 nm illumination and in the presence of the weak organic acid H-PhTMG + (2-phenyl-1,1,3,3- tetramethylguanidinium), the catalyst facilitates H 2 evolution with Faradaic efficiencies up to 90% and overpotentials as low as –90 mV. As a result, mechanistic studies enabled the construction of graphical “maps” that can guide the choice of catalyst and conditions to minimize overpotential, applicable to both photoelectrocatalysts and dark electrocatalysts.

14 SOLAR ENERGY↗

TiO 2 Nanoparticles Enabling Photocatalytic Desulfurization for C–C Coupling Reaction Using Visible Light

We uncover a surface-mediated mechanism for visible-light-driven photocatalytic desulfurization of thiols on titanium dioxide (TiO 2 ) nanoparticles under LED illumination (>420 nm). Combined diffuse reflectance, infrared, and cryogenic electron paramagnetic resonance spectroscopy verify the formation of a surface thiol–Ti charge-transfer complex that introduces visible-light absorption into an otherwise UV-active reaction system. Selective photoexcitation of this surface complex drives interfacial charge transfer, directly activating the adsorbed thiol and inducing C–S bond scission without the band gap excitation of TiO 2 . The resulting carbon-centered radical intermediates are intercepted by styrene to form C–C coupling products. These results uncover a distinct surface-mediated photocatalytic mechanism in which charge-transfer complex formation governs both light absorption and chemical activation, providing a mechanistic framework for extending visible-light reactivity to wide-bandgap metal oxide photocatalysts.

addition of carbon radical to vinyl groups↗

Chemical tuning of a honeycomb magnet through a critical point

BaCo 2⁢ (AsO 4 )⁢ 2 (BCAO) has seen extensive study since its initial identification as a proximate Kitaev quantum spin liquid candidate. Thought to be described by the highly anisotropic XXZ-J 1 -J 3 model, the ease with which magnetic order is suppressed in the system indicates proximity to a spin liquid phase. Upon chemical tuning via partial arsenic substitution with vanadium, we show an initial suppression of long-range incommensurate order in the BCAO system to T ≈ 3.0 K, followed by increased spin freezing at higher substitution levels. Between these two regions, at around 10% substitution, the system is shown to pass through a critical point where the competing J 1 /J 3 exchange interactions become more balanced, producing a more complex magnetic ground state, likely stabilized by quantum fluctuations. Here, this state shows how slight compositional change in magnetically frustrated systems may be leveraged to tune ground state degeneracies and potentially realize a quantum spin liquid state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pele: An Exascale-Ready Suite of Combustion Codes

High fidelity simulations of realistic combustion devices are extremely demanding computationally because of the requirements to capture complex fuel chemical decomposition, its intricate interactions with turbulent, often multiphase, flows, and the wide separation of space and time scales between the thin flame and the device boundaries. Software required to carry out such computations tends to be extremely complex, particularly when designed to exploit hardware accelerators, and can be difficult to port and maintain. We present Pele, a performance portable suite of tools for the simulation of combustion systems, including codes to evolve reactive multiphase configurations in the low Mach number and compressible flow regimes, along with a set of inter-compatible post processing and in situ analysis tools. The Pele suite of tools is built on top of the AMReX framework for block-structured adaptive mesh refinement, which provides efficient data structures and algorithms that enable the development of a wide variety of efficient mesh and particle based PDE integration schemes. A hierarchical MPI+X parallelism scheme supports CPU-only and accelerated architectures, where X can be OpenMP, CUDA, and HIP based approaches for intra-node computational work distribution. The algorithms and data structures underlying the Pele simulation and analysis tools are highly scalable and performant across a wide variety of high-performance computing platforms, including DOEs newest exascale-class machines, Frontier and Aurora. The simulation and analysis tools are fully documented and freely distributed as open source via GitHub. We present key algorithmic and software challenges, solution strategies, performance and resulting set of capabilities.

AMReX↗

Local Spin Density Approximation Strongly Improved by a Better-Informed Local Scaling of Its Self-Interaction Correction

The Perdew−Zunger self-interaction correction (PZSIC) makes density functional approximations (DFAs) exact for all one-electron densities. However, it overcorrects in manyelectron regions, introducing errors for the uniform-density limit, where uncorrected DFAs are exact. The locally scaled PZSIC (LSIC), based on the iso-orbital indicator zσ [which distinguishes single-orbital and slowly varying density regions and is used with the local spin density approximation (LSDA)], restores the uniform-density limit and significantly improves results for many properties, including chemical reaction barrier heights, atomization energies, and ionization potentials. Yet, LSIC performs poorly for weakly bonded systems, leaving many unbound, due to limitations of its iso-orbital indicator. To correct this, in this work we propose a new local scaling, LSIC-α, based on the iso-orbital indicator ασ (which additionally identifies regions of overlapping density tails). A two-parameter scaling function of ασ is fitted to a subset of the nonbonded appropriate norms for the SCAN and r2SCAN meta- GGAs, and tested on many properties of main-group atoms, molecules, and molecular complexes. LSIC-α greatly improves the interaction energies of weakly bonded systems in the S22 data set while retaining LSIC’s accuracy for other properties. This work shows that the errors of LSDA (and presumably of higher-level DFAs) can be largely but not entirely repaired by a proper “do no harm” self-interaction correction.

Approximation↗

A First-Principles Approach to Modeling Surface Site Stabilities on Multimetallic Catalysts

The study of multimetallic alloys and the multitude of possible surface compositions have sparked a tremendous interest in engineering low-cost materials with high activity and selectivity in heterogeneous catalysis. Multimetallic systems provide complementary functionalities and an unprecedented tunability when designing catalyst formulations. However, due to their immense structural and compositional complexity, the investigation and identification of an optimal catalyst is a tedious and time-consuming process, both experimentally and theoretically. Therefore, theoretical design principles are highly desirable to accelerate the screening of catalyst structures across the vast compositional space. In this paper, we introduce a simple and general model for predicting the site stability of multimetallic surfaces and nanoparticles, which is based on physical principles. The model requires only a small set of density functional theory (DFT) calculations of metal atom binding energies on monometallic and dilute alloy surface slabs to optimize the parameters in the simple model. The resulting model allows for the quantification of the stability of any particular atom site in any conceivable chemical environment across a wide range of morphologies, sizes, and arrangements by interpolating the derived parameters from a monometallic system to a completely diluted alloyed system. Herein, we demonstrate the robustness of the model across an extensive data set of transition metal alloy surfaces and 147-atoms cuboctahedral nanoparticles (NPs) composed of IrRhRu and PtPdRu. In conclusion, our approach yields mean absolute errors of ≈0.15 (IrRhRu), 0.20 (PtPdRu), 0.19 (IrRhRu NP), and 0.26 (PtPdRu NP) eV relative to site binding energies calculated using DFT.

36 MATERIALS SCIENCE↗

Tailoring Molecular Space to Navigate Phase Complexity in Cs-Based Quasi-2D Perovskites via Gated-Gaussian-Driven High-Throughput Discovery

Cesium-based quasi-2D halide perovskites (HPs) offer promising functionalities and low-temperature manufacturability, suited to stable tandem photovoltaics. However, the chemical interplays between the molecular spacers and the inorganic building blocks during crystallization cause substantial phase complexities in the resulting matrices. To successfully optimize and implement the quasi-2D HP functionalities, a systematic understanding of spacer chemistry, along with the seamless navigation of the inherently discrete molecular space, is necessary. Herein, by utilizing high-throughput automated experimentation, the phase complexities in the molecular space of quasi-2D HPs are explored, thus identifying the chemical roles of the spacer cations on the synthesis and functionalities of the complex materials. Furthermore, a novel active machine learning algorithm leveraging a two-stage decision-making process, called gated Gaussian process Bayesian optimization is introduced, to navigate the discrete ternary chemical space defined with two distinctive spacer molecules. Through simultaneous optimization of photoluminescence intensity and stability that “tailors” the chemistry in the molecular space, a ternary-compositional quasi-2D HP film realizing excellent optoelectronic functionalities is demonstrated. Finally, this work not only provides a pathway for the rational and bespoke design of complex HP materials but also sets the stage for accelerated materials discovery in other multifunctional systems.

36 MATERIALS SCIENCE↗