SEARCH · Engineering Papers

Results for “Butterflies”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS2018

An Engergetic Electron Flux Dropout Due to Magnetopause Shadowing on 1 June 2013

We examine the mechanisms responsible for the dropout of energetic electron flux during 31 May to 1 June 2013 using Van Allen Probe (Radiation Belt Storm Probes (RBSP)) electron flux data and simulations with the Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. During the storm main phase, L-shells at RBSP locations are greater than ~8, which are connected to open drift shells. Consequently, diminished electron fluxes were observed over a wide range of energies. The combination of drift shell splitting, magnetopause shadowing, and drift loss all results in butterfly electron pitch angle distributions (PADs) at the nightside. During storm sudden commencement, RBSP observations display electron butterfly PADs over a wide range of energies. However, it is difficult to determine whether there are butterfly PADs during the storm main phase since the maximum observable equatorial pitch angle from RBSP is not larger than ~40° during this period. To investigate the causes of the dropout, the CIMI model is used as a global 4-D kinetic inner magnetosphere model. The CIMI model reproduces the dropout with very similar timing and flux levels and PADs along the RBSP trajectory for 593 keV. Furthermore, the CIMI simulation shows butterfly PADs for 593 keV during the storm main phase. Based on comparison of observations and simulations, we suggest that the dropout during this event mainly results from magnetopause shadowing.

Kang, Suk-Bin↗

DOE OSTI2025

Spectroscopy of the fractal Hofstadter energy spectrum

Hofstadter’s butterfly, the predicted energy spectrum for non-interacting electrons confined to a two-dimensional lattice in a magnetic field, is one of the most remarkable fractal structures in nature. At rational ratios of magnetic flux quanta per lattice unit cell, this spectrum shows self-similar distributions of energy levels that reflect its recursive construction. For most materials, Hofstadter’s butterfly is predicted under experimental conditions that are unachievable using laboratory-scale magnetic fields. More recently, electrical transport studies have provided evidence for Hofstadter’s butterfly in materials engineered to have artificially large lattice constants, such as those with moiré superlattices. Yet, so far, direct spectroscopy of the fractal energy spectrum predicted by Hofstadter nearly 50 years ago has remained out of reach. Here, in this work, we use high-resolution scanning tunnelling microscopy/spectroscopy (STM/STS) to investigate the flat electronic bands in twisted bilayer graphene (TBG) near the predicted second magic angle, an ideal setting for spectroscopic studies of Hofstadter’s spectrum. Our study shows the fractionalization of flat moiré bands into discrete Hofstadter subbands and discerns experimental signatures of self-similarity of this spectrum. Moreover, our measurements uncover a spectrum that evolves dynamically with electron density, showing phenomena beyond that of Hofstadter’s original model owing to the combined effects of strong correlations, Coulomb interactions and the quantum degeneracy of electrons in TBG.

electronic properties and materials↗

DOE OSTI2025

Topological Phenomena in Artificial Quantum Materials Revealed by Local Chern Markers

A striking example of frustration in physics is Hofstadter’s butterfly, a fractal structure that emerges from the competition between a crystal’s lattice periodicity and the magnetic length of an applied field. Current methods for predicting the topological invariants associated with Hofstadter’s butterfly are challenging or impossible to apply to a range of materials, including those that are disordered or lack a bulk spectral gap. Here, in this work, we demonstrate a framework for predicting a material’s local Chern markers using its position-space description and validate it against experimental observations of quantum transport in artificial graphene in a semiconductor heterostructure, inherently accounting for fabrication disorder strong enough to close the bulk spectral gap. By resolving local changes in the system’s topology, we reveal the topological origins of antidot-localized states that appear in artificial graphene in the presence of a magnetic field. Moreover, we show the breadth of this framework by simulating how Hofstadter’s butterfly emerges from an initially unpatterned 2D electron gas as the system’s potential strength is increased and predict that artificial graphene becomes a topological insulator at the critical magnetic field. Overall, we anticipate that a position-space approach to determine a material’s Chern invariant without requiring prior knowledge of its occupied states or bulk spectral gaps will enable a broad array of fundamental inquiries and provide a novel route to material discovery, especially in metallic, aperiodic, and disordered systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1987

Magnetic field drift shell splitting - Cause of unusual dayside particle pitch angle distributions during storms and substorms

This paper presents a magnetic field drift shell-splitting model for the unusual butterfly and head-and-shoulder energetic (E greater than 25 keV) particle pitch angle distributions (PADs) which appear deep within the dayside magnetosphere during the course of storms and substorms. Drift shell splitting separates the high and low pitch angle particles in nightside injections as they move to the dayside magnetosphere, so that the higher pitch angle particles move radially away from earth. Consequently, butterfly PADs with a surplus of low pitch angle particles form on the inner edge of the injection, but head-and-shoulder PADs with a surplus of high pitch angle particles form on the outer edge. A similar process removes high pitch angle particles from the inner dayside magnetosphere during storms, leaving the remaining lower pitch angle particles to form butterfly PADs on the inner edge of the ring current. A detailed case and statistical study of Charge Composition Explorer/Medium-energy Particle Analyzer observations, as well as a review of previous work, shows most examples of unusual PADs to be consistent with the model.

Sibeck, D. G.↗

NASA NTRS1992

Aerodynamic Shutoff Valve

Aerodynamic flow achieved by adding fixed fairings to butterfly valve. When valve fully open, fairings align with butterfly and reduce wake. Butterfly free to turn, so valve can be closed, while fairings remain fixed. Design reduces turbulence in flow of air in internal suction system. Valve aids in development of improved porous-surface boundary-layer control system to reduce aerodynamic drag. Applications primarily aerospace. System adapted to boundary-layer control on high-speed land vehicles.

Horstman, Raymond H.↗

NASA NTRS1992

A Benes-like theorem for the shuffle-exchange graph

One of the first theorems on permutation routing, proved by V. E. Beness (1965), shows that given a set of source-destination pairs in an N-node butterfly network with at most a constant number of sources or destinations in each column of the butterfly, there exists a set of paths of lengths O(log N) connecting each pair such that the total congestion is constant. An analogous theorem yielding constant-congestion paths for off-line routing in the shuffle-exchange graph is proved here. The necklaces of the shuffle-exchange graph play the same structural role as the columns of the butterfly in Beness' theorem.

Schwabe, Eric J.↗

NASA NTRS2011

High-Throughput, Adaptive FFT Architecture for FPGA-Based Spaceborne Data Processors

Exponential growth in microelectronics technology such as field-programmable gate arrays (FPGAs) has enabled high-performance spaceborne instruments with increasing onboard data processing capabilities. As a commonly used digital signal processing (DSP) building block, fast Fourier transform (FFT) has been of great interest in onboard data processing applications, which needs to strike a reasonable balance between high-performance (throughput, block size, etc.) and low resource usage (power, silicon footprint, etc.). It is also desirable to be designed so that a single design can be reused and adapted into instruments with different requirements. The Multi-Pass Wide Kernel FFT (MPWK-FFT) architecture was developed, in which the high-throughput benefits of the parallel FFT structure and the low resource usage of Singleton s single butterfly method is exploited. The result is a wide-kernel, multipass, adaptive FFT architecture. The 32K-point MPWK-FFT architecture includes 32 radix-2 butterflies, 64 FIFOs to store the real inputs, 64 FIFOs to store the imaginary inputs, complex twiddle factor storage, and FIFO logic to route the outputs to the correct FIFO. The inputs are stored in sequential fashion into the FIFOs, and the outputs of each butterfly are sequentially written first into the even FIFO, then the odd FIFO. Because of the order of the outputs written into the FIFOs, the depth of the even FIFOs, which are 768 each, are 1.5 times larger than the odd FIFOs, which are 512 each. The total memory needed for data storage, assuming that each sample is 36 bits, is 2.95 Mbits. The twiddle factors are stored in internal ROM inside the FPGA for fast access time. The total memory size to store the twiddle factors is 589.9Kbits. This FFT structure combines the benefits of high throughput from the parallel FFT kernels and low resource usage from the multi-pass FFT kernels with desired adaptability. Space instrument missions that need onboard FFT capabilities such as the proposed DESDynl, SWOT (Surface Water Ocean Topography), and Europa sounding radar missions would greatly benefit from this technology with significant reductions in non-recurring cost and risk.

NguyenKobayashi, Kayla↗

NASA NTRSDate not supplied

Quantifying Leak Rates and Vacuum Impact on Granular Soil Column Pneumatic Seal

Introduction: All In-Situ Resource Utilization (ISRU) regolith reactor systems that aim to extract gases or vapors must maintain a seal along the flow path of the regolith from the inlet hopper to where extracted gases or vapors are generated. Using a column of regolith at the inlet has been proposed as an alternative to systems using series of mechanical valves to reduce mass and complexity. Preliminary investigations into the use of regolith itself as a barrier to escaping fluids have been conducted, which this research seeks to expand upon. Different vertical regolith simulant plug diameters and heights were tested under vacuum with dry nitrogen to observe leak rates and approximate steady state pressure gradient equilibria. The effect of absolute ambient pressure was also studied by repeating tests at incremental pressures ranging from 10e-5 Torr up to 1 atm to corroborate ambient bench tests from previous work. The impact of regolith dynamics was addressed to emulate continuous-feed regolith reaction processes by measuring leak rates through continuously draining simulant columns. OB-1 lunar highlands regolith simulant was used. Experimental Setup: Regolith simulant columns were constructed using a motorized butterfly valve separating upper acrylic tubes of various diameters from a larger, lower aluminum pipe plenum (See Figure 1). The butterfly valve was cracked open such that gas could freely pass through whereas regolith would bridge and stay suspended as a cohesive plug above in the acrylic portion. These assemblies were loaded into a vacuum chamber, and a compressed dry nitrogen line was routed into the lower plenum. Pressure transducers measured the absolute pressure of the lower plenum, which was compared to absolute chamber pressure to monitor the pressure gradient across the regolith simulant. A precision mass flow controller enabled fine control of plenum pressure. The plenum pressure was quasi-statically increased to set points informed by previous work [1] before shutting off the gas supply to allow the now pressurized gas volume to dissipate through the simulant into the chamber in what was referred to as a ‘leak test.’ In a separate set of tests, the pressure was allowed to quasi-statically build until visual disturbances in the soil were observed via remote camera feeds. Initial leak rate tests were repeated on the disturbed soil to observe changes in leak rates between initial as-poured and post-bubbling soil conditions. This test schedule was applied to simulant columns of three different diameters (approximately 2, 4, and 6 inches) as well as two different heights (around 8 and 16 inches) and was repeated at least three times for each configuration. Plenums were also pressurized while the butterfly valves were opened slightly further to facilitate slow, controlled, constant flowrate draining of the simulant to observe the impact that continuous-feed systems might have on leak rate through the regolith column. Finally, the regolith columns were reset and refilled and static leak rate tests were repeated approximately every decade in ambient pressure from 10e-5 Torr up to ambient conditions to observe the impact of transitioning out of the molecular flow regime on regolith column seal performance. Analysis: Trends with respect to each factor examined were constructed and reported alongside commentary on the most likely responsible physical phenomena. Confidence intervals were applied to the predictive findings of the data by quantifying instrumentation error and accounting for probabilistic impacts of curve fitting the generic exponential function onto pressure decay data. The leak rate is assumed directly proportional to the exponential term coefficient, while the steady-state static pressure gradient is assumed to be the asymptotic limit of each curve fit. Finally, a generic model developed by Ogino et al. [2] and later tuned for ambient testing was re-tuned using vacuum data to estimate leak rates on the lunar surface. Results: Pressure decay rates decreased up to 27-fold under vacuum when compared to ambient pressure, as shown in Figure 2. Simulant column diameter was found to exponentially correlate to both regolith column leak rate and steady-state equilibrium pressure delta. Height was found to be linearly related to leak rate –doubling the regolith column height consistently yielded around a 40% reduction in leak rate. Dynamic regolith draining was found to slightly increase leak rate, although to a lower degree than expected. Conclusion: Because of the transition into molecular gas flow through regolith, vacuum testing is required to properly estimate performance. Regolith columns should be designed as narrow as possible to avoid bridging and minimize leak rate, and as tall as is feasible to both further reduce leak rate and to help prevent bubbling, which has been shown to degrade the quality of the regolith seal. Disrupted soil was found to exhibit higher leak rates than pristine, as-poured soil, indicating the need to consider the impact of recovery from depressurization events. Regolith columns appear equally feasible for both dosed and continuous-feed reactors.

ISRU↗

DOE OSTI2026

A bioinspired approach for adaptive solid-solid phase change material coatings with optimized surface features for passive thermal regulation

The necessity to reduce global energy consumption calls for innovative strategies in building thermal management. Passive thermal regulation, particularly through bio-inspired designs, offers a promising avenue by mimicking nature's efficient control of optical properties. This research introduces a novel, climate-responsive coating that integrates optimized bio-inspired surface features with a solid-solid phase change material (SS-PCM) to dynamically manage solar absorptivity without adding additional thickness, enabling both heating and cooling as needed. Drawing on the photonic architectures of the Saharan silver ant and Morpho Didius butterfly, we employed a modeling and multi-objective optimization framework to tailor these surface features. Simulations reveal that surface texture, rather than the intrinsic phase transition of the SS PCM, dominates optical control. Relative to a flat SS PCM coating, optimized isotropic random roughness and broader range features yielded the highest passive heating power increase of about 144 % and 319 % respectively suitable for cold climates. Saharan ant-inspired features enhanced passive cooling for hot climates, achieving a 21.8 % improvement. For moderate climates, Butterfly-wing-inspired surface features provided a balanced enhancement of 19 % for heating and 7 % for cooling. Across all cases, the optimized surface features reduced combined heating and cooling energy demand more effectively than the baseline coating, while preserving material thickness. These findings demonstrate that climate-adaptive, optimized bio-inspired surface features can unlock the full potential of SS PCM coatings, providing a versatile pathway to significant energy savings in buildings and other applications. The methodology establishes a framework for designing next-generation adaptive envelopes that leverage natural photonic principles for high-impact, low-cost thermal regulation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Inhomogeneity in electronic phase and flat band in magnetic kagome metal Co 3 Sn 2 S 2

Co 3 Sn 2 S 2 has been reported to be a Weyl semimetal with c-axis ferromagnetism below a Curie temperature of 177 K. Despite the large interest in Co 3 Sn 2 S 2 , the magnetic structure is still unclear. Recent studies have challenged the magnetic phase diagram of Co 3 Sn 2 S 2 by reporting unusual magnetic phases including the presence of exchange bias. Here we show, using X-ray Magnetic Circular Dichroism, a shift in the magnetization hysteresis loop, reminiscent of exchange bias and establish that the magnetic moment in Co arises from the spin, with negligible orbital moment. At 6 K, using spatially-resolved angle-resolved photoemission spectroscopy, we detect a butterfly-shaped electronic band structure at small regions of the sample distinct from the known ferromagnetic band structure. Our density functional theory calculations suggest that the butterfly bands correspond to an antiferromagnetic phase. Separately, we detect a sharp flat band at the Fermi level at some regions in the sample, which we attribute to a surface state. These different electronic states found in a stoichiometric intermetallic invite further efforts to explore the origin and nature of the electronic inhomogeneity associated to magnetism on the mesoscale.

Electronic properties and materials↗

NASA NTRS1973

Electron pitch angle distributions throughout the magnetosphere as observed on Ogo 5.

A survey of the equatorial pitch angle distributions of energetic electrons is provided for all local times out to radial distances of 20 earth radii on the night side of the earth and to the magnetopause on the day side of the earth. In much of the inner magnetosphere and in the outer magnetosphere on the day side of the earth, the normal loss cone distribution prevails. The effects of drift shell splitting - i.e., the appearance of pitch angle distributions with minimums at 90 deg, called butterfly distributions - become apparent in the early afternoon magnetosphere at extended distances, and the distribution is observed in to 5.5 earth radii in the nighttime magnetosphere. Inside about 9 earth radii the pitch angle effects are quite energy-dependent. Beyond about 9 earth radii in the premidnight magnetosphere during quiet times the butterfly distribution is often observed. It is shown that these electrons cannot survive a drift to dawn without being considerably modified. The role of substorm activity in modifying these distributions is identified.

West, H. I., Jr.↗

NASA NTRS2001

Compressor Stall Recovery Through Tip Injection Assessed

Aerodynamic stability is a fundamental limit in the compressor design process. The development of robust techniques for increasing stability has several benefits: enabling higher loading and fewer blades, increasing safety throughout a mission, increasing tolerance to stage mismatch during part-speed operation and speed transients, and providing an opportunity to match stages at the compressor maximum efficiency point, thus reducing fuel burn. Mass injection upstream of the tip of a high-speed axial compressor rotor is a stability enhancement approach known to be effective in suppressing stall in tip-critical rotors if the injection is activated before stall occurs. This approach to stall suppression requires that a reliable stall warning system be available. Tests have recently been performed to assess whether steady injection can also be used to recover from fully developed stall. If mass injection is effective in recovering from stall quickly enough to avoid structural damage or loss of engine power, then a stall warning system may not be required. The stall recovery tests were performed on a transonic compressor rotor at its design tip speed of 1475 ft/sec using four injectors evenly spaced around the compressor case upstream of the rotor. The injectors were connected to an external air source. In an actual engine application, the injected air would be supplied with compressor bleed air. The injectors were isolated from the air source by a fast-acting butterfly valve. With the injectors turned off, the compressor was throttled into stall. Air injection was then activated with no change in throttle setting by opening the butterfly valve. The compressor recovered from stall at a fixed throttle setting with the aid of tip injection. The unsteady operating characteristic of the rotor was measured during these tests using high-response pressure sensors located upstream and downstream of the rotor. The figure shows the results, where the unsteady pressure and mass flow are superimposed on the steady operating characteristic. The total injected mass flow was equal to 1.3 percent of the compressor flow. The solid line with no solid squares on it denotes the operating point during the beginning of throttle closure and the initial drop into stall. The gray traces denote the operating point during an additional throttle closure that occurred over the next 1200 rotor revolutions (4 sec). The dashed line denotes the recovery from stall that occurred during 90 rotor revolutions (0.3 sec) after the injectors were activated with no change in throttle setting. Tip injection not only recovers the compressor from stall, but also restores the compressor to its pre-stall level of pressure rise. In contrast, standard stall recovery schemes such as compressor bleed, stator vane actuation, or engine throttle modulation result in a loss of pressure rise across the compressor, which results in a loss of engine power.

Suder, Ken L.↗

DOE OSTI2024

Comparative Study of Neutral and Cationic Sn 2 H 2 : Toward Laboratory Detection of the Cation

Group 14 M 2 H 2 isomers (M = Si, Ge, Sn, and Pb) have attracted interest due to their radically differing electronic structures from acetylene. To better understand the Sn–H interactions of the neutral and cationic Sn 2 H 2 structures, we present the most rigorous study of these systems to date. CCSD(T)/cc-pwCVTZ harmonic frequencies are presented as the first predictions for the neutral and cationic species to date. CCSDT(Q)/CBS relative energies are reported using the focal point approach, confirming the butterfly isomer as the global minimum on the potential energy surface for both the neutral and cationic species. In all, there exist 7 minima and 15 transition states. NBO analysis is also performed to elucidate the changes in bond order going from neutral to cation across all isomers of Sn 2 H 2 . Our results provide insights into the important Sn–H interaction and provide guidance for future work that may detect Sn 2 H$^{+}_{2}$ in the laboratory for the first time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Shape-persistent ladder molecules exhibit nanogap-independent conductance in single-molecule junctions

Molecular electronic devices require precise control over the flow of current in single molecules. However, the electron transport properties of single molecules critically depend on dynamic molecular conformations in nanoscale junctions. Here, in this work, we report a unique strategy for controlling molecular conductance using shape-persistent molecules. Chemically diverse, charged ladder molecules, synthesized via a one-pot multicomponent ladderization strategy, show a molecular conductance (d[log( G/G 0 )]/d x ≈ -0.1 nm -1 ) that is nearly independent of junction displacement, in stark contrast to the nanogap-dependent conductance (d[log( G/G 0 )]/d x ≈ -7 nm -1 ) observed for non-ladder analogues. Ladder molecules show an unusually narrow distribution of molecular conductance during dynamic junction displacement, which is attributed to the shape-persistent backbone and restricted rotation of terminal anchor groups. These principles are further extended to a butterfly-like molecule, thereby demonstrating the strategy's generality for achieving gap-independent conductance. Overall, our work provides important avenues for controlling molecular conductance using shape-persistent molecules. Achieving robust and controllable conductance in single-molecule junctions is challenging due to the dynamic nature of molecular conformations that fluctuate over operational timescales. A strategy using shape-persistent molecules has now been developed that demonstrates nearly junction-displacement-independent conductance, providing a stable solution for single-molecule electronic properties.

molecular electronics↗

DOE OSTI2025

Structural switching dynamically controls the doubly pseudoknotted Rous sarcoma virus–programmed ribosomal frameshifting element

A hallmark of retrovirus replication is the translation of two different polyproteins from one RNA through programmed –1 frameshifting. This is a mechanism in which the actively translating ribosome is induced to slip in the 5′ direction at a defined codon and then continues translating in the new reading frame. Programmed frameshifting controls the stoichiometry of viral proteins and is therefore under stringent evolutionary selection. Forty years ago, the first frameshifting stimulatory element was discovered in the Rous sarcoma virus. The ~120 nt RNA segment was predicted to contain a pseudoknot, but its 3D structure has remained elusive. Now, we have determined cryoEM and X-ray crystallographic structures of this classic retroviral element, finding that it adopts a butterfly-like double-pseudoknot fold. One “wing” contains a dynamic pyrimidine-rich helix, observed crystallographically in two conformations and in a third conformation via cryoEM. The other wing encompasses the predicted pseudoknot, which interacts with a second unexpected pseudoknot through a toggle residue, A2546. This key purine switches conformations between structural states and tunes the stability of interacting residues in the two wings. We find that its mutation can modulate frameshifting by as much as 50-fold, likely by altering the relative abundance of different structural states in the conformational ensemble of the RNA. Taken together, our structure–function analyses reveal how a dynamic double pseudoknot junction stimulates frameshifting by taking advantage of conformational heterogeneity, supporting a multistate model in which high Shannon entropy enhances frameshifting efficiency.

Science & Technology - Other Topics↗

DOE OSTI2024

Evolution of the SLATE linear algebra library

SLATE (Software for Linear Algebra Targeting Exascale) is a distributed, dense linear algebra library targeting both CPU-only and GPU-accelerated systems, developed over the course of the Exascale Computing Project (ECP). While it began with several documents setting out its initial design, significant design changes occurred throughout its development. In some cases, these were anticipated: an early version used a simple consistency flag that was later replaced with a full-featured consistency protocol. In other cases, performance limitations and software and hardware changes prompted a redesign. Sequential communication tasks were parallelized; host-to-host MPI calls were replaced with GPU device-to-device MPI calls; more advanced algorithms such as Communication Avoiding LU and the Random Butterfly Transform (RBT) were introduced. Early choices that turned out to be cumbersome, error prone, or inflexible have been replaced with simpler, more intuitive, or more flexible designs. Applications have been a driving force, prompting a lighter weight queue class, nonuniform tile sizes, and more flexible MPI process grids. Of paramount importance has been building a portable library that works across several different GPU architectures – AMD, Intel, and NVIDIA – while keeping a clean and maintainable codebase. Here we explore the evolving design choices and their effects, both in terms of performance and software sustainability.

Gates, Mark↗

DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Angle-resolved polarized Raman spectroscopy study of phosphorene nanoribbons

We present a systematic angle-resolved polarized Raman spectroscopy (ARPRS) study of black phosphorus (BP) nanostructures formed via electrochemical sodium- and lithium-ion intercalation. Sodium intercalation leads to bundles of densely packed, highly uniform phosphorene nanoribbons (PNRs) separated by parallel amorphous channels, whereas lithium intercalation results in shorter, irregular nanoribbon-like segments with lower aspect ratios. In both cases, six additional Raman peaks (P1–P6) appear alongside the three primary Raman-active modes of BP (A$^{1}_{g}$, B 2g , and A$^{2}_{g}$). These peaks are attributed to the amorphous regions, as confirmed by their isotropic angular dependence in ARPRS measurements. The three BP modes show pronounced angular variations that differ significantly between the two intercalated samples. In sodium-intercalated BP, A$^{1}_{g}$ and A$^{2}_{g}$ modes retain a dumbbell-like angular dependence under parallel polarization with enhanced anisotropy and reduced symmetry under crossed polarization. At the same time, the B 2g mode transitions from four-lobed (cloverleaf) polar plot to a butterfly-like one. In contrast, lithium-intercalated BP exhibits weaker anisotropy and less distinct angular polar plots for all three modes. These differences reflect the sensitivity of phonon behavior to underlying nanostructure morphology. The vibrational frequencies density of states (FDOS) calculations attribute the B 2g mode transformation to phonon band folding and mode mixing in PNRs. This study demonstrates the power of ARPRS in probing phonon-structure relationships and highlights the influence of edge geometry and quantum confinement on phonon dispersion in PNRs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗