Optimal control and numerical methods for hybrid stochastic SIS models

This work focuses on optimal controls of a class of stochastic SIS epidemic models under regime switching. By assuming that a decision maker can influence the infectivity period, our aim is to minimize the expected discounted cost due to illness, medical treatment, and the adverse effect on the society. In addition, a model with the incorporation of vaccination is proposed. Numerical schemes are developed by approximating the continuous-time dynamics using Markov chain approximation methods. It is demonstrated that the approximation schemes converge to the optimal strategy as the mesh size goes to zero. Numerical examples are provided to illustrate our results.     

Information percolation and cutoff for the random‐cluster model

We consider the random‐cluster model (RCM) on with parameters p∈(0,1) and q ≥ 1. This is a generalization of the standard bond percolation (with edges open independently with probability p) which is biased by a factor q raised to the number of connected components. We study the well‐known Fortuin‐Kasteleyn (FK)‐dynamics on this model where the update at an edge depends on the global geometry of the system unlike the Glauber heat‐bath dynamics for spin systems, and prove that for all small enough p (depending on the dimension) and any q>1, the FK‐dynamics exhibits the cutoff phenomenon at with a window size , where λ_∞ is the large n limit of the spectral gap of the process. Our proof extends the information percolation framework of Lubetzky and Sly to the RCM and also relies on the arguments of Blanca and Sinclair who proved a sharp mixing time bound for the planar version. A key aspect of our proof is the analysis of the effect of a sequence of dependent (across time) Bernoulli percolations extracted from the graphical construction of the dynamics, on how information propagates.     

An inequality connecting entropy distance,Fisher Information and large deviations

In this paper we introduce a new generalisation of the relative Fisher Information for Markov jump processes on a finite or countable state space, and prove an inequality which connects this object with the relative entropy and a large deviation rate functional. In addition to possessing various favourable properties, we show that this generalised Fisher Information converges to the classical Fisher Information in an appropriate limit. We then use this generalised Fisher Information and the aforementioned inequality to qualitatively study coarse-graining problems for jump processes on discrete spaces.     

Theory of chemical bonds in metalloenzymes XXII: a concerted bond-switching mechanism for the oxygen–oxygen bond formation coupled with one electron transfer for water oxidation in the oxygen-evolving complex of photosystem II

ABSTRACT

QM(UB3LYP)/MM(AMBER) calculations were performed for the locations of the transition structure (TS) of the oxygen–oxygen (O–O) bond formation in the S₄ state of the oxygen-evolving complex (OEC) of photosystem II (PSII). The natural orbital (NO) analysis of the broken-symmetry (BS) solutions was also performed to elucidate the nature of the chemical bonds at TS on the basis of several chemical indices defined by the occupation numbers of NO. The computational results revealed a concerted bond switching (CBS) mechanism for the oxygen–oxygen bond formation coupled with the one-electron transfer (OET) for water oxidation in OEC of PSII. The orbital interaction between the σ-HOMO of the Mn(IV)₄–O₍₅₎ bond and the π*-LUMO of the Mn(V)₁=O₍₆₎ bond plays an important role for the concerted O–O bond formation for water oxidation in the CaMn₄O₆ cluster of OEC of PSII. One electron transfer (OET) from the π-HOMO of the Mn(V)₁=O₍₆₎ bond to the σ*-LUMO of the Mn(IV)₄–O₍₅₎ bond occurs for the formation of electron transfer diradical, where the generated anion radical [Mn(IV)₄–O₍₅₎]^-? part is relaxed to the ?Mn(III)₄?…?O₍₅₎^- structure and the cation radical [O₍₆₎=Mn(V)₁]⁺ ? part is relaxed to the ⁺O₍₆₎–Mn(IV)₁? structure because of the charge-spin separation for the electron-and hole-doped Mn–oxo bonds. Therefore, the local spins are responsible for the one-electron reductions of Mn(IV)₄->Mn(III)₄ and Mn(V)₁->Mn(IV)₁. On the other hand, the O₍₅₎^- and O₍₆₎⁺ sites generated undergo the O–O bond formation in the CaMn₄O₆ cluster. The Ca(II) ion in the cubane- skeleton of the CaMn₄O₆ cluster assists the above orbital interactions by the lowering of the orbital energy levels of π*-LUMO of Mn(V)₁=O₍₆₎ and σ*-LUMO of Mn(IV)₄–O₍₅₎, indicating an important role of its Lewis acidity. Present CBS mechanism for the O–O bond formation coupled with one electron reductions of the high-valent Mn ions is different from the conventional radical coupling (RC) and acid-base (AB) mechanisms for water oxidation in artificial and native photosynthesis systems. The proton-coupled electron transfer (PC-OET) mechanism for the O–O bond formation is also touched in relation to the CBS-OET mechanism.     

Mixing time of an unaligned Gibbs sampler on the square

The paper concerns a particular example of the Gibbs sampler and its mixing efficiency. Coordinates of a point are rerandomized in the unit square ${[0,1]}^2$ to approach a stationary distribution with density proportional to $exp(?A^2{(u?v)}^2)$ for $(u,v)\in {[0,1]}^2$ with some large parameter $A$.Diaconis conjectured the mixing time of this process to be $O\left(A^2\right)$ which we confirm in this paper. This improves on the currently known $O(exp\left(A^2\right))$ estimate.