Search NASA⌕ Search

SEARCH · Search NASA

Results for “data parallel”

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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Hyper-Rayleigh scattering in centrosymmetric systems

Hyper-Rayleigh scattering (HRS) is an incoherent mechanism for optical second harmonic generation. The frequency-doubled light that emerges from this mechanism is not emitted in a laser-like manner, in the forward direction; it is scattered in all directions. The underlying theory for this effect involves terms that are quadratic in the incident field and involves an even-order optical susceptibility (for a molecule, its associated hyperpolarizability). In consequence, HRS is often regarded as formally forbidden in centrosymmetric media. However, for the fundamental three-photon interaction, theory based on the standard electric dipole approximation, representable as E1{sup 3}, does not account for all experimental observations. The relevant results emerge upon extending the theory to include E1{sup 2}M1 and E1{sup 2}E2 contributions, incorporating one magnetic dipolar or electric quadrupolar interaction, respectively, to a consistent level of multipolar expansion. Both additional interactions require the deployment of higher orders in the multipole expansion, with the E1{sup 2}E2 interaction analogous in rank and parity to a four-wave susceptibility. To elicit the correct form of response from fluid or disordered media invites a tensor representation which does not oversimplify the molecular components, yet which can produce results to facilitate the interpretation of experimental observations. The detailed derivation in this work leads to results which are summarized for the following: perpendicular detection of polarization components both parallel and perpendicular to the pump radiation, leading to distinct polarization ratio results, as well as a reversal ratio for forward scattered circular polarizations. The results provide a route to handling data with direct physical interpretation, to enable the more sophisticated design of molecules with sought nonlinear optical properties.

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

Anomalous resistivity and electron heating by lower hybrid drift waves inside reconnecting current sheets

Inside an electron diffusion region of laboratory reconnection experiments, the quasi-electrostatic lower hybrid drift wave (ES-LHDW) is observed when a significant guide field component is present. Through direct measurement of the anomalous drag term and quasilinear analysis, it is shown that ES-LHDW can account for approximately 20% of the mean reconnection electric field in a case with moderate guide field. This value exceeds the contribution from classical resistivity, which is around 10%. The effects of the Lorentz force term, often neglected for electrostatic waves, are crucial for the observed correlation between electric field and density fluctuations. Anomalous electron heating by the perturbed current and resistivity (2.6 MW/m 3 ) also surpasses the classical Ohmic heating, which is about 2.0 MW/m 3 . For the case with a high guide field, significantly higher local electron temperatures were observed during periods of strong ES-LHDW activity. A statistical analysis further supports electron heating by LHDW, showing a larger increase in electron temperature with a high guide field. Finally, data from the Magnetospheric Multiscale mission provide evidence of Landau damping of ES-LHDW, suggesting that ES-LHDW may contribute to the generation of nonthermal electrons along the direction parallel to the magnetic field.

Magnetic reconnection↗
DOE OSTI2014

Alignments of the galaxies in and around the Virgo cluster with the local velocity shear

Observational evidence is presented for the alignment between the cosmic sheet and the principal axis of the velocity shear field at the position of the Virgo cluster. The galaxies in and around the Virgo cluster from the Extended Virgo Cluster Catalog that was recently constructed by Kim et al. are used to determine the direction of the local sheet. The peculiar velocity field reconstructed from the Sloan Digital Sky Survey Data Release 7 is analyzed to estimate the local velocity shear tensor at the Virgo center. Showing first that the minor principal axis of the local velocity shear tensor is almost parallel to the direction of the line of sight, we detect a clear signal of alignment between the positions of the Virgo satellites and the intermediate principal axis of the local velocity shear projected onto the plane of the sky. Furthermore, the dwarf satellites are found to appear more strongly aligned than their normal counterparts, which is interpreted as an indication of the following. (1) The normal satellites and the dwarf satellites fall in the Virgo cluster preferentially along the local filament and the local sheet, respectively. (2) The local filament is aligned with the minor principal axis of the local velocity shear while the local sheet is parallel to the plane spanned by the minor and intermediate principal axes. Our result is consistent with the recent numerical claim that the velocity shear is a good tracer of the cosmic web.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The skeleton of the Milky Way

Recently, Goodman et al. argued that the very long, very thin infrared dark cloud “Nessie” lies directly in the Galactic midplane and runs along the Scutum–Centaurus Arm in position–position–velocity (p–p–v) space as traced by lower-density CO and higher-density NH{sub 3} gas. Nessie was presented as the first “bone” of the Milky Way, an extraordinarily long, thin, high-contrast filament that can be used to map our Galaxy’s “skeleton.” Here we present evidence for additional bones in the Milky Way, arguing that Nessie is not a curiosity but one of several filaments that could potentially trace Galactic structure. Our 10 bone candidates are all long, filamentary, mid-infrared extinction features that lie parallel to, and no more than 20 pc from, the physical Galactic mid-plane. We use CO, N{sub 2}H{sup +}, HCO{sup +}, and NH{sub 3} radial velocity data to establish the three-dimensional location of the candidates in p–p–v space. Of the 10 candidates, 6 also have a projected aspect ratio of ≥50:1; run along, or extremely close to, the Scutum–Centaurus Arm in p–p–v space; and exhibit no abrupt shifts in velocity. The evidence presented here suggests that these candidates mark the locations of significant spiral features, with the bone called filament 5 (“BC-18.88-0.09”) being a close analog to Nessie in the northern sky. As molecular spectral-line and extinction maps cover more of the sky at increasing resolution and sensitivity, it should be possible to find more bones in future studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Toward active disruption avoidance via real-time estimation of the safe operating region and disruption proximity in tokamaks

This paper describes a real-time capable algorithm for identifying the safe operating region around a tokamak operating point. The region is defined by a convex set of linear constraints, from which the distance of a point from a disruptive boundary can be calculated. The disruptivity of points is calculated from an empirical machine learning predictor that generates the likelihood of disruption. While the likelihood generated by such empirical models can be compared to a threshold to trigger a disruption mitigation system, the safe operating region calculation enables active optimization of the operating point to maintain a safe margin from disruptive boundaries. The proposed algorithm is tested using a random forest disruption predictor fit on data from DIII-D. The safe operating region identification algorithm is applied to historical data from DIII-D showing the evolution of disruptive boundaries and the potential impact of optimization of the operating point. Real-time relevant execution times are made possible by parallelizing many of the calculation steps and implementing the algorithm on a graphics processing unit. Lastly, a real-time capable algorithm for optimizing the target operating point within the identified constraints is also proposed and simulated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Measurements and modeling of type-I and type-II ELMs heat flux to the DIII-D divertor

Type-I and type-II edge-localized-modes (ELMs) heat flux profiles measured at the DIII-D divertor feature a peak in the vicinity of the strike-point and a plateau in the scrape-off-layer (SOL), which extends to the first wall. The plateau is present in attached and detached divertors and it is found to originate with plasma bursts upstream in the SOL. The integrated ELM heat flux is distributed at ~65% in the peak and ~35% in this plateau. The parallel loss model, currently used at ITER to predict power loads to the walls, is benchmarked using these results in the primary and secondary divertors with unprecedented constraints using experimental input data for ELM size, radial velocity, energy, electron temperature and density, heat flux footprints and number of filaments. The model can reproduce the experimental near-SOL peak within ~20%, but cannot match the SOL plateau. Employing a two-component approach for the ELM radial velocity, as guided by intermittent data, the full radial heat flux profile can be well matched. The ELM-averaged radial velocity at the separatrix, which explains profile widening, increases from ~0.2 km s –1 in attached to ~0.8 km s –1 in detached scenarios, as the ELM filaments' path becomes electrically disconnected from the sheath at the target. The results presented here indicate filaments fragmentation as a possible mechanism for ELM transport to the far-SOL and provide evidence on the beneficial role of detachment to mitigate ELM flux in the divertor far-SOL. However, these findings imply that wall regions far from the strike points in future machines should be designed to withstand significant heat flux, even for small-ELM regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

WIRE: Resource-efficient Scaling with Online Prediction for DAG-based Workflows

This paper introduces WIRE that manages resources for the DAG-based workflows on IaaS clouds. WIRE predicts and plans resources over the MAPE (Monitor-Analyze-Plan-Execute) loops to: 1) Estimate task performance with online data, 2) Conduct simulations to predict the upcoming loads based on online estimates and workflow DAGs, 3) Apply a resource-steering policy to size cloud instance pools for the maximal parallelism that is consistent with low cost. We implement WIRE on Pegasus WMS/HTCondor and evaluate its performance on the ExoGENI network cloud. The results show that WIRE attains low resource cost with the performance that is typically within a factor of two of optimal.

Xie, Bing↗
DOE OSTI2017

Faster than Real-Time Dynamic Simulation for Large-Size Power System with Detailed Dynamic Models using High-Performance Computing Platform

This paper presents a faster-than-real-time dynamic simulation software package that is designed for large-size power system dynamic simulation. It was developed on the GridPACKTM high-performance computing (HPC) framework. The key features of the developed software package include (1) faster-than-real-time dynamic simulation for a WECC system (17,000 buses) with different types of detailed generator, controller, and relay dynamic models, (2) a decoupled parallel dynamic simulation algorithm with optimized computation architecture to better leverage HPC resources and technologies, (3) options for HPC-based linear and iterative solvers, (4) hidden HPC details, such as data communication and distribution, to enable development centered on mathematical models and algorithms rather than on computational details for power system researchers, and (5) easy integration of new dynamic models and related algorithms into the software package.

Huang, Renke↗
DOE OSTI2021

DETECTING FIRE WITH MACHINE LEARNING-ENABLED VISUAL MONITORING FOR NUCLEAR POWER PLANT ENVIRONMENTS

Nuclear power plants are experiencing significant cost challenges to remain competitive with other energy-generation utilities. Unlike other industries, the cost of operation and maintenance activities is mostly attributed to workforce costs. To mitigate this, nuclear power plant stakeholders are increasingly interested in the development and deployment of machine learning methods to potentially automate or augment manually intensive tasks to reduce costs, especially for monitoring activities. One monitoring function that is visually demanding and that can occur frequently to meet the requirements of a fire protection program is visually monitoring an area for fire occurrence. Currently, fire watch activities consist of a worker physically stationed at a given location with the sole responsibility of observing a given area to ensure a fire is detected and mitigated promptly. This effort focused on the development and evaluation of a suitable deep convolutional neural network to classify individual video frames at a sub-second frequency for the occurrence of “fire” and “no fire” in varying industrial environments similar to nuclear power plants. It is believed that a trained neural network model could be integrated with existing facility video surveillance camera feeds to generate alerts when fire inferences occur in individual frames captured at sub-second temporal resolutions. Extensive effort was dedicated to identifying and curating suitable imagery training data representing varying environments and scene settings with and without flame features to maximize generalization in nuclear power plant environments. The data collection effort resulted in the aggregation of a large, labeled image library exceeding 12,000 images to support model training for diverse industrial environments. A deep neural network model incorporating parallel multi-scale capabilities was developed and trained to support accurate image-based detection of flame incidents of varying sizes and spectral feature properties within heterogeneous scenes. Analysis results show that the trained model can achieve high inference accuracy despite heterogeneous scene environments and components. Testing accuracy exceeded 95.0 percent with very low false positive and false negative inferences.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Raman Lidar MERGE Value-Added Product

The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Climate Research Facility Raman lidars (RLs) are semi-autonomous, land-based, laser remote sensing systems that provide height- and time-resolved measurements of water vapor mixing ratio, temperature, aerosol backscatter, extinction, and linear depolarization ratio from about 200 m to greater than 10 km AGL. These systems transmit at a wavelength of 355 nm with 300 mJ, ~5 ns pulses, and a pulse repetition frequency of 30 Hz. The receiver incorporates nine detection channels, including two water vapor channels at 408 nm, two nitrogen channels at 387 nm, three elastic channels, and two rotational Raman channels for temperature profiling at 354 and 353 nm. Figure 1 illustrates the layout of the ARM RL receiver system. Backscattered light from the atmosphere enters the telescope and is directed into the receiver system (i.e., aft optics). This signal is then split between a narrow-field-of-view radiometer (NFOV) path (blue) and a wide-field-of-view zenith radiometer (WFOV) path (red). The WFOV (2 mrad) path contains three channels (water vapor, nitrogen, and unpolarized elastic), and the NFOV (0.3 mrad) path contains six channels (water vapor, nitrogen, parallel and perpendicular elastic, and two rotational Raman). All nine detection channels use Electron Tubes 9954B photomultiplier tubes (PMTs). The signals from each of the nine PMTs are acquired using transient data recorders from Licel GbR (Berlin, Germany). The Licel data recorders provide simultaneous measurements of both analog photomultiplier current and photon counts at height resolution of 7.5 m and a time resolution of 10 s. The analog signal provides good linearity in the strong signal regime, but poor sensitivity at low signal levels. Conversely, the photo counting signal provides good sensitivity in the weak signal regime, but is strongly nonlinear at higher signal levels. The advantage in recording both signals is that they can be combined (or merged) into a single signal with improved dynamic range. The process of combining the analog and photon counting data has become known as “gluing” (Whiteman et al., 2006).

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

FY26 Progress Report on the Operation of the Activated Materials Laboratory at the Advanced Photon Source as a Nuclear Science User Facilities Partner Facility

The Activated Materials Laboratory (AML), located in the Long Beamline Building of the Advanced Photon Source (APS) at Argonne National Laboratory (ANL), provides a centralized radiological capability for preparing, handling, and supporting synchrotron experiments on activated materials. As a Nuclear Science User Facilities (NSUF) partner facility, the AML enables the receipt of radioactive shipments, open-form sample handling, encapsulation, transport of specimens to and from APS beamlines, and experimentation with dedicated equipment. The partnership includes the APS 1-ID and 20-ID beamlines, which provide high-energy x-ray scattering, tomography, and diffraction microscopy techniques for ex-situ, in-situ, and grainresolved three-dimensional (3D) characterization. All samples supported through the AML and partner beamlines must meet radiological limits of less than 100 mrem/h at 30 cm. This FY26 progress report summarizes the first full year of AML operations as an NSUF partner facility and highlights progress in both user support and capability development. By the end of FY26, the AML had received a total of 10 NSUF-awarded projects, including 2 Consolidated Innovative Nuclear Research (CINR), 7 Rapid Turnaround Experiment (RTE), and 1 Super RTE projects. Beamtime was fully delivered for 3 RTE projects and partially delivered for 1 CINR project, demonstrating successful workflows for receipt, encapsulation, beamline transfer, and radiological experiment execution. These efforts demonstrated safe radiological experiments at APS beamlines for samples with dose rates above the historical 5 mrem/h threshold and now up to 100 mrem/h at 30 cm, marking an important milestone for neutron-irradiated materials research. During FY26, the AML also expanded its experimental capabilities. A Psylotech xTS load frame with in-grip rotation was deployed for room-temperature mechanical testing with threedimensional x-ray characterization during interrupted loading. A Linkam TS1500V vacuum heater was commissioned for thermally driven studies, and temperature calibration experiments were performed to establish specimen-relevant thermal profiles. In parallel, a customized split-tube furnace for high-temperature mechanical testing entered commissioning, and initial work began on robotic sample handling to reduce worker dose and improve operational efficiency. Data workflow and beamline operations continued to mature. Together, these developments demonstrate that the AML is becoming a unique national resource for safe, efficient, and scientifically advanced characterization of activated materials at the APS.

Zhang, Xuan↗
DOE OSTI2022

Development, construction and qualification tests of the mechanical structures of the electromagnetic calorimeter of the Mu2e experiment at Fermilab

The “muon-to-electron conversion” (Mu2e) experiment at Fermilab will search for the Charged Lepton Flavour Violating neutrino-less coherent conversion of a muon into an electron in the field of an aluminum nucleus. The observation of this process would be the unambiguous evidence of physics beyond the Standard Model. Mu2e detectors comprise a straw-tracker, an electromagnetic calorimeter and an external veto for cosmic rays. The calorimeter provides excellent electron identification, complementary information to aid pattern recognition and track reconstruction, and a fast calorimetric online trigger. The detector has been designed as a state-of-the-art crystal calorimeter and employs 1340 pure Cesium Iodide (CsI) crystals readout by UV-extended silicon photosensors and fast front-end and digitization electronics. A design consisting of two identical annular matrices (named “disks”) positioned at the relative distance of 70 cm downstream the aluminum target along the muon beamline satisfies the Mu2e physics requirements.The hostile Mu2e operational conditions, in terms of radiation levels (total ionizing dose of 12 krad and a neutron fluence of 5x1010 n/cm2 @ 1 MeVeq (Si)/y), magnetic field intensity (1 T) and vacuum level (10$^{-4}$ Torr) have posed tight constraints on the design of the detector mechanical structures and materials choice. The support structure of the two 670 crystal matrices employs two aluminum hollow rings and parts made of open-cell vacuum-compatible carbon fiber. The photosensors and service front-end electronics for each crystal are assembled in a unique mechanical unit inserted in a machined copper holder. The 670 units are supported by a machined plate made of vacuum-compatible plastic material. The plate also integrates the cooling system made of a network of copper lines flowing a low temperature radiation-hard fluid and placed in thermal contact with the copper holders to constitute a low resistance thermal bridge. The data acquisition electronics is hosted in aluminum custom crates positioned on the external lateral surface of the two disks. The crates also integrate the electronics cooling system as lines running in parallel to the front-end system.The constraints on the calorimeter mechanical structures design, the development from the conceptual design to the specifications of all the structural components, the status of components production and the components quality assurance tests are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Adiabatic Quantum Computation Applied to Deep Learning Networks

Training deep learning networks is a difficult task due to computational complexity, and this is traditionally handled by simplifying network topology to enable parallel computation on graphical processing units (GPUs). However, the emergence of quantum devices allows reconsideration of complex topologies. We illustrate a particular network topology that can be trained to classify MNIST data (an image dataset of handwritten digits) and neutrino detection data using a restricted form of adiabatic quantum computation known as quantum annealing performed by a D-Wave processor. We provide a brief description of the hardware and how it solves Ising models, how we translate our data into the corresponding Ising models, and how we use available expanded topology options to explore potential performance improvements. Although we focus on the application of quantum annealing in this article, the work discussed here is just one of three approaches we explored as part of a larger project that considers alternative means for training deep learning networks. The other approaches involve using a high performance computing (HPC) environment to automatically find network topologies with good performance and using neuromorphic computing to find a low-power solution for training deep learning networks. Our results show that our quantum approach can find good network parameters in a reasonable time despite increased network topology complexity; that HPC can find good parameters for traditional, simplified network topologies; and that neuromorphic computers can use low power memristive hardware to represent complex topologies and parameters derived from other architecture choices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Regulation of Solar Wind Electron Temperature Anisotropy by Collisions and Instabilities

Abstract Typical solar wind electrons are modeled as being composed of a dense but less energetic thermal “core” population plus a tenuous but energetic “halo” population with varying degrees of temperature anisotropies for both species. In this paper, we seek a fundamental explanation of how these solar wind core and halo electron temperature anisotropies are regulated by combined effects of collisions and instability excitations. The observed solar wind core/halo electron data in ( β ∥ , T ⊥ / T ∥ ) phase space show that their respective occurrence distributions are confined within an area enclosed by outer boundaries. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures and β ∥ is the ratio of parallel thermal energy to background magnetic field energy. While it is known that the boundary on the high- β ∥ side is constrained by the temperature anisotropy-driven plasma instability threshold conditions, the low- β ∥ boundary remains largely unexplained. The present paper provides a baseline explanation for the low- β ∥ boundary based upon the collisional relaxation process. By combining the instability and collisional dynamics it is shown that the observed distribution of the solar wind electrons in the ( β ∥ , T ⊥ / T ∥ ) phase space is adequately explained, both for the “core” and “halo” components.

Yoon, Peter H. (ORCID:0000000181343790)↗
DOE OSTI2016

Progress In Developing An In-Pile Acoustically Telemetered Sensor Infrastructure

A salient grand challenge for a number of Department of Energy programs such as Fuels Cycle Research and Development ( includes Accident Tolerant Fuel research and the Transient Reactor Test Facility Restart experiments), Light Water Sustainability, and Advanced Reactor Technologies is to enhance our fundamental understanding of fuel and materials behavior under irradiation. Robust and accurate in-pile measurements will be instrumental to develop and validate a computationally predictive multi-scale understanding of nuclear fuel and materials. This sensing technology will enable the linking of fundamental micro-structural evolution mechanisms to the macroscopic degradation of fuels and materials. The in situ sensors and measurement systems will monitor local environmental parameters as well as characterize microstructure evolution during irradiation. One of the major road blocks in developing practical robust, and cost effective in-pile sensor systems, are instrument leads. If a wireless telemetry infrastructure can be developed for in-pile use, in-core measurements would become more attractive and effective. Thus to be successful in accomplishing effective in-pile sensing and microstructure characterization an interdisciplinary measurement infrastructure needs to be developed in parallel with key sensing technology. For the discussion in this research, infrastructure is defined as systems, technology, techniques, and algorithms that may be necessary in the delivery of beneficial and robust data from in-pile devices. The architecture of a system’s infrastructure determines how well it operates and how flexible it is to meet future requirements. The limiting path for the effective deployment of the salient sensing technology will not be the sensors themselves but the infrastructure that is necessary to communicate data from in-pile to the outside world in a non-intrusive and reliable manner. This article gives a high level overview of a promising telemetry infrastructure based on acoustic wireless transmission of data that is being developed and tested by the INL, Penn State and Westinghouse.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

BAWLD-CH 4 : a comprehensive dataset of methane fluxes from boreal and arctic ecosystems

Methane (CH 4 ) emissions from the boreal and arctic region are globally significant and highly sensitive to climate change. There is currently a wide range in estimates of high-latitude annual CH 4 fluxes, where estimates based on land cover inventories and empirical CH 4 flux data or process models (bottom-up approaches) generally are greater than atmospheric inversions (top-down approaches). A limitation of bottom-up approaches has been the lack of harmonization between inventories of site-level CH 4 flux data and the land cover classes present in high-latitude spatial datasets. Here we present a comprehensive dataset of small-scale, surface CH 4 flux data from 540 terrestrial sites (wetland and non-wetland) and 1247 aquatic sites (lakes and ponds), compiled from 189 studies. The Boreal–Arctic Wetland and Lake Methane Dataset (BAWLD-CH 4 ) was constructed in parallel with a compatible land cover dataset, sharing the same land cover classes to enable refined bottom-up assessments. BAWLD-CH 4 includes information on site-level CH 4 fluxes but also on study design (measurement method, timing, and frequency) and site characteristics (vegetation, climate, hydrology, soil, and sediment types, permafrost conditions, lake size and depth, and our determination of land cover class). The different land cover classes had distinct CH 4 fluxes, resulting from definitions that were either based on or co-varied with key environmental controls. Fluxes of CH 4 from terrestrial ecosystems were primarily influenced by water table position, soil temperature, and vegetation composition, while CH 4 fluxes from aquatic ecosystems were primarily influenced by water temperature, lake size, and lake genesis. Models could explain more of the between-site variability in CH 4 fluxes for terrestrial than aquatic ecosystems, likely due to both less precise assessments of lake CH 4 fluxes and fewer consistently reported lake site characteristics. Analysis of BAWLD-CH 4 identified both land cover classes and regions within the boreal and arctic domain, where future studies should be focused, alongside methodological approaches. Overall, BAWLD-CH 4 provides a comprehensive dataset of CH 4 emissions from high-latitude ecosystems that are useful for identifying research opportunities, for comparison against new field data, and model parameterization or validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

MultiPEM Toolbox: User Manual [Rev. 2]

This document explains use of the Multi-Phenomenology Explosion Monitoring (Multi PEM) Toolbox, a collection of R scripts for estimating the unknown device parameters of a new event with uncertainty quantification. The methodology and application used for illustration in this user manual are fully documented in a Los Alamos National Laboratory technical report hereafter designated “WPA” for reference. Additional details on the application are found in a recent journal article. Two assessment types are available: rapid and complete. Rapid assessments are conducted in two stages, as described in Section 2. In the first stage, calibration data are used to estimate forward and error model parameters (WPA, §5.1) and (if relevant) errors-in-variables yield values for calibration sources (WPA, §3, Equation (3)). In the second stage, new event data are used to estimate the unknown new event device parameters (WPA, §5.2) with uncertainty quantification. Two options for treating the inferred first stage parameters in second stage Bayesian analysis are available: fixing them at their maximum likelihood estimate (default), or multiple imputation. Multiple imputation involves utilizing several posterior samples (imputations) of the first stage parameters as fixed values in the second stage posterior sampling of the new event device parameters. Second stage sampling is conducted across imputations in parallel to improve computational efficiency. This method produces improved uncertainty quantification of the new event device parameters compared with the default treatment of the first stage parameters, at the expense of additional computation. Complete assessments are conducted in a single stage, as described in Section 3. Calibration and (if relevant) new event data are used simultaneously to estimate all forward model, error model, and (if relevant) new event device parameters with uncertainty quantification on the latter. As the name suggests, rapid assessments generally run substantially faster than complete assessments (even with multiple imputation), because the results of first stage analysis can be stored and incorporated into estimating a relatively low-dimensional space of new event device parameters whenever relevant new event data becomes available. On the other hand, complete assessments must be run on the full set of model and device parameters with calibration and new event data every time the latter becomes available.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Parallel hybrid quantum-classical machine learning for kernelized time-series classification

Supervised time-series classification garners widespread interest because of its applicability throughout a broad application domain including finance, astronomy, biosensors, and many others. Here, in this work, we tackle this problem with hybrid quantum-classical machine learning, deducing pairwise temporal relationships between time-series instances using a timeseries Hamiltonian kernel (TSHK). A TSHK is constructed with a sum of inner products generated by quantum states evolved using a parameterized time evolution operator. This sum is then optimally weighted using techniques derived from multiple kernel learning. Because we treat the kernel weighting step as a differentiable convex optimization problem, our method can be regarded as an end-to-end learnable hybrid quantum-classical-convex neural network, or QCC-net, whose output is a data set-generalized kernel function suitable for use in any kernelized machine learning technique such as the support vector machine (SVM). Using our TSHK as input to a SVM, we classify univariate and multivariate time-series using quantum circuit simulators and demonstrate the efficient parallel deployment of the algorithm to 127-qubit superconducting quantum processors using quantum multi-programming.

97 MATHEMATICS AND COMPUTING↗