Many bacterial species that cannot sporulate, such as the model bacterium Escherichia coli, can nevertheless survive for years, following exhaustion of external resources, in a state termed long-term stationary phase (LTSP). Here we describe the dynamics of E. coli adaptation during the first three years spent under LTSP. We show that during this time E. coli continuously adapts genetically, through the accumulation of mutations. For non-mutator clones, the majority of mutations accumulated appear to be adaptive under LTSP, reflected in an extremely convergent pattern of mutation accumulation. Despite the rapid and convergent manner in which populations adapt under LTSP, they continue to harbor extensive genetic variation. The dynamics of evolution of mutation rates under LTSP are particularly interesting. The emergence of mutators affects overall mutation accumulation rates as well as the mutational spectra and the ultimate spectrum of adaptive alleles acquired under LTSP. With time, mutators can evolve even higher mutation rates, through the acquisition of additional mutation-rate enhancing mutations. Different mutator and non-mutator clones within a single population and time point can display extreme variation in their mutation rates, resulting in differences in both the dynamics of adaptation and their associated deleterious burdens. Despite these differences, clones that vary greatly in their mutation rates tend to co-exist within their populations for many years, under LTSP.

Generalized-ensemble Monte Carlo simulations such as the multicanonical method and similar techniques are among the most efficient approaches for simulations of systems undergoing discontinuous phase transitions or with rugged freeenergy landscapes. As Markov chain methods, they are inherently serial computationally. It was demonstrated recently, however, that a combination of independent simulations that communicate weight updates at variable intervals allows for the efficient utilization of parallel computational resources for multicanonical simulations. Implementing this approach for the many-thread architecture provided by current generations of graphics processing units (GPUs), we show how it can be efficiently employed with of the order of 10 4 parallel walkers and beyond, thus constituting a versatile tool for Monte Carlo simulations in the era of massively parallel computing. We provide the fully documented source code for the approach applied to the paradigmatic example of the two-dimensional Ising model as starting point and reference for practitioners in the field. Solution method: The code uses a parallel variant of the multicanonical method employing many parallel walkers that accumulate a common histogram. The resulting histogram is used to determine the weight function for the next iteration. Once the iteration has converged, simulations visit all possible energies with the same probability. Restrictions: The system size and size of the population of replicas are limited depending on the memory of the GPU device used. Code repository at https://github.com/CQT-Leipzig/ cudamuca. Running time: Depends on system size (approx. 40 s for L=32 on a Tesla K20m)

We systematically investigate how the range of interaction between non-bonded monomers influences the formation of structural phases of elastic, flexible polymers. Massively parallel replica-exchange simulations of a generic, coarse-grained model, performed partly on graphics processing units and in multiple-gaussian modified ensembles, pave the way for the construction of the structural phase diagram, parametrized by interaction range and temperature. Conformational transitions between gas-like, liquid, and diverse solid (pseudo) phases are identified by microcanonical statistical inflection-point analysis. We find evidence for finite-size effects that cause the crossover of "collapse" and "freezing" transitions for very short interaction ranges.

No abstract

We report on a combined theoretical and experimental characterization of isolated Poly(3-hexylthiophene) (P3HT) chains weakly adsorbed on a reconstructed Au(001) surface. The local chain conformations of in situ deposited P3HT molecules were investigated by means of scanning tunneling microscopy. For comparison, Monte Carlo simulations of the system were performed up to a maximum chain length of 60 monomer units. The dependence of the end-to-end distance and the radius of gyration on the polymer chain length shows a good agreement between experiment and Monte Carlo simulations using simple updates for short chains.