Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space methods”

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 631 records · Page 35

Multi-Step Crystallization of Self-Organized Spiral Eutectics

A method for the solidification of metallic alloys involving spiral self-organization is presented as a new strategy for producing large-area chiral patterns with emergent structural and optical properties, with attention to the underlying mechanism and dynamics. This study reports the discovery of a new growth mode for metastable, two-phase spiral patterns from a liquid metal. Crystallization proceeds via a non-classical, two-step pathway consisting of the initial formation of a polytetrahedral seed crystal, followed by ordering of two solid phases that nucleate heterogeneously on the seed and grow in a strongly coupled fashion. Crystallographic defects within the seed provide a template for spiral self-organization. These observations demonstrate the ubiquity of defect-mediated growth in multi-phase materials and establish a pathway toward bottom-up synthesis of chiral materials with an inter-phase spacing comparable to the wavelength of infrared light. Given that liquids often possess polytetrahedral short-range order, our results are applicable to many systems undergoing multi-step crystallization.

36 MATERIALS SCIENCE↗

Quantum Thermodynamics of Nonequilibrium Processes in Lattice Gauge Theories

A key objective in nuclear and high-energy physics is to describe nonequilibrium dynamics of matter, e.g., in the early Universe and in particle colliders, starting from the standard model of particle physics. Classical computing methods, via the framework of lattice gauge theory, have experienced limited success in this mission. Quantum simulation of lattice gauge theories holds promise for overcoming computational limitations. Because of local constraints (Gauss’s laws), lattice gauge theories have an intricate Hilbert-space structure. This structure complicates the definition of thermodynamic properties of systems coupled to reservoirs during equilibrium and nonequilibrium processes. We show how to define thermodynamic quantities such as work and heat using strong-coupling thermodynamics, a framework that has recently burgeoned within the field of quantum thermodynamics. Our definitions suit instantaneous quenches, simple nonequilibrium processes undertaken in quantum simulators. To illustrate our framework, we compute the work and heat exchanged during a quench in a Z 2 lattice gauge theory coupled to matter in 1+1 dimensions. Here, the thermodynamic quantities, as functions of the quench parameter, evidence a phase transition. For general thermal states, we derive a simple relation between a quantum many-body system’s entanglement Hamiltonian, measurable with quantum-information-processing tools, and the Hamiltonian of mean force, used to define strong-coupling thermodynamic quantities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Entanglement in Gravity and Quantum Field Theory (Final Report)

It is becoming increasingly clear that ideas from quantum information theory, particularly the notion of quantum entanglement, play a fundamental role in some of the deepest aspects of our modern theories of quantum fields and gravity. The aim of this research was to explore the role that quantum entanglement plays in quantum field theories and in the nature of space-time and gravity. Building on a variety of earlier results obtained in these regards at the University of Illinois, we explored the constraints on the dynamical content of quantum field theories that follow from their entanglement properties. Topological field theories are important examples of particularly simple quantum field theories whose patterns of entanglement make connections between high energy physics, condensed matter physics and mathematics. These theories are directly relevant to low energy properties of certain materials. The study of such theories allowed us to investigate ideas that are relevant to quantum information research, such as new notions of entanglement between multiple parties and the quantum properties of interfaces between different phases of such materials. In addition, we employed new results in mathematics which strengthen monotonicity constraints on relative entropy to study their ramifications in quantum field theories, and we used quantum information methods to study the emergence of quantum gravity and string theory in holographic quantum field theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Guided search for desired functional responses via Bayesian optimization of generative model: Hysteresis loop shape engineering in ferroelectrics

Advances in theoretical modeling across multiple disciplines have yielded generative models capable of high veracity in predicting macroscopic functional responses of materials emerging as a result of complex non-local interactions. Correspondingly, of interest is the inverse problem of finding the model parameter that will yield desired macroscopic responses, such as stress–strain curves, ferroelectric hysteresis loops, etc. Here, we suggest and implement Gaussian process based methods that allow to effectively sample the degenerate parameter space of a complex non-local model to output regions of parameter space which yield desired functionalities. In this work, we discuss the specific adaptation of the acquisition function and sampling function to make the process efficient and balance the efficient exploration of parameter space for multiple possible minima and exploitation to densely sample the regions of interest where target behaviors are optimized. Additionally, this approach is illustrated via the hysteresis loop engineering in ferroelectric materials but can be adapted to other functionalities and generative models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quark and gluon two-loop beam functions for leading-jet p T and slicing at NNLO

We compute the complete set of two-loop beam functions for the transverse momentum distribution of the leading jet produced in association with an arbitrary colour-singlet system. Our results constitute the last missing ingredient for the calculation of the jet-vetoed cross section at small veto scales at the next-to-next-to-leading order, as well as an important ingredient for its resummation to next-to-next-to-next-to-leading logarithmic order. Our calculation is performed in the soft-collinear effective theory framework with a suitable regularisation of the rapidity divergences occurring in the phase-space integrals. We discuss the occurrence of soft-collinear mixing terms that might violate the factorisation theorem, and demonstrate that they are naturally absorbed into the beam functions at two loops in the exponential rapidity regularisation scheme when performing a multipole expansion of the measurement function. As in our recent computation of the two-loop soft function, we present the results as a Laurent expansion in the jet radius R. We provide analytic expressions for all flavour channels in x space with the exception of a set of R-independent non-logarithmic terms that are given as numerical grids. We also perform a fully numerical calculation with exact R dependence, and find that it agrees with our analytic expansion at the permyriad level or better. Our calculation allows us to define a next-to-next-to-leading order slicing method using the leading-jet p T as a slicing variable. As a check of our results, we carry out a calculation of the Higgs and Z boson total production cross sections at the next-to-next-to-leading order in QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theoretical search for possible Li–Ni–B crystal structures using an adaptive genetic algorithm

The structural diversity of rare-earth and transition metal borides indicates that alkali-transition metal borides (A-T-B) show tremendous promise in exhibiting a variety of crystal structures with different dimensionalities of T-B frameworks. On the other hand, the A-T-B ternary systems are severely underexplored because of the synthetic challenges associated with their preparation. Accurate and efficient computational predictions of low-energy stable and metastable phases can identify the optimal compositions of the hypothetical compounds in the A-T-B systems to guide the synthesis. As such, in this work, we have computationally discovered several new phases in the Li–Ni–B ternary system. The newly discovered LiNiB, Li 2 Ni 3 B, and Li 2 NiB phases expand the existing theoretical database, and the convex-hull surface of Li–Ni–B has been re-constructed. The lowest energy structure of the LiNiB compound has been found by an adaptive genetic algorithm with layered motif, which matches with the experimentally determined structure. According to our electrochemical calculations, LiNiB and another predicted layered Li 2 NiB compounds have great potential as anode materials for lithium batteries. The Li 2 Ni 3 B compound with the space group P4 3 32 was predicted to crystallize in a cubic structure composed of distorted octahedral units of BNi 6 , which is isostructural to two noncentrosymmetric superconductors Li 2 Pd 3 B and Li 2 Pt 3 B. While we were unable to experimentally confirm the Li 2 Ni 3 B compound utilizing the hydride synthetic route, attempts to synthesize this compound by alternate methods remain highly desirable, considering its potential superconducting properties.

36 MATERIALS SCIENCE↗

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis↗

Identifying crystallographic faces of the fluorites urania and thoria with rotational polarized Raman spectroscopy

Crystallite orientation identification is invaluable, but is often limited to small area identification or requires a large area sample. Nondestructive optical methods such as polarized Raman spectroscopy, in contrast, can be used to completely map a variety of sample sizes, but their potential is not yet fully realized. Here, we report a systematic study of polarized Raman scattering of high-quality, hydrothermally grown, single crystals of urania and thoria. The peak intensity variations for as-grown major crystal planes, post-growth polished crystal planes, and a post-growth polished non-crystallographic plane are directly linked to crystallographic orientation and crystal rotation, and agree with computed models. In particular, the parallel polarized peak intensity results are directly correlated with metal–oxygen–metal chains in the fluorite structure and can be used to determine both orientation and rotational alignment of a given crystal face if sufficiently small rotational steps are applied. These results are structure based, being applicable to the larger fluorite phase space, which is useful for optical, semiconductor, nuclear, and solid oxide fuel cell industries. Further, these results suggest that Raman spectroscopy can identify non-crystallographic orientations that are not discernable by traditional means.

36 MATERIALS SCIENCE↗

Scaling relations in the phase-space structure of dark matter haloes

We present new scaling relations for the isotropic phase-space distribution functions (DFs) and energy distributions of simulated dark matter haloes. These relations are inspired by those for the singular isothermal sphere with density profile $\rho (r)\propto r^{-2}$, for which the DF satisfies $f(E) \propto r_{\mathrm{max}}^{-2}(E)$ and the energy distribution satisfies $\mathrm{ d}M/\mathrm{ d}E \propto r_{\mathrm{max}}(E)$, with $r_{\mathrm{max}}(E)$ being the radius where the gravitational potential equals energy $E$. For the simulated haloes, we find $f(E)\propto r_{\mathrm{max}}^{-2.08}(E)$ and $\mathrm{ d}M/\mathrm{ d}E \propto r_{\mathrm{max}}(E)$ across broad energy ranges. In addition, the proportionality coefficients depend on the gravitational constant and the parameters of the best-fitting Navarro–Frenk–White density profile. These scaling relations are satisfied by haloes over a wide mass range and provide an efficient method to approximate their DFs and energy distributions. Understanding the origin of these relations may shed more light on halo formation.

galaxies: haloes↗

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗

Isochronous and period-doubling diagrams for symplectic maps of the plane

Symplectic mappings of the plane serve as key models for exploring the fundamental nature of complex behavior in nonlinear systems. Central to this exploration is the effective visualization of stability regimes, which enables the interpretation of how systems evolve under varying conditions. While the area-preserving quadratic Hénon map has received significant theoretical attention, a comprehensive description of its mixed parameter-space dynamics remain lacking. This limitation arises from early attempts to reduce the full two-dimensional phase space to a one-dimensional projection, a simplification that resulted in the loss of important dynamical features. Consequently, there is a clear need for a more thorough understanding of the underlying qualitative aspects. This paper aims to address this gap by revisiting the foundational concepts of reversibility and associated symmetries, first explored in the early works of G.D. Birkhoff. We extend the original framework proposed by Hénon by adding a period-doubling diagram to his isochronous diagram, which allows to represents the system’s bifurcations and the groups of symmetric periodic orbits that emerge in typical bifurcations of the fixed point. A qualitative and quantitative explanation of the main features of the region of parameters with bounded motion is provided, along with the application of this technique to other symplectic mappings, including cases of multiple reversibility. Modern chaos indicators, such as the Reversibility Error Method (REM) and the Generalized Alignment Index (GALI), are employed to distinguish between various dynamical regimes in the mixed space of variables and parameters. These tools prove effective in differentiating regular and chaotic dynamics, as well as in identifying twistless orbits and their associated bifurcations. Additionally, we discuss the application of these methods to real-world problems, such as visualizing dynamic aperture in accelerator physics, where our findings have direct relevance.

43 PARTICLE ACCELERATORS↗

Integrated control of edge localized modes and divertor flux using mixed toroidal harmonic resonant magnetic perturbations in EAST

Mixed harmonic resonant magnetic perturbations (RMPs) for integrated edge localized modes (ELMs) and divertor flux control are demonstrated in EAST target plasmas of low input torque and normalized beta β N ~1.7–1.9, which are close to the equivalent value in ITER high Q operation. The applied RMPs are designed to combine a static harmonic of the toroidal mode number n = 3 with a static or rotating harmonic of n = 2. ELM suppression is achieved without a drop of plasma energy confinement, and tungsten concentration is effectively reduced during the application of RMPs. With mixed harmonics, the toroidal varying steady state heat and particle fluxes on the divertor target can be modified with the rotating n = 2 harmonic, which agrees with the numerical modeling of three-dimensional magnetic topology, with plasma responses being taken into account. ELM suppression correlates with the times of larger n = 3 response with mixed n = 2 and n = 3 RMPs. The mixture of harmonics and the rotating n = 2 harmonic does not require additional coil current because the variation is only in the upper-lower coil current phase space. Furthermore, these results further affirm the effectiveness of integrated ELM and divertor flux control using RMPs with mixed harmonics and improve the understanding of the role of plasma responses in ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural Evolution and Photoluminescence Quenching across the FASnI 3– x Br x ( x = 0–3) Perovskites

One of the primary methods for band gap tuning in metal halide perovskites has been halide (I/Br) mixing. Despite widespread usage of this type of chemical substitution in perovskite photovoltaics, there is still little understanding of the structural impacts of halide alloying, with the assumption being the formation of ideal solid solutions. The FASnI 3–x Br x (x = 0–3) family of compounds provides the first example where the assumption breaks down, as the composition space is broken into two unique regimes (x = 0–2.9; x = 2.9–3) based on their average structure with the former having a 3D and the latter having an extended 3D (pseudo 0D) structure. Pair distribution function (PDF) analyses further suggest a dynamic 5s 2 lone pair expression resulting in increasing levels of off-centering of the central Sn as the Br concentration is increased. These antiferroelectric distortions indicate that even the x = 0–2.9 phase space behaves as a nonideal solid-solution on a more local scale. Solid-state NMR confirms the difference in local structure yielding greater insight into the chemical nature and local distributions of the FA + cation. In contrast to the FAPbI 3–x Br x series, a drastic photoluminescence (PL) quenching is observed with x ≥ 1.9 compounds having no observable PL. In conclusion, our detailed studies attribute this quenching to structural transitions induced by the distortions of the [SnBr 6 ] octahedra in response to stereochemically expressed lone pairs of electrons. This is confirmed through density functional theory, having a direct impact on the electronic structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating Covalent Bonding in f-elements using Gas-phase Ion Chemistry

Introduction In the reprocessing of f-elements present in used nuclear fuels, a variety of diglycolamides (DGA’s) are used as extractants for actinide partitioning. In particular, the Actinide-Lanthanide Separation (ALSEP) process typically utilizes either the N,N,N’,N’-tetraoctyl diglycolamide (TODGA) or N,N,N',N'-tetra-2-ethylhexyl diglycolamide (T2EHDGA) extractant ligands following the partitioning of uranium and plutonium from used nuclear fuel. To better understand fundamental interactions in these processes, covalent bonding of several f-elements with diglycolamides, primarily TODGA, is investigated in the gas phase using nanospray ionization and a quadrupole time-of-flight mass spectrometer. Further, analysis of the identity and relative strength of the cluster is enabled by MS2 isolation and collision induced dissociation. Methods Metal ion cluster analysis was completed using a Bruker mircOTOF-Q II quadrupole time-of-flight mass spectrometer equipped with a CaptiveSpray nanospray ion source. Metal:ligand solutions were prepared as 30 µM europium nitrate, samarium nitrate, cerium nitrate, or holmium nitrate and 3 µM DGA in acetonitrile or a 50:50 mixture of acetonitrile: isopropanol. Cluster mass spectra and collision-induced dissociation experiments were conducted in positive mode. Preliminary data To examine the patterns and relative strength of lanthanide: DGA interactions, MS2 experiments were completed with each lanthanide species listed above. Preliminary analyses of samarium and europium TODGA clusters suggest several combinations of TODGA and nitrate forming. The samarium cluster experiments yielded Sm(TODGA)x clusters with a samarium:TODGA ratio of up to 1:7 able to be isolated and evidence of greater ratios present in the mass spectrum. This is surprising, as metal clusters are not expected to have a coordination space able to accommodate this many ligands as large as TODGA. MS2 experiments show that, at higher ratios and with sufficient collision energy, entire TODGA ligands are removed instead of being fragmented. These experiments show that a lower collision energy is required to remove ligands as the number of bound TODGA’s increases, suggesting that in larger clusters, ligands are more delicately complexed to the metal. In addition to Sm(TODGA)x, several clusters were observed with nitrate ions bound to the metal in addition to TODGA. With a single nitrate ion, clusters with up to six TODGA’s were able to be isolated. In a similar pattern to the samarium clusters with only TODGA, less collision energy is required to eliminate one or more TODGA’s with increasing size. MS2 experiments suggest clusters with one nitrate appear to be of an equivalent or greater stability to clusters which replace the nitrate with a TODGA, as more collision energy is required to remove a TODGA ligand. These species with one nitrate are also in a higher abundance than the equivalent TODGA only cluster. With two nitrate ions, only clusters with a single TODGA were able to be isolated. Analogous europium experiments resulted in very similar clusters. Ratios of up to 1:7 Eu:TODGA were able to be isolated, and clusters with one nitrate and up to five TODGAs were isolated. In clusters with two nitrate ions, one or two TODGA’s could also be bound to the metal. MS2 experiments suggested, similarly to samarium, that larger clusters required less collision energy to eliminate TODGA. Europium clusters with one nitrate are in greater abundance and are stronger than the equivalent cluster which replaces the nitrate with TODGA. Similar analysis with cerium and holmium is ongoing, as well as analysis with other DGA ligands to compare relative strengths of the lanthanide metals with various extractant ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coherent correlation imaging for resolving fluctuating states of matter

Fluctuations and stochastic transitions are ubiquitous in nanometre-scale systems, especially in the presence of disorder. However, their direct observation has so far been impeded by a seemingly fundamental, signal-limited compromise between spatial and temporal resolution. Here we develop coherent correlation imaging (CCI) to overcome this dilemma. Our method begins by classifying recorded camera frames in Fourier space. Contrast and spatial resolution emerge by averaging selectively over same-state frames. Temporal resolution down to the acquisition time of a single frame arises independently from an exceptionally low misclassification rate, which we achieve by combining a correlation-based similarity metric with a modified, iterative hierarchical clustering algorithm. We apply CCI to study previously inaccessible magnetic fluctuations in a highly degenerate magnetic stripe domain state with nanometre-scale resolution. We uncover an intricate network of transitions between more than 30 discrete states. Our spatiotemporal data enable us to reconstruct the pinning energy landscape and to thereby explain the dynamics observed on a microscopic level. CCI massively expands the potential of emerging high-coherence X-ray sources and paves the way for addressing large fundamental questions such as the contribution of pinning and topology in phase transitions and the role of spin and charge order fluctuations in high-temperature superconductivity.

36 MATERIALS SCIENCE↗

Path Properties of Atmospheric Transitions: Illustration with a Low-Order Sudden Stratospheric Warming Model

Many rare weather events, including hurricanes, droughts, and floods, dramatically impact human life. To accurately forecast these events and characterize their climatology requires specialized mathematical techniques to fully leverage the limited data that are available. Here we describe transition path theory (TPT), a framework originally developed for molecular simulation, and argue that it is a useful paradigm for developing mechanistic understanding of rare climate events. TPT provides a method to calculate statistical properties of the paths into the event. As an initial demonstration of the utility of TPT, we analyze a low-order model of sudden stratospheric warming (SSW), a dramatic disturbance to the polar vortex that can induce extreme cold spells at the surface in the midlatitudes. SSW events pose a major challenge for seasonal weather prediction because of their rapid, complex onset and development. Climate models struggle to capture the long-term statistics of SSW, owing to their diversity and intermittent nature. We use a stochastically forced Holton–Mass-type model with two stable states, corresponding to radiative equilibrium and a vacillating SSW-like regime. In this stochastic bistable setting, from certain probabilistic forecasts TPT facilitates estimation of dominant transition pathways and return times of transitions. These “dynamical statistics” are obtained by solving partial differential equations in the model’s phase space. With future application to more complex models, TPT and its constituent quantities promise to improve the predictability of extreme weather events through both generation and principled evaluation of forecasts.

54 ENVIRONMENTAL SCIENCES↗

Phase control of nonlinear Breit-Wheeler pair creation

Electron-positron pair creation occurs throughout the universe in the environments of extreme astrophysical objects, such as pulsar magnetospheres and black hole accretion disks. The difficulty of emulating these environments in the laboratory has motivated the use of ultrahigh-intensity laser pulses for pair creation. Here we show that the phase offset between a laser pulse and its second harmonic can be used to control the relative transverse motion of electrons and positrons created in the nonlinear Breit-Wheeler process. Analytic theory and particle-in-cell simulations of a head-on collision between a two-color laser pulse and electron beam predict that with an appropriate phase offset, the electrons will drift in one direction and the positrons in the other. The resulting current may provide a collective signature of nonlinear Breit-Wheeler, while the spatial separation resulting from the relative motion may facilitate isolation of positrons for subsequent applications or detection. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗