Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical mechanics”

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 919 records · Page 51

Monopoles, exotic states and muon $$g-2$$ in TeV scale trinification

Abstract We study the low energy implications of a trinification model based on the gauge symmetry $$G= SU(3)_c \times SU(3)_L \times SU(3)_R$$ G = S U ( 3 ) c × S U ( 3 ) L × S U ( 3 ) R , without imposing gauge coupling unification. A minimal model requires two Higgs multiplets that reside in the bi-fundamental representation of G , and this is shown to be adequate for accommodating the Standard Model (SM) fermion masses and generate, via loop corrections and seesaw mechanism, suitable masses for the heavy neutral leptons as well as the observed SM neutrinos. We estimate a lower bound of around 15 TeV for the masses of the new down- type quarks that are required by the $$SU(3)_L \times SU(3)_R$$ S U ( 3 ) L × S U ( 3 ) R symmetry. We examine the resonant production at the LHC of the new gauge bosons, which leads to a lower bound of 16 TeV for the symmetry breaking scale of G . We also show how the muon $$g-2$$ g - 2 anomaly can be resolved in the presence of these new gauge bosons and the heavy charged leptons present in the model. Finally, the model predicts the presence of a topologically stable monopole carrying three quanta $$(6 \pi /e)$$ ( 6 π / e ) of Dirac magnetic charge and mass $$\gtrsim 160$$ ≳ 160 TeV. If new matter fields lying in the fundamental representations of G are included, the model predicts the presence of exotic leptons, mesons and baryons carrying fractional electric charges such as $$\pm e/3$$ ± e / 3 and $$\pm 2e/3$$ ± 2 e / 3 , fully compatible with the Dirac quantization condition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Robust finite-temperature many-body scarring on a quantum computer

Mechanisms for suppressing thermalization in disorder-free many-body systems, such as Hilbert space fragmentation and quantum many-body scars, have recently attracted much interest in foundations of quantum statistical physics and potential quantum information processing applications. However, their sensitivity to realistic effects such as finite temperature remains largely unexplored. Here, we have utilized IBM's Kolkata quantum processor to demonstrate an unexpected robustness of quantum many-body scars at finite temperatures when the system is prepared in a thermal Gibbs ensemble. We identify such robustness in the PXP model, which describes quantum many-body scars in experimental systems of Rydberg atom arrays and ultracold atoms in tilted Bose-Hubbard optical lattices. By contrast, other theoretical models which host exact quantum many-body scars are found to lack such robustness and their scarring properties quickly decay with temperature. Our study sheds light on the important differences between scarred models in terms of their algebraic structures, which impacts their resilience to finite temperature. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulation of decelerating streamers in inhomogeneous atmosphere with implications for runaway electron generation

The dynamics of positive and negative streamers is numerically simulated in atmospheric pressure air in the range of parameters corresponding to the streamer deceleration and termination in the middle of the discharge gap. A detailed comparison with experiments in air at constant and variable density demonstrates good agreement between the 2D simulation results and the observations. It is shown that positive and negative streamers behave in radically different ways when decelerating and stopping. When the head potential drops, the negative streamer transits to the mode in which the propagation is due to the forward electron drift. In this case, the radius of the ionization wave front increases, whereas the electric field at the streamer head decreases further and the streamer stops. Its head diameter continues to increase due to the slow drift of free electrons in the residual under-breakdown field. On the contrary, the only advancement mechanism for a positive streamer with a decreasing head potential is a decrease in the effective radius of the ionization wave, leading to a local increase in the electric field. This mechanism makes it possible to compensate for the reduction in the efficiency of gas photoionization at small head diameters. A qualitative 1D model is suggested to describe streamer deceleration and stopping for different discharge polarities. Estimates show that, during positive streamer stopping, the local electric field at the streamer head can exceed the threshold corresponding to the transition of electrons to the runaway mode when the head potential (relative to the surrounding space) decreases to~1.2 kV in atmospheric pressure air. In this case, pulsed generation of a beam of runaway electrons directed into the channel of a stopping positive streamer can occur. The energy of the formed pulsed electron beam depends on the intensity of photoionization in front of the streamer head. This energy can vary from 700 V (when increasing the photoionization rate by a factor of 10 with respect to the value in atmospheric pressure air) to 2.6 kV (when decreasing the photoionization rate by a factor of 1000). It is possible that this behavior of decelerating positive streamers can explain the observed bursts of x-ray radiation during the streamer propagation in long air gaps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

AENET–LAMMPS and AENET–TINKER : Interfaces for accurate and efficient molecular dynamics simulations with machine learning potentials

Machine-learning potentials (MLPs) trained on data from quantum-mechanics based first-principles methods can approach the accuracy of the reference method at a fraction of the computational cost. To facilitate efficient MLP-based molecular dynamics and Monte Carlo simulations, an integration of the MLPs with sampling software is needed. Here, we develop two interfaces that link the atomic energy network (ænet) MLP package with the popular sampling packages TINKER and LAMMPS. The three packages, ænet, TINKER, and LAMMPS, are free and open-source software that enable, in combination, accurate simulations of large and complex systems with low computational cost that scales linearly with the number of atoms. Scaling tests show that the parallel efficiency of the ænet–TINKER interface is nearly optimal but is limited to shared-memory systems. The ænet–LAMMPS interface achieves excellent parallel efficiency on highly parallel distributed memory systems and benefits from the highly optimized neighbor list implemented in LAMMPS. We demonstrate the utility of the two MLP interfaces for two relevant example applications: the investigation of diffusion phenomena in liquid water and the equilibration of nanostructured amorphous battery materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigation of the electrical conduction mechanisms in P-type amorphous germanium electrical contacts for germanium detectors in searching for rare-event physics

For the first time, electrical conduction mechanisms in the disordered material system is experimentally studied for p-type amorphous germanium (a-Ge) used for high-purity Ge detector contacts. The localization length and the hopping parameters in a-Ge are determined using the surface leakage current measured from three high-purity planar Ge detectors. The temperature dependent hopping distance and hopping energy are obtained for a-Ge fabricated as the electrical contact materials for high-purity Ge planar detectors. As a result, we find that the hopping energy in a-Ge increases as temperature increases while the hopping distance in a-Ge decreases as temperature increases. The localization length of a-Ge is on the order of 2.13$^{-0.05}_{+0.07}$A° to 5.07$^{0.83}_{+2.58}$A°,depending on the density of states near the Fermi energy level within bandgap. Using these parameters, we predict that the surface leakage current from a Ge detector with a-Ge contacts can be much smaller than one yocto amp (yA) at helium temperature, suitable for rare-event physics searches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Durability of plasma sprayed Thermal Barrier Coatings with controlled properties part Ⅱ: Effects of geometrical curvature

The present investigation elucidates the effects of substrate curvature on the durability of Air Plasma Sprayed (APS) Thermal Barrier Coatings (TBCs). Traditionally, planar disk specimens are utilized in Furnace Cycle Testing (FCT). In most cases, delamination is initiated at the disk's free edge and then propagates along the TBC-bond coat interface. However, in turbine components (e.g., blades/vanes), they lack significant free edges and the coatings are deposited on non-flat surfaces of varying radii of curvature. These geometrical discrepancies imply significant differences in the stress states as compared to planar disk specimens. Therefore, the part geometry inevitably affects the TBC failure mechanisms in the turbine system, and consequently influences TBC durability. Nevertheless, these effects have not been fully explored in the past, despite the anecdotal knowledge which suggests the leading edge of turbine blades to be one of the most common failure/spallation locations. Here in this study, representative TBC systems were deposited onto superalloy disks with flat substrate and rods with curved substrate. The experimental results from FCT suggest TBC durability on rods is significantly lower than when sprayed on disks. Furthermore, the durability of TBCs on rods appears to increase with higher porosity, which is contradictory to reported trends for TBCs on disks. In addition, the effects of various bond coats are not consistent with those observed in disks. To clarify the underlying effects of curved geometries, detailed stress analyses were performed. They revealed a state of thermal stresses which is unique to curved geometries. It was found that TBC on rods experience tensile radial stresses as well as tensile hoop stresses during the cooling. These stresses result in a complex interplay in failure processes in rods, which explains the different observed trends in the durability from those obtained with disks.

36 MATERIALS SCIENCE↗

X-ray linear dichroic ptychography

Biominerals such as seashells, coral skeletons, bone, and tooth enamel are optically anisotropic crystalline materials with unique nanoscale and microscale organization that translates into exceptional macroscopic mechanical properties, providing inspiration for engineering new and superior biomimetic structures. Using Seriatopora aculeata coral skeleton as a model, here, we experimentally demonstrate X-ray linear dichroic ptychography and map the c-axis orientations of the aragonite (CaCO 3 ) crystals. Linear dichroic phase imaging at the oxygen K-edge energy shows strong polarization-dependent contrast and reveals the presence of both narrow (<35°) and wide (>35°) c-axis angular spread in the coral samples. These X-ray ptychography results are corroborated by four-dimensional (4D) scanning transmission electron microscopy (STEM) on the same samples. Evidence of co-oriented, but disconnected, corallite subdomains indicates jagged crystal boundaries consistent with formation by amorphous nanoparticle attachment. We expect that the combination of X-ray linear dichroic ptychography and 4D STEM could be an important multimodal tool to study nano-crystallites, interfaces, nucleation, and mineral growth of optically anisotropic materials at multiple length scales.

4D scanning transmission electron microscopy↗

Optical current generation in graphene: CEP control vs. ω + 2 ω control

The injection of directional currents in solids with strong optical fields has attracted tremendous attention as a route to realize ultrafast electronics based on the quantum-mechanical nature of electrons at femto- to attosecond timescales. Such currents are usually the result of an asymmetric population distribution imprinted by the temporal symmetry of the driving field. Here we compare two experimental schemes that allow control over the amplitude and direction of light-field-driven currents excited in graphene. Both schemes rely on shaping the incident laser field with one parameter only: either the carrier-envelope phase (CEP) of a single laser pulse or the relative phase between pulses oscillating at angular frequencies ω and 2ω, both for comparable laser parameters. We observe that the efficiency in generating a current via two-color-control exceeds that of CEP control by more than two orders of magnitude (7 nA vs. 18 pA), as the ω + 2ω field exhibits significantly more asymmetry in its temporal shape. We support this finding with numerical simulations that clearly show that two-color current control in graphene is superior, even down to single-cycle pulse durations. We expect our results to be relevant to experimentally access fundamental properties of any solid at ultrafast timescales, as well as for the emerging field of petahertz electronics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evolution of 1 / f Flux Noise in Superconducting Qubits with Weak Magnetic Fields

The microscopic description of 1/f magnetic flux noise in superconducting circuits has remained an open question for several decades despite extensive experimental and theoretical investigation. Recent progress in superconducting devices for quantum information has highlighted the need to mitigate sources of qubit decoherence, driving a renewed interest in understanding the underlying noise mechanism(s). Though a consensus has emerged attributing flux noise to surface spins, their identity and interaction mechanisms remain unclear, prompting further study. Here, we apply weak in-plane magnetic fields to a capacitively shunted flux qubit (where the Zeeman splitting of surface spins lies below the device temperature) and study the flux-noise-limited qubit dephasing, revealing previously unexplored trends that may shed light on the dynamics behind the emergent 1/f noise. Notably, we observe an enhancement (suppression) of the spin-echo (Ramsey) pure-dephasing time in fields up to B = 100 G. With direct noise spectroscopy, we further observe a transition from a 1/f to approximately Lorentzian frequency dependence below 10 Hz and a reduction of the noise above 1 MHz with increasing magnetic field. Here, we suggest that these trends are qualitatively consistent with an increase of spin cluster sizes with magnetic field. These results should help to inform a complete microscopic theory of 1/f flux noise in superconducting circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗

Ferroelectric incommensurate spin crystals

Ferroics, especially ferromagnets, can form complextopological spin structures such as vortices and skyrmions when subjected to particular electrical and mechanical boundary conditions. Simple vortex-like, electric-dipole-based topological structures have been observed in dedicated ferroelectric systems, especially ferroelectric-insulator superlattices such as PbTiO 3 /SrTiO 3 , which was later shown to be a model system owing to its high depolarizing field. To date, the electric dipole equivalent of ordered magnetic spin lattices driven by the Dzyaloshinskii-Moriya interaction (DMi) has not been experimentally observed. Here we examine a domain structure in a single PbTiO 3 epitaxial layer sandwiched between SrRuO 3 electrodes. In this work, we observe periodic clockwise and anticlockwise ferroelectric vortices that are modulated by a second ordering along their toroidal core. The resulting topology, supported by calculations, is a labyrinth-like pattern with two orthogonal periodic modulations that form an incommensurate polar crystal that provides a ferroelectric analogue to the recently discovered incommensurate spin crystals in ferromagnetic materials. These findings further blur the border between emergent ferromagnetic and ferroelectrictopologies, clearing the way for experimental realization of further electric counterparts of magnetic DMi-driven phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fluid and hybrid simulations of the ionization instabilities in Hall thruster

Low-frequency axial oscillations in the range of 5–50 kHz stand out as a pervasive feature observed in many types of Hall thrusters. While it is widely recognized that the ionization effects play the central role in this mode, as manifested via the large-scale oscillations of neutral and plasma density, the exact mechanism(s) of the instabilities remain unclear. To gain further insight into the physics of the breathing mode and evaluate the role of kinetic effects, a one-dimensional time-dependent full nonlinear low-frequency model describing neutral atoms, ions, and electrons is developed in full fluid formulation and compared to the hybrid model in which the ions and neutrals are kinetic. Both models are quasi-neutral and share the same electron fluid equations that include the electron diffusion, mobility across the magnetic field, and the electron energy evolution. The ionization models are also similar in both approaches. Further, the predictions of fluid and hybrid simulations are compared for different test cases. Two main regimes are identified in both models: one with pure low-frequency behavior and the other one, where the low-frequency oscillations coexist with high-frequency oscillations in the range of 100–200 kHz, with the characteristic time scale of the ion channel fly-by time, 100–200 kHz. The other test case demonstrates the effect of a finite temperature of injected neutral atoms, which has a substantial suppression effect on the oscillation amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigating boosted decision trees as a guide for inertial confinement fusion design

Inertial confined fusion experiments at the National Ignition Facility have recently entered a new regime approaching ignition. Improved modeling and exploration of the experimental parameter space were essential to deepening our understanding of the mechanisms that degrade and amplify the neutron yield. The growing prevalence of machine learning in fusion studies opens a new avenue for investigation. Here in this paper, we have applied the Gradient-Boosted Decision Tree machine-learning architecture to further explore the parameter space and find correlations with the neutron yield, a key performance indicator. We find reasonable agreement between the measured and predicted yield, with a mean absolute percentage error on a randomly assigned test set of 35.5%. This model finds the characteristics of the laser pulse to be the most influential in prediction, as well as the hohlraum laser entrance hole diameter and an enhanced capsule fabrication technique. We used the trained model to scan over the design space of experiments from three different campaigns to evaluate the potential of this technique to provide design changes that could improve the resulting neutron yield. While these data-driven model cannot predict ignition without examples of ignited shots in the training set, it can be used to indicate that an unseen shot design will at least be in the upper range of previously observed neutron yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nuclear–electronic orbital approach to quantization of protons in periodic electronic structure calculations

The nuclear–electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born–Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems, such as two-dimensional materials and liquid–semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide–water interface. The zero-point energy effects of the protons as well as electron–proton correlation are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Locally commensurate charge-density wave with three-unit-cell periodicity in YBa 2 Cu 3 O y

In order to identify the mechanism responsible for the formation of charge-density waves (CDW) in cuprate superconductors, it is important to understand which aspects of the CDW’s microscopic structure are generic and which are material-dependent. Here, we show that, at the local scale probed by NMR, long-range CDW order in YBa 2 Cu 3 O y is unidirectional with a commensurate period of three unit cells ( λ = 3 b ), implying that the incommensurability found in X-ray scattering is ensured by phase slips (discommensurations). Furthermore, NMR spectra reveal a predominant oxygen character of the CDW with an out-of-phase relationship between certain lattice sites but no specific signature of a secondary CDW with λ = 6 b associated with a putative pair-density wave. These results shed light on universal aspects of the cuprate CDW. In particular, its spatial profile appears to generically result from the interplay between an incommensurate tendency at long length scales, possibly related to properties of the Fermi surface, and local commensuration effects, due to electron-electron interactions or lock-in to the lattice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Beyond quantum cluster theories: multiscale approaches for strongly correlated systems

The degrees of freedom that confer to strongly correlated systems their many intriguing properties also render them fairly intractable through typical perturbative treatments. For this reason, the mechanisms responsible for their technologically promising properties remain mostly elusive. Computational approaches have played a major role in efforts to fill this void. In particular, dynamical mean field theory and its cluster extension, the dynamical cluster approximation have allowed significant progress. However, despite all the insightful results of these embedding schemes, computational constraints, such as the minus sign problem in quantum Monte Carlo (QMC), and the exponential growth of the Hilbert space in exact diagonalization (ED) methods, still limit the length scale within which correlations can be treated exactly in the formalism. A recent advance aiming to overcome these difficulties is the development of multiscale many body approaches whereby this challenge is addressed by introducing an intermediate length scale between the short length scale where correlations are treated exactly using a cluster solver such QMC or ED, and the long length scale where correlations are treated in a mean field manner. At this intermediate length scale correlations can be treated perturbatively. This is the essence of multiscale many-body methods. Furthermore, we will review various implementations of these multiscale many-body approaches, the results they have produced, and the outstanding challenges that should be addressed for further advances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗