Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Glass dynamics”

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 199 records · Page 11

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl 2 from T m + 85 K to T m + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl 2 . We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. Furthermore, this result may contribute to new understandings of the structural relaxation of other network liquids.

36 MATERIALS SCIENCE↗

Memory in aging colloidal gels with time-varying attraction

Here, we report a combined rheology, x-ray photon correlation spectroscopy, and modeling study of gel formation and aging in suspensions of nanocolloidal spheres with volume fractions of 0.20 and 0.43 and with a short-range attraction whose strength is tuned by changing temperature. Following a quench from high temperature, where the colloids are essentially hard spheres, to a temperature below the gel point, the suspensions form gels that undergo aging characterized by a steadily increasing elastic shear modulus and slowing, increasingly constrained microscopic dynamics. The aging proceeds at a faster rate for stronger attraction strength. When the attraction strength is suddenly lowered during aging, the gel properties evolve non-monotonically in a manner resembling the Kovacs effect in glasses, in which the modulus decreases and the microscopic dynamics become less constrained for a period before more conventional aging resumes. Eventually, the properties of the gel following the decrease in attraction strength converge to those of a gel that has undergone aging at the lower attraction strength throughout. The time scale of this convergence increases as a power law with the age at which the attraction strength is decreased and decreases exponentially with the magnitude of the change in attraction. A model for gel aging in which particles attach and detach from the gel at rates that depend on their contact number reproduces these trends and reveals that the non-monotonic behavior results from the dispersion in the rates that the populations of particles with different contact number adjust to the new attraction strength.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-pressure annealing driven nanocrystal formation in Zr 50 Cu 40 Al 10 metallic glass and strength increase

Pressure-induced structural changes in metallic glasses have been of great interest as they are expected to open new ways to synthesize novel materials with unexpected properties. Here, we investigated the effect of simultaneous high-pressure and high-temperature treatment on the structure and properties of a Zr 50 Cu 40 Al 10 metallic glass by in situ X-ray structure measurement and property analysis of the final material. We found the unusual formation of Cu-rich nanocrystals at high pressure and temperature, accompanied by significant strength and hardness enhancement. Based on reverse Monte Carlo modeling and molecular dynamics simulations, the structure of the metallic glass changed to a densely packed, chemically uniform configuration with high short-range and medium-range ordering at high pressure and temperature. These results show that high-pressure annealing processes provide a new way to improve and control properties without changing their composition.

36 MATERIALS SCIENCE↗

Numerical modeling of Joule heated ceramic melter

Computational fluid dynamics and heat transfer models of waste glass melters are being developed to support the vitrification campaign at the Waste Treatment and Immobilization Plant (WTP) as part of the United States nuclear waste stewardship and environmental protection programs. The WTP mission aims to achieve the stable and durable processing of vitrified waste in the safest, most efficient, and economical manner possible. Hanford’s tank waste is far more complex and varied than the waste treated by other vitrification facilities and both startup and long-term operational challenges are likely to be even greater than prior experience has shown. There are thousands of different waste compositions, which makes processing the waste a complex and challenging task. Models of test melters, pilot-scale, and full-scale melters have been developed to provide understanding and insight into the various physical and chemical phenomena occurring during vitrification. So the results of these modeling activities will be used to support the WTP operations to process legacy nuclear waste into a stable form for disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Short- and medium-range orders in Al90Tb10 glass and their relation to the structures of competing crystalline phases

Molecular dynamics simulations using an interatomic potential developed by artificial neural network deep machine learning are performed to study the local structural order in Al90Tb10 metallic glass. We show that more than 80% of the Tb-centered clusters in Al90Tb10 glass have short-range order (SRO) with their 17 first coordination shell atoms stacked in a ‘3661’ or ‘15551’ sequence. Medium-range order (MRO) in Bergman-type packing extended out to the second and third coordination shells is also clearly observed. Analysis of the network formed by the ‘3661’ and ‘15551’ clusters show that ~82% of such SRO units share their faces or vertexes, while only ~6% of neighboring SRO pairs are interpenetrating. Such a network topology is consistent with the Bergman-type MRO around the Tb-centers. Moreover, crystal structure searches using genetic algorithm and the neural network interatomic potential reveal several low-energy metastable crystalline structures in the composition range close to Al 90 Tb 10 . Some of these crystalline structures have the ‘3661’ SRO while others have the ‘15551’ SRO. While the crystalline structures with the ‘3661’ SRO also exhibit the MRO very similar to that observed in the glass, the ones with the ‘15551’ SRO have very different atomic packing in the second and third shells around the Tb centers from that of the Bergman-type MRO observed in the glassy phase.

36 MATERIALS SCIENCE↗

Dynamics and Structure of Molten Neodymium Titanate

Binary rare earth titanates can be vitrified using levitation melt quenching, and their glasses are promising for optical applications. To better understand their thermophysical properties and dynamics, we have used molecular dynamics (MD) and empirical potential structure refinement (EPSR) to simulate molten neodymium titanate, 83TiO 2 –17Nd 2 O 3 . The models are compared with high-energy synchrotron X-ray diffraction of the melt structure over 1060–1960 K and thermophysical property measurements of electrostatically levitated droplets over 1500–2200 K. As temperature increases, the network periodicity increases and the correlation length decreases, both in monotonic trends that suggest continuous structural evolution through the melt and supercooled liquid regimes. At 1900 K, the mean Ti–O and Nd–O coordination numbers are 5.8(5) and 6.7(5), compared to the glass, 5.72(6) and 7.70(26). The equilibrium melting point is measured as 1860(20) K, over 100 K higher than previous estimates. Viscosity data from MD are in good agreement with experimental measurements. The dynamic propensity of Nd is similar to that of the network-forming Ti, suggesting that Nd does not participate as a conventional modifier. The viscosity of this fragile liquid is likely controlled by the breaking of cation–oxygen bonds and rearrangements of Ti–O polyhedra.

molecular dynamics↗

Microscopic Theory of Long-Time Center-of-Mass Self-Diffusion and Anomalous Transport in Ring Polymer Liquids

We construct a microscopic theory at the level of segment-scale correlated space–time intermolecular forces for the long-time center-of-mass (CM) diffusion constant and intermediate-time non-Fickian transport in dense solutions and melts of ring polymers. The approach combines ideas of polymer, colloid, and liquid-state statistical mechanics to quantify how the multifractal intra-ring conformational structure and inter-ring packing correlations determine dynamic caging constraints and time-dependent friction. Breakdown of Rouse theory is predicted to occur due to length scale-dependent temporal correlation of forces exerted on pairs of tagged ring segments from surrounding polymers. At large enough degrees of polymerization (N), a stronger scaling of the CM diffusion constant (D ∝ N –2 ) is predicted, with a crossover N D proportional to the product of the system-specific macromolecular volume fraction and dimensionless compressibility. In analogy with caging effects in glass-forming fluids, the theory does appear to begin to fail at sufficiently high N/N D for the center-of-mass diffusivity, likely due to another crossover to an even slower activated transport regime. However, use of N/N D with the theoretically predicted N D , in conjunction with dynamic blob scaling ideas for local Rouse friction, collapses semidilute and concentrated solution simulation data onto a master curve. Based on the same physical ideas employed to predict the diffusion constant, a generalized Langevin equation description is formulated for intermediate-time CM transport. It predicts two subdiffusive regimes that emerge due to the self-similar nature of internal ring structure. Analytic and numerical predictions for the apparent non-Fickian exponent as a function of time, N/N D and dimensionless compressibility, the evolution of the maximum degree of subdiffusive motion with N/N D , and the time scale for recovering Fickian diffusion are made, all of which are in good accordance with melt simulations. Furthermore, the present work sets the stage to address activated dynamics and glass formation on the macromolecular scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic and Kinetic Transitions of Liquids in Nanoconfinement

Core principles of chemistry are ubiquitously invoked to shed light on the nature of molecular level interactions in nanoconfined fluids, which play a pivotal role in a wide range of processes in geochemistry, biology, and engineering. A detailed understanding of the physicochemical processes involved in the flow, structural transitions, and freezing or melting behavior of fluids confined within nanometer-sized pores of solid materials is thus of enormous importance for both basic research and technological applications. This Account provides a perspective on new insights into the thermodynamic and kinetic transitions of nanoconfined fluids in their stable and metastable forms. After briefly introducing the unique properties of mesoporous silicas from the SBA, MCM, and FDU families that serve as the confinement matrices, combining highly ordered single and bimodal mesopore architectures with tunable pore sizes in the ~2-15 nm range and narrow size distributions, recent studies on melting/freezing behavior of water confined in these host matrices are reviewed. While differential scanning calorimetry (DSC) reveals a linear relationship between melting point depression and pore size (independent of the pore shape), as predicted by the Gibbs-Thomson relation, variable temperature 2 H wide-line nuclear magnetic resonance (NMR) spectroscopy studies confirm the core-shell model of water and give evidence for a layer-by-layer freezing mechanism, which gives rise to an apparent fragile-to-strong transition in the solidification dynamics. In contrast to the freezing/melting behavior of water, the effect of nanoconfinement on the glass transition of supercooled liquids is nonuniversal and the glass transition temperature T g can either increase or decrease with the dimensionality and extent of confinement. This nonuniversal behavior is exemplified by the two glass-forming molecular liquids, glycerol and ortho-terphenyl (OTP). While glycerol shows an increase in T g and a pronounced slowdown of the rotational dynamics of the constituent molecules due to a change in the molecular packing between the bulk and the confined liquid, OTP displays a linear and confining-media-dependent depression of T g with increased confinement that is strongly influenced by the pore-liquid interface characteristics. This Account concludes with a focus on recent experimental evidence of extreme spatial and dynamical heterogeneity in both freezing and glass transition processes. This discovery was enabled by the unique mesoporous structures of SBA-16 and FDU-5, possessing bimodal architectures with two interconnected pore types of different size and shape (spherical and cylindrical). For the very first time, two melting points for water and two glass transitions for supercooled OTP, corresponding to a specific pore type, were observed. Collectively, these observations strongly suggest a close mechanistic connection between the local fluctuations in the structure and dynamics of nanoconfined liquids. Finally, while the findings reviewed in this Account provide new insights into thermodynamic and kinetic transitions of fluids, there remain many unanswered questions regarding the effects of nanoconfinement on the fundamental properties of fluids, which offer exciting future opportunities in chemical research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-speed synchrotron X-ray phase-contrast imaging for evaluating microscale damage mechanisms and tracking cracking behaviors inside cross-ply GFRCs

In this study, we integrated the high-speed synchrotron X-ray phase-contrast imaging (PCI) with a modified Kolsky compression bar loading platform to visualize the dynamic failure processes of cross-ply glass fiber reinforced composites (GFRCs). Four S-2 glass/SC-15 composite specimens were prepared, having similar thicknesses but different stacking sequences, namely as [0 12 /90 12 ], [90 12 /0 12 ], [0 8 /90 8 /0 8 ], and [90 8 /0 8 /90 8 ]. Three-dimensional synchrotron X-ray computed tomography and scanning electron microscopy (SEM) were employed to examine the microstructures and quantify the fiber volume fractions. Each specimen was notched and subjected to a dynamic three-point flexural loading. The onset of cracking close to the notch tip, crack propagation in 0° or 90° plies and their interface, crack opening, and ultimately failure of the specimen were captured by high-speed synchrotron X-ray PCI. Additional dynamic experiments were performed to determine the average time when the stress wave propagated through the specimen and correlate the X-ray images with the specimen's force-deflection response. Finally, the surface morphology of each specimen after the dynamic loading was imaged by SEM. Comparison between real-time X-ray images and post-fracture SEM images demonstrated the capability of the X-ray method to record damage evolution inside composites. Furthermore, the high spatial and temporal resolutions of the X-ray setup and edge enhancement by PCI enabled the identification of microscale damage features within 1 μs. Two damaging processes were identified, crack growth was quantified, and fracture toughness of the composites was evaluated. The method is deemed useful to reveal microscale damage mechanisms and track cracking behaviors inside cross-ply GFRCs under dynamic loading in real time.

42 ENGINEERING↗

Rapid CFD Using Machine Learning Algorithms (CRADA Final Report)

This is a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (“LLNL”) and Guardian Glass, LLC ("Guardian Glass") to develop a fast-running emulator of the reactive Computational Fluid Dynamics (“CFD”) simulations needed to understand the complex reactions and flows in the glass melting, fining, and forming subprocesses. This CRADA project is sponsored under the High-Performance Computing for Manufacturing (“HPC4Mfg”) Program of the Department of Energy’s Advanced Manufacturing Office (“AMO”) within the Energy Efficiency and Renewable Energy (“EERE”) Office.

36 MATERIALS SCIENCE↗

Linking the Dynamic Response of the Carbon Dioxide-Concentrating Mechanism to Carbon Assimilation Behavior in Fremyella diplosiphon

ABSTRACT Cyanobacteria use a carbon dioxide (CO 2 )-concentrating mechanism (CCM) that enhances their carbon fixation efficiency and is regulated by many environmental factors that impact photosynthesis, including carbon availability, light levels, and nutrient access. Efforts to connect the regulation of the CCM by these factors to functional effects on carbon assimilation rates have been complicated by the aqueous nature of cyanobacteria. Here, we describe the use of cyanobacteria in a semiwet state on glass fiber filtration discs—cyanobacterial discs—to establish dynamic carbon assimilation behavior using gas exchange analysis. In combination with quantitative PCR (qPCR) and transmission electron microscopy (TEM) analyses, we linked the regulation of CCM components to corresponding carbon assimilation behavior in the freshwater, filamentous cyanobacterium Fremyella diplosiphon . Inorganic carbon (C i ) levels, light quantity, and light quality have all been shown to influence carbon assimilation behavior in F. diplosiphon . Our results suggest a biphasic model of cyanobacterial carbon fixation. While behavior at low levels of CO 2 is driven mainly by the C i uptake ability of the cyanobacterium, at higher CO 2 levels, carbon assimilation behavior is multifaceted and depends on C i availability, carboxysome morphology, linear electron flow, and cell shape. Carbon response curves (CRCs) generated via gas exchange analysis enable rapid examination of CO 2 assimilation behavior in cyanobacteria and can be used for cells grown under distinct conditions to provide insight into how CO 2 assimilation correlates with the regulation of critical cellular functions, such as the environmental control of the CCM and downstream photosynthetic capacity. IMPORTANCE Environmental regulation of photosynthesis in cyanobacteria enhances organismal fitness, light capture, and associated carbon fixation under dynamic conditions. Concentration of carbon dioxide (CO 2 ) near the carbon-fixing enzyme RubisCO occurs via the CO 2 -concentrating mechanism (CCM). The CCM is also tuned in response to carbon availability, light quality or levels, or nutrient access—cues that also impact photosynthesis. We adapted dynamic gas exchange methods generally used with plants to investigate environmental regulation of the CCM and carbon fixation capacity using glass fiber-filtered cells of the cyanobacterium Fremyella diplosiphon . We describe a breakthrough in measuring real-time carbon uptake and associated assimilation capacity for cells grown in distinct conditions (i.e., light quality, light quantity, or carbon status). These measurements demonstrate that the CCM modulates carbon uptake and assimilation under low-C i conditions and that light-dependent regulation of pigmentation, cell shape, and downstream stages of carbon fixation are critical for tuning carbon uptake and assimilation.

59 BASIC BIOLOGICAL SCIENCES↗

Mechanical properties of heterogeneous metallic glasses: Insights from brick-and-mortar designs

The brick-and-mortar architectural paradigm is commonly applied in structural design to achieve an exceptional combination of strength and plasticity due to its highly tunable mechanical attributes. This study applies molecular dynamics simulations of tensile loading to investigate the mechanical properties and failure mechanisms of brick-and-mortar metallic glasses (BMMGs). The focus is primarily on the implications of the aspect ratio of the bricks and interlayer thickness on the strength, plasticity, and deformation mechanisms. Results indicate that the failure mode is typified by multiple shear bands localized within the softer mortar regions, creating a staggered network under smaller aspect ratios that significantly enhances plastic deformation. However, a transition to a single dominant SB occurs when the aspect ratio exceeds 4.1. In addition, results indicate that a concurrent improvement in strength and plasticity can be achieved by modulating the distance between brick layers in BMMGs. This enhancement originates from the extensive generation and interaction of shear transformation zones. Furthermore, this study highlights that a desirable balance between strength and plasticity can be obtained in BMMGs with appropriate brick aspect ratios and interlayer thicknesses, providing a potential design strategy for the advancement of novel metallic glasses with superior properties.

36 MATERIALS SCIENCE↗

Synthesis of metallic glass nanoparticles by inert gas condensation

Inert Gas Condensation (IGC) is a versatile method for the synthesis of ultrafine nanoparticles. While it has been mostly used in the synthesis of nanocrystalline particles it has also been employed in the synthesis of metallic glass nanoparticles, which can be further used in the generation of nanostructured metallic glasses, known as nanoglasses. In this work, we use molecular dynamics simulations to investigate the IGC processes leading to the synthesis of Cu 64 Zr 36 amorphous nanoparticles. We investigate the growth of nanoparticles up to 7.3 nm diameter using long simulations of up to 130 ns, under relatively low pressures in the range of 4 to 10 bars. Results indicate a significant effect of pressure on the final structure of the amorphous nanoparticles, yet no clear influence on the segregation of Cu atoms to its surface. Cu segregation occurring at all pressures results in a thin Cu rich layer, about 2.5 Å thick, that coats all nanoparticles. Here, statistics of atomic Voronoi polyhedra shows that the fraction of full icosahedra in the nanoparticles produced at 4 bars is equivalent to that of bulk metallic glasses produced by melt and quench at 10 9 K/s. Higher pressures result in a lower fraction of full icosahedra.

36 MATERIALS SCIENCE↗

Tailoring Cu-Zr gradient nanoglass structures: Influence of nanoparticle size and cooling rates on glass-glass interfaces

The study of gradient nanoglasses (GNGs) has gained attention due to their unique mechanical properties and potential applications in advanced materials. This study employs molecular dynamics simulations to synthesize a GNG using Cu-Zr metallic glass nanoparticles (NPs) sized from 3 to 15 nm. The NPs were produced by melting and quenching metallic clusters at a relatively slow quench rate of 10 9 K/s. The synthesis of GNG is elucidated along with the characterization of its heterogeneous metallic glass nanostructure. A seamless GNG structure is formed through cold compression of Cu 64 Zr 36 amorphous NPs of varying sizes. The influence of NP size on the GNG structure is investigated, utilizing deeply relaxed NPs, which exhibit a characteristic Cu segregation pattern on their surfaces. The results highlight an increase in structural heterogeneity due to heterogeneous mass transport and the development of local composition and density variations caused by Cu segregation at glass-glass interfaces (GGIs). A reduction in NP size is correlated with decreased Cu atomic displacements and local density at GGIs, suggesting that larger NPs may produce stronger GGIs. This research presents a novel methodology for synthesizing heterogeneous metallic glasses, demonstrating the capacity to control and customize nanostructure heterogeneity through the manipulation of NP sizes and cooling rates. Furthermore, these findings enhance our understanding of structural evolution during nanoglass synthesis and lay the foundation for further exploration in nanomaterial synthesis and characterization.

36 MATERIALS SCIENCE↗

An asymptotically compatible treatment of traction loading in linearly elastic peridynamic fracture

Meshfree discretizations of state-based peridynamic models are attractive due to their ability to naturally describe fracture of general materials. However, two factors conspire to prevent meshfree discretizations of state-based peridynamics from converging to corresponding local solutions as resolution is increased: quadrature error prevents an accurate prediction of bulk mechanics, and the lack of an explicit boundary representation presents challenges when applying traction loads. Herein, we develop a reformulation of the linear peridynamic solid (LPS) model to address these shortcomings, using improved meshfree quadrature, a reformulation of the nonlocal dilatation, and a consistent handling of the nonlocal traction condition to construct a model with rigorous accuracy guarantees. In particular, these improvements are designed to enforce discrete consistency in the presence of evolving fractures, whose a priori unknown location render consistent treatment difficult. In the absence of fracture, when a corresponding classical continuum mechanics model exists, our improvements provide asymptotically compatible convergence to corresponding local solutions, eliminating surface effects and issues with traction loading which have historically plagued peridynamic discretizations. When fracture occurs, our formulation automatically provides a sharp representation of the fracture surface by breaking bonds, avoiding the loss of mass. We provide rigorous error analysis and demonstrate convergence for a number of benchmarks, including manufactured solutions, free-surface, nonhomogeneous traction loading, and composite material problems. Finally, we validate simulations of brittle fracture against a recent experiment of dynamic crack branching in soda-lime glass, providing evidence that the scheme yields accurate predictions for practical engineering problems.

42 ENGINEERING↗

Real-time damage characterization for GFRCs using high-speed synchrotron X-ray phase contrast imaging

Here, we report the application of high-speed synchrotron X-ray phase contrast imaging (PCI) in real-time damage characterization for glass fiber reinforced composites (GFRCs) subjected to dynamic loading. Dynamic single-edge notched bending (DSENB) experiments on pre-notched S-2 GFRCs were performed on a modified Kolsky compression bar. During loading, the synchrotron X-ray beam penetrated through the composite specimen from the side to detect damage evolution inside the material. Entire dynamic events were recorded by a high-speed camera as image sequences.0° and 90° unidirectional and cross-ply composites were investigated. An optical imaging technique was also employed to capture similar dynamic events in comparison with the radiographic imaging. It is demonstrated that high-speed X-ray PCI had sufficient phase contrast to characterize a crack initiation at a 20-μm spatial resolution within 920 ns and track the crack geometry during propagation, thereby providing reliable data to quantify the dynamic damage resistance of GFRCs. Furthermore, being capable of recognizing microscopic damage-related features at a sub-10-μm resolution, high-speed X-ray PCI provided fundamental material failure mechanisms to reveal the essential of macroscale structural failure of composites. It can also track the damage evolution inside and between individual plies of laminated composites. However, current high-speed X-ray PCI technique only supports in-situ observation and the high timing and spatial resolutions are limited within a field of view of ~2.5 mm in square, preventing its application in the three-dimensional and larger-area damage detection for GFRC structures.

42 ENGINEERING↗

Coarsening dynamics of ternary polymer solutions with mobility and viscosity contrasts

Using phase-field simulations, we investigate the bulk coarsening dynamics of ternary polymer solutions undergoing a glass transition for two models of phase separation: diffusion only and with hydrodynamics. The glass transition is incorporated in both models by imposing mobility and viscosity contrasts between the polymer-rich and polymer-poor phases of the evolving microstructure. For microstructures composed of polymer-poor clusters in a polymer-rich matrix, the mobility and viscosity contrasts significantly hinder coarsening, effectively leading to structural arrest. For microstructures composed of polymer-rich clusters in a polymer-poor matrix, the mobility and viscosity contrasts do not impede domain growth; rather, they change the transient concentration of the polymer-rich phase, altering the shape of the discrete domains. This effect introduces several complexities to the coarsening process, including percolation inversion of the polymer-rich and polymer-poor phases—a phenomenon normally attributed to viscoelastic phase separation.

Chemistry↗

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

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 .

59 BASIC BIOLOGICAL SCIENCES↗