Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cluster analysis”

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,099 records · Page 61

Subsurface microbial communities as a tool for characterizing regional-scale groundwater flow

Subsurface microbial community distribution patterns are influenced by biogeochemical and groundwater fluxes and may inform hydraulic connections along groundwater-flow paths. This study examined the regional-scale microbial community of the Death Valley Regional Flow System and evaluated whether subsurface communities can be used to identify groundwater-flow paths between recharge and discharge areas. Samples were collected from 36 sites in three groundwater basins: Pahute Mesa–Oasis Valley (PMOV), Ash Meadows (AM), and Alkali Flat–Furnace Creek Ranch (AFFCR). Microbial diversity within and between communities varied by location, and communities were separated into two overall groups that affiliated with the AM and PMOV/AFFCR basins. Network analysis revealed patterns between clusters of common microbes that represented groundwaters with similar geochemical conditions and largely corroborated hydraulic connections between recharge and discharge areas. Null model analyses identified deterministic and stochastic ecological processes contributing to microbial community assemblages. Most communities were more different than expected and governed by dispersal limitation, geochemical differences, or undominating processes. However, certain communities from sites located within or near the Nevada National Security Site were more similar than expected and dominated by homogeneous dispersal or selection. Overall, the (dis)similarities between the microbial communities of DVRFS recharge and discharge areas supported previously documented hydraulic connections between: (1) Spring Mountains and Ash Meadows; (2) Frenchman and Yucca Flat and Amargosa Desert; and (3) Amargosa Desert and Death Valley. However, only a portion of the flow path between Pahute Mesa and Oasis Valley could be supported by microbial community analyses, likely due to well-associated artifacts in samples from the two Oasis Valley sites. This study demonstrates the utility of combining microbial data with hydrologic, geologic, and water-chemistry information to comprehensively characterize groundwater systems, highlighting both strengths and limitations of this approach.

54 ENVIRONMENTAL SCIENCES↗

Thiohalorhabdus methylotropha sp. nov., an extremely halophilic autotrophic methylotiotroph from hypersaline lakes

So far, there have been no reports of trimethylamine (TMA)-utilizing extremely halophilic microorganisms in hypersaline habitats. Our aerobic enrichments at 4 M total Na + with 5 mM TMA inoculated with surface sediments from hypersaline soda (at pH 9.5) or chloride-sulfate (at pH 7) lakes in southwestern Siberia were successful only for the latter. The initial enrichment included both bacteria and haloarchaea but only the bacterial component was able to grow as a pure culture with TMA. Strain Cl-TMA forms a new-species lineage within the genus Thiohalorhabdus which includes extremely halophilic and obligate lithoautotrophic sulfur-oxidizing gammaproteobacteria. Cl-TMA can grow methyloautotrophically utilizing TMA, dimethylamine (DMA) and methanol (MeOH) as the electron donors or chemolithoautotrophically with thiosulfate. Mixotrophic growth was also observed with the three methyl compounds and thiosulfate. Carbon is assimilated autotrophically via the Calvin-Benson-Basham pathway. Unlike the type species of Thiohalorhabdus, T. denitrificans , Cl-TMA was incapable of anaerobic growth via denitrification. The isolate belongs to extreme halophiles growing between 2.5 and 5 M NaCl with an optimum at 3–3.5 M. Genome analysis identified two gene clusters coding for PQQ-dependent methanol dehydrogenases (MxaFI and XoxF), four homologues of the formaldehyde activating enzymes (Faes), a TMA/DMA oxidation locus, and two cluster of genes encoding an N-methylglutamate dehydrogenase pathway (NMGP) for methylamine oxidation. The first steps of C 1 -subtrate conversions are followed by the tetrahydrofolate (THF)-linked and tetrahydromethanopterin (H4MPT)-linked formaldehyde oxidation pathways and two formate dehydrogenases. All of those signatures of methylotrophy were absent in T. denitrificans . In contrast, genes for two key sulfur oxidation enzymes, thiosulfate dehydrogenase TsdAB and sulfide dehydrogenase FccAB, that are present in the type species are missing in Cl-TMA. Thiosulfate is oxidized to sulfate by a combination of an incomplete Sox cycle and an sHdr system. Strain Cl-TMA T (JCM 35977 = UQM 41915) is proposed to be classified as Thiohalorhabdus methylotrophus sp. nov.

59 BASIC BIOLOGICAL SCIENCES↗

The miR-183 family cluster alters zinc homeostasis in benign prostate cells, organoids and prostate cancer xenografts

The miR-183 cluster, which is comprised of paralogous miRs-183, -96 and -182, is overexpressed in many cancers, including prostate adenocarcinoma (PCa). Prior studies showed that overexpression of individual pre-miRs-182, -96 and -183 in prostate cells decreased zinc import, which is a characteristic feature of PCa tumours. Zinc is concentrated in healthy prostate 10-fold higher than any other tissue, and an >80% decrease in zinc is observed in PCa specimens. Here, we studied the effect of overexpression of the entire 4.8 kb miR-183 family cluster, including the intergenic region which contains highly conserved genomic regions, in prostate cells. This resulted in overexpression of mature miR-183 family miRs at levels that mimic cancer-related changes. Overexpression of the miR-183 cluster reduced zinc transporter and intracellular zinc levels in benign prostate cells, PCa xenografts and fresh prostate epithelial organoids. Microarray analysis of miR-183 family cluster overexpression in prostate cells showed an enrichment for cancer-related pathways including adhesion, migration and wound healing. An active secondary transcription start site was identified within the intergenic region of the miR-183 cluster, which may regulate expression of miR-182. Taken together, this study shows that physiologically relevant expression of the miR-183 family regulates zinc levels and carcinogenic pathways in prostate cells.

Science & Technology - Other Topics↗

Crystal genes in a marginal glass-forming system of Ni 50 Zr 50

Glass-forming motifs with B2 traits are found. A perfect Ni-centered B33 motif deteriorates the glass-forming ability of Ni 50 Zr 50 . The marginal glass-forming ability (GFA) of binary Ni-Zr system is an issue to be explained considering the numerous bulk metallic glasses (BMGs) found in the Cu-Zr system. Using molecular dynamics, the structures and dynamics of Ni 50 Zr 50 metallic liquid and glass are investigated at the atomistic level. To achieve a well-relaxed glassy sample, sub-T g annealing method is applied and the final sample is closer to the experiments than the models prepared by continuous cooling. With the state-of-the-art structural analysis tools such as cluster alignment and pair-wise alignment methods, two glass-forming motifs with some mixed traits of the metastable B2 crystalline phase and the crystalline Ni-centered B33 motif are found to be dominant in the undercooled liquid and glass samples. A new chemical order characterization on each short-range order (SRO) structure is accomplished based on the cluster alignment method. The significant amount of the crystalline motif and the few icosahedra in the glassy sample deteriorate the GFA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Communication: Newton homotopies for sampling stationary points of potential energy landscapes

One of the most challenging and frequently arising problems in many areas of science is to find solutions of a system of multivariate nonlinear equations. There are several numerical methods that can find many (or all if the system is small enough) solutions but they all exhibit characteristic problems. Moreover, traditional methods can break down if the system contains singular solutions. Here, we propose an efficient implementation of Newton homotopies, which can sample a large number of the stationary points of complicated many-body potentials. We demonstrate how the procedure works by applying it to the nearest-neighbor ϕ{sup 4} model and atomic clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Delocalization and stretch-bend mixing of the HOH bend in liquid water

Liquid water’s rich sub-picosecond vibrational dynamics arise from the interplay of different high- and low-frequency modes evolving in a strong yet fluctuating hydrogen bond network. Recent studies of the OH stretching excitations of H 2 O indicate that they are delocalized over several molecules, raising questions about whether the bending vibrations are similarly delocalized. In this paper, we take advantage of an improved 50 fs time-resolution and broadband infrared (IR) spectroscopy to interrogate the 2D IR lineshape and spectral dynamics of the HOH bending vibration of liquid H2O. Indications of strong bend-stretch coupling are observed in early time 2D IR spectra through a broad excited state absorption that extends from 1500 cm -1 to beyond 1900 cm -1 , which corresponds to transitions from the bend to the bend overtone and OH stretching band between 3150 and 3550 cm -1 . Pump-probe measurements reveal a fast 180 fs vibrational relaxation time, which results in a hot-ground state spectrum that is the same as observed for water IR excitation at any other frequency. The fastest dynamical time scale is 80 fs for the polarization anisotropy decay, providing evidence for the delocalized or excitonic character of the bend. Finally, normal mode analysis conducted on water clusters extracted from molecular dynamics simulations corroborate significant stretch-bend mixing and indicate delocalization of δ HOH on 2-7 water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cooling rate dependence of structural order in Ni 62 Nb 38 metallic glass

In this article, molecular dynamics (MD) simulations are performed to study the structure of Ni 62 Nb 38 bulk metallic glass at the atomistic level. Structural analysis based on the cluster alignment method is carried out and a new Ni-centered distorted-icosahedra (DISICO) motif is excavated. We show that the short-range order and medium-range order in the glass are enhanced with lower cooling rate. Almost 50% of the clusters around the Ni atoms in the well-annealed Ni 62 Nb 38 glass sample from our MD simulations can be classified as DISICO. It is revealed that the structural distortion with respect to the perfect icosahedra is driven by chemical ordering in the distorted region of the DISICO motif. The relationship between the structure, energy, and dynamics in this glass-forming alloy during the cooling and annealing processes is also established.

36 MATERIALS SCIENCE↗

STAR FORMATION ACROSS THE W3 COMPLEX

We present a multi-wavelength analysis of the history of star formation in the W3 complex. Using deep, near-infrared ground-based images combined with images obtained with Spitzer and Chandra observatories, we identified and classified young embedded sources. We identified the principal clusters in the complex and determined their structure and extension. We constructed extinction-limited samples for five principal clusters and constructed K-band luminosity functions that we compare with those of artificial clusters with varying ages. This analysis provided mean ages and possible age spreads for the clusters. We found that IC 1795, the centermost cluster of the complex, still hosts a large fraction of young sources with circumstellar disks. This indicates that star formation was active in IC 1795 as recently as 2 Myr ago, simultaneous to the star-forming activity in the flanking embedded clusters, W3-Main and W3(OH). A comparison with carbon monoxide emission maps indicates strong velocity gradients in the gas clumps hosting W3-Main and W3(OH) and shows small receding clumps of gas at IC 1795, suggestive of rapid gas removal (faster than the T Tauri timescale) in the cluster-forming regions. We discuss one possible scenario for the progression of cluster formation in the W3 complex. We propose that early processes of gas collapse in the main structure of the complex could have defined the progression of cluster formation across the complex with relatively small age differences from one group to another. However, triggering effects could act as catalysts for enhanced efficiency of formation at a local level, in agreement with previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Bias due to neutrinos must not uncorrect'd go

It is a well known fact that galaxies are biased tracers of the distribution of matter in the Universe. The galaxy bias is usually factored as a function of redshift and scale, and approximated as being scale-independent on large, linear scales. In cosmologies with massive neutrinos, the galaxy bias defined with respect to the total matter field (cold dark matter, baryons, and non-relativistic neutrinos) also depends on the sum of the neutrino masses M ν , and becomes scale-dependent even on large scales. This effect has been usually neglected given the sensitivity of current surveys. Yet, it becomes a severe systematic for future surveys aiming to provide the first detection of non-zero M ν . The effect can be corrected for by defining the bias with respect to the density field of cold dark matter and baryons, rather than the total matter field. In this work, we provide a simple prescription for correctly mitigating the neutrino-induced scale-dependent bias effect in a practical way. We clarify a number of subtleties regarding how to properly implement this correction in the presence of redshift-space distortions and non-linear evolution of perturbations. We perform a Markov Chain Monte Carlo analysis on simulated galaxy clustering data that match the expected sensitivity of the Euclid survey. We reveal that the neutrino-induced scale-dependent bias can lead to important shifts in both the inferred mean value of M ν , as well as its uncertainty, and provide an analytical explanation for the magnitude of the shifts. We show how these shifts propagate to the inferred values of other cosmological parameters correlated with M ν , such as the cold dark matter physical density Ω cdm h 2 and the scalar spectral index n s . In the end, we find that correctly accounting for the neutrino-induced scale-dependent bias will be of crucial importance for future galaxy clustering analyses. We encourage the cosmology community to correctly account for this effect using the simple prescription we present in our work. The tools necessary to easily correct for the neutrino-induced scale-dependent bias will be made publicly available in an upcoming release of the Boltzmann solver CLASS.

79 ASTRONOMY AND ASTROPHYSICS↗

Hybrid deep learning architecture for general disruption prediction across tokamaks

In this paper, we present a new deep learning disruption prediction algorithm based on important findings from explorative data analysis which effectively allows knowledge transfer from existing devices to new ones, thereby predicting disruptions using very limited disruptive data from the new devices. Here, the explorative data analysis conducted via unsupervised clustering techniques confirms that time-sequence data are much better separators of disruptive and non-disruptive behavior than the instantaneous plasma state data with further advantageous implications for a sequence-based predictor. Based on such important findings, we have designed a new algorithm for multi-machine disruption prediction that achieves high predictive accuracy on the C-Mod (AUC=0.801), DIII-D (AUC=0.947) and EAST (AUC=0.973) tokamaks with limited hyperparameter tuning. Through numerical experiments, we show that boosted accuracy (AUC=0.959) is achieved on EAST predictions by including in the training only 20 disruptive discharges, thousands of non-disruptive discharges from EAST, and combining this with more than a thousand discharges from DIII-D and C-Mod. The improvement of predictive ability obtained by combining disruptive data from other devices is found to be true for all permutations of the three devices. Furthermore, by comparing the predictive performance of each individual numerical experiment, we find that non-disruptive data are machine-specific while disruptive data from multiple devices contain device-independent knowledge that can be used to inform predictions for disruptions occurring on a new device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A HIGHLY ELONGATED PROMINENT LENS AT z = 0.87: FIRST STRONG-LENSING ANALYSIS OF EL GORDO

We present the first strong-lensing (SL) analysis of the galaxy cluster ACT-CL J0102-4915 (El Gordo), in recent HST/ACS images, revealing a prominent strong lens at a redshift of z = 0.87. This finding adds to the already-established unique properties of El Gordo: it is the most massive, hot, X-ray luminous, and bright Sunyaev-Zeldovich effect cluster at z {approx}> 0.6, and the only {sup b}ullet{sup -}like merging cluster known at these redshifts. The lens consists of two merging massive clumps, where, for a source redshift of z{sub s} {approx} 2, each clump exhibits only a small, separate critical area, with a total area of 0.69 {+-} 0.11{open_square}' over the two clumps. For a higher source redshift, z{sub s} {approx} 4, the critical curves of the two clumps merge together into one bigger and very elongated lens (axis ratio {approx_equal} 5.5), enclosing an effective area of 1.44 {+-} 0.22{open_square}'. The critical curves continue expanding with increasing redshift so that for high-redshift sources (z{sub s} {approx}> 9) they enclose an area of {approx}1.91 {+-} 0.30{open_square}' (effective {theta}{sub e} {approx_equal} 46.''8 {+-} 3.''7) and a mass of 6.09 {+-} 1.04 Multiplication-Sign 10{sup 14} M{sub Sun }. According to our model, the area of high magnification ({mu} > 10) for such high-redshift sources is {approx_equal}1.2{open_square}', and the area with {mu} > 5 is {approx_equal}2.3{open_square}', making El Gordo a compelling target for studying the high-redshift universe. We obtain a strong lower limit on the total mass of El Gordo, {approx}> 1.7 Multiplication-Sign 10{sup 15} M{sub Sun} from the SL regime alone, suggesting a total mass of roughly M{sub 200} {approx} 2.3 Multiplication-Sign 10{sup 15} M{sub Sun }. Our results should be revisited when additional spectroscopic and HST imaging data are available.

79 ASTRONOMY AND ASTROPHYSICS↗

HH 1158: THE LOWEST LUMINOSITY EXTERNALLY IRRADIATED HERBIG–HARO JET

We have identified a new externally irradiated Herbig–Haro (HH) jet, HH 1158, within ∼2 pc of the massive OB type stars in the σ Orionis cluster. At an L{sub bol} ∼ 0.1 L{sub ⊙}, HH 1158 is the lowest luminosity irradiated HH jet identified to date in any cluster. Results from the analysis of high-resolution optical spectra indicate asymmetries in the brightness, morphology, electron density, velocity, and the mass outflow rates for the blue and redshifted lobes. We constrain the position angle of the HH 1158 jet at 102° ± 5°. The mass outflow rate and the mean accretion rate for HH 1158 using multiple diagnostics are estimated to be (5.2 ± 2.6) × 10{sup −10} M{sub ⊙} yr{sup −1} and (3.0 ± 1.0) × 10{sup −10} M{sub ⊙} yr{sup −1}, respectively. The properties for HH 1158 are notably similar to the externally irradiated HH 444–HH 447 jets previously identified in σ Orionis. In particular, the morphology is such that the weaker jet beam is tilted toward the massive stars, indicating a higher extent of photo-evaporation. The high value for the Hα/[S ii] ratio is also consistent with the ratios measured in other irradiated jets, including HH 444–HH 447. The presence of an extended collimated jet that is bipolar and the evidence of shocked emission knots make HH 1158 the first unique case of irradiated HH jets at the very low-luminosity end, and provides an opportunity to learn the physical properties of very faint HH jet sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Microbial and viral-like rhodopsins present in coastal marine sediments from four polar and subpolar regions

Rhodopsins are broadly distributed. In this work, we analyzed 23 metagenomes corresponding to marine sediment samples from four regions that share cold climate conditions (Norway; Sweden; Argentina and Antarctica). In order to investigate the genes evolution of viral rhodopsins, an initial set of 6224 bacterial rhodopsin sequences according to COG5524 were retrieved from the 23 metagenomes. After selection by the presence of transmembrane domains and alignment, 123 viral (51) and non-viral (72) sequences (>50 amino acids) were finally included in further analysis. Viral rhodopsin genes were homologs of Phaeocystis globosa virus and Organic lake Phycodnavirus. Non-viral microbial rhodopsin genes were ascribed to Bacteroidetes, Planctomycetes, Firmicutes, Actinobacteria, Cyanobacteria, Proteobacteria, Deinococcus-Thermus and Cryptophyta and Fungi. A rescreening using Blastp, using as queries the viral sequences previously described, retrieved 30 sequences (>100 amino acids). Phylogeographic analysis revealed a geographical clustering of the sequences affiliated to the viral group. This clustering was not observed for the microbial non-viral sequences. The phylogenetic reconstruction allowed us to propose the existence of a putative ancestor of viral rhodopsin genes related to Actinobacteria and Chloroflexi. This is the first report about the existence of a phylogeographic association of the viral rhodopsin sequences from marine sediments.

marine sediment↗

AutoEnRichness: A hybrid empirical and analytical approach for estimating the richness of galaxy clusters

ABSTRACT We introduce AutoEnRichness, a hybrid approach that combines empirical and analytical strategies to determine the richness of galaxy clusters (in the redshift range of 0.1 ≤ z ≤ 0.35) using photometry data from the Sloan Digital Sky Survey Data Release 16, where cluster richness can be used as a proxy for cluster mass. In order to reliably estimate cluster richness, it is vital that the background subtraction is as accurate as possible when distinguishing cluster and field galaxies to mitigate severe contamination. AutoEnRichness is comprised of a multistage machine learning algorithm that performs background subtraction of interloping field galaxies along the cluster line of sight and a conventional luminosity distribution fitting approach that estimates cluster richness based only on the number of galaxies within a magnitude range and search area. In this proof-of-concept study, we obtain a balanced accuracy of 83.20 per cent when distinguishing between cluster and field galaxies as well as a median absolute percentage error of 33.50 per cent between our estimated cluster richnesses and known cluster richnesses within r200. In the future, we aim for AutoEnRichness to be applied on upcoming large-scale optical surveys, such as the Legacy Survey of Space and Time and Euclid, to estimate the richness of a large sample of galaxy groups and clusters from across the halo mass function. This would advance our overall understanding of galaxy evolution within overdense environments as well as enable cosmological parameters to be further constrained.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimating cluster masses from SDSS multiband images with transfer learning

ABSTRACT The total masses of galaxy clusters characterize many aspects of astrophysics and the underlying cosmology. It is crucial to obtain reliable and accurate mass estimates for numerous galaxy clusters over a wide range of redshifts and mass scales. We present a transfer-learning approach to estimate cluster masses using the ugriz-band images in the SDSS Data Release 12. The target masses are derived from X-ray or SZ measurements that are only available for a small subset of the clusters. We designed a semisupervised deep learning model consisting of two convolutional neural networks. In the first network, a feature extractor is trained to classify the SDSS photometric bands. The second network takes the previously trained features as inputs to estimate their total masses. The training and testing processes in this work depend purely on real observational data. Our algorithm reaches a mean absolute error (MAE) of 0.232 dex on average and 0.214 dex for the best fold. The performance is comparable to that given by redMaPPer, 0.192 dex. We have further applied a joint integrated gradient and class activation mapping method to interpret such a two-step neural network. The performance of our algorithm is likely to improve as the size of training data set increases. This proof-of-concept experiment demonstrates the potential of deep learning in maximizing the scientific return of the current and future large cluster surveys.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Cosmology with the Roman Space Telescope – Synergies with CMB lensing

ABSTRACT We explore synergies between the Nancy Grace Roman Space Telescope and CMB lensing data to constrain dark energy and modified gravity scenarios. A simulated likelihood analysis of the galaxy clustering and weak lensing data from the Roman Space Telescope High Latitude Survey combined with CMB lensing data from the Simons Observatory is undertaken, marginalizing over important astrophysical effects and calibration uncertainties. Included in the modelling are the effects of baryons on small-scale clustering, scale-dependent growth suppression by neutrinos, as well as uncertainties in the galaxy clustering biases, in the intrinsic alignment contributions to the lensing signal, in the redshift distributions, and in the galaxy shape calibration. The addition of CMB lensing roughly doubles the dark energy figure-of-merit from Roman photometric survey data alone, varying from a factor of 1.7 to 2.4 improvement depending on the particular Roman survey configuration. Alternatively, the inclusion of CMB lensing information can compensate for uncertainties in the Roman galaxy shape calibration if it falls below the design goals. Furthermore, we report the first forecast of Roman constraints on a model-independent structure growth, parametrized by σ8(z), and on the Hu-Sawicki f(R) gravity as well as an improved forecast of the phenomenological (Σ0, μ0) model. We find that CMB lensing plays a crucial role in constraining σ8(z) at z > 2, with percent-level constraints forecasted out to z = 4. CMB lensing information does not improve constraints on the f(R) model substantially. It does, however, increase the (Σ0, μ0) figure-of-merit by a factor of about 1.5.

79 ASTRONOMY AND ASTROPHYSICS↗

A measurement of CMB cluster lensing with SPT and DES year 1 data

Clusters of galaxies gravitationally lens the cosmic microwave background (CMB) radiation, resulting in a distinct imprint in the CMB on arcminute scales. Measurement of this effect offers a promising way to constrain the masses of galaxy clusters, particularly those at high redshift. We use CMB maps from the South Pole Telescope Sunyaev-Zel'dovich (SZ) survey to measure the CMB lensing signal around galaxy clusters identified in optical imaging from first year observations of the Dark Energy Survey. The cluster catalog used in this analysis contains 3697 members with mean redshift of $\bar{z} = 0.45$. We detect lensing of the CMB by the galaxy clusters at $8.1\sigma$ significance. Using the measured lensing signal, we constrain the amplitude of the relation between cluster mass and optical richness to roughly $17\%$ precision, finding good agreement with recent constraints obtained with galaxy lensing. The error budget is dominated by statistical noise but includes significant contributions from systematic biases due to the thermal SZ effect and cluster miscentering.

79 ASTRONOMY AND ASTROPHYSICS↗

The parallelism motifs of genomic data analysis

Genomic datasets are growing dramatically as the cost of sequencing continues to decline and small sequencing devices become available. Enormous community databases store and share these data with the research community, but some of these genomic data analysis problems require large-scale computational platforms to meet both the memory and computational requirements. These applications differ from scientific simulations that dominate the workload on high-end parallel systems today and place different requirements on programming support, software libraries and parallel architectural design. For example, they involve irregular communication patterns such as asynchronous updates to shared data structures. We consider several problems in high-performance genomics analysis, including alignment, profiling, clustering and assembly for both single genomes and metagenomes. We identify some of the common computational patterns or ‘motifs’ that help inform parallelization strategies and compare our motifs to some of the established lists, arguing that at least two key patterns, sorting and hashing, are missing. This article is part of a discussion meeting issue ‘Numerical algorithms for high-performance computational science’.

59 BASIC BIOLOGICAL SCIENCES↗