Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Phylogenetic analysis of anaerobic psychrophilic enrichment cultures obtained from a greenland glacier ice core

The examination of microorganisms in glacial ice cores allows the phylogenetic relationships of organisms frozen for thousands of years to be compared with those of current isolates. We developed a method for aseptically sampling a sediment-containing portion of a Greenland ice core that had remained at -9 degrees C for over 100,000 years. Epifluorescence microscopy and flow cytometry results showed that the ice sample contained over 6 x 10(7) cells/ml. Anaerobic enrichment cultures inoculated with melted ice were grown and maintained at -2 degrees C. Genomic DNA extracted from these enrichments was used for the PCR amplification of 16S rRNA genes with bacterial and archaeal primers and the preparation of clone libraries. Approximately 60 bacterial inserts were screened by restriction endonuclease analysis and grouped into 27 unique restriction fragment length polymorphism types, and 24 representative sequences were compared phylogenetically. Diverse sequences representing major phylogenetic groups including alpha, beta, and gamma Proteobacteria as well as relatives of the Thermus, Bacteroides, Eubacterium, and Clostridium groups were found. Sixteen clone sequences were closely related to those from known organisms, with four possibly representing new species. Seven sequences may reflect new genera and were most closely related to sequences obtained only by PCR amplification. One sequence was over 12% distant from its closest relative and may represent a novel order or family. These results show that phylogenetically diverse microorganisms have remained viable within the Greenland ice core for at least 100,000 years.

Evaluation Studies↗
DOE OSTI2021

Comprehensive analysis of differentially expressed mRNA and circRNA in Ankylosing spondylitis patients’ platelets

Highlights: • 4996 mRNAs and 2942 circRNAs were identified differentially expressed in platelets of AS patients. • Several platelet-derived immune mediators were upregulated to regulate the interaction between platelets and immune cells. • Two downregulated circRNA were identified, and the corresponding circRNA-miRNA-mRNA regulatory network was constructed. • mRNAs and circRNAs found herein could be used as diagnostic biomarkers or new therapeutic targets. Ankylosing spondylitis (AS) is a chronic inflammatory disease significantly decreasing the quality of life. Platelets play an important and active role in the development of AS. Accumulating evidence demonstrated platelets contain diverse RNA repository inherited from megakaryocytes or microvesicles. Platelet RNAs are dynamically affected by pathological conditions and could be used as diagnostic or prognostic biomarkers. However, the role of the platelet RNAs in AS is elusive. In this study, we compared mRNA and circRNA profiles in platelets between AS patients and healthy controls using RNA sequencing and bioinformatic analysis, and found 4996 mRNAs and 2942 circRNAs were differently expressed. The significantly over-expressed mRNAs in AS patients are involved in platelet activity, gap junction, focal adhesion, rap1 and toll and Imd signaling pathway. The previous identified platelet-derived immune mediators such as P2Y1, P2Y12, PF4, GPIbα, CD40L, ICAM2, CCL5 (RANTES), TGF-β (TGF-β1 and TGF-β2) and PDGF (PDGFB and PDGFA) are also included in these over expressed mRNAs, implying these factors may trigger inflammatory cascades and promote the development of AS. Additionally, we found two down-regulated circRNA (circPTPN22 and circFCHSD2) from the intersection analyses of platelets and spinal ligament tissues of AS patients. The circRNA-miRNA-mRNA regulatory network of these two circRNAs was constructed, and the target mRNAs were enriched in Th17 cell differentiation, inflammatory bowel disease, cell adhesion molecules, cytokine-cytokine receptor interaction, Jak-STAT and Wnt signaling pathway, all these pathways participate in the bone remodeling and pro/anti-inflammatory immune regulation in AS. Then, qRT-PCR was performed to validate the expression of selected key mRNAs and circRNAs and the results demonstrated that the expression levels of P2Y12, GPIbα, circPTPN22 and circFCHSD2 were consistent with the sequencing analysis. In addition, the high expression of five predicted miRNAs interacting with circPTPN22 and circFCHSD2 were also detected in AS by qRT-PCR. Taken together, our study presents a comprehensive overview of mRNAs and circRNAs in platelets in AS patients and offers new insight into the mechanisms of platelet involving in the pathogenesis of AS. The mRNAs and circRNAs identified in this study may serve as candidates for diagnosis and targeted treatment of AS.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

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

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

Shao, Nan↗
DOE OSTI2026

Data for High Yield Production of 3-Hydroxypropionic Acid Using Issatchenkia orientalis

Biomanufacturing provides a more sustainable alternative to fossil-based chemical manufacturing. 3-Hydroxypropionic acid (3HP) is a top Department of Energy value-added chemical and precursor to bioplastics, yet cost-effective microbial production remains elusive. Here, we establish the acid-tolerant yeast Issatchenkia orientalis as a robust host for low-pH 3HP biosynthesis. Genome-scale modeling identifies the β-alanine pathway as optimal, offering the highest theoretical yield and lowest oxygen requirement. Thermodynamic analysis confirms its favorability under acidic conditions. Using sequence similarity network analysis, we discover highly active aspartate 1-decarboxylase (PAND), β-alanine-pyruvate aminotransferase (BAPAT), and 3HP dehydrogenase (YDFG), which significantly improve the pathway efficiency. Next, to further elevate the production, pathway optimization through multi-copy PAND integration, byproduct elimination (knockouts of pyruvate decarboxylase and glycerol-3-phosphate dehydrogenase), and reinforcement of aspartate flux by overexpression of pyruvate carboxylase and aspartate amino transferase improves the titer to 29 g/L in shake flasks. Fed-batch fermentation at pH 4 with low-cost corn steep liquor medium further increases the production to 92 g/L with 0.7 g/g yield and 0.55 g/L/h productivity. Techno-economic analysis indicates that such performance could potentially enable a financially viable process for sustainable acrylic acid production. This work establishes I. orientalis as a next-generation platform for cost-effective 3HP production and paves the way toward industrial commercialization.

Bioproducts↗
DOE OSTI2024

Sequence, structure prediction, and epitope analysis of the polymorphic membrane protein family in Chlamydia trachomatis

The polymorphic membrane proteins (Pmps) are a family of autotransporters that play an important role in infection, adhesion and immunity in Chlamydia trachomatis. Here we show that the characteristic GGA(I,L,V) and FxxN tetrapeptide repeats fit into a larger repeat sequence, which correspond to the coils of a large beta-helical domain in high quality structure predictions. Analysis of the protein using structure prediction algorithms provided novel insight to the chlamydial Pmp family of proteins. While the tetrapeptide motifs themselves are predicted to play a structural role in folding and close stacking of the beta-helical backbone of the passenger domain, we found many of the interesting features of Pmps are localized to the side loops jutting out from the beta helix including protease cleavage, host cell adhesion, and B-cell epitopes; while T-cell epitopes are predominantly found in the beta-helix itself. This analysis more accurately defines the Pmp family of Chlamydia and may better inform rational vaccine design and functional studies.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2001

Ecological physiology of Synechococcus sp. strain SH-94-5, a naturally occurring cyanobacterium deficient in nitrate assimilation

Synechococcus sp. strain SH-94-5 is a nitrate assimilation-deficient cyanobacterium which was isolated from an ammonium-replete hot spring in central Oregon. While this clone could grow on ammonium and some forms of organic nitrogen as sole nitrogen sources, it could not grow on either nitrate or nitrite, even under conditions favoring passive diffusion. It was determined that this clone does not express functional nitrate reductase or nitrite reductase and that the lack of activity of either enzyme is not due to inactivation of the cyanobacterial nitrogen control protein NtcA. A few other naturally occurring cyanobacterial strains are also nitrate assimilation deficient, and phylogenetic analyses indicated that the ability to utilize nitrate has been independently lost at least four times during the evolutionary history of the cyanobacteria. This phenotype is associated with the presence of environmental ammonium, a negative regulator of nitrate assimilation gene expression, which may indicate that natural selection to maintain functional copies of nitrate assimilation genes has been relaxed in these habitats. These results suggest how the evolutionary fates of conditionally expressed genes might differ between environments and thereby effect ecological divergence and biogeographical structure in the microbial world.

Nitrates/metabolism↗
NASA NTRS2004

Molecular evidence for a terrestrial origin of snakes

Biologists have debated the origin of snakes since the nineteenth century. One hypothesis suggests that snakes are most closely related to terrestrial lizards, and reduced their limbs on land. An alternative hypothesis proposes that snakes are most closely related to Cretaceous marine lizards, such as mosasaurs, and reduced their limbs in water. A presumed close relationship between living monitor lizards, believed to be close relatives of the extinct mosasaurs, and snakes has bolstered the marine origin hypothesis. Here, we show that DNA sequence evidence does not support a close relationship between snakes and monitor lizards, and thus supports a terrestrial origin of snakes.

Environment↗
DOE OSTI2021

Physical, resource supply, and biological controls on nutrient processing along the river continuum

Nutrient impairment has led to damages to US surface and groundwater systems in excess of 100 billion dollars per year. Therefore, there is a strong need to develop methods to predict the transport, uptake, and export of nutrients along fluvial networks. We present results that are based on a data-driven mechanistic understanding of three factors that largely control nutrient uptake and export: 1) interactions between transport-related processes (mass transfer to metabolically active zones), 2) resource supply dynamics (nutrient concentration, stoichiometric constraints, etc.), and 3) biological controls (microbial community structure and function). Our results were generated from column experiments conducted along the Jemez River-Rio Grande continuum, which spans four orders of magnitude in mean annual discharge, more than 2000 m in altitude, and more than 500 km of stream longitude. Two resource supply injections were performed on each of the columns, i.e., a nitrate only addition, followed by a stoichiometrically ‘balanced’ 106Carbon:16Nitrogen:1Phosphorus addition. We quantified NO3-N uptake kinetics while constraining three variables: stream order, sediment type and type of injection (N vs stoichiometrically ‘balanced’ C:N:P). Following the laboratory nutrient uptake experiments, the columns were destructively sampled and the contents were homogenized to collect subsamples for DNA sequencing. Amplicon analysis was carried out as described by the Earth Microbiome Protocol for 16s and ITS sequencing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Sequence Design of Random Heteropolymers as Protein Mimics

Random heteropolymers (RHPs) have been computationally designed and experimentally shown to recapitulate protein-like phase behavior and function. However, unlike proteins, RHP sequences are only statistically defined and cannot be sequenced. Recent developments in reversible-deactivation radical polymerization allowed simulated polymer sequences based on the well-established Mayo–Lewis equation to more accurately reflect ground-truth sequences that are experimentally synthesized. This led to opportunities to perform bioinformatics-inspired analysis on simulated sequences to guide the design, synthesis, and interpretation of RHPs. We compared batches on the order of 10000 simulated RHP sequences that vary by synthetically controllable and measurable RHP characteristics such as chemical heterogeneity and average degree of polymerization. Our analysis spans across 3 levels: segments along a single chain, sequences within a batch, and batch-averaged statistics. We discuss simulator fidelity and highlight the importance of robust segment definition. Examples are presented that demonstrate the use of simulated sequence analysis for in-silico iterative design to mimic protein hydrophobic/hydrophilic segment distributions in RHPs and compare RHP and protein sequence segments to explain experimental results of RHPs that mimic protein function. To facilitate the community use of this workflow, the simulator and analysis modules have been made available through an open source toolkit, the RHPapp.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

cWINNOWER algorithm for finding fuzzy dna motifs

The cWINNOWER algorithm detects fuzzy motifs in DNA sequences rich in protein-binding signals. A signal is defined as any short nucleotide pattern having up to d mutations differing from a motif of length l. The algorithm finds such motifs if a clique consisting of a sufficiently large number of mutated copies of the motif (i.e., the signals) is present in the DNA sequence. The cWINNOWER algorithm substantially improves the sensitivity of the winnower method of Pevzner and Sze by imposing a consensus constraint, enabling it to detect much weaker signals. We studied the minimum detectable clique size qc as a function of sequence length N for random sequences. We found that qc increases linearly with N for a fast version of the algorithm based on counting three-member sub-cliques. Imposing consensus constraints reduces qc by a factor of three in this case, which makes the algorithm dramatically more sensitive. Our most sensitive algorithm, which counts four-member sub-cliques, needs a minimum of only 13 signals to detect motifs in a sequence of length N = 12,000 for (l, d) = (15, 4). Copyright Imperial College Press.

Evaluation Studies↗
DOE OSTI2022

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts—twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

transcription start and end sites↗
DOE OSTI2019

Genome sequence of the model rice variety KitaakeX

Abstract Background The availability of thousands of complete rice genome sequences from diverse varieties and accessions has laid the foundation for in-depth exploration of the rice genome. One drawback to these collections is that most of these rice varieties have long life cycles, and/or low transformation efficiencies, which limits their usefulness as model organisms for functional genomics studies. In contrast, the rice variety Kitaake has a rapid life cycle (9 weeks seed to seed) and is easy to transform and propagate. For these reasons, Kitaake has emerged as a model for studies of diverse monocotyledonous species. Results Here, we report the de novo genome sequencing and analysis of Oryza sativa ssp. japonica variety KitaakeX, a Kitaake plant carrying the rice XA21 immune receptor. Our KitaakeX sequence assembly contains 377.6 Mb, consisting of 33 scaffolds (476 contigs) with a contig N50 of 1.4 Mb. Complementing the assembly are detailed gene annotations of 35,594 protein coding genes. We identified 331,335 genomic variations between KitaakeX and Nipponbare (ssp. japonica ), and 2,785,991 variations between KitaakeX and Zhenshan97 (ssp. indica ). We also compared Kitaake resequencing reads to the KitaakeX assembly and identified 219 small variations. The high-quality genome of the model rice plant KitaakeX will accelerate rice functional genomics. Conclusions The high quality, de novo assembly of the KitaakeX genome will serve as a useful reference genome for rice and will accelerate functional genomics studies of rice and other species.

Nipponbare↗
NASA NTRS2004

Variation in sulfide tolerance of photosystem II in phylogenetically diverse cyanobacteria from sulfidic habitats

Physiological and molecular phylogenetic approaches were used to investigate variation among 12 cyanobacterial strains in their tolerance of sulfide, an inhibitor of oxygenic photosynthesis. Cyanobacteria from sulfidic habitats were found to be phylogenetically diverse and exhibited an approximately 50-fold variation in photosystem II performance in the presence of sulfide. Whereas the degree of tolerance was positively correlated with sulfide levels in the environment, a strain's phenotype could not be predicted from the tolerance of its closest relatives. These observations suggest that sulfide tolerance is a dynamic trait primarily shaped by environmental variation. Despite differences in absolute tolerance, similarities among strains in the effects of sulfide on chlorophyll fluorescence induction indicated a common mode of toxicity. Based on similarities with treatments known to disrupt the oxygen-evolving complex, it was concluded that sulfide toxicity resulted from inhibition of the donor side of photosystem II.

Phylogeny↗
NASA NTRS1995

Mutational analysis of photosystem I polypeptides in the cyanobacterium Synechocystis sp. PCC 6803. Targeted inactivation of psaI reveals the function of psaI in the structural organization of psaL

We cloned, characterized, and inactivated the psaI gene encoding a 4-kDa hydrophobic subunit of photosystem I from the cyanobacterium Synechocystis sp. PCC 6803. The psaI gene is located 90 base pairs downstream from psaL, and is transcribed on 0.94- and 0.32-kilobase transcripts. To identify the function of PsaI, we generated a cyanobacterial strain in which psaI has been interrupted by a gene for chloramphenicol resistance. The wild-type and the mutant cells showed comparable rates of photoautotrophic growth at 25 degrees C. However, the mutant cells grew slower and contained less chlorophyll than the wild-type cells, when grown at 40 degrees C. The PsaI-less membranes from cells grown at either temperature showed a small decrease in NADP+ photoreduction rate when compared to the wild-type membranes. Inactivation of psaI led to an 80% decrease in the PsaL level in the photosynthetic membranes and to a complete loss of PsaL in the purified photosystem I preparations, but had little effect on the accumulation of other photosystem I subunits. Upon solubilization with nonionic detergents, photosystem I trimers could be obtained from the wild-type, but not from the PsaI-less membranes. The PsaI-less photosystem I monomers did not contain detectable levels of PsaL. Therefore, a structural interaction between PsaL and PsaI may stabilize the association of PsaL with the photosystem I core. PsaL in the wild-type and PsaI-less membranes showed equal resistance to removal by chaotropic agents. However, PsaL in the PsaI-less strain exhibited an increased susceptibility to proteolysis. From these data, we conclude that PsaI has a crucial role in aiding normal structural organization of PsaL within the photosystem I complex and the absence of PsaI alters PsaL organization, leading to a small, but physiologically significant, defect in photosystem I function.

NASA Discipline Cell Biology↗
NASA NTRS2002

The Arabidopsis SKU5 gene encodes an extracellular glycosyl phosphatidylinositol-anchored glycoprotein involved in directional root growth

To investigate how roots respond to directional cues, we characterized a T-DNA-tagged Arabidopsis mutant named sku5 in which the roots skewed and looped away from the normal downward direction of growth on inclined agar surfaces. sku5 roots and etiolated hypocotyls were slightly shorter than normal and exhibited a counterclockwise (left-handed) axial rotation bias. The surface-dependent skewing phenotype disappeared when the roots penetrated the agar surface, but the axial rotation defect persisted, revealing that these two directional growth processes are separable. The SKU5 gene belongs to a 19-member gene family designated SKS (SKU5 Similar) that is related structurally to the multiple-copper oxidases ascorbate oxidase and laccase. However, the SKS proteins lack several of the conserved copper binding motifs characteristic of copper oxidases, and no enzymatic function could be assigned to the SKU5 protein. Analysis of plants expressing SKU5 reporter constructs and protein gel blot analysis showed that SKU5 was expressed most strongly in expanding tissues. SKU5 was glycosylated and modified by glycosyl phosphatidylinositol and localized to both the plasma membrane and the cell wall. Our observations suggest that SKU5 affects two directional growth processes, possibly by participating in cell wall expansion.

NASA Discipline Plant Biology↗
NASA NTRS1989

Mass analysis for the Space Station ECLSS using the balance spreadsheet method

The balance spreadsheet method is applied to mass analysis of the Environmental Control and Life Support System (ECLSS). The spreadsheet layout reduces the complexity of the ECLSS analysis by concisely defining the sources, sinks, and net changes in mass for each fluid. The analysis method is illustrated by using information from the latest Space Station ECLSS Architectural Control Documents and a given Space Station assembly sequence. The analysis results are plotted and discussed.

Chu, Wen-Ho↗
NASA NTRS2003

Implied alignment: a synapomorphy-based multiple-sequence alignment method and its use in cladogram search

A method to align sequence data based on parsimonious synapomorphy schemes generated by direct optimization (DO; earlier termed optimization alignment) is proposed. DO directly diagnoses sequence data on cladograms without an intervening multiple-alignment step, thereby creating topology-specific, dynamic homology statements. Hence, no multiple-alignment is required to generate cladograms. Unlike general and globally optimal multiple-alignment procedures, the method described here, implied alignment (IA), takes these dynamic homologies and traces them back through a single cladogram, linking the unaligned sequence positions in the terminal taxa via DO transformation series. These "lines of correspondence" link ancestor-descendent states and, when displayed as linearly arrayed columns without hypothetical ancestors, are largely indistinguishable from standard multiple alignment. Since this method is based on synapomorphy, the treatment of certain classes of insertion-deletion (indel) events may be different from that of other alignment procedures. As with all alignment methods, results are dependent on parameter assumptions such as indel cost and transversion:transition ratios. Such an IA could be used as a basis for phylogenetic search, but this would be questionable since the homologies derived from the implied alignment depend on its natal cladogram and any variance, between DO and IA + Search, due to heuristic approach. The utility of this procedure in heuristic cladogram searches using DO and the improvement of heuristic cladogram cost calculations are discussed. c2003 The Willi Hennig Society. Published by Elsevier Science (USA). All rights reserved.

Non-NASA Center↗
DOE OSTI2024

A General Materials Data Science Framework for Quantitative 2D Analysis of Particle Growth from Image Sequences

Abstract Phase transformations are a challenging problem in materials science, which lead to changes in properties and may impact performance of material systems in various applications. We introduce a general framework for the analysis of particle growth kinetics by utilizing concepts from machine learning and graph theory. As a model system, we use image sequences of atomic force microscopy showing the crystallization of an amorphous fluoroelastomer film. To identify crystalline particles in an amorphous matrix and track the temporal evolution of the particle dispersion, we have developed quantitative methods of 2D analysis. 700 image sequences were analyzed using a neural network architecture, achieving 0.97 pixel-wise classification accuracy as a measure of the correctly classified pixels. The growth kinetics of isolated and impinged particles were tracked throughout time using these image sequences. The relationship between image sequences and spatiotemporal graph representations was explored to identify the proximity of crystallites from each other. The framework enables the analysis of all image sequences without the requirement of sampling for specific particles or timesteps for various materials systems.

36 MATERIALS SCIENCE↗