Search NASA⌕ Search

SEARCH · Search NASA

Results for “Factor Analysis, Statistical”

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 127 records · Page 7

General-Purpose Unsupervised Cyber Anomaly Detection via Non-Negative Tensor Factorization

Distinguishing malicious anomalous activities from unusual but benign activities is a fundamental challenge for cyber defenders. Prior studies have shown that statistical user behavior analysis yields accurate detections by learning behavior profiles from observed user activity. These unsupervised models are able to generalize to unseen types of attacks by detecting deviations from normal behavior, without knowledge of specific attack signatures. However, approaches proposed to date based on probabilistic matrix factorization are limited by the information conveyed in a two-dimensional space. Non-negative tensor factorization, on the other hand, is a powerful unsupervised machine learning method that naturally models multi-dimensional data, capturing complex and multi-faceted details of behavior profiles. Herein, our new unsupervised statistical anomaly detection methodology matches or surpasses state-of-the-art supervised learning baselines across several challenging and diverse cyber application areas, including detection of compromised user credentials, botnets, spam e-mails, and fraudulent credit card transactions.

97 MATHEMATICS AND COMPUTING↗

Residential Building Energy Efficiency Field Studies: Low-Rise Multifamily

In recent years, the U.S. Department of Energy (DOE) has conducted a series of research studies to validate energy efficient building technologies in the field. Much of the work has focused on single-family construction, and some has also addressed commercial energy codes. The work detailed in this DOE-funded study (EE0007616) focuses on low-rise multifamily buildings (three stories or fewer above grade) in various regions of the United States, and reports on how state-level building codes are being implemented, both in terms of observed characteristics and also in terms of estimated energy impacts. Nearly 100 buildings across four states—Illinois, Minnesota, Oregon, and Washington—were sampled, which represent a range of climate types from mild temperature to very cold continental. Both common entry and outdoor entry buildings were included, and a parallel research project evaluated envelope air tightness and current still-evolving air tightness testing methods. Finally, a set of structured interviews of building designers and other relevant professionals was carried to out to gain more insight into this market. To the greatest extent possible, the methodology developed under the project for low-rise multifamily buildings mirrored the approach established by Pacific Northwest National Laboratory (PNNL) for single-family residential buildings (https://www.energy.gov/eere/buildings/downloads/residential-building-energy-code-field-study). This included the general approach to sampling, recruitment, and data collection, as well as data analysis and presentation. The range of permitting dates for the sites encompassed two energy code cycles in most regions. All states in the study had adopted a variation of the International Energy Conservation Code (IECC) for the structure of their state code. The low-rise multifamily occupancy presents a hybrid building type: most of the building’s conditioned floor area was covered by the residential chapter of the code while portions of the building (such as corridors and common spaces) fell under the commercial code chapter. The key items assessed in this work were: Building Shell—exterior wall insulation, ceiling insulation, foundation insulation, windows. Common Areas—HVAC and lighting. Living Units—lighting, ventilation. A few items were not assessed in detail, given their relative paucity in this occupancy type; these included duct leakage, pipe insulation, and hot water circulation controls. Building characteristics were collected via a combination of architectural, mechanical, electrical, and plumbing plan reviews and field inspections, and entered into a spreadsheet-based tool that was later queried to build a database. Data went through quality control both upon arrival and via a later semi-automated review and assurance process. Most of the data are presented graphically so that the reader can quickly assess compliance with the applicable energy codes (both by state and by code year). As a final step, EnergyPlus™ simulations were created for all buildings in the study to estimate both the as-found energy use intensity (EUI) and the energy and CO 2 that could be saved if features that were found to not meet code minimums were brought up to code. The savings estimates were tabulated for each of the four states in the study. The research team found that the single-family approach was largely applicable to low-rise multifamily buildings. This applies to both the data collection and the prototype EUI analysis. Most of the occupied space is living units and falls under residential energy codes, and many characteristics use similar envelope construction and relatively straightforward mechanical systems and lighting. One of the most challenging aspects of this work was to build an effective spreadsheet-based data collection instrument that could allow efficient collection of both building plan and field data. The research team is of the view that other methods could be equally effective if the work is done carefully with diligent quality control. The primary findings for the work center around the thermal envelope and mechanical systems and lighting at the sites: For thermal envelope components, the majority of buildings met or were better than the prescriptive code.This suggests that building designers and builders are aware of code requirements. In some cases, surveyed buildings were designed to qualify for energy efficiency certification programs. These buildings made up at least 20% of sampled buildings in each state. Almost all buildings met mechanical system efficiency requirements (for both living units and common areas). In some cases, sites employed systems that were considerably more efficient than required by the applicable energy code. Dwelling units had a majority of high-efficacy lighting, often in excess of the state’s residential code requirements. While high-efficacy fixtures were also typical in common areas (corridors and stairwells), lighting power densities (LPDs) in these areas were sometimes higher than levels dictated by the applicable part of the state commercial energy code. The simulation models run on a series of low-rise multifamily prototypes, informed by a composite of the field data collected, calculated annual EUIs of between 20 and 50 kBtu/ft2-yr, with the range representing the effects of both building characteristics and building location (climate zone). A detailed process (based on simulations of prototype buildings) was used to estimate the amount of avoided energy use that would occur if 100% adherence to energy codes were attained. The results indicated modest savings are attainable for items such as window thermal performance and common area lighting. The result is overall only a modest potential for additional energy savings, averaging about 10% of EUI.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Multiple Component Event-Related Potential (mcERP) Estimation

We show how model-based estimation of the neural sources responsible for transient neuroelectric signals can be improved by the analysis of single trial data. Previously, we showed that a multiple component event-related potential (mcERP) algorithm can extract the responses of individual sources from recordings of a mixture of multiple, possibly interacting, neural ensembles. McERP also estimated single-trial amplitudes and onset latencies, thus allowing more accurate estimation of ongoing neural activity during an experimental trial. The mcERP algorithm is related to informax independent component analysis (ICA); however, the underlying signal model is more physiologically realistic in that a component is modeled as a stereotypic waveshape varying both in amplitude and onset latency from trial to trial. The result is a model that reflects quantities of interest to the neuroscientist. Here we demonstrate that the mcERP algorithm provides more accurate results than more traditional methods such as factor analysis and the more recent ICA. Whereas factor analysis assumes the sources are orthogonal and ICA assumes the sources are statistically independent, the mcERP algorithm makes no such assumptions thus allowing investigators to examine interactions among components by estimating the properties of single-trial responses.

Knuth, K. H.↗

Radiation-induced hypothyroidism in patients with breast cancer: a retrospective analysis of 243 cases

This study aims to estimate the incidence of hypothyroidism (HT) and to evaluate the predictors affecting the development of HT after radiotherapy (RT) for breast cancer, with a focus on radiation dose-volumetric parameters. Between 2009 and 2015, 243 patients undergoing RT for breast cancer were retrospectively analyzed. Free triiodothyronine (FT3), free thyroxin (FT4), and thyrotropin (TSH) were monitored before and after RT. The relation between the doses to thyroid gland (D{sub mean}, D{sub max}, D{sub min}), percentage of thyroid volume receiving > 10 Gy, 20 Gy, 30 Gy, 40 Gy, and 50 Gy (V10 to V50), absolute thyroid volume, and HT were analyzed. The risk of HT according to radiation fields and the other clinic factors were also evaluated. The median follow-up was 41 (range; 6 to 130) months. Sixty-seven percent of the patients received RT to the breast/chest wall and ipsilateral supraclavicular fossa. Of 243 patients, 51 (21%) were diagnosed with HT. The median time to the onset of HT was 27 (range; 5 to 64) months. There were no significant relationships observed between D{sub min} or V10 to V50 and HT. The surgery type, clinical stage, nodal status, RT field, D{sub mean}, and D{sub max} were statistically significant predictors for HT in univariate analysis. The D{sub mean} was the only prognostic factor predicting HT in multivariate analysis, and D{sub mean} > 21 Gy was a threshold value for the evolvement of HT. In this study, we present evidence that postoperative irradiation of patients with breast cancer may frequently lead to HT. Patients who have received RT for breast cancer, especially irradiation on the supraclavicular region, may require thyroid function screening after RT.

61 RADIATION PROTECTION AND DOSIMETRY↗

A novel measurement of marginal Alfven Eigenmode stability during high power auxiliary heating in JET

The interaction of Alfven Eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25-250 kHz and toroidal mode numbers |n| < 20. During the 2019-2020 deuterium campaign, almost 7500 resonances and their frequencies f0, net damping rates \gamma < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient's effect on accessibility and ELM-related noise's impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized Ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ICRH and NBI) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f0, \gamma, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A novel measurement of marginal Alfvén eigenmode stability during high power auxiliary heating in JET

The interaction of Alfvén eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25–250 kHz and toroidal mode numbers |n| < 20. During the 2019–2020 deuterium campaign, almost 7500 resonances and their frequencies f 0 , net damping rates γ < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient’s effect on accessibility and ELM-related noise’s impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ion cyclotron resonance heating and neutral beam injection) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f 0 , γ, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

W Boson Mass Measurement at CDF

This is the closeout report for the grant for experimental research at the energy frontier in high energy physics. The report describes the precise measurement of the W boson mass at the CDF experiment at Fermilab, with an uncertainty of ≈ 12 MeV, using the full dataset of ≈ 9 fb -1 collected by the experiment up to the shutdown of the Tevatron in 2011. In this analysis, the statistical and most of the experimental systematic uncertainties have been reduced by a factor of two compared to the previous measurement with 2.2 fb -1 of CDF data. This research has been the culmination of the PI's track record of producing world-leading measurements of the W boson mass from the Tevatron. The PI performed the first and only measurement to date of the W boson mass using high-rapidity leptons using the D0 endcap calorimeters in Run 1. He has led this measurement in Run 2 at CDF, publishing two world-leading measurements in 2007 and 2012 with total uncertainties of 48 MeV and 19 MeV respectively. The analysis of the final dataset is currently under internal review in CDF. Upon approval of the internal review, the result will be available for public release.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Understanding and improving anode performance in an alkaline membrane electrolyzer using statistical design of experiments

The behavior of the oxygen-evolving positive electrode (i.e. anode) in the anion exchange membrane water electrolyzer (AEMEL) is complex and influenced by several factors. Very few studies have been performed to understand oxygen evolution reaction (OER) behavior by optimizing the individual factors that influence performance. Here, this study highlights the effects of catalyst loading, catalyst selection, porous transport layer (PTL) type and conductive additive content. The influence of each factor is elucidated through a design of experiments (DoE) approach with a full statistical analysis. Electrochemical data, alongside Pareto charts, parametric trends and their mutual interactions are discussed. This DoE approach is also helpful in making useful predictions and discovering new combinations to be tested. The end result was a high-performance AEMEL able to operate at a current density of 1.0 A/cm 2 at 1.80 V with IrOx OER and PtNi hydrogen evolution reaction (HER) catalysts using 0.3 M KOH fed to the anode. Even lower operating voltage was observed with PbRuOx catalyst at the anode, 1.64 V @ 1.0 A/cm 2 , though the cell decay rate was higher. Lastly, a IrOx/PtNi cell was stably operated continuously for 30 days (720 h) at 1.0 A/cm 2 . This study can serve as a guide for optimal electrode design with insights into component-performance compromises, which can aid in making design choices and performing techno-economic analyses.

08 HYDROGEN↗

Identification of Critical Surface Parameters Driving Lectin-Mediated Capture of Bacteria from Solution

Lectin-functional interfaces are useful for isolation of bacteria from solution because they are low-cost and allow nondestructive, reversible capture. This study provides a systematic investigation of physical and chemical surface parameters that influence bacteria capture over lectin-functionalized polymer interfaces and then applies these findings to construct surfaces with significantly enhanced bacteria capture. The designer block copolymer poly(glycidyl methacrylate)-block-poly(vinyldimethyl azlactone) was used as a lectin attachment layer, and lectin coupling into the polymer film through azlactone–lectin coupling reactions was first characterized. Here, experimental parameters including polymer areal chain density, lectin molecular weight, and lectin coupling buffer were systematically varied to identify parameters driving highest azlactone conversions and corresponding lectin surface densities. To introduce physical nanostructures into the attachment layer, nanopillar arrays (NPAs) of varied heights (300 and 2100 nm) were then used to provide an underlying surface template for the functional polymer layer. Capture of Escherichia coli on lectin–polymer surfaces coated over both flat and NPA surfaces was then investigated. For flat polymer interfaces, bacteria were detected on the surface after incubation at a solution concentration of 10 3 cfu/mL, and a corresponding detection limit of 1.7 × 10 3 cfu/mL was quantified. This detection limit was 1 order of magnitude lower than control lectin surfaces functionalized with standard, carbodiimide coupling chemistry. NPA surfaces containing 300 nm tall pillars further improved the detection limit to 2.1 × 10 2 cfu/mL, but also reduced the viability of captured cells. Finally, to investigate the impact of cell surface parameters on capture, we used Agrobacterium tumefaciens cells genetically modified to allow manipulation of exopolysaccharide adhesin production levels. Statistical analysis of surface capture levels revealed that lectin surface density was the primary factor driving capture, as opposed to exopolysaccharide adhesin expression. These findings emphasize the critical importance of the synthetic interface and the development of surfaces that combine high lectin densities with tailored physical features to drive high levels of capture. These insights will aid in design of biofunctional interfaces with physicochemical surface properties favorable for capture and isolation of bacteria cells from solutions.

59 BASIC BIOLOGICAL SCIENCES↗

General Bayesian Inference over the Stiefel Manifold via the Givens Transform

We introduce here an approach based on the Givens representation for posterior inference in statistical models with orthogonal matrix parameters, such as factor models and probabilistic principal component analysis (PPCA). We show how the Givens representation can be used to develop practical methods for transforming densities over the Stiefel manifold into densities over subsets of Euclidean space. We demonstrate how to deal with issues arising from the topology of the Stiefel manifold and how to inexpensively compute the change-of-measure terms. We introduce an auxiliary parameter approach that limits the impact of topological issues. We provide both analysis of our methods and numerical examples demonstrating the effectiveness of the approach. We also discuss how our Givens representation can be used to define general classes of distributions over the space of orthogonal matrices. We then give demonstrations on several examples showing how the Givens approach performs in practice in comparison with other methods.

97 MATHEMATICS AND COMPUTING↗

Effects of miso- and mesoscale obstructions on PAM winds obtained during project NIMROD

Data from 27 PAM (Portable Automated Mesonet) stations, operational as a phase of project NIMROD (Northern Illinois Meteorological Research on Downburst), are presented. It was found that PAM-measured winds are influenced by the mesoscale obstruction of the Chicago metropolitan area, as well as by the misoscale obstruction of identified trees and buildings. The mesoscale obstruction was estimated within the range of near zero to 50%, increasing toward the city limits, while the misoscale obstruction was estimated as being as large as 58% near obstructing trees which were empirically calculated to cause a wind speed deficit 50-80 times their height. Despite a statistical analysis based on one-million PAM winds, wind speed and stability transmission factors could not be accurately calculated; thus, in order to calculate the airflow free from obstacle, PAM-measured winds must be corrected.

Fujita, T. T.↗

Radiation Dose Escalation in Esophageal Cancer Revisited: A Contemporary Analysis of the National Cancer Data Base, 2004 to 2012

Purpose: To evaluate the effect of radiation dose escalation on overall survival (OS) for patients with nonmetastatic esophageal cancer treated with concurrent radiation and chemotherapy. Methods and Materials: Patients diagnosed with stage I to III esophageal cancer treated from 2004 to 2012 were identified from the National Cancer Data Base. Patients who received concurrent radiation and chemotherapy with radiation doses of ≥50 Gy and did not undergo surgery were included. OS was compared using Cox proportional hazards regression and propensity score matching. Results: A total of 6854 patients were included; 3821 (55.7%) received 50 to 50.4 Gy and 3033 (44.3%) received doses >50.4 Gy. Univariate analysis revealed no significant difference in OS between patients receiving 50 to 50.4 Gy and those receiving >50.4 Gy (P=.53). The dose analysis, binned as 50 to 50.4, 51 to 54, 55 to 60, and >60 Gy, revealed no appreciable difference in OS within any group compared with 50 to 50.4 Gy. Subgroup analyses investigating the effect of dose escalation by histologic type and in the setting of intensity modulated radiation therapy also failed to reveal a benefit. Propensity score matching confirmed the absence of a statistically significant difference in OS among the dose levels. The factors associated with improved OS on multivariable analysis included female sex, lower Charlson-Deyo comorbidity score, private insurance, cervical/upper esophagus location, squamous cell histologic type, lower T stage, and node-negative status (P<.01 for all analyses). Conclusions: In this large national cohort, dose escalation >50.4 Gy did not result in improved OS among patients with stage I to III esophageal cancer treated with definitive concurrent radiation and chemotherapy. These data suggest that despite advanced contemporary treatment techniques, OS for patients with esophageal cancer remains unaltered by escalation of radiation dose >50.4 Gy, consistent with the results of the INT-0123 trial. Furthermore, these data highlight that many radiation oncologists have not embraced the concept that dose escalation does not improve OS. Although local control, not investigated in the present study, might benefit from dose escalation, novel therapies are needed to improve the OS of patients with esophageal cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Power spectrum and form factor in random diagonal matrices and integrable billiards

Triggered by a controversy surrounding a universal behavior of the power spectrum in quantum systems exhibiting regular classical dynamics, we focus on a model of random diagonal matrices (RDM), often associated with the Poisson spectral universality class, and examine how the power spectrum and the form factor get affected by two-sided truncations of RDM spectra. Having developed a nonperturbative description of both statistics, we perform their detailed asymptotic analysis to demonstrate explicitly how a traditional assumption (lying at the heart of the controversy) – that the power spectrum is merely determined by the spectral form factor – breaks down for truncated spectra. This observation has important consequences as we further argue that bounded quantum systems with integrable classical dynamics are described by heavily truncated rather than complete RDM spectra. High-precision numerical simulations of semicircular and irrational rectangular billiards lend independent support to these conclusions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hierarchical semi-Markov models with duration-aware dynamics for activity sequences

Residential electricity demand at granular scales is driven by what people do and for how long. Accurately forecasting this demand for applications like microgrid management and demand response therefore requires generative models for activities that can produce realistic daily activity sequences, capturing both the timing and duration of human behavior. This paper develops a generative model of human activity sequences using nationally representative time-use diaries at a 10-min resolution. We use this model to quantify which demographic factors are most critical for improving predictive performance. We propose a hierarchical semi-Markov framework that addresses two key modeling challenges. First, a time-inhomogeneous Markov router learns the patterns of “which activity comes next.” Second, a semi-Markov hazard component explicitly models activity durations, capturing “how long” activities realistically last. To ensure statistical stability when data are sparse, the model pools information across related demographic groups and time blocks. The entire framework is trained and evaluated using survey design weights to ensure our findings are representative of the U.S. population. On a held-out test set, we demonstrate that explicitly modeling durations with the hazard component provides a substantial and statistically significant improvement over purely Markovian models. Furthermore, our analysis reveals a clear hierarchy of demographic factors: Sex, Day-Type, and Household Size provide the largest predictive gains, while Region and Season, though important for energy calculations, contribute little to predicting the activity sequence itself. The result is an interpretable and robust generator of synthetic activity traces, providing a high-fidelity foundation for downstream energy systems modeling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Environmental risk factors associated with low birth weight: The case study of the Haifa Bay Area in Israel

Highlights: • Environmental factors affecting low birth weight (LBW) rate are investigated. • Analysis is conducted by alternative techniques of geo-statistical analysis. • Proximity to petrochemical industries is found to increase LBW rate. • Seashore proximity and elevation above the sea level reduce the LBW rate. • We attribute the latter to climate-moderating effect of elevated seashore locations. Low birth weight (LBW) is known to be associated with infant mortality and postnatal health complications. Previous studies revealed strong relationships between LBW rate and several socio-demographic factors, including ethnicity, maternal age, and family income. However, studies of association between LBW rate and environmental risk factors remain infrequent.

54 ENVIRONMENTAL SCIENCES↗

The role of social support on midwestern farmers’ willingness to grow perennial bioenergy crops

The lack of farmers' willingness to grow perennial bioenergy crops (PBCs) presents a critical barrier to the emergence of cellulosic biofuel production. The willingness relies on a complex network of economic, environmental, and social drivers, among which the influence of social factors (e.g., the influence of neighborhood, community, and communication) is less understood. This study addresses this knowledge gap via a survey analysis of midwestern farmers. The survey data are analyzed through ordinary least square regression and structural equation model, which together investigate the individual and interactive impacts of multiple factors on farmers' decisions to adopt PBCs. Based on a farm-scale analysis, six statistically significant predictors of farmer willingness to grow PBCs are identified: perception of PBCs' environment benefits, education level, willingness to take risks, familiarity with PBCs, portion of peers already growing PBCs, and support of biorefineries locating in the local community. Among these, the latter three predictors are social support variables. It is found that familiarity with the crops is the most significant predictor of willingness; familiarity is also an important intermediate variable that mediates the influence of many other predictors. In addition, peer adoption can both directly and indirectly affect willingness via its influence on familiarity. Furthermore, these findings suggest that it is a pressing need to improve farmers’ knowledge of PBCs to promote the adoption of such crops.

09 BIOMASS FUELS↗

Probing cosmic velocities with the pairwise kinematic Sunyaev-Zel’dovich signal in DESI Bright Galaxy Sample DR1 and ACT DR6

We present a measurement of the pairwise kinematic Sunyaev-Zel’dovich (kSZ) signal using the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Sample (BGS) Data Release 1 (DR1) galaxy sample overlapping with the Atacama Cosmology Telescope (ACT) CMB temperature map. Our analysis makes use of 1.6 million galaxies with stellar masses log⁡ 𝑀 ⋆ /𝑀 ⊙ >10, and we explore measurements across a range of aperture sizes (2.1′ <𝜃 ap <3.5′) and stellar mass selections. This statistic directly probes the velocity field of the large-scale structure, a unique observable of cosmic dynamics and modified gravity. In particular, at low redshifts, this quantity is especially interesting, as deviations from General Relativity are expected to be largest. Notably, our result represents the highest-significance low-redshift (𝑧 ∼ 0.3) detection of the kSZ pairwise effect yet. In our most optimal configuration (𝜃 ap =3.3′, log⁡ 𝑀 ⋆ >11), we achieve a 5⁢𝜎 detection. Assuming that an estimate of the optical depth and galaxy bias of the sample exists via e.g., external observables, this measurement constrains the fundamental cosmological combination 𝐻 0 ⁡𝑓⁡𝜎$^2_8$. A key challenge is the degeneracy with the galaxy optical depth. We address this by combining CMB lensing, which allows us to infer the halo mass and galaxy population properties, with hydrodynamical simulation estimates of the mean optical depth, $\bar{𝜏}$ . We stress that this is a proof-of-concept analysis; with BGS DR2 data we expect to improve the statistical precision by roughly a factor of two, paving the way toward robust tests of modified gravity with kSZ-informed velocity-field measurements at low redshift.

Hadzhiyska, Boryana [Institute of Astronomy; Kavli↗

A Survey of Singular Value Decomposition Methods for Distributed Tall/Skinny Data

The Singular Value Decomposition (SVD) is one of the most important matrix factorizations, enjoying a wide variety of applications across numerous application domains. In statistics and data analysis, the common applications of SVD inclue Principal Components Analysis (PCA) and regression. Usually these applications arise on data that has far more rows than columns, so-called "tall/skinny" matrices. In the big data analytics context, this may take the form of hundreds of millions to billions of rows with only a few hundred columns. There is a need, therefore, for fast, accurate, and scalable tall/skinny SVD implementations which can fully utilize modern computing resources. To that end, we present a survey of three different algorithms for computing the SVD for these kinds of tall/skinny data layouts using MPI for communication. We contextualize these with common big data analytics techniques. Finally, we present both CPU and GPU timing results from the Summit supercomputer, and discuss possible alternative approaches.

Schmidt, Drew↗