SEARCH · Search NASA

Results for “Long Time Integration”

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 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Ka-band linearizer for the Ultra-Compact X-ray free-electron laser at UCLA

There is a strong demand for accelerating structures able to achieve higher gradients and more compact dimensions for the next generation of linear accelerators for research, industrial and medical applications. Notably innovative technologies will permit compact and affordable advanced accelerators as the linear collider and X-ray free-electron lasers (XFELs) with accelerating gradients over twice the value achieved with current technologies. In particular XFELs are able to produce coherent X-ray pulses with peak brightness 10 orders of magnitude greater than preceding approaches, which has revolutionized numerous research fields through imaging of the nanoscopic world at the time and length scale of atom-based systems, that is of femtosecond and Angstrom. There is a strong interest for combining these two fields, to form a proper tool with the goal of producing a very compact XFEL in order to investigate multi-disciplinary topics in chemistry, biology, materials science, medicine and physics. In the framework of the Ultra-Compact XFEL project under study at the University of California, Los Angeles, a high gradient radio-frequency accelerating structure for the longitudinal phase-space linearization with an integrated voltage of at least 15 MV working on 6th harmonic of the main Linac frequency is required. We here present the electromagnetic design of a cryogenic normal-conducting 8 cm long Ka-band standing-wave linearizer working on π mode with a target accelerating gradient beyond 100 MV/m. The studies have been performed analytically and numerically to investigate the beam dynamics and electromagnetic issues.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Differential Impact of Climate Change on the Hydropower Economies of Karakoram and Himalayan River Basins

Water stored in the form of snow and glaciers in the High Mountain Asia (HMA) region regulates the water supply, and resultant water-based economies, that support the livelihoods of millions of people. Trends in the seasonal and long-term melting of snow and glaciers, governed by initial ice reserves, meteorological factors and geographic features, vary across sub-basins in the HMA region. We examined the economic impacts of climate-led changes in river flow in two drainage basins, one each from the Karakoram and Central Himalaya region. We used an integrated assessment framework to estimate the changes in economic value of the hydropower generation from hydropower plants on rivers fed by snow and glacier melt in the two sub-basins. The framework, developed under a NASA High Mountain Asia project, coupled biophysical models (a suite of climate models, snow/glacier-hydrology, and hydropower model) with economic analysis. We compared the differences in estimated river flow over historic and future time using the water balance model in sixteen scenarios (eight climate models and two emissions scenarios) for rivers upstream of hydropower plants in each sub-basin. Using the hydropower model we developed, we estimated the changes in hydropower generation at the Naltar IV hydropower plant, with an 18 MW capacity, located in Hunza, Karakoram, and the Trishuli hydropower plant, with a 19.6 MW capacity, in Trishuli, Central Himalaya. When compared to their baselines, the estimated impact of climate change and temporal variability were higher for the Naltar plant than for the Trishuli plant. Our sensitivity analysis shows that hydropower plants with water storage facilities help reduce the impact of changes, but the estimated impacts are higher for the higher capacity plants. This study provides an example of the differential impacts of climate change on hydropower plants located in rivers fed by varying amounts of snow and glacier melt at different decades in this century. This type of integrated assessment of climate change impact will support the scientific understanding of hydrologic flow and its impacts on a hydropower economy under various climate scenarios, as well as generate information about water resource management in a changing climate.

climate change↗

DOE OSTI2023

Face mask integrated with flexible and wearable manganite oxide respiration sensor

Face masks are key personal protective equipment for reducing exposure to viruses and other environmental hazards such as air pollution. Integrating flexible and wearable sensors into face masks can provide valuable insights into personal and public health. The advantages that a breath-monitoring face mask requires, including multi-functional sensing ability and continuous, long-term dynamic breathing process monitoring, have been underdeveloped to date. Here, we design an effective human breath monitoring face mask based on a flexible La 0.7 Sr 0.3 MnO 3 (LSMO)/Mica respiration sensor. The sensor’s capabilities and systematic measurements are investigated under two application scenes, namely clinical monitoring mode and daily monitoring mode, to monitor, recognise, and analyse different human breath status, i.e., cough, normal breath, and deep breath. This sensing system exhibits super-stability and multi-modal capabilities in continuous and long-time monitoring of the human breath. We determine that during monitoring human breath, thermal diffusion in LSMO is responsible for the change of resistance in flexible LSMO/Mica sensor. Both simulated and experimental results demonstrate good discernibility of the flexible LSMO/Mica sensor operating at different breath status. Our work opens a route for the design of novel flexible and wearable electronic devices.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Reparameterizing Litter Decomposition Using a Simplified Monte Carlo Method Improves Litter Decay Simulated by a Microbial Model and Alters Bioenergy Soil Carbon Estimates

Litter decomposition determines soil organic matter (SOM) formation and plant-available nutrient cycles. Therefore, accurate model representation of litter decomposition is critical to improving soil carbon (C) projections of bioenergy feedstocks. Soil C models that simulate microbial physiology (i.e., microbial models) are new to bioenergy agriculture, and their parameterization is often based on small datasets or manual calibration to reach benchmarks. Here, we reparameterized litter decomposition in a microbial soil C model (CORPSE - Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment) using the continental-scale Long-term Inter-site Decomposition Experiment Team (LIDET) dataset which documents decomposition across a range of litter qualities over a decade. Here, we conducted a simplified Monte Carlo simulation that constrained parameter values to reduce computational costs. The LIDET-derived parameters improved modeled C and nitrogen (N) remaining, decomposition rates, and litter mean residence times as compared to Baseline parameters. We applied the LIDET litter decomposition parameters to a microbial bioenergy model (Fixation and Uptake of Nitrogen – Bioenergy Carbon, Rhizosphere, Organisms, and Protection) to examine soil C estimates generated by Baseline and LIDET parameters. LIDET parameters increased estimated soil C in bioenergy feedstocks, with even greater increases under elevated plant inputs (i.e., by increasing residue, N fertilization). This was due to the integrated effects of plant litter quantity, quality, and agricultural practices (tillage, fertilization). Collectively, we developed a simple framework for using large-scale datasets to inform the parameterization of microbial models that impacts projections of soil C for bioenergy feedstocks.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Mono-Ether and Alcohol Bioblendstocks to Reduce the Fuel Penalty of Mixing Controlled Compression Ignition (MCCI) Engine Aftertreatment

An integrated approach utilizing catalysis experiments, process systems engineering, fuel property modeling, and engine testing was utilized in this project to optimize the production process and composition of a #2 diesel bioblendstock produced from ethanol consisting primarily of long-chain mono-ethers. The primary objective of the project was to determine the composition and to design the production process for a bioblendstock for #2 diesel fuel with > 50% reduction in greenhouse gas emissions relative to conventional diesel fuel. The desired bioblendstock needed to be blendable with #2 diesel fuel at > 5 vol. % while still meeting ASTM D975 diesel fuel specification properties and achieving improvements in fuel properties: increased cetane number, decreased sooting, and reduced pour point and cloud point temperatures. At the same time, it needed to reduce the fuel energy penalty associated with MCCI engine aftertreatment resulting in improved system efficiency. The findings of the current project demonstrate that primary objective of the work has been met, i.e., to determine the composition and to design the production process for a bioblendstock for #2 diesel fuel with > 50% reduction in greenhouse gas emissions relative to conventional diesel fuel. Additionally, the results indicate that the property objectives (increased cetane number, reduced pour and cloud points, and reduced sooting propensity) for the designed bioblendstock composition have also been met. Engine testing performed has also confirmed that the increased reactivity of the bioblendstock can be used to improve catalyst heating operation and to reduce the fuel penalty associated with this operation mode. The fuel property results demonstrate that >5 vol.% blending is easily achieved while meeting the ASTM D975 #2 diesel fuel property specifications tested in this work, as this was achieved for a blend with 43 vol. % of the bioblendstock. The results provide a foundation for future work to scaleup the catalytic production process designed in this work, with many of the challenges and areas for improvement being identified in this work to enable economic production with low GHG lifecycle emissions.

09 BIOMASS FUELS↗

DOE OSTI2023

Real-time evolution of Anderson impurity models via tensor network influence functionals

In this work, we present and analyze two tensor network-based influence functional approaches for simulating the real-time dynamics of quantum impurity models such as the Anderson model. Via comparison with recent numerically exact simulations, we show that such methods accurately capture the long-time nonequilibrium quench dynamics. The two parameters that must be controlled in these tensor network influence functional approaches are a time discretization (Trotter) error and a bond dimension (tensor network truncation) error. We show that the actual numerical uncertainties are controlled by an intricate interplay of these two approximations, which we demonstrate in different regimes. Our work opens the door to using these tensor network influence functional methods as general impurity solvers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Tritium Aging Effects on Fracture Toughness of Stainless Steel Weldments

The long-term embrittlement effects of tritium and decay helium on the structural properties of stainless steels have been studied for years at Savannah River to provide required data for establishing safe operating conditions and lifetimes of pressure vessels used to contain tritium gas. In this study, the fracture toughness properties of the longest-aged tritium-precharged stainless-steel base metals and weldments tested at Savannah River were measured and compared to earlier results. The fracture toughness values were the lowest recorded here for tritium-exposed stainless steel. As-forged and as-welded specimens were thermally precharging with tritium gas at 34.5 MPa and 623 K, then aged for up to 17 years to buildin decay helium prior to testing. ASTM J-Integral fracture mechanics analyses, transmission electron microscopy (TEM), and small angle neutron scattering (SANS) examinations were conducted to characterize the effects of tritium and its radioactive decay product, helium-3. Results show that the fracture toughness values were reduced to less than 2-4% of the as-forged values for specimens with more than 1300 appm helium from tritium decay. The trend of decreasing fracture toughness values with increasing helium content was consistent with earlier observations, and the data show that Type 304L stainless steel is more resistant to tritium-induced cracking than Type 21-6-9 stainless at similar decay helium levels. The fracture toughness properties of long-aged weldments were also affected, but the reductions were not as severe over time because weldments did not retain as much tritium as did the base metals. TEM observations were used to characterize the effects of decay helium bubbles on the deformation substructures, but nanometer-sized helium bubbles were not easily resolved because of high dislocation densities within the forged microstructures. SANS results are presented that suggest the technique can provide information on decay helium bubble size, spacing, and distribution in these steels.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Refueling infrastructure planning in intercity networks considering route choice and travel time delay for mixed fleet of electric and conventional vehicles

The range anxiety has been a major factor that affects the market acceptance of electric vehicles. Even with the recent development of battery technologies, a lack of charging stations and range anxiety are still significant concerns, specifically for intercity trips. This calls for more investments in building charging stations and advancing battery technologies to increase the market share of electric vehicles and improve sustainability. This study suggests a configuration for plug-in electric vehicle charging infrastructure to support long-distance intercity trips of electric vehicles at the network level. A model is proposed to minimize the total system cost including infrastructure investment (building charging stations/spots) and travel time delays (charging time, waiting time in the queue, and detour time to access charging stations). This study fills existing gaps in the literature by capturing realistic patterns of travel demand and considering flow-dependent charging delays at charging stations. Furthermore, the proposed model, which is formulated as a mixed-integer program with nonlinear constraints, solves the optimization problem at the network level. At the network level, impacts of charging station locations on the traffic assignment problem with a mixed fleet of electric and conventional vehicles need to be considered. To this end, a traffic assignment module is integrated with a simulated annealing algorithm. The numerical experiments show a satisfactory application of the model for a full-scale case study (intercity network in Michigan). The solution quality and efficiency of the proposed solution algorithm are evaluated against those of an enumeration approach for a small case study. The results suggest that even for the current market share and charging stations’ setting, a significant investment is needed to support intercity trips without range anxiety issues and with acceptable delays. Additionally, through sensitivity analyses, the required infrastructure and battery investments to support intercity trips with acceptable delays are established for hypothetical increased market shares and battery size in the future.

42 ENGINEERING↗

DOE OSTI2022

The root-zone soil moisture spectrum in a mediterranean ecosystem

Storage of water within soil pores of the root zone introduce memory effects in the dynamics of soil moisture that are considerably longer than the integral timescale of many atmospheric processes. Thus, hydro-climatic states can be “sustained” through land-surface heat and water vapor fluxes primarily because they can “feed off” on this long-term soil moisture memory. Root-zone soil moisture memory is only but one feature characterizing the spectrum of soil moisture dynamics, which is analyzed here using a combination of long-term measurements and models. In particular, the spectrum of root-zone soil moisture content in a Mediterranean ecosystem is examined using 14-years of half-hourly measurements. A distinguishing hydro-climatic feature in such ecosystems is that sources (mainly rainfall) and sinks (mainly evapotranspiration) of soil moisture are roughly out of phase with each other. For over 4 decades of time scales and 7 decades of energy, the canonical shape of the measured soil moisture spectrum is shown to be approximately Lorentzian determined by the soil moisture variance and its memory but with two exceptions: the occurrences of a peak at diurnal-to-daily time scales and a weaker peak at near annual time scales. Model calculations and spectral analysis demonstrate that diurnal and seasonal variations in hydroclimate forcing responsible for variability in evapotranspiration had minor impact on the normalized shape of the soil moisture spectrum. However, their impact was captured by adjustments in the temporal variance. These findings indicate that precipitation and not evapotranspiration variability dominates the multi-scaling properties of soil moisture variability consistent with prior climate model simulations. Furthermore, the soil moisture memory inferred by the annual peak of soil moisture (340 d) is consistent with climate model simulations, while the memory evaluated from the loss function of a linearized mass balance approach leads to a smaller value (50 d), highlighting the effect of weak non-stationarity on soil moisture variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Experimental Study of the 3 He and Neutron Spin Structure at Low Q 2 and The Upgrade of the JLab Polarized 3 He Target

In the scattering process off a nuclear or nucleon target, the Gerasimov-Drell-Hearn (GDH) sum rule for real photons (Q 2 =0 where Q 2 &equiv; -q 2 with q the photon's 4-momentum) relates static properties of the target particle's ground state to dynamic properties of all its excited states. On the other side of the Q 2 spectrum, the Bjorken sum rule holds in the Bjorken limit Q 2 → ∞. Bjorken sum rule relates the final structure functions of the proton and neutron to the nucleon axial coupling constant in weak decay. These two sum rules belong to domains where calculations are achievable but use different degrees of freedom: hadronic degrees of freedom at low Q 2 versus partonic degrees of freedom at intermediate Q 2 . Meanwhile, different methods have been used to connect the two sum rules at finite Q 2 values: Chiral Perturbation Theory is used to expand the GDH sum rule while Operator Product Expansion is used to expand the Bjorken sum rule. In recent decades, improvements in polarized beam and polarized target techniques have made it possible to test theoretical predictions in the intermediate Q 2 region. During the Jefferson Lab (JLab) Hall A E97110 experiment, a precise measurement of polarized cross sections was performed at 0.02 < Q 2 < 0.3 GeV 2 using a polarized 3 He target as an effective polarized neutron target. The measured data allowed us to test predictions of Chiral Perturbation Theory at very low Q 2 . Furthermore, an extrapolation to the real photon point Q 2 =0 tests the GDH sum rule on the neutron. In order to reach the small angles necessary for the low Q 2 range, a new septum magnet was installed in Hall A for this experiment. Unfortunately, the magnet was mis-wired during initial running. There were therefore two periods for this experiment: the first period had the defective magnet due to mis-wiring; while in the second period, the magnet had been fixed and was working properly. The cross sections and the asymmetries for the first period must be extracted using a difficult and unusual method employing focal plane variables in the spectrometer. In the work described in this thesis, the target-plane and focal-plane method will be established first and confirmed by elastic scattering cross section measurements from carbon foil and 3 He targets. Cross sections and asymmetries of inelastic scattering were then extracted using the same method, which constitute the main physics results of this thesis. Preliminary results on ITT , the integral of the polarized cross section, were also extracted and are presented in this thesis. In addition to the data analysis of E97110, this thesis includes instrumentation work performed on the polarized 3 He target at JLab. For the 12 GeV program of JLab, the polarized 3 He target is being upgraded to satisfy new experimental requirements. The new target is a convection-based cell with two transfer tubes rather than a single transfer tube as in the former 6 GeV design. The first stage of the target upgrade aims to produce a 40 cm long, 10 amg target that can withstand 30 μA of electron beam current with an in-beam polarization of 55%, doubling the figure-of-merit from the 6 GeV time. Several studies were conducted at the Test Lab of JLab on the new convection cell. Work and improvement done on the Pulse Nuclear Magnetic Resonance (NMR) polarimetry and its test results will be presented. As this thesis is written, the target is being installed in Hall C of JLab in preparation for a 6-month long running from November 2019 to May 2020.

Ton, Nguyen↗

DOE OSTI2020

Detection of DoS Attacks Using ARFIMA Modeling of GOOSE Communication in IEC 61850 Substations

Integration of Information and Communication Technology (ICT) in modern smart grids (SGs) offers many advantages including the use of renewables and an effective way to protect, control and monitor the energy transmission and distribution. To reach an optimal operation of future energy systems, availability, integrity and confidentiality of data should be guaranteed. Research on the cyber-physical security of electrical substations based on IEC 61850 is still at an early stage. In the present work, we first model the network traffic data in electrical substations, then, we present a statistical Anomaly Detection (AD) method to detect Denial of Service (DoS) attacks against the Generic Object Oriented Substation Event (GOOSE) network communication. According to interpretations on the self-similarity and the Long-Range Dependency (LRD) of the data, an Auto-Regressive Fractionally Integrated Moving Average (ARFIMA) model was shown to describe well the GOOSE communication in the substation process network. Based on this ARFIMA-model and in view of cyber-physical security, an effective model-based AD method is developed and analyzed. Two variants of the statistical AD considering statistical hypothesis testing based on the Generalized Likelihood Ratio Test (GLRT) and the cumulative sum (CUSUM) are presented to detect flooding attacks that might affect the availability of the data. Our work presents a novel AD method, with two different variants, tailored to the specific features of the GOOSE traffic in IEC 61850 substations. The statistical AD is capable of detecting anomalies at unknown change times under the realistic assumption of unknown model parameters. The performance of both variants of the AD method is validated and assessed using data collected from a simulation case study. We perform several Monte-Carlo simulations under different noise variances. The detection delay is provided for each detector and it represents the number of discrete time samples after which an anomaly is detected. In fact, our statistical AD method with both variants (CUSUM and GLRT) has around half the false positive rate and a smaller detection delay when compared with two of the closest works found in the literature. Our AD approach based on the GLRT detector has the smallest false positive rate among all considered approaches. Whereas, our AD approach based on the CUSUM test has the lowest false negative rate thus the best detection rate. Depending on the requirements as well as the costs of false alarms or missed anomalies, both variants of our statistical detection method can be used and are further analyzed using composite detection metrics.

IEC 61850 electrical substations↗

DOE OSTI2023

Demonstration of a quantum-classical coprocessing protocol for simulating nuclear reactions

Quantum computers hold great promise for exact simulations of nuclear dynamical processes (e.g., scattering and reactions), which are paramount to the study of nuclear matter at the limit of stability and in the formation of chemical elements in stars. However, quantum simulations of the unitary (real) time dynamics of fermionic many-body systems require a currently prohibitive number of reliable and long-lived qubits. Here we propose a co-processing algorithm for the simulation of real-time dynamics in which the time evolution of the spatial coordinates is carried out on a classical processor, while the evolution of the spin degrees of freedom is carried out on quantum hardware. We demonstrate this hybrid scheme with the simulation of two neutrons scattering at the Lawrence Berkeley National Laboratory's Advanced Quantum Testbed. After implementing error mitigation strategies to improve the accuracy of the algorithm in addition to a combination of circuit compression techniques and tomography as methods to elucidate the onset of decoherence, our results validate the principle of the proposed co-processing scheme. A generalization of this present scheme will open the way for (real-time) path integral simulations of nuclear scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Microwave-Based Quantum Control and Coherence Protection of Tin-Vacancy Spin Qubits in a Strain-Tuned Diamond-Membrane Heterostructure

Robust spin-photon interfaces in solids are essential components in quantum networking and sensing technologies. Ideally, these interfaces combine a long-lived spin memory, coherent optical transitions, fast and high-fidelity spin manipulation, and straightforward device integration and scaling. The tin- vacancy center (SnV) in diamond is a promising spin-photon interface with desirable optical and spin properties at 1.7 K. However, the SnV spin lacks efficient microwave control, and its spin coherence degrades with higher temperature. In this work, we introduce a new platform that overcomes these challenges—SnV centers in uniformly strained thin diamond membranes. The controlled generation of crystal strain introduces orbital mixing that allows microwave control of the spin state with 99.36(9)% gate fidelity and spin coherence protection beyond a millisecond. Moreover, the presence of crystal strain suppresses temperature-dependent dephasing processes, leading to a considerable improvement of the coherence time up to 223(10) μs at 4 K, a widely accessible temperature in common cryogenic systems. Critically, the coherence of optical transitions is unaffected by the elevated temperature, exhibiting nearly lifetime-limited optical linewidths. Combined with the compatibility of diamond membranes with device integration, the demonstrated platform is an ideal spin-photon interface for future quantum technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics↗

DOE OSTI2025

Electrodermal activity as a proxy for sweat rate monitoring during physical and mental activities

Electrodermal activity has long been used for mental activity monitoring by measuring skin conductance at specific locations, such as fingertips, with high sweat gland density. However, electrodermal activity has not been considered useful for physical activity monitoring, where large sweat volumes are generated, resulting in the accumulation of sweat at the skin–electrode interface and, thus, preventing further dynamic response to sweating events. Here we show that electrodermal activity can be used as a proxy for sweat loss measurement under both low and high physical activity levels. We use wearable sweat sensors that consist of water-permeable electrodes and microfluidic-based sweat analysers, and show that skin conductance is proportional to the instantaneous sweat loss. We demonstrate that sweat loss during exercise can be estimated by integrating skin conductance over time, which can be applied to assess the body hydration status of exercisers. From multisite measurements of skin conductance, we show that the wrist, forearm and upper arm are reflective of physical activity levels, whereas the finger is indicative of mental activity. Simultaneous measurement of two different sites selectively decouples mental and physical activities.

Biomedical engineering↗

DOE OSTI2024

Sequential Inference of Hospitalization Electronic Health Records Using Probabilistic Models

In the dynamic hospital setting, decision support can be a valuable tool for improving patient outcomes. Data-driven inference of future outcomes is challenging in this dynamic setting, where long sequences such as laboratory tests and medications are updated frequently. This is due in part to heterogeneity of data types and mixed-sequence types contained in variable length sequences. In this work we design a probabilistic unsupervised model for multiple arbitrary-length sequences contained in hospitalization Electronic Health Record (EHR) data. The model uses a latent variable structure and captures complex relationships between medications, diagnoses, laboratory tests, neurological assessments, and medications. It can be trained on original data, without requiring any lossy transformations or time binning. Inference algorithms are derived that use partial data to infer properties of the complete sequences, including their length and presence of specific values. We train this model on data from subjects receiving medical care in the Kaiser Permanente Northern California integrated healthcare delivery system. The results are evaluated against held-out data for predicting the length of sequences and presence of Intensive Care Unit (ICU) in hospitalization bed sequences. Our method outperforms a baseline approach, showing that in these experiments the trained model captures information in the sequences that is informative of their future values.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Quadrature Based Neural Network Learning of Stochastic Hamiltonian Systems

Hamiltonian Neural Networks (HNNs) provide structure-preserving learning of Hamiltonian systems. In this paper, we extend HNNs to structure-preserving inversion of stochastic Hamiltonian systems (SHSs) from observational data. We propose the quadrature-based models according to the integral form of the SHSs’ solutions, where we denoise the loss-by-moment calculations of the solutions. The integral pattern of the models transforms the source of the essential learning error from the discrepancy between the modified Hamiltonian and the true Hamiltonian in the classical HNN models into that between the integrals and their quadrature approximations. This transforms the challenging task of deriving the relation between the modified and the true Hamiltonians from the (stochastic) Hamilton–Jacobi PDEs, into the one that only requires invoking results from the numerical quadrature theory. Meanwhile, denoising via moments calculations gives a simpler data fitting method than, e.g., via probability density fitting, which may imply better generalization ability in certain circumstances. Numerical experiments validate the proposed learning strategy on several concrete Hamiltonian systems. The experimental results show that both the learned Hamiltonian function and the predicted solution of our quadrature-based model are more accurate than that of the corrected symplectic HNN method on a harmonic oscillator, and the three-point Gaussian quadrature-based model produces higher accuracy in long-time prediction than the Kramers–Moyal method and the numerics-informed likelihood method on the stochastic Kubo oscillator as well as other two stochastic systems with non-polynomial Hamiltonian functions. Moreover, the Hamiltonian learning error εH arising from the Gaussian quadrature-based model is lower than that from Simpson’s quadrature-based model. These demonstrate the superiority of our approach in learning accuracy and long-time prediction ability compared to certain existing methods and exhibit its potential to improve learning accuracy via applying precise quadrature formulae.

Mathematics↗

DOE OSTI2022

Search for High-energy Neutrino Emission from Galactic X-Ray Binaries with IceCube

We present the first comprehensive search for high-energy neutrino emission from high- and low-mass X-ray binaries conducted by IceCube. Galactic X-ray binaries are long-standing candidates for the source of Galactic hadronic cosmic rays and neutrinos. The compact object in these systems can be the site of cosmic-ray acceleration, and neutrinos can be produced by interactions of cosmic rays with radiation or gas, in the jet of a microquasar, in the stellar wind, or in the atmosphere of the companion star. We study X-ray binaries using 7.5 yr of IceCube data with three separate analyses. In the first, we search for periodic neutrino emission from 55 binaries in the Northern Sky with known orbital periods. In the second, the X-ray light curves of 102 binaries across the entire sky are used as templates to search for time-dependent neutrino emission. Finally, we search for time-integrated emission of neutrinos for a list of 4 notable binaries identified as microquasars. In the absence of a significant excess, we place upper limits on the neutrino flux for each hypothesis and compare our results with theoretical predictions for several binaries. In addition, we evaluate the sensitivity of the next generation neutrino telescope at the South Pole, IceCube-Gen2, and demonstrate its power to identify potential neutrino emission from these binary sources in the Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗