Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Timing Analysis”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ultra-high Bandwidth, Ultra-high Dynamic Range X-Ray Shock Characterization via Photonic Integrated Circuit

Nuclear weapon component assessment tests at the Z Machine rely on accurate X-ray yield measurements for model validation, design and analysis of component survivability, and source optimization. Yield measurement devices currently used do not provide the requisite certainty to enable accurate and efficient data analysis of test results, which leads to an increase in Z-shots required for evaluation, longer device development times, and higher operation costs. Additionally, current X-ray flux detectors lack high temporal resolution at the tails. Similarly, shock measurement techniques used for component assessment often lack the ability to spatially resolve wave behavior thereby limiting the ability to measure shock propagation dynamics needed to design the next generation of ND components. Here we present the development of a photonic micro-calorimeter and shock sensor consisting of meter-long waveguide spirals that are optomechanically coupled to X-ray absorbing layers to characterize yield and shock propagation at nanosecond timescales.

42 ENGINEERING↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2020

In situ , in vivo , and in operando imaging and spectroscopy of liquids using microfluidics in vacuum

This review offers a succinct overview of the development of a vacuum-compatible microfluidic reactor system for analysis at the liquid vacuum interface (SALVI), and its diverse applications in in situ, in vivo, and in operando imaging of liquid surfaces as well as the air-liquid (a-l), liquid-liquid (l-l), and solid-liquid (s-l) interfaces in the past decade. SALVI is one of the first microfluidics-based reactors that has enabled direct analysis of volatile liquids in vacuum surface tools such as scanning electron microscopy (SEM) and time-of-flight secondary ion mass spectrometry (ToF-SIMS). Its integration into ambient and vacuum spectroscopy and microscopy is illustrated. Several applications are highlighted including (1) imaging nanoparticles in liquid using in situ SEM; (2) mapping the evolving l-l interface using in situ x-ray absorption spectroscopy and ToF-SIMS; (3) following complex a-l interfacial oxidation reaction products using in situ ToF-SIMS; (4) capturing biological interfaces of cells and microbes via in vivo multimodal and correlative imaging; and (5) monitoring the dynamic solid electrode and liquid electrolyte interface using in operando molecular imaging. Finally, outlook and recommendations are presented. Besides showing the holistic information volume obtained by real-time multiplexed imaging, this review intends to convey the importance of tool development in revolutionizing surface and interface analysis using vacuum platforms previously limited to solid surfaces. Microfluidics is manifested to be not limited to ambient conditions in many examples in this review. Moreover, fundamental interfacial phenomena underpinning mass and charge transfer can now be pursued in real time via innovated chemical imaging and spectroscopy.

42 ENGINEERING↗
DOE OSTI2016

Improving Mixed-phase Cloud Parameterization in Climate Model with the ACRF Measurements

Mixed-phase cloud microphysical and dynamical processes are still poorly understood, and their representation in GCMs is a major source of uncertainties in overall cloud feedback in GCMs. Thus improving mixed-phase cloud parameterizations in climate models is critical to reducing the climate forecast uncertainties. This study aims at providing improved knowledge of mixed-phase cloud properties from the long-term ACRF observations and improving mixed-phase clouds simulations in the NCAR Community Atmosphere Model version 5 (CAM5). The key accomplishments are: 1) An improved retrieval algorithm was developed to provide liquid droplet concentration for drizzling or mixed-phase stratiform clouds. 2) A new ice concentration retrieval algorithm for stratiform mixed-phase clouds was developed. 3) A strong seasonal aerosol impact on ice generation in Arctic mixed-phase clouds was identified, which is mainly attributed to the high dust occurrence during the spring season. 4) A suite of multi-senor algorithms was applied to long-term ARM observations at the Barrow site to provide a complete dataset (LWC and effective radius profile for liquid phase, and IWC, Dge profiles and ice concentration for ice phase) to characterize Arctic stratiform mixed-phase clouds. This multi-year stratiform mixed-phase cloud dataset provides necessary information to study related processes, evaluate model stratiform mixed-phase cloud simulations, and improve model stratiform mixed-phase cloud parameterization. 5). A new in situ data analysis method was developed to quantify liquid mass partition in convective mixed-phase clouds. For the first time, we reliably compared liquid mass partitions in stratiform and convective mixed-phase clouds. Due to the different dynamics in stratiform and convective mixed-phase clouds, the temperature dependencies of liquid mass partitions are significantly different due to much higher ice concentrations in convective mixed phase clouds. 6) Systematic evaluations of mixed-phase cloud simulations by CAM5 were performed. Measurement results indicate that ice concentrations control stratiform mixed-phase cloud properties. The improvement of ice concentration parameterization in the CAM5 was done in close collaboration with Dr. Xiaohong Liu, PNNL (now at University of Wyoming).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Perspectives on Microfluidics for the Study of Asphaltenes in Upstream Hydrocarbon Production: A Minireview

The utilization of microfluidics has generated deep insights into asphaltene precipitation mechanisms and oil–water emulsion stabilization. Agglomeration and precipitation of asphaltenes can cause flow assurance problems during the extraction and transportation of crude oil. Change in temperature, pressure, reservoir conditions, and solvents can change the local environment, leading to asphaltene precipitation. Understanding asphaltene properties and precipitation pathways becomes critical in devising mitigation methods, demulsifiers, and suitable conditions during hydrocarbon processing. Microfluidics has helped in high throughput measurement studies, understanding critical processing conditions, fast demulsifier screening, and the effect of solvent concentration on deposition, generating useful information for utilization at the point of resource extraction facilitating improved resource management. It has become possible to capture the porous, complex nature of reservoir formations and the interaction of chemicals during precipitation through integrated analytics and visualization studies available only through microfluidics. The use of droplet microfluidics, with optical microscopy and high-speed imaging to study the oil–water interface, has resulted in greater understanding of the role of asphaltenes in interfacial properties and emulsion stabilization. Here, this minireview highlights the crucial aspects of microfluidics that have been used to understand physicochemical behavior and dynamics of asphaltene deposition. Some of the unique devices have been presented focusing on the key elements of microfluidics design, fabrication, and analysis, as the insight obtained from microfluidics strongly depends on the device design and the controllability of the experimental parameters. Successful implementation of microfluidics for efficient and controlled experiments, short analysis time scales and rapid screening, and generation of high-quality, reliable data that convey asphaltene deposition issues and interface behavior in emulsions shows the importance of microsystems for advancing knowledge in hydrocarbon production and processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Modeling Isothermal Reduction of Iron Ore Pellet Using Finite Element Analysis Method: Experiments & Validation

Iron ore pellet reduction experiments were performed with pure hydrogen (H2) and mixtures with carbon monoxide (CO) at different ratios. For direct reduction processes that switch dynamically between reformed natural gas and hydrogen as the reductant, it is important to understand the effects of the transition on the oxide reduction kinetics to optimize the residence time of iron ore pellets in a shaft reactor. Hence, the reduction rates were studied by varying experimental parameters such as the temperature (800, 850 & 900 °C), reactant gas flow rate (100, 150 & 200 cm3/min), pellet size and composition of the reactant gas mixture. The rate of reduction was observed to increase with an increase in temperature and reactant gas flow rate, but it decreased with an increase in pellet size. SEM greyscale analysis was performed to analyze the porosity and phase composition of partially reduced pellets. The porosity of the pellets was observed to increase from 0.3 for unreacted pellet to 0.42 for a completely reduced pellet. Energy-dispersive X-ray spectroscopy (EDAX) analysis was performed to identify the phases observed in the SEM images. The fraction of iron phase was observed to increase from the shell region of the pellet to the core region with an increase in the degree of reduction. A 2D-axisymmetric numerical model was developed on COMSOL Multiphysics, and it was validated using the conversion (X) vs. time curves obtained from each experiment. The model was able to accurately predict the total time needed for the complete conversion of a single iron ore pellet for multiple experiments. Effects of changes in the porosity and tortuosity of the pellet on the model were also studied and the rate of reduction was observed to be sensitive to changes in both porosity and tortuosity. The SEM analysis and the model results show that tortuosity is higher for pellets reduced with H2 than for pellets reduced with H2-CO gas mixtures.

08 HYDROGEN↗
DOE OSTI2019

Effects of Inflow Spatiotemporal Discretization on Wake Meandering and Turbine Structural Response using FAST.Farm

FAST.Farm is a newly developed multiphysics, midfidelity engineering tool that can be used to predict turbine power and structural loads of wind turbines in a wind farm with minimal computational expense. Previous studies have shown similarities and differences in wind turbine performance between FAST.Farm and high-fidelity large-eddy simulations (LES) using both LES precursor-generated inflow and TurbSim-generated synthetic inflow. While conservative resolutions have been used to date, no formal spatial or temporal discretization study has been performed for the FAST.Farm model wind domains. This work aims to study the effects of varying the spatial and temporal discretization of the wind domains on wake meandering (low-resolution domain) and turbine structural response (high-resolution domains) and resulting wake and load calculations. The purpose of this study is to establish convergence criteria and recommendations for discretization values in terms of rotor diameter (D, expressed in meters) that will maximize computational efficiency. To ensure a percent error of ≤ 1% in standard deviation of wake center position relative to the finest included resolution, a low-resolution time step of 0.024D s, 0.016D s, and 0.0079D s and a low-resolution spatial discretization of 0.079D m, 0.16D m, and 0.24D m is recommended for 8 m/s, 12 m/s, and 18 m/s mean ambient wind speed at hub height, respectively. These guidelines are likely applicable to other implementations of the dynamic wake meandering model. To ensure a percent error of ≤ 1% in standard deviation of all considered structural response outputs, a high-resolution time step that captures the highest influential excitation and natural frequencies in the system and a high-resolution spatial discretization comparable to the maximum airfoil cord length is recommended. Lastly, these guidelines are likely applicable to any aeroelastic analysis.

17 WIND ENERGY↗
DOE OSTI2021

Comparative Analysis of ML Techniques for Data-Driven Anomaly Detection, Classification and Localization in Distribution System

High penetration of Distributed Energy Resources (DERs), fundamental load behavior changes, controllable loads, and significant increase in Electrical Vehicles (EVs) lead to complex dynamic behavior of the electric distribution system. Increasing number of components also means more measurements, more data and more data anomalies. Detecting, classifying and localizing these anomalies are important for situational awareness, and at the same time, very challenging given increasing complexity of the system. Highly accurate and high-resolution analytical techniques are needed to support anomaly detection, classification and localization (AD-C-L) for monitoring, root cause analysis and decision making. This paper provides comprehensive review and analysis of the existing spatio-temporal AD-C-L techniques within the distribution system. Challenges for specific problems in AD-C- L have been also discussed in this paper. Existing AD-C- L techniques have been categorized and synthesized for specific merits and limitations of multiple Machine Learning (ML) methodologies using common developed metrics of performance. The comparative analysis is summarized and presented with the open research challenges and path forward for future research needs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Recurrent neural networks for short-term and long-term prediction of geothermal reservoirs

Accurate prediction of geothermal reservoir responses to alternative energy production scenarios is critical for optimizing the development of the underlying resources. While the conventional physics-based models offer a comprehensive prediction tool, data-driven models provide an efficient alternative to build fit-for-purpose predictive models by extracting and using the statistical patterns in the collected data to make predictions. The recurrent neural network (RNN) is a data-driven model that is commonly applied to predict time series sequences. This paper presents a variant of RNN that also utilizes the efficiency of convolutional neural networks (CNN) for the prediction of energy production from geothermal reservoirs. Specifically, a CNN–RNN architecture is developed that takes historical well controls as input (features) and their corresponding production response data as output (labels) to learn an input-output mapping that can predict the future well production responses/performance for any given future well control inputs. The model is paired with a labeling scheme to handle real field disturbances that create data gaps. In addition to the model structure, we introduce a thorough workflow for applying the model, which includes data pre-processing, feature selection, as well as different training strategies for short-term and long-term prediction. Finally, the performance and accuracy of the model are evaluated by applying it to multiple datasets, including a field reservoir model.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Identification and Online Updating of Dynamic Models for Demand Response of an Industrial Air Separation Unit

Demand-response operation of air separation units requires frequent changes in production rate(s), and scheduling calculations must explicitly consider process dynamics to ensure feasibility of the solutions. To this end, scale-bridging models (SBMs) approximate the scheduling-relevant dynamics of a process and its controller in a low-order representation. In contrast to previous works that have employed nonlinear SBMs, this paper proposes linear SBMs, developed using time-series analysis, to facilitate online scheduling computations. Using a year-long industrial dataset, we nd that compact linear SBMs are suitable approximations over typical scheduling horizons, but that their accuracies are unpredictable over time. We introduce a strategy for online updating of the SBMs, based on Kalman ltering schemes for online parameter estimation. The approach greatly improves the accuracy of SBM predictions and will enable the use of linear SBM-based demand-response scheduling in the future.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Web-based wide-area monitoring platform for ringdown and clustering analytics in power systems

This paper introduces an open-source research platform for monitoring the Mexican interconnected power grid, allowing real-time processing and information extraction of the grid’s dynamic condition. Moreover, the platform is a Python-based development that embeds different ringdown and clustering analytics tools. In the case of ringdown analysis, the modal information can be extracted using some of the most known algorithms, i.e., Prony analysis, eigensystem realization algorithm (ERA), and matrix pencil (MP). For clustering analysis, the coherent behaviour of generator and non-generator buses is provided by applying recent state-of-the-art techniques such as affinity propagation, K-means, hierarchical agglomerative clustering, and typicality data analysis. The results of up to 93 PMUs show that this open-source platform suits researchers’ and engineers’ power system dynamic analysis requirements.

Clustering↗
DOE OSTI2025

Cyanobacteria dynamically regulate phycobilisome-to-photosystem excitation energy transfer

In cyanobacteria and red algae, the phycobilisome (PBS) absorbs light and transfers its energy to the chlorophylls in photosystems II (PSII) and I (PSI). With the help of target analysis of time-resolved emission spectra measured at room temperature (RT) and at 77 K, we establish a general kinetic scheme for excitation energy transfer (EET) and trapping based upon a PBS-PSII-PSI megacomplex. At RT it is found that in the dark-adapted cells (State II), the terminal emitter of PBS, allophycocyanin APC680, transfers energy to PSII and PSI with equal rates of ≈50 ns −1 , and that spillover from PSII to PSI is present with rate ≈6 ns −1 . At 77 K, upon transition from State I to State II the EET rate from APC680 to PSII is constant, whereas the rate to PSI increases by 67%. This indicates that a structural change in EET distance in the PBS-PSII-PSI megacomplex underlies the state transition.

Science & Technology - Other Topics↗
DOE OSTI2018

Progress on Complex Langevin simulations of a finite density matrix model for QCD

We study the Stephanov model, which is an RMT model for QCD at finite density, using the Complex Langevin algorithm. Naive implementation of the algorithm shows convergence towards the phase quenched or quenched theory rather than to intended theory with dynamical quarks. A detailed analysis of this issue and a potential resolution of the failure of this algorithm are discussed. We study the effect of gauge cooling on the Dirac eigenvalue distribution and time evolution of the norm for various cooling norms, which were specifically designed to remove the pathologies of the complex Langevin evolution. The cooling is further supplemented with a shifted representation for the random matrices. Unfortunately, none of these modifications generate a substantial improvement on the complex Langevin evolution and the final results still do not agree with the analytical predictions.

Bloch, Jacques↗
DOE OSTI2015

Recurrent dynamics in an epidemic model due to stimulated bifurcation crossovers

Epidemics are known to persist in the form of recurrence cycles. Despite intervention efforts through vaccination and targeted social distancing, peaks of activity for infectious diseases like influenza reappear over time. Analysis of a stochastic model is here undertaken to explore a proposed cycle-generating mechanism – the bifurcation crossover. Time series from simulations of the model exhibit oscillations similar to the temporal signature of influenza activity. Power-spectral density indicates a resonant frequency, which corresponds to the annual seasonality of influenza in temperate zones. The study finds that intervention actions influence the extinguishability of epidemic activity. Asymptotic solution to a backward Kolmogorov equation corresponds to a mean extinction time that is a function of both intervention efficacy and population size. Intervention efficacy must be greater than a certain threshold to increase the chances of extinguishing the epidemic. Agreement of the model with several phenomenological features of epidemic cycles lends to it a tractability that may serve as early warning of imminent outbreaks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Characterization of a new ultra-high pressure shock tube facility for combustion and propulsion studies

A new shock tube facility has been designed, constructed, and characterized at the University of Central Florida. This facility is capable of withstanding pressures of up to 1000 atm, allowing for combustion diagnostics of extreme conditions, such as in rocket combustion chambers or in novel power conversion cycles. For studies with toxic gas impurities, the high initial pressures required the development of a gas delivery system to ensure the longevity of the facility and the safety of the personnel. Data acquisition and experimental propagation were implemented with remote access to ensure safety, paired with a LabVIEW- and Python-based user interface. Thus far, test pressures of 270 atm, blast pressures of 730 atm, and temperatures approaching 10 000 K have been achieved. The extreme limitations of this facility allow for emission spectroscopy to be performed during the oxidation of fuel mixtures, e.g., alkanes diluted in argon and carbon dioxide. Ignition delay times were determined and compared to simulations using chemical kinetic mechanisms. The design, experimental procedures, processes of analysis, and uncertainty determination are outlined, and typical pressure profiles are compared with a new gas dynamics solver and empirical correlations developed across multiple shock tube facilities. Preliminary reactive mixture analyses are included with further investigation of the mixtures outlined.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

First observations from the Kr multi-monochromatic x-ray imager for time and spatially resolved diagnosis of hot implosion cores

This paper presents initial findings from the recently deployed Kr multi-monochromatic x-ray imager (Kr MMI) at the Omega laser facility. The experiment focuses on exploring implosion dynamics in exploding pusher capsules at three distinct initial gas fill densities. Utilizing time-gated and spatially integrated measurements, core size, electron temperature (Te), and electron densities (ne) are extracted through the analysis of the spectral region encompassing the Kr He α and its satellite lines. A comprehensive spectral database, incorporating atomic kinetics, spectroscopic-quality radiation transport, and Stark broadened line shapes, has been developed for rigorous data analysis. These measurements underscore the utility of the new Kr MMI instrument, which, combined with sophisticated analysis techniques, enables the diagnosis of plasma conditions at Te>2000 eV, thereby extending the capabilities beyond the prior Ar MMI design. This is an important stepping stone for achieving time-gated and space-resolved diagnostics of electron temperature, electron density, and heat transport in high temperature implosion cores.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Reversible motions and disordered structure of soft particles in amorphous solids

In amorphous solids, soft vibrational modes derived from normal mode analysis can be utilized to identify the soft particles that are prone to irreversible rearrangements. However, the normal mode analysis approach cannot explain why the spatial distributions of clustered soft particles do not change over time. We define a softness parameter based on the vibrational density of states calculated directly from molecular dynamics simulations with both the harmonic vibrations and anharmonic relaxations being properly captured at finite temperatures. This parameter spontaneously correlates with the real space atomic motions and the dynamics heterogeneity. Using the softness parameter, we show that the softest particles are confined within rigid cages. These particles keep rearranging reversibly without long-range diffusion. The moderately soft particles rearrange irreversibly, and the hard particles mainly participate in vibrations without rearrangement. Here, we also show that the soft particles form locally disordered structures, while the hard particles present strong ordering. These findings confirm the defective nature of soft particles, and provide insights on the nature of softness as the ability to rearrange, but not necessarily irreversibly.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spectral Denoising for Accelerated Analysis of Correlated Ionic Transport

Computation of correlated ionic transport properties from molecular dynamics in the Green-Kubo formalism is expensive, as one cannot rely on the affordable mean square displacement approach. We use spectral decomposition of the short-time ionic displacement covariance to learn a set of diffusion eigenmodes that encode the correlation structure and form a basis for analyzing the ionic trajectories. This allows systematic reduction of the uncertainty and accelerate computations of ionic conductivity in systems with a steady-state correlation structure. We provide mathematical and numerical proofs of the method’s robustness and demonstrate it on realistic electrolyte materials.

74 ATOMIC AND MOLECULAR PHYSICS↗