Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear dynamical systems”

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

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

At least 595 records · Page 33

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗

FARM User Guidance and Instructions

The FARM (Feasible Actuator Range Modifier) module is a component of the RAVEN-based FORCE framework for the analysis of Integrated Energy Systems (IES). FARM aids HERON in the solution of the power dispatch problem by evaluating feasible set-point signals to be issued to the control systems of the different IES unit components. Set-points need to satisfy limits on both production variables (i.e., the variables to be optimized such as the electrical power, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To enforce all these limits, a two-stage approach is adopted. First, the power dispatcher algorithm in HERON module estimates set-points that meet the constraints on the production variables, e.g., power levels and power ramp rate limits. These constraints are called explicit constraints. Then, if necessary, FARM adjusts these set-points to ensure the respect of the limits on the process variables of interest, given the knowledge of the system dynamics acquired through machine learning algorithms. These constraints are called implicit constraints. From this standpoint, FARM constitutes a bridge between the HERON power dispatcher that adopts a simplified description of the IES unit (low-resolution physics) and the HYBRID high-fidelity models (high-resolution physics). In this report, the general workflow and the software requirements of FARM module are summarized, and the detailed instructions for installing FARM software, running built-in example cases, deriving Linear Parameter-Varying (LPV) state-space models, and using FARM for user-defined power dispatch problems are provided.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Relationships between Vehicle Pricing and Features: Data Driven Analysis of the Chinese Vehicle Market

A full-scale understanding of the dynamics of the Chinese vehicle market can benefit stakeholders with respect to rational decision-making and effective long-term investment. This study attempts to discover the common vehicle pricing patterns in the Chinese market by quantifying statistical correlations among critical vehicle features from intrinsic powertrain systems to extrinsic market positioning. The data samples involve almost all passenger vehicle models sold in 2013 to 2019. After comparing multiple statistical methodologies, a log-transformation variant of the multinomial linear regression model was found to be the best one, and the goodness of fit shows that this model can offer stable estimates, which were validated using 2019 market data. The insights achieved are: (1) The price and major performance features of SUVs/crossovers are similar to those of sedans; (2) If all other explicit features remain the same, the price of a Japanese midsize sedan is 62% higher than that of a Chinese midsize sedan, and European midsize vehicles have the highest prices overall. (3) The incremental price of fuel consumption varies by vehicle class and fuel economy. For example, from 30 to 50 MPG, the vehicle price increases by $119 for a Chinese brand sedan vehicle, by $69 for a Chinese brand SUV.

33 ADVANCED PROPULSION SYSTEMS↗

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

Leveraging Groundwater Dynamics to Improve Predictions of Summer Low-Flow Discharges

Summer streamflow predictions are critical for managing water resources; however, warming-induced shifts from snow to rain regimes impact low-flow predictive models. Additionally, reductions in snowpack drive earlier peak flows and lower summer flows across the western United States increasing reliance on groundwater for maintaining summer streamflow. However, it remains poorly understood how groundwater contributions vary interannually. Here, we quantify recession limb groundwater (RLGW), defined as the proportional groundwater contribution to the stream during the period between peak stream flow and low flow, to predict summer low flows across three diverse western US watersheds. We ask (a) how do snow and rain dynamics influence interannual variations of RLGW contributions and summer low flows?; (b) which watershed attributes impact the effectiveness of RLGW as a predictor of summer low flows? Linear models reveal that RLGW is a strong predictor of low flows across all sites and drastically improves low-flow prediction compared to snow metrics at a rain-dominated site. Results suggest that strength of RLGW control on summer low flows may be mediated by subsurface storage. Subsurface storage can be divided into dynamic (i.e., variability saturated) and deep (i.e., permanently saturated) components, and we hypothesize that interannual variability in dynamic storage contribution to streamflow drives RLGW variability. In systems with a higher proportion of dynamic storage, RLGW is a better predictor of summer low flow because the stream is more responsive to dynamic storage contributions compared to deep-storage-dominated systems. Overall, including RLGW improved low-flow prediction across diverse watersheds.

54 ENVIRONMENTAL SCIENCES↗

Topological resonance in graphene-like materials

Abstract In topological materials, interacting with short and strong optical pulses, electrons can accumulate a topological phase during the pulse. Such phase can compensate the dynamic phase resulting in topological resonance, which is visible as a large inter-band transfer of electron population. We study theoretically the topological resonance in materials of the gapped multilayer graphene type. We show that the resonance can be observed only in the systems with finite bandgap. For graphene monolayer the topological resonance can occur only in the field of an elliptically polarized pulse, while for graphene systems with many layers the topological resonance can be also realized in a linearly polarized pulse.

Physics↗

Casimir Light in Dispersive Nanophotonics

Time-varying optical media, whose dielectric properties are actively modulated in time, introduce a host of novel effects in the classical propagation of light, and are of intense current interest. In the quantum domain, time-dependent media can be used to convert vacuum fluctuations (virtual photons) into pairs of real photons. Here, we refer to these processes broadly as “dynamical vacuum effects” (DVEs). Despite interest for their potential applications as sources of quantum light, DVEs are generally very weak, presenting many opportunities for enhancement through modern techniques in nanophotonics, such as using media which support excitations such as plasmon and phonon polaritons. Here, we present a theory of weakly modulated DVEs in arbitrary nanostructured, dispersive, and dissipative systems. A key element of our framework is the simultaneous incorporation of time-modulation and “dispersion” through time-translation-breaking linear response theory. As an example, we use our approach to propose a highly efficient scheme for generating entangled surface polaritons based on time-modulation of the optical phonon frequency of a polar insulator.

42 ENGINEERING↗

Power Analysis of an ePump Applied to the Linear Functions of an Agricultural Planter

Like many other industries, the agricultural industry has recently ex-perienced pressure to reduce vehicle emissions while improving productivity. Electric actuation is perceived as a viable solution to replace or augment hydrau-lic and mechanical actuation. However, electrification presents challenges with regards to linear functions, where hydraulic actuators have advantages in terms of compactness, tolerance to contamination and resistance to shocks. A combined electro-hydraulic actuation architecture can leverage the benefits of both electric and hydraulic actuation, while reducing the drawbacks of both approaches. This work investigates the potential of a centralized electric driven pump (ePump) system powering the pressure-controlled linear functions of an agricul-tural planter, with the goal of improving the operating point of the main supply pump. In this application, the rotary functions are hydraulically actuated, alt-hough such a solution could be applied also to electric rotary actuation. This is accomplished by setting the ePump to boost the pressure supplied by the tractor to the level required by the linear functions. An accumulator is used to stabilize the flow requirements of the linear functions, and control the pressure supplied to the actuators. Two control schemes are proposed for the regulation of the ac-cumulator pressure, one favoring an efficient operating point for the ePump, the other favoring stable steady state operation. A simulation model of the baseline system and proposed system is developed and validated using experimental data from a full-scale machine. Using the vali-dated simulation, both control architectures are then evaluated for improvement in system power consumption and dynamic requirements on the ePump to assess their effectiveness. Both systems demonstrate significant improvement in power consumption over the baseline system, with the best solution improving effi-ciency by 64%.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Dynamic Response of a Semiactive Suspension System with Hysteretic Nonlinear Energy Sink Based on Random Excitation by means of Computer Simulation

This paper aims to investigate the property and behavior of the hysteretic nonlinear energy sink (HNES) coupled to a half vehicle system which is a nine-degree-of-freedom, nonlinear, and semiactive suspension system in order to improve the ride comfort and increase the stability in shock mitigation by using the computer simulation method. The HNES model is a semiactive suspension device, which comprises the famous Bouc–Wen (B-W) model employed to describe the force produced by both the purely hysteretic spring and linear elastic spring of potentially negative stiffness connected in parallel, for the half vehicle system. Nine nonlinear motion equations of the half vehicle system are derived in terms of the seven displacements and the two dimensionless hysteretic variables, which are integrated numerically by employing the direct time integration method for studying both the variables of vertical displacements, velocities, accelerations, chassis pitch angle, and the ride comfort and driver safety, respectively, based on the bump and random road inputs of the pseudoexcitation method as excitation signal. Simulation results show that, compared with the HNES model and the magnetorheological (MR) model coupled to the half vehicle system, the ride comfort and stability have been evidently improved. A successful validation process has been performed, which indicated that both the ride comfort and driver safety properties of the HNES model coupled to half vehicle significantly improved.

Chen, Hui↗

Spin wave Hamiltonian and anomalous scattering in NiPS 3

Here, we report a comprehensive spin wave analysis of the semiconducting honeycomb van der Waal antiferromagnet NiPS 3 . Using single-crystal inelastic neutron scattering, we map out the full Brillouin zone and fit the observed modes to a spin wave model with rigorously defined uncertainty. We find that the third-nearest-neighbor exchange J 3 dominates the Hamiltonian, a feature which we fully account for by ab initio density functional theory calculations. We also quantify the degree to which the threefold rotation symmetry is broken and account for the Q = 0 excitations observed in other measurements, yielding a spin exchange model which is consistent across multiple experimental probes. We also identify a strongly reduced static ordered moment and reduced low-energy intensity relative to the linear spin wave calculations, signaling unexplained features in the magnetism which requires going beyond the linear spin wave approximation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Polymer Macrocycles: A novel topology to control dynamics of rubbery materials (Final Report)

One of the current challenges in polymer physics and engineering is to understand the behavior of circular macromolecules. Though such molecules have been made by dilute solution ring closure methods in the past, such procedures have been limited to relatively small, unentangled rings as well as to small sample quantities. Initial studies were also hampered by issues surrounding purity of the rings, though this problem has been somewhat overcome by the use of liquid chromatography at critical conditions (LCCC) methods that separate rings and linear chains. Furthermore, in spite of the continued progress a full understanding of ring dynamics remains elusive. The present work has addressed multiple aspects of the behavior of ring molecules by using a novel synthetic route: reversible radical recombination redox polymerization (R3P) to produce poly(3,6-dioxa-1,8-octanedithiol) (PolyDODT) molecules. The synthesis conditions were chosen to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT) in order to investigate dilute solution and rheological responses of the materials. Importantly, the R3P synthesis makes larger quantities and higher molecular weights than do the dilute solution ring closure methods and also has the potential to be scaled to industrial quantities, thus becoming commercializable. We emphasize two important outcomes of the work. The first is that the rheological investigation shows that very large molecular weight the PolyDODT rings permitted the investigation of ring dynamics from glassy towards terminal flow for molecular sizes corresponding to approximately 300 entanglements of the linear counterpart. At the same time, these very large sized rings could be diluted in a solvent in order to reduce the molecular entanglement density to less than unity. The results demonstrate that ring viscosity follows a Rouse-like linear dependence on molecular weight or entanglement number to approximately 15 entanglements. This range of "Rouse-like" response is greater than what was achieved in all prior works where synthesis conditions have limited the molecular sizes achievable. After the entanglement number $Z_w$ increases beyond approximately 15, we find that the viscosity-molecular weight scaling undergoes a rapid change from the Rouse-like $η\sim{Z}^1_w$ to a very strong power-law scaling that is approximately $η\sim{Z}^{5.8}_w$. This aspect of the work strongly suggests that prior investigations on highly pure rings have invariably been performed on unentangled systems simply because ring molecules have very high entanglement thresholds. The results are further supported by plateau modulus measurements that show that in the "Rouse-like" viscosity regime, for the RDODT samples there is a corresponding lack of a rubbery entanglement plateau and once entanglement coupling is seen, the rubbery plateau also appears in the dynamic modulus data. Furthermore, upon appropriate scaling, we show that the unentangled polyDODT rings (i.e., the diluted high molecular weight rings) show virtually identical dynamics (upon appropriate scaling to account for the chain dilutions) to the dynamics of similarly unentangled LCCC fractionated polystyrene rings.

36 MATERIALS SCIENCE↗

Extracted Water Induces Concentration Fluctuations in Model Ternary Liquid–Liquid Extraction System

Organic phase aggregation and phase splitting are two important (and connected) phenomena in liquid–liquid extraction of metal ions. Previously, we demonstrated how structure in binary mixtures of extractant and diluent originates from concentration fluctuations associated with the so-called “third-phase formation” liquid–liquid phase instability. For more complex systems, we sought to understand how extracted aqueous components affect these concentration fluctuations. Here, in this study, we use small-angle X-ray scattering and molecular dynamics simulations to investigate organic phase structural changes upon contact with water in a model ternary system of water/alkane/extractant. Critical exponents were determined from scattering with temperature variation and were consistent with mean field values. Then, by controlling the water content in the organic phase through direct dissolution of fixed quantities at constant temperature, we relate the observed increase in the fluctuation length scale to changes to the spinodal temperature of the third-phase formation phase boundary. We find a roughly linear relationship between water concentration and the spinodal temperature and a similar slope in that relationship for different linear alkane diluents. This suggests that the mechanism by which water extraction impacts the organic phase structure is to enhance concentration fluctuations through the introduction of additional polar extractant-water complexes, which increases the spinodal temperature and thereby decreases the reduced temperature, i.e., the distance to the instability. This explanation for how solutes induce aggregation in liquid–liquid extraction organic phases may extend to more complex systems, explaining the large increases in correlation lengths reported under high acid and metal loading.

Bonnett, Brittany L. [Argonne National Laboratory ↗

Molecular diversity of dissolved organic matter reflects macroecological patterns in river networks

Deciphering dissolved organic matter (DOM) molecular complexity is crucial for understanding ecosystem function. Using the continental-scale Worldwide Hydrobiogeochemistry Observation Network for Dynamic Rivers Systems (WHONDRS) Fourier-transform ion cyclotron resonance mass spectrometry (FTICR-MS) dataset, we reveal fundamental scaling patterns of DOM chemodiversity with watershed characteristics. Analysis of 54 river sites shows local and regional watershed features significantly influence DOM chemodiversity (2500–8718 unique formulae), exhibiting consistent scaling patterns across compound classes and a novel latitudinal gradient (decreasing diversity with increasing latitude). Scaling relationships for DOM composition vary by compound class. Crucially, the scaling parameters (B, baseline chemodiversity; Z, sensitivity) are linearly interrelated. This B–Z relationship is most robust for potentially bio-labile carbohydrates (coefficient of determination R 2 ≈ 0.85), diminishing for recalcitrant, plant-derived molecules (such as lignin), and indicates (potential) biolability-dependent coupling between baseline diversity and environmental responsiveness. These quantitative scaling relationships, with scaling exponents ranging from − 2.1 to 2.2 across compound classes, enable prediction of DOM composition across watersheds, offering a framework to understand ecosystem responses to environmental change. This research bridges biogeochemistry and ecology, providing tools to anticipate molecular transformations across scales.

59 BASIC BIOLOGICAL SCIENCES↗

Evaluation of silicon and indium gallium arsenide photodiodes as direct timing detectors for pulsed x-ray systems

Benchtop pulsed x-ray systems are commonly used to record dynamic material data on the order of nanoseconds, but pulse timing is often difficult to accurately determine. This study demonstrates that commercially available photodiodes can be used effectively for direct x-ray pulse detection without the need for visible light scintillators. X-ray pulses from four commercially available flash x-ray systems were quantified using one silicon and two indium gallium arsenide (InGaAs) photodiodes. The measured InGaAs pulse durations were strongly dependent on radiation dose in the non-linear operating regime, so the photodiodes were shielded to operate below the 2.5 V non-linear regime threshold. The average pulse duration and pulse arrival time jitter of the photodiodes for each x-ray source were within several nanoseconds with the exception of two sets of measurements that were affected by low instrument sensitivity and electrical noise. These results show that InGaAs photodiodes can be used as effective and repeatable stand-alone timing diagnostics for x-ray pulses as short as 20 ns or less.

47 OTHER INSTRUMENTATION↗

Stability and Performance of an Undersea Kite Operating in a Turbulent Flow Field

In this paper, we examine the effects of flow disturbances resulting from turbulence on the dynamic behavior of an underwater energy-harvesting kite system that executes periodic figure-8 flight. Due to the periodic nature of the kite’s operation, we begin by assessing orbital stability using the Floquet analysis and stroboscopic intersection analysis of a Poincaré section, with the former analysis performed on a simplified “unifoil” model and the latter performed on a six-degree-of-freedom (6-DOF)/flexible tether model. With periodic stability established, a frequency-domain analysis based on a linearization about the kite’s path is used to predict the quality of flight path tracking as a function of the turbulence frequency. To validate the accuracy of these simulation-based predictions under flow disturbances, we compare the predictions of the kite’s behavior against the results of small-scale tow testing experiments performed in a controlled pool environment.

13 HYDRO ENERGY↗

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗

Link statistics of dislocation network during strain hardening

Dislocations are line defects in crystals that multiply and self-organize into a complex network during strain hardening. The length of dislocation links, connecting neighboring nodes within this network, contains crucial information about the evolving dislocation microstructure. By analyzing data from Discrete Dislocation Dynamics (DDD) simulations in face-centered cubic (fcc) Cu, we characterize the statistical distribution of link lengths of dislocation networks during strain hardening on individual slip systems. Here, our analysis reveals that link lengths on active slip systems follow a double-exponential distribution, while those on inactive slip systems conform to a single-exponential distribution. The distinctive long tail observed in the double-exponential distribution is attributed to the stress-induced bowing out of long links on active slip systems, a feature that disappears upon removal of the applied stress. We further demonstrate that both observed link length distributions can be explained by extending a one-dimensional Poisson process to include different growth functions. Specifically, the double-exponential distribution emerges when the growth rate for links exceeding a critical length becomes super-linear, which aligns with the physical phenomenon of long links bowing out under stress. This work advances our understanding of dislocation microstructure evolution during strain hardening and elucidates the underlying physical mechanisms governing its formation.

Crystal plasticity↗