Researcher profile

Todd R. Gingrich

Todd R. Gingrich contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
4topics
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

10 published item(s)

preprint2026arXiv

Chemical master equation parameter exploration using DMRG

Well-mixed chemical reaction networks (CRNs) contain many distinct chemical species with copy numbers that fluctuate in correlated ways. While those correlations are typically monitored via Monte Carlo sampling of stochastic trajectories, there is interest in systematically approximating the joint distribution over the exponentially large number of possible microstates using tensor networks or tensor trains. We exploit the tensor network strategy to determine when the steady state of a seven-species gene toggle switch CRN model supports bistability as a function of two decomposition rates, both parameters of the kinetic model. We highlight how the tensor network solution captures the effects of stochastic fluctuations, going beyond mean field and indeed deviating meaningfully from a mean-field analysis. The work furthermore develops and demonstrates several technical advances that will allow steady-states of broad classes of CRNs to be computed in a manner conducive to parameter exploration. We show that the steady-state distributions can be computed via the ordinary density matrix renormalization group (DMRG) algorithm despite having a non-Hermitian rate operator with a small spectral gap, we illustrate how that steady-state distribution can be efficiently projected to an order parameter that identifies bimodality, and we employ excited-state DMRG to calculate a relaxation timescale for the bistability.

preprint2026arXiv

It Takes Two to Make a Thing Go Right: Boosting Current in Coupled Motors

Catalysis-driven synthetic molecular motors operate in a loose mechanochemical coupling regime, one in which a decomposition of a fuel molecule does not reliably produce a forward step. In that regime, stochastic backward steps can significantly degrade the motor's current, prompting us to ask whether mechanically coupling multiple such motors can boost their averaged current. By simulating rotaxane-based motors with two classes of models--particle-based nonequilibrium molecular dynamics and jump-diffusion models--we show that current boosts are physically achievable. Our observed boosts, which amplify current by single-digit factors, emerge when coupling between motors can increase the activity, speeding up the rate of both forward and backward steps. In doing so, the bias for preferring forward steps actually degrades, but the lost bias can be largely recovered by raising the fuel concentration, demonstrating a general design strategy: amplify activity through coupling and restore bias through stronger driving.

preprint2024arXiv

From high-dimensional committors to reactive insights

Transition path theory (TPT) offers a powerful formalism for extracting the rate and mechanism of rare dynamical transitions between metastable states. Most applications of TPT either focus on systems with modestly sized state spaces or use collective variables to try to tame the curse of dimensionality. Increasingly, expressive function approximators like neural networks and tensor networks have shown promise in computing the central object of TPT, the committor function, even in very high dimensional systems. That progress prompts our consideration of how one could use such a high dimensional function to extract mechanistic insight. Here, we present and illustrate a straightforward but powerful way to track how individual dynamical coordinates evolve during a reactive event. The strategy, which involves marginalizing the reactive ensemble, naturally captures the evolution of the dynamical coordinate's distribution, not just its mean reactive behavior.

preprint2023arXiv

Limits on the Precision of Catenane Molecular Motors: Insights from Thermodynamics and Molecular Dynamics Simulations

Thermodynamic uncertainty relations (TURs) relate precision to the dissipation rate, yet the inequalities can be far from saturation. Indeed, in catenane molecular motor simulations, we record precision far below the TUR limit. We further show that this inefficiency can be anticipated by four physical parameters: the thermodynamic driving force, fuel decomposition rate, coupling between fuel decomposition and motor motion, and rate of undriven motor motion. The physical insights might assist in designing molecular motors in the future.

preprint2022arXiv

Computing time-periodic steady-state currents via the time evolution of tensor network states

We present an approach based upon binary tree tensor network (BTTN) states for computing steady-state current statistics for a many-particle 1D ratchet subject to volume exclusion interactions. The ratcheted particles, which move on a lattice with periodic boundary conditions subject to a time-periodic drive, can be stochastically evolved in time to sample representative trajectories via a Gillespie method. In lieu of generating realizations of trajectories, a BTTN state can variationally approximate a distribution over the vast number of many-body configurations. We apply the density matrix renormalization group (DMRG) algorithm to initialize BTTN states, which are then propagated in time via the time-dependent variational principle (TDVP) algorithm to yield the steady-state behavior, including the effects of both typical and rare trajectories. The application of the methods to ratchet currents is highlighted in a companion letter, but the approach extends naturally to other interacting lattice models with time-dependent driving. Though trajectory sampling is conceptually and computationally simpler, we discuss situations for which the BTTN TDVP strategy could be more favorable.

preprint2022arXiv

Using tensor network states for multi-particle Brownian ratchets

The study of Brownian ratchets has taught how time-periodic driving supports a time-periodic steady state that generates nonequilibrium transport. When a single particle is transported in one dimension, it is possible to rationalize the current in terms of the potential, but experimental efforts have ventured beyond that single-body case to systems with many interacting carriers. Working with a lattice model of volume-excluding particles in one dimension, we analyze the impact of interactions on a flashing ratchet's current. To surmount the many-body problem, we employ the time-dependent variational principle with a binary tree tensor network, methods discussed at length in a companion paper. Rather than propagating individual trajectories, the tensor network approach propagates a distribution over many-body configurations via a controllable variational approximation. The calculations, which reproduce Gillespie trajectory sampling, identify and explain a shift in the frequency of maximum current to higher driving frequency as the lattice occupancy increases.

preprint2021arXiv

Simulating a Chemically-Fueled Molecular Motor with Nonequilibrium Molecular Dynamics

Most computer simulations of molecular dynamics take place under equilibrium conditions--in a closed, isolated system, or perhaps one held at constant temperature or pressure. Sometimes, extra tensions, shears, or temperature gradients are introduced to those simulations to probe one type of nonequilibrium response to external forces. Catalysts and molecular motors, however, function based on the nonequilibrium dynamics induced by a chemical reaction's thermodynamic driving force. In this scenario, simulations require chemostats capable of preserving the chemical concentrations of the nonequilibrium steady state. We develop such a dynamic scheme and use it to observe cycles of a new particle-based classical model of a catenane-like molecular motor. Molecular motors are frequently modeled with detailed-balance-breaking Markov models, and we explicitly construct such a picture by coarse graining the microscopic dynamics of our simulations in order to extract rates. This work identifies inter-particle interactions that tune those rates to create a functional motor, thereby yielding a computational playground to investigate the interplay between directional bias, current generation, and coupling strength in molecular information ratchets.

preprint2020arXiv

Current inversion in a periodically driven two-dimensional Brownian ratchet

It is well-known that Brownian ratchets can exhibit current reversals, wherein the sign of the current switches as a function of the driving frequency. We introduce a spatial discretization of such a two-dimensional Brownian ratchet to enable spectral methods that efficiently compute those currents. These discrete-space models provide a convenient way to study the Markovian dynamics conditioned upon generating particular values of the currents. By studying such conditioned processes, we demonstrate that low-frequency negative values of current arise from typical events and high-frequency positive values of current arises from rare events. We demonstrate how these observations can inform the sculpting of time-dependent potential landscapes with a specific frequency response.

preprint2020arXiv

Isometric Uncertainty Relations

We generalize the link between fluctuation theorems and thermodynamic uncertainty relations by deriving a bound on the variance of fluxes that satisfy an isometric fluctuation theorem. The resulting bound, which depends on the system's dimension $d$, naturally interpolates between two known bounds. The bound derived from the entropy production fluctuation theorem is recovered for $d=1$, and the original entropy production thermodynamic uncertainty relation is obtained in the $d \to \infty$ limit. We show that our result can be generalized to order parameters in equilibrium systems, and we illustrate the results on a Heisenberg spin chain.

preprint2020arXiv

Universal thermodynamic bounds on nonequilibrium response with biochemical applications

Diverse physical systems are characterized by their response to small perturbations. Near thermodynamic equilibrium, the fluctuation-dissipation theorem provides a powerful theoretical and experimental tool to determine the nature of response by observing spontaneous equilibrium fluctuations. In this spirit, we derive here a collection of equalities and inequalities valid arbitrarily far from equilibrium that constrain the response of nonequilibrium steady states in terms of the strength of nonequilibrium driving. Our work opens new avenues for characterizing nonequilibrium response. As illustrations, we show how our results rationalize the energetic requirements of two common biochemical motifs.