Twice Continuously Differentiable Extension of a Transition-Rate Family

Let $l$ and $m$ be natural numbers with $l\ge2$ and $m\ge1$, let $\beta$ be a transition-rate family on $l$ states with control dimension $m$ and rate bound $B$, let $\Delta^l\subset\mathbb{R}^l$ be the probability simplex, and let $K\ge0$ be a real number. Points of $\mathbb{R}^l\times\mathbb{R}^m$ are written $x=(\Sigma,\alpha)$ and identified with points of Euclidean space $\mathbb{R}^{l+m}$, with coordinates $x_1,\dots,x_{l+m}$, so that $x_\gamma=\Sigma^\gamma$ for $\gamma\le l$ and $x_{l+j}=\alpha^j$ for $j\le m$. For a real-valued function $f$ on an open subset of $\mathbb{R}^{l+m}$ we write $\partial_i f$ for the partial derivative of $f$ with respect to the coordinate $x_i$, and $\partial_j\partial_i f$ for $\partial_j$ applied to the function $\partial_i f$.

A pair $(U,\bar{\beta})$ is a twice continuously differentiable extension of the transition-rate family $\beta$ with derivative bound $K$ if it consists of an open set $U\subseteq\mathbb{R}^l$ with $\Delta^l\subset U$ and a family of functions $\bar{\beta}(\sigma,\gamma,\cdot,\cdot):U\times\mathbb{R}^m\to\mathbb{R}$, indexed by the ordered pairs $(\sigma,\gamma)$ with $\sigma,\gamma\in\{1,\dots,l\}$ and $\sigma\neq\gamma$, such that for every such pair:

1. (Extension.) $\bar{\beta}(\sigma,\gamma,\Sigma,\alpha)=\beta(\sigma,\gamma,\Sigma,\alpha)$ for all $(\Sigma,\alpha)\in\Delta^l\times\mathbb{R}^m$.

2. (Regularity.) $\bar{\beta}(\sigma,\gamma,\cdot,\cdot)$ is a $C^1$ map on the open set $U\times\mathbb{R}^m$, and for every $i\in\{1,\dots,l+m\}$ the partial derivative $\partial_i\bar{\beta}(\sigma,\gamma,\cdot,\cdot)$ is again a $C^1$ map on $U\times\mathbb{R}^m$.

3. (Derivative bounds.) $|\partial_i\bar{\beta}(\sigma,\gamma,x)|\le K$ and $|\partial_j\partial_i\bar{\beta}(\sigma,\gamma,x)|\le K$ for all $i,j\in\{1,\dots,l+m\}$ and all $x\in U\times\mathbb{R}^m$.

4. (Uniform continuity of second derivatives.) For every real $\varepsilon>0$ there is a real $\delta>0$ such that $|\partial_j\partial_i\bar{\beta}(\sigma,\gamma,x)-\partial_j\partial_i\bar{\beta}(\sigma,\gamma,y)|\le\varepsilon$ for all $i,j\in\{1,\dots,l+m\}$ and all $x,y\in U\times\mathbb{R}^m$ whose Euclidean distance satisfies $d(x,y)\le\delta$.