Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global Optimization”

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 235 records · Page 13

Accelerating the Structure Exploration of Diverse Bi–Pt Nanoclusters via Physics‐Informed Machine Learning Potential and Particle Swarm Optimization

Bimetallic Bi–Pt nanoclusters exhibit diverse structural motifs, including core-shell, Janus, and mixed alloy configurations, due to the unique bonding characteristics between Bi and Pt atoms. Using density functional theory refinements from ChIMES physically machine-learned potential and CALYPSO particle swarm optimization global searches, 34 Bi20-Pt20 nanoclusters are systematically classified. The results reveal that Bi atoms predominantly occupy surface sites, driven by charge transfer effects. Cohesive energy trends alone prove insufficient for structure differentiation, necessitating a data-driven approach employing principal component analysis and K-means clustering. Furthermore, vibrational, electronic, and infrared spectral analyses provide additional insights into structure-property relationships. The findings offer an original framework for the automated classification and analysis of bimetallic nanoclusters, enhancing the understanding of their stability and functional properties.

bimetallic nanoparticles↗

Bulk‐Boundary Correspondence of Semimetal Ru 3 Sn 7 and Topological Surface States on Chemically Realistic Terminations

Ru 3 Sn 7 is experimentally demonstrated as an efficient catalyst, with potential utilization of topological surface states for hydrogen evolution reaction. Despite its promising catalytic performance, the topological nature of Ru 3 Sn 7 remains uncertain. Particularly, the bulk-boundary correspondence has not yet been established, hence hindering a rigorous justification of its topologically-protected surface states. In this work, the bulk topology of Ru 3 Sn 7 is detailed using first-principles calculations and the topological quantum chemistry formalism. Ru3Sn7 turns out to be an enforced semimetal possessing symmetry-protected crossings within a set of bands near the Fermi level, which are enforced and prescribed by the violations of symmetry-prescribed compatibility relations. Moreover, the surface states and the associated origin from the same set of entangled bands are identified, thereby establishing the bulk-boundary correspondence. To evaluate the effects of chemical modifications, the response of topological surface states to various surface terminations, stoichiometry, and oxidation is examined. The surface structures are globally optimized, and the phase diagrams for various experimental conditions are built. It is shown that, due to changes in the local chemical environment, the original surface states are significantly altered. Modified surface bands can be observed near the Fermi level on surface terminations that preserve the C 4v symmetry.

36 MATERIALS SCIENCE↗

Parameter estimation from spontaneous imbibition into volcanic tuff

Two-phase fluid flow properties underlie quantitative prediction of water and gas movement, but constraining these properties typically requires multiple time-consuming laboratory methods. The estimation of two-phase flow properties (van Genuchten parameters, porosity, and intrinsic permeability) is illustrated in cores of vitric nonwelded volcanic tuff using Bayesian parameter estimation that fits numerical models to observations from spontaneous imbibition experiments. The uniqueness and correlation of the estimated parameters is explored using different modeling assumptions and subsets of the observed data. The resulting estimation process is sensitive to both moisture retention and relative permeability functions, thereby offering a comprehensive method for constraining both functions. The data collected during this relatively simple laboratory experiment, used in conjunction with a numerical model and a global optimizer, result in a viable approach for augmenting more traditional capillary pressure data obtained from hanging water column, membrane plate extractor, or mercury intrusion methods. This method may be useful when imbibition rather than drainage parameters are sought, when larger samples (e.g., including heterogeneity or fractures) need to be tested that cannot be accommodated in more traditional methods, or when in educational laboratory settings.

58 GEOSCIENCES↗

McCormick envelopes in mixed-integer PDE-constrained optimization

McCormick envelopes are a standard tool for deriving convex relaxations of optimization problems that involve polynomial terms. Such McCormick relaxations provide lower bounds, for example, in branch-and-bound procedures for mixed-integer nonlinear programs but have not gained much attention in PDE-constrained optimization so far. This lack of attention may be due to the distributed nature of such problems, which on the one hand leads to infinitely many linear constraints (generally state constraints that may be difficult to handle) in addition to the state equation for a pointwise formulation of the McCormick envelopes and renders bound-tightening procedures that successively improve the resulting convex relaxations computationally intractable. We analyze McCormick envelopes for a model problem class that is governed by a semilinear PDE involving a bilinearity and integrality constraints. We approximate the nonlinearity and in turn the McCormick envelopes by averaging the involved terms over the cells of a partition of the computational domain on which the PDE is defined. This yields convex relaxations that underestimate the original problem up to an a priori error estimate that depends on the mesh size of the discretization. These approximate McCormick relaxations can be improved by means of an optimization-based bound-tightening procedure. We show that their minimizers converge to minimizers to a limit problem with a pointwise formulation of the McCormick envelopes when driving the mesh size to zero. We provide a computational example, for which we certify all of our imposed assumptions. The results point to both the potential of the methodology and the gaps in the research that need to be closed. Our methodology provides a framework first for obtaining pointwise underestimators for nonconvexities and second for approximating them with finitely many linear inequalities in an infinite-dimensional setting.

Approximations and Expansions↗

Revisiting point defect thermodynamics in group IVB and VB transition metal carbides

We present a comprehensive re-examination of point defect thermodynamics in group IVB and VB transition metal carbides (TMCs) with the rocksalt structure using a combination of density functional theory (DFT) calculations and a statistical mechanical Wagner-Schottky model within the canonical ensemble. The most stable configurations of point defects were discovered using basin-hopping global optimization, driven by either a machine learning interatomic potential (MLIP) or DFT. A key finding is the identification of previously unreported dicarbon antisites—a C–C dimer occupying a metal site—as the structural (constitutional) defects on the carbon-rich side of stoichiometry in all group IVB and VB TMCs except TaC. Furthermore, dicarbon antisite-containing thermal defect complexes, such as quadruple and interbranch defects, can dominate in TMCs under specific stoichiometric and temperature conditions. In conclusion, by incorporating dicarbon antisites into the defect landscape, this work provides a revised understanding of the thermodynamics of point defects in TMCs.

Carbides↗

Hyperplane decision trees as piecewise linear surrogate models for chemical process design

Recent trends in chemical engineering research point towards an increasing reliance on data-driven modeling approaches. Neural networks, for instance, have proven to be accurate when data is plentiful and high-dimensional, but in many cases, they require computationally-intensive training procedures. Here, in this work, we describe hyperplane decision trees (HT) as a highly expressive and low-compute machine learning model architecture. These models are locally linear and have linear decision boundaries, resulting in a piecewise linear model of the data. This property allows them to be converted into mixed-integer linear constraints which can be globally optimized. Our open-source PyTorch implementation of this method is a fast, flexible, and accessible way to build accurate piecewise linear models of data.

Decision trees↗

Designing a drone delivery network with automated battery swapping machines

Drones are projected to alter last-mile delivery, but their short travel range is a concern. In this study, we propose a drone delivery network design using automated battery swapping machines (ABSMs) to extend ranges. The design minimizes the long-term delivery costs, including ABSM investment, drone ownership, and cost of the delivery time, and locates ABSMs to serve a set of customers. We build a mixed-integer nonlinear program that captures the nonlinear waiting time of drones at ABSMs. To solve the problem, we create an exact solution algorithm that finds the globally optimal solution using a derivative-supported cutting-plane method. To validate the applicability of our program, we conduct a case study on the Chicago Metropolitan area using cost data from leading ABSM manufacturer and geographical data from the planning and operations language for agent-based regional integrated simulation (more commonly known as POLARIS). A sensitivity analysis identifies that ABSM service times and costs are the key parameters impacting the long-term adoption of drone delivery.

25 ENERGY STORAGE↗

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗

Navigating Transition-Metal Chemical Space: Artificial Intelligence for First-Principles Design

Conspectus The variability of chemical bonding in open-shell transition-metal complexes not only motivates their study as functional materials and catalysts but also challenges conventional computational modeling tools. Here, tailoring ligand chemistry can alter preferred spin or oxidation states as well as electronic structure properties and reactivity, creating vast regions of chemical space to explore when designing new materials atom by atom. Although first-principles density functional theory (DFT) remains the workhorse of computational chemistry in mechanism deduction and property prediction, it is of limited use here. DFT is both far too computationally costly for widespread exploration of transition-metal chemical space and also prone to inaccuracies that limit its predictive performance for localized d electrons in transition-metal complexes. These challenges starkly contrast with the well-trodden regions of small-organic-molecule chemical space, where the analytical forms of molecular mechanics force fields and semiempirical theories have for decades accelerated the discovery of new molecules, accurate DFT functional performance has been demonstrated, and gold-standard methods from correlated wavefunction theory can predict experimental results to chemical accuracy. The combined promise of transition-metal chemical space exploration and lack of established tools has mandated a distinct approach. In this Account, we outline the path we charted in exploration of transition-metal chemical space starting from the first machine learning (ML) models (i.e., artificial neural network and kernel ridge regression) and representations for the prediction of open-shell transition-metal complex properties. The distinct importance of the immediate coordination environment of the metal center as well as the lack of low-level methods to accurately predict structural properties in this coordination environment first motivated and then benefited from these ML models and representations. Once developed, the recipe for prediction of geometric, spin state, and redox potential properties was straightforwardly extended to a diverse range of other properties, including in catalysis, computational “feasibility”, and the gas separation properties of periodic metal–organic frameworks. Interpretation of selected features most important for model prediction revealed new ways to encapsulate design rules and confirmed that models were robustly mapping essential structure–property relationships. Encountering the special challenge of ensuring that good model performance could generalize to new discovery targets motivated investigation of how to best carry out model uncertainty quantification. Distance-based approaches, whether in model latent space or in carefully engineered feature space, provided intuitive measures of the domain of applicability. With all of these pieces together, ML can be harnessed as an engine to tackle the large-scale exploration of transition-metal chemical space needed to satisfy multiple objectives using efficient global optimization methods. In practical terms, bringing these artificial intelligence tools to bear on the problems of transition-metal chemical space exploration has resulted in ML-model assessments of large, multimillion compound spaces in minutes and validated new design leads in weeks instead of decades.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Functionalization of Electrodes with Tunable [EMIM] x [Cl] x +1 – Ionic Liquid Clusters for Electrochemical Separations

Functionalization of electrodes with clusters of hydrophobic molecules may improve the energy efficiency and selectivity of electrochemical separations by modulating the desolvation process occurring at the interface. Ionic liquids (IL), which have a broad range of compositions and properties, are potential candidates for controlling the transport, desolvation, and adsorption of target ions at electrochemical interfaces. We report a joint experimental and theoretical investigation of the structure, stability, and selective adsorption properties of the IL clusters 1-ethyl-3-methylimidazolium chloride [EMIM] x [Cl] x+1 - (x = 1 – 10) and demonstrate their ability to adsorb and separate ions from solution. The structure and stability of the IL clusters are determined experimentally using high-mass-resolution electrospray ionization mass spectrometry, collision-induced dissociation, and negative ion photoelectron spectroscopy. Global optimization theory and ab initio molecular dynamics simulations provide molecular-level insight into the bonding and structural fluxionality of these species. Ion soft landing is used to selectively functionalize the surface of highly oriented pyrolytic graphite (HOPG) working electrodes with [EMIM] 1 [Cl] 2 - , [EMIM] 3 [Cl] 4 - , and [EMIM] 5 [Cl] 6 - clusters. Kelvin probe microscopy provides insight into the relative stability of the clusters on HOPG and their effect on the work function of IL-functionalized electrodes. Cyclic voltammetry measurements reveal irreversible adsorption of Fe(CN) 6 4-/3- anions during redox cycling, while electrochemical impedance spectroscopy indicates a substantial decrease in the electron transfer resistance of the IL-functionalized electrodes due to adsorption of Fe(CN) 6 4-/3- . Overall, our findings demonstrate that IL clusters with different size and stoichiometry may be used to increase the efficiency of electrochemical separations, opening new horizons in selective electrode functionalization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interpreting the Operando X-ray Absorption Near-Edge Structure of Supported Cu and CuPd Clusters in Conditions of Oxidative Dehydrogenation of Propane: Dynamic Changes in Composition and Size

Supported subnano-cluster catalysts are highly dynamic, developing true active sites only under the pressures and temperatures of reaction conditions. Operando X-ray absorption near-edge structure (XANES) spectroscopy can track changes in the oxidation state and the local environment of cluster atoms, providing insight into the development of these active sites. While bulk metal, oxide, and hydroxide standards are often used for fitting experimental XANES spectra to obtain average oxidation states, we recently showed that computed cluster standards of relevant compositions are a more suitable basis, producing more accurate fits. Here, we theoretically interpret the operando XANES of supported Cu 3 Pd and Cu 4 clusters during temperature-programmed reaction (TPRx) of oxidative dehydrogenation of propane. We use an expanded basis set including both globally optimized computed clusters and bulk standards. Not only can we track reversible composition/oxidation state change with temperature, but also the irreversible growth of the bulk fraction upon heating, which we attribute to cluster sintering. This has important implications for the mechanism of the catalyzed reaction and the nature of the available active sites. Here, we propose that operando XANES provides most significant insight into the nature of supported cluster catalysts in reaction conditions when interpreted using mixed computed cluster and bulk standards.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanism of Stoichiometrically Governed Titanium Oxide Brownian Tree Formation on Stepped Au(111)

Previously observed formation of substoichiometric titanium oxide dendritic structures across terraces of Au(111) is computationally studied and shown to follow the classical fractal formation mechanism of diffusion-limited aggregation (DLA). Globally optimized gas-phase oxide cluster structures are sampled in a variety of landing formations on gold surfaces and shown to favor isomers driving polymerization to Brownian tree fractal structures. Mobility of Ti 3 O 5 monomers is shown to be extremely high, with diffusion barriers of 0.21 eV or less. Through bonding stabilization, polymerization of these monomers is energetically favorable and irreversible on the 111 terrace but geometrically impossible to propagate along the step edge. Simulated scanning tunneling microscopy (STM) images show strong similarity to experiment. By contrast, observation of Ti 3 O 6 aggregating as wires along step edges is explained by the affinity of oxygen to step edges and statistical arguments for aggregation entropy at the step, in addition to low barriers for monomer diffusion and polymerization.

36 MATERIALS SCIENCE↗

Structure and Stability of the Ionic Liquid Clusters [EMIM]n[BF4]n+1- (n = 1 – 9): Implications for Electrochemical Separations

Precise functionalization of electrodes with size-selected ionic liquid (IL) clusters may improve the application of ILs in electrochemical separations. Herein, we report our combined experimental and theoretical investigation of the IL clusters [EMIM]????[BF4]????+1- (n = 1 – 9) and demonstrate their selectivity and efficiency towards targeted adsorption of molecules from solution. The structures and energies of the IL clusters, predicted with global optimization, agree with and help interpret the ion abundances and stabilities measured by high mass resolution electrospray ionization mass spectrometry and collision induced dissociation experiments. The [EMIM][BF4]2- cluster, which was identified as the most stable IL cluster, was selectively soft landed onto a working electrode. Electrochemical impedance spectroscopy revealed a lower charge transfer resistance on the soft-landed electrode containing [EMIM] [BF_4 ]_2^- compared to an electrode prepared by drop casting of IL solution, which contains the full range of IL clusters. Our findings indicate that specific IL clusters may be used to increase the efficiency of electrochemical separations by lowering the overpotentials involved.

Zhang, Jun↗

CO (2) Electroreduction on Borated Copper Surfaces: Boron Active Sites, Not Copper

Boron-doped copper has recently emerged as an active and stable catalyst for the electrochemical reduction of CO 2 to value-added C 2 products. Here, in this work, we develop a realistic model of CO electroreduction on surface borides of copper under operational conditions, taking into account the effects of electrode potential, electrolyte environment, and pH. We study the possible reconstruction of the electrocatalyst surface using grand canonical DFT and global optimization to obtain a potential-dependent grand canonical ensemble description of metastable, hydrogen-covered catalyst surfaces. Two key surface configurations, low H coverage (LC, −0.6 V SHE ) and high H coverage (HC, −0.8 V SHE ), dominate this ensemble, with the former being kinetically persistent and C 2 selective under strongly reducing conditions. Nonmetallic boron sites on the surface copper boride are found to bind CO more strongly than copper sites, and mechanistic investigation of CO electroreduction pathways presents a surprisingly unconventional case of boron-centered reactivity in contrast to typical copper-centered reactivity. Neighboring boron sites present along boron chains on the surface copper boride are found to facilitate C−C coupling, thereby driving the high C 2 selectivity of this electrocatalyst.

boron-doped copper↗

Electrocatalytic Hydrogen Evolution at Full Atomic Utilization over ITO-Supported Sub-nano-Pt n Clusters: High, Size-Dependent Activity Controlled by Fluxional Pt Hydride Species

A combination of density functional theory (DFT) and experiments with atomically size-selected Pt n clusters deposited on indium-tin oxide (ITO) electrodes was used to examine the effects of applied potential and Pt n size on the electrocatalytic activity of Pt n (n = 1, 4, 7, 8) for the hydrogen evolution reaction (HER). Activity is found to be negligible for isolated Pt atoms on ITO, increasing rapidly with Pt n size, such that Pt 7 /ITO and Pt 8 /ITO have roughly double the activity per Pt atom compared to atoms in the surface layer of polycrystalline Pt. Both DFT and experiment find that hydrogen under-potential deposition (H upd ) results in Pt n /ITO (n = 4, 7, 8) adsorbing ~2 H atoms/Pt atom at the HER threshold potential, equal to ca. double the Hupd observed for Pt bulk or nanoparticles. Here, the cluster catalysts under electrocatalytic conditions are hence best described as a Pt hydride compound, significantly departing from a metallic Pt cluster. The exception is Pt 1 /ITO, where H adsorption at the HER threshold potential is energetically unfavorable. Theory combines global optimization with grand canonical approaches for the influence of potential, uncovering that several metastable structures contribute to HER, changing with the applied potential. It is hence critical to include reactions of the ensemble of energetically accessible Pt n H x /ITO structures to correctly predict the activity vs. Pt n size and applied potential. For the small clusters, spillover of H ads from the clusters to the ITO support is significant, resulting in a competing channel for loss of H ads , particularly at slow potential scan rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Off-Stoichiometric Restructuring and Sliding Dynamics of Hexagonal Boron Nitride Edges in Conditions of Oxidative Dehydrogenation of Propane

Boron-containing materials, such as hexagonal boron nitride (h-BN), recently shown to be active and selective catalysts for the oxidative dehydrogenation of propane (ODHP), have been shown to undergo significant surface oxyfunctionalization and restructuring. Although experimental ex situ studies have probed the change in chemical environment on the surface, the structural evolution of it under varying reaction conditions has not been established. Herein, we perform global optimization structure search with a grand canonical genetic algorithm to explore the chemical space of off-stoichiometric restructuring of the h-BN surface under ambient as well as ODHP-relevant conditions. A grand canonical ensemble representation of the surface is established, and the predicted 11B solid-state NMR spectra are consistent with previous experimental reports. In addition, we investigated the relative sliding of h-BN sheets and how it influences the surface chemistry with ab initio molecular dynamics simulations. Furthermore, the B–O linkages on the edges are found to be significantly strained during the sliding, causing the metastable sliding configurations to have higher reactivity toward the activation of propane and water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spontaneous dynamical disordering of borophenes in MgB2 and related metal borides

Abstract Layered boron compounds have attracted significant interest in applications from energy storage to electronic materials to device applications, owing in part to a diversity of surface properties tied to specific arrangements of boron atoms. Here we report the energy landscape for surface atomic configurations of MgB 2 by combining first-principles calculations, global optimization, material synthesis and characterization. We demonstrate that contrary to previous assumptions, multiple disordered reconstructions are thermodynamically preferred and kinetically accessible within exposed B surfaces in MgB 2 and other layered metal diborides at low boron chemical potentials. Such a dynamic environment and intrinsic disordering of the B surface atoms present new opportunities to realize a diverse set of 2D boron structures. We validated the predicted surface disorder by characterizing exfoliated boron-terminated MgB 2 nanosheets. We further discuss application-relevant implications, with a particular view towards understanding the impact of boron surface heterogeneity on hydrogen storage performance.

36 MATERIALS SCIENCE↗

Cobalt-based magnetic Weyl semimetals with high-thermodynamic stabilities

Abstract Recent experiments identified Co 3 Sn 2 S 2 as the first magnetic Weyl semimetal (MWSM). Using first-principles calculation with a global optimization approach, we explore the structural stabilities and topological electronic properties of cobalt (Co)-based shandite and alloys, Co 3 MM’X 2 (M/M’ = Ge, Sn, Pb, X = S, Se, Te), and identify stable structures with different Weyl phases. Using a tight-binding model, for the first time, we reveal that the physical origin of the nodal lines of a Co-based shandite structure is the interlayer coupling between Co atoms in different Kagome layers, while the number of Weyl points and their types are mainly governed by the interaction between Co and the metal atoms, Sn, Ge, and Pb. The Co 3 SnPbS 2 alloy exhibits two distinguished topological phases, depending on the relative positions of the Sn and Pb atoms: a three-dimensional quantum anomalous Hall metal, and a MWSM phase with anomalous Hall conductivity (~1290 Ω −1 cm −1 ) that is larger than that of Co 2 Sn 2 S 2 . Our work reveals the physical mechanism of the origination of Weyl fermions in Co-based shandite structures and proposes topological quantum states with high thermal stability.

36 MATERIALS SCIENCE↗