Search NASA⌕ Search

SEARCH · Search NASA

Results for “time-dependent basis”

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 19 records

Hybrid Basis and Multi-Center Grid Method for Strong-Field Processes

We present a time-dependent framework that combines a hybrid basis, consisting of Gaussian-type orbitals (GTOs) and finite-element discrete-variable representation (FEDVR) functions, with a multicenter grid to simulate strong-field and attosecond dynamics in atoms and molecules. The method incorporates the construction of the orthonormal hybrid basis, the evaluation of electronic integrals, a unitary time-propagation scheme, and the extraction of optical and photoelectron observables. Its accuracy and robustness are benchmarked on one-electron systems such as atomic hydrogen and the dihydrogen cation (H$^+_2$) through comparisons with essentially-exact reference results for bound-state energies, high-harmonic generation spectra, photoionization cross sections, and photoelectron momentum distributions. This work establishes the groundwork for its integration with quantum-chemistry methods, which is already operational but will be detailed in future work, thereby enabling ab initio simulations of correlated polyatomic systems in intense ultrafast laser fields.

74 ATOMIC AND MOLECULAR PHYSICS↗

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

FAST-1.2.2: A Computer Code for Thermal-Mechanical Nuclear Fuel Analysis under Steady-state and Transients

Fuel Analysis under Steady-state and Transients (FAST) is the U.S. Nuclear Regulatory Commission (NRC)’s computer code that calculates the steady-state and transient response of nuclear reactor fuel rods during long-term in-reactor burnup, anticipated operational occurrences (AOOs), design basis accidents (DBAs), and dry storage conditions. The code calculates the temperature, pressure, and deformation of a fuel rod as functions of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by the code include heat conduction through the fuel and other materials, heat transfer from the cladding-to-coolant, cladding elastic and plastic deformation (including creep), fuel-cladding mechanical interaction, fission gas release from the fuel, rod internal pressure, void volume, and cladding oxidation. The code contains necessary material and coolant properties, as well as clad-to-coolant heat transfer correlations, for normal operation through postulated accidents and AOOs for today’s U.S.-based light water reactor (LWR) fuel designs. FAST-1.2.2 also contains preliminary materials and models for new LWR fuel concepts, such as accident tolerant fuel (ATF), and non-LWR fuel concepts such as metallic fuels for sodium fast reactors (SFRs). FAST has been developed for use on Windows and Linux operating systems. This document describes FAST-1.2.2 and is one of a series of documents on the code; the other documents detail the material properties used by FAST as well as its integral assessment to experiments and commercial data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An Optimization-Based Coupling of Reduced Order Models with an Efficient Reduced Adjoint Basis Generation Approach

Optimization-based coupling (OBC) is an attractive alternative to traditional Lagrange multiplier approaches in multiple modeling and simulation contexts. However, application of OBC to time-dependent problems has been hindered by the computational cost of finding the stationary points of the associated Lagrangian, which requires primal and adjoint solves. This issue can be mitigated by using OBC in conjunction with computationally efficient reduced order models (ROMs). To demonstrate the potential of this combination, in this paper, we develop an optimization-based ROM-ROM coupling for a transient advection-diffusion transmission problem. We pursue the “optimize-then-reduce” path toward solving the minimization problem at each time step and solve reduced space adjoint system of equations, where the main challenge in this formulation is the generation of adjoint snapshots and reduced bases for the adjoint systems required by the optimizer. One of the main contributions of the paper is a new technique for an efficient adjoint snapshot collection for gradient-based optimizers in the context of optimization-based ROM-ROM couplings. In conclusion, we present numerical studies demonstrating the accuracy of the approach along with comparison between various approaches for selecting a reduced order basis for the adjoint systems, including decay of snapshot energy, average iteration counts, and timings.

coupled problems↗

KTaO 3 -Based Supercurrent Diode

The supercurrent diode effect (SDE), characterized by nonreciprocal critical currents, represents a promising building block for future dissipationless electronics and quantum circuits. Realizing SDE requires breaking both time-reversal and inversion symmetry in the device. Here we use conductive atomic force microscope (c-AFM) lithography to pattern reconfigurable superconducting weak links (WLs) at the LaAlO 3 /KTaO 3 (LAO/KTO) interface. By deliberately engineering the WL geometry at the nanoscale, we realize SDE in these devices in the presence of modest out-of-plane magnetic fields. The SDE polarity can be reversed by simply changing the WL position, and the rectification efficiency reaches up to 13% under optimal magnetic field conditions. Time-dependent Ginzburg–Landau simulations reveal that the observed SDE originates from asymmetric vortex motion in the inversion-symmetry-breaking device geometry. This demonstration of SDE in the LAO/KTO system establishes a versatile platform for investigating and engineering vortex dynamics, forming the basis for engineered quantum circuit elements.

C-AFM lithography↗

Leveraging operator learning to accelerate convergence of the preconditioned conjugate gradient method

We propose a new deflation strategy to accelerate the convergence of the preconditioned conjugate gradient (PCG) method for solving parametric large-scale linear systems of equations. Unlike traditional deflation techniques that rely on eigenvector approximations or recycled Krylov subspaces, we generate the deflation subspaces using operator learning, specifically the Deep Operator Network (DeepONet). To this aim, we introduce two complementary approaches for assembling the deflation operators. The first approach approximates near-null space vectors of the discrete PDE operator using the basis functions learned by the DeepONet. The second approach directly leverages solutions predicted by the DeepONet. To further enhance convergence, we also propose several strategies for prescribing the sparsity pattern of the deflation operator. Here, a comprehensive set of numerical experiments encompassing steady-state, time-dependent, scalar, and vector-valued problems posed on both structured and unstructured geometries is presented and demonstrates the effectiveness of the proposed DeepONet-based deflated PCG method, as well as its generalization across a wide range of model parameters and problem resolutions.

Deflation↗

When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu 14 and Au 20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps

In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.

Chemical calculations↗

Conformation-specific synthetic intrabodies modulate mTOR signaling with subcellular spatial resolution

Subcellular compartmentalization is integral to the spatial regulation of mechanistic target of rapamycin (mTOR) signaling. However, the biological outputs associated with location-specific mTOR signaling events are poorly understood and challenging to decouple. Here, we engineered synthetic intracellular antibodies (intrabodies) that are capable of modulating mTOR signaling with genetically programmable spatial resolution. Epitope-directed phage display was exploited to generate high affinity synthetic antibody fragments (Fabs) against the FKBP12–Rapamycin binding site of mTOR (mTOR FRB ). We determined high-resolution crystal structures of two unique Fabs that discriminate distinct conformational states of mTOR FRB through recognition of its substrate recruitment interface. By leveraging these conformation-specific binders as intracellular probes, we uncovered the structural basis for an allosteric mechanism governing mTOR complex 1 (mTORC1) stability mediated by subtle structural adjustments within mTOR FRB . Furthermore, our results demonstrated that synthetic binders emulate natural substrates by employing divergent yet complementary hydrophobic residues at defined positions, underscoring the broad molecular recognition capability of mTOR FRB . Intracellular signaling studies showed differential time-dependent inhibition of S6 kinase 1 and Akt phosphorylation by genetically encoded intrabodies, thus supporting a mechanism of inhibition analogous to the natural product rapamycin. Finally, we implemented a feasible approach to selectively modulate mTOR signaling in the nucleus through spatially programmed intrabody expression. These findings establish intrabodies as versatile tools for dissecting the conformational regulation of mTORC1 and should be useful to explore how location-specific mTOR signaling influences disease progression.

Science & Technology - Other Topics↗

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗

Plasmon-Assisted Electrochemical Epoxidation using Water as an Oxidant

Olefin epoxidation, an important industrial reaction, often uses hazardous oxidants, causing challenges in waste disposal. Here, we demonstrate the use of water as an oxidant by a plasmon-assisted electrochemical strategy. The electrocatalyst is comprised of a hybrid of a water oxidation catalyst, manganese oxide, and plasmonic gold nanoparticles. Visible-light irradiation of the electrocatalyst enhanced the epoxidation of 4-styrenesulfonate 5-fold as compared to dark conditions at the same temperature. From electrochemical analyses conducted under plasmon excitation conditions, complemented by real-time time-dependent density functional tight binding simulations, it is found that the plasmonic boost of the electrochemical styrene epoxidation is due to energetic holes generated by the excitation of localized surface plasmon resonances of gold nanoparticles. These photogenerated holes activate adsorbed water for oxidation and enhance the binding of 4-styrenesulfonate at interfacial sites. Here, this work demonstrates a proof of concept and establishes the mechanistic basis for plasmon-assisted activation of water as an O atom source for electrochemical epoxidations.

Gold↗

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗

Overview of the physics basis for the ARC fusion power plant

Commonwealth Fusion Systems plans to build ARC as the first fusion power plant at a site in Chesterfield County, Virginia, USA by the early 2030s. We present an overview of analysis comprising the physics basis of the ARC V3A design, a high-magnetic-field tokamak with 𝐵 0 =11.4 T, 𝐼 𝑝 =12.0 MA, 𝑅 0 =4.62 m, 𝑎 =1.18 m. ARC V3A is designed to produce 𝑃 𝑓⁡𝑢⁢𝑠 ≈1.13 GW DT fusion power and deliver ⩾400 MW net electric power to the grid. This overview includes quantitative analysis of fundamental issues for design of and operational plasma scenarios for a tokamak power plant, and lays out the design targets and strategic choices for ARC, including empirical fusion performance projections, assessment of H-mode access, ion cyclotron resonance heating simulations, alpha particle physics and time-dependent full-pulse simulations. This is complemented by topical papers on fusion performance and transport, disruption physics, boundary physics and magnetohydrodynamic stability. Critically, these studies identify key model uncertainties and physics risks to be retired through SPARC operation. Due to the modular nature of ARC, early results from SPARC can be incorporated into the design of the first ARC as well as subsequent replacements of the ARC vacuum vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Statistical Uncertainty of Inhalation Dose Coefficients: Impact of Particle Deposition in ICRP 66 Human Respiratory Tract Model

Inhaled radioactive materials can pose a long-term health concern, as the material can be incorporated into the body’s metabolic pathways and remain in organs and tissues for extended durations. During the retention period, the radioactive material may localize in a source organ and irradiate adjacent target organs and tissues. Distribution of these materials changes over time, requiring biokinetic modeling to evaluate their movement through various tissues and organs. The evolving distribution depends on multiple inputs characterizing the inhaled material, such as particle size and size distribution, particle density, aspect ratio, specific radionuclide, the chemical form, and solubility. In addition, biological parameters such as breathing rate, breathing type (nasal or nasal/oral), respiratory system morphometry, tidal volume, functional residual capacity, and anatomical dead space all influence material transport. These aerosol properties and physiological characteristics of the respiratory tract jointly define a range of initial conditions that influence the time-dependent distribution of radioactive material. To evaluate both uncertainty in the initial conditions of inhalation exposure and the final output (committed effective dose) from biokinetic models, a Python-based software tool, Radiological Exposure Dose Calculator (REDCAL), was developed to propagate uncertainty within the human respiratory tract model. Focusing on deposition fraction uncertainty, the primary objective was to characterize the initial activity distribution across respiratory regions as a function of anticipated particle sizes and distributions. The impact of the deposition fraction uncertainty was propagated to committed effective dose coefficients for selected radionuclides in a companion publication. For each particle size, a lognormal distribution, characterized by its geometric mean as defined within ICRP Publication 66, serves as the basis for introducing uncertainty into the physical processes governing deposition in various lung regions. Finally, this study addresses the deposition process and examines how uncertainty in deposition mechanisms affects activity distribution in the airways, ultimately presenting the expected range and standard deviation of deposited activity as a function of particle size.

International Commission on Radiological Protectio↗

Representativity error scaling of models of the high temperature test facility

Error scaling is a critical step toward the validation of modeling capabilities for high-temperature gas-cooled reactors (HTGRs). Extensive effort is being made to bring HTGRs into the validation basis of numerous thermal-hydraulics codes. Here, this paper demonstrates how systems-level codes can be leveraged to perform error scaling analyses between an experimental facility and a plant-to-be facility. Specifically, we focus on two conduction cooldown experiments from the High Temperature Test Facility (HTTF) and the General Atomics 350 MW th modular high-temperature gas-cooled reactor (MHTGR-350). The error scaling methodology employed in this study is representativity, which in the context of this work was used to quantify how well experiments captured the physics of the plant facility by comparing sensitivity vectors between the two facilities. In addition, two different RELAP5–3D models of the HTTF were included in the comparison to determine if modeling methodologies notably impact the representativity results. The key figures of merit are the maximum block temperature and the coolant outlet temperature. The time-dependent maximum block temperature had a low representativity of below 0.1 for both models across both transients. The coolant outlet temperature had a higher representativity of around 0.6 for both models during the pressurized conduction cooldown, but it was below 0.2 for the depressurized conduction cooldown. Overall, the HTTF experiments were not representative of the transients in the MHTGR-350. However, these results can play a significant role in informing future potential experiments for HTGR systems that iterate on what the HTTF accomplished.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

MOSAiC studies of long-lasting mixed-phase cloud events and analysis of the liquid-phase properties of Arctic clouds

Vertically resolved observations of the temporal evolution of mixed-phase clouds (MPCs) were performed over the central Arctic during the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition, which lasted from October 2019 to September 2020. The research icebreaker Polarstern , drifting with the pack ice for more than 7 months, mostly at latitudes > 85° N, served as a platform for state-of-the-art remote sensing of aerosols and clouds. The use of the recently introduced dual field-of-view (FOV) polarization lidar technique in combination with the well-established lidar-radar retrieval technique provided, for the first time, a robust instrumental basis to monitor the evolution of the liquid and the ice phase of MPCs and the interplay between the two phases. Two long-lasting Arctic MPC events observed close to the North Pole in mid-winter (December 2019) and late summer (September 2020) are discussed to provide new insight into Arctic MPC evolution processes. In the second part of the article, cloud statistics, covering all seasons of a year, are presented. The focus is on the optical and microphysical properties of the liquid phase. These results are solely derived from the dual-FOV lidar observations. The key findings of the study can be summarized as follows: persistent activation of aerosol particles to form water droplets is of great importance for the longevity of MPCs. The observations confirm that ice formation occurs predominantly via immersion freezing. The field studies suggest that the free tropospheric reservoirs of cloud condensation nuclei (CCN) and of ice-nucleating particles (INPs) were always well filled, i.e., the clouds did not exhaust their supply of activatable and activated particles. The observation of long-lasting MPC events, low ice production rates, and a sufficiently large INP reservoir leads to the recommendation to use a time-dependent immersion freezing parameterization in MPC modeling efforts.

Jimenez, Cristofer [Leibniz Inst. for Tropospheric↗

Radiological Releases from Novel Fuel Forms in Advanced Reactors During Severe Accidents for Consequence Analyses

Various advanced reactor developers are exploring the potential for reductions in the size of physical security forces and emergency planning zones. These reductions are based on robust fuel forms and inherently safe reactor designs. However, such reductions in physical protection measures could increase the risk of sabotage. To assess the possibility of reducing these measures, sabotage-induced radiological consequence analyses were carried out. These analyses considered accident scenarios that were beyond design basis accidents and overly conservative (Shah, 2025a; Shah, 2025b; Shah and Hartanto, 2026), yielding very large release fractions. These fractions, which can be used to evaluate physical protection and emergency planning requirements, have been crudely determined and applied as demonstrations for a sodium-cooled fast reactor (SFR) (Shah and Hartanto, 2025a), a high-temperature gas-cooled reactor (HTGR) (Shah and Hartanto, 2025b), a heat pipe–cooled reactor (HPR) (Shah and Hartanto, 2025c), and a molten salt–cooled reactor (MSR) (Shah et al., 2026). A Sandia National Laboratories (SNL) team used MELCOR—a fully integrated severe accident analysis code—to demonstrate the code’s capability to analyze advanced (i.e., not light water–cooled) reactors (including a fluoride salt–cooled high-temperature reactor [FHR]) and calculate radiological releases to the environment during severe accidents (Wagner et al., 2022a, 2022b, 2022c, 2023a, and 2023b). Although the analyses were carried out to demonstrate MELCOR’s growing capability, the release source terms were estimated for advanced reactors, providing valuable insights into the accident progression and radiological releases. These findings from prior SNL studies, including estimated source terms and related sensitivity studies, were leveraged to derive source terms for postulated sabotage-induced accidents. Insights from these sensitivity studies informed the scaling of SNL’s estimated source terms for the defined accident scenarios. The derived release fractions for the severe accident scenarios for the respective reactor designs can be used to perform more nuanced dose consequence analyses to evaluate the reactors’ physical protection and emergency planning zone requirements. These analyses are in accordance with the risk-informed, performance-based approach proposed under 10 CFR Part 53. This study builds on the prior source term analyses and associated sensitivity studies by SNL to derive time-dependent and design-informed release fractions. Section 2 describes the diverse advanced reactor designs analyzed by the SNL team. Section 3 discusses the severe accident analyses, the release fractions calculated, and the limitations and assumptions of the demonstration project. Section 4 presents the release percentages derived for the hypothetical sabotage-induced severe accidents at the advanced reactors. Section 5 summarizes the study’s findings and conclusions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Systematic multi-machine analysis of the exhaust time-dependent behavior in tokamaks

The understanding of the time-scales and associated transient behavior of fusion exhaust plasmas plays a crucial role in its dynamic modeling and its control. This work presents an overview of experimental investigations of the exhaust dynamics in TCV, MAST-U, ASDEX-Upgrade, WEST, DIII-D, and JET. From the presented experiments, a clear picture arises on properties of the exhaust dynamics across machines. Particularly, we observe that the scrape-off layer equilibrates on fast time-scales ($>$ 70 Hz) and that exhaust dynamics measured in response to gas valve modulations mostly behave smoothly and linearly, with similarities across devices, across scenarios (H-mode, L-mode), injected species, and injection locations. The measurements presented have formed the basis for systematic exhaust control on the considered devices. We now present this database for the essential validation of dynamic exhaust models for reactor design and control.

control↗

Floquet Engineering of Interactions and Entanglement in Periodically Driven Rydberg Chains

Neutral atom arrays driven into Rydberg states constitute a promising approach for realizing programmable quantum systems. Enabled by strong interactions associated with Rydberg blockade, they allow for simulation of complex spin models and quantum dynamics. We introduce a new Floquet engineering technique for systems in the blockade regime that provides control over novel forms of interactions and entanglement dynamics in such systems. Our approach is based on time-dependent control of Rydberg laser detuning and leverages perturbations around periodic many-body trajectories as resources for operator spreading. These time-evolved operators are utilized as a basis for engineering interactions in the effective Hamiltonian describing the stroboscopic evolution. As an example, we show how our method can be used to engineer strong spin exchange, consistent with the blockade, in a one-dimensional chain, enabling the exploration of gapless Luttinger liquid phases. In addition, we demonstrate that combining gapless excitations with Rydberg blockade can lead to dynamic generation of large-scale multipartite entanglement. Experimental feasibility and possible generalizations are discussed.

Floquet systems↗