Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cluster Computing”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

3-D Simulations of earthquakes rupture jumps: 1. Homogeneous pre-stress conditions

SUMMARY Observational and modelling studies indicate that earthquake ruptures can jump between fault sections as large as ∼3 and ∼5 km for compressional and extensional offsets, respectively. Here, we compare characteristics of the rupture jump process on parallel but offset fault sections from traditional 3-D dynamic rupture simulations governed by slip weakening friction using the finite element code, FaultMod, to those from quasi-dynamic simulations governed by rate- and state-dependent friction (rate-state friction) using the code RSQSim. These simulations use spatially uniform initial stresses. For a variety of measures the rupture renucleation position on the offset fault, the rate-state friction and slip weakening friction models produce very similar results. The principal difference is the additional occurrence of delayed rupture jumps that arise from the time- and stress-dependent nucleation that is characteristic of rate-state friction. For immediate rupture jumps, models with slip weakening friction span greater offsets than those with rate-state friction. However the jump distances are nearly identical when delayed rupture jumps are included in the comparisons. We propose that delayed rupture jumps are the likely mechanism for adjacent large-earthquake pairs and clusters. Based on the similarity of renucleation positions with both dynamic and quasi-dynamic models, we conclude that the renucleation positions for rupture initiation on the receiver fault (separated by less than ∼3 km from the source fault) are primarily controlled by static stress changes induced by slip on the initiating fault. However, in light of the slightly greater maximum jump distances (>3 km) seen with the dynamic slip weakening friction model, dynamic stress changes from seismic waves play an increasingly important role as offset distances increase.

, RSQSim↗
DOE OSTI2019

DCA++: A software framework to solve correlated electron problems with modern quantum cluster methods

We present the first open release of the DCA++ project, a high-performance research software framework to solve quantum many-body problems with cutting edge quantum cluster algorithms. DCA++ implements the dynamical cluster approximation (DCA) and its DCA extension with a continuous self-energy. The algorithms capture nonlocal correlations in strongly correlated electron systems, thereby giving insight into high-$T_c$ superconductivity. The code’s scalability allows efficient usage of systems at all scales, from workstations to leadership computers. With regard to the increasing heterogeneity of modern computing machines, DCA++ provides portable performance on conventional and emerging new architectures, such as hybrid CPU–GPU, sustaining multiple petaflops on ORNL’s Titan and CSCS’ Piz Daint supercomputers. Moreover, we show how sustainable and scalable development of the code base has been achieved by adopting standard techniques of the software industry. These include employing a distributed version control system, applying test-driven development and following continuous integration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Angular-momentum projection in coupled-cluster theory: Structure of 34 Mg

Single- reference coupled-cluster theory is an accurate and affordable computational method for the nuclear many-body problem. For open-shell nuclei, the reference state typically breaks rotational invariance and angular momentum must be restored as a good quantum number. We perform angular-momentum projection after variation and employ the disentangled coupled-cluster formalism and a Hermitian approach. We compare our results with benchmarks for 8 Be and 20 Ne using a two-nucleon interaction from chiral effective field theory and for pf-shell nuclei within the traditional shell model. We compute the rotational band in the exotic nucleus 34 Mg and find agreement with data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

DUST IN CLUSTERS: SEPARATING THE CONTRIBUTION OF GALAXIES AND INTRACLUSTER MEDIA

We have analyzed a sample of 327 clusters of galaxies, spanning the range of 0.06–0.70 in redshift. Strong constraints on their mean intracluster emission of dust have been obtained using maps and catalogs from the Herschel MerMES project; within a radius of 5 arcmin centered in each cluster, the 95% C.L. limits obtained are 86.6, 48.2, and 30.9 mJy at the observed frequencies of 250, 350, and 500 μ m. From these restrictions, and assuming physical parameters typical of interstellar media in the Milky Way, we have obtained tight upper limits on the visual extinction of background galaxies due to the intracluster media (ICM): A {sub V} (95% C.L.) ≲ 10{sup −3} mag. Strong constraints are also obtained for the mass of such dust; for instance, using the data at 350 μ m we establish a 95% upper limit of <10{sup 9} M {sub ⊙} within a circle with a radius of 5 arcmin centered in the clusters. This corresponds to a fraction of the total mass of the clusters of 9.5 × 10{sup −6}, and indicates a deficiency in the gas-to-dust ratio in the ICM by about three orders of magnitude in relation to the value found in the Milky Way. Computing the total infrared luminosity of the clusters in three ranges of redshift (0.05–0.24, 0.24–0.42, and 0.42–0.71) and two ranges of mass (<10{sup 14} and >10{sup 14} M {sub ⊙}), respectively, a strong evolution of luminosity in redshift ( L ∼ z {sup 1.5}) for both ranges of masses is found. The results indicate a strong declining in star formation rate with time in the last ∼6 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Coupled-cluster representation of Green function employing modified spectral resolutions of similarity transformed Hamiltonians

In this paper we discuss a new formalism for producing an analytic coupled-cluster (CC) Green's function for an N-electron system by shifting the poles of similarity transformed Hamiltonians represented in N − 1 and N + 1 electron Hilbert spaces. Simple criteria are derived for the states in N − 1 and N + 1 electron spaces that are then corrected in the spectral resolution of the corresponding matrix representations of the similarity transformed Hamiltonian. The accurate description of excited state processes within a Green's function formalism would be of significant importance to a number of scientific communities ranging from physics and chemistry to engineering and the biological sciences. This is because the Green's function methodology provides a direct path for not only calculating properties whose underlying origins come from coupled many-body interactions but also provides a straightforward path for calculating electron transport, response, and correlation functions that allows for a direct link with experiment. As a special case of this general formulation, we discuss the application of this technique for Green's function defined by the CC with singles and doubles representation of the ground-state wave function.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Artificial Neural Network Representation of Dynamic Criticality Excursion Experiment Data

This paper outlines time dependent considerations in applying the artificial neural network (ANN) paradigm to dynamic nuclear criticality excursion experiment data. The ANN amongst a variety of uses, can apply to function fitting of multiple input-multiple output data sets, pattern recognition, cluster analysis, and time series prediction. The ANN is a component of Machine Learning (ML) data analysis. The practical application of ML to criticality safety has been funded through the Y-12 Plant Directed Research and Demonstration (PDRD) program to support criticality safety and emergency planning efforts. The initial effort considered function fitting of criticality excursion data. Recent efforts successfully extended the ANN to criticality safety code validation. The initial ANN work applied to criticality excursions centered around multivariate function fitting of static (tabular) data taken from the French Consequences Radiologiques d'un Accident de Criticite (CRAC) and SILENE reactor campaigns. The former report (LeCorche) summarized important tabular data for the first 44 CRAC experiments. The initial ANN was trained to specific time-independent inputs, then applied to predict peak fission rate and total fissions for a wide variety of fissile solutions process conditions. This work extends that effort to apply the ANN paradigm to dynamic prediction of time-dependent criticality excursion experiment data. A dynamic criticality excursion in fissile solution represents a complex multiphysics phenomenon. Factors include among other things, coupled neutronic, multiphase fluid flow, temperature, pressure, acoustic shock, boiling heat transfer, reactivity feedback, and radiolytic gas considerations. There have been many efforts to model (over time) as coupled ordinary and (space) as partial differential equations (ODE, PDE). The ANN representation noted here does not rely on solving explicit ODE or PDEs. Rather the intent is to demonstrate that as a time series, the criticality excursion can be predicted by training specific ANN architectures to experimental input and output time series exclusively. The rationale for extending this work to the time-dependent regime is important in developing a 'criticality excursion emulator' solely on the basis of learning characteristics from experiment dynamic data. An emulator is different than a computer code simulator driven by the solution of ODEs and PDEs given select 'tuning parameters' that match outputs to experiments. A criticality excursion emulator, derived exclusively from an ANN or ML, could predict most likely excursion characteristics, trained from historic dynamic experiment data, presented with input values within the range of those experiments. The work noted shows that generic time-series ANN architectures can be satisfactorily applied to the prediction of dynamic, extremely chaotic criticality excursion power given an initial time dependent input series. In reality, an excursion emulator could be driven by user-supplied initial (t=0), time independent or constant values to generate time dependent outputs. Thus ongoing work includes developing a custom ANN such that initial dynamic conditions would be predicted. Those results then serve as inputs to a customized ANN architecture currently in development. (author)

42 ENGINEERING↗
DOE OSTI2015

A two-stage method for microcalcification cluster segmentation in mammography by deformable models

Purpose: Segmentation of microcalcification (MC) clusters in x-ray mammography is a difficult task for radiologists. Accurate segmentation is prerequisite for quantitative image analysis of MC clusters and subsequent feature extraction and classification in computer-aided diagnosis schemes. Methods: In this study, a two-stage semiautomated segmentation method of MC clusters is investigated. The first stage is targeted to accurate and time efficient segmentation of the majority of the particles of a MC cluster, by means of a level set method. The second stage is targeted to shape refinement of selected individual MCs, by means of an active contour model. Both methods are applied in the framework of a rich scale-space representation, provided by the wavelet transform at integer scales. Segmentation reliability of the proposed method in terms of inter and intraobserver agreements was evaluated in a case sample of 80 MC clusters originating from the digital database for screening mammography, corresponding to 4 morphology types (punctate: 22, fine linear branching: 16, pleomorphic: 18, and amorphous: 24) of MC clusters, assessing radiologists’ segmentations quantitatively by two distance metrics (Hausdorff distance—HDIST{sub cluster}, average of minimum distance—AMINDIST{sub cluster}) and the area overlap measure (AOM{sub cluster}). The effect of the proposed segmentation method on MC cluster characterization accuracy was evaluated in a case sample of 162 pleomorphic MC clusters (72 malignant and 90 benign). Ten MC cluster features, targeted to capture morphologic properties of individual MCs in a cluster (area, major length, perimeter, compactness, and spread), were extracted and a correlation-based feature selection method yielded a feature subset to feed in a support vector machine classifier. Classification performance of the MC cluster features was estimated by means of the area under receiver operating characteristic curve (Az ± Standard Error) utilizing tenfold cross-validation methodology. A previously developed B-spline active rays segmentation method was also considered for comparison purposes. Results: Interobserver and intraobserver segmentation agreements (median and [25%, 75%] quartile range) were substantial with respect to the distance metrics HDIST{sub cluster} (2.3 [1.8, 2.9] and 2.5 [2.1, 3.2] pixels) and AMINDIST{sub cluster} (0.8 [0.6, 1.0] and 1.0 [0.8, 1.2] pixels), while moderate with respect to AOM{sub cluster} (0.64 [0.55, 0.71] and 0.59 [0.52, 0.66]). The proposed segmentation method outperformed (0.80 ± 0.04) statistically significantly (Mann-Whitney U-test, p < 0.05) the B-spline active rays segmentation method (0.69 ± 0.04), suggesting the significance of the proposed semiautomated method. Conclusions: Results indicate a reliable semiautomated segmentation method for MC clusters offered by deformable models, which could be utilized in MC cluster quantitative image analysis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Pair extended coupled cluster doubles

The accurate and efficient description of strongly correlated systems remains an important challenge for computational methods. Doubly occupied configuration interaction (DOCI), in which all electrons are paired and no correlations which break these pairs are permitted, can in many cases provide an accurate account of strong correlations, albeit at combinatorial computational cost. Recently, there has been significant interest in a method we refer to as pair coupled cluster doubles (pCCD), a variant of coupled cluster doubles in which the electrons are paired. This is simply because pCCD provides energies nearly identical to those of DOCI, but at mean-field computational cost (disregarding the cost of the two-electron integral transformation). Here, we introduce the more complete pair extended coupled cluster doubles (pECCD) approach which, like pCCD, has mean-field cost and reproduces DOCI energetically. We show that unlike pCCD, pECCD also reproduces the DOCI wave function with high accuracy. Moreover, pECCD yields sensible albeit inexact results even for attractive interactions where pCCD breaks down.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Tuning Catalytic Reactivity via Wetting Control through Oxygen Vacancies: Ru Clusters on Anatase TiO 2 and CeO 2 Supports

The shape of supported metal particles regulates their catalytic reactivity and is determined by the degree of wetting between the metal particle and the support surface. Flattened particles that wet support surfaces were reported in various catalytic systems, particularly in the subnanometer size regime. Such consequential metal–support wetting phenomena are poorly understood, and methods to study them on powder catalysts under realistic conditions are lacking. Here, we investigate the size-dependent wetting behaviors of Ru particles on two reducible-oxide supports, anatase TiO 2 (TiO 2 -A) and CeO 2 , under reducing catalytic conditions. X-ray absorption spectroscopy (XAS), low-energy ion scattering (LEIS), and density functional theory (DFT) are combined to determine the shape of Ru particles. Ru particles remain three-dimensional without wetting the TiO 2 -A support within the coverage range studied (0.06–0.98 Ru nm –2 ). In contrast, at low coverages (<0.25 Ru nm –2 ), Ru wets the CeO 2 support to form flat, disordered structures. The higher wettability of CeO 2 than TiO 2 -A is attributed to oxygen vacancies in the near-surface region. The shape difference between small Ru particles or clusters on the two supports leads to drastically contrasting catalytic reactivities in polyolefin hydrogenolysis, despite similar diameters. As a result, this work highlights the implications of metal–support wetting, or cluster shape, on catalytic behaviors of small metal clusters, while establishing the foundation for future systematic studies of such a phenomenon in realistic systems, by delivering a multitechnique methodology and revealing governing fundamental principles.

42 ENGINEERING↗
DOE OSTI2015

Status on the Development of a Modeling and Simulation Framework for the Economic Assessment of Nuclear Hybrid Energy Systems

An effort to design and build a modeling and simulation framework to assess the economic viability of Nuclear Hybrid Energy Systems (NHES) was undertaken in fiscal year 2015 (FY15). The purpose of this report is to document the various tasks associated with the development of such a framework and to provide a status on its progress. Several tasks have been accomplished. First, starting from a simulation strategy, a rigorous mathematical formulation has been achieved in which the economic optimization of a Nuclear Hybrid Energy System is presented as a constrained robust (under uncertainty) optimization problem. Some possible algorithms for the solution of the optimization problem are presented. A variation of the Simultaneous Perturbation Stochastic Approximation algorithm has been implemented in RAVEN and preliminary tests have been performed. The development of the software infrastructure to support the simulation of the whole NHES has also moved forward. The coupling between RAVEN and an implementation of the Modelica language (OpenModelica) has been implemented, migrated under several operating systems and tested using an adapted model of a desalination plant. In particular, this exercise was focused on testing the coupling of the different code systems; testing parallel, computationally expensive simulations on the INL cluster; and providing a proof of concept for the possibility of using surrogate models to represent the different NHES subsystems. Another important step was the porting of the RAVEN code under the Windows™ operating system. This accomplishment makes RAVEN compatible with the development environment that is being used for dynamic simulation of NHES components. A very simplified model of a NHES on the electric market has been built in RAVEN to confirm expectations on the analysis capability of RAVEN to provide insight into system economics and to test the capability of RAVEN to identify limit surfaces even for stochastic constraints. This capability will be needed in the future to enforce the stochastic constraints on the electric demand coverage from the NHES. The development team gained experience with many of the tools that are currently envisioned for use in the economic analysis of NHES and completed several important steps. Given the complexity of the project, preference has been given to a structural approach in which several independent efforts have been used to build the cornerstone of the simulation framework. While this is good approach in establishing such a complex framework, it may delay reaching more complete results on the performance of analyzed system configurations. The integration of the previously reported exergy analysis approach was initially proposed as part of this milestone. However, in reality, the exergy-based apportioning of cost will take place only in a second stage of the implementation since it will be used to properly allocate cost among the different NHES subsystems. Therefore, exergy does not appear at the level of the main drivers in the analysis framework; the latter development of the base framework is the focus of this report.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2018

The Physics and Chemistry of Cluster-based Catalyst Systems

In this project we used first principles calculations to investigate the physical and chemical properties of atomic clusters, small clumps of matter containing up to a few hundred atoms. The calculations are based on density functional theory, the most efficient and accurate method available for computing and thereby predicting the properties of molecules, clusters, and solids. An initial focus of our work involved studying the reactivity of small Pd clusters. We demonstrated that H2 molecules readily dissociate on Pd clusters and that multiple dissociations can occur on clusters larger than a few Pd atoms. For example, up to eight H2 molecules will dissociate on Pd13 at room temperature. A second focus of the project was to study the evolution of the electrical behavior of clusters, in particular, the evolution of metallic behavior in clusters of metal atoms. We developed a technique for using the electric dipole polarizability of the clusters – the response of the cluster to an applied electric field -- as a probe of this behavior. By decomposing the polarizability into so-called charge transfer (or metallic) and local (or insulating) contributions, we obtained the remarkable result that the per atom charge transfer contributions in sodium clusters attain the same value as occurs in bulk sodium metal for clusters containing as few as 10 atoms. Our calculated total polarizabilities matched experimentally measured trends essentially perfectly, validating the reliability of the calculations and supporting our conclusion that tiny clusters already display bulk-like metallic character. Similar results were obtained for other metal atom clusters, including potassium, aluminum, and zinc. We discovered that the evolution toward bulk metallic behavior for each of these cluster species is essentially identical, when judged using these polarizability-based descriptors. This suggests a universality in the approach of metal clusters to bulk-like metallicity. This discovery is a major outcome of our project. We also used our polarizability-based metrics to investigate the evolution of the properties of Si and NaCl clusters. These systems are non-metallic in the bulk limit. For Si, we find the behavior of small clusters to be metal-like and we attribute this to bonding arrangements in the clusters that are very different from the open structure of crystalline silicon. As the clusters get larger, the bonding in the cluster interior becomes more and more like that of the bulk material. We expect a transition from metal-like to insulator-like behavior to accompany this change is structure. Our results point to a transition at sizes beyond 146 atoms, the largest size treated in our study. We hope to study this transition further in the future. A third thrust of our project involved work on a new approach to eliminating unphysical electron self-interaction from density functional theory calculations. Two electrons exert repulsive forces on each other, but a single electron should not interact with itself. Because of approximations made in density functional theory, individual electrons effectively interact with themselves. The effect of this self-interaction is often small, especially when atoms are near their equilibrium positions in molecules. But the effects can be pronounced when bonds are stretched and also in some situations involving transition metals even when the atoms are near equilibrium. Our research involved the further development of a method for removing self-interaction called the Fermi-Löwdin orbital self-interaction correction or FLO-SIC. We conducted the first systematic FLO-SIC calculations for the individual atoms ranging from Li – Kr. We also carried out a critical comparison of the results of FLO-SIC with an older method known as the Perdew-Zunger self-interaction correction (PZ-SIC). FLO-SIC can be thought of as a modified approach to PZ-SIC that is expected to be computationally more efficient. We demonstrated that the two methods are distinct, but give essentially the same results for quantities such as molecular binding energies or electron removal energies. We are continuing to develop the FLO-SIC methodology with a goal of creating a computationally efficient method that is free of self-interaction effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Accelerating Simulation for High-Fidelity PV Inverter System Reliability Assessment with High-Performance Computing

The overall cost of photovoltaic (PV) systems has shown a downward trend during the last decade; however, PV inverter failures account for the highest cost of operation and maintenance. To address this, reliability tools with powerful computation and better accuracy are required for the lifetime prediction and degradation evaluation of PV inverters. This paper proposes an event-driven parallel computing-based simulator. The proposed simulator applies high-performance computing techniques and other accessory optimization techniques-including cluster merging, adaptive model updates, and steady-state identification-to make reliability assessments for PV inverters under given input mission profiles and operating conditions with high efficiency and high fidelity. The main idea of the simulator and its workflow are introduced. Then, a demo PV inverter system simulator is implemented, and the speedup of the total simulations of the switching model reaches 123.03 times.

high-performance computing↗
DOE OSTI2023

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Assessment of low-scaling approximations to the equation of motion coupled-cluster singles and doubles equations

Methods for fast and reliable computation of electronic excitation energies are in short supply, and little is known about their systematic performance. This work reports a comparison of several low-scaling approximations to the equation of motion coupled cluster singles and doubles (EOM–CCSD) and linear-response coupled cluster singles and doubles (LR–CCSD) equations with other single reference methods for computing the vertical electronic transition energies of 11 small organic molecules. The methods, including second order equation-of-motion many-body perturbation theory (EOM–MBPT2) and its partitioned variant, are compared to several valence and Rydberg singlet states. We find that the EOM–MBPT2 method was rarely more than a tenth of an eV from EOM–CCSD calculated energies, yet demonstrates a performance gain of nearly 30%. The partitioned equation-of-motion approach, P–EOM–MBPT2, which is an order of magnitude faster than EOM–CCSD, outperforms the CIS(D) and CC2 in the description of Rydberg states. CC2, on the other hand, excels at describing valence states where P–EOM–MBPT2 does not. The difference between the CC2 and P–EOM–MBPT2 can ultimately be traced back to how each method approximates EOM–CCSD and LR–CCSD. The results suggest that CC2 and P–EOM–MBPT2 are complementary: CC2 is best suited for the description of valence states while P–EOM–MBPT2 proves to be a superior O(N{sup 5}) method for the description of Rydberg states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Structural and Dynamical Design Principles for Iron-Sulfur Clusters in Sequence Defined Polymers

This integrated synthetic and computational effort uses sequence defined polymers such as peptoids, and triazines, capable of mimicking a flexible protein environment to control the positioning of side chains around redox cofactors. This strategy will allow modulation of redox properties through both direct and indirect mechanisms. To this end the aim is defined to test the underlying hypothesis: that principles and critical functionalities governing the binding and tuning of naturally occurring electron carrier co-factors can be identified and used to predictably design and synthesize maquettes based on non-natural polymers resulting in new modular biomimetic scaffolds for Electron Transfer Materials. This will be achieved through design and synthesis of ligands that can stabilize iron-sulfur clusters and tune the redox potential based on sequence.

36 MATERIALS SCIENCE↗
DOE OSTI2015

CHANDRA and XMM-NEWTON observations of the bimodal PLANCK SZ -detected clustered CKG345.40-39.34 (A3716) with high and low entropy subcluster cores

Here, we present results from Chandra, XMM-Newton, and ROSAT observations of the Planck SZ-detected cluster A3716 (PLCKG345.40-39.34-G345). We show that G345 is, in fact, two subclusters separated on the sky by 400 kpc. We measure the subclusters' gas temperatures (~2–3 keV), total (~1–2 × 10 14 ${{M}_{\odot }}$) and gas (~1–2 × 10 13 ${{M}_{\odot }}$) masses, gas mass fraction within r500, entropy profiles, and X-ray luminosities (~10 43 erg s -1 ). Using the gas density and temperature profiles for both subclusters, we show that there is good (0.8σ) agreement between the expected Sunyaev–Zel'dovich signal predicted from the X-ray data and that measured from the Planck mission, and better agreement within 0.6σ when we re-computed the Planck value assuming a two component cluster model, with relative amplitudes fixed based on the X-ray data. Dynamical analysis shows that the two galaxy subclusters are very likely ($\gt 97\%$ probability) gravitationally bound, and in the most likely scenario, the subclusters will undergo core passage in 500 ± 200 Myr. The northern subcluster is centrally peaked and has a low entropy core, while the southern subcluster has a high central entropy. Finally, the high central entropy in the southern subcluster can be explained either by the mergers of several groups, as suggested by the presence of five giant ellipticals or by active galactic nucleus energy injection, as suggested by the presence of a strong radio source in one of its massive elliptical galaxies, or by a combination of both processes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

MVP: a modular viromics pipeline to identify, filter, cluster, annotate, and bin viruses from metagenomes

While numerous computational frameworks and workflows are available for recovering prokaryote and eukaryote genomes from metagenome data, only a limited number of pipelines are designed specifically for viromics analysis. With many viromics tools developed in the last few years alone, it can be challenging for scientists with limited bioinformatics experience to easily recover, evaluate quality, annotate genes, dereplicate, assign taxonomy, and calculate relative abundance and coverage of viral genomes using state-of-the-art methods and standards. Here, we describe Modular Viromics Pipeline (MVP) v.1.0, a user-friendly pipeline written in Python and providing a simple framework to perform standard viromics analyses. MVP combines multiple tools to enable viral genome identification, characterization of genome quality, filtering, clustering, taxonomic and functional annotation, genome binning, and comprehensive summaries of results that can be used for downstream ecological analyses. Overall, MVP provides a standardized and reproducible pipeline for both extensive and robust characterization of viruses from large-scale sequencing data including metagenomes, metatranscriptomes, viromes, and isolate genomes. As a typical use case, we show how the entire MVP pipeline can be applied to a set of 20 metagenomes from wetland sediments using only 10 modules executed via command lines, leading to the identification of 11,656 viral contigs and 8,145 viral operational taxonomic units (vOTUs) displaying a clear beta-diversity pattern. Further, acting as a dynamic wrapper, MVP is designed to continuously incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field of viromics. MVP is available at https://gitlab.com/ccoclet/mvp and as versioned packages in PyPi and Conda.

59 BASIC BIOLOGICAL SCIENCES↗