Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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

At least 937 records · Page 52

Discrete generative diffusion models without stochastic differential equations: A tensor network approach

Diffusion models (DMs) are a class of generative machine learning methods that sample a target distribution by transforming samples of a trivial (often Gaussian) distribution using a learned stochastic differential equation. In standard DMs, this is done by learning a “score function” that reverses the effect of adding diffusive noise to the distribution of interest. Here we consider the generalisation of DMs to lattice systems with discrete degrees of freedom, and where noise is added via Markov chain jump dynamics. We show how to use tensor networks (TNs) to efficiently define and sample such “discrete diffusion models” (DDMs) without explicitly having to solve a stochastic differential equation. We show the following: (i) by parametrising the data and evolution operators as TNs, the denoising dynamics can be represented exactly; (ii) the auto-regressive nature of TNs allows to generate samples efficiently and without bias; (iii) for sampling Boltzmann-like distributions, TNs allow to construct an efficient learning scheme that integrates well with Monte Carlo. We illustrate this approach to study the equilibrium of two models with non-trivial thermodynamics, the d = 1 constrained Fredkin chain and the d = 2 Ising model. Published by the American Physical Society 2025

Causer, Luke (ORCID:0000000194243473)↗

Symmetrizing relativistic three-body partial wave amplitudes

S matrix principles and symmetries impose constraints on three-particle scattering amplitudes, which can be formulated as a class of integral equations for their partial wave projections. However, these amplitudes are typically expressed in an asymmetric basis, where one of the initial and final state particles is singled out, and all quantum numbers are defined relative to this spectator. In this work, we show how to construct symmetric partial wave amplitudes, which have been symmetrized over all possible spectator combinations, using their asymmetric counterparts and sets of recoupling coefficients. We derive these recoupling coefficients for arbitrary angular momentum and isospin for arbitrary systems of spinless particles with SU(2) flavor symmetry. We propose a simple intensity observable suitable for visualizing the structure of three-body dynamics in Dalitz distributions. Finally, we provide some numerical examples of Dalitz distributions relevant for future lattice QCD calculations of 3⁢𝜋 systems and provide numerical evidence that the symmetrization procedure is consistent with expected symmetries of Dalitz plots.

Lattice QCD↗

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↗

Towards Improving luminosity using optics tuning and data-driven methods

The results of Run 24 experiments at Relativistic Heavy Ion Collider (RHIC) for improving luminosity using optics tuning are presented in this study. In the first experiment, MADx matching was used to output magnet strengths corresponding to specific s star movements around Interaction Region 8 (IR8). The corresponding Zero Degree Calorimeter (ZDC) signal was measured in place of luminosity, and Bayesian Optimization aids search of optimal movements. It was found that values retrieved from matching were inaccurate, resulting in negative feedback loops. The second experiment focused on calculating accurate s star movements. The matching method was replaced with a linear sensitivity matrix, directly relating optics to power supply, and its null space was used to fit constraints such as hysteresis effects. At the experiment, beam losses were observed at collimators around boundary of IR8, which were fixed for the third experiment. Dynamic mode decomposition was also introduced to improve quality of turn-by-turn (TBT) data as well as accuracy and consistency of optics measurements at IR8. These improvements will be tested in the experiment of next RHIC run for luminosity optimization.

Accelerator Physics↗

Enabling Predictive Scale-Bridging Simulations through Active Learning (Institutional Computing Annual Report (Project w21_alscalebridging)) [Slides]

The goal of this project was to develop, demonstrate, and provide a new capability to achieve greater physical fidelity in large-scale simulations, rather than the usual brute-force increases in the number of mesh elements or particles. Transport in nanoporous media, critical to hydrocarbon extraction from tight shale formations, is affected by molecular-level interactions. Several coarse-scale Lattice Boltzmann model (LBM) parameters cannot be directly computed, so instead we calibrate them to molecular dynamics (MD) simulations, by building emulators that mimic MD and LBM behavior and then training an upscaler. The resulting machine learning (ML)-based scale-bridging framework is up to 7 orders of magnitude faster than direct MD.

74 ATOMIC AND MOLECULAR PHYSICS↗

Varied phenomenology of models displaying dynamical large-deviation singularities

Singularities of dynamical large-deviation functions are often interpreted as the signal of a dynamical phase transition and the coexistence of distinct dynamical phases, by analogy with the correspondence between singularities of free energies and equilibrium phase behavior. Here we study models of driven random walkers on a lattice. These models display large-deviation singularities in the limit of large lattice size, but the extent to which each model's phenomenology resembles a phase transition depends on the details of the driving. We also compare the behavior of ergodic and nonergodic models that present large-deviation singularities. Finally, we argue that dynamical large-deviation singularities indicate the divergence of a model timescale, but not necessarily one associated with cooperative behavior or the existence of distinct phases.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dynamically reconfigurable topological routing in nonlinear photonic systems

The propagation path of topologically protected states is bound to the interface between regions with different topology, and as such, the functionality of linear photonic devices leveraging these states is fixed during fabrication. Here, we propose a mechanism for dynamic control over a driven dissipative system’s local topology, yielding reconfigurable topological interfaces and thus tunable paths for protected routing. We illustrate our approach in non-resonantly pumped polariton lattices, where the nonlinear interaction between the polaritons and the exciton reservoir due to non-resonant pumping can yield picosecond-scale changes in the propagation paths of the chiral edge states. To analytically confirm the numerically observed topological dynamics, we generalize the spectral localizer framework to non-linear non-Hermitian Chern materials and apply this framework to a continuous model of the polariton system based on a driven-dissipative Gross-Pitaevskii equation. In doing so, we show that the local changes in the polariton lattice’s topology are captured by a local Chern marker. Looking forward, we anticipate such reconfigurable topological routing will enable the realization of novel classes of topological photonic devices.

Wong, Stephan [Sandia National Laboratories (SNL-N↗

Modeling of standing wave RF cavities for tracking through multi-pass energy recovery linac

Short bunches, high current and multiple linac pass are all characteristics of Energy Recovery Linacs (ERLs), which may result in collective effects. They in turn, may affect the beam, degrading its quality, or even yield to instabilities causing a beam loss. To study and mitigate these effects one needs a numerical simulation code, that can take into account both the collective effects, as well as, particular ERL features, such as a multi-turn design that does not reach a steady state or the multiple passages of the beam through Radio-Frequency (RF) cavities at different energies. CODAL [1], a code developed by SOLEIL in collaboration with IJCLab, enables such studies. It is a 6 dimensional (6D) tracking code applying 'kicks' based on the integration of the local Hamiltonian for each element of the lattice. It is also capable of simulating space charge, wakefields and coherent synchrotron radiation. However, to correctly take into account the ERL dynamics, an upgrade had to be made to include the effect of a standing wave RF cavity in 6D. In this paper, we will concentrate on the implementation and benchmarking (with DESY’s tracking code ASTRA [2]) of both the longitudinal and the transverse models (by J.B. Rosenzweig and L. Serafini [3]), which we use to carry out tracking of fully analytical 6D RF cavity.

Bogacz, A.↗

Quantifying Dynamic Signal Spread in Real-Time High-Energy X-ray Diffraction

Abstract Measured intensity in high-energy monochromatic X-ray diffraction (HEXD) experiments provides information regarding the microstructure of the crystalline material under study. The location of intensity on an areal detector is determined by the lattice spacing and orientation of crystals so that changes in the heterogeneity of these quantities are reflected in the spreading of diffraction peaks over time. High temporal resolution of such dynamics can now be experimentally observed using technologies such as the mixed-mode pixel array detector (MM-PAD) which facilitates in situ dynamic HEXD experiments to study plasticity and its underlying mechanisms. In this paper, we define and demonstrate a feature computed directly from such diffraction time series data quantifying signal spread in a manner that is correlated with plastic deformation of the sample. A distinguishing characteristic of the analysis is the capability to describe the evolution from the distinct diffraction peaks of an undeformed alloy sample through to the non-uniform Debye–Scherrer rings developed upon significant plastic deformation. We build on our previous work modeling data using an overcomplete dictionary by treating temporal measurements jointly to improve signal spread recovery. We demonstrate our approach in simulations and on experimental HEXD measurements captured using the MM-PAD. Our method for characterizing the temporal evolution of signal spread is shown to provide an informative means of data analysis that adds to the capabilities of existing methods. Our work draws on ideas from convolutional sparse coding and requires solving a coupled convex optimization problem based on the alternating direction method of multipliers.

Engineering↗

Ultrafast infrared nano-imaging of far-from-equilibrium carrier and vibrational dynamics

Ultrafast infrared nano-imaging has demonstrated access to ultrafast carrier dynamics on the nanoscale in semiconductor, correlated-electron, or polaritonic materials. However, mostly limited to short-lived transient states, the contrast obtained has remained insufficient to probe important long-lived excitations, which arise from many-body interactions induced by strong perturbation among carriers, lattice phonons, or molecular vibrations. Here, we demonstrate ultrafast infrared nano-imaging based on excitation modulation and sideband detection to characterize electron and vibration dynamics with nano- to micro-second lifetimes. As an exemplary application to quantum materials, in phase-resolved ultrafast nano-imaging of the photoinduced insulator-to-metal transition in vanadium dioxide, a distinct transient nano-domain behavior is quantified. In another application to lead halide perovskites, transient vibrational nano-FTIR spatially resolves the excited-state polaron-cation coupling underlying the photovoltaic response. These examples show how heterodyne pump-probe nano-spectroscopy with low-repetition excitation extends ultrafast infrared nano-imaging to probe elementary processes in quantum and molecular materials in space and time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Preparations for quantum simulations of quantum chromodynamics in 1+1 dimensions. I. Axial gauge

Here, tools necessary for quantum simulations of 1+1 dimensional quantum chromodynamics are developed. When formulated in axial gauge and with two flavors of quarks, this system requires 12 qubits per spatial site with the gauge fields included via nonlocal interactions. Classical computations and D-wave’s quantum annealer advantage are used to determine the hadronic spectrum, enabling a decomposition of the masses and a study of quark entanglement. Color “edge” states confined within a screening length of the end of the lattice are found. IBM’s seven-qubit quantum computers, ibmq_jakarta and ibm_perth, are used to compute dynamics from the trivial vacuum in one-flavor QCD with one spatial site. More generally, the Hamiltonian and quantum circuits for time evolution of 1+1 dimensional SU(N c ) gauge theory with Nf flavors of quarks are developed, and the resource requirements for large-scale quantum simulations are estimated.

1+1 dimensions↗

Driving Force Dependent Photoinduced Charge Transfer Dynamics in Polymer-Wrapped Semiconducting Single-Walled Carbon Nanotubes

Here, we investigate the thermodynamic driving-force dependences of photoinduced charge separation (CS) and subsequent charge transfer dynamics in single-walled carbon nanotube (SWNT)–perylenediimide (PDI) donor–acceptor (D–A) superstructures. Pump–probe spectroscopy reveals that [SWNT (•+)n ]-(PDI –• ) n CS states form on an ∼100 fs time scale following photoexcitation; these dynamics are invariant across an ∼400 mV driving force range, indicating that SWNT hole polaron formation time scales are determined by nanotube lattice and solvent relaxation. These CS states feature SWNT hole polarons adjacent to (geminate) and nearby (nongeminate) PDI radical anions. Analysis of the free energy dependence for charge recombination (CR) of [SWNT •+ ] geminate -(PDI –• ) CS states highlights an ∼2 meV value for D–A electronic coupling (HAB) and ∼0.93 eV for the total reorganization energy (λ T ). A corresponding driving force dependence of the CR dynamics for [SWNT •+ ] nongeminate -(PDI –• ) CS states indicates a diminished H AB value (∼0.6 meV) and a larger λ T (∼1.1 eV), consistent with larger transfer distances. SWNT excitons that persist following photoinduced CS drive photooxidation of PDI –• components of [SWNT (•+)n ]-(PDI –• ) n CS states ( 1 SWNT* + PDI •– → PDI + SWNT •– ); this reaction manifests a significantly reduced λ T value (∼0.67 eV) as the initially prepared SWNT reduced state bears the character of a conduction band injected electron ([SWNT •– ] CB ). This reaction thus gives rise to relaxed, nongeminate SWNT electron and hole polarons on the same nanotube; these polarons react on a 10 2 ps time scale independent of the electronic structure of these SWNT-PDI superstructures.

anions↗

Averaged invariants in storage rings with synchrotron motion

In an ideal accelerator, the single-particle dynamics can be decoupled into transverse motion—the betatron oscillations— and longitudinal motion—the synchrotron oscillations. Chromatic and dispersive effects introduce a coupling between these dynamics, the so-called synchro-betatron coupling. We present an analysis of the fully coupled dynamics over a single synchrotron oscillation that leads to an averaged invariant with synchro-betatron coupling in a generic lattice. Here, we apply this analysis to two problems: first, a toy lattice where the computations are analytically tractable, then a design for a rapid cycling synchrotron built using the integrable optics described by Danilov and Nagaitsev, showing that although there is fairly complex behavior over the course of a synchrotron oscillation, the Danilov-Nagaitsev invariants are nevertheless periodic with the synchrotron motion.

43 PARTICLE ACCELERATORS↗

Dynamic Tuning of Moiré Excitons in a WSe 2 /WS 2 Heterostructure via Mechanical Deformation

Moiré superlattices in van der Waals (vdW) heterostructures form by stacking atomically thin layers on top of one another with a twist angle or lattice mismatch. The resulting moiré potential leads to a strong modification of the band structure, which can give rise to exotic quantum phenomena ranging from correlated insulators and superconductors to moiré excitons and Wigner crystals. Here, in this study, we demonstrate the dynamic tuning of moiré potential in a WSe 2 /WS 2 heterostructure at cryogenic temperature. We utilize the optical fiber tip of a cryogenic scanning near-field optical microscope (SNOM) to locally deform the heterostructure and measure its near-field optical response simultaneously. The deformation of the heterostructure increases the moiré potential, which leads to a red shift of the moiré exciton resonances. We observe the interlayer exciton resonance shifts up to 20 meV, while the intralayer exciton resonances shift up to 17 meV.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Nonperturbative and perturbative dynamics of a light QCD axion: Dark matter and the strong 𝐶⁢𝑃 problem

Considerable theoretical efforts have gone into expanding the reach of the quantum chromodynamics (QCD) axion beyond its canonical mass–decay-constant relation. The 𝑍 𝒩 QCD axion model reduces the QCD axion mass naturally, by invoking a discrete 𝑍 𝒩 symmetry through which the axion field is coupled to 𝒩 copies of the Standard Model. Before the QCD phase transition at temperature 𝑇 QCD , the 𝑍 𝒩 potential has a minimum at misalignment angle 𝜃 = 𝜋. At 𝑇 QCD , 𝜃 = 𝜋 becomes a maximum; the axion potential becomes exponentially suppressed and develops 𝒩 minima—only one of which actually solves the strong 𝐶⁢𝑃 problem. Before 𝑇 QCD , 𝜃 relaxes toward 𝜋. After 𝑇 QCD , the axion field starts from around the hilltop and may have sufficient kinetic energy to overcome the newly suppressed potential barriers. Such a field evolution leads to nonperturbative effects via the self-interactions near the hilltop, which can cause the exponential growth of fluctuations and backreaction on the coherent motion. This behavior can influence the relic density of the field and the minimum in which it settles. We conduct the first lattice simulations of the 𝑍 𝒩 QCD axion using 𝒞osmoℒattice to accurately calculate dark matter abundances and find nonperturbative dynamics reduce the abundance by up to a factor of two. We furthermore find that the probability of solving the strong 𝐶⁢𝑃 problem tends to diverge considerably from the naïve expectation of 1/𝒩.

Axions↗

Quantum simulation of nonequilibrium dynamics and thermalization in the Schwinger model

In this work, we present simulations of nonequilibrium dynamics of quantum field theories on digital quantum computers. As a representative example, we consider the Schwinger model, a (1+1)-dimensional U(1) gauge theory, coupled through a Yukawa-type interaction to a thermal environment described by a scalar field theory. We use the Hamiltonian formulation of the Schwinger model discretized on a spatial lattice. With the thermal scalar fields traced out, the Schwinger model can be treated as an open quantum system and its real-time dynamics are governed by a Lindblad equation in the Markovian limit. The interaction with the environment ultimately drives the system to thermal equilibrium. In the quantum Brownian motion limit, the Lindblad equation is related to a field theoretical Caldeira-Leggett equation. By using the Stinespring dilation theorem with ancillary qubits, we perform studies of both the nonequilibrium dynamics and the preparation of a thermal state in the Schwinger model using IBM's simulator and quantum devices. The real-time dynamics of field theories as open quantum systems and the thermal state preparation studied here are relevant for a variety of applications in nuclear and particle physics, quantum information and cosmology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Gauge Picture of Quantum Dynamics

Although local Hamiltonians exhibit local time dynamics, this locality is not explicit in the Schrödinger picture in the sense that the wavefunction amplitudes do not obey a local equation of motion. We show that geometric locality can be achieved explicitly in the equations of motion by "gauging" the global unitary invariance of quantum mechanics into a local gauge invariance. That is, expectation values $\langle \psi|A|\psi \rangle$ are invariant under a global unitary transformation acting on the wavefunction $|\psi\rangle \to U |\psi\rangle$ and operators $A \to U A U^\dagger$, and we show that it is possible to gauge this global invariance into a local gauge invariance. To do this, we replace the wavefunction with a collection of local wavefunctions $|\psi_J\rangle$, one for each patch of space $J$. The collection of spatial patches is chosen to cover the space; e.g. we could choose the patches to be single qubits or nearest-neighbor sites on a lattice. Local wavefunctions associated with neighboring pairs of spatial patches $I$ and $J$ are related to each other by dynamical unitary transformations $U_{IJ}$. The local wavefunctions are local in the sense that their dynamics are local. That is, the equations of motion for the local wavefunctions $|\psi_J\rangle$ and connections $U_{IJ}$ are explicitly local in space and only depend on nearby Hamiltonian terms. (The local wavefunctions are many-body wavefunctions and have the same Hilbert space dimension as the usual wavefunction.) We call this picture of quantum dynamics the gauge picture since it exhibits a local gauge invariance. The local dynamics of a single spatial patch is related to the interaction picture, where the interaction Hamiltonian consists of only nearby Hamiltonian terms. We can also generalize the explicit locality to include locality in local charge and energy densities.

Physics↗

Dataset_for_Molecular_Motion_Below_the_Glass_Transition_A_Solid-State_NMR_Study_of_Siloxane_Polymer_Dynamics Study

This dataset contains solid-state 1H and 13C NMR relaxometry data, differential scanning calorimetry (DSC) data, and size exclusion chromatography (SEC/GPC) data supporting the study of sub-glass-transition (sub-Tg) molecular dynamics in a composition- and sequence-controlled series of diphenyl-substituted polysiloxanes (PDMS, 14Ph, 33Ph, 50Ph, 67Ph, and 100Ph; 0–100% diphenylsiloxane content by mole).All solid-state NMR data were acquired on a 200 MHz Bruker Avance III HD spectrometer using a static 7 mm HX probe or a 4 mm HX probe under 4 kHz magic-angle spinning. Raw Bruker TopSpin experiment folders are included for: (1) variable-temperature 1H lineshape measurements used to determine linewidth (FWHM) as a function of temperature across the glass transition; (2) 1H T1 (saturation recovery with solid-echo detection), probing nanosecond-scale dynamics near the 1H Larmor frequency; (3) 1H T1rho (direct spin-lock, 62.5 kHz), probing microsecond-scale segmental dynamics; (4) 13C-detected Lee–Goldburg cross-polarization 1H T1rho (LGCPH T1rho) for 33Ph and 50Ph, resolving aromatic and aliphatic proton environments; and (5) 13C T1 relaxation for 33Ph and 50Ph. Differential scanning calorimetry data (TA Instruments DSC 25, −150 to +120 °C, up to +300 °C for 100Ph, 10 °C/min) are included for all six compositions and support the glass-transition temperatures in Table 1 and Figure 1. Size exclusion chromatography data (Agilent 1200 Series, PL-Gel 300 mixed-C column, THF mobile phase, polystyrene calibration standards) are included for the three synthesized copolymers (33Ph, 50Ph, 67Ph) and support the number-average molecular weights in Table 1. Processed data include per-composition relaxation-time summaries (Excel), curve-fitting and Bloembergen-Purcell-Pound (BPP) model analysis notebooks (Jupyter/Python), and Igor Pro (.pxp) master files used to generate the manuscript's figures.

Bloembergen-Purcell-Pound theory↗