Source author record

Hao Ge

Hao Ge 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

28works
15topics
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

28 published item(s)

preprint2022arXiv

Direct measurement of acoustic spectral density and fractional topological charge

Local density-of-states (LDOS) is a fundamental spectral property that plays a central role in various physical phenomena such as wave-matter interactions. Here, we report on the direct measurement of the LDOS of acoustic systems and derive from which the fractional topological number in an acoustic Su-Schrieffer-Heeger system. The acoustic LDOS is quantified here with a state-of-the-art technique through the measurement of the volume flow rate and the acoustic pressure with a local excitation-probe configuration. Based on this method, we study the acoustic Purcell effect and establish experimentally the important relation between the near-field LDOS and the far-field acoustic emission power. Moreover, we detect the LDOS in the one-dimensional acoustic Su-Schrieffer-Heeger model and observe the fractional topological number of the system. Our work unveils the important role of the LDOS in acoustic phenomena and paves the way toward characterizing and tailoring the LDOS in topological systems.

preprint2022arXiv

Sample-path large deviation principle for a 2-D stochastic interacting vortex dynamics with singular kernel

We consider a stochastic interacting vortex system of $N$ particles, approximating the vorticity formulation of 2-D Navier-Stokes equation on torus. The singular interaction kernel is given by the Biot-Savart law. We only require the initial state to have finite energy, and obtain a sample-path large deviation principle for the empirical measure when the number of vortices goes to infinity. The rate function is characterized by an explicit formula supporting on sample paths with finite energy and finite integral of $L^2$ norms over time. The proof utilizes a symmetrization technique for the representation of singular kernel, together with a detailed regularity analysis of the sample path with finite rate function. The key step is to prove that the singular term after symmetrization can be bounded by the integral of $L^2$ norms along sample paths.

preprint2022arXiv

The Nonequilibrium Mechanism of Noise Enhancer synergizing with Activator in HIV Latency Reactivation

Noise-modulating chemicals can synergize with transcriptional activators in reactivating latent HIV to eliminate latent HIV reservoirs. To understand the underlying biomolecular mechanism, we investigate a previous two-gene-state model and identify two necessary conditions for the synergy: an assumption of inhibition effect of transcription activators on noise enhancers; and frequent transitions to the gene non-transcription-permissive state. We then develop a loop-four-gene-state model with Tat transcription/translation and find that drug synergy is mainly determined by the magnitude and direction of energy input into the genetic regulatory kinetics of the HIV promoter. The inhibition effect of transcription activators is actually a phenomenon of energy dissipation in the nonequilibrium gene transition system. Overall, the loop-four-state model demonstrates that energy dissipation plays a crucial role in HIV latency reactivation, which might be useful for improving drug effects and identifying other synergies on lentivirus latency reactivation.

preprint2020arXiv

NN-PARS: A Parallelized Neural Network Based Circuit Simulation Framework

The shrinking of transistor geometries as well as the increasing complexity of integrated circuits, significantly aggravate nonlinear design behavior. This demands accurate and fast circuit simulation to meet the design quality and time-to-market constraints. The existing circuit simulators which utilize lookup tables and/or closed-form expressions are either slow or inaccurate in analyzing the nonlinear behavior of designs with billions of transistors. To address these shortcomings, we present NN-PARS, a neural network (NN) based and parallelized circuit simulation framework with optimized event-driven scheduling of simulation tasks to maximize concurrency, according to the underlying GPU parallel processing capabilities. NN-PARS replaces the required memory queries in traditional techniques with parallelized NN-based computation tasks. Experimental results show that compared to a state-of-the-art current-based simulation method, NN-PARS reduces the simulation time by over two orders of magnitude in large circuits. NN-PARS also provides high accuracy levels in signal waveform calculations, with less than $2\%$ error compared to HSPICE.

preprint2020arXiv

What's the relationship between CNNs and communication systems?

The interpretability of Convolutional Neural Networks (CNNs) is an important topic in the field of computer vision. In recent years, works in this field generally adopt a mature model to reveal the internal mechanism of CNNs, helping to understand CNNs thoroughly. In this paper, we argue the working mechanism of CNNs can be revealed through a totally different interpretation, by comparing the communication systems and CNNs. This paper successfully obtained the corresponding relationship between the modules of the two, and verified the rationality of the corresponding relationship with experiments. Finally, through the analysis of some cutting-edge research on neural networks, we find the inherent relation between these two tasks can be of help in explaining these researches reasonably, as well as helping us discover the correct research direction of neural networks.

preprint2019arXiv

Defending from adversarial examples with a two-stream architecture

In recent years, deep learning has shown impressive performance on many tasks. However, recent researches showed that deep learning systems are vulnerable to small, specially crafted perturbations that are imperceptible to humans. Images with such perturbations are the so called adversarial examples, which have proven to be an indisputable threat to the DNN based applications. The lack of better understanding of the DNNs has prevented the development of efficient defenses against adversarial examples. In this paper, we propose a two-stream architecture to protect CNN from attacking by adversarial examples. Our model draws on the idea of "two-stream" which commonly used in the security field, and successfully defends different kinds of attack methods by the differences of "high-resolution" and "low-resolution" networks in feature extraction. We provide a reasonable interpretation on why our two-stream architecture is difficult to defeat, and show experimentally that our method is hard to defeat with state-of-the-art attacks. We demonstrate that our two-stream architecture is robust to adversarial examples built by currently known attacking algorithms.

preprint2019arXiv

Universal Relation Between Thermodynamic Driving Force and One-Way Fluxes in a Nonequilibrium Chemical Reaction with Complex Mechanism

In nonequilibrium chemical reaction systems, a fundamental relationship between unbalanced kinetic one-way fluxes and thermodynamic chemical driving forces is believed to exists. However this relation has been rigorously demonstrated only in a few cases in which one-way fluxes are well defined. In terms of its stochastic kinetic representation, we formulate the one-way fluxes for a general chemical reaction far from equilibrium, with arbitrary complex mechanisms, multiple intermediates, and internal kinetic cycles. For each kinetic cycle, the logarithm of the ratio of the steady-state forward and backward one-way fluxes is equal to the free energy difference between the reactants and products along the cycle. This fundamental relation is further established for general chemical reaction networks with multiple input and output complexes. Our result not only provides an equivalent definition of free energy difference in nonequilibrium chemical reaction networks, it also unifies the stochastic and macroscopic nonequilibrium chemical thermodynamics in a very broad sense.

preprint2016arXiv

Mathematical Formalism of Nonequilibrium Thermodynamics for Nonlinear Chemical Reaction Systems with General Rate Law

This paper studies a mathematical formalism of nonequilibrium thermodynamics for chemical reaction models with $N$ species, $M$ reactions, and general rate law. We establish a mathematical basis for J. W. Gibbs' macroscopic chemical thermodynamics under G. N. Lewis' kinetic law of entire equilibrium (detailed balance in nonlinear chemistry kinetics). In doing so, the equilibrium thermodynamics is then naturally generalized to nonequilibrium settings without detailed balance. The kinetic models are represented by a Markovian jumping process. A generalized macroscopic chemical free energy function and its associated balance equation with nonnegative source and sink are the major discoveries. The proof is based on the large deviation principle of this type of Markov processes. A general fluctuation dissipation theorem for stochastic reaction kinetics is also proved. The mathematical theory illustrates how a novel macroscopic dynamic law can emerges from the mesoscopic kinetics in a multi-scale system.

preprint2016arXiv

Mesoscopic Kinetic Basis of Macroscopic Chemical Thermodynamics: A Mathematical Theory

From a mathematical model that describes a complex chemical kinetic system of $N$ species and $M$ elementrary reactions in a rapidly stirred vessel of size $V$ as a Markov process, we show that a macroscopic chemical thermodynamics emerges as $V\rightarrow\infty$. The theory is applicable to linear and nonlinear reactions, closed systems reaching chemical equilibrium, or open, driven systems approaching to nonequilibrium steady states. A generalized mesoscopic free energy gives rise to a macroscopic chemical energy function $φ^{ss}(\vx)$ where $\vx=(x_1,\cdots,x_N)$ are the concentrations of the $N$ chemical species. The macroscopic chemical dynamics $\vx(t)$ satisfies two emergent laws: (1) $(\rd/\rd t)φ^{ss}[\vx(t)]\le 0$, and (2)$(\rd/\rd t)φ^{ss}[\vx(t)]=\text{cmf}(\vx)-σ(\vx)$ where entropy production rate $σ\ge 0$ represents the sink for the chemical energy, and chemical motive force $\text{cmf}\ge 0$ is non-zero if the system is driven under a sustained nonequilibrium chemostat. For systems with detailed balance $\text{cmf}=0$, and if one assumes the law of mass action,$φ^{ss}(\vx)$ is precisely the Gibbs' function $\sum_{i=1}^N x_i\big[μ_i^o+\ln x_i\big]$ for ideal solutions. For a class of kinetic systems called complex balanced, which include many nonlinear systems as well as many simple open, driven chemical systems, the $φ^{ss}(\vx)$, with global minimum at $\vx^*$, has the generic form $\sum_{i=1}^N x_i\big[\ln(x_i/x_i^*)-x_i+x_i^*\big]$,which has been known in chemical kinetic literature.Macroscopic emergent "laws" are independent of the details of the underlying kinetics. This theory provides a concrete example from chemistry showing how a dynamic macroscopic law can emerge from the kinetics at a level below.

preprint2016arXiv

Nonequilibrium steady state of biochemical cycle kinetics under non-isothermal conditions

Nonequilibrium steady state of isothermal biochemical cycle kinetics has been extensively studied, but much less investigated under non-isothermal conditions. However, once the heat exchange between subsystems is rather slow, the isothermal assumption of the whole system meets great challenge, which is indeed the case inside many kinds of living organisms. Here we generalize the nonequilibrium steady-state theory of isothermal biochemical cycle kinetics, in the master-equation models, to the situation in which the temperatures of subsystems can be far from uniform. We first obtain a new thermodynamic relation between the chemical reaction rates and thermodynamic potentials under such a non-isothermal circumstances, which immediately implies simply applying the isothermal transition-state rate formula for each chemical reaction in terms of only the reactants' temperature, is not thermodynamically consistent. Therefore, we mathematically derive several revised reaction-rate formulas which not only obey the new thermodynamic relation but also approximate the exact reaction rate better than the rate formula under isothermal condition. The new thermodynamic relation also predicts that in the transporter system with different temperatures inside and outside the membrane, the net flux of the transported molecules can possibly even go against the temperature gradient in the absence of the chemical driving force.

preprint2016arXiv

Nonequilibrium Thermodynamic Formalism of Nonlinear Chemical Reaction Systems with Waage-Guldberg's Law of Mass Action

Macroscopic entropy production $σ^{(tot)}$ in the general nonlinear isothermal chemical reaction system with mass action kinetics is decomposed into a free energy dissipation and a house-keeping heat: $σ^{(tot)}=σ^{(fd)}+σ^{(hk)}$; $σ^{(fd)}=-\rd A/\rd t$, where $A$ is a generalized free energy function. This yields a novel nonequilibrium free energy balance equation $\rd A/\rd t=-σ^{(tot)}+σ^{(hk)}$, which is on a par with celebrated entropy balance equation $\rd S/\rd t=σ^{(tot)}+η^{(ex)}$ where $η^{(ex)}$ is the rate of entropy exchange with the environment.For kinetic systems with complex balance, $σ^{(fd)}$ and $σ^{(hk)}$ are the macroscopic limits of stochastic free energy dissipation and house-keeping heat, which are both nonnegative, in the Delbrück-Gillespie description of the stochastic chemical kinetics.Therefore, we show that a full kinetic and thermodynamic theory of chemical reaction systems that transcends mesoscopic and macroscopic levels emerges.

preprint2015arXiv

Acoustic topological insulator and robust one-way sound transport

Discovery of novel topological orders of condensed matters is of a significant interest in both fundamental and applied physics due to the associated quantum conductance behaviors and unique symmetry-protected backscattering-immune propagation against defects, which inspired similar fantastic effects in classical waves system, leading to the revolution of the manipulation of wave propagation. To date, however, only few theoretical models were proposed to realize acoustic topological states. Here, we theoretically and experimentally demonstrate a two dimensional acoustic topological insulators with acoustic analogue of quantum spin Hall Effect. Due to the band inversion mechanism near the double Dirac cones, acoustic one-way pseudospin dependent propagating edge states, corresponding to spin-plus and spin-minus, can be observed at the interface between two graphene-like acoustic crystals. We have also experimentally verified the associated topological immunity of such one-way edge states against the different lattice defects and disorders, which can always lead to inherent propagation loss and noise. We show that this unique acoustic topological phenomenon can offer a new promising application platform for the design of novel acoustic devices, such as one-way sound isolators, acoustic mode switchers, splitters, filters etc.

preprint2015arXiv

Macroscopic fluxes and local reciprocal relation in second-order stochastic processes far from equilibrium

Stochastic process is an essential tool for the investigation of the physical and life sciences at nanoscale. In the first-order stochastic processes widely used in chemistry and biology, only the flux of mass rather than that of heat can be well defined. Here we investigate the two macroscopic fluxes in second-order stochastic processes driven by position-dependent forces and temperature gradient. We prove that the thermodynamic equilibrium defined through the vanishing of macroscopic fluxes is equivalent to that defined via time reversibility at mesoscopic scale. In the small noise limit, we find that the entropy production rate, which has previously been defined by the mesoscopic irreversible fluxes on the phase space, matches the classic macroscopic expression as the sum of the products of macroscopic fluxes and their associated thermodynamic forces. Further we show that the two pairs of forces and fluxes in such a limit follow a linear phenomenonical relation and the associated scalar coefficients always satisfy the reciprocal relation for both transient and steady states. The scalar coefficient is proportional to the square of local temperature divided by the local frictional coefficient and originated from the second moment of velocity distribution along each dimension. This result suggests the very close connection between Soret effect (thermal diffusion) and Dufour effect at nano scale even far from equilibrium.

preprint2015arXiv

Stochastic phenotype transition of a single cell in an intermediate region of gene-state switching

Multiple phenotypic states often arise in a single cell with different gene-expression states that undergo transcription regulation with positive feedback. Recent experiments have shown that at least in E. coli, the gene state switching can be neither extremely slow nor exceedingly rapid as many previous theoretical treatments assumed. Rather it is in the intermediate region which is difficult to handle mathematically.Under this condition, from a full chemical-master-equation description we derive a model in which the protein copy-number, for a given gene state, follow a deterministic mean-field description while the protein synthesis rates fluctuate due to stochastic gene-state switching. The simplified kinetics yields a nonequilibrium landscape function, which, similar to the energy function for equilibrium fluctuation, provides the leading orders of fluctuations around each phenotypic state, as well as the transition rates between the two phenotypic states. This rate formula is analogous to Kramers theory for chemical reactions. The resulting behaviors are significantly different from the two limiting cases studied previously.

preprint2015arXiv

Stochastic robustness and relative stability of multiple pathways in biological networks

Multiple dynamic pathways always exist in biological networks, but their robustness against internal fluctuations and relative stability have not been well recognized and carefully analyzed yet. Here we try to address these issues through an illustrative example, namely the Siah-1/beta-catenin/p14/19 ARF loop of protein p53 dynamics. Its deterministic Boolean network model predicts that two parallel pathways with comparable magnitudes of attractive basins should exist after the protein p53 is activated when a cell becomes harmfully disturbed. Once the low but non-neglectable intrinsic fluctuations are incorporated into the model, we show that a phase transition phenomenon is emerged: in one parameter region the probability weights of the normal pathway, reported in experimental literature, are comparable with the other pathway which is seemingly abnormal with the unknown functions, whereas, in some other parameter regions, the probability weight of the abnormal pathway can even dominate and become globally attractive. The theory of exponentially perturbed Markov chains is applied and further generalized in order to quantitatively explain such a phase transition phenomenon, in which the nonequilibrium "activation energy barriers" along each transiting trajectory between the parallel pathways and the number of "optimal transition paths" play a central part. Our theory can also determine how the transition time and the number of optimal transition paths between the parallel pathways depend on each interaction's strength, and help to identify those possibly more crucial interactions in the biological network.

preprint2015arXiv

The Large Deviation Principle and Steady-state Fluctuation Theorem for the Entropy Production Rate of a Stochastic Process in Magnetic Fields

Fluctuation theorem is one of the major achievements in the field of nonequilibrium statistical mechanics during the past two decades. Steady-state fluctuation theorem of sample entropy production rate in terms of large deviation principle for diffusion processes have not been rigorously proved yet due to technical difficulties. Here we give a proof for the steady-state fluctuation theorem of a diffusion process in magnetic fields, with explicit expressions of the free energy function and rate function. The proof is based on the Karhunen-Loéve expansion of complex-valued Ornstein-Uhlenbeck process.

preprint2015arXiv

Time reversibility and nonequilibrium thermodynamics of second-order stochastic processes

Nonequilibrium thermodynamics of a general second-order stochastic system is investigated. We prove that at steady state, under inversion of velocities, the condition of time-reversibility over the phase space is equivalent to the antisymmetry of spatial flux and the symmetry of velocity flux. Then we show that the condition of time-reversibility alone could not always guarantee the Maxwell-Boltzmann distribution. Comparing the two conditions together, we found that the frictional force naturally emerges as the unique odd term of the total force at thermodynamic equilibrium, and is followed by the Einstein relation. The two conditions respectively correspond to two previously reported different entropy production rates. In case that the external force is only position-dependent, the two entropy production rates become one. We prove that such an entropy production rate can be decomposed into two nonnegative terms, expressed respectively by the conditional mean and variance of the thermodynamic force associated with the irreversible velocity flux at any given spatial coordinate. In the small inertia limit, the former term becomes the entropy production rate of the overdamped dynamics; while the anomalous entropy production rate originated from the latter term. Furthermore, regarding the connection between the First Law and Second Law, we found that in the steady state of such a limit, the anomalous entropy production rate is also the leading order of the Boltzmann-factor weighted difference between the spatial heat dissipation densities of the underdamped and overdamped dynamics, while their unweighted difference always tends to vanish.

preprint2013arXiv

Dissipation, Generalized Free Energy, and a Self-consistent Nonequilibrium Thermodynamics of Chemically Driven Open Subsystems

Nonequilibrium thermodynamics of a system situated in a sustained environment with influx and efflux is usually treated as a subsystem in a larger, closed "universe". It remains a question what the minimally required description for the surrounding of such an open driven system is, so that its nonequilibrium thermodynamics can be established solely based on the internal stochastic kinetics. We provide a solution to this problem using insights from studies of molecular motors in a chemical nonequilibrium steady state (NESS) with sustained external drive through a regenerating system, or in a quasi-steady state (QSS) with an excess amount of ATP, ADP, and Pi. We introduce the key notion of {\em minimal work} that is needed, $W_{min}$, for the external regenerating system to sustain a NESS ({\em e.g.}, maintaining constant concentrations of ATP, ADP and Pi for a molecular motor). Using a Markov (master-equation) description of a motor protein, we illustrate that the NESS and QSS have identical kinetics as well as the Second Law in terms of a same positive entropy production rate. The difference between the heat dissipation of a NESS and its corresponding QSS is exactly the $W_{min}$. This provides a justification for introducing an {\em ideal external regenerating system} and yields a {\em free energy balance equation} between the net free energy input $F_{in}$ and total dissipation $F_{dis}$ in an NESS: $F_{in}$ consists of chemical input minus mechanical output; $F_{dis}$ consists of dissipative heat; and the amount of useful energy becoming heat is the NESS entropy production. Furthermore, we show that for non-stationary systems, the $F_{dis}$ and $F_{in}$ correspond to the entropy production rate and housekeeping heat in stochastic thermodynamics, and identify a relative entropy $H$ as a generalized free energy.

preprint2012arXiv

Analytical Mechanics in Stochastic Dynamics: Most Probable Path, Large-Deviation Rate Function and Hamilton-Jacobi Equation

Analytical (rational) mechanics is the mathematical structure of Newtonian deterministic dynamics developed by D'Alembert, Langrange, Hamilton, Jacobi, and many other luminaries of applied mathematics. Diffusion as a stochastic process of an overdamped individual particle immersed in a fluid, initiated by Einstein, Smoluchowski, Langevin and Wiener, has no momentum since its path is nowhere differentiable. In this exposition, we illustrate how analytical mechanics arises in stochastic dynamics from a randomly perturbed ordinary differential equation $dX_t=b(X_t)dt+εdW_t$ where $W_t$ is a Brownian motion. In the limit of vanishingly small $ε$, the solution to the stochastic differential equation other than $\dot{x}=b(x)$ are all rare events. However, conditioned on an occurence of such an event, the most probable trajectory of the stochastic motion is the solution to Lagrangian mechanics with $\mathcal{L}=\|\dot{q}-b(q)\|^2/4$ and Hamiltonian equations with $H(p,q)=\|p\|^2+b(q)\cdot p$. Hamiltonian conservation law implies that the most probable trajectory for a "rare" event has a uniform "excess kinetic energy" along its path. Rare events can also be characterized by the principle of large deviations which expresses the probability density function for $X_t$ as $f(x,t)=e^{-u(x,t)/ε}$, where $u(x,t)$ is called a large-deviation rate function which satisfies the corresponding Hamilton-Jacobi equation. An irreversible diffusion process with $\nabla\times b\neq 0$ corresponds to a Newtonian system with a Lorentz force $\ddot{q}=(\nabla\times b)\times \dot{q}+1/2\nabla\|b\|^2$. The connection between stochastic motion and analytical mechanics can be explored in terms of various techniques of applied mathematics, for example, singular perturbations, viscosity solutions, and integrable systems.

preprint2012arXiv

Landscapes of Non-gradient Dynamics Without Detailed Balance: Stable Limit Cycles and Multiple Attractors

Landscape is one of the key notions in literature on biological processes and physics of complex systems with both deterministic and stochastic dynamics. The large deviation theory (LDT) provides a possible mathematical basis for the scientists' intuition. In terms of Freidlin-Wentzell's LDT, we discuss explicitly two issues in singularly perturbed stationary diffusion processes arisen from nonlinear differential equations: (1) For a process whose corresponding ordinary differential equation has a stable limit cycle, the stationary solution exhibits a clear separation of time scales: an exponential terms and an algebraic prefactor. The large deviation rate function attains its minimum zero on the entire stable limit cycle, while the leading term of the prefactor is inversely proportional to the velocity of the non-uniform periodic oscillation on the cycle. (2) For dynamics with multiple stable fixed points and saddles, there is in general a breakdown of detailed balance among the corresponding attractors. Two landscapes, a local and a global, arise in LDT, and a Markov jumping process with cycle flux emerges in the low-noise limit. A local landscape is pertinent to the transition rates between neighboring stable fixed points; and the global landscape defines a nonequilibrium steady state. There would be nondifferentiable points in the latter for a stationary dynamics with cycle flux. LDT serving as the mathematical foundation for emergent landscapes deserves further investigations.

preprint2012arXiv

Mesoscopic Biochemical Basis of Isogenetic Inheritance and Canalization: Stochasticity, Nonlinearity, and Emergent Landscape

Biochemical reaction systems in mesoscopic volume, under sustained environmental chemical gradient(s), can have multiple stochastic attractors. Two distinct mechanisms are known for their origins: ($a$) Stochastic single-molecule events, such as gene expression, with slow gene on-off dynamics; and ($b$) nonlinear networks with feedbacks. These two mechanisms yield different volume dependence for the sojourn time of an attractor. As in the classic Arrhenius theory for temperature dependent transition rates, a landscape perspective provides a natural framework for the system's behavior. However, due to the nonequilibrium nature of the open chemical systems, the landscape, and the attractors it represents, are all themselves {\em emergent properties} of complex, mesoscopic dynamics. In terms of the landscape, we show a generalization of Kramers' approach is possible to provide a rate theory. The emergence of attractors is a form of self-organization in the mesoscopic system; stochastic attractors in biochemical systems such as gene regulation and cellular signaling are naturally inheritable via cell division. Delbrück-Gillespie's mesoscopic reaction system theory, therefore, provides a biochemical basis for spontaneous isogenetic switching and canalization.

preprint2012arXiv

One-Dimensional Birth-Death Process and Delbrück-Gillespie Theory of Mesoscopic Nonlinear Chemical Reactions

As a mathematical theory for the stochasstic, nonlinear dynamics of individuals within a population, Delbrück-Gillespie process (DGP) $n(t)\in\mathbb{Z}^N$, is a birth-death system with state-dependent rates which contain the system size $V$ as a natural parameter. For large $V$, it is intimately related to an autonomous, nonlinear ordinary differential equation as well as a diffusion process. For nonlinear dynamical systems with multiple attractors, the quasi-stationary and stationary behavior of such a birth-death process can be underestood in terms of a separation of time scales by a $T^*\sim e^{αV}$ $(α>0)$: a relatively fast, intra-basin diffusion for $t\ll T^*$ and a much slower inter-basin Markov jump process for $t\gg T^*$. In the present paper for one-dimensional systems, we study both stationary behavior ($t=\infty$) in terms of invariant distribution $p_n^{ss}(V)$, and finite time dynamics in terms of the mean first passsage time (MFPT) $T_{n_1\rightarrow n_2}(V)$. We obtain an asymptotic expression of MFPT in terms of the "stochastic potential" $Φ(x,V)=-(1/V)\ln p^{ss}_{xV}(V)$. We show in general no continuous diffusion process can provide asymptotically accurate representations for both the MFPT and the $p_n^{ss}(V)$ for a DGP. When $n_1$ and $n_2$ belong to two different basins of attraction, the MFPT yields the $T^*(V)$ in terms of $Φ(x,V)\approx ϕ_0(x)+(1/V)ϕ_1(x)$. For systems with a saddle-node bifurcation and catastrophe, discontinuous "phase transition" emerges, which can be characterized by $Φ(x,V)$ in the limit of $V\rightarrow\infty$. In terms of time scale separation, the relation between deterministic, local nonlinear bifurcations and stochastic global phase transition is discussed. The one-dimensional theory is a pedagogic first step toward a general theory of DGP.

preprint2011arXiv

Markov processes follow from the principle of Maximum Caliber

Markov models are widely used to describe processes of stochastic dynamics. Here, we show that Markov models are a natural consequence of the dynamical principle of Maximum Caliber. First, we show that when there are different possible dynamical trajectories in a time-homogeneous process, then the only type of process that maximizes the path entropy, for any given singlet statistics, is a sequence of identical, independently distributed (i.i.d.) random variables, which is the simplest Markov process. If the data is in the form of sequentially pairwise statistics, then maximizing the caliber dictates that the process is Markovian with a uniform initial distribution. Furthermore, if an initial non-uniform dynamical distribution is known, or multiple trajectories are conditioned on an initial state, then the Markov process is still the only one that maximizes the caliber. Second, given a model, MaxCal can be used to compute the parameters of that model. We show that this procedure is equivalent to the maximum-likelihood method of inference in the theory of statistics.

preprint2011arXiv

Maximum Entropy Principle, Equal Probability a Priori and Gibbs Paradox

We prove that information-theoretic maximum entropy (MaxEnt) approach to canonical ensemble is mathematically equivalent to the classic approach of Boltzmann, Gibbs and Darwin-Fowler. The two approaches, however, "interpret" a same mathematical theorem differently; most notably observing mean-energy in the former and energy conservation in the latter. However, applying the same MaxEnt method to grand canonical ensemble fails; while carefully following the classic approach based on Boltzmann's microcanonical {\em equal probability a priori} produces the correct statistics: One does not need to invoke quantum mechanics; and there is no Gibbs paradox. MaxEnt and related minimum relative entropy principle are based on the mathematical theorem concerning large deviations of rare fluctuations. As a scientific method, it requires classic mechanics, or some other assumptions, to provide meaningful {\em prior distributions} for the expected-value based statistical inference. A naive assumption of uniform prior is not valid in statistical mechanics.

preprint2011arXiv

Thermodynamics and Geometry of Reversible and Irreversible Markov Processes

Master equation with microscopic reversibility ($q_{ij}\neq 0$ iff $q_{ji}\neq 0$) has a {\em thermodynamic superstructure} in terms of two state functions $S$, entropy, and $F$, free energy: It is discovered recently that entropy production rate $e_p=-dF/dt+Q_{hk}$ with both $-dF/dt=f_d, Q_{hk} \ge 0$. The free energy dissipation $f_d\ge 0$ reflects irreversibility in spontaneous self-organization; house-keeping heat $Q_{hk}\ge 0$ reveals broken time-symmetry in open system driven away from equilibrium. In a Riemannian geometric space, the master equation is a geodesic flow when $Q_{hk}=0$; here we show that the $e_p$ decomposition is orthogonal: $e_p$, $f_d$, $Q_{hk}$ forms a pythagorean triples. Gradient flow means {\em maximum dissipation principle} outside Onsager's regime. The presence of $Q_{hk}$ makses gradient flow no longer generally true. Thermodynamics of stochastic physics requires a new geometric perspective.

preprint2010arXiv

van't Hoff-Arrhenius Analysis of Mesoscopic and Macroscopic Dynamics of Simple Biochemical Systems: Stochastic vs. Nonlinear Bistabilities

Multistability of mesoscopic, driven biochemical reaction systems has implications to a wide range of cellular processes. Using several simple models, we show that one class of bistable chemical systems has a deterministic counterpart in the nonlinear dynamics based on the Law of Mass Action, while another class, widely known as noise-induced stochastic bistability, does not. Observing the system's volume ($V$) playing a similar role as the inverse temperature ($β$) in classical rate theory, an van't Hoff-Arrhenius like analysis is introduced. In one-dimensional systems, a transition rate between two states, represented in terms of a barrier in the landscape for the dynamics $Φ(x,V)$, $k\propto\exp\{-VΔΦ^‡(V)\}$, can be understood from a decomposition $ΔΦ^‡(V) \approxΔϕ_0^‡ Δϕ_1^‡/V$. Nonlinear bistability means $Δϕ_0^‡>0$ while stochastic bistability has $Δϕ_0^‡<0$ but $Δϕ_1^‡>0$. Stochastic bistabilities can be viewed as remants (or "ghosts) of nonlinear bifurcations or extinction phenomenon, and $Δϕ_0^‡$ and $Δϕ_1^‡$ as "enthalpic" and "entropic" barriers to a transition.

preprint2009arXiv

The Physical Origins of Entropy Production, Free Energy Dissipation and their Mathematical Representations

A complete mathematical theory of nonequilibrium thermodynamics of stochastic systems in terms of master equations is presented. As generalizations of isothermal entropy and free energy, two functions of states play central roles: the Gibbs entropy $S$ and the relative entropy $F$, which are related via the stationary distribution of the stochastic dynamics. $S$ satisfies the fundamental entropy balance equation $dS/dt=e_p-h_d/T$ with entropy production rate $e_p\ge 0$ and heat dissipation rate $h_d$, while $dF/dt=-f_d\le 0$. For closed systems that satisfy detailed balance: $Te_p(t)=f_d(t)$. For open system one has $Te_p(t)=f_d(t)+Q_{hk}(t)$ where the housekeeping heat $Q_{hk}\ge 0$ was first introduced in the phenomenological nonequilibrium steady state thermodynamics. Entropy production $e_p$ consists of free energy dissipation associated with spontaneous relaxation, $f_d$, and active energy pumping that sustains the open system $Q_{hk}$. The amount of excess heat involved in the relaxation $Q_{ex}=h_d-Q_{hk} = f_d-T(dS/dt)$.

preprint2009arXiv

Thermodynamic Limit of a Nonequilibrium Steady-State: Maxwell-Type Construction for a Bistable Biochemical System

We show that the thermodynamic limit of a bistable phosphorylation-dephosphorylation cycle has a selection rule for the "more stable" macroscopic steady state. The analysis is akin to the Maxwell construction. Based on the chemical master equation approach, it is shown that, except at a critical point, bistability disappears in the stochastic model when fluctuation is sufficiently low but unneglectable. Onsager's Gaussian fluctuation theory applies to the unique macroscopic steady state. With initial state in the basin of attraction of the "less stable" steady state, the deterministic dynamics obtained by the Law of Mass Action is a metastable phenomenon. Stability and robustness in cell biology are stochastic concepts.