Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “collective 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 397 records · Page 22

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Construction and characterization of a genome-scale ordered mutant collection of Bacteroides thetaiotaomicron

Ordered transposon-insertion collections, in which specific transposon-insertion mutants are stored as monocultures in a genome-scale collection, represent a promising tool for genetic dissection of human gut microbiota members. However, publicly available collections are scarce and the construction methodology remains in early stages of development. Here, we describe the assembly of a genome-scale ordered collection of transposon-insertion mutants in the model gut anaerobe Bacteroides thetaiotaomicron VPI-5482 that we created as a resource for the research community. We used flow cytometry to sort single cells from a pooled library, located mutants within this initial progenitor collection by applying a pooling strategy with barcode sequencing, and re-arrayed specific mutants to create a condensed collection with single-insertion strains covering >2500 genes. To demonstrate the potential of the condensed collection for phenotypic screening, we analyzed growth dynamics and cell morphology. We identified both growth defects and altered cell shape in mutants disrupting sphingolipid synthesis and thiamine scavenging. Finally, we analyzed the process of assembling the B. theta condensed collection to identify inefficiencies that limited coverage. We demonstrate as part of this analysis that the process of assembling an ordered collection can be accurately modeled using barcode sequencing data. We expect that utilization of this ordered collection will accelerate research into B. theta physiology and that lessons learned while assembling the collection will inform future efforts to assemble ordered mutant collections for an increasing number of gut microbiota members.

59 BASIC BIOLOGICAL SCIENCES↗

Ice dynamics at the mouth of ice stream B, Antarctica

Data collected in the region of the mouth of ice stream B, West Antarctica, during three field seasons are presented. The physical characteristics of the mouth of ice stream B are described, and the dynamics in the vicinity of the DNB network are discussed. The dynamics of ice stream B from DNB to the grounding line is briefly considered, and a force analysis of the grounding line region is made. The results demonstrate that the dynamic situation of the region at the mouth of ice stream B is distinctly different from either the greater portion of the ice stream upstream or the Ross ice shelf downstream.

Bindschadler, R. A.↗

Capture and Decomposition of the Nerve Agent Simulant, DMCP, Using the Zeolitic Imidazolate Framework (ZIF-8)

Understanding mechanisms of decontamination of chemical warfare agents (CWA) is an area of intense research aimed at developing new filtration materials to protect soldiers and civilians in case of state-sponsored or terrorist attack. In this study, we employed complementary structural, chemical, and dynamic probes and in situ data collection, to elucidate the complex chemistry, capture, and decomposition of the CWA simulant, dimethyl chlorophosphonate (DMCP). Our work reveals key details of the reactive adsorption of DMCP and demonstrates the versatility of zeolitic imidazolate framework (ZIF-8) as a plausible material for CWA capture and decomposition. Here, the in situ synchrotron-based powder X-ray diffraction (PXRD) and pair distribution function (PDF) studies, combined with Fourier transform infrared (FT-IR) spectroscopy, X-ray photoelectron spectroscopy (XPS), zinc K-edge X-ray absorption near edge structure (XANES), and Raman spectroscopies, showed that the unique structure, chemical state, and topology of ZIF-8 enable accessibility, adsorption, and hydrolysis of DMCP into the pores and revealed the importance of linker chemistry and Zn 2+ sites for nerve agent decomposition. DMCP decontamination and decomposition product(s) formation were observed by thermogravimetric analysis, FT-IR spectroscopy, and phosphorus (P) K-edge XANES studies. Differential PDF analysis indicated that the average structure of ZIF-8 (at the 30 Å scale) remains unchanged after DMCP dosing and provided information on the dynamics of interactions of DMCP with the ZIF-8 framework. Using in situ PXRD and diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), we showed that nearly 90% regeneration of the ZIF-8 structure and complete liberation of DMCP and decomposition products occur upon heating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Statistical analysis of correlations in the x-ray induced Coulomb explosion of iodopyridine

Coulomb explosion imaging is a promising experimental tool to study individual molecules. Here, in this work simulation data for the x-ray induced Coulomb explosion of 2-iodopyridine is analyzed and the involved fragmentation dynamics are described. It is found that particular final ion momenta show correlations that reflect a collision of two atoms during the explosion. Variations of the proton kinetic energies can be associated with variations in the charge build-up time. Covariances of the forces during the explosion can be utilized to simplify the description of the dynamics in reduced dimensionality using only four collective coordinates.

74 ATOMIC AND MOLECULAR PHYSICS↗

Coordinating Multi-Rover Systems: Evaluation Functions for Dynamic and Noisy Environments

This paper addresses the evolution of control strategies for a collective: a set of entities that collectively strives to maximize a global evaluation function that rates the performance of the full system. Directly addressing such problems by having a population of collectives and applying the evolutionary algorithm to that population is appealing, but the search space is prohibitively large in most cases. Instead, we focus on evolving control policies for each member of the collective. The fundamental issue in this approach is how to create an evaluation function for each member of the collective that is both aligned with the global evaluation function and is sensitive to the fitness changes of the member, while relatively insensitive to the fitness changes of other members. We show how to construct evaluation functions in dynamic, noisy and communication-limited collective environments. On a rover coordination problem, a control policy evolved using aligned and member-sensitive evaluations outperfoms global evaluation methods by up to 400%. More notably, in the presence of a larger number of rovers or rovers with noisy and communication limited sensors, the proposed method outperforms global evaluation by a higher percentage than in noise-free conditions with a small number of rovers.

Turner, Kagan↗

Randomized Algorithms for Scientific Computing (RASC)

Randomized algorithms have propelled advances in artificial intelligence (AI) and represent a foundational research area in advancing AI for Science. Future advancements in DOE Office of Science priority areas such as climate science, astrophysics, fusion, advanced materials, combustion, and quantum computing all require randomized algorithms for surmounting challenges of complexity, robustness, and scalability. Advances in data collection and numerical simulation have changed the dynamics of scientific research and motivate the need for randomized algorithms. For instance, advances in imaging technologies such as X-ray ptychography, electron microscopy, electron energy loss spectroscopy, or adaptive optics lattice light-sheet microscopy collect hyperspectral imaging and scattering data in terabytes, at breakneck speed enabled by state-of-the-art detectors. The data collection is exceptionally fast compared with its analysis. Likewise, advances in high-performance architectures have made exascale computing a reality and changed the economies of scientific computing in the process. Floating-point operations that create data are essentially free in comparison with data movement. Thus far, most approaches have focused on creating faster hardware. Ironically, this faster hardware has exacerbated the problem by making data still easier to create. Under such an onslaught, scientists often resort to heuristic deterministic sampling schemes (e.g., low-precision arithmetic, sampling every nth element) and sacrifice potentially valuable accuracy. Dramatically better results can be achieved via randomized algorithms, reducing the data size as much as or more than naive deterministic subsampling can achieve, while retaining the high accuracy of computing on the full data set. By randomized algorithms we mean those algorithms that employ some form of randomness in internal algorithmic decisions to accelerate time to solution, increase scalability, or improve reliability. Examples include matrix sketching for solving large-scale least-squares problems (see Figure 1) and stochastic gradient descent for training machine learning models. We are not recommending heuristic methods but rather randomized algorithms that have certificates of correctness and probabilistic guarantees of optimality and near-optimality. Such approaches can be useful beyond acceleration, for example, in understanding how to avoid measure zero worst-case scenarios that plague methods such as QR matrix factorization.

97 MATHEMATICS AND COMPUTING↗

Relaxation and Aging of Nanosphere Assemblies at a Water–Oil Interface

The relaxation and aging of an assembly of spherical nanoparticles (NPs) at a water-oil interface are characterized in situ by grazing incidence X-ray photon correlation spectroscopy. The dynamics of the interfacial assembly is measured while the interface saturates with NPs. Weak attractions between NPs lead to gel-like structures in the assembly, where the in-plane ordering is inhibited by the broad size distribution of the NPs. Structural rearrangements on the length scale of the NP-NP center-to-center distances proceed by intermittent fluctuations instead of continuous cooperative motions. Here, the coexistence of rapid and slow NP populations is confirmed, as commonly observed in soft glass-forming materials. Dynamics are increasingly slowed as the NPs initially segregate to the locally clustered interface. The structural relaxation of the NPs in these localized clusters is 5 orders of magnitude slower than that of free particles in the bulk. When the interface is nearly saturated, the time for relaxation increases suddenly due to the onset of local jamming, and the dynamics slow exponentially afterward until the system reaches collective jamming by cooperative rearrangements. This investigation provides insights into structural relaxations near the glass transition and the evolution of the structure and dynamics of the assemblies as they transition from an isotropic liquid to a dense disordered film.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-consistent calculations of the electric charge, ion drag force, and the drift velocity of spherical grains using Langevin dynamics and comparisons against canonical experiments

We present trajectory simulation-based modeling to capture the interactions between ions and charged grains in dusty or complex plasmas. Our study is motivated by the need for a self-consistent and experimentally validated approach for accurately calculating the ion drag force and grain charge that determine grain collective behavior in plasmas. We implement Langevin dynamics in a computationally efficient predictor–corrector approach to capture multiscale ion and grain dynamics. Predictions of grain velocity, grain charge, and ion drag force are compared with prior measurements to assess our approach. The comparisons reveal excellent agreement to within ±20% between predicted and measured grain velocities [Yaroshenko et al., Phys. Plasmas 12, 093503 (2005) and Khrapak et al., Europhys. Lett. 97, 35001 (2012)] for 0.64, 1.25 μm grains at ∼20−500 Pa. Comparisons with the measured grain charge [Khrapak et al., Phys. Rev. E 72, 016406 (2005)] under similar conditions reveal agreement to within ∼20% as well. Measurements of the ion drag force [Hirt et al., Phys. Plasmas 11, 5690 (2004); IEEE Trans. Plasma Sci. 32, 582 (2004)] are used to assess the viability of the presented approach to calculate the ion drag force experienced by grains exposed to ion beams of well-defined energy. Excellent agreement between calculations and measurements is obtained for beam energies >10 eV, and the overprediction below 10 eV is attributed to the neglect of charge exchange collisions in our modeling. Along with critical assessments of our approach, suggestions for future experimental design to probe charging of and momentum transfer onto grains that capture the effect of space charge concentration and external fields are outlined.

Physics↗

Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of 48 Ca

Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. In this work, we use our scheme to compute the matrix element that governs the neutrinoless double beta decay of 48 Ca to 48 Ti, and find it to have the value 0.61, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as 76 Ge, 130 Te, and 136 Xe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Orientational Pair Correlations and Local Structure of Benzonitrile from Molecular Dynamics Simulations with Comparisons to Experiments [plus supplemental information]

Here, we present an experimentally parametrized molecular dynamics study of single-molecule and collective orientational relaxation in neat benzonitrile through the analysis of the reorientational anisotropy and polarizability anisotropy time-correlation function (PA-TCF). The simulations show that the PA-TCF is dominated by collective reorientation after 20 ps. Collective reorientation is found to be slower than single-molecule reorientation by a factor of 1.67, consistent with recent experiments. The simulations provide direct evidence of local antiparallel benzonitrile configurations. These structures, which have been the center of some debate, are responsible for the slower rate of collective vs. single-molecule reorientation in the liquid. Further structural analysis indicates that significant Coulombic interactions between the nitrile group and hydrogen atoms on adjacent molecules play a role in the formation of the antiparallel structures. The single-molecule dynamics reflected in the anisotropy are complex and consist of a ballistic regime, restricted angular diffusion, and spatially anisotropic free diffusion. The principal components of the rotational diffusion tensor are independently obtained and shown to reproduce the free diffusion regime of the anisotropy for each principal axis according to the predictions of a previous theory.

14 SOLAR ENERGY↗

New Perspectives on Vibrational Energy Transfer in Energetic Materials: Insights from Pressure-Tuned Ultrafast Spectroscopy

Understanding the manner in which vibrational energy flows between molecular and lattice vibrations is of great interest in physical chemistry due to its central role in reactivity and energy dissipation in molecular materials. Here, in this feature article, we highlight our recent efforts employing ultrafast broadband infrared spectroscopy toward understanding the interplay between molecular and lattice vibrations in energetic materials, motivated by the open questions surrounding the role of vibrational energy transfer (VET) in reaction initiation in these materials. Our work addresses the ongoing debate on the participation of doorway modes in VET. We further present new results from high-pressure ultrafast experiments on RDX, a hydrogen-bonded material, and BNFF, a hydrogen-free material, to explore how intermolecular interaction strength governs VET pathways and time scales. Collectively, our findings reveal that vibrational dynamics in these systems occurs across three distinct time regimes, with VET being incomplete out to hundreds of picoseconds, suggesting the importance of considering nonstatistical reactions in the modeling of these materials. These time scales vary as intermolecular interaction strength is indirectly modified by application of static pressure, indicating dramatic changes to the vibrational structure of these materials under shock-relevant conditions. Our results thus shed light on how intermolecular interactions shape vibrational energy redistribution in molecular materials, and highlight the need for further theoretical and experimental investigation.

Carlson, Daniel Ryan [Sandia National Laboratories↗

Gas puff imaging on the TCV tokamak

We present the design and operation of a suite of Gas Puff Imaging (GPI) diagnostic systems installed on the Tokamak à Configuration Variable (TCV) for the study of turbulence in the plasma edge and Scrape-Off-Layer (SOL). These systems provide the unique ability to simultaneously collect poloidal 2D images of plasma dynamics at the outboard midplane, around the X-point, in both the High-Field Side (HFS) and Low-Field Side (LFS) SOL, and in the divertor region. We describe and characterize an innovative control system for deuterium and helium gas injection, which is becoming the default standard for the other gas injections at TCV. Extensive pre-design studies and the different detection systems are presented, including an array of avalanche photodiodes and a high-speed CMOS camera. First results with spatial and time resolutions of up to ≈2 mm and 0.5 µs, respectively, are described, and future upgrades of the GPI diagnostics for TCV are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Demonstration of a two-dimensional divertor Thomson scattering system on DIII-D

Recent upgrades providing two-dimensional divertor Thomson scattering (DTS-2D) measurements of T e and n e during a DIII-D plasma shot and a thorough description of system components and their functionality are presented. This system expands the capabilities of the existing single divertor Floor measurement location by introducing seven additional laser beam path options in the poloidal plane, spanning major radii from 1.062 to 1.335 m. The system redirects ~1 J, 50 Hz Nd:YAG laser pulses to the new beam paths within 20 ms, stepping through each path on the divertor Floor every 200 ms during a plasma shot. The laser is redirected using an ex-vessel, fast-steering mirror to one of eight in-vessel beam paths oriented underneath the vessel tiles. Up to eleven measurement positions per beam path, from –1.35 to –1.13 m below the machine midplane, are available by dynamically refocusing the ex-vessel collection fiber array using a high-speed linear stage. As a result, current measurement positions above the divertor Shelf are retained via a hole in the fixed, in-vessel mirror, allowing laser pass through.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of new and existing methods for time-domain simulations of turbulence and loads

We seek to obtain a second-by-second match between the simulated and measured structural loads of a utility-scale wind turbine. To obtain the one-to-one load simulations, we start with the furthest upstream component of the modeling chain: the turbulent inflow. We consider new and existing methods to generate constrained-turbulence flow fields. The new method is based on large-eddy simulations (LES) and machine learning (ML). The existing methods include Kaimal-based TurbSim and the superstatistical wind field model. The inflow measurements used to constrain these simulations are obtained with a nacelle-mounted scanning lidar. We compare the flow fields for the different inflow simulation approaches and validate their associated load predictions against measurements collected in the Rotor Aero-dynamics, Aeroelastics, and Wake (RAAW) field campaign. We find that the rotor-position control developed for this study is key in enabling the time match between measurements and simulations. When this control approach is used, the load simulation performance tracks with the inflow simulation fidelity, with LES+ML yielding errors ≤ 4% for the damage-equivalent loads of flapwise bending moment, and tower fore-aft bending moments.

17 WIND ENERGY↗

Stability of trapped solutions of a nonlinear Schrödinger equation with a nonlocal nonlinear self-interaction potential

This work focuses on the study of the stability of trapped soliton-like solutions of a (1 + 1)-dimensional nonlinear Schrödinger equation (NLSE) in a nonlocal, nonlinear, self-interaction potential of the form [|Ψ(x,t)| 2 +|Ψ(-x,t)| 2 ] κ where κ is an arbitrary nonlinearity parameter. Although the system with κ = 1 (i.e. fully integrable case) was first reported by Yang (2018 Phys. Rev. E 98 042202), here in the present work, we extend this model to the one in which κ is arbitrary. This allows us to compare the stability properties of the now trapped solutions to previously found solutions of the more usual NLSE with κ ≠ 1 which are moving soliton solutions. We show that there is a simple, one-component, nonlocal Lagrangian and corresponding action governing the dynamics of the system. Using a collective coordinate method derived from the action as well as assuming the validity of Derrick's theorem, we find that these trapped solutions are stable for 0 < κ < 2 and unstable when κ > 2. At the critical value of κ, i.e. κ = 2, the solution can either collapse or blowup linearly in time when q 0 = 0, where q 0 is the center of the initial density ρ(x, t = 0) = ψ*ψ of the solution. For q 0 ≠ 0 the displaced solution collapses. When κ > 2 initial small displacements from the origin also lead to collapse of the wave function. This phenomenon is not seen in the usual NLSE.

collective coordinates↗

Damped Dirac magnon in the metallic kagome antiferromagnet FeSn

The kagome lattice is a fertile platform to explore topological excitations with both Fermi-Dirac and Bose-Einstein statistics. While relativistic Dirac fermions and flat bands have been discovered in the electronic structure of kagome metals, the spin excitations have received less attention. Here, we report inelastic neutron scattering studies of the prototypical kagome magnetic metal FeSn. The spectra display well-defined spin waves extending to 120 meV. Above this energy, the spin waves become progressively broadened, reflecting interactions with the Stoner continuum. Using linear spin-wave theory, we determine an effective spin Hamiltonian that reproduces the measured dispersion. This analysis indicates that the Dirac magnon at the K point remarkably occurs on the brink of a region where well-defined spin waves become unobservable. Furthermore, our results emphasize the influential role of itinerant carriers on the topological spin excitations of metallic kagome magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗