Search NASA⌕ Search

SEARCH · Search NASA

Results for “direct counts”

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 55 records · Page 3

Quantum dynamics simulation of the advection-diffusion equation

The advection-diffusion equation is simulated via several quantum algorithms. Three formulations are considered: (1) Trotterization, (2) variational quantum time evolution (VarQTE), and (3) adaptive variational quantum dynamics simulation (AVQDS). These schemes were originally developed for the Hamiltonian simulation of many-body quantum systems. The finite-difference discretized operator of the transport equation is formulated as a Hamiltonian and solved without the need for ancillary qubits. Computations are conducted on a quantum simulator (IBM Qiskit Aer) and a superconducting quantum hardware (IBM Fez). The former emulates the latter without the noise. The actual hardware implementation experiences significant noise. The results of the quantum simulator are compared with data from direct numerical simulation (DNS) with infidelities of the order 10 −5 . In the quantum simulator, Trotterization is observed to have the lowest infidelity and is suitable for fault-tolerant computation. The AVQDS algorithm requires the lowest gate count and circuit depth. The VarQTE algorithm is the next best in terms of gate counts, but the number of its optimization variables is directly proportional to the number of qubits. Due to current hardware limitations, Trotterization cannot be implemented, as it has an overwhelmingly large number of operations. Meanwhile, AVQDS and VarQTE can be executed at the hardware level. These algorithms present a new paradigm for computational transport phenomena on quantum computers.

Alipanah, Hirad [Univ. of Pittsburgh, PA (United S↗

Characteristics of Veterans With Depression Who Use the Veterans Choice Program of the Veterans Health Administration

Objective: The Veterans Choice Program (VCP) of the Veterans Health Administration (VHA) allowed eligible veterans to use their benefits with participating providers outside the VHA. The authors aimed to identify characteristics of veterans with depression who used or did not use mental health care through the VCP. Methods: In this cross-sectional study, the authors analyzed secondary data from the national VHA Corporate Data Warehouse. VHA administrative data were linked with VCP claims to examine characteristics of VCP-eligible veterans with depression. Further, the study sample included 595,943 unique veterans who were enrolled in the VHA before 2013, were eligible for the VCP in 2016, were alive in 2018, and had an assessed Patient Health Questionnaire–9 (PHQ-9) score or depressive disorder diagnosis documented in the VHA between 2016 and 2018. Results: Veterans who used the VCP had lower medical comorbidity scores and lived in less socioeconomically disadvantaged counties, compared with veterans who received only VHA care. VCP veterans were also more likely to have a PHQ-9 score assessment and to have higher mean depression scores. Mean counts of annual mental health visits per 1,000 veterans were markedly higher for direct VHA care than for care provided via the VCP. As a percentage of the total counts of visits per 1,000 veterans across the VCP and VHA, residential programs and outpatient procedures were the services that were most frequently delivered through the VCP. Conclusions: Between 2016 and 2018, the VCP was used primarily to augment mental health care provided by the VHA, rather than to fill a gap in care.

60 APPLIED LIFE SCIENCES↗

Revealing the landscape of globally color-dual multi-loop integrands

Abstract We report on progress in understanding how to construct color-dual multi-loop amplitudes. First we identify a cubic theory, semi-abelian Yang-Mills, that unifies many of the color-dual theories studied in the literature, and provides a prescriptive approach for constructingD-dimensional color-dual numerators through one-loop directly from Feynman rules. By a simple weight counting argument, this approach does not further generalize to two-loops. As a first step in understanding the two-loop challenge, we use aD-dimensional color-dual bootstrap to successfully construct globally color-duallocaltwo-loop four-point nonlinear sigma model (NLSM) numerators. The double-copy of these NLSM numerators with themselves, pure Yang-Mills, and$$\mathcal{N}$$= 4 super-Yang-Mills correctly reproduce the known unitarity constructed integrands of special Galileons, Born-Infeld theory, and Dirac-Born-Infeld-Volkov-Akulov theory, respectively. Applying our bootstrap to two-loop four-point pure Yang-Mills, we exhaustively search the space of local numerators and find that it fails to satisfy global color-kinematics duality, completing a search previously initiated in the literature. We pinpoint the failure to the bowtie unitarity cut, and discuss a path forward towards non-local construction of color-dual integrands at generic loop order.

Physics↗

Rapid neutron and gamma-ray source localization using machine learning

Rapid localization of radiation sources is critical for applications including nuclear emergency response, safeguards, and security. However, conventional imaging systems such as neutron scatter cameras and Compton cameras depend on rare coincidence events, which often result in long acquisition times. In this work, we address the challenge of rapid source localization by developing a machine learning approach to predict the direction of a single radiation source using only count rates from an array of neutron and gamma-ray detectors. The proposed model is a fully connected neural network (FCNN) trained using Monte Carlo simulation data from a 252 Cf source. The model hyperparameters are optimized with a small set of routine 252 Cf measurements. We benchmarked the performance of the trained and optimized machine learning model using additional 252 Cf , 137 Cs , and PuBe measurements under laboratory conditions with varying source-detector configurations. For these measurements, the machine learning model achieved a mean localization error smaller than 30° with 3 x 10 3 system counts, corresponding to 8 s measurement time for the imaging system used in this work. In this low-statistics regime, the method outperformed traditional scatter-based imaging by more than 75% in localization accuracy for the evaluated measurement configurations. These results demonstrate that a machine learning-based approach can significantly reduce the time required for accurate single-source localization, providing a robust and computationally efficient alternative to traditional imaging systems in time-critical nuclear security and emergency response scenarios.

Gamma-ray imaging↗

Comprehensive Physical Activity Assessment During U.S. Army Basic Combat Training

Abstract Alemany, JA, Pierce, JR, Bornstein, DB, Grier, TL, Jones, BH, and Glover, SH. Comprehensive physical activity assessment during U.S. Army Basic Combat Training. J Strength Cond Res 36(12): 3505–3512, 2022—Physical activity (PA) volume, intensity, and qualitative contextual information regarding activity type and loads carried are limited during U.S. Army Basic Combat Training (BCT). The purpose of this study was to characterize daily (05:00–20:00 hours) PA during BCT using a comprehensive approach. During 2 10-week BCT cycles ( n = 40 trainees per cycle), pedometers, accelerometers, and direct observation were used to estimate daily step count, PA volume, and intensity. Physical activity intensity was categorized by metabolic equivalents (METs) such as “sedentary” (1–2 METs), “light” (2–3 METs), “moderate” (3–6 METs), or “vigorous” (≥6 METs). Daily PA data were analyzed longitudinally using linear mixed models, with significance set at p ≤ 0.05. The mean daily step count was 13,459 ± 4,376 steps, and the mean daily accelerometer-assessed PA volume and intensity were as follows: sedentary: 505 ± 98 minutes, light: 190 ± 78 minutes, moderate: 168 ± 51 minutes, and vigorous: 14 ± 14 minutes, with no differences between cycles for all measures ( p > 0.50). Cumulative time on feet (∼50%) and sitting (20–25%) accounted for most daily activity types during both cycles. Trainees, on average, carried between 3 and 9 kg, and ≥9 kg, for 60% and 10% of the monitored day, respectively. Basic Combat Training's physical demands are high, where trainees achieved 1.7 to 2.7 times greater daily ambulation and 6 times the recommended weekly moderate-to-vigorous PA compared with civilian counterparts and performed weight-bearing load carriage for nearly half of the day. Basic Combat Training-associated PA may increase injury risk among trainees unaccustomed to arduous PA and exercise. Implementing national PA policies to improve physical fitness and facilitate acclimatization to BCT's high physical demands could reduce public health burdens and military nonreadiness.

Sport Sciences↗

Data Analysis of X-ray Panels Used in X-ray CT Scientific Instruments

This report is a follow-up study on M. Skeate’s study on System Drift Detection for Health Monitoring (Skeate, 2023). The following report highlights analyses of the panel’s behavior with extended use. This includes effects of burn-in on the panel from over-use and bad pixels (defined in Methodology.) Using data-forward statistical analyses across a sequence of scans, and given that the other conditions present in the data can be replicated, this study show that overtime use of panels does not affect bad pixel count. Also, I present conclusive evidence of a direct relationship between the regions of the panel that are exposed to radiation and burn-in damage to that region overtime. Additionally, this study also highlights the effect of the duty cycle on the dark current change in the panel.

42 ENGINEERING↗

Towards a Unified Low-Cost Flow Plate, Flow-Field, PTL Solution for Proton Exchange Membrane Electrolyzers

Proton exchange membrane (PEM) water electrolysis is a highly efficient method for hydrogen production. Research cells typically consist of one proton exchange membrane, two catalyst layers, two porous transport layers, two flow-field plates, and two endplates. In commercial systems, the machined flow-field plates that are employed in research cells are typically replaced by stamped parts or open mesh material solutions to reduce manufacturing cost at scale. Nonetheless, the cell contains about 8 total interfaces: bipolar plate / flow plate material / porous transport medium / electrode / membrane / electrode / porous transport medium / flow plate material / bipolar plate. All these materials and interfaces need to be optimized for maximum performance and efficiency. Reducing the amount of interfaces by combining individual cell components directly benefits the fabrication cost (by reducing the parts count and the needs for surface coatings) and the electrochemical performance (by reducing ohmic losses). We have designed a novel PEM electrolysis cell with a piece of channeled titanium felt functioning as both the anode flow-field and the PTL, referred to as the channeled diffusion layer (CDL). The pores of the felt facilitate both in-plane and through-plane diffusion, ensuring maximum catalyst utilization while also minimizing mass transport loss. The titanium felt can be mass manufactured with existing stamping and forming methods and is therefore a promising candidate to reduce the capital cost of PEM electrolyzers whilst improving hydrogen production efficiency. Experiments conducted with 3mg IrOx/cm2 loading MEAs have shown a approximately 40% boost in peak current by implementing the CDL design. Low catalyst-loading MEAs are being tested in ongoing experiments and their results will be discussed and compared.

08 HYDROGEN↗

Hierarchical Control and Stability Analysis for a Nonisolated Grid-Tied DC Energy Router Integrating Energy Storage and Partial Distributed Generation

This article proposes a nonisolated dc converter-based energy router (dc-ER) and its operating strategy. Here, the intent is to integrate energy storage (ES), distributed generation (DG), local load, and dc power grid in an autonomous and more efficient way. The ES/DG power ports are coupled with each other in partial-type connections, which obtains higher voltage supply gain and direct input–output power transmission. High-voltage supply gain allows a wide solar power tracking range and the possibility of optimal battery charging/discharging. Direct input–output power transmission improves the energy conversion efficiency. In this article, the operating modes of dc-ER are first analyzed, followed by the optimal hardware design counting the DG current ripple minimization and the limitation caused by the power flow direction. Second, the mathematical model of dc-ER is deduced. The hierarchical control is then proposed to manage port energy in a flexible manner. The stability is analyzed by using impedance modeling. Finally, experimental and simulation results verify the superiorities of the proposed topology in terms of flexible operating mode transition and high-power conversion efficiency.

25 ENERGY STORAGE↗

A Novel Measurement of the Anomalous Muon Spin Precession Frequency in the Muon $g-2$ Experiment at Fermilab

The Muon $g-2$ Experiment operated at Fermi National Accelerator Laboratory (FNAL, or Fermilab) between 2018 and 2023 to produce the world's most precise measurement of the muon's \textit{anomalous magnetic moment}, $a_\mu = \frac{g_\mu - 2}{2}$, which expresses the relative deviation in the muon's $g$-factor from a baseline theoretical expectation that $g_\mu = 2$. In the Standard Model of particle physics, $g_\mu > 2$ and hence $a_\mu > 0$ by a calculable amount that depends on all possible interactions between the muon and all other fundamental particles, including any potentially undiscovered interactions beyond the Standard Model. For this reason, measurements of the electron anomaly $a_e$ and later the muon anomaly $a_\mu$ have helped guide the development of the Standard Model since the inception of quantum field theory, and the measured value of $a_\mu$ provides a valuable constraint for new hypotheses that extend the Standard Model. As of 2006, the leading measurement and Standard Model prediction for $a_\mu$ exhibited tension at the level of about three standard deviations, motivating an improved measurement at Fermilab that could test the tension more precisely. The experiment functions by storing a polarized beam of $\mu^+$ in a uniform magnetic field, which simultaneously induces circular motion and spin precession. As the stored muons undergo the Michel decay $\mu^+ \to e^+ + \nu_e + \bar{\nu}_\mu$, mediated by the parity-violating weak interaction, the rest-frame $e^+$ emission direction is correlated with the parent $\mu^+$ spin orientation. Boosting into the laboratory frame encodes this correlation in the decay $e^+$ energy, which is higher when the emission (i.e. $\mu^+$ spin direction) is aligned with the $\mu^+$ momentum, and lower when opposite. Detectors then count the rate of high-energy decay $e^+$, which modulates at the difference between the $\mu^+$ revolution and spin precession frequencies. This observed frequency, called the \textit{anomalous spin precession frequency} $\omega_a$, is directly proportional to $a_\mu$. The extraction of $\omega_a$ proceeds by fitting the time spectrum of detected $e^+$, which requires precise modeling of the $\omega_a$ oscillation as well as any perturbations from beam dynamics and detector acceptance. Using the $\omega_a$ analysis presented in this work, based on Runs 4 -- 6 of the Muon $g-2$ Experiment at Fermilab, we find that $a_\mu = 0.001\,165\,920\,738(162)$ with a relative uncertainty of 139 parts per billion.

Barrett, Tyler [Cornell U.]↗

Search for New Physics with the Compact Muon Solenoid Experiment and QIS-enabled Technology

Understanding the fundamental nature of dark matter (DM)---its cosmological origin, constituents, and interactions---is one of the most important questions in fundamental science today. In this thesis, I present two novel and highly complementary approaches to cover the gaps in sensitivity of current DM searches. The searches are enabled by a first-of-its-kind reconstruction technique to search for hidden-sector particles using the Compact Muon Solenoid (CMS) and by new advances in quantum sensing technology to search for axions and hidden-sector DM. In the first part of this thesis, I present a search for long-lived hidden sector particles, predicted by many extensions of the SM, using a novel technique to reconstruct decays of long-lived particles (LLPs) in the CMS muon detector. The innovative LLP reconstruction technique is sensitive to a broad range of LLP decays and to LLP masses below GeV. The search yields competitive sensitivity for proper lifetime 0.1--1000 m with the full Run 2 dataset recorded at the LHC between 2016--2018 at $\sqrt{s} = 13~$TeV. To extend the physics reach of this novel muon detector shower (MDS) signature, I present the model-independence of MDS and the reinterpretation of the search to a large number of LLP models, demonstrating its complementarity with proposed and existing dedicated LLP experiments. Finally, I present a new dedicated MDS trigger that improves the trigger efficiency by at least an order of magnitude and was deployed in 2022, at the start of Run 3 of the LHC operations. In the second part of the thesis, I present for the first time, the use of a novel quantum sensor, the low-noise and single-photon sensitive superconducting nanowire single photon detectors (SNSPDs), to directly detect dark matter. The low detection threshold and ultra-low dark count rate of SNSPDs can close the gap in DM discovery reach due to the current limitations in detector sensitivity. I will present my work on the development and characterization of SNSPDs for two entirely new experiments to directly detect axions via absorption and hidden-sector DM via electron scattering. The search for axions employs a novel broadband reflector technique with the Broadband Reflector Experiment for Axion Detection (BREAD). A unique parabolic mirror is then used to focus axion-converted photons to the SNSPDs, extending the reach to axion masses of 0.04--1 eV. On the other hand, by coupling the SNSPDs with gallium arsenide, a bright cryogenic scintillator well matched to SNSPD detection, a prototype sensing system can be built as a basis of new direct DM detection experiments capable of extending the discovery to DM masses as low as 1 MeV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A zone-level occupancy counting system for commercial office spaces using low-resolution time-of-flight sensors

Understanding the locations of occupants in a commercial built environment is critical for realizing energy savings by delivering lighting, heating, and cooling only where it is needed. In this paper, we present an indoor occupancy counting system using a sparse array of inexpensive, low-resolution, and privacy-preserving time-of-flight sensors. We develop and validate an algorithm for zonal occupancy counting that can deal with multiple people walking underneath the sensors in arbitrary directions, and evaluate the system both in realistic simulations of office spaces and in a real-world installation. Finally, we found that our system has an error rate of around 0.4%, resulting in highly accurate person localization and zone counting using only a few sensors per space.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Stopping criteria for ending autonomous, single detector radiological source searches

While the localization of radiological sources has traditionally been handled with statistical algorithms, such a task can be augmented with advanced machine learning methodologies. The combination of deep and reinforcement learning has provided learning-based navigation to autonomous, single-detector, mobile systems. However, these approaches lacked the capacity to terminate a surveying/search task without outside influence of an operator or perfect knowledge of source location (defeating the purpose of such a system). Two stopping criteria are investigated in this work for a machine learning navigated system: one based upon Bayesian and maximum likelihood estimation (MLE) strategies commonly used in source localization, and a second providing the navigational machine learning network with a “stop search” action. A convolutional neural network was trained via reinforcement learning in a 10 m × 10 m simulated environment to navigate a randomly placed detector-agent to a randomly placed source of varied strength (stopping with perfect knowledge during training). The network agent could move in one of four directions (up, down, left, right) after taking a 1 s count measurement at the current location. During testing, the stopping criteria for this navigational algorithm was based upon a Bayesian likelihood estimation technique of source presence, updating this likelihood after each step, and terminating once the confidence of the source being in a single location exceeded 0.9. A second network was trained and tested with similar architecture as the previous but which contained a fifth action: for self-stopping. The accuracy and speed of localization with set detector and source initializations were compared over 50 trials of MLE-Bayesian approach and 1000 trials of the CNN with self-stopping. The statistical stopping condition yielded a median localization error of ~1.41 m and median localization speed of 12 steps. The machine learning stopping condition yielded a median localization error of 0 m and median localization speed of 17 steps. This work demonstrated two stopping criteria available to a machine learning guided, source localization system.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Scalable Thin Light-Emitting Diode (LED) Light Sheet Platform

The goal of the Scalable Thin Light-Emitting Diode (LED) Light Sheet Platform project is to derisk the manufacturing scalability for a high chip count heterogeneous smart lighting platform, with the potential to enable significant energy savings through dynamic directional light control and custom. Specifically, the project derisks the scalability of a disruptive new computer controlled microassembly fabrication process SRI International is developing, to address fundamental cost barriers to mass production of high chip count systems. This is key for enabling mass adoption and thus maximum societal energy savings impact. The LED light sheet technology is a smart illumination platform that can be applied to lighting, signage, and display. It has features and a form factor similar to bendable or conformal OLED light sheets but uses more efficient LEDs with a remote phosphor layer that is very close to the LED to deliver a luminous efficacy that exceeds 125 lumen/Watt and enable > 50% Lighting Application Efficiency (LAE) energy savings. During Budget Period 1 (BP1), all of the BP1 milestones were successfully completed: (M2.0.1) Design demonstrator details that will meet the final project goals, (M6.1.1) Show lighting output model can predict illuminance, spatial, spectral distribution for specific light-sheet designs and use case, (M3.2.1) Automated loading supports 200 chip arrays, and (M3.3.1) Finish first interconnect process run for 200 chip array. The BP1 Go/No-Go Decision Point (G/NG 1) 200 chip demonstrator sample was in the process of final assembly and test at the end of BP1 on December 31, 2023. The objective was successfully achieved on January 16, 2024, when a wired 200-chip sample underwent confirmation testing that demonstrated a 93.75% first pass electrical yield that successfully exceeded the 75% first-pass electrical yield BP1 G/NG requirement. In addition, all 93.75% of the sample devices lit up, further demonstrating the ability to assemble and transfer small LED chips without damaging their functionality.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Time-Resolved Line Shapes of Single Quantum Emitters via Machine Learned Photon Correlations

Solid-state single-photon emitters (SPEs) are quantum light sources that combine atomlike optical properties with solid-state integration and fabrication capabilities. SPEs are hindered by spectral diffusion, where the emitter’s surrounding environment induces random energy fluctuations. Timescales of spectral diffusion span nanoseconds to minutes and require probing single emitters to remove ensemble averaging. Photon correlation Fourier spectroscopy (PCFS) can be used to measure time-resolved single emitter line shapes, but is hindered by poor signal-to-noise ratio in the measured correlation functions at early times due to low photon counts. Here, we develop a framework to simulate PCFS correlation functions directly from diffusing spectra that match well with experimental data for single colloidal quantum dots. We use these simulated datasets to train a deep ensemble autoencoder machine learning model that outputs accurate, noiseless, and probabilistic reconstructions of the noisy correlations. Using this model, we obtain reconstructed time-resolved single dot emission line shapes at timescales as low as 10 ns, which are otherwise completely obscured by noise. This enables PCFS to extract optical coherence times on the same timescales as Hong-Ou-Mandel two-photon interference, but with the advantage of providing spectral information in addition to estimates of photon indistinguishability. Further, our machine learning approach is broadly applicable to different photon correlation spectroscopy techniques and SPE systems, offering an enhanced tool for probing single emitter line shapes on previously inaccessible timescales.

74 ATOMIC AND MOLECULAR PHYSICS↗

Automated Bacterial Identification and Morphological Feature Analysis in Low‐Dose Cryo‐EM Using YOLOv11

Bacteria rapidly adapt to environmental cues through morphological and ultrastructural changes that correlate with physiology and behavior. Cryogenic transmission electron microscopy (cryo‐TEM) can capture these phenotypic changes in near‐native, vitrified states, but manual analysis of low‐dose micrographs is labor intensive and limits throughput. Here, we present an end‐to‐end workflow that combines low‐dose cryo‐TEM imaging with a YOLOv11‐based instance‐segmentation model to automatically identify bacteria and quantify key structural features directly from the micrographs. This workflow enables (i) robust bacterial localization and counting from low‐magnification atlas/montage images, (ii) automated measurements of cell‐envelope (outer–inner membrane) thickness and anisotropy from higher‐magnification views, and (iii) detection and quantification of bacteria–flagella interactions, including overlap length and curvature metrics for interacting versus noninteracting flagella. Using Pantoea sp. YR343 grown under distinct media conditions, we show that the automated measurements agree with manual annotations while substantially reducing analysis time. Together, these tools provide a practical framework for scalable bacterial identification and quantitative phenotyping in low‐dose cryo‐TEM datasets and establish a foundation for extending cryo‐TEM image analysis toward higher‐throughput studies of microbial heterogeneity and biointerfaces.

YOLOv11↗

Integration of the Kromek D3S Detector and Spot Robot For Secondary Inspections

Inspecting vehicles and containers for the presence of nuclear material is a challenging task for border control and security. When performed manually by inspectors, this task also has an associated risk of exposing the inspectors to unknown radiation. With the advent of agile, easy-to-program, quadruped robots like the Boston Dynamics Spot, automation of secondary inspection can improve the efficiency of the inspection process and alleviates the radiation risks to inspectors. In this project, Brookhaven National Laboratory and the University of Massachussetts at Lowell explored how to automate a simple secondary inspection mission. The Spot robot comes with its own software development kit (SDK) that allows clients/users to write custom code in the Python programming language to control the robot. Spot also has a payload computer called Spot-CORE, which runs the Ubuntu Linux operating system and allows users to integrate external sensors, such as a radiation detector, with Spot. In this study, the Kromek D3S detector has been integrated with Spot via the Spot-CORE, allowing Spot to capture gamma spectra and neutron counts for a specified acquisition period. Two custom routines, search and confirmation, have been developed and executed in this specified order. The search routine directs Spot to go around the nearest obstacle, e.g., vehicle and container, in a preset distance and step to collect gamma and neutron gross counts with the D3S detector. The radiation data and the robot location corresponding to each step are stored and fed to the confirmation routine at the end of the search. The confirmation routine then navigates Spot to the locations of the highest gamma or neutron counts to perform a long, e.g., one minute, measurement, and gives the operators the signature gamma spectra and neutron counts at the hotspots. This paper presents a detailed description of this automated system along with results of the preliminary tests in identifying the location and signature of a 137Cs radiation source in a vehicle.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Scaled ILU Smoothers for Navier-Stokes Pressure Projection

Incomplete LU (ILU) smoothers are effective in the algebraic multigrid (AMG) V-cycle for reducing high-frequency components of the error. However, the requisite direct triangular solves are comparatively slow on GPUs. Previous work has demonstrated the advantages of Jacobi iteration as an alternative to direct solution of these systems. Depending on the threshold and fill-level parameters chosen, the factors can be highly nonnormal and Jacobi is unlikely to converge in a low number of iterations. We demonstrate that row scaling can reduce the departure from normality, allowing us to replace the inherently sequential solve with a rapidly converging Richardson iteration. There are several advantages beyond the lower compute time. Scaling is performed locally for a diagonal block of the global matrix because it is applied directly to the factor. Further, an ILUT Schur complement smoother maintains a constant GMRES iteration count as the number of MPI ranks increases, and thus parallel strong-scaling is improved. Our algorithms have been incorporated into hypre, and we demonstrate improved time to solution for linear systems arising in the Nalu-Wind and PeleLM pressure solvers. For large problem sizes, GMRES+AMG executes at least five times faster when using iterative triangular solves compared with direct solves on massively parallel GPUs.

algebraic multigrid↗

Congruity of genomic and epidemiological data in modelling of local cholera outbreaks

Cholera continues to be a global health threat. Understanding how cholera spreads between locations is fundamental to the rational, evidence-based design of intervention and control efforts. Traditionally, cholera transmission models have used cholera case-count data. More recently, whole-genome sequence data have qualitatively described cholera transmission. Integrating these data streams may provide much more accurate models of cholera spread; however, no systematic analyses have been performed so far to compare traditional case-count models to the phylodynamic models from genomic data for cholera transmission. Here, we use high-fidelity case-count and whole-genome sequencing data from the 1991 to 1998 cholera epidemic in Argentina to directly compare the epidemiological model parameters estimated from these two data sources. We find that phylodynamic methods applied to cholera genomics data provide comparable estimates that are in line with established methods. Our methodology represents a critical step in building a framework for integrating case-count and genomic data sources for cholera epidemiology and other bacterial pathogens.

59 BASIC BIOLOGICAL SCIENCES↗