A preprint reports that a model's lowest energy level changed when its negative-drift setting was changed. At -0.8, the calculation gave 1.255; at -2.5, it gave 0.399, a reported shift of -0.856. Under the stronger setting, the calculated probability-amplitude peak sat near coordinate 0.8. These are outputs from an asymmetric Ornstein-Uhlenbeck model with a stepped potential used to approximate molecular Brownian-motor behavior, not measurements from a biological system.
The question behind the calculation
The paper asks how path-space thermodynamics changes when two trajectory quantities are considered together: residence time in the positive half-plane and an integral current. In ordinary language, it describes the statistical behavior of complete paths while tracking both how long the process occupies one side of the origin and how much current accumulates. The focus is the joint structure of those quantities.
For the ion-channel application, positions below zero are treated as closed and positions at or above zero as open. The model assigns the restoring force to protein elasticity. This mapping is part of the mathematical setup, as is the use of the stepped potential to approximate molecular Brownian-motor behavior.
How the model handles rare paths
To solve the model, the authors use a Feynman-Kac gauge transformation and turn the problem into a stationary Schrödinger equation with a piecewise potential. They match parabolic-cylinder, or Weber, functions at the interface. For the numerical work, they report Mathcad, a fourth-order Runge-Kutta integrator and an iterative secant method for the matching equation.
The stepped penalty gives the calculated wave function a kink at the origin. The paper interprets that interface behavior as a first-order dynamic phase transition along individual trajectories. That interpretation is a feature of the model, not experimental evidence of a biological phase transition.
The same energy spectrum is used to build a joint large-deviation rate function for residence time and current. A large-deviation rate function assigns relative weight to unusual combinations of those quantities. The paper connects it to the spectrum through a Legendre-Fenchel transform, making the calculated eigenvalue the generating potential for the model's rare fluctuations.
A test of opposite currents
The model also states a symmetry for the spectrum. If the first control force is reversed and the second is reversed and shifted by twice the asymmetry parameter, the spectrum is unchanged. That symmetry produces a Gallavotti-Cohen-type relation between the rate functions for opposite currents at the same residence time: the rate function for the negative current minus the rate function for the positive current equals the asymmetry parameter multiplied by the current. This is a theoretical relation specific to the stated model.
Another output is chi12, a cross-sensitivity coefficient. It is defined from the mixed second derivative of the energy spectrum with respect to the two control forces, evaluated when both forces are zero, and is interpreted as the asymptotic covariance between residence time and integrated conductance. In practical terms, it measures how the two responses vary together in the model. The reciprocity result is reported only within the Gaussian fluctuation regime and at second-order perturbation theory.
The model identifies its strongest sensitivity when baseline asymmetry U0 is zero. As baseline asymmetry increases, chi12 is exponentially suppressed. In the zero-diffusion limit, the paper reports chi12 driven to zero for nonzero U0, while increasing diffusion reduces that suppression. These are parameter trends from the model, not calibrated predictions for a biological channel or motor.
The edge of the claim
The calculation also marks the limits of its Gaussian picture. As the residence-time functional approaches zero or one, corresponding to the fully closed or fully open boundaries, the appendix says the Gaussian approximation breaks down and nonlinear large-deviation effects appear. Those boundary cases are therefore outside the same approximation used for the cross-sensitivity result.
Taken together, the work extends a mathematical framework for trajectory fluctuations and offers a way to relate open-time-like behavior to current-like behavior in the model. Its biological examples remain applications or analogies of the asymmetric Ornstein-Uhlenbeck construction. The reported interface transition, fluctuation relation and sensitivity window should therefore be read as results of that construction, not as evidence that real ion channels or molecular motors obey them.
Paper data and sources
Original title: Thermodynamics of the space of trajectories governed by a combination of two additive boundary functionals
Authors: V. V. Ryazanov
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-09-29
DOI: Not available
Original paper · Full text