Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “clusters”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ion Clusters Reveal the Sources, Impacts, and Drivers of Freshwater Salinization

Population growth, land use change, climate change, and natural resource extraction are driving the salinization of freshwater resources worldwide. Reversing these trends will require data-centric approaches that identify salt sources, environmental drivers, and ecosystem responses. In this study, we applied principal component analysis and hierarchical clustering to identify ion covariance patterns, or “ion clusters,” in Broad Run, an urban stream in the Mid-Atlantic United States. These clusters correspond to distinct hydrologic regimes and reveal specific salinization risks: (1) phosphorus pollution mobilized during summer storms (Cluster 1); (2) elevated concentrations of sulfate and bicarbonate during baseflow (Cluster 2), likely reflecting groundwater discharge; and (3) elevated specific conductance and sodium, chloride, and potassium ion concentrations during snowmelt and rain-on-snow events (Cluster 3), driven by deicer and anti-icer wash-off. These ion fingerprints offer a transferable framework for diagnosing salt sources, assessing ecological risk, and identifying management targets. Our findings underscore the need for next-generation stormwater infrastructure and smart growth policies to protect aquatic life in rapidly urbanizing watersheds.

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

Microcarbonation of Naphthalene: An Experimental and Computational Study of Photoionization in Naphthalene-Carbon Dioxide Clusters

The photoionization of naphthalene (N)-carbon dioxide (CO 2 ) clusters was studied using tunable vacuum ultraviolet (VUV) radiation from a synchrotron in the photon range of 8.0 to 13.7 eV, in combination with time-of-flight mass spectrometry. Clusters of monomer, dimer, and trimer naphthalene with CO 2 (N­(CO 2 ) 0–6 , N 2 ­(CO 2 ) 0–3 , N3) were observed. The lowest-energy conformers were obtained via a conformer search, followed by geometry optimizations at the ωB97X-V2/aug-cc-pVTZ (monomer) and ωB97X-V2/aug-cc-pVDZ (dimer) levels of theory. Carbon dioxide was found to preferentially cluster on top of the naphthalene molecule (in an out-of-plane configuration). From the mass spectra, photoionization intensity curves (PICs) were constructed, and appearance energies (AEs) were determined. No substantial trend in AE was observed with increasing size of the naphthalene-carbon dioxide clusters; rather, AE oscillations around the value for pure naphthalene were observed. These AE oscillations are also observed in a recently studied naphthalene-water cluster system, though in this system a slight downward trend in AE for pure naphthalene clusters (N1/2/3/4) was observed. The differences between the two systems are attributed to differing interaction strengths. In conclusion, understanding these differences may aid in determining how photoprocessing can proceed differently depending on the dominant matrix component in interstellar ices.

Wannenmacher, Anna [Lawrence Berkeley National Lab↗
DOE OSTI2025

Ion Clusters in Multicomponent Solutions Determined from X-ray PDF and SAXS Analysis: The NaNO 2 –NaOH–H 2 O System

Here, this study explores ion cluster formation in the NaNO 2 –NaOH–H 2 O system to understand how ion cluster formation is influenced by the composition in multicomponent mixtures. X-ray pair distribution function (PDF) and small-angle X-ray scattering (SAXS) identified complex ion clusters in concentrated NaNO 2 and NaOH solutions as well as their mixtures. PDF analysis showed that the Na–O distance depended primarily on total Na + concentrations rather than anion composition, whereas the nitrite-water oxygen distance stayed the same regardless of concentration or composition. This result indicates that the nitrite ion hydration was relatively independent of composition. SAXS confirms local fluctuations and coherent clusters with notable size differences between NaOH and NaNO 2 solutions. SAXS analysis of mixed solutions shows that their clusters are an average of the individual solutions, indicating mixed clusters.

Reynolds, Jacob G. [Central Plateau Cleanup Compan↗
DOE OSTI2020

Supramolecular Assembly of U(IV) Clusters and Superatoms with Unconventional Countercations

Superatoms are nanometer-sized molecules or particles that form ordered lattices, mimicking their atomic counterparts. Hierarchical assembly of superatoms gives rise to emergent properties in lattices of quantum dots, p-block clusters, and fullerenes. In this report, we introduce a family of uranium-oxysulfate cluster anions whose hierarchical assembly in water is controlled by two parameters: acidity and the lanthanide or transition-metal countercation. In acid, larger Ln III (Ln = La–Ho) link hexamer (U 6 ) oxoclusters into body-centered cubic frameworks, while smaller Ln III (Ln = Er–Lu and Y) promote linking of 14 U 6 clusters into hollow superclusters (U 84 superatoms). U 84 assembles into superlattices including cubic-closest packed, body-centered cubic, and interpenetrating networks, bridged by interstitial countercations and U 6 clusters. Divalent transition metals (TM = Mn II and Zn II ) charge-balance and promote the fusion of 10 U 6 and 10 U monomers into a wheel-shaped cluster (U 70 ). Dissolution of U 70 in organic media shows (by small-angle X-ray scattering) that differing supramolecular assemblies are accessed, controlled by TM II -linking of U 70 clusters. Magnetic measurements of these assemblies reveal Curie–Weiss behavior at high temperatures, without pairing of the 5f 2 -electrons down to 2 K.

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

Examining Clustered MCSs and Their Precipitation Significance Over Global Land MCS Hotspots

Multiple independent mesoscale convective systems (MCSs) can occur closely in space and time, subjecting common geographic regions to multiple precipitation events over short periods. Motivated by an earlier study highlighting significant hydrological impacts of clustered MCSs in the US, we extend the analysis to seven global continental MCS hotspots using a 20‐year global MCS data set. Our findings reveal that (a) 40%–80% of MCSs occur in clusters, depending on the region; (b) clustered MCSs exhibit larger rainfall areas and higher mean areal rainfall, with the highest values observed in Argentina and India; and (c) clustered MCSs contribute increasing fractions to more extreme 24‐ and 72‐hr precipitation events across all regions. Our results underscore the critical role of clustered MCSs in driving hydrological extremes over global land regions and highlight the need to better understand how their formation mechanisms differ from non‐clustered events across the globe.

Hu, Huancui [Pacific Northwest National Laboratory↗
DOE OSTI2025

Computational ranking identifies Plexin-B2 in circulating tumor cell clustering with monocytes in breast cancer metastasis

Abstract Multicellular circulating tumor cell (CTC) clusters can be up to 50 times more efficient than single CTCs in mediating viable metastasis. Here, combining computational ranking and functional determination, we identify the transmembrane protein Plexin-B2 (PLXNB2) as one of the top molecular targets associated with unfavorable distant metastasis-free survival, showing enriched expression in CTC clusters versus single CTCs from patients with advanced breast cancer (mostly female). Loss of PLXNB2 (Plxnb2) reduces the formation of homotypic tumor cell clusters and heterotypic tumor-myeloid cell clusters, reducing spontaneous metastases in female mice bearing human (mouse) breast cancer. Interactions of PLXNB2 with its ligands SEMA4C on tumor cells and SEMA4A on myeloid cells (monocytes) promote homotypic and heterotypic CTC cluster formation, respectively, thereby driving lung metastasis. Global proteomic analysis reveals downstream effectors of the PLXNB2 pathway associated with tumor cell clustering. Thus, PLXNB2 is a therapeutic target for preventing new metastasis in breast cancer.

Science & Technology - Other Topics↗
DOE OSTI2022

Attosecond spectroscopy of size-resolved water clusters

Electron dynamics in water are of fundamental importance for a broad range of phenomena, but their real-time study faces numerous conceptual and methodological challenges. Here, in this paper, we introduce attosecond size-resolved cluster spectroscopy and build up a molecular-level understanding of the attosecond electron dynamics in water. We measure the effect that the addition of single water molecules has on the photoionization time delays of water clusters. We find a continuous increase of the delay for clusters containing up to four to five molecules and little change towards larger clusters. We show that these delays are proportional to the spatial extension of the created electron hole, which first increases with cluster size and then partially localizes through the onset of structural disorder that is characteristic of large clusters and bulk liquid water. These results indicate a previously unknown sensitivity of photoionization delays to electron-hole delocalization and indicate a direct link between electronic structure and attosecond photoionization dynamics. Our results offer new perspectives for studying electron-hole delocalization and its attosecond dynamics.

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

The minimal SUF system is not required for Fe–S cluster biogenesis in the methanogenic archaeon Methanosarcina acetivorans

Abstract Iron–sulfur (Fe–S) proteins are essential for the ability of methanogens to carry out methanogenesis and biological nitrogen fixation (diazotrophy). Nonetheless, the factors involved in Fe–S cluster biogenesis in methanogens remain largely unknown. The minimal SUF Fe–S cluster biogenesis system (i.e., SufBC) is postulated to serve as the primary system in methanogens. Here, the role of SufBC in Methanosarcina acetivorans , which contains two sufCB gene clusters, was investigated. The CRISPRi-dCas9 and CRISPR-Cas9 systems were utilized to repress or delete sufC1B1 and sufC2B2 , respectively. Neither the dual repression of sufC1B1 and sufC2B2 nor the deletion of both sufC1B1 and sufC2B2 affected the growth of M. acetivorans under any conditions tested, including diazotrophy. Interestingly, deletion of only sufC1B1 led to a delayed-growth phenotype under all growth conditions, suggesting that the deletion of sufC2B2 acts as a suppressor mutation in the absence of sufC1B1 . In addition, the deletion of sufC1B1 and/or sufC2B2 did not affect the total Fe–S cluster content in M. acetivorans cells. Overall, these results reveal that the minimal SUF system is not required for Fe–S cluster biogenesis in M. acetivorans and challenge the universal role of SufBC in Fe–S cluster biogenesis in methanogens.

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

Exciton energy transfer reveals spectral signatures of excited states in clusters

Electronic excitation and concomitant energy transfer leading to Penning ionization in argon–acetylene clusters generated in a supersonic expansion are investigated with synchrotron-based photoionization mass spectrometry and electronic structure calculations. Spectral features in the photoionization efficiency of the mixed argon–acetylene clusters reveal a blue shift from the 2 P 1/2 and 2 P 3/2 excited states of atomic argon. Analysis of this feature suggests that excited states of argon clusters transfer energy to acetylene, resulting in its ionization and successive evaporation of argon. Theoretically calculated Ar n (n = 2–6) cluster spectra are in excellent agreement with experimental observations, and provide insight into the structure and ionization dynamics of the clusters. A comparison between argon–acetylene and argon–water clusters reveals that argon solvates water better, allowing for higher-order excitons and Rydberg states to be populated. Furthermore, these results are explained by theoretical calculations of respective binding energies and structures.

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

Advances in the electron diffraction characterization of atomic clusters and nanoparticles

Nanoparticles and metallic clusters continue to make a remarkable impact on novel and emerging technologies. In recent years, there have been impressive advances in the controlled synthesis of clusters and their advanced characterization. One of the most common ways to determine the structures of nanoparticles and clusters is by means of X-ray diffraction methods. However, this requires the clusters to crystallize in a similar way to those used in protein studies, which is not possible in many cases. Novel methods based on electron diffraction have been used to efficiently study individual nanoparticles and clusters and these can overcome the obstacles commonly encountered during X-ray diffraction methods without the need for large crystals. These novel methodologies have improved with advances in electron microscopy instrumentation and electron detection. Here, we review advanced methodologies for characterizing metallic nanoparticles and clusters using a variety of electron diffraction procedures. These include selected area electron diffraction, nanobeam diffraction, coherent electron diffraction, precession electron diffraction, scanning transmission electron microcopy diffraction, and high throughput data analytics, which leverage deep learning to reduce the propensity for data errors and translate nanometer and atomic scale measurements into material data.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Size-dependent femtosecond proton transfer in protonated methanol clusters

Methanol clusters (CH 3 OH) n are prototypical hydrogen-bonded systems that bridge isolated molecules and bulk liquid methanol. Strong-field ionization of these clusters can trigger ultrafast charge migration, proton transfer, and isomerization; however, the onset timescales of these pathways remain largely unexplored despite extensive steady-state mass-spectrometric evidence for protonation and chemical rearrangement. Here we use femtosecond time-resolved strong-field ionization and disruptive probing to directly track the early-time dynamics of methanol clusters. Transient mass spectra yield size-dependent timescales for the formation of protonated clusters H + (CH 3 OH) n (n = 1–3) and an associated CHO˙ + radical cation. We find that both proton-transfer and subsequent stabilization become markedly faster with increasing cluster size. These measurements establish methanol clusters as a complementary benchmark to water for understanding ionization-initiated proton motion and chemical reorganization in hydrogen-bonded liquids.

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

The first high-redshift cavity power measurements of cool-core galaxy clusters with the International LOFAR Telescope

Radio-mode feedback associated with the active galactic nuclei (AGNs) at the cores of galaxy clusters injects a large amount of energy into the intracluster medium (ICM), offsetting radiative losses through X-ray emission. This mechanism prevents the ICM from rapidly cooling down and fueling extreme starburst activity as it accretes onto the central galaxies, and it is therefore a key ingredient in the evolution of galaxy clusters. However, the influence and mode of feedback at high redshifts (z ∼ 1) remains largely unknown. Low-frequency sub-arcsecond-resolution radio observations taken with the International LOFAR Telescope have demonstrated their ability to assist X-ray observations with constraining the energy output from the AGNs (or “cavity power”) in galaxy clusters, thereby enabling research at higher redshifts than before. In this pilot project, we tested this hybrid method on a high-redshift (0.6 < z < 1.3) sample of 13 galaxy clusters for the first time with the aim of verifying the performance of this method at these redshifts and providing the first estimates of the cavity power associated with the central AGN for a sample of distant clusters. We were able to detect clear radio lobes in three out of 13 galaxy clusters at redshifts of 0.7 < z < 0.9, and we used these detections in combination with ICM pressures surrounding the radio lobes obtained from standard profiles to calculate the corresponding cavity powers of the AGNs. Combining our results with the literature, the current data appear to suggest that the average cavity power peaked at a redshift ofz ∼ 0.4 and slowly decreases toward higher redshifts. However, we require more and tighter constraints on the cavity volume and a better understanding of our observational systematics to confirm any deviation of the cavity power trend from a constant level.

Astronomy & Astrophysics↗
DOE OSTI2024

Observation of a super-tetrahedral cluster of acetonitrile-solvated dodecaborate dianion via dihydrogen bonding

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B 12 H 12 2- ·nCH 3 CN (n = 1–4). The electron binding energies of B 12 H 12 2- ·nCH 3 CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B 12 H 12 2- via a threefold dihydrogen bond (DHB) B3–H3 ⋮⋮⋮ H3C–CN unit, in which three adjacent nucleophilic H atoms in B 12 H 12 2- interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B 12 H 12 2- solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B 12 H 12 2- and CH 3 CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B 12 H 12 2- ·CH 3 CN cluster. The binding motif in B 12 H 12 2- ·CH 3 CN is largely retained for the whole halogenated series B 12 X 12 2- ·CH 3 CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions.

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

Cosmological analysis of three-dimensional BOSS galaxy clustering and Planck CMB lensing cross correlations via Lagrangian perturbation theory

In this work, we present a formalism for jointly fitting pre- and post-reconstruction redshift-space clustering (RSD) and baryon acoustic oscillations (BAO) plus gravitational lensing (of the CMB) that works directly with the observed 2-point statistics. The formalism is based upon (effective) Lagrangian perturbation theory and a Lagrangian bias expansion, which models RSD, BAO and galaxy-lensing cross correlations within a consistent dynamical framework. As an example we present an analysis of clustering measured by the Baryon Oscillation Spectroscopic Survey in combination with CMB lensing measured by Planck. The post-reconstruction BAO strongly constrains the distance-redshift relation, the full-shape redshift-space clustering constrains the matter density and growth rate, and CMB lensing constrains the clustering amplitude. Using only the redshift space data we obtain Ω m = 0.303 ± 0.008, $H_0$ = 69.21 ± 0.78 and $σ_8$ = 0.743 ± 0.043. The addition of lensing information, even when restricted to the Northern Galactic Cap, improves constraints to Ω m = 0.303 ± 0.008, $H_0$ = 69.21 ± 0.77 and $σ_8$ = 0.707 ± 0.035, in tension with CMB and cosmic shear constraints. The combination of Ω m and $H_0$ are consistent with Planck, though their constraints derive mostly from redshift-space clustering. The low $σ_8$ value are driven by cross correlations with CMB lensing in the low redshift bin ($\textit{z}$ ≃ 0.38) and at large angular scales, which show a 20% deficit compared to expectations from galaxy clustering alone. We conduct several systematics tests on the data and find none that could fully explain these tensions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Weak Lensing Mass Calibration of the ACT DR5 Galaxy Clusters with the DES Year 3 Weak Lensing Data

We use weak gravitational lensing measurements from Year 3 Dark Energy Survey data to calibrate the masses of 443 galaxy clusters selected via the Sunyaev-Zel'dovich effect from Atacama Cosmology Telescope Data Release 5 maps of the cosmic microwave background. We incorporate redshift and SZ measurements for individual clusters into a hierarchical model for the stacked lensing signals and perform Bayesian analyses to constrain the hydrostatic mass bias of the clusters. Our treatment of systematic uncertainties includes a prescription for measuring and accounting for the weak lensing boost factor, consideration of a miscentering effect, as well as marginalization over uncertainties in the source galaxy photometric redshift distributions and shear calibration. The resultant constraints on the normalization of the mass-observable relation have a precision of approximately 7%, with the mean WL halo mass of M $_{500c}$ = 5.4 × 10$^{14}$ M $_{⊙}$. We measure the bias between the true cluster mass and the mass estimated from the SZ signal based on an X-ray-calibrated scaling relation assuming hydrostatic equilibrium, to be 1 - b = 0.74$^{+0.06}$ $_{-0.05}$ over the full sample. When splitting the clusters into high (z = 0.43-0.70) and low (z = 0.15-0.43) redshift bins, we measure 1 - b = 0.58$^{+0.06}$ $_{-0.05}$ and 0.81$^{+0.08}$ $_{-0.06}$, respectively. When introducing additional freedom in redshift and mass to the hydrostatic bias model, we find that 1 - b decreases with redshift (with the power law of -1.8$^{+0.5}$ $_{-0.6}$, 99.95% confidence), consistent with findings from other recent studies, while we do not find any significant trend in mass. We also demonstrate that our result is robust against various systematics such as a scale cut, priors on baryonic and miscentering parameters, and degree of scatter in mass-observable relation. The weak-lensing mass calibration presented in this study will be a useful tool for using the ACT clusters as probes of astrophysics, and as a step towards using their abundance as a cosmological probe.

Shin, T. [Carnegie Mellon U.] (ORCID:0000000263895↗
DOE OSTI2022

Deep structural clustering for single-cell RNA-seq data jointly through autoencoder and graph neural network

Abstract Single-cell RNA sequencing (scRNA-seq) permits researchers to study the complex mechanisms of cell heterogeneity and diversity. Unsupervised clustering is of central importance for the analysis of the scRNA-seq data, as it can be used to identify putative cell types. However, due to noise impacts, high dimensionality and pervasive dropout events, clustering analysis of scRNA-seq data remains a computational challenge. Here, we propose a new deep structural clustering method for scRNA-seq data, named scDSC, which integrate the structural information into deep clustering of single cells. The proposed scDSC consists of a Zero-Inflated Negative Binomial (ZINB) model-based autoencoder, a graph neural network (GNN) module and a mutual-supervised module. To learn the data representation from the sparse and zero-inflated scRNA-seq data, we add a ZINB model to the basic autoencoder. The GNN module is introduced to capture the structural information among cells. By joining the ZINB-based autoencoder with the GNN module, the model transfers the data representation learned by autoencoder to the corresponding GNN layer. Furthermore, we adopt a mutual supervised strategy to unify these two different deep neural architectures and to guide the clustering task. Extensive experimental results on six real scRNA-seq datasets demonstrate that scDSC outperforms state-of-the-art methods in terms of clustering accuracy and scalability. Our method scDSC is implemented in Python using the Pytorch machine-learning library, and it is freely available at https://github.com/DHUDBlab/scDSC.

Gan, Yanglan↗
DOE OSTI2025

Statistical relationships across epigenomes using large-scale hierarchical clustering

Recent advances in genomics and sequencing platforms have revolutionized our ability to create immense data sets, particularly for studying epigenetic regulation of gene expression. However, the avalanche of epigenomic data is difficult to parse for biological interpretation given nonlinear complex patterns and relationships. This attractive challenge in epigenomic data lends itself to machine learning for discerning infectivity and susceptibility. In this study, we explore over 3000 epigenomes of uninfected individuals and provide a framework to characterize the relationships among epigenetic modifiers, their modifiers, genetic loci, and specific immune cell types across all chromosomes using hierarchical clustering. Hierarchical clustering of epigenomic data revealed consistent epigenetic patterns across chromosomes, demonstrating that variation due to epigenetic modifiers is greater than variation between cell types. Gene Ontology and KEGG pathway analyses indicated significant enrichment of genes involved in chromatin remodeling, mRNA splicing, immune responses, and the regulation of microRNAs and snoRNAs. Epigenetic modifiers frequently formed biologically relevant clusters, including the cohesin complex, RNA Polymerase II transcription factors, and PRC2 complex members. These clustering behaviors remained consistent across all chromosomes, supported by entropy analysis and high Adjusted Rand Index scores, indicating robust cross-chromosomal similarity. Co-occurrence analysis further revealed specific sets of modifiers that consistently appeared together within clusters, reflecting shared biological functions and interactions. Validation using another dataset confirmed the reproducibility of these clustering patterns and modifier co-occurrence relationships, underscoring the reliability and generalizability of the methodology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Poisson hurdle model-based method for clustering microbiome features

Abstract Motivation High-throughput sequencing technologies have greatly facilitated microbiome research and have generated a large volume of microbiome data with the potential to answer key questions regarding microbiome assembly, structure and function. Cluster analysis aims to group features that behave similarly across treatments, and such grouping helps to highlight the functional relationships among features and may provide biological insights into microbiome networks. However, clustering microbiome data are challenging due to the sparsity and high dimensionality. Results We propose a model-based clustering method based on Poisson hurdle models for sparse microbiome count data. We describe an expectation–maximization algorithm and a modified version using simulated annealing to conduct the cluster analysis. Moreover, we provide algorithms for initialization and choosing the number of clusters. Simulation results demonstrate that our proposed methods provide better clustering results than alternative methods under a variety of settings. We also apply the proposed method to a sorghum rhizosphere microbiome dataset that results in interesting biological findings. Availability and implementation R package is freely available for download at https://cran.r-project.org/package=PHclust. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗