1
|
Boffi NM, Guo Y, Rycroft CH, Amir A. How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution. eLife 2024; 12:RP87895. [PMID: 38376390 PMCID: PMC10942580 DOI: 10.7554/elife.87895] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024] Open
Abstract
The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population's mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn.
Collapse
Affiliation(s)
- Nicholas M Boffi
- Courant Institute of Mathematical Sciences, New York UniversityNew YorkUnited States
| | - Yipei Guo
- Janelia Research CampusAshburnUnited States
| | - Chris H Rycroft
- Department of Mathematics, University of Wisconsin–MadisonMadisonUnited States
- Mathematics Group, Lawrence Berkeley National LaboratoryBerkeleyUnited States
| | - Ariel Amir
- Weizmann Institute of ScienceRehovotIsrael
- John A. Paulson School of Engineering and Applied Sciences, Harvard UniversityCambridgeUnited States
| |
Collapse
|
2
|
Lin YL, Smith SN, Kanso E, Septer AN, Rycroft CH. A subcellular biochemical model for T6SS dynamics reveals winning competitive strategies. PNAS Nexus 2023; 2:pgad195. [PMID: 37441614 PMCID: PMC10335733 DOI: 10.1093/pnasnexus/pgad195] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/03/2022] [Revised: 05/30/2023] [Accepted: 06/01/2023] [Indexed: 07/15/2023]
Abstract
The type VI secretion system (T6SS) is a broadly distributed interbacterial weapon that can be used to eliminate competing bacterial populations. Although unarmed target populations are typically used to study T6SS function in vitro, bacteria most likely encounter other T6SS-armed competitors in nature. However, the connection between subcellular details of the T6SS and the outcomes of such mutually lethal battles is not well understood. Here, we incorporate biological data derived from natural competitors of Vibrio fischeri light organ symbionts to build a biochemical model for T6SS at the single-cell level, which we then integrate into an agent-based model (ABM). Using the ABM, we isolate and experiment with strain-specific physiological differences between competitors in ways not possible with biological samples to identify winning strategies for T6SS-armed populations. Through in vitro experiments, we discover that strain-specific differences exist in T6SS activation speed. ABM simulations corroborate that faster activation is dominant in determining survival during competition. Once competitors are fully activated, the energy required for T6SS creates a tipping point where increased weapon building and firing becomes too costly to be advantageous. Through ABM simulations, we identify the threshold where this transition occurs in the T6SS parameter space. We also find that competitive outcomes depend on the geometry of the battlefield: unarmed target cells survive at the edges of a range expansion where unlimited territory can be claimed. Alternatively, competitions within a confined space, much like the light organ crypts where natural V. fischeri compete, result in the rapid elimination of the unarmed population.
Collapse
Affiliation(s)
| | | | - Eva Kanso
- Department of Aerospace and Mechanical Engineering, University of Southern California, 3650 McClintock Ave, Los Angeles, CA 90089, USA
| | | | | |
Collapse
|
3
|
Gaucherand L, Iyer A, Gilabert I, Rycroft CH, Gaglia MM. Cut site preference allows influenza A virus PA-X to discriminate between host and viral mRNAs. Nat Microbiol 2023; 8:1304-1317. [PMID: 37349586 PMCID: PMC10690756 DOI: 10.1038/s41564-023-01409-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2022] [Accepted: 05/10/2023] [Indexed: 06/24/2023]
Abstract
Many viruses block host gene expression to take over the infected cell. This process, termed host shutoff, is thought to promote viral replication by preventing antiviral responses and redirecting cellular resources to viral processes. Several viruses from divergent families accomplish host shutoff through RNA degradation by endoribonucleases. However, viruses also need to ensure expression of their own genes. The influenza A virus endoribonuclease PA-X solves this problem by sparing viral mRNAs and some host RNAs necessary for viral replication. To understand how PA-X distinguishes between RNAs, we characterized PA-X cut sites transcriptome-wide using 5' rapid amplification of complementary DNA ends coupled to high-throughput sequencing. This analysis, along with RNA structure predictions and validation experiments using reporters, shows that PA-Xs from multiple influenza strains preferentially cleave RNAs at GCUG tetramers in hairpin loops. Importantly, GCUG tetramers are enriched in the human but not the influenza transcriptome. Moreover, optimal PA-X cut sites inserted in the influenza A virus genome are quickly selected against during viral replication in cells. This finding suggests that PA-X evolved these cleavage characteristics to preferentially target host over viral mRNAs in a manner reminiscent of cellular self versus non-self discrimination.
Collapse
Affiliation(s)
- Lea Gaucherand
- Program in Molecular Microbiology, Tufts University Graduate School of Biomedical Sciences, Boston, MA, USA
- Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, MA, USA
- Architecture et Réactivité de l'ARN, Institut de Biologie Moléculaire et Cellulaire du CNRS, Université de Strasbourg, Strasbourg, France
| | - Amrita Iyer
- Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, MA, USA
| | - Isabel Gilabert
- Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, MA, USA
- Faculty of Experimental Sciences, Universidad Francisco de Vitoria, Madrid, Spain
| | - Chris H Rycroft
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
- Department of Mathematics, University of Wisconsin-Madison, Madison, WI, USA
| | - Marta M Gaglia
- Program in Molecular Microbiology, Tufts University Graduate School of Biomedical Sciences, Boston, MA, USA.
- Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, MA, USA.
- Institute for Molecular Virology and Department of Medical Microbiology and Immunology, University of Wisconsin-Madison, Madison, WI, USA.
| |
Collapse
|
4
|
Leembruggen M, Andrejevic J, Kudrolli A, Rycroft CH. Computational model of twisted elastic ribbons. Phys Rev E 2023; 108:015003. [PMID: 37583198 DOI: 10.1103/physreve.108.015003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Accepted: 06/26/2023] [Indexed: 08/17/2023]
Abstract
We develop an irregular lattice mass-spring model to simulate and study the deformation modes of a thin elastic ribbon as a function of applied end-to-end twist and tension. Our simulations reproduce all reported experimentally observed modes, including transitions from helicoids to longitudinal wrinkles, creased helicoids and loops with self-contact, and transverse wrinkles to accordion self-folds. Our simulations also show that the twist angles at which the primary longitudinal and transverse wrinkles appear are well described by various analyses of the Föppl-von Kármán equations, but the characteristic wavelength of the longitudinal wrinkles has a more complex relationship to applied tension than previously estimated. The clamped edges are shown to suppress longitudinal wrinkling over a distance set by the applied tension and the ribbon width, but otherwise have no apparent effect on measured wavelength. Further, by analyzing the stress profile, we find that longitudinal wrinkling does not completely alleviate compression, but caps the magnitude of the compression. Nonetheless, the width over which wrinkles form is observed to be wider than the near-threshold analysis predictions: the width is more consistent with the predictions of far-from-threshold analysis. However, the end-to-end contraction of the ribbon as a function of twist is found to more closely follow the corresponding near-threshold prediction as tension in the ribbon is increased, in contrast to the expectations of far-from-threshold analysis. These results point to the need for further theoretical analysis of this rich thin elastic system, guided by our physically robust and intuitive simulation model.
Collapse
Affiliation(s)
- Madelyn Leembruggen
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Jovana Andrejevic
- Department of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Arshad Kudrolli
- Department of Physics, Clark University, Worcester, Massachusetts 01610, USA
| | - Chris H Rycroft
- Department of Mathematics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
- Computational Research Division, Lawrence Berkeley Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
5
|
Andrejevic N, Andrejevic J, Bernevig BA, Regnault N, Han F, Fabbris G, Nguyen T, Drucker NC, Rycroft CH, Li M. Machine-Learning Spectral Indicators of Topology. Adv Mater 2022; 34:e2204113. [PMID: 36193763 DOI: 10.1002/adma.202204113] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/06/2022] [Revised: 07/18/2022] [Indexed: 06/16/2023]
Abstract
Topological materials discovery has emerged as an important frontier in condensed matter physics. While theoretical classification frameworks have been used to identify thousands of candidate topological materials, experimental determination of materials' topology often poses significant technical challenges. X-ray absorption spectroscopy (XAS) is a widely used materials characterization technique sensitive to atoms' local symmetry and chemical bonding, which are intimately linked to band topology by the theory of topological quantum chemistry (TQC). Moreover, as a local structural probe, XAS is known to have high quantitative agreement between experiment and calculation, suggesting that insights from computational spectra can effectively inform experiments. In this work, computed X-ray absorption near-edge structure (XANES) spectra of more than 10 000 inorganic materials to train a neural network (NN) classifier that predicts topological class directly from XANES signatures, achieving F1 scores of 89% and 93% for topological and trivial classes, respectively is leveraged. Given the simplicity of the XAS setup and its compatibility with multimodal sample environments, the proposed machine-learning-augmented XAS topological indicator has the potential to discover broader categories of topological materials, such as non-cleavable compounds and amorphous materials, and may further inform field-driven phenomena in situ, such as magnetic field-driven topological phase transitions.
Collapse
Affiliation(s)
- Nina Andrejevic
- Center for Nanoscale Materials, Argonne National Laboratory, Lemont, IL, 60439, USA
- Quantum Measurement Group, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| | - Jovana Andrejevic
- Department of Physics, University of Pennsylvania, Philadelphia, PA, 19104, USA
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - B Andrei Bernevig
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
- Department of Physics, Princeton University, Princeton, NJ, 08544, USA
- Donostia International Physics Center, P. Manuel de Lardizabal 4, Donostia-San Sebastian, 20018, Spain
- IKERBASQUE, Basque Foundation for Science, Plaza Euskadi 5, Bilbao, 48009, Spain
| | - Nicolas Regnault
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
- Department of Physics, Princeton University, Princeton, NJ, 08544, USA
| | - Fei Han
- Center for Nanoscale Materials, Argonne National Laboratory, Lemont, IL, 60439, USA
- Quantum Measurement Group, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
- IKERBASQUE, Basque Foundation for Science, Plaza Euskadi 5, Bilbao, 48009, Spain
- Department of Nuclear Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| | - Gilberto Fabbris
- Department of Nuclear Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
- Advanced Photon Source, Argonne National Laboratory, Lemont, IL, 60439, USA
| | - Thanh Nguyen
- Center for Nanoscale Materials, Argonne National Laboratory, Lemont, IL, 60439, USA
- Quantum Measurement Group, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
- IKERBASQUE, Basque Foundation for Science, Plaza Euskadi 5, Bilbao, 48009, Spain
- Department of Nuclear Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| | - Nathan C Drucker
- Center for Nanoscale Materials, Argonne National Laboratory, Lemont, IL, 60439, USA
- Quantum Measurement Group, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
- Department of Physics, University of Pennsylvania, Philadelphia, PA, 19104, USA
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Chris H Rycroft
- Department of Physics, University of Pennsylvania, Philadelphia, PA, 19104, USA
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
- Advanced Photon Source, Argonne National Laboratory, Lemont, IL, 60439, USA
- Department of Mathematics, University of Wisconsin-Madison, Madison, WI, 53706, USA
- Computational Research Division, Lawrence Berkeley Laboratory, Berkeley, CA, 94720, USA
| | - Mingda Li
- Center for Nanoscale Materials, Argonne National Laboratory, Lemont, IL, 60439, USA
- Quantum Measurement Group, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
- IKERBASQUE, Basque Foundation for Science, Plaza Euskadi 5, Bilbao, 48009, Spain
- Department of Nuclear Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| |
Collapse
|
6
|
Deng B, Zareei A, Ding X, Weaver JC, Rycroft CH, Bertoldi K. Inverse Design of Mechanical Metamaterials with Target Nonlinear Response via a Neural Accelerated Evolution Strategy. Adv Mater 2022; 34:e2206238. [PMID: 36103610 DOI: 10.1002/adma.202206238] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/09/2022] [Indexed: 06/15/2023]
Abstract
Materials with target nonlinear mechanical response can support the design of innovative soft robots, wearable devices, footwear, and energy-absorbing systems, yet it is challenging to realize them. Here, mechanical metamaterials based on hinged quadrilaterals are used as a platform to realize target nonlinear mechanical responses. It is first shown that by changing the shape of the quadrilaterals, the amount of internal rotations induced by the applied compression can be tuned, and a wide range of mechanical responses is achieved. Next, a neural network is introduced that provides a computationally inexpensive relationship between the parameters describing the geometry and the corresponding stress-strain response. Finally, it is shown that by combining the neural network with an evolution strategy, one can efficiently identify geometries resulting in a wide range of target nonlinear mechanical responses and design optimized energy-absorbing systems, soft robots, and morphing structures.
Collapse
Affiliation(s)
- Bolei Deng
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
- Computer Science and Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
- Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| | - Ahmad Zareei
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Xiaoxiao Ding
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - James C Weaver
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Chris H Rycroft
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
- Computational Research Division, Lawrence Berkeley Laboratory, Berkeley, CA, 94720, USA
| | - Katia Bertoldi
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| |
Collapse
|
7
|
Andrejevic J, Lee LM, Rubinstein SM, Rycroft CH. Author Correction: A model for the fragmentation kinetics of crumpled thin sheets. Nat Commun 2021; 12:2362. [PMID: 33859200 PMCID: PMC8050261 DOI: 10.1038/s41467-021-22791-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-22791-z
Collapse
Affiliation(s)
- Jovana Andrejevic
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
| | - Lisa M Lee
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
| | - Shmuel M Rubinstein
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA.,The Racah Institute of Physics, The Hebrew University of Jerusalem, Jerusalem, 91904, Israel
| | - Chris H Rycroft
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA. .,Computational Research Division, Lawrence Berkeley Laboratory, Berkeley, CA, USA.
| |
Collapse
|
8
|
Abstract
As a confined thin sheet crumples, it spontaneously segments into flat facets delimited by a network of ridges. Despite the apparent disorder of this process, statistical properties of crumpled sheets exhibit striking reproducibility. Experiments have shown that the total crease length accrues logarithmically when repeatedly compacting and unfolding a sheet of paper. Here, we offer insight to this unexpected result by exploring the correspondence between crumpling and fragmentation processes. We identify a physical model for the evolution of facet area and ridge length distributions of crumpled sheets, and propose a mechanism for re-fragmentation driven by geometric frustration. This mechanism establishes a feedback loop in which the facet size distribution informs the subsequent rate of fragmentation under repeated confinement, thereby producing a new size distribution. We then demonstrate the capacity of this model to reproduce the characteristic logarithmic scaling of total crease length, thereby supplying a missing physical basis for the observed phenomenon.
Collapse
Affiliation(s)
- Jovana Andrejevic
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
| | - Lisa M Lee
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
| | - Shmuel M Rubinstein
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
- The Racah Institute of Physics, The Hebrew University of Jerusalem, Jerusalem, 91904, Israel
| | - Chris H Rycroft
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA.
- Computational Research Division, Lawrence Berkeley Laboratory, Berkeley, CA, USA.
| |
Collapse
|
9
|
Derr NJ, Fronk DC, Weber CA, Mahadevan A, Rycroft CH, Mahadevan L. Flow-Driven Branching in a Frangible Porous Medium. Phys Rev Lett 2020; 125:158002. [PMID: 33095596 DOI: 10.1103/physrevlett.125.158002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/21/2020] [Accepted: 08/18/2020] [Indexed: 06/11/2023]
Abstract
Channel formation and branching is widely seen in physical systems where movement of fluid through a porous structure causes the spatiotemporal evolution of the medium. We provide a simple theoretical framework that embodies this feedback mechanism in a multiphase model for flow through a frangible porous medium with a dynamic permeability. Numerical simulations of the model show the emergence of branched networks whose topology is determined by the geometry of external flow forcing. This allows us to delineate the conditions under which splitting and/or coalescing branched network formation is favored, with potential implications for both understanding and controlling branching in soft frangible media.
Collapse
Affiliation(s)
- Nicholas J Derr
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
| | - David C Fronk
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Christoph A Weber
- Max Planck Institute for the Physics of Complex Systems, Dresden 01187, Germany
| | - Amala Mahadevan
- Woods Hole Oceanographic Institution, Woods Hole, Massachusetts 02450, USA
| | - Chris H Rycroft
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - L Mahadevan
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts 02138, USA
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA
| |
Collapse
|
10
|
Boffi NM, Rycroft CH. Coordinate transformation methodology for simulating quasistatic elastoplastic solids. Phys Rev E 2020; 101:053304. [PMID: 32575210 DOI: 10.1103/physreve.101.053304] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2019] [Accepted: 02/03/2020] [Indexed: 11/07/2022]
Abstract
Molecular dynamics simulations frequently employ periodic boundary conditions where the positions of the periodic images are manipulated in order to apply deformation to the material sample. For example, Lees-Edwards conditions use moving periodic images to apply simple shear. Here, we examine the problem of precisely comparing this type of simulation to continuum solid mechanics. We employ a hypoelastoplastic mechanical model, and develop a projection method to enforce quasistatic equilibrium. We introduce a simulation framework that uses a fixed Cartesian computational grid on a reference domain, and which imposes deformation via a time-dependent coordinate transformation to the physical domain. As a test case for our method, we consider the evolution of shear bands in a bulk metallic glass using the shear transformation zone theory of amorphous plasticity. We examine the growth of shear bands in simple shear and pure shear conditions as a function of the initial preparation of the bulk metallic glass.
Collapse
Affiliation(s)
- Nicholas M Boffi
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Chris H Rycroft
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA.,Computational Research Division, Lawrence Berkeley Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
11
|
Gaucherand L, Porter BK, Levene RE, Price EL, Schmaling SK, Rycroft CH, Kevorkian Y, McCormick C, Khaperskyy DA, Gaglia MM. The Influenza A Virus Endoribonuclease PA-X Usurps Host mRNA Processing Machinery to Limit Host Gene Expression. Cell Rep 2020; 27:776-792.e7. [PMID: 30995476 PMCID: PMC6499400 DOI: 10.1016/j.celrep.2019.03.063] [Citation(s) in RCA: 58] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2018] [Revised: 02/08/2019] [Accepted: 03/15/2019] [Indexed: 12/29/2022] Open
Abstract
Many viruses shut off host gene expression to inhibit antiviral responses. Viral proteins and host proteins required for viral replication are typically spared in this process, but the mechanisms of target selectivity during host shutoff remain poorly understood. Using transcriptome-wide and targeted reporter experiments, we demonstrate that the influenza A virus endoribonuclease PA-X usurps RNA splicing to selectively target host RNAs for destruction. Proximity-labeling proteomics reveals that PA-X interacts with cellular RNA processing proteins, some of which are partially required for host shutoff. Thus, PA-X taps into host nuclear pre-mRNA processing mechanisms to destroy nascent mRNAs shortly after their synthesis. This mechanism sets PA-X apart from other viral host shutoff proteins that target actively translating mRNAs in the cytoplasm. Our study reveals a unique mechanism of host shutoff that helps us understand how influenza viruses suppress host gene expression. Influenza A virus PA-X targets the majority of host mRNAs for destruction Downregulation by PA-X correlates with the number of splice sites in a transcript Splicing renders RNAs susceptible to PA-X The cellular CFIm complex interacts with PA-X and contributes to PA-X activity
Collapse
Affiliation(s)
- Lea Gaucherand
- Graduate Program in Molecular Microbiology, Sackler School of Graduate Biomedical Sciences, Tufts University, Boston, MA, USA; Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, MA, USA
| | - Brittany K Porter
- Department of Microbiology & Immunology, Dalhousie University, Halifax, NS, Canada
| | - Rachel E Levene
- Graduate Program in Molecular Microbiology, Sackler School of Graduate Biomedical Sciences, Tufts University, Boston, MA, USA; Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, MA, USA
| | - Emma L Price
- Department of Microbiology & Immunology, Dalhousie University, Halifax, NS, Canada
| | - Summer K Schmaling
- Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, MA, USA
| | - Chris H Rycroft
- Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA; Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Yuzo Kevorkian
- Graduate Program in Molecular Microbiology, Sackler School of Graduate Biomedical Sciences, Tufts University, Boston, MA, USA; Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, MA, USA
| | - Craig McCormick
- Department of Microbiology & Immunology, Dalhousie University, Halifax, NS, Canada.
| | - Denys A Khaperskyy
- Department of Microbiology & Immunology, Dalhousie University, Halifax, NS, Canada.
| | - Marta M Gaglia
- Graduate Program in Molecular Microbiology, Sackler School of Graduate Biomedical Sciences, Tufts University, Boston, MA, USA; Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, MA, USA.
| |
Collapse
|
12
|
Choi GPT, Chiu B, Rycroft CH. Area-Preserving Mapping of 3D Carotid Ultrasound Images Using Density-Equalizing Reference Map. IEEE Trans Biomed Eng 2020; 67:2507-2517. [PMID: 31905128 DOI: 10.1109/tbme.2019.2963783] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Atherosclerotic plaques are focal and tend to occur at arterial bends and bifurcations. To quantitatively monitor the local changes in the carotid vessel-wall-plus-plaque thickness (VWT) and compare the VWT distributions for different patients or for the same patients at different ultrasound scanning sessions, a mapping technique is required to adjust for the geometric variability of different carotid artery models. In this work, we propose a novel method called density-equalizing reference map (DERM) for mapping 3D carotid surfaces to a standardized 2D carotid template, with an emphasis on preserving the local geometry of the carotid surface by minimizing the local area distortion. The initial map was generated by a previously described arc-length scaling (ALS) mapping method, which projects a 3D carotid surface onto a 2D non-convex L-shaped domain. A smooth and area-preserving flattened map was subsequently constructed by deforming the ALS map using the proposed algorithm that combines the density-equalizing map and the reference map techniques. This combination allows, for the first time, one-to-one mapping from a 3D surface to a standardized non-convex planar domain in an area-preserving manner. Evaluations using 20 carotid surface models show that the proposed method reduced the area distortion of the flattening maps by over 80% as compared to the ALS mapping method. The proposed method is capable of improving the accuracy of area estimation for plaque regions without compromising inter-scan reproducibility.
Collapse
|
13
|
Tang Y, Zhang L, Guo Q, Xia B, Yin Z, Cao J, Tong J, Rycroft CH. Analysis of the pebble burnup profile in a pebble-bed nuclear reactor. Nuclear Engineering and Design 2019. [DOI: 10.1016/j.nucengdes.2019.01.030] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
14
|
Hoffmann J, Bar-Sinai Y, Lee LM, Andrejevic J, Mishra S, Rubinstein SM, Rycroft CH. Machine learning in a data-limited regime: Augmenting experiments with synthetic data uncovers order in crumpled sheets. Sci Adv 2019; 5:eaau6792. [PMID: 31032399 PMCID: PMC6486215 DOI: 10.1126/sciadv.aau6792] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Accepted: 03/06/2019] [Indexed: 05/09/2023]
Abstract
Machine learning has gained widespread attention as a powerful tool to identify structure in complex, high-dimensional data. However, these techniques are ostensibly inapplicable for experimental systems where data are scarce or expensive to obtain. Here, we introduce a strategy to resolve this impasse by augmenting the experimental dataset with synthetically generated data of a much simpler sister system. Specifically, we study spontaneously emerging local order in crease networks of crumpled thin sheets, a paradigmatic example of spatial complexity, and show that machine learning techniques can be effective even in a data-limited regime. This is achieved by augmenting the scarce experimental dataset with inexhaustible amounts of simulated data of rigid flat-folded sheets, which are simple to simulate and share common statistical properties. This considerably improves the predictive power in a test problem of pattern completion and demonstrates the usefulness of machine learning in bench-top experiments where data are good but scarce.
Collapse
Affiliation(s)
- Jordan Hoffmann
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
| | - Yohai Bar-Sinai
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
- Corresponding author. (Y.B.-S.); (S.M.R.)
| | - Lisa M. Lee
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
| | - Jovana Andrejevic
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
| | - Shruti Mishra
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
| | - Shmuel M. Rubinstein
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
- Corresponding author. (Y.B.-S.); (S.M.R.)
| | - Chris H. Rycroft
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
- Computational Research Division, Lawrence Berkeley Laboratory, Berkeley, CA 94720, USA
| |
Collapse
|
15
|
Hoffmann J, Donoughe S, Li K, Salcedo MK, Rycroft CH. A simple developmental model recapitulates complex insect wing venation patterns. Proc Natl Acad Sci U S A 2018; 115:9905-9910. [PMID: 30224459 PMCID: PMC6176563 DOI: 10.1073/pnas.1721248115] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
Insect wings are typically supported by thickened struts called veins. These veins form diverse geometric patterns across insects. For many insect species, even the left and right wings from the same individual have veins with unique topological arrangements, and little is known about how these patterns form. We present a large-scale quantitative study of the fingerprint-like "secondary veins." We compile a dataset of wings from 232 species and 17 families from the order Odonata (dragonflies and damselflies), a group with particularly elaborate vein patterns. We characterize the geometric arrangements of veins and develop a simple model of secondary vein patterning. We show that our model is capable of recapitulating the vein geometries of species from other, distantly related winged insect clades.
Collapse
Affiliation(s)
- Jordan Hoffmann
- Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138
| | - Seth Donoughe
- Department of Molecular Genetics and Cell Biology, University of Chicago, Chicago, IL 60637;
| | - Kathy Li
- Applied Physics and Applied Mathematics Department, Columbia University, New York, NY 10027
| | - Mary K Salcedo
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138
| | - Chris H Rycroft
- Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138;
- Computational Research Division, Lawrence Berkeley Laboratory, Berkeley, CA 94720
| |
Collapse
|
16
|
Weber CA, Rycroft CH, Mahadevan L. Differential Activity-Driven Instabilities in Biphasic Active Matter. Phys Rev Lett 2018; 120:248003. [PMID: 29956948 DOI: 10.1103/physrevlett.120.248003] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2017] [Revised: 02/25/2018] [Indexed: 06/08/2023]
Abstract
Active stresses can cause instabilities in contractile gels and living tissues. Here we provide a generic hydrodynamic theory that treats these systems as a mixture of two phases of varying activity and different mechanical properties. We find that differential activity between the phases causes a uniform mixture to undergo a demixing instability. We follow the nonlinear evolution of the instability and characterize a phase diagram of the resulting patterns. Our study complements other instability mechanisms in mixtures driven by differential adhesion, differential diffusion, differential growth, and differential motion.
Collapse
Affiliation(s)
- Christoph A Weber
- Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Chris H Rycroft
- Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
- Mathematics Group, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - L Mahadevan
- Paulson School of Engineering and Applied Sciences, Department of Physics, Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| |
Collapse
|
17
|
Hinkle AR, Rycroft CH, Shields MD, Falk ML. Coarse graining atomistic simulations of plastically deforming amorphous solids. Phys Rev E 2017; 95:053001. [PMID: 28618619 DOI: 10.1103/physreve.95.053001] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2017] [Indexed: 06/07/2023]
Abstract
The primary mode of failure in disordered solids results from the formation and persistence of highly localized regions of large plastic strains known as shear bands. Continuum-level field theories capable of predicting this mechanical response rely upon an accurate representation of the initial and evolving states of the amorphous structure. We perform molecular dynamics simulations of a metallic glass and propose a methodology for coarse graining discrete, atomistic quantities, such as the potential energies of the elemental constituents. A strain criterion is established and used to distinguish the coarse-grained degrees-of-freedom inside the emerging shear band from those of the surrounding material. A signal-to-noise ratio provides a means of evaluating the strength of the signal of the shear band as a function of the coarse graining. Finally, we investigate the effect of different coarse graining length scales by comparing a two-dimensional, numerical implementation of the effective-temperature description in the shear transformation zone (STZ) theory with direct molecular dynamics simulations. These comparisons indicate the coarse graining length scale has a lower bound, above which there is a high level of agreement between the atomistics and the STZ theory, and below which the concept of effective temperature breaks down.
Collapse
Affiliation(s)
- Adam R Hinkle
- Department of Materials Science & Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Chris H Rycroft
- Paulson School of Engineering & Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Michael D Shields
- Department of Materials Science & Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Department of Civil Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Michael L Falk
- Department of Materials Science & Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Department of Physics & Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
18
|
Abstract
An advection-diffusion-limited dissolution model of an object being eroded by a two-dimensional potential flow is presented. By taking advantage of the conformal invariance of the model, a numerical method is introduced that tracks the evolution of the object boundary in terms of a time-dependent Laurent series. Simulations of a variety of dissolving objects are shown, which shrink and collapse to a single point in finite time. The simulations reveal a surprising exact relationship, whereby the collapse point is the root of a non-analytic function given in terms of the flow velocity and the Laurent series coefficients describing the initial shape. This result is subsequently derived using residue calculus. The structure of the non-analytic function is examined for three different test cases, and a practical approach to determine the collapse point using a generalized Newton-Raphson root-finding algorithm is outlined. These examples also illustrate the possibility that the model breaks down in finite time prior to complete collapse, due to a topological singularity, as the dissolving boundary overlaps itself rather than breaking up into multiple domains (analogous to droplet pinch-off in fluid mechanics). The model raises fundamental mathematical questions about broken symmetries in finite-time singularities of both continuous and stochastic dynamical systems.
Collapse
Affiliation(s)
- Chris H Rycroft
- Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA; Department of Mathematics, Lawrence Berkeley Laboratory, Berkeley, CA 94720, USA
| | - Martin Z Bazant
- Departments of Chemical Engineering and Mathematics , Massachusetts Institute of Technology , Cambridge, MA 02139, USA
| |
Collapse
|
19
|
Gaglia MM, Rycroft CH, Glaunsinger BA. Transcriptome-Wide Cleavage Site Mapping on Cellular mRNAs Reveals Features Underlying Sequence-Specific Cleavage by the Viral Ribonuclease SOX. PLoS Pathog 2015; 11:e1005305. [PMID: 26646420 PMCID: PMC4672902 DOI: 10.1371/journal.ppat.1005305] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2015] [Accepted: 11/03/2015] [Indexed: 11/18/2022] Open
Abstract
Many viruses express factors that reduce host gene expression through widespread degradation of cellular mRNA. An example of this class of proteins is the mRNA-targeting endoribonuclease SOX from the gamma-herpesvirus Kaposi's sarcoma-associated herpesvirus (KSHV). Previous studies indicated that cleavage of messenger RNAs (mRNA) by SOX occurs at specific locations defined by the sequence of the target RNA, which is at odds with the down-regulation of a large portion of cellular transcripts. In this study, we address this paradox by using high-throughput sequencing of cleavage intermediates combined with a custom bioinformatics-based analysis pipeline to identify SOX cleavage sites across the mRNA transcriptome. These data, coupled with targeted mutagenesis, reveal that while cleavage sites are specific and reproducible, they are defined by a degenerate sequence motif containing a small number of conserved residues rather than a strong consensus sequence. This degenerate element is well represented in both human and KSHV mRNA, and its presence correlates with RNA destabilization by SOX. This represents a new endonuclease targeting strategy, in which use of a degenerate targeting element enables RNA cleavage at specific locations without restricting the range of targets. Furthermore, it shows that strong target selectivity can be achieved without a high degree of sequence specificity.
Collapse
Affiliation(s)
- Marta Maria Gaglia
- Program in Molecular Microbiology, Sackler School of Graduate Biomedical Sciences, Tufts University, Boston, Massachusetts, United States of America
- Department of Molecular Biology and Microbiology, Tufts University School of Medicine, Boston, Massachusetts, United States of America
- * E-mail: (MMG); (BAG)
| | - Chris H. Rycroft
- School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts, United States of America
- Department of Mathematics, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America
| | - Britt A. Glaunsinger
- Howard Hughes Medical Institute, University of California, Berkeley, Berkeley, California, United States of America
- Department of Plant and Microbial Biology, University of California, Berkeley, Berkeley, California, United States of America
- * E-mail: (MMG); (BAG)
| |
Collapse
|
20
|
Rycroft CH, Dehbi A, Lind T, Güntay S. Granular flow in pebble-bed nuclear reactors: Scaling, dust generation, and stress. Nuclear Engineering and Design 2013. [DOI: 10.1016/j.nucengdes.2013.07.010] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
21
|
|
22
|
Abstract
Quantitative understanding of the fracture toughness of metallic glasses, including the associated ductile-to-brittle (embrittlement) transitions, is not yet available. Here, we use a simple model of plastic deformation in glasses, coupled to an advanced Eulerian level set formulation for solving complex free-boundary problems, to calculate the fracture toughness of metallic glasses as a function of the degree of structural relaxation corresponding to different annealing times near the glass temperature. Our main result indicates the existence of an elastoplastic crack tip instability for sufficiently relaxed glasses, resulting in a marked drop in the toughness, which we interpret as annealing-induced embrittlement transition similar to experimental observations.
Collapse
Affiliation(s)
- Chris H Rycroft
- Department of Mathematics, University of California, Berkeley, California 94720, USA
| | | |
Collapse
|
23
|
Lin LC, Berger AH, Martin RL, Kim J, Swisher JA, Jariwala K, Rycroft CH, Bhown AS, Deem MW, Haranczyk M, Smit B. In silico screening of carbon-capture materials. Nat Mater 2012; 11:633-41. [PMID: 22635045 DOI: 10.1038/nmat3336] [Citation(s) in RCA: 267] [Impact Index Per Article: 22.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/11/2011] [Accepted: 04/18/2012] [Indexed: 05/26/2023]
Abstract
One of the main bottlenecks to deploying large-scale carbon dioxide capture and storage (CCS) in power plants is the energy required to separate the CO(2) from flue gas. For example, near-term CCS technology applied to coal-fired power plants is projected to reduce the net output of the plant by some 30% and to increase the cost of electricity by 60-80%. Developing capture materials and processes that reduce the parasitic energy imposed by CCS is therefore an important area of research. We have developed a computational approach to rank adsorbents for their performance in CCS. Using this analysis, we have screened hundreds of thousands of zeolite and zeolitic imidazolate framework structures and identified many different structures that have the potential to reduce the parasitic energy of CCS by 30-40% compared with near-term technologies.
Collapse
Affiliation(s)
- Li-Chiang Lin
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720-1462, USA
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
24
|
|
25
|
Affiliation(s)
- Chris H Rycroft
- Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA
| |
Collapse
|
26
|
Rycroft CH, Orpe AV, Kudrolli A. Physical test of a particle simulation model in a sheared granular system. Phys Rev E Stat Nonlin Soft Matter Phys 2009; 80:031305. [PMID: 19905108 DOI: 10.1103/physreve.80.031305] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/21/2009] [Revised: 06/16/2009] [Indexed: 05/28/2023]
Abstract
We report a detailed comparison of a slow gravity-driven sheared granular flow with a discrete-element simulation performed in the same geometry. In the experiments, grains flow inside a silo with a rectangular cross section and are sheared by a rough boundary on one side and smooth boundaries on the other sides. Individual grain position and motion are measured using a particle index-matching imaging technique where a fluorescent dye is added to the interstitial liquid which has the same refractive index as the glass beads. The simulations use a Cundall-Strack contact model between the grains using contact parameters that have been used in many other previous studies and ignore the hydrodynamic effects of the interstitial liquid. Computations are performed to understand the effect of particle coefficient of friction, elasticity, contact model, and polydispersity on mean flow properties. We then perform a detailed comparison of the particle fluctuation properties as measured by the displacement probability distribution function and the mean square displacement. All in all, our study suggests a high level of quantitative agreement between the simulations and experiments.
Collapse
Affiliation(s)
- Chris H Rycroft
- Department of Mathematics, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA.
| | | | | |
Collapse
|
27
|
Rycroft CH, Grest GS, Landry JW, Bazant MZ. Analysis of granular flow in a pebble-bed nuclear reactor. Phys Rev E Stat Nonlin Soft Matter Phys 2006; 74:021306. [PMID: 17025418 DOI: 10.1103/physreve.74.021306] [Citation(s) in RCA: 95] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2006] [Indexed: 05/12/2023]
Abstract
Pebble-bed nuclear reactor technology, which is currently being revived around the world, raises fundamental questions about dense granular flow in silos. A typical reactor core is composed of graphite fuel pebbles, which drain very slowly in a continuous refueling process. Pebble flow is poorly understood and not easily accessible to experiments, and yet it has a major impact on reactor physics. To address this problem, we perform full-scale, discrete-element simulations in realistic geometries, with up to 440,000 frictional, viscoelastic 6-cm-diam spheres draining in a cylindrical vessel of diameter 3.5m and height 10 m with bottom funnels angled at 30 degrees or 60 degrees. We also simulate a bidisperse core with a dynamic central column of smaller graphite moderator pebbles and show that little mixing occurs down to a 1:2 diameter ratio. We analyze the mean velocity, diffusion and mixing, local ordering and porosity (from Voronoi volumes), the residence-time distribution, and the effects of wall friction and discuss implications for reactor design and the basic physics of granular flow.
Collapse
Affiliation(s)
- Chris H Rycroft
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | | | | | | |
Collapse
|
28
|
Rycroft CH, Bazant MZ, Grest GS, Landry JW. Dynamics of random packings in granular flow. Phys Rev E Stat Nonlin Soft Matter Phys 2006; 73:051306. [PMID: 16802932 DOI: 10.1103/physreve.73.051306] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2005] [Indexed: 05/10/2023]
Abstract
We present a multiscale simulation algorithm for amorphous materials, which we illustrate and validate in a canonical case of dense granular flow. Our algorithm is based on the recently proposed spot model, where particles in a dense random packing undergo chainlike collective displacements in response to diffusing "spots" of influence, carrying a slight excess of interstitial free volume. We reconstruct the microscopic dynamics of particles from the "coarse grained" dynamics of spots by introducing a localized particle relaxation step after each spot-induced block displacement, simply to enforce packing constraints with a (fairly arbitrary) soft-core repulsion. To test the model, we study to what extent it can describe the dynamics of up to 135,000 frictional, viscoelastic spheres in granular drainage simulated by the discrete-element method (DEM). With only five fitting parameters (the radius, volume, diffusivity, drift velocity, and injection rate of spots), we find that the spot simulations are able to largely reproduce not only the mean flow and diffusion, but also some subtle statistics of the flowing packings, such as spatial velocity correlations and many-body structural correlations. The spot simulations run over 100 times faster than the DEM and demonstrate the possibility of multiscale modeling for amorphous materials, whenever a suitable model can be devised for the coarse-grained spot dynamics.
Collapse
Affiliation(s)
- Chris H Rycroft
- Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA.
| | | | | | | |
Collapse
|