BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David A. Sivak

David A. Sivak appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechBiological Physicsphysics.chem-phphysics.comp-phQuantitative MethodsBiomoleculescond-mat.softGenomicsInformation Theorymath.ITSubcellular Processes

Catalog footprint

What is connected

20works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

20 published item(s)

preprint2026arXiv

Global thermodynamic manifold for conservative control of stochastic systems

Optimal control of stochastic systems plays a central role in nonequilibrium physics, with applications in the study of biological molecular motors and the design of single-molecule experiments. While exact analytic solutions to optimization problems are rare, under slow driving conditions, the problem can be reformulated geometrically solely in terms of equilibrium properties. In this framework, minimum-work protocols are geodesics on a thermodynamic manifold whose metric is a generalized friction tensor. Here, we introduce a new foundation for this friction-tensor formalism for conservatively driven systems. Under complete control of the potential energy, a global thermodynamic manifold (on which points are identified with instantaneous energy landscapes) has as its metric a full-control friction tensor. Arbitrary partial-control friction tensors arise naturally as inherited metrics on submanifolds of this global manifold. Leveraging a simple mathematical relationship between system dynamics and the geometry of the global manifold, we derive new expressions for the friction tensor that offer powerful tools for interpretation and computation of friction tensors and minimum-work protocols. Our results elucidate a connection between relaxation and dissipation in slowly driven systems and suggest optimization heuristics. We demonstrate the utility of these developments in three illustrative examples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamic and thermodynamic bounds for collective motor-driven transport

Molecular motors work collectively to transport cargo within cells, with anywhere from one to several hundred motors towing a single cargo. For a broad class of collective-transport systems, we use tools from stochastic thermodynamics to derive a new lower bound for the entropy production rate which is tighter than the second law. This implies new bounds on the velocity, efficiency, and precision of general transport systems and a set of analytic Pareto frontiers for identical motors. In a specific model, we identify conditions for saturation of these Pareto frontiers.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Efficient two-dimensional control of barrier crossing

Driven barrier crossings are pervasive in optical-trapping experiments and steered molecular-dynamics simulations. Despite the high fidelity of control, the freedom in the choice of driving protocol is rarely exploited to improve efficiency. We design protocols that reduce dissipation for rapidly driven barrier crossing under two-dimensional control of a harmonic trapping potential, controlling both trap center and stiffness. For fast driving, the minimum-dissipation protocol jumps halfway between the control-parameter endpoints. For slow driving, the minimum-dissipation protocol generically slows down and tightens the trap as it crosses the barrier, resulting in both significant energy savings and increased flux compared to naive and one-dimensional protocols (that only change trap center). Combining fast and slow results, we design protocols that improve performance at all speeds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Information Thermodynamics of the Transition-Path Ensemble

The reaction coordinate describing a transition between reactant and product is a fundamental concept in the theory of chemical reactions. Within transition-path theory, a quantitative definition of the reaction coordinate is found in the committor, which is the probability that a trajectory initiated from a given microstate first reaches the product before the reactant. Here we develop an information-theoretic origin for the committor and show how selecting transition paths from the equilibrium ensemble induces entropy production which exactly equals the information that system dynamics provide about the reactivity of trajectories. This equality of entropy production and dynamical information generation also holds at the level of arbitrary individual coordinates, providing parallel measures of the coordinate's relevance to the reaction, each of which is maximized by the committor.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multidimensional minimum-work control of a 2D Ising model

A system's configurational state can be manipulated using dynamic variation of control parameters, such as temperature, pressure, or magnetic field; for finite-duration driving, excess work is required above the equilibrium free-energy change. Minimum-work protocols in multidimensional control-parameter space have potential to significantly reduce work relative to one-dimensional control. By numerically minimizing a linear-response approximation to the excess work, we design protocols in control-parameter spaces of a 2D Ising model that efficiently drive the system from the all-down to all-up configuration. We find that such designed multidimensional protocols take advantage of more flexible control to avoid control-parameter regions of high system resistance, heterogeneously input and extract work to make use of system relaxation, and flatten the energy landscape, making accessible many configurations that would otherwise have prohibitively high energy and thus decreasing spin correlations. Relative to one-dimensional protocols, this speeds up the rate-limiting spin-inversion reaction, thereby keeping the system significantly closer to equilibrium for a wide range of protocol durations, and significantly reducing resistance and hence work.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimal control with a strong harmonic trap

Quadratic trapping potentials are widely used to experimentally probe biopolymers and molecular machines and drive transitions in steered molecular-dynamics simulations. Approximating energy landscapes as locally quadratic, we design multidimensional trapping protocols that minimize dissipation. The designed protocols are easily solvable and applicable to a wide range of systems. The approximation does not rely on either fast or slow limits and is valid for any duration provided the trapping potential is sufficiently strong. We demonstrate the utility of the designed protocols with a simple model of a periodically driven rotary motor. Our results elucidate principles of effective single-molecule manipulation and efficient nonequilibrium free-energy estimation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Performance scaling and trade-offs for collective motor-driven transport

Motor-driven intracellular transport of organelles, vesicles, and other molecular cargo is a highly collective process. An individual cargo is often pulled by a team of transport motors, with numbers ranging from only a few to several hundred. We explore the behavior of these systems using a stochastic model for transport of molecular cargo by an arbitrary number N of motors obeying linear Langevin dynamics, finding analytic solutions for the N-dependence of the velocity, precision of forward progress, energy flows between different system components, and efficiency. In two opposing regimes, we show that these properties obey simple scaling laws with N. Finally, we explore trade-offs between performance metrics as N is varied, providing insight into how different numbers of motors might be well-matched to distinct contexts where different performance metrics are prioritized.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Therapeutic algebra of immunomodulatory drug responses at single-cell resolution

Therapeutic modulation of immune states is central to the treatment of human disease. However, how drugs and drug combinations impact the diverse cell types in the human immune system remains poorly understood at the transcriptome scale. Here, we apply single-cell mRNA-seq to profile the response of human immune cells to 502 immunomodulatory drugs alone and in combination. We develop a unified mathematical model that quantitatively describes the transcriptome scale response of myeloid and lymphoid cell types to individual drugs and drug combinations through a single inferred regulatory network. The mathematical model reveals how drug combinations generate novel, macrophage and T-cell states by recruiting combinations of gene expression programs through both additive and non-additive drug interactions. A simplified drug response algebra allows us to predict the continuous modulation of immune cell populations between activated, resting and hyper-inhibited states through combinatorial drug dose titrations. Our results suggest that transcriptome-scale mathematical models could enable the design of therapeutic strategies for programming the human immune system using combinations of therapeutics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Skewed Thermodynamic Geometry and Optimal Free Energy Estimation

Free energy differences are a central quantity of interest in physics, chemistry, and biology. We develop design principles that improve the precision and accuracy of free energy estimators, which has potential applications to screening for targeted drug discovery. Specifically, by exploiting the connection between the work statistics of time-reversed protocol pairs, we develop near-equilibrium approximations for moments of the excess work and analyze the dominant contributions to the precision and accuracy of standard nonequilibrium free-energy estimators. Within linear response, minimum-dissipation protocols follow geodesics of the Riemannian metric induced by the Stokes' friction tensor. We find the next-order contribution arises from the rank-3 supra-Stokes' tensor that skews the geometric structure such that minimum-dissipation protocols follow geodesics of a generalized cubic Finsler metric. Thus, near equilibrium the supra-Stokes' tensor determines the leading-order contribution to the bias of bidirectional free-energy estimators.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modeling work-speed-accuracy trade-offs in a stochastic rotary machine

Molecular machines are stochastic systems that catalyze the energetic processes keeping living cells alive and structured. Inspired by the examples of F1-ATP synthase and the bacterial flagellum, we present a minimal model of an externally driven stochastic rotary machine. We explore the trade-offs of work, driving speed, and driving accuracy when changing driving strength, speed, and the underlying system dynamics. We find an upper bound on accuracy and work for a particular speed. Our results favor slow driving when tasked with minimizing the work-accuracy ratio and maximizing the rate of successful cycles. Finally, in the parameter regime mapping to the dynamics of F1-ATP synthase, we find a significant decay of driving accuracy at physiological rotation rates, raising questions about how ATP synthase achieves reasonable or even remarkable efficiency in vivo.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal control of protein copy number

Cell-cell communication is often achieved by secreted signaling molecules that bind membrane-bound receptors. A common class of such receptors are G-protein coupled receptors, where extracellular binding induces changes on the membrane affinity near the receptor for certain cytosolic proteins, effectively altering their chemical potential. We analyze the minimum-dissipation schedules for dynamically changing chemical potential to induce steady-state changes in protein copy-number distributions, and illustrate with analytic solutions for linear chemical reaction networks. Protocols that change chemical potential on biologically relevant timescales are experimentally accessible using optogenetic manipulations, and our framework provides non-trivial predictions about functional dynamical cell-cell interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermodynamic geometry of minimum-dissipation driven barrier crossing

We explore the thermodynamic geometry of a simple system that models the bistable dynamics of nucleic acid hairpins in single molecule force-extension experiments. Near equilibrium, optimal (minimum-dissipation) driving protocols are governed by a generalized linear response friction coefficient. Our analysis and simulations demonstrate that the friction coefficient of the driving protocols is sharply peaked at the interface between metastable regions, which leads to minimum-dissipation protocols that drive rapidly within a metastable basin, but then linger longest at the interface, giving thermal fluctuations maximal time to kick the system over the barrier. Intuitively, the same principle applies generically in free energy estimation (both in steered molecular dynamics simulations and in single-molecule experiments), provides a design principle for the construction of thermodynamically efficient coupling between stochastic objects, and makes a prediction regarding the construction of evolved biomolecular motors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Time step rescaling recovers continuous-time dynamical properties for discrete-time Langevin integration of nonequilibrium systems

When simulating molecular systems using deterministic equations of motion (e.g., Newtonian dynamics), such equations are generally numerically integrated according to a well-developed set of algorithms that share commonly agreed-upon desirable properties. However, for stochastic equations of motion (e.g., Langevin dynamics), there is still broad disagreement over which integration algorithms are most appropriate. While multiple desiderata have been proposed throughout the literature, consensus on which criteria are important is absent, and no published integration scheme satisfies all desiderata simultaneously. Additional nontrivial complications stem from simulating systems driven out of equilibrium using existing stochastic integration schemes in conjunction with recently-developed nonequilibrium fluctuation theorems. Here, we examine a family of discrete time integration schemes for Langevin dynamics, assessing how each member satisfies a variety of desiderata that have been enumerated in prior efforts to construct suitable Langevin integrators. We show that the incorporation of a novel time step rescaling in the deterministic updates of position and velocity can correct a number of dynamical defects in these integrators. Finally, we identify a particular splitting that has essentially universally appropriate properties for the simulation of Langevin dynamics for molecular systems in equilibrium, nonequilibrium, and path sampling contexts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Optimal control of transitions between nonequilibrium steady states

Biological systems fundamentally exist out of equilibrium in order to preserve organized structures and processes. Many changing cellular conditions can be represented as transitions between nonequilibrium steady states, and organisms have an interest in optimizing such transitions. Using the Hatano-Sasa Y-value, we extend a recently developed geometrical framework for determining optimal protocols so that it can be applied to systems driven from nonequilibrium steady states. We calculate and numerically verify optimal protocols for a colloidal particle dragged through solution by a translating optical trap with two controllable parameters. We offer experimental predictions, specifically that optimal protocols are significantly less costly than naive ones. Optimal protocols similar to these may ultimately point to design principles for biological energy transduction systems and guide the design of artificial molecular machines.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Using nonequilibrium fluctuation theorems to understand and correct errors in equilibrium and nonequilibrium discrete Langevin dynamics simulations

Common algorithms for computationally simulating Langevin dynamics must discretize the stochastic differential equations of motion. These resulting finite time step integrators necessarily have several practical issues in common: Microscopic reversibility is violated, the sampled stationary distribution differs from the desired equilibrium distribution, and the work accumulated in nonequilibrium simulations is not directly usable in estimators based on nonequilibrium work theorems. Here, we show that even with a time-independent Hamiltonian, finite time step Langevin integrators can be thought of as a driven, nonequilibrium physical process. Once an appropriate work-like quantity is defined -- here called the shadow work -- recently developed nonequilibrium fluctuation theorems can be used to measure or correct for the errors introduced by the use of finite time steps. In particular, we demonstrate that amending estimators based on nonequilibrium work theorems to include this shadow work removes the time step dependent error from estimates of free energies. We also quantify, for the first time, the magnitude of deviations between the sampled stationary distribution and the desired equilibrium distribution for equilibrium Langevin simulations of solvated systems of varying size. While these deviations can be large, they can be eliminated altogether by Metropolization or greatly diminished by small reductions in the time step. Through this connection with driven processes, further developments in nonequilibrium fluctuation theorems can provide additional analytical tools for dealing with errors in finite time step integrators.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Near-equilibrium measurements of nonequilibrium free energy

A central endeavor of thermodynamics is the measurement of free energy changes. Regrettably, although we can measure the free energy of a system in thermodynamic equilibrium, typically all we can say about the free energy of a non-equilibrium ensemble is that it is larger than that of the same system at equilibrium. Herein, we derive a formally exact expression for the probability distribution of a driven system, which involves path ensemble averages of the work over trajectories of the time-reversed system. From this we find a simple near-equilibrium approximation for the free energy in terms of an excess mean time-reversed work, which can be experimentally measured on real systems. With analysis and computer simulation, we demonstrate the accuracy of our approximations for several simple models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The geometry of thermodynamic control

A deeper understanding of nonequilibrium phenomena is needed to reveal the principles governing natural and synthetic molecular machines. Recent work has shown that when a thermodynamic system is driven from equilibrium then, in the linear response regime, the space of controllable parameters has a Riemannian geometry induced by a generalized friction tensor. We exploit this geometric insight to construct closed-form expressions for minimal-dissipation protocols for a particle diffusing in a one dimensional harmonic potential, where the spring constant, inverse temperature, and trap location are adjusted simultaneously. These optimal protocols are geodesics on the Riemannian manifold, and reveal that this simple model has a surprisingly rich geometry. We test these optimal protocols via a numerical implementation of the Fokker-Planck equation and demonstrate that the friction tensor arises naturally from a first order expansion in temporal derivatives of the control parameters, without appealing directly to linear response theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The thermodynamics of prediction

A system responding to a stochastic driving signal can be interpreted as computing, by means of its dynamics, an implicit model of the environmental variables. The system's state retains information about past environmental fluctuations, and a fraction of this information is predictive of future ones. The remaining nonpredictive information reflects model complexity that does not improve predictive power, and thus represents the ineffectiveness of the model. We expose the fundamental equivalence between this model inefficiency and thermodynamic inefficiency, measured by dissipation. Our results hold arbitrarily far from thermodynamic equilibrium and are applicable to a wide range of systems, including biomolecular machines. They highlight a profound connection between the effective use of information and efficient thermodynamic operation: any system constructed to keep memory about its environment and to operate with maximal energetic efficiency has to be predictive.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermodynamic metrics and optimal paths

A fundamental problem in modern thermodynamics is how a molecular-scale machine performs useful work, while operating away from thermal equilibrium without excessive dissipation. To this end, we derive a friction tensor that induces a Riemannian manifold on the space of thermodynamic states. Within the linear-response regime, this metric structure controls the dissipation of finite-time transformations, and bestows optimal protocols with many useful properties. We discuss the connection to the existing thermodynamic length formalism, and demonstrate the utility of this metric by solving for optimal control parameter protocols in a simple nonequilibrium model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Consequences of local inter-strand dehybridization for large-amplitude bending fluctuations of double-stranded DNA

The wormlike chain (WLC) model of DNA bending accurately reproduces single-molecule force-extension profiles of long (kilobase) chains. These bending statistics over large scales do not, however, establish a unique microscopic model for elasticity at the 1-10 bp scale, which holds particular interest in biological contexts. Here we examine a class of microscopic models which allow for disruption of base pairing (i.e., a `melt' or `kink', generically an `excitation') and consequently enhanced local flexibility. We first analyze the effect on the excitation free energy of integrating out the spatial degrees of freedom in a wormlike chain. Based on this analysis, we present a formulation of these models that ensures consistency with the well-established thermodynamics of melting in long chains. Using a new method to calculate cyclization statistics of short chains from enhanced-sampling Monte Carlo simulations, we compute J-factors of a meltable wormlike chain (MWLC) over a broad range of chain lengths, including very short molecules (30 bp) that have not yet been explored experimentally. For chains longer than about 120 bp, including most molecules studied to date in the laboratory, we find that melting excitations have little impact on cyclization kinetics. Strong signatures of melting, which might be resolved within typical experimental scatter, emerge only for shorter chains.

0 citations0 reviewsNo code linkedOpen access