Results for “dynamical mean field theory” — Page 14

DOE OSTI2022

Stochastic exciton-scattering theory of optical line shapes: Renormalized many-body contributions

Spectral line shapes provide a window into the local environment coupled to a quantum transition in the condensed phase. In this paper, we build upon a stochastic model to account for non-stationary background processes produced by broad-band pulsed laser stimulation, as distinguished from those for stationary phonon bath. In particular, we consider the contribution of pair-fluctuations arising from the full bosonic many-body Hamiltonian within a mean-field approximation, treating the coupling to the system as a stochastic noise term. Herein, using the Itô transformation, we consider two limiting cases for our model, which lead to a connection between the observed spectral fluctuations and the spectral density of the environment. In the first case, we consider a Brownian environment and show that this produces spectral dynamics that relax to form dressed excitonic states and recover an Anderson–Kubo-like form for the spectral correlations. In the second case, we assume that the spectrum is Anderson–Kubo like and invert to determine the corresponding background. Using the Jensen inequality, we obtain an upper limit for the spectral density for the background. The results presented here provide the technical tools for applying the stochastic model to a broad range of problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1987

Unsteady aerodynamics of blade rows

The requirements placed on an unsteady aerodynamic theory intended for turbomachinery aeroelastic applications are discussed along with a brief description of the various theoretical models that are available to address these requirements. The main emphasis is placed on the description of a linearized inviscid theory which fully accounts for the effects of a nonuniform mean or steady flow on unsteady aerodynamic response. Although this theory has been developed primarily for blade flutter prediction, more general equations are presented which account for unsteady excitations due to incident external aerodynamic disturbances as well as those due to blade motions. The resulting equations consist of a system of three field equations along with conditions imposed at blade, wake and shock surfaces and in the far field. These equations can be solved to determine the fluctuations in all fluid dynamic properties throughout the required solution domain. Example solutions are presented to demonstrate several effects associated with nonuniform steady flows on the linearized unsteady aerodynamic response to prescribed blade motions.

Verdon, J. M.↗

NASA NTRS1986

Adiabaticity and viscosity in deep mantle convection

A study has been conducted of steady convection with adiabatic and viscous heating for variable viscosity in the Boussinesq limit using the mean-field theory. A strong nonlinear coupling is found between the thermodynamic constants governing adiabatic heating and the rheological parameters. The range of rheological values for which adiabaticity would occur throughout the mantle has been established. Too large an activation volume, greater than 6 cu cm/mol for the cases examined, would produce unreasonably high temperature at the bottom of the mantle (greater than 6000 K) and superadiabatic gradients, especially in the lower mantle. Radiogenic heating plays a profound role in controlling dynamically mantle temperatures. Present values for the averaged mantle heat production would yield objectionably high temperatures in the lower mantle.

Quareni, F.↗

DOE OSTI2020

Momentum-resolved analysis of condensate dynamic and Higgs oscillations in quenched superconductors with time-resolved ARPES

Higgs oscillations in nonequilibrium superconductors provide a unique tool to obtain information about the underlying order parameter. Several properties like the absolute value, existence of multiple gaps, and the symmetry of the order parameter can be encoded in the Higgs oscillation spectrum. Studying Higgs oscillations with time-resolved angle-resolved photoemission spectroscopy (ARPES) has the advantage over optical measurements that a momentum-resolved analysis of the condensate dynamic is possible. In this paper, we investigate the time-resolved spectral function measured in ARPES for different quench protocols. We find that analyzing amplitude oscillations of the ARPES intensity in the whole Brillouin zone allows to understand how the condensate dynamic contributes to the emerging of collective Higgs oscillations. Furthermore, by evaluating the phase of these oscillations, the symmetry deformation dynamic of the condensate can be revealed, which gives insight about the ground state symmetry of the system. With such an analysis, time-resolved ARPES experiments might be used in the future as a powerful tool in the field of Higgs spectroscopy.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Towards excitations and dynamical quantities in correlated lattices with density matrix embedding theory

Density matrix embedding theory (DMET) provides a framework to describe ground-state expectation values in strongly correlated systems, but its extension to dynamical quantities is still an open problem. We show one route to obtaining excitations and dynamical spectral functions by using the techniques of DMET to approximate the matrix elements that arise in a single-mode inspired excitation ansatz. We demonstrate this approach in the one-dimensional Hubbard model, comparing the neutral excitations, single-particle density of states, charge, and spin dynamical structure factors to benchmarks from the Bethe ansatz and density matrix renormalization group. Finally, our work highlights the potential of these ideas in building computationally efficient approaches for dynamical quantities.

1-dimensional systems↗

NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗

NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗

NASA NTRS1990

Measurement realities of current collection in dynamic space plasma environments

Theories which describe currents collected by conducting and non-conducting bodies immersed in plasmas have many of their concepts based upon the fundamentals of sheath-potential distributions and charged-particle behavior in superimposed electric and magnetic fields. Those current-collecting bodies (or electrodes) may be Langmuir probes, electric field detectors, aperture plates on ion mass spectrometers and retarding potential analyzers, or spacecraft and their rigid and tethered appendages. Often the models are incomplete in representing the conditions under which the current-voltage characteristics of the electrode and its system are to be measured. In such cases, the experimenter must carefully take into account magnetic field effects and particle anisotropies, perturbations caused by the current collection process itself and contamination on electrode surfaces, the complexities of non-Maxwellian plasma distributions, and the temporal variability of the local plasma density, temperature, composition and fields. This set of variables is by no means all-inclusive, but it represents a collection of circumstances guaranteed to accompany experiments involving energetic particle beams, plasma discharges, chemical releases, wave injection and various events of controlled and uncontrolled spacecraft charging. Here, an attempt is made to synopsize these diagnostic challenges and frame them within a perspective that focuses on the physics under investigation and the requirements on the parameters to be measured. Examples include laboratory and spaceborne applications, with specific interest in dynamic and unstable plasma environments.

Szuszczewicz, Edward P.↗

DOE OSTI2026

Time-resolved unsteady shock motion and spectral analysis of a hypersonic shock-wave/boundary-layer interaction over a double-cone

Direct experimental measurements of double-cone shock-wave/boundary-layer interaction dynamics at Mach 7.2 were collected using a schlieren-based shock-detection algorithm at the University of Texas at San Antonio. The points of flow separation and reattachment were spatially and spectrally evaluated at 300 kHz for time-resolved characterization of the unsteadiness of this flow field. The separation bubble was found to rapidly expand and contract about a large mean length. Spectral analysis of the velocity fluctuations of this motion revealed this behavior to be a turbulent phenomenon, with spectral energy trends consistent with the Kolmogorov theory for the turbulent energy cascade.

Glasby, Ryan [ORNL]↗

DOE OSTI2020

Low-temperature spin dynamics in the TmFeO 3 orthoferrite with a non-Kramers ion

We investigate the magnetic dynamics of the orthorhombic perovskite TmFeO 3 at low temperatures, below the spin reorientation transition at TSR≈80 K, by means of time-of-flight neutron spectroscopy. We find that the magnetic excitation spectrum combines two emergent collective modes associated with different magnetic sublattices. The Fe subsystem orders below TN~632 K into a canted antiferromagnetic structure and exhibits sharp, high-energy magnon excitations. We describe them using linear spin-wave theory, and reveal a pronounced anisotropy between in- and out-of-plane exchange interactions, which was mainly neglected in previous reports on the spin dynamics in orthoferrites. At lower energies, we find two crystalline electrical field (CEF) excitations of Tm3+ ions at energies of ~2 and 5 meV. In contrast to the sister compound YbFeO 3 , where the Yb 3+ ions form quasi-one-dimensional chains along the c axis, the Tm excitations show dispersion along both directions in the (0KL) scattering plane. Analysis of the neutron scattering polarization factor reveals a longitudinal polarization of the 2 meV excitation. To evaluate the effect of the CEF on the Tm 3+ ions, we perform point-charge model calculations, and their results quantitatively capture the main features of Tm single-ion physics, such as energies, intensities, and polarization of the CEF transitions, and the type of magnetic anisotropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1978

Three-dimensional magnetohydrodynamic turbulence in cylindrical geometry

A theory is developed in which a perfectly conducting cylindrical wall is used to describe incompressible MHD turbulence. Dynamic equations of incompressible MHD and its associated fields are expanded in a system of vector eigenfunctions of the curl. The total energy, magnetic helicity, and cross helicity have simple quadratic expressions in the expansion coefficients. The invariants are temporally constant for a truncation at a large, but finite, number of coefficients. The invariants are used to construct canonical distributions. It is noted that for all initial conditions, including quiescent ones, the mean-square velocity fields have finite values. The stability problem is transformed into a search for values of the integral variants which minimize the mean-square velocity fields. This leads to an extremal helicity principle.

Montgomery, D.↗

NASA NTRS2010

Downscaling Soil Moisture in the Southern Great Plains Through a Calibrated Multifractal Model for Land Surface Modeling Applications

Accounting for small-scale spatial heterogeneity of soil moisture (theta) is required to enhance the predictive skill of land surface models. In this paper, we present the results of the development, calibration, and performance evaluation of a downscaling model based on multifractal theory using aircraft!based (800 m) theta estimates collected during the southern Great Plains experiment in 1997 (SGP97).We first demonstrate the presence of scale invariance and multifractality in theta fields of nine square domains of size 25.6 x 25.6 sq km, approximately a satellite footprint. Then, we estimate the downscaling model parameters and evaluate the model performance using a set of different calibration approaches. Results reveal that small-scale theta distributions are adequately reproduced across the entire region when coarse predictors include a dynamic component (i.e., the spatial mean soil moisture ) and a stationary contribution accounting for static features (i.e., topography, soil texture, vegetation). For wet conditions, we found similar multifractal properties of soil moisture across all domains, which we ascribe to the signature of rainfall spatial variability. For drier states, the theta fields in the northern domains are more intermittent than in southern domains, likely because of differences in the distribution of vegetation coverage. Through our analyses, we propose a regional downscaling relation for coarse, satellite-based soil moisture estimates, based on ancillary information (static and dynamic landscape features), which can be used in the study area to characterize statistical properties of small-scale theta distribution required by land surface models and data assimilation systems.

Mascaro, Giuseppe↗

NASA NTRS1991

Particle motion in the tail current sheet

Theory of particle motion in current sheets is reviewed. For small, approximately constant normal magnetic field, Bz, particles oscillate about the current sheet and 'live' within the sheet for one-half gyroperiod based on Bz. This lifetime replaces the mean collision time in the Lorentzian conductivity and thus gives rise to the concept of an inertial (or gyro-) conductivity. A substorm model by Coroniti utilizes this conductivity to allow reconnection to proceed without anomalous processes, due to wave-particle interactions. Chaotic particle orbits may at times be important to the dynamics, depending on parameters such as particle energy, current sheet thickness, and field line curvature. A current sheet model with neutral line predicts a ridge structure and asymmetries in the distribution function. Ion distributions near the plasma sheet boundary layer, during the CDAW 6 interval, are consistent with the model predictions. In recent studies by Mitchell et al. and Williams et al., the major current carriers during the growth phase of a substorm were found to be adiabatic electrons not more than 1 keV, but just before a current disruption event, the tail current was mainly carried by energetic ions undergoing current sheet oscillation.

Speiser, T. W.↗

NASA NTRS1994

Toward a theory of interstellar turbulence. 1: Weak Alfvenic turbulence

We study weak Alfvenic turbulence of an incompressible, magnetized fluid in some detail, with a view to developing a firm theoretical basis for the dynamics of small-scale turbulence in the interstellar medium. We prove that resonant 3-wave interactions are absent. We also show that the Iroshnikov-Kraichnan theory of incompressible, magnetohydrodynamic turbulence -- which is widely accepted -- describes weak 3-wave turbulence; consequently, it is incorrect. Physical arguments, as well as detailed calculations of the coupling coefficients are used to demonstrate that these interactions are empty. We then examine resonant 4-wave interactions, and show that the resonance relations forbid energy transport to small spatial scales along the direction of the mean magnetic field, for both the shear Alfven wave and the pseudo Alfven wave. The three-dimensional inertial-range energy spectrum of 4-wave shear Alfven turbulence guessed from physical arguments reads E(k(sub z), k(sub perpendicular)) approximately V(sub A)v(sub L)L(exp -1/3)k(sub perpendicular)(exp -10/3), where V(sub A) is the Alfven speed, and v(sub L) is the velocity difference across the outer scale L. Given this spectrum, the velocity difference across lambda(sub perpendicular) approximately k(sub perpendicular exp -1) is v(sub lambda (sub perpendicular)) is approximately v(sub L)(lambda(sub perpendicular)/L)(exp 2/3). We derive a kinetic equation, and prove that this energy spectrum is a stationary solution and that it implies a positive flux of energy in k-space, along directions perpendicular to the mean magnetic field. Using this energy spectrum, we deduce that 4-wave interactions strengthen as the energy cascades to small, perpendicular spatial scales; beyond an upper bound in perpendicular wavenumber, k(sub perpendicular)L is approximately (V(sub A)/v(sub L))(exp 3/2), weak turbulence theory ceases to be valid. Energy excitation amplitudes must be very small for the 4-wave inertial-range to be substantial. When the excitation is strong, the width of the 4-wave inertial-range shrinks to zero. This seems likely to be the case in the interstellar medium.

Sridhar, S.↗

NASA NTRS1997

Effects of Buoyancy on Laminar and Turbulent Premixed V-Flame

Turbulent combustion occurs naturally in almost all combustion systems and involves complex dynamic coupling of chemical and fluid mechanical processes. It is considered as one of the most challenging combustion research problems today. Though buoyancy has little effect on power generating systems operating under high pressures (e.g., IC engines and turbines), flames in atmospheric burners and the operation of small to medium furnaces and boilers are profoundly affected by buoyancy. Changes in burner orientation impacts on their blow-off, flash-back and extinction limits, and their range of operation, burning rate, heat transfer, and emissions. Theoretically, buoyancy is often neglected in turbulent combustion models. Yet the modeling results are routinely compared with experiments of open laboratory flames that are obviously affected by buoyancy. This inconsistency is an obstacle to reconciling experiments and theories. Consequently, a fundamental understanding of the coupling between turbulent flames and buoyancy is significant to both turbulent combustion science and applications. The overall effect of buoyancy relates to the dynamic interaction between the flame and its surrounding, i.e., the so-called elliptical problem. The overall flame shape, its flowfield, stability, and mean and local burning rates are dictated by both upstream and downstream boundary conditions. In steady propagating premixed flames, buoyancy affects the products region downstream of the flame zone. These effects are manifested upstream through the mean and fluctuating pressure fields to influence flame stretch and flame wrinkling. Intuitively, the effects buoyancy should diminish with increasing flow momentum. This is the justification for excluding buoyancy in turbulent combustion models that treats high Reynolds number flows. The objectives of our experimental research program is to elucidate flame-buoyancy coupling processes in laminar and turbulent premixed flames, and to characterize microgravity (micro g) premixed flames. The results are used to derive appropriate scaling parameters for guiding the development of theoretical models to include the effects of buoyancy. Knowledge gain from the analysis will also contribute to further understanding of the elliptical nature of premixed flames. Our current emphasis is to examine the momentum limit above which the effects of buoyancy would become insignificant. This is accomplished by comparing the flowfields and the mean properties of normal gravity flames (+g), and reversed gravity flames (-g, up-side-down flames) at different flow velocities and turbulence intensities. Microgravity (micro g) flames experiments provide the key reference data to reconcile the differences between flames in +g and -g. As flame configuration has significant impact on premixed flames characteristics we have studied axi-symmetric conical flames and plane-symmetric rod-stabilized v-flames. The two configurations produce distinct features that dictates how the flames couple with buoyancy. In a conical flame, the hot products plume completely envelopes the flame cone and shields the flame from direct interaction with the ambient air. The plume originates at the burner rim and generates a divergent flowfield. In comparison, the products region of v-flames forms between the twin flame sheets and it is convergent towards the center-plane. Interaction with ambient air is limited to the two end regions of the stabilized rod and beyond the flame sheets.

Cheng, Robert K.↗

NASA NTRS1991

A theory of neutral gas emissions from a plasma contactor and its effect on electrodynamic tether performance

An existing anisotropic anodic plasma contactor model is extended to include neutral gas emissions. Within the framework of the contactor model, the inclusion of external ionization leads to an integrodifferential equation. By choosing the ratio of the contactor radius to the neutral gas-electron ionization mean free path at the contactor exit as a suitable expansion parameter, an analytical expression for the current-voltage characteristic can be found. The performance of this contactor model is examined via a dynamic simulation model developed for electrodynamic tethers. Detailed models of the geomagnetic field and the ionosphere are used to create a realistic environment. The highly nonlinear contactor voltagle-current characteristic is incorporated into a circuit equation, which includes radiation impedance, and is solved along the tether's orbit. Results show that it is more effective to fully ionize the contactor gas internally, than to partially ionize it externally, and that based on the specific power for the tether system, the optimum gas to use is argon. In addition, the effect of the radiation impedance on tether system performance is examined.

Roy, R. I. S.↗

DOE OSTI2022

Feasibility study of Mg storage in a bilayer silicene anode via application of an external electric field

With the goal of developing a Si-based anode for Mg-ion batteries (MIBs) that is both efficient and compatible with the current semiconductor industry, the current research utilized classical Molecular Dynamics (MD) simulation in investigating the intercalation of a Mg 2+ ion under an external electric field (E-field) in a 2D bilayer silicene anode (BSA). First principles density functional theory calculations were used to validate the implemented EDIP potentials. Our simulation shows that there exists an optimum E-field value in the range of 0.2–0.4 V Å –1 for Mg 2+ intercalation in BSA. To study the effect of the E-field on Mg 2+ ions, an exhaustive spread of investigations was carried out under different boundary conditions, including calculations of mean square displacement (MSD), interaction energy, radial distribution function (RDF), and trajectory of ions. Our results show that the Mg 2+ ions form a stable bond with Si in BSA. The effects of E-field direction and operating temperature were also investigated. In the X–Y plane in the 0°–45° range, 15° from the X-direction was found to be the optimum direction for intercalation. The results of this work also suggest that BSA does not undergo drastic structural changes during the charging cycles with the highest operating temperature being ~300 K

25 ENERGY STORAGE↗

DOE OSTI2024

Coarse-graining Hamiltonian systems using WSINDy

Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.

97 MATHEMATICS AND COMPUTING↗