Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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 325 records · Page 18

Quasinormal g-modes of neutron stars with quarks

Quasinormal oscillation modes of neutron stars provide a means to probe their interior composition using gravitational wave astronomy. We compute the frequencies and damping times of composition-dependent core 𝑔-modes of neutron stars containing quark matter employing linearized perturbative equations of general relativity. We find that ignoring background metric perturbations due to the oscillating fluid, as in the Cowling approximation, underestimates the 𝑔-mode frequency by up to 10% for higher mass stars, depending on the parameters of the nuclear equation of state and how the mixed phase is constructed. The 𝑔-mode frequencies are well described by a linear scaling with the central lepton (or combined lepton and quark) fraction for nucleonic (hybrid) stars. Furthermore, our findings suggest that neutron stars with and without quarks are manifestly different with regards to their quasinormal 𝑔-mode spectrum, and may thus be distinguished from one another in future observations of gravitational waves from merging neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Primitive quantum gates for dihedral gauge theories

We describe the simulation of dihedral gauge theories on digital quantum computers. The nonabelian discrete gauge group $D_N$-- the dihedral group -- serves as an approximation to $U(1)\times\mathbb{Z}_2$ lattice gauge theory. In order to carry out such a lattice simulation, we detail the construction of efficient quantum circuits to realize basic primitives including the nonabelian Fourier transform over $D_N$, the trace operation, and the group multiplication and inversion operations. For each case the required quantum resources scale linearly or as low-degree polynomials in $n=\log N$. We experimentally benchmark our gates on the Rigetti Aspen-9 quantum processor for the case of $D_4$. The fidelity of all $D_4$ gates was found to exceed $80\%$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fragile Topology and Flat-Band Superconductivity in the Strong-Coupling Regime

In flat bands, superconductivity can lead to surprising transport effects. The superfluid “mobility”, in the form of the superfluid weight D s , does not draw from the curvature of the band but has a purely band-geometric origin. In a mean-field description, a nonzero Chern number or fragile topology sets a lower bound for D s , which, via the Berezinskii-Kosterlitz-Thouless mechanism, might explain the relatively high superconducting transition temperature measured in magic-angle twisted bilayer graphene (MATBG). For fragile topology, relevant for the bilayer system, the fate of this bound for finite temperature and beyond the mean-field approximation remained, however, unclear. Here, we numerically use exact Monte Carlo simulations to study an attractive Hubbard model in flat bands with topological properties akin to those of MATBG. We find a superconducting phase transition with a critical temperature that scales linearly with the interaction strength. Then, we investigate the robustness of the superconducting state to the addition of trivial bands that may or may not trivialize the fragile topology. Our results substantiate the validity of the topological bound beyond the mean-field regime and further stress the importance of fragile topology for flat-band superconductivity.

2-dimensional systems↗

Optimality and Noise Resilience of Critical Quantum Sensing

We compare critical quantum sensing to passive quantum strategies to perform frequency estimation, in the case of single-mode quadratic Hamiltonians. We show that, while in the unitary case both strategies achieve precision scaling quadratic with the number of photons, in the presence of dissipation this is true only for critical strategies. We also establish that working at the exceptional point or beyond threshold provides suboptimal performance. This critical enhancement is due to the emergence of a transient regime in the open critical dynamics, and is invariant to temperature changes. When considering both time and system size as resources, for both strategies the precision scales linearly with the product of the total time and the number of photons, in accordance with fundamental bounds. However, we show that critical protocols outperform optimal passive strategies if preparation and measurement times are not negligible. Our results are applicable to a broad variety of critical sensors whose phenomenology can be reduced to that of a single-mode quadratic Hamiltonian, including systems described by finite-component and fully connected models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measuring the angular momentum of a neutron using Earth's rotation

A coupling between Earth's rotation and orbital angular momentum (OAM), known as the Sagnac effect, is observed in entangled neutrons produced using a spin-echo interferometer. After correction for instrument systematics the measured coupling is within 5% of theory, with an uncertainty of 7.2%. The OAM in our setup is transverse to the propagation direction and scales linearly with neutron wavelength (4–12.75 Å), so the Sagnac coupling can be varied without mechanically rotating the device, which avoids systematic errors present in previous experiments. The detected transverse OAM of our beam corresponds to 4098 ± 295 ℏ Å − 1 , 10 5 times lower than in the previous neutron experiments. This demonstrates the feasibility of using the Sagnac effect to definitively measure neutron OAM and paves the way towards a future observation of the quantum Sagnac effect. Published by the American Physical Society 2025

Geerits, Niels (ORCID:0000000248072084)↗

Scattering neutrinos, spin models, and permutations

We consider a class of Heisenberg all-to-all coupled spin models inspired by neutrino interactions in a supernova with N degrees of freedom. These models are characterized by a coupling matrix that is relatively simple in the sense that there are only a few, relative to N , nontrivial eigenvalues, in distinction to the classic Heisenberg spin-glass models, leading to distinct behavior in both the high-temperature and low-temperature regimes. When the momenta of the neutrinos are uniform and random in directions, we can calculate the large- N partition function for the quantum Heisenberg model. In particular, the high-temperature partition function predicts a non-Gaussian density of states, providing interesting counterexamples showing the limits of general theorems on the density of states for quantum spin models. We can repeat the same argument for classical Heisenberg models, also known as rotor models, and we find the high-temperature expansion is completely controlled by the eigenvalues of the coupling matrix, and again predicts non-Gaussian behavior for the density of states as long as the number of eigenvalues does not scale linearly with N . Indeed, we derive the amusing fact that these partition functions are essentially the generating function for counting permutations in the high-temperature regime. Finally, for the case relevant to neutrinos in a supernova, we identify the low-temperature phase as a unique state with the direction of the momenta of the neutrino dictating its coherent state in flavor-space, a state we dub the “flavor-momentum-locked” state. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamic, symmetry-preserving, and hardware-adaptable circuits for quantum computing many-body states and correlators of the Anderson impurity model

We present a hardware-reconfigurable ansatz on N q -qubits for the variational preparation of many-body states of the Anderson impurity model (AIM) with N imp + N bath = N q /2 sites, which conserves total charge and spin z component within each variational search subspace. The many-body ground state of the AIM is determined as the minimum over all minima of O(N$^2_ q$) distinct charge-spin sectors. Hamiltonian expectation values are shown to require ω(N q ) < N meas. $\leqslant$ O(N imp N bath ) symmetry-preserving, parallelizable measurement circuits, each amenable to postselection. To obtain the one-particle impurity Green’s function we show how initial Krylov vectors can be computed via midcircuit measurement and how Lanczos iterations can be computed using the symmetry-preserving ansatz. For a single-impurity Anderson model with a number of bath sites increasing from one to seven, we show using numerical emulation that the ease of variational ground-state preparation is suggestive of linear scaling in circuit depth and subquartic scaling in optimizer complexity. We therefore expect that, combined with time-dependent methods for Green’s function computation, our ansatz provides a useful tool to account for electronic correlations on early fault-tolerant processors. Finally, with a view towards computing real materials properties of interest like magnetic susceptibilities and electron-hole propagators, we provide a straightforward method to compute many-body, time-dependent correlation functions using a combination of time evolution, midcircuit measurement-conditioned operations, and the Hadamard test.

36 MATERIALS SCIENCE↗

Entanglement Witness for Indistinguishable Electrons Using Solid-State Spectroscopy

Characterizing entanglement in quantum materials is crucial for advancing next-generation quantum technologies. Despite recent strides in witnessing entanglement in magnetic materials with distinguishable spin modes, quantifying entanglement in systems formed by indistinguishable electrons remains a formidable challenge. To solve this problem, we introduce a method to extract various four-fermion correlations by analyzing the nonlinearity in resonant inelastic x-ray scattering spectra. These correlations constitute the primary components of the cumulant two-particle reduced density matrix. We further derive bounds for its eigenvalues and demonstrate the linear scaling with fermionic entanglement depth, providing a reliable witness for entanglement. Using the material-relevant strongly correlated models as examples, we show how this entanglement witness can efficiently quantify multipartite entanglement across different phase regions, highlighting its advantage over quantum Fisher information. Published by the American Physical Society 2025

Liu, Tongtong (ORCID:0000000295324061)↗

SpeckleNN: a unified embedding for real-time speckle pattern classification in X-ray single-particle imaging with limited labeled examples

With X-ray free-electron lasers (XFELs), it is possible to determine the three-dimensional structure of noncrystalline nanoscale particles using X-ray single-particle imaging (SPI) techniques at room temperature. Classifying SPI scattering patterns, or `speckles', to extract single-hits that are needed for real-time vetoing and three-dimensional reconstruction poses a challenge for high-data-rate facilities like the European XFEL and LCLS-II-HE. Here, we introduce SpeckleNN, a unified embedding model for real-time speckle pattern classification with limited labeled examples that can scale linearly with dataset size. Trained with twin neural networks, SpeckleNN maps speckle patterns to a unified embedding vector space, where similarity is measured by Euclidean distance. We highlight its few-shot classification capability on new never-seen samples and its robust performance despite having only tens of labels per classification category even in the presence of substantial missing detector areas. Without the need for excessive manual labeling or even a full detector image, our classification method offers a great solution for real-time high-throughput SPI experiments.

47 OTHER INSTRUMENTATION↗

Hvac: Removing I/O Bottleneck for Large-Scale Deep Learning Applications

Scientific communities are increasingly adopting deep learning (DL) models in their applications to accelerate scientific discovery processes. However, with rapid growth in the computing capabilities of HPC supercomputers, large-scale DL applications have to spend a significant portion of training time performing I/O to a parallel storage system. Previous research works have investigated optimization techniques such as prefetching and caching. Unfortunately, there exist non-trivial challenges to adopting the existing solutions on HPC supercomputers for large-scale DL training applications, which include non-performance and/or failures at extreme scale, lack of portability and generality in design, complex deployment methodology, and being limited to a specific application or dataset. To address these challenges, we propose High-Velocity AI Cache (HVAC), a distributed read-cache layer that targets and fully exploits the node-local storage or near node-local storage technology. HVAC seamlessly accelerates read I/O by aggregating node-local or near node-local storage, avoiding metadata lookups and file locking while preserving portability in the application code. We deploy and evaluate HVAC on 1,024 nodes (with over 6000 NVIDIA V100 GPUS) of the Summit supercomputer. In particular, we evaluate the scalability, efficiency, accuracy, and load distribution of HVAC compared to GPFS and XFS-on-NVMe. With four different DL applications, we observe an average 25 % performance improvement atop GPFS and 9% drop against XFS-on-NVMe, which scale linearly and are considered the performance upper bound. We envision HVAC as an important caching library for upcoming HPC supercomputers such as Frontier.

Khan, Awais↗

Parallel String Graph Construction and Transitive Reduction for De Novo Genome Assembly

One of the most computationally intensive tasks in computational biology is de novo genome assembly, the decoding of the sequence of an unknown genome from redundant and erroneous short sequences. A common assembly paradigm identifies overlapping sequences, simplifies their layout, and creates consensus. Despite many algorithms developed in the literature, the efficient assembly of large genomes is still an open problem. In this work, we introduce new distributed-memory parallel algorithms for overlap detection and layout simplification steps of de novo genome assembly, and implement them in the diBELLA 2D pipeline. Our distributed memory algorithms for both overlap detection and layout simplification are based on linear-algebra operations over semirings using 2D distributed sparse matrices. Our layout step consists of performing a transitive reduction from the overlap graph to a string graph. We provide a detailed communication analysis of the main stages of our new algorithms. diBELLA 2D achieves near linear scaling with over 80% parallel efficiency for the human genome, reducing the runtime for overlap detection by 1.2-1.3× for the human genome and 1.5-1.9× for C.elegans compared to the state-of-the-art. Our transitive reduction algorithm outperforms an existing distributed-memory implementation by 10.5-13.3× for the human genome and 18-29× for the C. elegans. Our work paves the way for efficient de novo assembly of large genomes using long reads in distributed memory.

59 BASIC BIOLOGICAL SCIENCES↗

Catalytic resonance of ammonia synthesis by simulated dynamic ruthenium crystal strain

Ammonia affords dense storage for renewable energy as a fungible liquid fuel, provided it can be efficiently synthesized from hydrogen and nitrogen. In this work, the catalysis of ammonia synthesis was computationally explored beyond the Sabatier limit by dynamically straining a ruthenium crystal (±4%) at the resonant frequencies (10 2 to 10 5+ Hz) of N 2 surface dissociation and hydrogenation. Density functional theory calculations at different strain conditions indicated that the energies of NH x surface intermediates and transition states scale linearly, allowing the description of ammonia synthesis at a continuum of strain conditions. A microkinetic model including multiple sites and surface diffusion between step and Ru(0001) terrace sites of varying ratios for nanoparticles of differing size revealed that dynamic strain yields catalytic ammonia synthesis conversion and turnover frequency comparable to industrial reactors (400°C, 200 atm) but at lower temperature (320°C) and an order of magnitude lower pressure (20 atm).

10 SYNTHETIC FUELS↗

Large magnetoelectric resistance in the topological Dirac semimetal α-Sn

The spin-momentum locking of surface states in topological materials can produce a resistance that scales linearly with magnetic and electric fields. Such a bilinear magnetoelectric resistance (BMER) effect offers a new approach for information reading and field sensing applications, but the effects demonstrated so far are too weak or for low temperatures. This article reports the first observation of BMER effects in topological Dirac semimetals; the BMER responses were measured at room temperature and were substantially stronger than those reported previously. The experiments used topological Dirac semimetal α-Sn thin films grown on silicon substrates. The films showed BMER responses that are 10 6 times larger than previously measured at room temperature and are also larger than those previously obtained at low temperatures. These results represent a major advance toward realistic BMER applications. Significantly, the data also yield the first characterization of three-dimensional Fermi-level spin texture of topological surface states in α-Sn.

36 MATERIALS SCIENCE↗

Review—Development of Highly Active and Durable Hybrid Compressive Platinum Lattice Catalysts for Polymer Electrolyte Membrane Fuel Cells: Mathematical Modeling and Experimental Work

This review provides a comprehensive overview on the development of highly active and durable platinum catalysts with ultra-low Pt loadings for polymer electrolyte membrane fuel cells (PEMFCs) through a combined mathematical modeling and experimental work. First, simulation techniques were applied to evaluate the validity of the Tafel approximation for the calculation of the mass activity (MA) and specific activity (SA). A one-dimensional agglomeration model was developed and solved to understand the effects of exchange current density, porosity, agglomerate size, Nafion ® film thickness, and Pt loading on the MA and SA. High porosity (> 60%) and agglomerations at high Pt loadings cause the loss of the Tafel approximation and consequently the decrease in MA and SA. A new structure parameter was introduced to estimate the real porous structure using the fractal theory. The volumetric catalyst density was corrected by the fractal dimension (measured by Hg porosimetry), which gave a good agreement with the experimental values. The loading-dependent Tafel equation was then derived, which contains both the utilization and the non-linear scaling factor. Second, activated carbon composite support (ACCS) with optimized surface area, porosity, pore size, and pore size distribution was developed. The hydrophilic/hydrophobic ratio, structural properties (amorphous/crystalline ratio), and the number of active sites were optimized through metal-catalyzed pyrolysis. Stability of ACCS and Pt/ACCS were evaluated using an accelerated stress test (AST). The results indicated that Pt/ACCS showed no significant loss of MA and power density after 5,000 cycles at 1.0–1.5 V, while the commercial Pt/C catalysts showed drastic losses of MA and power density. Finally, monolayers of compressed Pt (core–shell-type Pt 3 Co 1 ) catalysts were structured by diffusing Co atoms (previously embedded in ACCS) into Pt. Compressive Pt lattice (Pt * ) catalysts were synthesized through an annealing procedure developed at the University of South Carolina (USC). The Pt * /ACCS catalyst showed high initial power density (rated) of 0.174 g Pt kW −1 and high stability (24 mV loss) at 0.8 A cm −2 after 30,000 cycles (0.6–1.0 V). The outstanding performance of Pt * /ACCS is due to the synergistic effect of ACCS and compressive Pt * lattice.

25 ENERGY STORAGE↗

Enabling particle applications for exascale computing platforms

The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four “sub-motifs”: short-range particle–particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed particle hydrodynamics (SPH) methods), long-range particle–particle interactions (e.g., electrostatic MD and gravitational N-body), particle-in-cell (PIC) methods, and linear-scaling electronic structure and quantum molecular dynamics (QMD) algorithms. Our crosscutting co-designed technologies fall into two categories: proxy applications (or “apps”) and libraries. Proxy apps are vehicles used to evaluate the viability of incorporating various types of algorithms, data structures, and architecture-specific optimizations and the associated trade-offs; examples include ExaMiniMD, CabanaMD, CabanaPIC, and ExaSP2. Libraries are modular instantiations that multiple applications can utilize or be built upon; CoPA has developed the Cabana particle library, PROGRESS/BML libraries for QMD, and the SWFFT and fftMPI parallel FFT libraries. Success is measured by identifiable “lessons learned” that are translated either directly into parent production application codes or into libraries, with demonstrated performance and/or productivity improvement. The libraries and their use in CoPA’s ECP application partner codes are also addressed.

97 MATHEMATICS AND COMPUTING↗

Exploring limitations in the induced polarization versus surface conductivity relationship in the case of wetland soils

Recent induced polarization studies suggest that the real part of surface conductivity (σ' surf ) scales linearly with the imaginary conductivity (σ'' = σ'' surf ) or normalized chargeability (Mn) for a range of soil types. The coefficients of this relationship l and l_Mn ( l = σ''/σ' surf or l_Mn = Mn/σ' surf ) allow the separation of the surface and electrolytic conductivities from the bulk conductivity. However, the dependence of these constants on varying soil physicochemical properties, including under unsaturated conditions, is yet to be assessed. Using estimates of σ' surf from 18 undisturbed soil samples from a restored wetland and σ'' measured over a frequency range of 0.01 Hz to 10 kHz, the σ' surf and σ'' were compared with the laboratory measurements of soil properties. Also, l and l_Mn were calculated for each soil sample and regressed them against the soil properties. We find an apparent dependence of l on soil texture, bulk density, organic matter, and moisture contents, with coefficients of determination (R 2 ) ranging from 0.5 to 0.65 at low frequencies (e.g., 1 Hz) but not at high frequencies (e.g., 936 Hz). This dependence of l on soil texture results from the insensitivity of σ'' at low frequency to σ' surf and, by implication, to the soil properties controlling σ' surf . In contrast, l_Mn indicates no correlation with the soil properties because Mn is linearly correlated with σ' surf and correlated with the soil properties controlling σ' surf . Finally, our results call for caution on the application of σ'' at a single frequency as a proxy of σ' surf because σ'' is not necessarily correlated with σ' surf across all soil types. Although using l_Mn derived from multifrequency measurements overcomes this limitation, field acquisition of spectral information (e.g., up to 1000 Hz) remains a challenge.

54 ENVIRONMENTAL SCIENCES↗

Pattern formation on free surfaces via plastic buckling and periodic folding

Spontaneous formation of surface patterns such as folds and creases is widely observed in engineered systems, such as thin films and machined surfaces, as well as in natural systems such as floral petals and leaves. A recognized route for their formation is via surface elastic instabilities. Here we demonstrate an alternate route for their formation —plastic buckling and folding—that occurs when ductile materials are subject to large-strain (>1) deformation. Using in situ imaging of simple-shear deformation of polypropylene, we delineate the mechanics underlying formation of folds and the resulting surface morphology. Here, the fold pattern shows a single wavelength that scales linearly with the size of the plastic zone. An analytical model predicts onset of the plastic buckling instability and characteristics of the folds. Besides obvious applications for controllably generating surface features over large areas, our results suggest a contributing mechanism for origin of surface roughness in solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On-chip terahertz emission from Floquet-Bloch states [Invited]

Floquet engineering uses time-periodic electromagnetic fields to modify the electronic properties of quantum materials via the creation of Floquet-Bloch states. These photon-dressed states inherit features from both the material and the driving field, enabling the exploration and control of quantum phenomena in light-matter hybrid systems. In non-centrosymmetric materials, shift currents can arise from the quantum geometric properties of electronic wavefunctions. However, shift currents from Floquet-Bloch states remain experimentally unexplored. Here, we employ an on-chip optoelectronic circuit to detect intrinsic terahertz emission from Floquet-Bloch states in T d -WTe 2 under intense optical driving. We observe strong edge-localized terahertz emission that scales linearly with the driving field, consistent with the theoretical prediction for shift currents generated by Floquet-Bloch states. The results advance our understanding of strongly driven quantum materials and provide insights for developing efficient, bias-free terahertz sources for future optoelectronic technologies.

Li, Xinyu [Max Planck Institute for the Structure ↗