Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “clustering algorithm”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71

Real‐time XFEL data analysis at SLAC and NERSC: A trial run of nascent exascale experimental data analysis

X‐ray scattering experiments using free electron lasers (XFELs) are a powerful tool to determine the molecular structure and function of unknown samples (such as COVID‐19 viral proteins). XFEL experiments are a challenge to computing in two ways: (i) due to the high cost of running XFELs, a fast turnaround time from data acquisition to data analysis is essential to make informed decisions on experimental protocols; (ii) data‐collection rates are growing exponentially, requiring new scalable algorithms. Here we report our experiences analyzing data from two experiments at the Linac Coherent Light Source (LCLS) during September 2020. Raw data were analyzed on NERSC's Cori XC40 system, using the Superfacility paradigm: our workflow automatically moves raw data between LCLS and NERSC, where it is analyzed using the software package CCTBX. We achieved real time data analysis with a turnaround time from data acquisition to full molecular reconstruction in as little as 10 min—sufficient time for the experiment's operators to make informed decisions. By hosting the data analysis on Cori, and by automating LCLS‐NERSC interoperability, we achieved a data analysis rate which matches the data acquisition rate. Completing data analysis within 10 min is a first for XFEL experiments and an important milestone if we are to keep up with data‐collection trends.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Numerical experiments with model monophyletic and paraphyletic taxa

The problem of how accurately paraphyletic taxa versus monophyletic (i.e., holophyletic) groups (clades) capture underlying species patterns of diversity and extinction is explored with Monte Carlo simulations. Phylogenies are modeled as stochastic trees. Paraphyletic taxa are defined in an arbitrary manner by randomly choosing progenitors and clustering all descendants not belonging to other taxa. These taxa are then examined to determine which are clades, and the remaining paraphyletic groups are dissected to discover monophyletic subgroups. Comparisons of diversity patterns and extinction rates between modeled taxa and lineages indicate that paraphyletic groups can adequately capture lineage information under a variety of conditions of diversification and mass extinction. This suggests that these groups constitute more than mere "taxonomic noise" in this context. But, strictly monophyletic groups perform somewhat better, especially with regard to mass extinctions. However, when low levels of paleontologic sampling are simulated, the veracity of clades deteriorates, especially with respect to diversity, and modeled paraphyletic taxa often capture more information about underlying lineages. Thus, for studies of diversity and taxic evolution in the fossil record, traditional paleontologic genera and families need not be rejected in favor of cladistically-defined taxa.

Non-NASA Center↗

The Power of Many: An Ensemble Approach to Spectral Similarity

Quantifying the similarity between two mass spectra─a known reference mass spectrum and an unidentified sample mass spectrum─is at the heart of compound identification workflows in gas chromatography–mass spectrometry (GC-MS). The reference spectrum most like the sample is assigned as its identification (provided some quantitative similarity threshold is met, e.g., 80%) and thus accurately measuring similarity is essential. Significant research has gone toward developing metrics for this purpose, each of which has attempted to improve upon existing methods by incorporating GC-MS-specific information (e.g., peak ratios or retention times) or adopting various statistical and algorithmic frameworks. While this active development has led to a plethora of similarity metrics with demonstrated value across different contexts, the unfortunate consequence has been confusion surrounding which metric should be used as a global standard. No such metric is currently accepted as the standard method because different metrics have demonstrated optimal performance in different contexts. In this work, we propose an ensemble approach to spectral similarity scoring that combines the collective information from across existing similarity metrics to form an improved, globally representative similarity metric as a step toward establishing a global standard method. In conclusion, the resulting ensemble metrics are evaluated on over 88,000 spectra of varying complexity and demonstrate improved abilities to accurately rank the correct reference spectrum as the top-matching candidate for a sample relative to the rankings generated by individual similarity scores.

Carbohydrates↗

Acceleration of Graph Neural Network-Based Prediction Models in Chemistry via Co-Design Optimization on Intelligence Processing Units

Atomic structure prediction and associated property calculations are the bedrock of chemical physics. Since high-fidelity ab initio modeling techniques for computing the structure and properties can be prohibitively expensive, this motivates the development of machine-learning (ML) models that make these predictions more efficiently. Training graph neural networks over large atomistic databases introduces unique computational challenges such as the need to process millions of small graphs with variable size and support communication patterns that are distinct from learning over large graphs such as social networks. We demonstrate a novel hardware-software co-design approach to scale up the training of atomistic graph neural networks (GNN) for structure and property prediction. First, to eliminate redundant computation and memory associated with alternative padding techniques and to improve throughput via minimizing communication, we formulate the effective coalescing of the batches of variable-size atomistic graphs as the bin packing problem and introduce a hardware-agnostic algorithm to pack these batches. In addition, we propose hardware-specific optimizations including a planner and vectorization for the gather-scatter operations targeted for Graphcore’s Intelligence Processing Unit (IPU), as well as model-specific optimizations such as merged communication collectives and optimized softplus. Putting these all together, we demonstrate the effectiveness of the proposed co-design approach by providing an implementation of a well-established atomistic GNN on the Graphcore IPUs. We evaluate the training performance on multiple atomistic graph databases with varying degrees of graph counts, sizes and sparsity. Here, we demonstrate that such a co-design approach can reduce the training time of atomistic GNNs and can improve the performance by up to 1.5× compared to the baseline implementation of the model on the IPUs. Additionally, we compare our IPU implementation with a Nvidia GPU-based implementation and show that our atomistic GNN implementation on the IPUs can run 1.8× faster on average compared to the execution time on the GPUs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New constraints on sodium production in globular clusters from the Na 23 ( He 3 , d ) Mg 24 reaction

The star-to-star anticorrelation of sodium and oxygen is a defining feature of globular clusters, but, to date, the astrophysical site responsible for this unique chemical signature remains unknown. Sodium enrichment within these clusters depends sensitively on reaction rate of the sodium destroying reactions 23 Na(p, γ) and 23 Na(p,α). In this paper, we report the results of a 23 Na( 3 He,d) 24 Mg transfer reaction carried out at Triangle Universities Nuclear Laboratory using a 21 MeV 3 He beam. Astrophysically relevant states in 24 Mg between 11 < E x < 12 MeV were studied using high-resolution magnetic spectroscopy, thereby allowing the extraction of excitation energies and spectroscopic factors. Bayesian methods are combined with the distorted wave Born approximation to assign statistically meaningful uncertainties to the extracted spectroscopic factors. For the first time, these uncertainties are propagated through to the estimation of proton partial widths. Our experimental data are used to calculate the reaction rate. The impact of the new rates are investigated using asymptotic giant branch star models. Furthermore, it is found that while the astrophysical conditions still dominate the total uncertainty, intramodel variations on sodium production from the 23 Na(p, γ) and 23 Na(p,α) reaction channels are a lingering source of uncertainty.

20 ≤ A ≤ 38↗

Continuous Filament Network of the Local Universe

Simulated galaxy distributions are suitable for developing filament detection algorithms. However, samples of observed galaxies, being of limited size, cause difficulties that lead to a discontinuous distribution of filaments. We created a new galaxy filament catalog composed of a continuous cosmic web with no lone filaments. The core of our approach is a ridge filter used within the framework of image analysis. We considered galaxies from the HyperLeda database with redshifts 0.02 ≤ z ≤ 0.1, and in the solid angle 120° ≤ R.A. ≤ 240°, 0° ≤ decl. ≤ 60°. We divided the sample into 16 two-dimensional celestial projections with redshift bin Δ z = 0.005, and compared our continuous filament network with a similar recent catalog covering the same region of the sky. We tested our catalog on two application scenarios. First, we compared the distributions of the distances to the nearest filament of various astrophysical sources (Seyfert galaxies and other active galactic nuclei, radio galaxies, low-surface-brightness galaxies, and dwarf galaxies), and found that all source types trace the filaments well, with no systematic differences. Next, among the HyperLeda galaxies, we investigated the dependence of the g - r color distribution on the distance to the nearest filament, and confirmed that early-type galaxies are located on average further from the filaments than late-type ones.

79 ASTRONOMY AND ASTROPHYSICS↗

Fossil group origins

Aims: We analyse the large-scale structure out to 100 Mpc around a sample of 16 confirmed fossil systems using spectroscopic information from the Sloan Digital Sky Survey Data Release 16. Methods: We computed the distance between our fossil groups (FGs) and the centres of filaments and nodes from the literature. We also studied the density of bright galaxies, since this parameter is thought to be a good mass tracers, as well as the projected over-densities of galaxies. Finally, we applied a friends-of-friends (FoF) algorithm to detect virialised structures around our FGs and obtain an estimate of the mass available in their surroundings. Results: We find that FGs are mainly located close to filaments, with a mean distance of 3.7 ± 1.1 R 200 and a minimum distance of 0.05 R 200 . On the other hand, none of our FGs were found close to intersections, with a mean and minimum distance of 19.3 ± 3.6 and 6.1 R 200 , respectively. There is a correlation that indicates FGs at higher redshifts are found in denser regions, when we use bright galaxies as tracers of the mass. At the same time, FGs with the largest magnitude gaps (Δm 12 > 2.5) are found in less dense environments and tend to host (on average) smaller central galaxies. Conclusions: Our results suggest that FGs formed in a peculiar position within the cosmic web, close to filaments and far from nodes, whereby their interaction with the cosmic web itself may be limited. We deduce that FGs with brightest central galaxies (BCGs) that are relatively faint, high values of Δm 12 , and low redshifts could, in fact, be systems that are at the very last stage of their evolution. Moreover, we confirm theoretical predictions that systems with the largest magnitude gap are not massive.

79 ASTRONOMY AND ASTROPHYSICS↗

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Visualizing Temporal Topic Embeddings with a Compass

—Dynamic topic modeling is useful at discovering the development and change in latent topics over time. However, present methodology relies on algorithms that separate document and word representations. This prevents the creation of a meaningful embedding space where changes in word usage and documents can be directly analyzed in a temporal context. This paper proposes an expansion of the compass-aligned temporal Word2Vec methodology into dynamic topic modeling. Such a method allows for the direct comparison of word and document embeddings across time in dynamic topics. This enables the creation of visualizations that incorporate temporal word embeddings within the context of documents into topic visualizations. In experiments against the current state-of-the-art, our proposed method demonstrates overall competitive performance in topic relevancy and diversity across temporal datasets of varying size. Simultaneously, it provides insightful visualizations focused on temporal word embeddings while maintaining the insights provided by global topic evolution, advancing our understanding of how topics evolve over time.

Cluster analysis↗

Efficient computation of N -point correlation functions in D dimensions

We present efficient algorithms for computing the N-point correlation functions (NPCFs) of random fields in arbitrary D-dimensional homogeneous and isotropic spaces. Such statistics appear throughout the physical sciences and provide a natural tool to describe stochastic processes. Typically, algorithms for computing the NPCF components have $\mathscr O$(n N ) complexity (for a dataset containing n particles); their application is thus computationally infeasible unless N is small. By projecting the statistic onto a suitably defined angular basis, we show that the estimators can be written in a separable form, with complexity $\mathscr O$(n 2 ) or $\mathscr O$(n g log n g ) if evaluated using a Fast Fourier Transform on a grid of size n g . Our decomposition is built upon the D-dimensional hyperspherical harmonics; these form a complete basis on the (D – 1) sphere and are intrinsically related to angular momentum operators. Concatenation of (N – 1) such harmonics gives states of definite combined angular momentum, forming a natural separable basis for the NPCF. As N and D grow, the number of basis components quickly becomes large, providing a practical limitation to this (and all other) approaches: However, the dimensionality is greatly reduced in the presence of symmetries; for example, isotropic correlation functions require only states of zero combined angular momentum. We provide a Julia package implementing our estimators and show how they can be applied to a variety of scenarios within cosmology and fluid dynamics. The efficiency of such estimators will allow higher-order correlators to become a standard tool in the analysis of random fields.

97 MATHEMATICS AND COMPUTING↗

Development of a Ground Multi-Mission Low-Cost Optical Terminal (LCOT) for Free-Space Optical Communications

Once confined to the realm of laboratory experiments and theoretical papers, space-based laser communications (lasercomm) are on the verge of achieving mainstream status. Organizations from Facebook to NASA, and missions from cubesats to Orion are employing lasercomm to achieve gigabit communication speeds at mass and power requirements lower than that of traditional radio frequency (RF) methods. Since first demonstrating free-space optical communications services with Lunar Laser Communications Demonstration (LLCD) in 2013, NASA has invested in developing optical communications technologies and capabilities to enhancing its space communications networks. Along with evolving optical space terminals, NASA is also developing lasercomm grounds stations capable of meeting the rapidly increasing data volume demands of upcoming missions, from low-earth to lunar orbits and beyond and integrating these advanced capabilities into its Near-Space and Deep Space Networks. To meet this emerging need, the Low-Cost Optical Terminal (LCOT) project at NASA’s Goddard Space Flight Center (GSFC) is designing, building and validating a prototype for a flexible, multi-mission, and economical optical ground terminal that could be used as a blueprint for a global network of optical ground stations, capable of supporting a wide variety of missions. To date a major impediment to widespread adoption of laser communication has been the lack of an existing ground network infrastructure. A mission that wishes to take advantage of laser communication not only needs to invest in an optical space terminal, but it must also finance the creation of ground terminals to receive the downlink signal. This adds significant additional cost. Missions that do decide to incur the cost of financing a network of ground terminals end up building highly specialized optical receivers that are operable as receivers for that specific mission only. Significant Non-Recurring Engineering (NRE) cost is invested to build highly specialized one-of-a-kind ground terminals that go into storage after that particular mission is over. This is not an economical approach and does nothing to grow the number of optical ground stations available to future missions. In essence each mission that wants to take advantage of the benefits of lasercom has to start from scratch to provide a ground terminal network to support it. As long as this is the case, the cost for using laser communications will be too high for most missions to consider. LCOT intends to close this gap in technology by designing and developing a standard optical ground terminal design that is flexible enough to serve as a receiver for a wide range of future missions – a ground terminal that can be quickly reconfigured to receive downlinks at different wavelengths using different signal formats. Not only does LCOT have the industry-building objectives of utilizing commercial-off-the-shelf (COTS) components to the maximum extent possible, but also spurs the commercial development of other necessary lasercomm components not currently offered by industry. Finally, LCOT will give NASA scientists and engineers a facility where they can gain real-world experience with optical communications. It will give engineers a cost-effective way to try out new concepts and processes by providing the infrastructure for such testing. In this way it is hoped LCOT will serve as a stimulus for innovation in optical communications and speed its widespread adoption by future missions. The LCOT is comprised of five subsystems: Free-Space Optical, Transceiver, Amplifier, Monitor and Control, and Observatory Infrastructure. In August 2021, the LCOT team installed a 70 cm telescope, developed by Planewave Instruments that was optimized for optical communications. Free-Space Optical subsystem comprises of the telescope and its associated hardware, including a transmitter optical assembly, wide field cameras, two optical benches, and an adaptive optics subsystem. The transmit optical assembly, a unique concept design, is a cluster of four functionally independent transmit subassemblies located on the receive telescope. In addition to receiving optical signals and directing the expanded beam with high precision to the space terminal, it also performs tracking functions. The transmit optical assembly will support operations from Low Earth Orbit (LEO) through lunar and will be used as a template for industry manufacturing. The Optical Infrastructure subsystem is responsible for providing environmentally controlled shelters for LCOT equipment and various other systems. To maintain the safety and proper functionality of the telescope, a 16 ft Astrohaven clamshell dome procured which provides all-sky coverage without the need to rotate the dome. Additionally, Atmospheric Monitoring Assembly (AMA) will be part of optical infrastructure subsystem to ensure accurate performance of the LCOT. Like existing optical ground stations, LCOT will measure standard weather station parameters, infrared all sky image of cloud cover, and cloud height. LCOT, however, adds requirements for measuring night time seeing and, in the future, daytime seeing. Unlike other optical ground terminals, the LCOT is transceiver agnostic; user transceivers may be duplex transceivers, standalone receivers, or standalone transmitters with or without acquisition beacon functionality. As such, the LCOT project accommodates testing with external customer transceivers in a flexible manor, further complimenting its intended multi-mission goals. Another unique component of LCOT is the use of a new amplifier technology – the Very Large Mode Area (VLMA) amplifiers. This new technology allows LCOT to avoid the issues faced by previous laser communications ground terminals, gives users more flexibility and modular capability, and is capable of reaching an order of magnitude higher peak power than traditional High Power Optical Amplifiers (HPOA). One drawback of the VLMA HPOA approach is that the amplified light is output into free-space. The solution developed by LCOT is an optics train that couples the output of the VLMA amplifier into a short fiber for transport to the transmit telescopes with high efficiency. Like many of the LCOT components, a set of detailed manufacturing drawings have been created for the optics train to allow any machine shop with a multi-axis Computer Numerical Control (CNC) machine to fabricate the piece parts from commonly available materials. In line with the goals of LCOT, the monitor and control functions are developed as a modular and flexible system with the ability to support future hardware or algorithm changes, minimizing disruptions. A main priority of development in the Monitor and Control Subsystem (MCS) is the safety monitor system.

laser communication↗

Development of a Ground Multi-mission Low Cost Optical Terminal(LCOT) for Free-Space Optical Communication

Once confined to the realm of laboratory experiments and theoretical papers, space-based laser communications (lasercomm) are on the verge of achieving mainstream status. Organizations from Facebook to NASA, and missions from cubesats to Orion are employing lasercomm to achieve gigabit communication speeds at mass and power requirements lower than that of traditional radio frequency (RF) methods. Since first demonstrating free-space optical communications services with Lunar Laser Communications Demonstration (LLCD) in 2013, NASA has invested in developing optical communications technologies and capabilities to enhancing its space communications networks. Along with evolving optical space terminals, NASA is also developing lasercomm grounds stations capable of meeting the rapidly increasing data volume demands of upcoming missions, from low-earth to lunar orbits and beyond and integrating these advanced capabilities into its Near-Space and Deep Space Networks. To meet this emerging need, the Low-Cost Optical Terminal (LCOT) project at NASA’s Goddard Space Flight Center (GSFC) is designing, building and validating a prototype for a flexible, multi-mission, and economical optical ground terminal that could be used as a blueprint for a global network of optical ground stations, capable of supporting a wide variety of missions. To date a major impediment to widespread adoption of laser communication has been the lack of an existing ground network infrastructure. A mission that wishes to take advantage of laser communication not only needs to invest in an optical space terminal, but it must also finance the creation of ground terminals to receive the downlink signal. This adds significant additional cost. Missions that do decide to incur the cost of financing a network of ground terminals end up building highly specialized optical receivers that are operable as receivers for that specific mission only. Significant Non-Recurring Engineering (NRE) cost is invested to build highly specialized one-of-a-kind ground terminals that go into storage after that particular mission is over. This is not an economical approach and does nothing to grow the number of optical ground stations available to future missions. In essence each mission that wants to take advantage of the benefits of lasercom has to start from scratch to provide a ground terminal network to support it. As long as this is the case, the cost for using laser communications will be too high for most missions to consider. LCOT intends to close this gap in technology by designing and developing a standard optical ground terminal design that is flexible enough to serve as a receiver for a wide range of future missions – a ground terminal that can be quickly reconfigured to receive downlinks at different wavelengths using different signal formats. Not only does LCOT have the industry-building objectives of utilizing commercial-off-the-shelf (COTS) components to the maximum extent possible, but also spurs the commercial development of other necessary lasercomm components not currently offered by industry. Finally, LCOT will give NASA scientists and engineers a facility where they can gain real-world experience with optical communications. It will give engineers a cost-effective way to try out new concepts and processes by providing the infrastructure for such testing. In this way it is hoped LCOT will serve as a stimulus for innovation in optical communications and speed its widespread adoption by future missions. The LCOT is comprised of five subsystems: Free-Space Optical, Transceiver, Amplifier, Monitor and Control, and Observatory Infrastructure. In August 2021, the LCOT team installed a 70 cm telescope, developed by Planewave Instruments that was optimized for optical communications. Free-Space Optical subsystem comprises of the telescope and its associated hardware, including a transmitter optical assembly, wide field cameras, two optical benches, and an adaptive optics subsystem. The transmit optical assembly, a unique concept design, is a cluster of four functionally independent transmit subassemblies located on the receive telescope. In addition to receiving optical signals and directing the expanded beam with high precision to the space terminal, it also performs tracking functions. The transmit optical assembly will support operations from Low Earth Orbit (LEO) through lunar and will be used as a template for industry manufacturing. The Optical Infrastructure subsystem is responsible for providing environmentally controlled shelters for LCOT equipment and various other systems. To maintain the safety and proper functionality of the telescope, a 16 ft Astrohaven clamshell dome procured which provides all-sky coverage without the need to rotate the dome. Additionally, Atmospheric Monitoring Assembly (AMA) will be part of optical infrastructure subsystem to ensure accurate performance of the LCOT. Like existing optical ground stations, LCOT will measure standard weather station parameters, infrared all sky image of cloud cover, and cloud height. LCOT, however, adds requirements for measuring night time seeing and, in the future, daytime seeing. Unlike other optical ground terminals, the LCOT is transceiver agnostic; user transceivers may be duplex transceivers, standalone receivers, or standalone transmitters with or without acquisition beacon functionality. As such, the LCOT project accommodates testing with external customer transceivers in a flexible manor, further complimenting its intended multi-mission goals. Another unique component of LCOT is the use of a new amplifier technology – the Very Large Mode Area (VLMA) amplifiers. This new technology allows LCOT to avoid the issues faced by previous laser communications ground terminals, gives users more flexibility and modular capability, and is capable of reaching an order of magnitude higher peak power than traditional High Power Optical Amplifiers (HPOA). One drawback of the VLMA HPOA approach is that the amplified light is output into free-space. The solution developed by LCOT is an optics train that couples the output of the VLMA amplifier into a short fiber for transport to the transmit telescopes with high efficiency. Like many of the LCOT components, a set of detailed manufacturing drawings have been created for the optics train to allow any machine shop with a multi-axis Computer Numerical Control (CNC) machine to fabricate the piece parts from commonly available materials. In line with the goals of LCOT, the monitor and control functions are developed as a modular and flexible system with the ability to support future hardware or algorithm changes, minimizing disruptions. A main priority of development in the Monitor and Control Subsystem (MCS) is the safety monitor system.

Haleh Safavi↗

A new calibration method for charm jet identification validated with proton-proton collision events at $\sqrt{s}$ =13 TeV

Many measurements at the LHC require efficient identification of heavy-flavour jets, i.e. jets originating from bottom (b) or charm (c) quarks. An overview of the algorithms used to identify c jets is described and a novel method to calibrate them is presented. This new method adjusts the entire distributions of the outputs obtained when the algorithms are applied to jets of different flavours. It is based on an iterative approach exploiting three distinct control regions that are enriched with either b jets, c jets, or light-flavour and gluon jets. Results are presented in the form of correction factors evaluated using proton-proton collision data with an integrated luminosity of 41.5 fb$^{-1}$ at √s = 13 TeV, collected by the CMS experiment in 2017. The closure of the method is tested by applying the measured correction factors on simulated data sets and checking the agreement between the adjusted simulation and collision data. Furthermore, a validation is performed by testing the method on pseudodata, which emulate various mismodelling conditions. The calibrated results enable the use of the full distributions of heavy-flavour identification algorithm outputs, e.g. as inputs to machine-learning models. Thus, they are expected to increase the sensitivity of future physics analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SMC 2021 : Analyzing Resource Utilization and User Behavior on Titan Supercomputer

Resource utilization statistics of submitted jobs on a supercomputer can help us understand how users from various scientific domains use HPC platforms and better design a job scheduler. We explore to generate insight regarding workload distribution and usage pattern domains from job scheduler trace, GPU failure information, and project-specific information collected from Titan supercomputer. Furthermore, we want to know how the scheduler performance varies over time and how the users' scheduling behavior changes following a system failure. These observations have the potential to provide valuable insight, which is helpful to prepare for system failures. These practices will help us develop and apply novel machine learning algorithms in understanding system behavior, requirement, and better scheduling of HPC systems. There are two datasets, RUR and GPU. RUR: This dataset is the job scheduler traces collected from the Titan supercomputerfrom 01/01/2015 to 07/31/2019 (2015.csv - 2019.csv). These were collected usingResource Utilization Report (RUR), a Cray-developed resource-usage data collectionand reporting system. It contains the usage information of its critical resources (CPU,Memory, GPU, and I/O) of each running job on Titan during that period [2]. ProjectAreas: Every job is associated with a project ID. TheProjectAreas.csvdatasetprovides a mapping of the project ID to its domain science. GPU: There have been some hardware-related issues in the GPUs in Titan that caused some GPUs to fail, sometimes irrecoverably during some job runs. This dataset provides information regarding these failures during the execution of the submitted jobs. GPUs on Titan are uniquely identified by a serial number (SN), and they are installed in a location. A GPU can be installed in a location, then removed from that location following a failure, and then re-installed in a different location after fixing the problem. If the failure can't be recovered, the GPU might be removed entirely from Titan. There are two prominent types of failures that resulted in the removal of GPUs from Titan: Double Bit Error (DBE) and Out of the Bus (OTB). The dataset (gc_full.csv) has the following fields: 1. SN : Serial number of a GPU 2. location : The location where it is installed 3. insert : The time when it was inserted into that location 4. remove : The time when it was removed from that location 5. duration : Amount of time the GPU spent in this location 6. out : If the device was taken out entirely w/o a re-installment into a new location. 7. event : If the GPU was taken out entirely, the reason for its removal. To learn more about this dataset, please refer to the git repositoryhttps://github.com/olcf/TitanGPULifeand the related publication [1]. References [1] George Ostrouchov, Don Maxwell, Rizwan A Ashraf, Christian Engelmann, MallikarjunShankar, and James H Rogers. Gpu lifetimes on titan supercomputer: Survival analysisand reliability. InSC20: International Conference for High Performance Computing,Networking, Storage and Analysis, pages 1-14. IEEE, 2020. [2] Feiyi Wang, Sarp Oral, Satyabrata Sen, and Neena Imam. Learning from five-yearresource-utilization data of titan system. In2019 IEEE International Conference onCluster Computing (CLUSTER), pages 1-6. IEEE, 2019.

42 ENGINEERING↗

Classical Preoptimization Approach for ADAPT-VQE: Maximizing the Potential of High-Performance Computing Resources to Improve Quantum Simulation of Chemical Applications

The ADAPT-VQE algorithm is a promising method for generating a compact ansatz based on derivatives of the underlying cost function, and it yields accurate predictions of electronic energies for molecules. In this work, we report the implementation and performance of ADAPT-VQE with our recently developed sparse wave function circuit solver (SWCS) in terms of accuracy and efficiency for molecular systems with up to 52 spin orbitals. The SWCS can be tuned to balance computational cost and accuracy, which extends the application of ADAPT-VQE for molecular electronic structure calculations to larger basis sets and a larger number of qubits. Using this tunable feature of the SWCS, we propose an alternative optimization procedure for ADAPT-VQE to reduce the computational cost of the optimization. Furthermore, by preoptimizing a quantum simulation with a parametrized ansatz generated with ADAPT-VQE/SWCS, we aim to utilize the power of classical high-performance computing in order to minimize the work required on noisy intermediate-scale quantum hardware, which offers a promising path toward demonstrating quantum advantage for chemical applications.

ADAPT-VQE↗

Using graph neural networks to reconstruct charged pion showers in the CMS High Granularity Calorimeter

A novel method to reconstruct the energy of hadronic showersin the CMS High Granularity Calorimeter (HGCAL) is presented. TheHGCAL is a sampling calorimeter with very fine transverse andlongitudinal granularity. The active media are silicon sensors andscintillator tiles readout by SiPMs and the absorbers are acombination of lead and Cu/CuW in the electromagnetic section, andsteel in the hadronic section. The shower reconstruction method isbased on graph neural networks and it makes use of a dynamicreduction network architecture. It is shown that the algorithm isable to capture and mitigate the main effects that normally hinderthe reconstruction of hadronic showers using classicalreconstruction methods, by compensating for fluctuations in themultiplicity, energy, and spatial distributions of the shower'sconstituents. The performance of the algorithm is evaluated usingtest beam data collected in 2018 prototype of the CMS HGCALaccompanied by a section of the CALICE AHCAL prototype. Thecapability of the method to mitigate the impact of energy leakagefrom the calorimeter is also demonstrated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Semi-analytical covariance matrices for two-point correlation function for DESI 2024 data

We present an optimized way of producing the fast semi-analytical covariance matrices for the Legendre moments of the two-point correlation function, taking into account survey geometry and mimicking the non-Gaussian effects. We validate the approach on simulated (mock) catalogs for different galaxy types, representative of the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, used in 2024 analyses. We find only a few percent differences between the mock sample covariance matrix and our results, which can be expected given the approximate nature of the mocks, although we do identify discrepancies between the shot-noise properties of the DESI fiber assignment algorithm and the faster approximation (emulator) used in the mocks. Importantly, we find a close agreement (≤ 8% relative differences) in the projected errorbars for distance scale parameters for the baryon acoustic oscillation measurements. This confirms our method as an attractive alternative to simulation-based covariance matrices, especially for non-standard models or galaxy sample selections, making it particularly relevant to the broad current and future analyses of DESI data.

79 ASTRONOMY AND ASTROPHYSICS↗