
There is a growing awareness that catastrophic phenomena in biology and
medicine can be mathematically represented in terms of saddlenode
bifurcations. In particular, the term `tipping', or critical transition has in
recent years entered the discourse of the general public in relation to
ecology, medicine, and public health. The saddlenode bifurcation and its
associated theory of catastrophe as put forth by Thom and Zeeman has seen
applications in a wide range of fields including molecular biophysics,
mesoscopic physics, and climate science. In this paper, we investigate a simple
model of a nonautonomous system with a timedependent parameter $p(\tau)$ and
its corresponding `dynamic' (timedependent) saddlenode bifurcation by the
modern theory of nonautonomous dynamical systems. We show that the actual
point of no return for a system undergoing tipping can be significantly delayed
in comparison to the {\em breaking time} $\hat{\tau}$ at which the
corresponding autonomous system with a timeindependent parameter $p_{a}=
p(\hat{\tau})$ undergoes a bifurcation. A dimensionless parameter
$\alpha=\lambda p_0^3V^{2}$ is introduced, in which $\lambda$ is the curvature
of the autonomous saddlenode bifurcation according to parameter $p(\tau)$,
which has an initial value of $p_{0}$ and a constant rate of change $V$. We
find that the breaking time $\hat{\tau}$ is always less than the actual point
of no return $\tau^*$ after which the critical transition is irreversible;
specifically, the relation $\tau^*\hat{\tau}\simeq 2.338(\lambda
V)^{\frac{1}{3}}$ is analytically obtained. For a system with a small $\lambda
V$, there exists a significant window of opportunity $(\hat{\tau},\tau^*)$
during which rapid reversal of the environment can save the system from
catastrophe.

We study finite state random dynamical systems (RDS) and their induced Markov
chains (MC) as stochastic models for complex dynamics. The linear
representation of deterministic maps in RDS are matrixvalued random variables
whose expectations correspond to the transition matrix of the MC. The
instantaneous Gibbs entropy, ShannonKhinchin entropy per step, and the entropy
production rate of the MC are discussed. These three concepts as key anchor
points in stochastic dynamics, characterize respectively the uncertainties of
the system at instant time $t$, the randomness generated per step, and the
dynamical asymmetry with respect to time reversal. The entropy production rate,
expressed in terms of the cycle distributions, has found an expression in terms
of the probability of the deterministic maps with the single attractor in the
maximum entropy RDS. For finite RDS with invertible transformations, the
nonnegative entropy production rate of its MC is bounded above by the
KullbackLeibler divergence of the probability of the deterministic maps with
respect to its timereversal dual probability.

Inference in hidden Markov model has been challenging in terms of scalability
due to dependencies in the observation data. In this paper, we utilize the
inherent memory decay in hidden Markov models, such that the forward and
backward probabilities can be carried out with subsequences, enabling efficient
inference over long sequences of observations. We formulate this forward
filtering process in the setting of the random dynamical system and there exist
Lyapunov exponents in the i.i.d random matrices production. And the rate of the
memory decay is known as $\lambda_2\lambda_1$, the gap of the top two Lyapunov
exponents almost surely. An efficient and accurate algorithm is proposed to
numerically estimate the gap after the softmax parametrization. The length of
subsequences $B$ given the controlled error $\epsilon$ is
$B=\log(\epsilon)/(\lambda_2\lambda_1)$. We theoretically prove the validity
of the algorithm and demonstrate the effectiveness with numerical examples. The
method developed here can be applied to widely used algorithms, such as
minibatch stochastic gradient method. Moreover, the continuity of Lyapunov
spectrum ensures the estimated $B$ could be reused for the nearby parameter
during the inference.

We revisit the celebrated WilemskiFixman (WF) treatment for the looping time
of a freedraining polymer. The WF theory introduces a sink term into the
FokkerPlanck equation for the $3(N+1)$dimensional OrnsteinUhlenbeck process
of the polymer dynamics, which accounts for the appropriate boundary condition
due to the formation of a loop. The assumption for WF theory is considerably
relaxed. A perturbation method approach is developed that justifies and
generalizes the previous results using either a Delta sink or a Heaviside sink.
For both types of sinks, we show that under the condition of a small
dimensionless $\epsilon$, the ratio of capture radius to the Kuhn length, we
are able to systematically produce all known analytical and asymptotic results
obtained by other methods. This includes most notably the transition regime
between the $N^2$ scaling of Doi, and $N\sqrt{N}/\epsilon$ scaling of Szabo,
Schulten, and Schulten. The mathematical issue at play is the nonuniform
convergence of $\epsilon\to 0$ and $N\to\infty$, the latter being an inherent
part of the theory of a Gaussian polymer. Our analysis yields a novel term in
the analytical expression for the looping time with small $\epsilon$, which is
previously unknown. Monte Carlo numerical simulations corroborate the
analytical findings. The systematic method developed here can be applied to
other systems modeled by multidimensional Smoluchowski equations.