Petsev, Nikolai D.; Stillinger, Frank H.; Debenedetti, Pablo G.
Source code for our energy-conserving reformulation of the 4-site molecular model for chiral phenomena originally introduced by Latinwo et al. [F. Latinwo, F. H. Stillinger, and P. G. Debenedetti, Molecular Model for Chirality Phenomena, J. Chem. Phys. 145, 154503 (2016)]. The reformulation includes an additional 8-body force that arises from an explicit configuration-dependent term in the potential energy function, resulting in a coarse-grained energy-conserving force field for molecular dynamics simulations of chirality phenomena. In this model, the coarse-grained interaction energy between two tetramers depends on their respective chiralities, and is controlled by a parameter λ, where favors local configurations involving tetramers of opposite chirality, and gives energetic preference to configurations involving tetramers of the same chirality. The source code is for use with the LAMMPS simulation package.
Bhattacharjee, Tapomoy; Amchin, Daniel; Alert, Ricard; Ott, Jenna; Datta, Sujit
Collective migration -- the directed, coordinated motion of many self-propelled agents -- is a fascinating emergent behavior exhibited by active matter that has key functional implications for biological systems. Extensive studies have elucidated the different ways in which this phenomenon may arise. Nevertheless, how collective migration can persist when a population is confronted with perturbations, which inevitably arise in complex settings, is poorly understood. Here, by combining experiments and simulations, we describe a mechanism by which collectively migrating populations smooth out large-scale perturbations in their overall morphology, enabling their constituents to continue to migrate together. We focus on the canonical example of chemotactic migration of Escherichia coli, in which fronts of cells move via directed motion, or chemotaxis, in response to a self-generated nutrient gradient. We identify two distinct modes in which chemotaxis influences the morphology of the population: cells in different locations along a front migrate at different velocities due to spatial variations in (i) the local nutrient gradient and in (ii) the ability of cells to sense and respond to the local nutrient gradient. While the first mode is destabilizing, the second mode is stabilizing and dominates, ultimately driving smoothing of the overall population and enabling continued collective migration. This process is autonomous, arising without any external intervention; instead, it is a population-scale consequence of the manner in which individual cells transduce external signals. Our findings thus provide insights to predict, and potentially control, the collective migration and morphology of cell populations and diverse other forms of active matter.
Muniz, Maria Carolina; Gartner III, Thomas E.; Riera, Marc; Knight, Christopher; Yue, Shuwen; Paesani, Francesco; Panagiotopoulos, Athanassios Z.
This dataset contains all data (including input files, simulation trajectories as well as other data files and analysis scripts) related to the publication "Vapor-liquid equilibrium of water with the MB-pol many-body potential" by Muniz et al. in preparation (2021). In this work, we assessed the performance of the MB-pol many-body potential with respect to water's vapor-liquid equilibrium properties. Through the use of direct coexistence molecular dynamics, we calculated properties such as coexistence densities, surface tension, vapor pressures and enthalpy of vaporization. We found that MB-pol is able to predict these properties in good agreement with experimental data. The results attest to the chemical accuracy of MB-pol and its large range of application across water's phase diagram.