Source author record

Aaron R. Dinner

Aaron R. Dinner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2026arXiv

Composing diffusion priors with explicit physical context via generative Gibbs sampling

Pretrained diffusion models provide powerful learned priors, but in scientific sampling the target distribution often depends on physical context that is not fully represented by one generative model. We introduce Generative Gibbs for Physics-Aware Sampling (GG-PA), a training-free framework that formulates the composition of learned partial priors and explicit physical context as inference over a joint target distribution in an augmented state space. We derive a Gibbs sampler for this joint target, show that it is asymptotically exact as the diffusion time approaches zero, and prove that in settings with quadratic interactions it remains exact at finite diffusion times. We further introduce replica exchange over diffusion time to accelerate mixing. Experiments on a double-well system, a $φ^4$ lattice model, and atomistic peptide systems show that GG-PA recovers context-induced distribution shifts and emergent collective behavior in interacting systems using partial priors without retraining. These results demonstrate GG-PA as a practical approach for combining pretrained generative priors with explicit physical context.

preprint2026arXiv

In-context learning emerges in chemical reaction networks without attention

We investigate whether chemical processes can perform in-context learning (ICL), a mode of computation typically associated with transformer architectures. ICL allows a system to infer task-specific rules from a sequence of examples without relying solely on fixed parameters. Traditional ICL relies on a pairwise attention mechanism which is not obviously implementable in chemical systems. However, we show theoretically and numerically that chemical processes can achieve ICL through a mechanism we call subspace projection, in which the entire input vector is mapped onto comparison subspaces, with the dominant projection determining the computational output. We illustrate this mechanism analytically in small chemical systems and show numerically that performance is robust to input encoding and dynamical choices, with the number of tunable degrees of freedom in the input encoding as a key limitation. Our results provide a blueprint for realizing ICL in chemical or other physical media and suggest new directions for designing adaptive synthetic chemical systems and understanding possible biological computation in cells.

preprint2022arXiv

Computing transition path theory quantities with trajectory stratification

Transition path theory computes statistics from ensembles of reactive trajectories. A common strategy for sampling reactive trajectories is to control the branching and pruning of trajectories so as to enhance the sampling of low probability segments. However, it can be challenging to apply transition path theory to data from such methods because determining whether configurations and trajectory segments are part of reactive trajectories requires looking backward and forward in time. Here, we show how this issue can be overcome efficiently by introducing simple data structures. We illustrate the approach in the context of nonequilibrium umbrella sampling (NEUS), but the strategy is general and can be used to obtain transition path theory statistics from other methods that sample segments of unbiased trajectories.

preprint2020arXiv

Long-timescale predictions from short-trajectory data: A benchmark analysis of the trp-cage miniprotein

Elucidating physical mechanisms with statistical confidence from molecular dynamics simulations can be challenging owing to the many degrees of freedom that contribute to collective motions. To address this issue, we recently introduced a dynamical Galerkin approximation (DGA) [Thiede et al. J. Phys. Chem. 150, 244111 (2019)], in which chemical kinetic statistics that satisfy equations of dynamical operators are represented by a basis expansion. Here, we reformulate this approach, clarifying (and reducing) the dependence on the choice of lag time. We present a new projection of the reactive current onto collective variables and provide improved estimators for rates and committors. We also present simple procedures for constructing suitable smoothly varying basis functions from arbitrary molecular features. To evaluate estimators and basis sets numerically, we generate and carefully validate a dataset of short trajectories for the unfolding and folding of the trp-cage miniprotein, a well-studied system. Our analysis demonstrates a comprehensive strategy for characterizing reaction pathways quantitatively.

preprint2020arXiv

Machine learning force fields and coarse-grained variables in molecular dynamics: application to materials and biological systems

Machine learning encompasses a set of tools and algorithms which are now becoming popular in almost all scientific and technological fields. This is true for molecular dynamics as well, where machine learning offers promises of extracting valuable information from the enormous amounts of data generated by simulation of complex systems. We provide here a review of our current understanding of goals, benefits, and limitations of machine learning techniques for computational studies on atomistic systems, focusing on the construction of empirical force fields from ab-initio databases and the determination of reaction coordinates for free energy computation and enhanced sampling.

preprint2020arXiv

Stratification as a general variance reduction method for Markov chain Monte Carlo

The Eigenvector Method for Umbrella Sampling (EMUS) belongs to a popular class of methods in statistical mechanics which adapt the principle of stratified survey sampling to the computation of free energies. We develop a detailed theoretical analysis of EMUS. Based on this analysis, we show that EMUS is an efficient general method for computing averages over arbitrary target distributions. In particular, we show that EMUS can be dramatically more efficient than direct MCMC when the target distribution is multimodal or when the goal is to compute tail probabilities. To illustrate these theoretical results, we present a tutorial application of the method to a problem from Bayesian statistics.

preprint2016arXiv

A cycling state that can lead to glassy dynamics in intracellular transport

Power-law dwell times have been observed for molecular motors in living cells, but the origins of these trapped states are not known. We introduce a minimal model of motors moving on a two-dimensional network of filaments, and simulations of its dynamics exhibit statistics comparable to those observed experimentally. Analysis of the model trajectories, as well as experimental particle tracking data, reveals a state in which motors cycle unproductively at junctions of three or more filaments. We formulate a master equation for these junction dynamics and show that the time required to escape from this vortex-like state can account for the power-law dwell times. We identify trends in the dynamics with the motor valency for further experimental validation. We demonstrate that these trends exist in individual trajectories of myosin II on an actin network. We discuss how cells could regulate intracellular transport and, in turn, biological function, by controlling their cytoskeletal network structures locally.

preprint2016arXiv

Eigenvector method for umbrella sampling enables error analysis

Umbrella sampling efficiently yields equilibrium averages that depend on exploring rare states of a model by biasing simulations to windows of coordinate values and then combining the resulting data with physical weighting. Here, we introduce a mathematical framework that casts the step of combining the data as an eigenproblem. The advantage to this approach is that it facilitates error analysis. We discuss how the error scales with the number of windows. Then, we derive a central limit theorem for averages that are obtained from umbrella sampling. The central limit theorem suggests an estimator of the error contributions from individual windows, and we develop a simple and computationally inexpensive procedure for implementing it. We demonstrate this estimator for simulations of the alanine dipeptide and show that it emphasizes low free energy pathways between stable states in comparison to existing approaches for assessing error contributions. We discuss the possibility of using the estimator and, more generally, the eigenvector method for umbrella sampling to guide adaptation of the simulation parameters to accelerate convergence.

preprint2016arXiv

Shape dynamics of growing cell walls

We introduce a general theoretical framework to study the shape dynamics of actively growing and remodeling surfaces. Using this framework we develop a physical model for growing bacterial cell walls and study the interplay of cell shape with the dynamics of growth and constriction. The model allows us to derive constraints on cell wall mechanical energy based on the observed dynamics of cell shape. We predict that exponential growth in cell size requires a constant amount of cell wall energy to be dissipated per unit volume. We use the model to understand and contrast growth in bacteria with different shapes such as spherical, ellipsoidal, cylindrical and toroidal morphologies. Coupling growth to cell wall constriction, we predict a discontinuous shape transformation, from partial constriction to cell division, as a function of the chemical potential driving cell-wall synthesis. Our model for cell wall energy and shape dynamics relates growth kinetics with cell geometry, and provides a unified framework to describe the interplay between shape, growth and division in bacterial cells.

preprint2015arXiv

Intergenerational continuity of cell shape dynamics in Caulobacter crescentus

We investigate the intergenerational shape dynamics of single Caulobacter crescentus cells using a novel combination of imaging techniques and theoretical modeling. We determine the dynamics of cell pole-to-pole lengths, cross-sectional widths, and medial curvatures from high accuracy measurements of cell contours. Moreover, these shape parameters are determined for over 250 cells across approximately 10000 total generations, which affords high statistical precision. Our data and model show that constriction is initiated early in the cell cycle and that its dynamics are controlled by the time scale of exponential longitudinal growth. Based on our extensive and detailed growth and contour data, we develop a minimal mechanical model that quantitatively accounts for the cell shape dynamics and suggests that the asymmetric location of the division plane reflects the distinct mechanical properties of the stalked and swarmer poles. Furthermore, we find that the asymmetry in the division plane location is inherited from the previous generation. We interpret these results in terms of the current molecular understanding of shape, growth, and division of C. crescentus

preprint2014arXiv

Scaling laws governing stochastic growth and division of single bacterial cells

Uncovering the quantitative laws that govern the growth and division of single cells remains a major challenge. Using a unique combination of technologies that yields unprecedented statistical precision, we find that the sizes of individual Caulobacter crescentus cells increase exponentially in time. We also establish that they divide upon reaching a critical multiple ($\approx$1.8) of their initial sizes, rather than an absolute size. We show that when the temperature is varied, the growth and division timescales scale proportionally with each other over the physiological temperature range. Strikingly, the cell-size and division-time distributions can both be rescaled by their mean values such that the condition-specific distributions collapse to universal curves. We account for these observations with a minimal stochastic model that is based on an autocatalytic cycle. It predicts the scalings, as well as specific functional forms for the universal curves. Our experimental and theoretical analysis reveals a simple physical principle governing these complex biological processes: a single temperature-dependent scale of cellular time governs the stochastic dynamics of growth and division in balanced growth conditions.

preprint2014arXiv

Universality in stochastic exponential growth

Recent imaging data for single bacterial cells reveal that their mean sizes grow exponentially in time and that their size distributions collapse to a single curve when rescaled by their means. An analogous result holds for the division-time distributions. A model is needed to delineate the minimal requirements for these scaling behaviors. We formulate a microscopic theory of stochastic exponential growth as a Master Equation that accounts for these observations, in contrast to existing quantitative models of stochastic exponential growth (e.g., the Black-Scholes equation or geometric Brownian motion). Our model, the stochastic Hinshelwood cycle (SHC), is an autocatalytic reaction cycle in which each molecular species catalyzes the production of the next. By finding exact analytical solutions to the SHC and the corresponding first passage time problem, we uncover universal signatures of fluctuations in exponential growth and division. The model makes minimal assumptions, and we describe how more complex reaction networks can reduce to such a cycle. We thus expect similar scalings to be discovered in stochastic processes resulting in exponential growth that appear in diverse contexts such as cosmology, finance, technology, and population growth.

preprint2012arXiv

Phase resetting reveals network dynamics underlying a bacterial cell cycle

Genomic and proteomic methods yield networks of biological regulatory interactions but do not provide direct insight into how those interactions are organized into functional modules, or how information flows from one module to another. In this work we introduce an approach that provides this complementary information and apply it to the bacterium Caulobacter crescentus, a paradigm for cell-cycle control. Operationally, we use an inducible promoter to express the essential transcriptional regulatory gene ctrA in a periodic, pulsed fashion. This chemical perturbation causes the population of cells to divide synchronously, and we use the resulting advance or delay of the division times of single cells to construct a phase resetting curve. We find that delay is strongly favored over advance. This finding is surprising since it does not follow from the temporal expression profile of CtrA and, in turn, simulations of existing network models. We propose a phenomenological model that suggests that the cell-cycle network comprises two distinct functional modules that oscillate autonomously and couple in a highly asymmetric fashion. These features collectively provide a new mechanism for tight temporal control of the cell cycle in C. crescentus. We discuss how the procedure can serve as the basis for a general approach for probing network dynamics, which we term chemical perturbation spectroscopy (CPS).

preprint2012arXiv

Requirements for contractility in disordered cytoskeletal bundles

Actomyosin contractility is essential for biological force generation, and is well understood in highly organized structures such as striated muscle. Additionally, actomyosin bundles devoid of this organization are known to contract both in vivo and in vitro, which cannot be described by standard muscle models. To narrow down the search for possible contraction mechanisms in these systems, we investigate their microscopic symmetries. We show that contractile behavior requires non-identical motors that generate large enough forces to probe the nonlinear elastic behavior of F-actin. This suggests a role for filament buckling in the contraction of these bundles, consistent with recent experimental results on reconstituted actomyosin bundles.

preprint2012arXiv

Steered Transition Path Sampling

We introduce a path sampling method for obtaining statistical properties of an arbitrary stochastic dynamics. The method works by decomposing a trajectory in time, estimating the probability of satisfying a progress constraint, modifying the dynamics based on that probability, and then reweighting to calculate averages. Because the progress constraint can be formulated in terms of occurrences of events within time intervals, the method is particularly well suited for controlling the sampling of currents of dynamic events. We demonstrate the method for calculating transition probabilities in barrier crossing problems and survival probabilities in strongly diffusive systems with absorbing states, which are difficult to treat by shooting. We discuss the relation of the algorithm to other methods.

preprint2012arXiv

Structural Responses of Quasi-2D Colloid Fluids to Excitations Elicited by Nonequilibrium Perturbations

We investigate the response of a dense monodisperse quasi-two-dimensional (q2D) colloid suspension when a particle is dragged by a constant velocity optical trap. Consistent with microrheological studies of other geometries, the perturbation induces a leading density wave and trailing wake, and we use Stokesian Dynamics (SD) simulations to parse direct colloid-colloid and hydrodynamic interactions. We go on to analyze the underlying individual particle-particle collisions in the experimental images. The displacements of particles form chains reminiscent of stress propagation in sheared granular materials. From these data, we can reconstruct steady-state dipolar flow patterns that were predicted for dilute suspensions and previously observed in granular analogs to our system. The decay of this field differs, however, from point Stokeslet calculations, indicating that the finite size of the colloids is important. Moreover, there is a pronounced angular dependence that corresponds to the surrounding colloid structure, which evolves in response to the perturbation. Put together, our results show that the response of the complex fluid is highly anisotropic owing to the fact that the effects of the perturbation propagate through the structured medium via chains of colloid-colloid collisions.

preprint2011arXiv

Flow-dependent unfolding and refolding of an RNA by nonequilibrium umbrella sampling

Nonequilibrium experiments of single biomolecules such as force-induced unfolding reveal details about a few degrees of freedom of a complex system. Molecular dynamics simulations can provide complementary information, but exploration of the space of possible configurations is often hindered by large barriers in phase space that separate metastable regions. To solve this problem, enhanced sampling methods have been developed that divide a phase space into regions and integrate trajectory segments in each region. These methods boost the probability of passage over barriers, and facilitate parallelization since integration of the trajectory segments does not require communication, aside from their initialization and termination. Here we present a parallel version of an enhanced sampling method suitable for systems driven far from equilibrium: nonequilibrium umbrella sampling (NEUS). We apply this method to a coarse-grained model of a 262-nucleotide RNA molecule that unfolds and refolds in an explicit flow field modeled with stochastic rotation dynamics. Using NEUS we are able to observe extremely rare unfolding events that have mean first passage times as long as 1.4 s (3.4 E13 dynamics steps). We examine the unfolding process for a range of flow rates of the medium, and we describe two competing pathways in which different intramolecular contacts are broken.