SEARCH · Search NASA

Results for “fluctuation and dynamics”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Molecular basis for higher affinity of SARS-CoV-2 spike RBD for human ACE2 receptor

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused substantially more infections, deaths, and economic disruptions than the 2002-2003 SARSCoV. The key to understanding SARS-CoV-2's higher infectivity lies partly in its host receptor recognition mechanism. Experiments show that the human angiotensin converting enzyme 2 (ACE2) protein, which serves as the primary receptor for both CoVs, binds to the receptor binding domain (RBD) of CoV-2's spike protein stronger than SARS-CoV's spike RBD. The molecular basis for this difference in binding affinity, however, remains unexplained from X-ray structures. To go beyond insights gained from X-ray structures and investigate the role of thermal fluctuations in structure, we employ all-atom molecular dynamics simulations. Microseconds-long simulations reveal that while CoV and CoV-2 spike-ACE2 interfaces have similar conformational binding modes, CoV-2 spike interacts with ACE2 via a larger combinatorics of polar contacts, and on average, makes 45% more polar contacts. Correlation analysis and thermodynamic calculations indicate that these differences in the density and dynamics of polar contacts arise from differences in spatial arrangements of interfacial residues, and dynamical coupling between interfacial and non-interfacial residues. Furthermore, these results recommend that ongoing efforts to design spike-ACE2 peptide blockers will benefit from incorporating dynamical information as well as allosteric coupling effects.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Influence of Rigidity–Hydration Coupling on Size-Dependent Diffusion in Hydrated Polymer Membranes

Selective ion transport in polymer membranes depends critically on how penetrant motion couples to polymer dynamics and hydration. Yet, the mechanistic interplay between polymer rigidity, water content, and penetrant size remains poorly understood, especially in the regime where the penetrant diameter, polymer Kuhn length, and correlation length are comparable. Here, we employ coarse-grained molecular dynamics simulations to systematically investigate penetrant diffusion in hydrated polymer networks across a broad range of water volume fractions, chain rigidities, and penetrant sizes. The results reveal a transition from a decoupled regime, where small penetrants diffuse nearly independently of polymer relaxation, to a coupled regime in which large penetrants require cooperative polymer motion for transport. Increasing polymer rigidity amplifies the sensitivity of diffusivity to hydration, particularly at low water content, leading to pronounced deviations from Stokes−Einstein scaling. Comparison with scaling theories and free-volume models shows that classical nanoparticle-based frameworks fail to capture this intermediate regime. To address this gap, we extend the Yasuda model to incorporate polymer rigidity through a single parameter that quantifies the dynamic contribution of chain stiffness to free-volume fluctuations. The resulting model collapses diffusivity data across all sizes, water contents, and rigidities, providing a unified description of penetrant transport in hydrated polymer matrices. Furthermore, these findings establish polymer rigidity as a key, tunable determinant of diffusion and offer a framework for interpreting size-dependent transport in ion-selective membranes.

diffusion↗

DOE OSTI2026

Spectral Properties and Energy Injection in Mercury's Magnetotail Current Sheet

Mercury's magnetotail hosts a thin and highly dynamic current sheet (CS), where magnetic reconnection and strong fluctuations frequently occur. Here, we statistically analyze magnetic field power spectra across 370 magnetotail CSs observed by MESSENGER. About 20% of the events are quasi-laminar, showing single power-law spectra, whereas ∼80% are turbulent, exhibiting a spectral break separating inertial and kinetic ranges. A dawn–dusk asymmetry is identified: inertial-range slopes are systematically shallower on the dawnside, whereas kinetic-range slopes are steeper, indicating more developed turbulence there, consistent with the higher occurrence of reconnection-related processes on the dawnside. Component analysis shows that the transverse components, orthogonal to the tail-aligned principal field ( B X ), display shallow slopes near −1 in the inertial range, suggesting energy injection at ion scales rather than a classical inertial range. These results demonstrate that Mercury's unique plasma environment fundamentally reshapes the initiation of turbulence and the redistribution of energy in the magnetotail.

Mercury↗

DOE OSTI2024

ELM-free enhanced D α H-mode with near zero NBI torque injection in DIII-D tokamak

Enhanced D α H-mode (EDA H-mode), an ELM-free H-mode regime, and the concomitant quasi-coherent mode (QCM) are explored in neutral beam heated, lower single null plasmas with near zero external torque injection. This regime exhibits good energy confinement (H 98y2 ~ 1) with β N ~ 2, high density, regime access at low input power, and no ELMs. New time–resolved measurements of electron and ion density, temperature, plasma rotation, and radial electric field as well as turbulence, flows, thermal and particle transport, and linear instability calculations examining the dynamics of the QCM are presented. Measurements using multiple fluctuation diagnostics reveal the QCM to be a separatrix spanning mode, peaking just inside the separatrix, existing in a wide range of k ⊥ ρ s ~ 0.1–1.2 with multiple harmonics, and propagating with a very small phase velocity in the plasma frame, where k ⊥ is the binormal wavenumber and ρ s is the ion sound radius. Linear gyrokinetic simulations of this EDA H-mode discharge with CGYRO indicate that the trapped electron mode (TEM) and electron temperature gradient mode are the dominant instabilities in the radial region where QCM is also unstable. Qualitative analysis indicates that the properties of TEM are consistent with the experimentally observed characteristics of the QCM. These similarities suggest that the QCM is consistent with a TEM instability existing in the edge region of the EDA H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Characterizing non-Markovian and coherent errors in quantum simulation

Quantum simulation of many-body systems, particularly using ultracold atoms and trapped ions, presents a unique form of quantum control—it is a direct implementation of a multi-qubit gate generated by the Hamiltonian. As a consequence, it also faces a unique challenge in terms of benchmarking, because the well-established gate benchmarking techniques are unsuitable for this form of quantum control. Here we show that the symmetries of the target many-body Hamiltonian can be used not only to benchmark but to characterize experimental errors in the quantum simulation. We use our results to develop protocols to characterize these errors, which can be implemented using state-of-the-art technology. We consider two forms of errors: (i) unitary errors arising out of systematic errors in the applied Hamiltonian and (ii) canonical non-Markovian errors arising out of random shot-to-shot fluctuations in the applied Hamiltonian. We show that the dynamics of the expectation value of the target Hamiltonian itself, which is ideally constant in time, can be used to characterize these errors. In the presence of errors, the expectation value of the target Hamiltonian shows a characteristic thermalization dynamics, when it satisfies the operator thermalization hypothesis (OTH). That is, an oscillation in the short time followed by relaxation to a steady-state value in the long time limit. We show that while the steady-state value can be used to characterize the coherent errors, the amplitude of the oscillations can be used to estimate the non-Markovian errors. We prove a sandwich theorem to establish a linear relation between the amplitude of the oscillations and the magnitude of the non-Markovian errors. Moreover, by varying the initial state, we show that the steady state values can be used to completely construct the generator of the coherent errors. Using these results, we develop two experimental protocols to characterize the unitary errors based on these results, one of which requires single-qubit addressing and the other one doesn't. We also develop a protocol to characterize non-Markovian errors. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Diffraction-quality, ultraflexible protein single crystals engineered with DNA

DNA-functionalized colloidal nanoparticles assemble through flexible, nanoscale DNA hybridization interactions that limit atomic-level structural order. Here, we report a valence-centric strategy that enables DNA-bonded, protein single crystals with unconventional mechanical properties. An octameric enzyme, glutarate L-2-hydroxylase, was site- and number-selectively conjugated with eight self-complementary single-stranded DNA, yielding octavalent molecular bonds. The resulting conjugate assembled into the designed body-centered tetragonal crystals that diffracted to 1.42- to 2.61-angstrom resolution, with contacts mediated by B-form DNA helices spanning 17 to 25 angstroms. Increasing oligonucleotide length induces anisotropic lattice expansion while preserving atomic periodicity, even with partial DNA occupancy. Mechanistic studies suggest that the dynamic motion of unhybridized DNA facilitates crystallization, analogous to fluctuating electron clouds in atomic bonding. Compared with native protein crystals, DNA-hybridized crystals are 23-fold softer. These results challenge the assumption that flexibility is incompatible with structural order and establish a programmable framework for biomolecular crystallization and nanomaterials engineering with atomic precision.

Han, Zhenyu [Department of Chemistry, Northwestern↗

DOE OSTI2025

Quantum entanglement in nuclear fission

Nuclear fission presents a unique example of quantum entanglement in strongly interacting many-body systems. A heavy nucleus can split into hundreds of combinations of two complementary fragments in the fission process. The entanglement of fragment wave functions is persistent even after separation and impacts the partition of particles and energies between fragments. Based on microscopic dynamical calculations of the fission of 240 Pu, this work finds that dynamical quantum entanglement is indispensable in the appearance of sawtooth distributions of average excitation energies of fragments and thus neutron multiplicities, but not in average neutron excess of fragments. Both sawtooth slopes from particle-number projections are found to be steep – a feature which can be alleviated by random fluctuations. The persistent entanglement is mainly due to non-adiabatic dynamics since the final splitting is so fast that the non-localization of wave functions is kept during the separation. These findings may impact the understanding of quantum entanglement more broadly in mesoscopic systems

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Data-Driven Learning for the Mori--Zwanzig Formalism: A Generalization of the Koopman Learning Framework

A theoretical framework which unifies the conventional Mori--Zwanzig formalism and the approximate Koopman learning of deterministic dynamical systems from noiseless observation is presented. In this framework, the Mori--Zwanzig formalism, developed in statistical mechanics to tackle the hard problem of construction of reduced-order dynamics for high-dimensional dynamical systems, can be considered as a natural generalization of the Koopman description of the dynamical system. We next show that, similar to the approximate Koopman learning methods, data-driven methods can be developed for the Mori--Zwanzig formalism with Mori's linear projection operator. We have developed two algorithms to extract the key operators, the Markov and the memory kernel, using time series of a reduced set of observables in a dynamical system. We have adopted the Lorenz `96 system as a test problem and solved for the above operators. These operators exhibit complex behaviors, which are unlikely to be captured by traditional modeling approaches in Mori--Zwanzig analysis. The nontrivial generalized fluctuation-dissipation relationship, which relates the memory kernel with the two-time correlation statistics of the orthogonal dynamics, was numerically verified as a validation of the solved operators. Here we present numerical evidence that the generalized Langevin equation, a key construct in the Mori--Zwanzig formalism, is more advantageous in predicting the evolution of the reduced set of observables than the conventional approximate Koopman operators.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Construction of approximate invariants for non-integrable Hamiltonian systems

We present a method to construct high-order polynomial approximate invariants (AI) for non integrable Hamiltonian dynamical systems, and apply it to a modern ring-based particle accelerator. Taking advantage of a special property of one-turn transformation maps in the form of a square matrix, AIs can be constructed order-by-order iteratively. Evaluating AI with simulation data, we observe that AI’s fluctuation is actually a measure of chaos. Through minimizing the fluctuations, the stable region of long-term motions, i.e., the dynamic aperture of the accelerator, could be enlarged.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

PV-Fed Smart Inverters for Mitigation of Voltage and Frequency Fluctuations in Islanded Microgrids

The presence of low-inertia distributed generation (DG) units makes islanded microgrids vulnerable to voltage and frequency variations. Besides, the considerable difference between the inertia of synchronous-based and inverter-based DGs results in a power imbalance between supply and demand during abnormal conditions. As a result, both voltage and frequency of microgrid ac-bus start oscillating which might have severe consequences such as blackouts. This paper deploys the traditional controller of photovoltaic (PV) units to improve the dynamic behavior of islanded microgrids by suppressing the voltage and frequency fluctuations. To this end, an adaptive piecewise droop (APD) characteristic is proposed and employed in PV units to attain a faster balance between generation and consumption during transients, leading to an enhanced frequency response. Besides, the reactive-power control loop is equipped with a droop characteristic which enables the PV units to inject/absorb reactive power during transients and participate in voltage-profile enhancement of the system. Case study results are presented using PSCAD/EMTDC to confirm the validity of proposed method in improving the dynamics of islanded microgrids.

Pilehvar, Mohsen S.↗

DOE OSTI2020

Lone Pair Rotational Dynamics in Solids

Traditional classifications of crystalline phases focus on nuclear degrees of freedom. Through examination of both electronic and nuclear structure, we introduce the concept of an electronic plastic crystal. Such a material is classified by crystalline nuclear structure, while localized electronic degrees of freedom — here lone pairs — exhibit orientational motion at finite temperatures. This orientational motion is an emergent phenomenon arising from the coupling between electronic structure and polarization fluctuations generated by collective motions, such as phonons. Using ab initio molecular dynamics simulations, we predict the existence of electronic plastic crystal motion in halogen crystals and halide perovskites, and suggest that such motion may be found in a broad range of solids with lone pair electrons. Finally, such fluctuations in the charge density should be observable, in principle via synchrotron scattering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

An Information Theoretic Approach to Identify Dominant Voltage Influencers for Unbalanced Distribution Systems

Smart distribution grid with multiple renewable energy sources can experience random voltage fluctuations due to variable generation, which may result in voltage violations. Traditional voltage control algorithms are inadequate to handle fast voltage variations. Therefore, new dynamic control methods are being developed that can significantly benefit from the knowledge of dominant voltage influencer (DVI) nodes. DVI nodes for a particular node of interest refer to nodes that have a relatively high impact on the voltage fluctuations at that node. Conventional power flow-based algorithms to identify DVI nodes are computationally complex, which limits their use in real-time applications. This paper proposes a novel information theoretic voltage influencing score (VIS) that quantifies the voltage influencing capacity of nodes with DERs/active loads in a three phase unbalanced distribution system. VIS is then employed to rank the nodes and identify the DVI set. VIS is derived analytically in a computationally efficient manner and its efficacy to identify DVI nodes is validated using the IEEE 37-node test system. It is shown through experiments that KL divergence and Bhattacharyya distance are effective indicators of DVI nodes with an identifying accuracy of more than 90%. Additionally, the computation burden is also reduced by an order of 5, thus providing the foundation for efficient voltage control.

42 ENGINEERING↗

DOE OSTI2024

Studying Wind Loading on CSP Collectors to Improve Performance and Reliability

Electricity generation through Concentrating Solar-thermal Power (CSP) adds the advantage of thermal energy storage and heat production. CSP collector costs constitute nearly one-third of total plant costs. Managing wind loading, especially dynamic wind loading caused by the turbulent wind flow, is a significant design challenge. Traditional collector designs rely on wind tunnel experiments and numerical simulations, which do not fully capture dynamic effects. Here, we provide an overview of NREL's activities on wind loading on CSP collectors, aiming to improve their reliability and cost-efficiency. We present insights from field measurements at operational parabolic trough and power-tower CSP plants. Over two years, we measured turbulent wind fields and resulting structural loads in a parabolic trough plant, revealing wind, turbulence, and collector field interactions. Similar measurements are ongoing at the Crescent Dunes heliostat field. These campaigns show how atmospheric turbulence, wind direction, and collector orientation affect dynamic wind loading. Upwind collector structures can generate turbulent structures, causing fluctuating loads at downstream collectors, impacting fatigue lifetime and optical efficiency. Preliminary results indicate that interior collectors experience higher turning moments than those at the edges, suggesting higher drive wear. At Crescent Dunes, we use instrumented heliostats at the field's edge and interior to study the translation of turbulent wind to dynamic structural loads. The campaign is still ongoing, and we will present first findings. Another focus of our work is tying the experimental wind loading findings to optical performance of the collectors, using simulation tools. For example, we derive torsional errors of the parabolic troughs - the angular offset from the sun position along a collector row. This error is impacted by wind-induced mirror deformations. Future work will focus on continuing the measurement efforts and combining them with computational fluid dynamic simulations. By providing detailed field measurements and validated models, our efforts aim to enhance understanding of wind loading on CSP collectors, improving their structural integrity and optical performance.

CSP↗

DOE OSTI2021

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

From One- to Two-Magnon Excitations in the S = 3 / 2 Magnet β - CaCr 2 O 4

Here, we apply neutron spectroscopy to measure the magnetic dynamics in the S = 3 / 2 magnet β-CaCr 2 O 4 ( T N = 21 K ). The low-energy fluctuations, in the ordered state, resemble large- S linear spin waves from the incommensurate ground state. However, at higher energy transfers, these semiclassical and harmonic dynamics are replaced by an energy and momentum broadened continuum of excitations. Applying kinematic constraints required for energy and momentum conservation, sum rules of neutron scattering, and comparison against exact diagonalization calculations, we show that the dynamics at high-energy transfers resemble low- S one-dimensional quantum fluctuations. β-CaCr 2 O 4 represents an example of a magnet at the border between classical Néel and quantum phases, displaying dual characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Solid-state platform for cooperative quantum dynamics driven by correlated emission

While traditionally regarded as an obstacle to quantum coherence, recent breakthroughs in quantum optics have shown that the dissipative interaction of a qubit with its environment can be leveraged to protect quantum states and synthesize many-body entanglement. Inspired by this progress, here we set the stage for the—as yet uncharted—exploration of analogous cooperative phenomena in hybrid solid-state platforms. We develop a comprehensive formalism for the quantum many-body dynamics of an ensemble of solid-state spin defects interacting with the magnetic field fluctuations of a common solid-state reservoir. Our framework applies to any solid-state reservoir whose fluctuating spin, pseudospin, or charge degrees of freedom generate magnetic fields. To understand whether correlations induced by dissipative processes can play a relevant role in a realistic experimental setup, we apply our model to a qubit array interacting via the spin fluctuations of a ferromagnetic bath. Our results show that the low-temperature collective relaxation rates of the qubit ensemble can display clear signatures of super- and subradiance, i.e., forms of cooperative dynamics traditionally achieved in atomic ensembles. We find that the solid-state analog of these cooperative phenomena is robust against spatial disorder in the qubit ensemble and thermal fluctuations of the magnetic reservoir, providing a route for their feasibility in near-term experiments. Finally, our work lays the foundation for a multiqubit approach to quantum sensing of solid-state systems and the direct generation of many-body entanglement in spin-defect ensembles. Furthermore, we discuss how the tunability of solid-state reservoirs opens up novel pathways for exploring cooperative phenomena in regimes beyond the reach of conventional quantum optics setups.

NV centers↗

DOE OSTI2024

River–aquifer interactions enhancing evapotranspiration in a semiarid riparian zone: A modelling study

The hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three-dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high-frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia-Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river-stage fluctuations. Adding 6.64 × 10 6 m 3 year –1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5-m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Finally, our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

On a Critical Acceleration Scale of Dark Matter in ΛCDM and Dynamical Dark Energy

Abstract Universal acceleration a 0 emerges in various empirical laws, yet its fundamental nature remains unclear. Using Illustris and Virgo N -body simulations, we focus on the velocity and acceleration fluctuations in collisionless dark matter involving long-range gravity. For comparison, in the kinetic theory of gases, molecules undergo random elastic collisions involving short-range interactions, where only velocity fluctuations are relevant. Hierarchical structure formation proceeds through the merging of smaller halos to form larger halos, which facilitates a continuous energy cascade from small to large halos at a constant rate ε u ≈ −10 −7 m 2 s −3 . Velocity fluctuations involve a critical velocity u c ∝ (1 + z ) −3/4 . Acceleration fluctuations involve a critical acceleration a c ∝ (1 + z ) 3/4 . Two critical quantities are related by the rate of energy cascade ε u ≈ − a c u c /[2(3 π ) 2 ], where factor 3 π is from the angle of incidence during merging. With critical velocity u c on the order of 300 km s −1 at z = 0, the critical acceleration is determined to be a c 0 ≡ a c ( z = 0) ≈ 10 −10 m s −2 , suggesting a c might explain the universal acceleration a 0 ≈ 10 −10 m s −2 in the empirical Tully–Fisher relation or modified Newtonian dynamics. The redshift evolution a c ∝ (1 + z ) 3/4 is in good agreement with Magneticum and EAGLE simulations and in reasonable agreement with limited observations. This suggests a larger a 0 at a higher redshift such that galaxies of fixed mass rotate faster at a higher redshift. Note that for dark energy (DE) density ρ DE 0 ≈ a c 0 2 / G = 1 0 − 10 J m −3 , we postulate an entropic origin of the DE from acceleration fluctuations of dark matter, analogous to the gas pressure from velocity fluctuations. This leads to a dynamical DE coupled to the structure evolution involving a relatively constant DE density followed by a slow weakening phase, suggesting possible deviations from the standard ΛCDM paradigm.

N-body simulations↗