Search NASA⌕ Search

SEARCH · Search NASA

Results for “State Estimation”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Variational Quantum Fidelity Estimation

Computing quantum state fidelity will be important to verify and characterize states prepared on a quantum computer. In this work, we propose novel lower and upper bounds for the fidelity F ( ρ , σ ) based on the ``truncated fidelity'' F ( ρ m , σ ) , which is evaluated for a state ρ m obtained by projecting ρ onto its m -largest eigenvalues. Our bounds can be refined, i.e., they tighten monotonically with m . To compute our bounds, we introduce a hybrid quantum-classical algorithm, called Variational Quantum Fidelity Estimation, that involves three steps: (1) variationally diagonalize ρ , (2) compute matrix elements of σ in the eigenbasis of ρ , and (3) combine these matrix elements to compute our bounds. Our algorithm is aimed at the case where σ is arbitrary and ρ is low rank, which we call low-rank fidelity estimation, and we prove that no classical algorithm can efficiently solve this problem under reasonable assumptions. Finally, we demonstrate that our bounds can detect quantum phase transitions and are often tighter than previously known computable bounds for realistic situations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Large uncertainty in permafrost carbon stocks due to hillslope soil deposits

Here, northern circumpolar permafrost soils contain more than a third of the global Soil Organic Carbon pool (SOC). The sensitivity of this carbon pool to a changing climate is a primary source of uncertainty in simulationbased climate projections. These projections, however, do not account for the accumulation of soil deposits at the base of hillslopes (hill-toes), and the influence of this accumulation on the distribution, sequestration, and decomposition of SOC in landscapes affected by permafrost. Here we combine topographic models with soil-profile data and topographic analysis to evaluate the quantity and uncertainty of SOC mass stored in perennially frozen hill-toe soil deposits. We show that in Alaska this SOC mass introduces an uncertainty that is > 200% than state-wide estimates of SOC stocks (77 PgC), and that a similarly large uncertainty may also pertain at a circumpolar scale. Soil sampling and geophysical-imaging efforts that target hill-toe deposits can help constrain this large uncertainty.

58 GEOSCIENCES↗
DOE OSTI2017

Global Sensitivity of Simulated Water Balance Indicators Under Future Climate Change in the Colorado Basin

The Colorado River basin is a fundamentally important river for society, ecology and energy in the United States. Streamflow estimates are often provided using modeling tools which rely on uncertain parameters; sensitivity analysis can help determine which parameters impact model results. Despite the fact that simulated flows respond to changing climate and vegetation in the basin, parameter sensitivity of the simulations under climate change has rarely been considered. In this study, we conduct a global sensitivity analysis to relate changes in runoff, evapotranspiration, snow water equivalent and soil moisture to model parameters in the Variable Infiltration Capacity (VIC) hydrologic model. Here, we combine global sensitivity analysis with a space-filling Latin Hypercube sampling of the model parameter space and statistical emulation of the VIC model to examine sensitivities to uncertainties in 46 model parameters following a variance-based approach.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Memristive Behavior Enabled by Amorphous–Crystalline 2D Oxide Heterostructure

The emergence of memristive behavior in amorphous–crystalline 2D oxide heterostructures, which are synthesized by atomic layer deposition (ALD) of a few-nanometer amorphous Al 2 O 3 layers onto atomically thin single-crystalline ZnO nanosheets, is demonstrated. The conduction mechanism is identified based on classic oxygen vacancy conductive channels. ZnO nanosheets provide a 2D host for oxygen vacancies, while the amorphous Al 2 O 3 facilitates the generation and stabilization of the oxygen vacancies. The conduction mechanism in the high-resistance state follows Poole–Frenkel emission, and in the the low-resistance state is fitted by the Mott–Gurney law. From the slope of the fitting curve, the mobility in the low-resistance state is estimated to be ≈2400 cm 2 V -1 s -1 , which is the highest value reported in semiconductor oxides. When annealed at high temperature to eliminate oxygen vacancies, Al is doped into the ZnO nanosheet, and the memristive behavior disappears, further confirming the oxygen vacancies as being responsible for the memristive behavior. The 2D heterointerface offers opportunities for new design of high-performance memristor devices.

2D heterostructures↗
DOE OSTI2013

Radiative Heating of the ISCCP Upper Level Cloud Regimes and its Impact on the Large-scale Tropical Circulation

Radiative heating profiles of the International Satellite Cloud Climatology Project (ISCCP) cloud regimes (or weather states) were estimated by matching ISCCP observations with radiative properties derived from cloud radar and lidar measurements from the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) sites at Manus, Papua New Guinea, and Darwin, Australia. Focus was placed on the ISCCP cloud regimes containing the majority of upper level clouds in the tropics, i.e., mesoscale convective systems (MCSs), deep cumulonimbus with cirrus, mixed shallow and deep convection, and thin cirrus. At upper levels, these regimes have average maximum cloud occurrences ranging from 30% to 55% near 12 km with variations depending on the location and cloud regime. The resulting radiative heating profiles have maxima of approximately 1 K/day near 12 km, with equal heating contributions from the longwave and shortwave components. Upper level minima occur near 15 km, with the MCS regime showing the strongest cooling of 0.2 K/day and the thin cirrus showing no cooling. The gradient of upper level heating ranges from 0.2 to 0.4 K/(day∙km), with the most convectively active regimes (i.e., MCSs and deep cumulonimbus with cirrus) having the largest gradient. When the above heating profiles were applied to the 25-year ISCCP data set, the tropics-wide average profile has a radiative heating maximum of 0.45Kday-1 near 250 hPa. Column-integrated radiative heating of upper level cloud accounts for about 20% of the latent heating estimated by the Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR). The ISCCP radiative heating of tropical upper level cloud only slightly modifies the response of an idealized primitive equation model forced with the tropics-wide TRMM PR latent heating, which suggests that the impact of upper level cloud is more important to large-scale tropical circulation variations because of convective feedbacks rather than direct forcing by the cloud radiative heating profiles. However, the height of the radiative heating maxima and gradient of the heating profiles are important to determine the sign and patterns of the horizontal circulation anomaly driven by radiative heating at upper levels.

circulation↗
DOE OSTI2018

High-precision control of tracked field robots in the presence of unknown traction coefficients

Accurate steering through crop rows that avoids crop damage is one of the most important tasks for agricultural robots utilized in various field operations, such as monitoring, mechanical weeding, or spraying. In practice, varying soil conditions can result in off-track navigation due to unknown traction coefficients so that it can cause crop damage. To address this problem, this work presents the development, application, and experimental results of a real-time receding horizon estimation and control (RHEC) framework applied to a fully autonomous mobile robotic platform to increase its steering accuracy. Recent advances in cheap and fast microprocessors, as well as advances in solution methods for nonlinear optimization problems, have made nonlinear receding horizon control (RHC) and receding horizon estimation (RHE) methods suitable for field robots that require high-frequency (milliseconds) updates. A real-time RHEC framework is developed and applied to a fully autonomous mobile robotic platform designed by the authors for in-field phenotyping applications in sorghum fields. Nonlinear RHE is used to estimate constrained states and parameters, and nonlinear RHC is designed based on an adaptive system model that contains time-varying parameters. The capabilities of the real-time RHEC framework are verified experimentally, and the results show an accurate tracking performance on a bumpy and wet soil field. Finally, the mean values of the Euclidean error and required computation time of the RHEC framework are equal to 0.0423 m and 0.88 ms, respectively.

42 ENGINEERING↗
DOE OSTI2025

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Setting a baseline for Integrated Landscape Design: Cost and risk assessment in herbaceous feedstock supply chains

The literature on bioenergy has pointed out that financial risks and limited performance benchmarks are leading barriers to wide-spread industry expansion. The purpose of this paper is to document a baseline cost and risk assessment from the perspective of biorefinery management for two herbaceous feedstocks, switchgrass and corn stover, in an application of Integrated Landscape Design (ILD). Integrating energy crops into row crops has attracted considerable attention leading to a variety of ways to implement ILD. However, variation in production conditions, different methods of implementing unit operations, and technology performance variability are among many factors leading to financial risk facing biorefineries from feedstock supply based on ILD. Utilizing data from field-based trials on switchgrass and corn stover in Virginia and Iowa in the United States, we estimate feedstock logistics cost and characterize it by a distribution of potential outcomes. The results show that, across fields included in the study, logistics cost for switchgrass vary from 73 per Mg to 171 per Mg. For corn stover, the inter-field variation in cost ranges from 83 per Mg to 105 per Mg. Results for intra-field uncertainty show that the cost for switchgrass ranges from 58 per Mg to 74 per Mg, and for corn stover they range from 50 per Mg to 65 per Mg. Furthermore, these results inform a baseline from which improvements in ILD can be compared, and provide biorefinery management with a starting point for cost and risk estimates of feedstock from ILD.

09 BIOMASS FUELS↗
DOE OSTI2022

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗
DOE OSTI2018

Fundamental catalytic challenges to design improved biomass conversion technologies

During the past ten years, there has been significant interest and investment in the study of catalytic conversion of biomass-derived feedstocks into renewable fuels and chemicals. In the United States, an estimated $25 billion has been spent by venture capitalists, industry, and government agencies during this period of time to commercialize “renewable technologies” including solar energy, wind power, batteries, and biofuels (e.g., cellulosic ethanol). Four societal factors are driving these investments including: (1) the increased price of crude oil; (2) concerns about global warming; (3) the desire to improve rural economies where biomass is produced, and; (4) national goals to become energy self-sufficient. Furthermore, underlying these efforts is the realization that lignocellulosic biomass is the only realistic, near-term and non-food-competitive source of renewable organic carbon. To this end, legislative efforts, such as the US Renewable Fuel Standards, have been implemented to create subsidies, tax credits, mandates, and loan guarantees to help bring renewable fuel technologies to market. However, the representative body of industrial efforts to this end has faced significant challenges, with several startup companies having commercialized biomass conversion technologies but having struggled to reach commercial scale; in fact, many of these companies have filed for bankruptcy. This situation is likely due to the challenges associated with scaling up unproven pioneer processes, as well as neophyte investors not understanding the decade-long time frames and the sheer amount of funding that is often required to bring chemical process technologies to market. Nevertheless, several emerging catalytic technologies have either entered the market place, or are currently demonstrating their technologies in fully integrated pilot plants. Commercial and near commercial technologies for second generation biomass conversion technologies to-date include, among others: biomass-derived jet and diesel fuel from both waste vegetable oils and ethanol; small scale production of renewable jet and diesel from landfill gases via Fischer-Tropsch synthesis; catalytic conversion of carbohydrates into gasoline and aromatics; hydropyrolysis of biomass into gasoline and diesel, and; catalytic conversion of wood into aromatics. To be successful, biomass conversion technologies must ultimately be able to compete economically with petroleum technologies, which are already operating at large commercial scales and have been practiced for decades. To this end, various factors must be considered when evaluating the potential for new processes to compete with incumbent technologies; factors such as regional variations in feedstock quality and availability, government policy, subsidies and tax rates, and the proprietary positions of the ancillary technologies that might support the process. By any measure, however, a critical metric of the economic potential of a process is its efficiency with respect to the yield of products from raw materials, and this metric of performance is directly related to atom efficiency of the underlying chemical reactions. Therefore, while promising biomass conversion technologies continue to demonstrate progress towards commercialization, there remains an important need for the catalysis community to aid in this effort by: (1) designing more active, selective and stable catalysts, and (2) elucidating a more detailed understanding of the catalytic chemistries underlying these processes. Indeed, the study of catalytic biomass conversion has grown tremendously during the past decade; and new or emerging technologies that are in the laboratory stage have allowed for biomass to be converted into a wider variety of commodity chemicals and the full range of liquid fuels that are produced from petroleum. The objective of this perspective is to highlight some of the ongoing, fundamental challenges with respect to the catalytic conversion of biomass into renewable products, and to provide insight as to how the catalysis community can overcome several of these challenges. It is our belief that, as an international academic community of catalysis researchers, we can (and must) work together to address these challenges, and to drive progress toward the production of next-generation renewable chemicals and fuels from biomass.

09 BIOMASS FUELS↗
DOE OSTI2022

Kernel learning backward SDE filter for data assimilation

In this paper, we develop a kernel learning backward SDE filter method to estimate the state of a stochastic dynamical system based on its partial noisy observations. A system of forward backward stochastic differential equations is used to propagate the state of the target dynamical model, and Bayesian inference is applied to incorporate the observational information. Further, to characterize the dynamical model in the entire state space, we introduce a kernel learning method to learn a continuous global approximation for the conditional probability density function of the target state by using discrete approximated density values as training data. Numerical experiments demonstrate that the kernel learning backward SDE is highly effective.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Plasma-parameter dependence of ro-vibrational temperatures for H 2 in LHD divertor

We analyzed a thousand visible spectra of Fulcher-α band measured for divertor plasmas in Large Helical Device. With a coronal model and Bayesian inference, the population distribution of hydrogen molecule in the electronic ground state was estimated. The non-thermal population distribution was recovered with a two-temperature model, which has two sets of rotational and vibrational temperatures, as well as their mixture coefficient. The lower rotational temperature significantly changes according to the plasma parameters. Its nearly linear dependence on the electron density was found, which is consistent with previous works. The lower vibrational temperature also shows a small density dependence, as reported by a previous work. On the other hand, the higher rotational and vibrational temperatures as well as the mixture coefficient only show slight changes over the broad range of plasma parameters. These population parameters show a significant correlation; with higher electron density, all the temperatures and the fraction of the higher-temperature component increase simultaneously. In conclusion, this suggests that the electron impact plays an important role to determine the population distribution.

42 ENGINEERING↗
DOE OSTI2017

Sensitivity of polycrystal plasticity to slip system kinematic hardening laws for Al 7075-T6

The prediction of formation and early growth of microstructurally small fatigue cracks requires use of constitutive models that accurately estimate local states of stress, strain, and cyclic plastic strain. However, few research efforts have attempted to systematically consider the sensitivity of overall cyclic stress-strain hysteresis and higher order mean stress relaxation and plastic strain ratcheting responses introduced by the slip system back-stress formulation in crystal plasticity, even for face centered cubic (FCC) crystal systems. This paper explores the performance of two slip system level kinematic hardening models using a finite element crystal plasticity implementation as a User Material Subroutine (UMAT) within ABAQUS, with fully implicit numerical integration. The two kinematic hardening formulations aim to reproduce the cyclic deformation of polycrystalline Al 7075-T6 in terms of both macroscopic cyclic stress-strain hysteresis loop shape, as well as ratcheting and mean stress relaxation under strain- or stress-controlled loading with mean strain or stress, respectively. The first formulation is an Armstrong-Frederick type hardening-dynamic recovery law for evolution of the back stress. This approach is capable of reproducing observed deformation under completely reversed uniaxial loading conditions, but overpredicts the rate of cyclic ratcheting and associated mean stress relaxation. The second formulation corresponds to a multiple back stress Ohno-Wang type hardening law with nonlinear dynamic recovery. The adoption of this back stress evolution law greatly improves the capability to model experimental results for polycrystalline specimens subjected to cycling with mean stress or strain. As a result, the relation of such nonlinear dynamic recovery effects are related to slip system interactions with dislocation substructures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

California's harvested wood products: A time-dependent assessment of life cycle greenhouse gas emissions

Following life-cycle assessment (LCA) methodology, this study presents a state-level estimation of embodied carbon of wood products harvested in 2019 from California and subsequently processed, manufactured, transported, used, and disposed at the end-of-life (EoL). In a conventional static approach to LCA, all GHG emissions were aggregated and considered to occur at year 0 of the given time horizon (500 years in this study) and used a static characterization factor (CF). In dynamic LCA, GHG emissions occurring in different years were considered, and their global warming impact (GWI) was determined using a time-dependent CF over the selected time horizon of 500 years. Four scenarios were developed to examine the impact of EoL choices on GWI. It was found that dynamic GWI for all scenarios ranged from 0.27 to 0.93 million tonne CO₂e, which were 45–73 % lower than those estimated with static LCA approach, indicating that the static LCA approach could lead to an underestimation of the benefits of substituting wood for non-wood products, compared to those based on dynamic LCA approach. This analysis also demonstrated that the choice of EoL treatment option is a key factor affecting the estimated GWI as it directly determines the annual emission of GHGs released into atmosphere and subsequently their warming effect depending on the time harvested wood products (HWPs) spend in the horizon of assessment. Altogether, the dynamic LCA performed in this study enabled more robust interpretations of embodied carbon by including temporal boundaries associated with the HWPs life cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Physical Properties of Candidate X-ray Detector Material Rb 4 Ag 2 BiBr 9

Recently, metal halide perovskites have emerged as promising semiconductor candidates for sensitive X-ray photon detection due to their suitable band gap energies, excellent charge transport properties, and low material cost afforded by their low-temperature solution-processing preparation. Furthermore, we report an improved methodology for single crystal growth and thermal and electrical properties of a two-dimensional (2D) layered halide material Rb 4 Ag 2 BiBr 9 , which has been identified as a potential candidate for X-ray radiation detection applications. The measured heat capacity for Rb 4 Ag 2 BiBr 9 implies that there are no structural phase transitions upon cooling. Temperature dependence of thermal transport measurements further suggests remarkably low thermal conductivities of Rb 4 Ag 2 BiBr 9 that are comparable to the lowest reported in literature. The bulk crystal resistivity is determined to be 2.59 × 10 9 Ω·cm from the current–voltage (I–V) curve. Density of trap states is estimated to be ~10 10 cm –3 using the space-charge-limited-current measurements. The fabricated Rb 4 Ag 2 BiBr 9 -based X-ray detector shows good operational stability with no apparent current drift, which may be ascribed to the 2D crystal structure of Rb 4 Ag 2 BiBr 9 . Finally, by varying the X-ray tube current to change the corresponding dose rate, the Rb4Ag2BiBr9 X-ray detector sensitivity is determined to be 222.03 μC Gy –1 cm –2 (at an electric field of E = 24 V/mm).

36 MATERIALS SCIENCE↗
DOE OSTI2025

PySIDT: Subgraph Isomorphic Decision Trees for Molecular Property Prediction

Accurate molecular property prediction is important across all fields of chemistry. Deep neural networks (DNNs) have become increasingly popular due to their ability to train automatically, avoiding the incredibly tedious process of constructing and extending traditional property estimation schemes. However, DNNs require large amounts of training data, are challenging to interpret, require large amounts of memory to load even during inference, and have severe difficulties incorporating qualitative chemical knowledge, which are often desired for molecular property prediction tasks. Here, in this study, we present PySIDT (https://github.com/zadorlab/PySIDT), a software for training and running inference on Subgraph Isomorphic Decision Trees (SIDTs). SIDTs are graph-based decision trees made of nodes associated with molecular substructures. Inference is done by descending target molecular structures down the decision tree to nodes with matching subgraph isomorphic substructures and making predictions based on the final (most specific) nodes matched. SIDTs scale down well to dataset sizes much smaller than is feasible for DNNs. As trees of molecular substructures, SIDTs are inherently readable and easy to visualize, making them easy to analyze. They are also straightforward to extend and retrain, facilitate uncertainty estimation, and enable easy integration of expert knowledge. We demonstrate the SIDT approach discussing its application to a diverse range of molecular prediction tasks: rate coefficient estimation, diffusion coefficient estimation, thermochemistry estimation, transition state bond stretch prediction, p K a prediction, stability of molecular structures, stability of surface structures, and prediction of surface lateral interaction energetics. Additionally, we demonstrate the power of the SIDT algorithms in two direct learning curve vanilla comparisons with the popular DNN-based software Chemprop and the popular gradient boosted trees-based software XGBoost on enthalpy of formation and rate coefficient prediction tasks. In particular, in the enthalpy of formation case, vanilla PySIDT is able to outperform vanilla Chemprop and XGBoost across the full range of training/validation set sizes out to 11,560 data points.

Johnson, Matthew Sean [Sandia National Laboratorie↗
DOE OSTI2020

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗