Unlocking massively parallel spectral proper orthogonal decompositions in the PySPOD package
Not Available
SEARCH · Engineering Papers
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.
Not Available
Unusual successive order in NdB 4 with the frustrated Shastry-Sutherland lattice was investigated by means of inelastic neutron scattering. Here, we found that the pseudoquartet with |±5/2⟩ and |±7/2⟩ in the paramagnetic state splits into four levels in the magnetically ordered states. The observed spectra are in good agreement with the 𝑓-electron states obtained from the on-site Hamiltonian with the magnetic interaction ⟨𝐽 𝑥 ⟩ 𝐽 𝑥 , the coupling of 𝐽 𝑥 with quadrupole 𝑂 𝑧𝑥 , and the crystalline electric field potential. The 𝑓-electron states reproduce the temperature dependences of the in-plane moment 𝑚 ⊥ = ⟨𝐽 𝑥 ⟩, which is the primary order parameter, and a weakly induced 𝑚 𝑐 = ⟨𝐽 𝑧 ⟩ as the secondary order parameter. The magnetic-quadrupole coupling emerges and plays an essential role under the geometrical frustration, where the Heisenberg interactions are canceled out.
Performing a variety of numerical computations efficiently and, at the same time, in a portable fashion requires both an overarching design followed by a number of implementation strategies. All of these are exemplified below as we present transitioning the PLASMA numerical library from relying on dependence-driven large tasks to achieving utilization of fine grain tasking and offload to hardware accelerators while keeping its core dependence sets: OpenMP source code pragmas and runtime for most system-level functionality and basic low-level numerical kernels provided directly by hardware vendors or open source projects with vendor contributions. We also present new algorithmic methods and their efficient parallel implementations including fine grained tasking for eigen-spectrum slicing and offload for mixed-precision eigenvalue refinement. We provide performance, scaling, and numerical results showing sizable gains over the available solutions from either the open source and vendor-provided packages.
Explore the source record for details and available documents.
Systems theoretic process analysis (STPA) is becoming an increasingly popular technique to assess how complex digital software systems can fail. Rather than defining failures by their observable failure events, which may be sparse especially for safety rated nuclear digital instrumentation and control systems (DI&C), failures are defined as postulated unsafe actions under specific contextual conditions. This permits a top-down analysis of system hazards and identifies whether imposed constraints and requirements can sufficiently address undesirable hazards. However, STPA is a qualitative approach at identifying inadequacies in the development process and cannot currently be used to quantify unsafe action likelihoods for probabilistic risk assessment. Therefore, in this work, we examine the root causes of software failure and explore whether a consistent correlation can be linked to specific unsafe action classes. We implement Lbl2Vec, an unsupervised document classification and retrieval algorithm, on a database of 4,096 software defect reports acquired from various open-source software systems. By analyzing sentence structure, embedded labels, and word vectors, we show that certain defect types positively correlate to specific unsafe action classes over others. The correlations developed can be used to estimate the failure probability of safety intended DI&C systems which provides a licensing basis for nuclear plant modernization efforts.
Systems theoretic process analysis (STPA) is becoming an increasingly popular technique to assess how complex digital software systems can fail. Rather than defining failures by their observable failure events, which may be sparse especially for safety rated nuclear digital instrumentation and control systems (DI&C), failures are defined as postulated unsafe actions under specific contextual conditions. This permits a top-down analysis of system hazards and identifies whether imposed constraints and requirements can sufficiently address undesirable hazards. However, STPA is a qualitative approach at identifying inadequacies in the development process and cannot currently be used to quantify unsafe action likelihoods for probabilistic risk assessment. Therefore, in this work, we examine the root causes of software failure and explore whether a consistent correlation can be linked to specific unsafe action classes. We implement Lbl2Vec, an unsupervised document classification and retrieval algorithm, on a database of 4,096 software defect reports acquired from various open-source software systems. By analyzing sentence structure, embedded labels, and word vectors, we show that certain defect types positively correlate to specific unsafe action classes over others. The correlations developed can be used to estimate the failure probability of safety intended DI&C systems which provides a licensing basis for nuclear plant modernization efforts.
Efficiently learning expectation values of a quantum state using classical shadow tomography has become a fundamental task in quantum information theory. In a classical shadows protocol, one measures a state in a chosen basis $\mathcal{W}$ after it has evolved under a unitary transformation randomly sampled from a chosen distribution $\mathcal{U}$. In this work we study the case where $\mathcal{U}$ corresponds to either local or global orthogonal Clifford gates, and $\mathcal{W}$ consists of real-valued vectors. Our results show that for various situations of interest, this ‘real’ classical shadow protocol improves the sample complexity over the standard scheme based on general Clifford unitaries. For example, when one is interested in estimating the expectation values of arbitrary real-valued observables, global orthogonal Cliffords typically decrease the required number of samples by a factor of two. More dramatically, for k-local observables composed only of real-valued Pauli operators, sampling local orthogonal Cliffords leads to a reduction by an exponential-in-k factor in the sample complexity over local unitary Cliffords. Finally, we show that by measuring in a basis containing complex-valued vectors, orthogonal shadows can, in the limit of large system size, exactly reproduce the original unitary shadows protocol.
A combination of experiments and optical modeling provided insight into the mechanism of mesoscale woodpile formation in response to an orthogonal shift in polarization during photoelectrochemical deposition of Se–Te. Cathodic deposition of semiconducting Se–Te using spatially uniform, linearly polarized illumination produced arrays of lamellae that were aligned parallel to the optical E-field oscillation. Continued deposition in conjunction with an orthogonal shift in the polarization direction then produced aligned bridging features that spanned the void space between, and were orthogonal to, the preexisting lamellae. The height and pitch, respectively, in each layer of the woodpile were a function of the charge density and illumination wavelength during deposition. A Monte Carlo model, in which material addition was scaled by the absorption magnitude obtained from electromagnetic simulations, produced morphologies that were nominally identical to those observed experimentally. Here, the formation of mesoscale woodpiles is consistent with a mechanism that involves a series of spontaneously initiated, concerted light–matter interactions during the photoelectrochemical deposition process.
Distortion generated in machined, monolithic, thin-walled aerospace components due to residual stresses leads to significant material and economic waste in the manufacturing industry. Inherent residual stress (IRS) present in stock materials combines with machining-induced residual stress (MIRS) to influence the final machined part distortion. It is hypothesized that MIRS can be controlled, based on the part geometry, through deliberate cutting tool geometry and process parameter modifications to negate the effect of IRS on distortion, consequently resulting in distortion-free parts. A finite element (FE) orthogonal cutting model is developed to study how tool geometry and process parameters influence machining-induced residual stress (MIRS). Orthogonal cutting experiments are performed on Al 7075-T651 samples to measure cutting forces and MIRS. A cutting force dynamometer is used to measure forces during cutting and a novel digital image correlation (DIC) based hole drilling technique is employed to measure the near-surface residual stress (RS) in the cut samples. These data are subsequently utilized to validate the FE prediction model. Various cases of cutting simulations involving different depths of cut, tool tip radii, and rake angles are performed to study their effects on RS. Similar to prior literature, increasing the depth of cut, tool tip radius, or rake angle is found to promote the formation of near-surface tensile stresses. The competing effects of material plowing and temperature are shown to determine the type of RS at the end of the cut. Moreover, a window of variation of RS (up to ± 400 MPa) is estimated within the given range of conditions, allowing for the control of MIRS through tool and process modification.
Abstract Cosmic microwave background (CMB) experiments have deployed focal planes with $$\mathcal {O}(10^{4})$$ O ( 10 4 ) transition edge sensor (TES) bolometers cooled to sub-Kelvin temperatures by multiplexing the readout of many TES channels onto a single pair of wires. Digital Frequency-domain Multiplexing (DfMux) is a multiplexing technique used in many CMB polarization experiments, such as the Simons Array, SPT-3 G, and EBEX. The DfMux system studied here uses LC filters with resonant frequencies ranging from 1.5 to 4.5 MHz connected to an array of TESs. Each detector has an amplitude-modulated carrier tone at the resonant frequency of its accompanying LC resonator. The signal is recovered via quadrature demodulation where the in-phase (I) component of the demodulated current is in phase with the complex admittance of the circuit and the quadrature (Q) component is orthogonal to I. Observed excess current noise in the Q component is consistent with fluctuations in the resonant frequency. This noise has been shown to be non-orthogonal to the phase of the detector’s responsivity. We present a detailed analysis of the phase of responsivity of the TES and noise sources in our DfMux readout system. Further, we investigate how modifications to the TES operating resistance and bias frequency can affect the phase of noise relative to the phase of the detector responsivity, using data from Simons Array to evaluate our predictions. We find that both the phase of responsivity and phase of noise are functions of the two tuning parameters, which can be purposefully selected to maximize signal-to-noise (SNR) ratio.
We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.
The ability to uniquely identify a compound requires highly precise and orthogonal measurements. Here we describe a newly developed analytical platform that integrates high resolution ion mobility and cryogenic vibrational ion spectroscopy for high-precision structural characterizations. This platform allows for the temporal separation of isomeric/isobaric ions and provides a highly sensitive description of the ion’s adopted geometry in the gas phase. The combination of these orthogonal structural measurements yields precise descriptors that can be used to resolve between and confidently identify highly similar ions. The unique benefits of our instrument, which integrates a structures for lossless ion manipulations ion mobility (SLIM IM) device with messenger tagging infrared spectroscopy, include the ability to perform high-resolution ion mobility separations and to record the IR spectra of all ions simultaneously. The SLIM IM device, with its 13 m separation path length, allows for multipass experiments to be performed for increased resolution as needed. It is integrated with an Agilent qTOF MS where the collision cell was replaced with a cryogenically held (30 K) TW-SLIM module. The cryo-SLIM is operated in a novel manner that allows ions to be streamed through the device and collisionally cooled to a temperature where they can form noncovalently bound N 2 complexes that are maintained as they exit the device and are detected by the TOF mass analyzer. The instrument can be operated in two modes: IMS+IR where the IR spectra for mobility-selected ions can be recorded and IR-only mode where the IR spectra for all mass-resolved ions can be recorded. In IR-only mode, IR spectra (400 cm –1 spectral range) can be recorded in as short as 2 s for high throughput measurements. Further, this work details the construction of the instrument and modes of operation. It provides initial benchmarking of CCS and IR measurements to demonstrate the utility of this instrument for targeted and untargeted approaches.
Circularly polarized (CP) vibrations are the nuclear-motion analog of CP light, emerging in molecules with non-Abelian point-group symmetry that support orthogonal and degenerate vibrational normal modes. Here, in this work, we theoretically explore the optical addressability of CP vibrations, motivated by their potential as nanoscale angular momentum states suitable for information storage and manipulation. We investigate how symmetry-breaking chemical modifications affect a molecule’s capacity to support CP vibrations. Notably, we find cases where modes retain their orthogonality and degeneracy in spite of such modifications, broadening the opportunities afforded by CP vibrations beyond non-Abelian point-group molecules. Prospects for experimental implementations using CP pump–probe techniques are briefly discussed.
This study provides an in-depth investigation of the primary decomposition of polyethylene terephthalate (PET) via pyrolysis, employing an experimental-analytic workflow that integrates design of experiments (DoE), micropyrolysis coupled with comprehensive two-dimensional gas chromatography (GC×GC), and multivariate data analysis to verify intrinsic kinetic conditions and elucidate evolving product distributions for mapping key reaction pathways. Peaks that could not be identified using commercial spectral libraries were assigned using Mass Frontier simulations, enabling the identification of divinyl terephthalate, ethyl vinyl terephthalate, and 2-(benzoyloxy)ethyl vinyl terephthalate. A polar×polar (non-orthogonal) column set tailored for the detection of carboxylic acids enhanced the quantification of benzoic acid, 4-vinylbenzoic acid, 4-ethylbenzoic acid, and methylbenzoic acid by up to 6-fold relative to an orthogonal column combination (non-polar×mid-polar). Moreover, pyrolysis variables were systematically evaluated using a Box- Behnken design (BBD), encompassing pyrolysis temperature (500−600 °C), sample weight (50−150 μg), and carrier gas flow rate (100−300 mL min −1 ). Among these, pyrolysis temperature was the only statistically significant factor influencing product yields, ranging from 58.78 to 84.26 wt %. In contrast, neither the sample weight nor the carrier gas flow rate had a significant effect on product yields within the evaluated experimental space. At 600 °C, the major pyrolysis products were benzoic acid (up to 20.20 ± 1.46 wt %) and CO 2 (up to 21.28 ± 1.46 wt %), which can be produced through decarboxylation reactions. These findings underscore the critical importance of selecting appropriate analytical columns for the accurate quantification of heteroatomcontaining products such as carboxylic acids, which may otherwise be underestimated or undetected due to their reactivity with the stationary phase of non-polar and mid-polar columns, as well as other GC components. They also highlight the importance of selecting pyrolysis conditions for investigating the primary decomposition of PET under an isothermal kinetically limited regime.
Efficient depolymerization of polymers with all-carbon backbones under mild conditions would be valuable in chemically recycling commodity plastics. Poly(methyl methacrylate) (PMMA) is a commodity thermoplastic that is currently depolymerized under temperatures in excess of 400 °C. Herein, we lower the temperatures needed to achieve depolymerization of PMMA by performing radical generation and depropagation with orthogonal stimuli. This first demonstration of electrochemically initiated PMMA depolymerization relies on reduction of phthalimide esters that, upon subsequent decarboxylation, generate polymer-centered radicals. These radicals then spontaneously unzip the polymer back to its monomeric constituents at temperatures as low as 105 °C. We studied the mechanism and efficiency of this transformation as a function of phthalimide ester placement, incorporation density, and polymer molecular weight. We found that chain-end activation is effective for modest molecular weights but suffers diminished efficiency at higher degrees of polymerization. In contrast, pendent-group activation is more effective for depolymerizing higher molecular weight species. Integrating higher molar amounts of phthalimide ester pendants leads to more effective depolymerization, with >95% depolymerization in copolymers with 5 mol% phthalimide ester incorporation. We leveraged this understanding to create a custom electro-distillation apparatus that allowed us to simultaneously electrochemically depolymerize PMMA and directly distill methyl methacrylate in >22% yield, which could ultimately be repolymerized. These findings establish electrochemistry as a versatile and orthogonal stimulus for vinyl polymer depolymerization and provide a foundation for closed-loop electrochemical recycling of widely used plastics.
Chiral organic-inorganic hybrid perovskites offer a promising platform for developing non-linear chiro-optical applications and chiral-induced spin selectivity. Here, we show that achiral hybrid perovskites that have highly ordered ferroelectric domains with orthogonal polarization exhibit planar chirality, as manifested by second harmonic generation with strong circular dichroism. Interestingly, the handedness of the second harmonic generation circular dichroism response can be alternatingly switched between orthogonally polarized domains and domain walls. To correlate the origin of planar chirality in these domains, atom-resolved atomic force microscopy is used to reveal distinct arrangements of atomic rows in these ferroelectric crystals that indicate planar symmetry breaking. Remarkably, a large SHG anisotropy factor of 0.41 is measured at the domain wall. Additionally, spatially resolved two-photon photoluminescence excitation measurement provide evidence of larger Rashba spin-splitting in these chiral ferroelectric domains compared to domain-free areas. Our discoveries of domain-dependent planar chirality and spin texture hold great promise for advancement of domain-specific chiro-optical applications.