BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

David A. Sivak

David A. Sivak contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechBiological PhysicsGenomicsphysics.chem-phQuantitative MethodsSubcellular Processes

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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