SEARCH · Engineering Papers

Results for “Complex 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.

SourceResearch connectionFrom yearThrough year

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

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Shot noise detection in hBN-based tunnel junctions

High quality Au/hBN/Au tunnel devices are fabricated using transferred atomically thin hexagonal boron nitride as the tunneling barrier. All tunnel junctions show tunneling resistance on the order of several kX/lm2. Ohmic I-V curves at small bias with no signs of resonances indicate the sparsity of defects. Tunneling current shot noise is measured in these devices, and the excess shot noise shows consistency with theoretical expectations. Furthermore, these results show that atomically thin hBN is an excellent tunnel barrier, especially for the study of shot noise properties, and this can enable the study of the tunneling density of states and shot noise spectroscopy in more complex systems.

42 ENGINEERING↗

DOE OSTI2018

Multi-fidelity uncertainty quantification strategies for large-scale multiphysics applications: PSAAP II particle-based solar energy receiver

The study of complex systems is often based on computationally intensive, high-fidelity, simulations. To build confidence or improve the prediction accuracy of such simulations, the impact of uncertainties in model inputs, or even the structure of the models, on the quantities of interest must be measured. This, however, requires a computational budget that is a possibly large multiple of the cost of a single simulation, and thus may become infeasible for expensive simulation models featuring a large number of uncertain inputs and highly non-linear behavior. Therefore, this work explores multi-fidelity strategies to accelerate the estimation of the effect of uncertainties. The main idea behind multi-fidelity methods is to utilize cheaper, lower-fidelity models – than the intended high-fidelity, expensive model of the problem – to generate a baseline solution that together with relatively small number of high-fidelity simulations can lead to accurate predictions. The methods are briefly presented, and their performance assessed on an irradiated particle-laden turbulent flow case related to Stanford’s PSAAP II particle-based solar energy receiver.

Physics↗

DOE OSTI2019

Ultrathin transmissive metasurfaces for multi-wavelength optics in the visible

Metasurface-based optical elements have the ability to shape wavefronts by locally changing the properties of incident illumination. They hold great potential to promote a new generation of wearable devices and thin optical systems for imaging and sensing. However, due to the diffractive nature of these metasurfaces, chromatic aberration represents a critical challenge toward practical implementations, especially for the visible spectrum. Existing solutions for multi-wavelength meta-optics inevitably increase the device thickness and system complexity, similar to the conventional refractive and diffractive counterparts. They are also limited to polarization sensitive lenses with a small size and numerical aperture (NA). In this manuscript, we present the design of single-layer multi-wavelength metasurfaces using ultrathin dielectric resonators (thickness << wavelength) as structural elements. The electromagnetic response of these resonators can be tailored by the in-plane geometrical parameters and used to manipulate the transmitted light. We introduce a digital transmission modulation scheme which enables opportunities to design multi-wavelength optics beyond conventional binary diffractive optics. The design versatility of our approach is demonstrated by polarization-independent achromatic metalenses with a larger lens size and larger NA than existing approaches. Furthermore, the simplicity of the design and fabrication process makes these ultrathin metasurfaces good candidates for building flat optical elements that can be easily integrated onto CMOS electronics and MEMS devices.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Spiral waves in externally excited neuronal network: Solvable model with a monotonically differentiable magnetic flux

Spiral waves are particular spatiotemporal patterns connected to specific phase singularities representing topological wave dislocations or nodes of zero amplitude, witnessed in a wide range of complex systems such as neuronal networks. The appearance of these waves is linked to the network structure as well as the diffusion dynamics of its blocks. We report a novel form of the Hindmarsh-Rose neuron model utilized as a square neuronal network, showing the remarkable multistructure of dynamical patterns ranging from characteristic spiral wave domains of spatiotemporal phase coherence to regions of hyperchaos. The proposed model comprises a hyperbolic memductance function as the monotone differentiable magnetic flux. Hindmarsh-Rose neurons with an external electromagnetic excitation are considered in three different cases: no excitation, periodic excitation, and quasiperiodic excitation. We performed an extensive study of the neuronal dynamics including calculation of equilibrium points, bifurcation analysis, and Lyapunov spectrum. We have found the property of antimonotonicity in bifurcation scenarios with no excitation or periodic excitation and identified wide regions of hyperchaos in the case of quasiperiodic excitation. Furthermore, the formation and elimination of the spiral waves in each case of external excitation with respect to stimuli parameters are investigated. We have identified novel forms of Hindmarsh-Rose bursting dynamics. Our findings reveal multipartite spiral wave formations and symmetry breaking spatiotemporal dynamics of the neuronal model that may find broad practical applications. © 2019 Author(s).

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

Memory embedded non-intrusive reduced order modeling of non-ergodic flows

Generating a digital twin of any complex system requires modeling and computational approaches that are efficient, accurate, and modular. Traditional reduced order modeling techniques are targeted at only the first two, but the novel nonintrusive approach we present in this study is an attempt at taking all three into account effectively compared to their traditional counterparts. Based on dimensionality reduction using proper orthogonal decomposition (POD), we introduce a long short-term memory neural network architecture together with a principal interval decomposition (PID) framework as an enabler to account for localized modal deformation. As an effective partitioning tool for breaking the Kolmogorov barrier, our PID framework, therefore, can be considered a key element in the accurate reduced order modeling of convective flows. Our applications for convection-dominated systems governed by Burgers, Navier-Stokes, and Boussinesq equations demonstrate that the proposed approach yields significantly more accurate predictions than the POD-Galerkin method and could be a key enabler toward near real-time predictions of unsteady flows.

42 ENGINEERING↗

DOE OSTI2025

Filling data analysis gaps in time-resolved crystallography by machine learning

There is a growing understanding of the structural dynamics of biological molecules fueled by x-ray crystallography experiments. Time-resolved serial femtosecond crystallography (TR-SFX) with x-ray Free Electron Lasers allows the measurement of ultrafast structural changes in proteins. Nevertheless, this technique comes with some limitations. One major challenge is the quality of data from TR-SFX measurements, which often faces issues like data sparsity, partial recording of Bragg reflections, timing errors, and pixel noise. To overcome these difficulties, conventionally, large volumes of data are collected and grouped into a few temporal bins. The data in each bin are then averaged and paired with the mean of their corresponding jittered timestamps. This procedure provides one structure per bin, resulting in a limited number of averaged structures for the entire time interval spanned by the experiment. Therefore, the information on ultrafast structural dynamics at high temporal resolution is lost. This has initiated research for advanced methods of analyzing experimental TR-SFX data beyond the standard binning and averaging method. To address this problem, we use a machine learning algorithm called Nonlinear Laplacian Spectral Analysis (NLSA), which has emerged as a promising technique for studying the dynamics of complex systems. In this work, we demonstrate the power of this algorithm using synthetic x-ray diffraction snapshots from a protein with significant data incompleteness, timing uncertainties, and noise. Our study confirms that NLSA is a suitable approach that effectively mitigates the effects of these artifacts in TR-SFX data and recovers accurate structural dynamics information hidden in such data.

Trujillo, Justin (ORCID:0000000285505360)↗

DOE OSTI2021

Time-resolved photoelectron imaging of complex resonances in molecular nitrogen

We have used the FERMI free-electron laser to perform time-resolved photoelectron imaging experiments on a complex group of resonances near 15.38 eV in the absorption spectrum of molecular nitrogen, N-2, under jet-cooled conditions. The new data complement and extend the earlier work of Fushitani et al. [Opt. Express 27, 19702-19711 (2019)], who recorded time-resolved photoelectron spectra for this same group of resonances. Time-dependent oscillations are observed in both the photoelectron yields and the photoelectron angular distributions, providing insight into the interactions among the resonant intermediate states. In addition, for most states, we observe an exponential decay of the photoelectron yield that depends on the ionic final state. This observation can be rationalized by the different lifetimes for the intermediate states contributing to a particular ionization channel. Finally, although there are nine resonances within the group, we show that by detecting individual photoelectron final states and their angular dependence, we can identify and differentiate quantum pathways within this complex system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

AENET–LAMMPS and AENET–TINKER : Interfaces for accurate and efficient molecular dynamics simulations with machine learning potentials

Machine-learning potentials (MLPs) trained on data from quantum-mechanics based first-principles methods can approach the accuracy of the reference method at a fraction of the computational cost. To facilitate efficient MLP-based molecular dynamics and Monte Carlo simulations, an integration of the MLPs with sampling software is needed. Here, we develop two interfaces that link the atomic energy network (ænet) MLP package with the popular sampling packages TINKER and LAMMPS. The three packages, ænet, TINKER, and LAMMPS, are free and open-source software that enable, in combination, accurate simulations of large and complex systems with low computational cost that scales linearly with the number of atoms. Scaling tests show that the parallel efficiency of the ænet–TINKER interface is nearly optimal but is limited to shared-memory systems. The ænet–LAMMPS interface achieves excellent parallel efficiency on highly parallel distributed memory systems and benefits from the highly optimized neighbor list implemented in LAMMPS. We demonstrate the utility of the two MLP interfaces for two relevant example applications: the investigation of diffusion phenomena in liquid water and the equilibration of nanostructured amorphous battery materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Inexact iterative numerical linear algebra for neural network-based spectral estimation and rare-event prediction

Understanding dynamics in complex systems is challenging because there are many degrees of freedom, and those that are most important for describing events of interest are often not obvious. The leading eigenfunctions of the transition operator are useful for visualization, and they can provide an efficient basis for computing statistics, such as the likelihood and average time of events (predictions). Here, we develop inexact iterative linear algebra methods for computing these eigenfunctions (spectral estimation) and making predictions from a dataset of short trajectories sampled at finite intervals. We demonstrate the methods on a low-dimensional model that facilitates visualization and a high-dimensional model of a biomolecular system. Implications for the prediction problem in reinforcement learning are discussed.

Chemistry↗

DOE OSTI2024

Estimating the carbon footprint of computational fluid dynamics

In this study, computational resources have grown exponentially in the past few decades. These machines make possible research and design in fields as diverse as medicine, astronomy, and engineering. Despite ever-increasing computational capabilities, direct simulation of complex systems has remained challenging owing to the degrees of freedom involved. At the cusp of exascale computing, high-resolution simulation of practical problems with minimal model assumptions may soon experience a renaissance. However, growing reliance on modern computers comes at the cost of a growing carbon footprint. To illustrate this, we examine historic computations in fluid dynamics where larger computers have afforded the opportunity to simulate flows at increasingly relevant Reynolds numbers. Under a variety of flow configurations, the carbon footprint of such simulations is found to scale roughly with the fourth power of Reynolds number. This is primarily explained by the computation cost in core-hours, which is also described by similar scaling, though regional differences in renewable energy use also play a role. Using the established correlation, we examine a large database of simulations to develop estimates for the carbon footprint of computational fluid dynamics in a given year. Collectively, the analysis provides an additional benchmark for new computations where, in addition to balancing considerations of model fidelity, carbon footprint should also be considered.

42 ENGINEERING↗

DOE OSTI2024

Out of time order correlation of the Hubbard model with random local disorder

The out-of-time-order correlator (OTOC) serves as a powerful tool for investigating quantum information spreading and chaos in complex systems. We present a method employing non-equilibrium dynamical mean-field theory and coherent potential approximation combined with diagrammatic perturbation on the Schwinger–Keldysh contour to calculate the OTOC for correlated fermionic systems subjected to both random disorder and electron interaction. Furthermore, our key finding is that random disorder enhances the OTOC decay in the Hubbard model for the metallic phase in the weakly interacting limit. However, the current limitation of our perturbative solver restricts the applicability to weak interaction regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Extended dynamic mode decomposition for model reduction in fluid dynamics simulations

High computational cost and storage/memory requirements of fluid dynamics simulations constrain their usefulness as a predictive tool. Reduced-order models (ROMs) provide a viable solution to this challenge by extracting the key underlying dynamics of a complex system directly from data. We investigate the efficacy and robustness of an extended dynamic mode decomposition (xDMD) algorithm in constructing ROMs of three-dimensional cardiovascular computations. Focusing on the ROMs' accuracy in representation and interpolation, we relate these metrics to the truncation rank of singular value decomposition, which underpins xDMD and other approaches to ROM construction. Our key innovation is to relate the truncation rank to the singular values of the original flow problem. This result establishes a priori guidelines for the xDMD deployment and its likely success as a means of data compression and reconstruction of the system's dynamics from dominant spatiotemporal structures present in the data.

Mechanics↗

DOE OSTI2025

Reproducibility of fixed-node diffusion Monte Carlo across diverse community codes: The case of water–methane dimer

Fixed-node diffusion quantum Monte Carlo (FN-DMC) is a widely trusted many-body method for solving the Schrödinger equation, known for its reliable predictions of material and molecular properties. Furthermore, its excellent scalability with system complexity and near-perfect utilization of computational power make FN-DMC ideally positioned to leverage new advances in computing to address increasingly complex scientific problems. Even though the method is widely used as a computational gold standard, reproducibility across the numerous FN-DMC code implementations has yet to be demonstrated. This difficulty stems from the diverse array of DMC algorithms and trial wave functions, compounded by the method’s inherent stochastic nature. Here, this study represents a community-wide effort to assess the reproducibility of the method, affirming that yes, FN-DMC is reproducible (when handled with care). Using the water–methane dimer as the canonical test case, we compare results from eleven different FN-DMC codes and show that the approximations to treat the non-locality of pseudopotentials are the primary source of the discrepancies between them. In particular, we demonstrate that, for the same choice of determinantal component in the trial wave function, reliable and reproducible predictions can be achieved by employing the T-move, the determinant locality approximation, or the determinant T-move schemes, while the older locality approximation leads to considerable variability in results. These findings demonstrate that, with appropriate choices of algorithmic details, fixed-node DMC is reproducible across diverse community codes—highlighting the maturity and robustness of the method as a tool for open and reliable computational science.

Della Pia, Flaviano [Univ. of Cambridge (United Ki↗

DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗

DOE OSTI2025

Demonstration of a 1820 channel multiplexer for transition-edge sensor bolometers

The scalability of most transition-edge sensor arrays is limited by the multiplexing technology, which combines their signals over a reduced number of wires and amplifiers. Here, in this Letter, we present and demonstrate a multiplexer design optimized for transition-edge sensor bolometers with 1820 sensors per readout unit, a factor of two more than the previous state-of-the-art. The design is optimized for cosmic microwave background imaging applications, and it builds on previous microwave superconducting quantum interference device multiplexers by doubling the available readout bandwidth to the full 4–8 GHz octave. Evaluating the key performance metrics of yield, sensitivity, and crosstalk through laboratory testing, we find an end-to-end operable detector yield of 78%, a typical nearest-neighbor crosstalk amplitude of ∼0.4%, and a median white noise level of 83 pA/$\sqrt{\textrm{Hz}}$ due to the multiplexer, corresponding to an estimated contribution of 4% to the total system noise for a ground-based cosmic microwave background telescope. Additionally, we identify a possible path toward reducing resonator loss for future designs with reduced noise. We expect these developments to alleviate the system complexity, cryogenic requirements, and cost of future large arrays of low temperature detectors.

cosmic microwave background↗

DOE OSTI2015

Transversely diode-pumped alkali metal vapour laser

We have studied theoretically the operation of a transversely diode-pumped alkali metal vapour laser. For the case of high-intensity laser radiation, we have obtained an analytical solution to a complex system of differential equations describing the laser. This solution allows one to exhaustively determine all the energy characteristics of the laser and to find optimal parameters of the working medium and pump radiation (temperature, buffer gas pressure, and intensity and width of the pump spectrum). (lasers)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Data-driven parameterization of the generalized Langevin equation

Significance The generalized Langevin equation (GLE) provides a precise description of coarse-grained variable dynamics in reduced dimension models. However, computation of the memory kernel poses a major challenge to the practical use of the GLE. This paper presents a data-driven approach to compute the memory kernel, relying on a hierarchy of parameterized rational approximations in terms of the Laplace transform, which can be expanded to arbitrarily high order as necessary. This parameterization makes it convenient to represent the GLE via an extended stochastic model where the memory term is eliminated by properly introducing auxiliary variables. The present method is well-suited for constructing reduced models for nonequilibrium properties of complex systems such as biomolecules, chemical reaction networks, and climate simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Hidden domain boundary dynamics toward crystalline perfection

A central paradigm of nonequilibrium physics concerns the dynamics of heterogeneity and disorder, impacting processes ranging from the behavior of glasses to the emergent functionality of active matter. Understanding these complex mesoscopic systems requires probing the microscopic trajectories associated with irreversible processes, the role of fluctuations and entropy growth, and the timescales on which nonequilibrium responses are ultimately maintained. Approaches that illuminate these processes in model systems may enable a more general understanding of other heterogeneous nonequilibrium phenomena, and potentially define ultimate speed and energy cost limits for information processing technologies. Here, we apply ultrafast single-shot X-ray photon correlation spectroscopy to resolve the nonequilibrium, heterogeneous, and irreversible mesoscale dynamics during a light-induced phase transition in a (PbTiO 3 ) 16 /(SrTiO 3 ) 16 superlattice. Such ferroelectric superlattice systems are a useful platform to study phase transitions and topological dynamics due to their high degree of tunability. This provides an approach for capturing the nucleation of the light-induced phase, the formation of transient mesoscale defects at the boundaries of the nuclei, and the eventual annihilation of these defects, even in systems with complex polarization topologies. We identify a nonequilibrium correlation response spanning >10 orders of magnitude in timescales, with multistep behavior similar to the plateaus observed in supercooled liquids and glasses. We further show how the observed time-dependent long-time correlations can be understood in terms of stochastic and non-Markovian dynamics of domain walls, encoded in waiting-time distributions with power-law tails. This work defines possibilities for probing the nonequilibrium and correlated dynamics of disordered and heterogeneous media.

36 MATERIALS SCIENCE↗