Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maximum Likelihood”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Simulator for Neyer Tests of Explosives

Explosives and explosive devices such as detonators are typically tested by applying a range of stimuli such as voltage or mechanical shock, and recording binary “detonated/did not detonate” responses. These are analyzed using maximum likelihood or generalized linear models to provide estimates of quantities such as the all-fire and no-fire points. Given that the true threshold for detonation is unknown a priori , sequential design methods are typically used to optimize the set of test points. One popular method, implemented in commercial software, is Neyer’s algorithm. To support simulation and experimental design, we have developed code in the R programming language to duplicate the functions of the Neyer software. We provide code for the simulator along with a description and examples of usage.

42 ENGINEERING↗
DOE OSTI2026

Poisson-response Tensor-on-Tensor Regression and Applications

We introduce Poisson-response tensor-on-tensor regression (PToTR), a novel regression framework designed to handle tensor responses composed element-wise of random Poisson-distributed counts. Tensors, or multi-dimensional arrays, composed of counts are common data in fields such as inter national relations, social networks, epidemiology, and medical imaging, where events occur across multiple dimensions like time, location, and dyads. PToTR accommodates such tensor responses alongside tensor covariates, providing a versatile tool for multi dimensional data analysis. We propose algorithms for maximum likelihood estimation under a canonical polyadic (CP) structure on the regression coefficient tensor that satisfy the positivity of Poisson parameters and then provide an initial theoretical error analysis for PToTR estimators. We also demonstrate the utility of PToTR through three concrete applications: longitudinal data analysis of the Integrated Crisis Early Warning System database, positron emission tomography (PET) image reconstruction, and change-point detection of communication patterns in longitudinal dyadic data. These applications highlight the versatility of PToTR in addressing complex, structured count data across various domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Probabilistic Data Association for Orbital-Element Estimation Using Multistage Expectation–Maximization

Tracking space objects is important for managing space traffic and predicting collisions, but is difficult in part due to data association and orbit model uncertainty. Expectation–maximization (EM) is a commonly used tracking method that has not been widely considered for tracking space objects. The technique consists of iteratively computing data association probabilities with a set of current element estimates, and updating estimates of the elements by solving a nonlinear weighted least-squares regression problem where the weights are the data association probabilities. In this paper, we demonstrate the use of EM for probabilistic data association and orbital-element estimation by applying the technique to simulated data from two angles-only tracking scenarios. In both scenarios, EM provides correct data associations and accurate maximum likelihood estimates of orbital elements. One scenario considers tracking a single object in clutter and quantifies the improvement of the orbital-element estimates and data associations as the detection probability increases. However, standard application of EM requires knowing the number of objects or may fail when a large number of objects are present. To address these issues, this paper employs a multistage version of EM that is applicable when there are a large and possibly unknown number of objects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

On the uncertainty of long-period return values of extreme daily precipitation

Methods for calculating return values of extreme precipitation and their uncertainty are compared using daily precipitation rates over the Western U.S. and Southwestern Canada from a large ensemble of climate model simulations. The roles of return-value estimation procedures and sample size in uncertainty are evaluated for various return periods. We compare two different generalized extreme value (GEV) parameter estimation techniques, namely L-moments and maximum likelihood (MLE), as well as empirical techniques. Even for very large datasets, confidence intervals calculated using GEV techniques are narrower than those calculated using empirical methods. Furthermore, the more efficient L-moments parameter estimation techniques result in narrower confidence intervals than MLE parameter estimation techniques at small sample sizes, but similar best estimates. It should be noted that we do not claim that either parameter fitting technique is better calibrated than the other to estimate long period return values. While a non-stationary MLE methodology is readily available to estimate GEV parameters, it is not for the L-moments method. Comparison of uncertainty quantification methods are found to yield significantly different estimates for small sample sizes but converge to similar results as sample size increases. Finally, practical recommendations about the length and size of climate model ensemble simulations and the choice of statistical methods to robustly estimate long period return values of extreme daily precipitation statistics and quantify their uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Molecular evolution of human coronavirus-NL63, -229E, -HKU1 and -OC43 in hospitalized children in China

Human coronaviruses (HCoVs) HCoV-NL63, HCoV-229E, HCoV-HKU1 and HCoV-OC43 have been circulated in the human population worldwide, and they are associated with a broad range of respiratory diseases with varying severity. However, there are neither effective therapeutic drugs nor licensed vaccines available for the treatment and prevention of infections by the four HCoVs. In this study, we collected nasopharyngeal aspirates of children hospitalized for respiratory tract infection in China during 2014–2018 and conducted next-generation sequencing. Sequences of four HCoVs were then selected for an in-depth analysis. Genome sequences of 2 HCoV-NL63, 8 HCoV-229E, 2 HCoV-HKU1, and 6 HCoV-OC43 were obtained. Based on the full-length S gene, a strong temporal signal was found in HCoV-229E and the molecular evolutionary rate was 6 × 10 −4 substitutions/site/year. Based on the maximum-likelihood (ML) phylogenetic tree of complete S gene, we designated H78 as a new sub-genotype C2 of HCoV-HKU1, and the obtained P43 sequence was grouped into the reported novel genotype K of HCoV-OC43 circulating in Guangzhou, China. Based on the complete genome, potential recombination events were found to occur as two phenomena, namely intraspecies and interspecies. Moreover, we observed two amino acid substitutions in the S1 subunit of obtained HCoV-NL63 (G534V) and HCoV-HKU1 (H512R), while residues 534 and 512 are important for the binding of angiotensin-converting enzyme 2 and neutralizing antibodies, respectively. Our findings might provide a clue for the molecular evolution of the four HCoVs and help in the early diagnosis, treatment and prevention of broad-spectrum HCoV infection.

Shao, Nan↗
DOE OSTI2022

Identification of a new circulating recombinant form of human immunodeficiency virus type 1, CRF124_cpx involving subtypes A, G, H, and CRF27_cpx in Angola

Angola, located in Central Africa, has around 320,000 (270,000–380,000) people living with human immunodeficiency virus (HIV)/AIDS, equivalent to 1% of the country’s population at the end of 2021. A previous study conducted in 2012, using Angolan samples collected between 2008 and 2010 revealed a high prevalence of HIV-1 recombinants, around 42% of sequences, with 21% showing the same UH profile in partial pol region which were grouped into a monophyletic cluster with high bootstrap support. Thus, the objective of the present work was to obtain complete genomes of those sequences and characterize them, aiming at a description of a new circulating recombinant form (CRF). Whole blood from nine HIV-1 UH pol -infected individuals had their genomic DNA extracted, and nested PCR was used to amplify seven overlapping fragments targeting the full-length HIV-1 genome. The final classification was based on maximum likelihood trees, and recombination analyses were performed using a bootscan from the Simplot program. BLAST and Los Alamos Database inspections were used to search other similar H-like pol sequences. Complete genome amplification was possible for three samples, partial genomes were obtained for the other three, and only pol was available for the remaining three sequences. Bootscan analysis of the two whole-genome and three partial genome sequences retrieved from people living with HIV/AIDS (PLHIVA) without epidemiological linkage showed the same complex recombination profile involving HIV-1 subtypes A/G/H/CRF27_cpx, with a total of six recombinant breakpoints, aiming to classify a new HIV-1 CRF124_cpx. We found no other full-length HIV-1 genomes with the same mosaic profile; however, we identified 33 partial pol sequences, mainly sampled from Angola between 2001 to 2019, with the same H-like profile. Bayesian analysis of H and H-like pol sequences indicates that CRF124_cpx probably originated in Angola at mid-1970s, indicating that this CRF has been circulating in the country for a long time. In summary, our study describes a new CRF circulating principally in Angola and highlights the importance of continuing molecular surveillance studies, especially in countries with high molecular diversity of HIV.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

classLog: Logistic regression for the classification of genetic sequences

Introduction Sequencing and phylogenetic classification have become a common task in human and animal diagnostic laboratories. It is routine to sequence pathogens to identify genetic variations of diagnostic significance and to use these data in realtime genomic contact tracing and surveillance. Under this paradigm, unprecedented volumes of data are generated that require rapid analysis to provide meaningful inference. Methods We present a machine learning logistic regression pipeline that can assign classifications to genetic sequence data. The pipeline implements an intuitive and customizable approach to developing a trained prediction model that runs in linear time complexity, generating accurate output rapidly, even with incomplete data. Our approach was benchmarked against porcine respiratory and reproductive syndrome virus (PRRSv) and swine H1 influenza A virus (IAV) datasets. Trained classifiers were tested against sequences and simulated datasets that artificially degraded sequence quality at 0, 10, 20, 30, and 40%. Results When applied to a poor-quality sequence data, the classifier achieved between >85% to 95% accuracy for the PRRSv and the swine H1 IAV HA dataset and this increased to near perfect accuracy when using the full dataset. The model also identifies amino acid positions used to determine genetic clade identity through a feature selection ranking within the model. These positions can be mapped onto a maximum-likelihood phylogenetic tree, allowing for the inference of clade defining mutations. Discussion Our approach is implemented as a python package with code available at https://github.com/flu-crew/classLog .

Zeller, Michael A.↗
DOE OSTI2022

Transcriptome Sequencing of Agave angustifolia Reveals Conservation and Diversification in the Expression of Cinnamyl Alcohol Dehydrogenase Genes in Agave Species

Agave angustifolia is an important crassulacean acid metabolism plant species, with wide applications in beverage and sisal fiber production. In this study, we carried out a transcriptome analysis of A. angustifolia leaves, generating a total of 58,482,436 clean reads through Illumina paired-end sequencing. De novo transcriptome assembly generated 67,314 unigenes, with about half of them having homologs in four public databases. In the Nr database, Asparagus officinalis was shown to be most closely related to agave, and the early angiosperm Amborella trichopoda was selected as an outgroup species. We further identified five, two, seven, seven, seven, six and six CAD genes in asparagus, amborella, A. deserti, A. tequilana, A. americana, A. hybrid H11648 and A. angustifolia, respectively. The maximum likelihood phylogenetic tree revealed the species-specific expansion of CAD genes in arabidopsis, rice and agave. The expression analysis indicated the conservatively expressed CAD1/2/4/6, providing candidate targets for manipulation to improve lignin traits. The species-specific expression of CAD3/5/7 indicates the existence of different regulatory mechanisms controlling the expression of these genes in agave species. This study presents the first transcriptome dataset of A. angustifolia, facilitating future studies on lignin biosynthesis in agave.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Comparative Study of Wind Energy Potential Estimation Methods for Wind Sites in Togo and Benin (West Sub-Saharan Africa)

The characterization of wind speed distribution and the optimal assessment of wind energy potential are critical factors in selecting a suitable site for wind power plants (WPP). The Weibull distribution law has been used extensively to analyze the wind characteristics of candidate WPP sites, and to estimate the available and deliverable energy. This paper presents a comparative study of five wind energy resource assessment methods as they applied to the context of wind sites in West Sub-Saharan Africa. We investigated three numerical approaches, namely, the adaptive neuro-fuzzy inference system (ANFIS), the multilayer perceptron method (MLP), and support vector regression (SVR), to derive the distribution law of wind speeds and to optimally quantify the corresponding wind energy potential. Next, we compared these three approaches to two well-known Weibull distribution law-based methods: the empirical method of Justus (EMJ) and the maximum likelihood method (MLM). Case study results indicated that the neural network-based methods, ANFIS and MLP, yielded the most accurate distribution fits and wind energy potential estimates, and consequently, are the most recommended methods for the wind sites in Togo and Benin. The orders of magnitude of the root mean squared error (RMSE) in estimating the recoverable energy using ANFIS were, respectively, 10-4 and 10-5 for Lomé and Cotonou, while MLP achieved an RMSE order of magnitude of 10-3 for both sites.

17 WIND ENERGY↗
DOE OSTI2021

High-Resolution Scanning Coded-Mask-Based X-ray Multi-Contrast Imaging and Tomography

Near-field X-ray speckle tracking has been used in phase-contrast imaging and tomography as an emerging technique, providing higher contrast images than traditional absorption radiography. Most reported methods use sandpaper or membrane filters as speckle generators and digital image cross-correlation for phase reconstruction, which has either limited resolution or requires a large number of position scanning steps. Recently, we have proposed a novel coded-mask-based multi-contrast imaging (CMMI) technique for single-shot measurement with superior performance in efficiency and resolution compared with other single-shot methods. We present here a scanning CMMI method for the ultimate imaging resolution and phase sensitivity by using a coded mask as a high-contrast speckle generator, the flexible scanning mode, the adaption of advanced maximum-likelihood optimization to scanning data, and the multi-resolution analysis. Scanning CMMI can outperform other speckle-based imaging methods, such as X-ray speckle vector tracking, providing higher quality absorption, phase, and dark-field images with fewer scanning steps. Scanning CMMI is also successfully demonstrated in multi-contrast tomography, showing great potentials in high-resolution full-field imaging applications, such as in vivo biomedical imaging.

Qiao, Zhi (ORCID:0000000233928732)↗
DOE OSTI2021

Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances

Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.

Shayesteh Zadeh, Armin (ORCID:0000000299896710)↗
DOE OSTI2023

Comprehensive Genetic Characterization of Four Novel HIV-1 Circulating Recombinant Forms (CRF129_56G, CRF130_A1B, CRF131_A1B, and CRF138_cpx): Insights from Molecular Epidemiology in Cyprus

Molecular investigations of the HIV-1 pol region (2253–5250 in the HXB2 genome) were conducted on sequences obtained from 331 individuals infected with HIV-1 in Cyprus between 2017 and 2021. This study unveiled four distinct HIV-1 putative transmission clusters, encompassing 19 previously unidentified HIV-1 recombinants. These recombinants, each comprising eight, three, four, and four sequences, respectively, did not align with previously established Circulating Recombinant Forms (CRFs). To characterize these novel HIV-1 recombinants, near-full-length genome sequences were successfully obtained for 16 of the 19 recombinants (790–8795 in the HXB2 genome) using an in-house-developed RT-PCR assay. Phylogenetic analyses, employing MEGAX and Cluster-Picker, along with confirmatory neighbor-joining tree analyses of subregions, were conducted to identify distinct clusters and determine subtypes. The uniqueness of the HIV-1 recombinants was evident in their exclusive clustering within generated maximum likelihood trees. Recombination analyses highlighted the distinct chimeric nature of these recombinants, with consistent mosaic patterns observed across all sequences within each of the four putative transmission clusters. Conclusive genetic characterization identified four novel HIV-1 CRFs: CRF129_56G, CRF130_A1B, CRF131_A1B, and CRF138_cpx. CRF129_56G exhibited two recombination breakpoints and three fragments of subtypes CRF56_cpx and G. Both CRF130_A1B and CRF131_A1B featured seven recombination breakpoints and eight fragments of subtypes A1 and B. CRF138_cpx displayed five recombination breakpoints and six fragments of subtypes CRF22_01A1 and F2, along with an unclassified fragment. Additional BLAST analyses identified a Unique Recombinant Form (URF) of CRF138_cpx with three additional recombination sites, involving subtype F2, a fragment of unknown subtype origin, and CRF138_cpx. Post-identification, all putative transmission clusters remained active, with CRF130_A1B, CRF131_A1B, and CRF138_cpx clusters exhibiting further growth. Furthermore, international connections were identified through BLAST analyses, linking one sequence from the USA to the CRF130_A1B strain, and three sequences from Belgium and Cameroon to the CRF138_cpx strain. This study contributes valuable insights into the dynamic landscape of HIV-1 diversity and transmission patterns, emphasizing the need for ongoing molecular surveillance and global collaboration in tracking emerging viral variants.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Phylogeny and character evolution in the Dacrymycetes , and systematics of Unilacrymaceae and Dacryonaemataceae fam. nov.

We present a multilocus phylogeny of the class Dacrymycetes, based on data from the 18S, ITS, 28S, RPB1,RPB2,TEF-1α, 12S, andATP6DNA regions, with c. 90 species including the types of most currently accepted genera. A variety of methodological approaches was used to infer phylogenetic relationships among the Dacrymycetes, from a supermatrix strategy using maximum likelihood and Bayesian inference on a concatenated dataset, to coalescence-based calculations, such as quartet-based summary methods of independent single-locus trees, and Bayesian integration of single-locus trees into a species tree under the multispecies coalescent. We evaluate for the first time the taxonomic usefulness of some cytological phenotypic characters, i.e., vacuolar contents (vacuolar bodies and lipid bodies), number of nuclei of recently discharged basidiospores, and pigments, with especial emphasis on carotenoids. These characters, along with several others traditionally used for the taxonomy of this group (basidium shape, presence and morphology of clamp connections, morphology of the terminal cells of cortical/marginal hyphae, presence and degree of ramification of the hyphidia), are mapped on the resulting phylogenies and their evolution through the class Dacrymycetes discussed. Our analyses reveal five lineages that putatively represent five different families, four of which are accepted and named. Three out of these four lineages correspond to previously circumscribed and published families (Cerinomycetaceae, Dacrymycetaceae, and Unilacrymaceae), and one is proposed as the new family Dacryonaemataceae. Provisionally, only a single order, Dacrymycetales, is accepted with in the class. Furthermore, the systematics of the two smallest families, Dacryonaemataceae and Unilacrymaceae, are investigated to the species level, using coalescence-based species delimitation on multilocus DNA data, and a detailed morphological study including morphometric analyses of the basidiospores. Three species are accepted in Dacryonaema, the type, Da. rufum, the newly combined Da. Macnabbii (basionym Dacrymyces macnabbii), and a new species named Da. macrosporum. Two species are accepted in Unilacryma, the new U. bispora, and the type, U. unispora, the latter treated in a broad sense pending improved sampling across the Holarctic.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Galactic Chemical Evolution of Short-lived Radioactive Isotopes Produced by Explosive Nucleosynthesis: 60 Fe and 53 Mn

Several short-lived radionuclides (SLRs) are know to have existed in the early solar system (ESS). These species, which typically decay with half-lives of the order of a few million years, can be used to probe the timescales of events preceding the birth of the Sun. We investigate the ESS origin of 53 Mn, produced by core-collapse supernovae (CCSNe) and Type Ia supernovae (SNe Ia), and 60 Fe, produced exclusively by CCSNe. We model the evolution of the radioactive-to-stable abundance ratios of these SLRs with a Galactic chemical evolution (GCE) framework accounting for different SN yields, SN Ia delay times, and other Galactic features (K). A further set of models is calculated assuming that SNe Ia did not contribute any 53 Mn to the ESS. The predicted ratios are compared to meteoritic ratios to derive a distribution of solar isolation times that includes uncertainties due to stochastic chemical enrichment and measurements of the ESS values. The isolation times are then compared to those of 107 Pd and 182 Hf calculated in previous work. A self-consistent solution can be found within the current uncertainties, especially when using the GCE setups with K = 1.6 and 2.3, although the maximum likelihood for the 60 Fe distribution is typically ∼4–5 Myr shorter than for 53 Mn. The predicted 60 Fe/ 53 Mn ratio, instead, is completely inconsistent with the ESS value; this could be resolved using a larger fraction of faint CCSNe than usually considered in GCE models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

TARDIS. I. A Constrained Reconstruction Approach to Modeling the z ~ 2.5 Cosmic Web Probed by Lyα Forest Tomography

Recent Lyα forest tomography measurements of the intergalactic medium (IGM) have revealed a wealth of cosmic structures at high redshift (z ~ 2.5). In this work, we present the Tomographic Absorption Reconstruction and Density Inference Scheme (TARDIS), a new chronocosmographic analysis tool for understanding the formation and evolution of these observed structures. Here, we use maximum likelihood techniques with a fast nonlinear gravitational model to reconstruct the initial density field of the observed regions. We find that TARDIS allows accurate reconstruction of smaller-scale structures than standard Wiener-filtering techniques. Applying this technique to mock Lyα forest data sets that simulate ongoing and future surveys such as CLAMATO, Subaru PFS, or the ELTs, we are able to infer the underlying matter density field at observed redshift and classify the cosmic web structures. We find good agreement with the underlying truth in both the characteristic eigenvalues and eigenvectors of the pseudo-deformation tensor, with the eigenvalues inferred from 30 m class telescopes correlated at r = 0.95 relative to the truth. As an output of this method, we are able to further evolve the inferred structures to late time (z = 0) and also track the trajectories of coeval z = 2.5 galaxies to their z = 0 cosmic web environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

IceCube Search for Neutrinos Coincident with Gravitational Wave Events from LIGO/Virgo Run O3

Using data from the IceCube Neutrino Observatory, we searched for high-energy neutrino emission from the gravitational-wave events detected by the advanced LIGO and Virgo detectors during their third observing run. We did a low-latency follow-up on the public candidate events released during the detectors' third observing run and an archival search on the 80 confident events reported in the GWTC-2.1 and GWTC-3 catalogs. An extended search was also conducted for neutrino emission on longer timescales from neutron star containing mergers. Follow-up searches on the candidate optical counterpart of GW190521 were also conducted. We used two methods; an unbinned maximum likelihood analysis and a Bayesian analysis using astrophysical priors, both of which were previously used to search for high-energy neutrino emission from gravitational-wave events. No significant neutrino emission was observed by any analysis, and upper limits were placed on the time-integrated neutrino flux as well as the total isotropic equivalent energy emitted in high-energy neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Search for IceCube Sub-TeV Neutrinos Correlated with Gravitational-wave Events Detected By LIGO/Virgo

The LIGO/Virgo collaboration published the catalogs GWTC-1, GWTC-2.1, and GWTC-3 containing candidate gravitational-wave (GW) events detected during its runs O1, O2, and O3. These GW events can be possible sites of neutrino emission. In this paper, we present a search for neutrino counterparts of 90 GW candidates using IceCube DeepCore, the low-energy infill array of the IceCube Neutrino Observatory. The search is conducted using an unbinned maximum likelihood method, within a time window of 1000 s, and uses the spatial and timing information from the GW events. The neutrinos used for the search have energies ranging from a few GeV to several tens of TeV. We do not find any significant emission of neutrinos, and place upper limits on the flux and the isotropic-equivalent energy emitted in low-energy neutrinos. We also conduct a binomial test to search for source populations potentially contributing to neutrino emission. We report a nondetection of a significant neutrino-source population with this test.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗