Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “invariance”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

New insights into Mn 2+ and Mg 2+ inhibition of calcite growth

Impurity ion and isotope partitioning into carbonate minerals provide a window into the molecular processes occurring at the fluid-mineral interface during crystal growth. Here, we employ calcium isotope fractionation together with process-based modeling to elucidate the mechanisms by which two divalent cations with starkly contrasting compatibility, magnesium and manganese, inhibit calcite growth and incorporate into the mineral lattice. Calcite growth inhibition by Mg 2+ is log-linear and K Mg is on the order of 0.02–0.03 throughout the range of {Mg 2+ }/{Ca 2+ } studied here (0.01–2.6). Mn 2+ exhibits much stronger log-linear growth rate inhibition at low Mn 2+ concentrations (fluid {Mn 2+ }/{Ca 2+ } = 0.001–0.02). Mn 2+ is readily incorporated into the calcite lattice to form a calcite-rhodochrosite solid solution, with large partition coefficients (K Mn 4.6–15.6) inversely correlated to growth rate. For both Mn 2+ and Mg 2+ , calcium isotope fractionation is found to be invariant with {Me 2+ }/{Ca 2+ } despite more than an order of magnitude decline in growth rate. This invariant Δ 44/40 Ca suggests that the presence of Mn 2 + or Mg 2+ does not significantly change the relative rates of Ca 2+ attachment and detachment at kink sites during growth, indicative of a dominantly kink blocking inhibition mechanism. Because the partitioning behavior dictates that Mn 2+ must attach to the surface significantly faster than Ca 2+ , attachment of Mn 2+ is likely to be as a non-monomer species such as an ion pair or possibly a larger polynuclear cluster. We propose that calcite growth rate inhibition by Mn is determined by the kinetics of carbonate attachment at Mn-occupied kink sites, potentially due to slow re-orientation kinetics of carbonate ions that have formed an inner-sphere complex with Mn 2+ at the surface but must reorient to incorporate into the lattice. We demonstrate that patterns in Mg 2+ partitioning and inhibition behavior are broadly consistent with growth inhibition driven by slow Mg 2+ -aquo complex dehydration relative to Ca 2+ but argue that this mechanism likely represents one endmember scenario, seen in Mg-calcite growth at low supersaturations and net precipitation rates. During growth at faster net precipitation rates, some portion of Mg 2+ is likely incorporated as a partially hydrated or otherwise complexed species, but calcite growth remains significantly inhibited by the kinetics of CO 3 2- attachment at Mg 2+ kink sites. These findings suggest a hybrid classical/nonclassical growth mechanism whereby Ca 2+ incorporates largely as a free ion at kink sites while Mn 2+ and some portion of Mg 2+ are incorporated via non-monomer attachment. This pattern may be generalizable; trace constituent cations with aquo-complex desolvation rates significantly slower than the mineral growth rate preferentially incorporate as a non-monomer species during otherwise classical crystal growth.

58 GEOSCIENCES↗
DOE OSTI2021

A consistent and conservative model and its scheme for N -phase- M -component incompressible flows

Here, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Bicrystallography-informed Frenkel–Kontorova model for interlayer dislocations in strained 2D heterostructures

In recent years, van der Waals (vdW) heterostructures and homostructures, which consist of stacks of two-dimensional (2D) materials, have risen to prominence due to their association with exotic quantum phenomena originating from correlated electronic states harbored by them. Atomistic scale relaxation effects play an extremely important role in the electronic scale quantum physics of these systems, providing means of manipulation of these materials and allowing them to be tailored for emergent technologies. We investigate such structural relaxation effects in this work using atomistic and mesoscale models, within the context of twisted bilayer graphene — a well-known heterostructure system that features moiré patterns arising from the lattices of the two graphene layers. For small twist angles, atomic relaxation effects in this system are associated with the natural emergence of interface dislocations or strain solitons, which result from the cyclic nature of the generalized stacking fault energy (GSFE), that measures the interface energy based on the relative movement of the two layers. Here, in this work, we first demonstrate using atomistic simulations that atomic reconstruction in bilayer graphene under a large twist also results from interface dislocations, although the Burgers vectors of such dislocations are considerably smaller than those observed in small-twist systems. To reveal the translational invariance of the heterointerface responsible for the formation of such dislocations, we derive the translational symmetry of the GSFE of a 2D heterostructure using the notions of coincident site lattices (CSLs) and displacement shift complete lattices (DSCLs). The workhorse for this exercise is a recently developed Smith normal form bicrystallography framework. Next, we construct a bicrystallography-informed and frame-invariant Frenkel–Kontorova model, which can predict the formation of strain solitons in arbitrary 2D heterostructures, and apply it to study a heterostrained, large-twist bilayer graphene system. Our mesoscale model is found to produce results consistent with atomistic simulations. We anticipate that the model will be invaluable in predicting structural relaxation and for providing insights into various heterostructure systems, especially in cases where the fundamental unit cell is large and therefore, atomistic simulations are computationally expensive.

2D heterostructures↗
DOE OSTI2020

Higgs boson pair production via gluon fusion at N 3 LO in QCD

We present next-to-next-to-next-to-leading order ( N 3 LO ) QCD predictions for the Higgs boson pair production via gluon fusion at hadron colliders in the infinite top-quark mass limit. Besides the inclusive total cross sections at various collision energies, we also provide the invariant mass distribution of the Higgs boson pair. Our results show that the N 3 LO QCD corrections enhance the next-to-next-to-leading order cross section by 3.0% (2.7%) at s = 13 ( 100 ) TeV , while the scale uncertainty is reduced substantially below 3% (2%). We also find that a judicious scale choice can significantly improve the perturbative convergence. For the invariant mass distribution, our calculation demonstrates that the N 3 LO corrections improve the scale dependence but almost do not change the shape.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Conundrum for the free energy of a holonomous gluonic plasma at cubic order

We compute the term ~ g 3 in the free energy for a S U ( N ) gauge theory with nonzero holonomy at nonzero temperature. If the holonomy is generated kinematically by the introduction of gauge invariant sources coupled to Polyakov loops, the contribution of charged (off-diagonal) gluons to the free energy at order g 3 , F ( 3 : c . g . ) , is singular: F ( 3 : c . g . ) ≠ 0 without holonomy, but F ( 3 : c . g . ) = 0 when the holonomy is nonzero, even infinitesimally. We show that the absence of the charged gluon contribution is required by gauge invariance alone and is therefore a universal feature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Evidence for WW/WZ vector boson scattering in the decay channel ℓνqq produced in association with two jets in proton-proton collisions at s = 13 TeV

Evidence is reported for electroweak (EW) vector boson scattering in the decay channel $\ell\nu$qq of two weak vector bosons WV (V = W or Z), produced in association with two parton jets. The search uses a data set of proton-proton collisions at 13 TeV collected with the CMS detector during 2016-2018 with an integrated luminosity of 138 fb$^{-1}$. Events are selected requiring one lepton (electron or muon), moderate missing transverse momentum, two jets with a large pseudorapidity separation and a large dijet invariant mass, and a signature consistent with the hadronic decay of a W/Z boson. The cross section is computed in a fiducial phase space defined at parton level requiring all parton transverse momenta $p_\mathrm{T}$$\gt$ 10 GeV and at least one pair of outgoing partons with invariant mass $m_\mathrm{qq}$$\gt$ 100 GeV. The measured and expected EW WV production cross sections are 1.90 $^{+0.53}_{-0.46}$ pb and 2.23 $^{+0.08}_{-0.11}$ (scale) $\pm$ 0.05 (PDF) pb, respectively, where PDF is the parton distribution function. The observed EW signal strength is $m_\mathrm{EW}$ = 0.85 $\pm$ 0.12 (stat) $^{+0.19}_{-0.17}$ (syst), corresponding to a signal significance of 4.4 standard deviations with 5.1 expected, and it is measured keeping the quantum chromodynamics (QCD) associated diboson production fixed to the standard model prediction. This is the first evidence of vector boson scattering in the $\ell\nu$qq decay channel at LHC. The simultaneous measurement of the EW and QCD associated diboson production agrees with the standard model prediction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Measurement of the $t\bar{t}$ charge asymmetry in events with highly Lorentz-boosted top quarks in pp collisions at $\sqrt{s}$ = 13 TeV

The measurement of the charge asymmetry in top quark pair events with highly Lorentz-boosted top quarks decaying to a single lepton and jets is presented. The analysis is performed using proton-proton collisions at $\sqrt{s}$ = 13 TeV with the CMS detector at the LHC and corresponding to an integrated luminosity of 138 fb -1 . The selection is optimized for top quarks produced with large Lorentz boosts, resulting in nonisolated leptons and overlapping jets. The top quark charge asymmetry is measured for events with a $t\bar{t}$ invariant mass larger than 750 GeV and corrected for detector and acceptance effects using a binned maximum likelihood fit. The measured top quark charge asymmetry of ($0.42_{-0.69}^{+0.64}$) % is in good agreement with the standard model prediction at next-to-next-to-leading order in quantum chromodynamic perturbation theory with next-to-leading-order electroweak corrections. The result is also presented for two invariant mass ranges, 750–900 and > 900 GeV .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Direct-photon production in inelastic and high-multiplicity proton–proton collisions at $\sqrt{s}$ = 13 TeV

In this letter, we present the first measurement of direct photons at the transverse momentum of 1 < 𝑝 T < 6 GeV/𝑐 at midrapidity 𝜂 < 0.8 in inelastic and high-multiplicity proton–proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV. The fraction of virtual direct photons in the inclusive virtual photon spectrum is obtained from a fit to the dielectron invariant mass spectrum. In the limit of zero invariant mass, this fraction is equal to the relative contribution of real direct photons in the inclusive real photon spectrum. Contributions from decays of light- flavour neutral mesons are estimated using independent measurements in proton–proton collisions at the same energy and the same event class. For the first time at the LHC energies, a direct-photon signal is observed at low 𝑝T in both inelastic and high-multiplicity event classes, with a significance of 3.2𝜎 and 1.9𝜎 in terms of standard deviations, correspondingly. The yield of direct photons in inelastic pp collisions is compared to perturbative QCD calculations. The integrated photon yield is studied as a function of charged-particle multiplicity and is compared to the results from other experiments and theoretical calculations. The results show a significant increase of direct- photon yield with charged-particle multiplicity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Adsorption of ethane, ethene, and ethyne on reconstructed Fe 3 O 4 (001)

The adsorption of the C2 hydrocarbons including ethane, ethene, and ethyne are studied on magnetite Fe 3 O 4 (001) by a combination of molecular beam dosing, temperature programmed desorption, and X-ray photoelectron spectroscopy. The ethane desorption profile has a single temperature invariant peak at 100 K while ethene and ethyne exhibit additional peaks at ~120 K and ~135 K. An inversion analysis is used to extract coverage dependent desorption energies as well as coverage-averaged prefactors for each molecule. Ethene and ethyne exhibit moderate coverage dependent desorption energies decreasing from ~80 to ~30 kJ/mol at saturation while ethane shows a relatively coverage invariant desorption energy ~28 kJ/mol. The desorption energies of the C2 hydrocarbons increase in order of increasing bond order. This is likely due to the enhanced interaction of the ethene and ethyne p system to the coordinatively unsaturated octahedral Fe sites (Fe oct ) on the ¬oxide surface. The saturation coverages of each C2 hydrocarbon at 82 K were also determined. Furthermore, these match well with the 2-dimensional spacing determined from their liquid densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

On heat conduction in an irregular magnetic field. Part 1

Anisotropic heat conduction in a plasma embedded in a magnetic field with irregular, possibly chaotic, field lines is discussed. If the collisional mean free path exceeds the electron gyroradius, the heat conductivity is much larger along the field lines than across them, and this enhances the transport across a domain where good flux surfaces do not exist. Recognising that anisotropic heat conduction may be cast in a variational form, and by constructing increasingly sophisticated trial functions that are based on invariant and almost-invariant structures under the magnetic field-line flow, bounds are derived on this enhancement and on the temperature variation along the magnetic field. In this way, remarkably accurate approximations for the temperature can be rapidly constructed without solving the diffusion equation, even in the small perpendicular-diffusion limit when the solution for the temperature is dominated by the fractal structure the magnetic field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electron kinetics in a high- Z plasmoid

The problem of the electron dynamics on a closed magnetic field line passing through a high-Z plasmoid is considered. The electron kinetic equation is integrated over bounce motion and pitch angle, reducing the independent variables to a single adiabatic invariant plus time. Integration of the full Landau self-collision operator is carried out exactly, resulting in a nonlinear integro-differential operator in the new invariant. Conservation laws and the H theorem of the integrated self-collision operator are proven. Numerical solutions of the integrated kinetic equation are obtained with a self-consistent quasineutral electric potential, given the initial condition of a cold plasmoid immersed in a hot ambient plasma. The fact that cold electrons are deeply trapped in a potential with a parabolic peak leads to exactly 3/4 the usual rate of collisional heating by the ambient plasma, independent of any other parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Narrow operator models of stellarator equilibria in Fourier Zernike basis

Numerical computation of the ideal magnetohydrodynamic (MHD) equilibrium magnetic field is at the base of stellarator optimisation and provides the starting point for solving more sophisticated partial differential equations like transport or turbulence models. Conventional approaches solve for a single stationary point of the ideal MHD equations, which is fully defined by three invariants and the numerical scheme employed by the solver. We present the first numerical approach that can solve for a continuous distribution of equilibria with fixed boundary and rotational transform, varying only the pressure invariant. This approach minimises the force residual by optimising parameters of multilayer perceptrons that map from a scalar pressure multiplier to the Fourier Zernike basis as implemented in the modern stellarator equilibrium solver DESC.

fusion plasma↗
DOE OSTI2023

Topological Bismuth (1 1̅ 0) Facet for Efficient Oxygen Reduction Cathode in Fuel Cells

The performance of proton-exchange membrane fuel cells is strongly dependent on the efficiency of the electrochemical oxygen reduction reaction (ORR). Here, we report both two-dimensional (4 monolayer, ML) and three-dimensional Bi (16 ML) (11̅0) facets deliver a strong correlation between nontrivial topological surface states and low electrochemical overpotential toward ORR. Both topological Bi (11̅0) film and slab have two slightly gapped (9–45 meV) Dirac cones with the invariant Z 2 =1, derived from the parity products for the occupied states at four time-reversal invariant momentum points. By contrast, the Bi (110) film is stabilized as an ordinary insulator due to atomic buckling. The much lower thermodynamic ORR overpotential (η = 0.58 V) and higher H 2 O selectivity by four steps of H + /e – transfer are obtained on the topological Bi (11̅0) surface compared with the nontopological (110) facet (η = 0.84 V). This promotional effect can be attributed to the accumulation of topological surface bands (Bi 6p orbitals) near the Fermi energy, which provide stronger electronic transport channels along the surface Bi–Bi and Bi–O bonds. Finally, a 16 ML Bi (11̅0) slab is also constructed to simulate the topological (11̅0) surface of three-dimensional bulk Bi as a topological crystalline insulator (TCI). Furthermore, the consistent topological properties and slightly lower thermodynamic ORR overpotential (η = 0.53 V) are also observed on the thicker slab.

25 ENERGY STORAGE↗
DOE OSTI2021

Quantification of Mono- and Multivalent Counterion-Mediated Bridging in Polyelectrolyte Brushes

Multivalent counterion-induced bridging interactions have been identified as the key mechanism of drastic collapse of the height of polyelectrolyte (PE) brushes. In this article, we employ all-atom molecular dynamics (MD) simulations to quantify the bridging interactions in PE brushes for counterions of different sizes and valences. We identify that unlike the current notion, bridging interactions are not the sole function of the counterion valence. Rather the bridging interactions depend on the fraction of counterions (of a given type) that get physically condensed on the PE backbone as well as the size of the counterion solvation shell. These mechanisms ensure that certain monovalent counterions demonstrate much stronger bridging interactions than that witnessed for certain divalent and trivalent counterions, while certain counterions of identical valences show drastically different bridging. Here, we argue that these counterion-specific bridging interactions eventually enable not only the significant reduction of the PE brush height in presence of certain multivalent screening counterions, but may also give rise to scenarios where the brush height reduction for certain monovalent counterions is larger than certain divalent and trivalent counterions. This latter observation contradicts the experimental findings where the multivalent counterions invariably led to a larger decrease in the height of the PE brushes: we argue that this discrepancy stems from the fact that in our simulations we only consider densely grafted and short (and hence less flexible) PE brushes that hinder the formation of different laterally inhomogeneous structures (like pinned micelles and cylindrical bundles) that would have led to a larger brush height reduction (in experiments, which invariably consider longer and less densely grafted brushes, the formation of such inhomogeneous structures are primarily responsible for larger brush height reduction in presence of multivalent counterions). Finally, we also probe the dynamic properties of the counterions (i.e., their time-dependent displacements) and their bridging interactions (i.e., lifetime of bridging interactions).

42 ENGINEERING↗
DOE OSTI2022

Artificial Graphene Nanoribbons: A Test Bed for Topology and Low-Dimensional Dirac Physics

Here, we synthesize artificial graphene nanoribbons by positioning carbon monoxide molecules on a copper surface to confine its surface state electrons into artificial atoms positioned to emulate the low-energy electronic structure of graphene derivatives. We demonstrate that the dimensionality of artificial graphene can be reduced to one dimension with proper "edge" passivation, with the emergence of an effectively gapped one-dimensional nanoribbon structure. These one-dimensional structures show evidence of topological effects analogous to graphene nanoribbons. Guided by first-principles calculations, we spatially explore robust, zero-dimensional topological states by altering the topological invariants of quasi-one-dimensional artificial graphene nanostructures. The robustness and flexibility of our platform allow us to toggle the topological invariants between trivial and nontrivial on the same nanostructure. Ultimately, we spatially manipulate the states to understand fundamental coupling between adjacent topological states that are finely engineered and simulate complex Hamiltonians.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Observation of Persistent Zero Modes and Superconducting Vortex Doublets in UTe 2

Superconducting vortices can reveal electron pairing details and nucleate topologically protected states. Yet, vortices of bulk spin-triplet superconductors have never been visualized at the atomic scale. Recently, UTe 2 has emerged as a prime spin-triplet superconductor, but its superconducting order parameter is elusive, and whether time-reversal symmetry (TRS) is broken remains unsettled. Here, in this work, we visualize vortices on the (011) surface of ultraclean UTe 2 single crystals (T c = 2.1 K) using scanning tunneling microscopy (STM). We introduce $\frac{d^{2}𝐼}{d𝑉^{2}}$ imaging as an effective technique for vortex visualization in superconductors with substantial residual zero-energy density of states (DOS), as in UTe 2 . Anisotropic single-flux-quantum vortices, with coherence lengths of ∼12 nm (4 nm) parallel (perpendicular) to the a -axis, form a triangular vortex lattice (VL) under a small out-of-plane magnetic field. The invariance of vortex structures and VL under changes of field polarity and cooling history strongly supports time-reversal invariant superconductivity under zero field. At vortex cores (VCs), nonsplit, spectrally sharp zero-bias conductance peaks (ZBPs) persist up 8 T that are consistent with symmetry-protected Majorana zero modes (MZMs) in a topological vortex line. Close examination of vortex structures reveals a mirror-asymmetric doublet─one with ZBPs and another with an enhanced apparent gap, possibly originating from a field-induced multicomponent order parameter.

UTe2↗
DOE OSTI2023

Bias-Adjustment Methods for Future Subdaily Precipitation Extremes Consistent Across Durations

Model output from climate projections often requires bias-adjustment to compensate for systematic model errors. A bias-adjustment method for extreme precipitation intensity is proposed that preserves the scaling equation for different accumulation levels from hourly to daily, using intensity-duration-frequency (IDF) modeling. A validation is performed within a pseudo-reality setting, based on hourly precipitation from 28 regional climate model projections of the EURO-CORDEX ensemble over Belgium. The scaling-based adjustment methods improve upon previous methods, an optimal method is identified, and, analytical quantile mapping methods must be avoided due to three identified problems. The ensemble mean of the adjusted extreme precipitation intensity obeys the above-mentioned scale-invariance property, which is consistent with observed extreme intensities. We thus show that IDF modeling provides added value in the context of bias-adjustment, and, that the particular IDF model proposed balances well between accuracy and the preservation of desired properties such as scale invariance and consistency among rainfall durations.

54 ENVIRONMENTAL SCIENCES↗