Let l l l and m m m be natural numbers with l ≥ 2 l\ge2 l ≥ 2 and m ≥ 1 m\ge1 m ≥ 1 , let A \mathcal{A} A be a nonempty subset of Euclidean space R m \mathbb{R}^m R m , let β \beta β be a transition-rate family on l l l states with control set A \mathcal{A} A and rate bound B B B , let ( U , V , β ˉ ) (U,V,\bar{\beta}) ( U , V , β ˉ ) be a twice continuously differentiable extension of β \beta β with derivative bound K K K , and let b ˉ \bar{b} b ˉ be the extended aggregate state drift of ( U , V , β ˉ ) (U,V,\bar{\beta}) ( U , V , β ˉ ) . Adopt the coordinate and partial-derivative notation ∂ i \partial_i ∂ i , ∂ j ∂ i \partial_j\partial_i ∂ j ∂ i of the extension definition, write Δ l \Delta^l Δ l for the probability simplex , write d d d for the Euclidean distance , and for i ∈ { 1 , … , l + m } i\in\{1,\dots,l+m\} i ∈ { 1 , … , l + m } and σ ∈ { 1 , … , l } \sigma\in\{1,\dots,l\} σ ∈ { 1 , … , l } set δ i σ = 1 \delta_{i\sigma}=1 δ iσ = 1 if i = σ i=\sigma i = σ and δ i σ = 0 \delta_{i\sigma}=0 δ iσ = 0 otherwise. Then for every γ ∈ { 1 , … , l } \gamma\in\{1,\dots,l\} γ ∈ { 1 , … , l } :
(i) (Regularity, derivative formulas, and restriction.) b ˉ γ \bar{b}^\gamma b ˉ γ is of class C 2 C^2 C 2 on U × V U\times V U × V , which is an open subset of R l + m \mathbb{R}^{l+m} R l + m as noted in clause 2 of the extension definition; in particular each ∂ i b ˉ γ \partial_i\bar{b}^\gamma ∂ i b ˉ γ exists and is of class C 1 C^1 C 1 on U × V U\times V U × V , and for all i , j ∈ { 1 , … , l + m } i,j\in\{1,\dots,l+m\} i , j ∈ { 1 , … , l + m } and x = ( Σ , α ) ∈ U × V x=(\Sigma,\alpha)\in U\times V x = ( Σ , α ) ∈ U × V :
∂ i b ˉ γ ( x ) = ∑ σ : σ ≠ γ ( δ i σ β ˉ ( σ , γ , x ) + Σ σ ∂ i β ˉ ( σ , γ , x ) − δ i γ β ˉ ( γ , σ , x ) − Σ γ ∂ i β ˉ ( γ , σ , x ) ) , \partial_i\bar{b}^\gamma(x)=\sum_{\sigma:\sigma\neq\gamma}\Big(\delta_{i\sigma}\,\bar{\beta}(\sigma,\gamma,x)+\Sigma^\sigma\,\partial_i\bar{\beta}(\sigma,\gamma,x)-\delta_{i\gamma}\,\bar{\beta}(\gamma,\sigma,x)-\Sigma^\gamma\,\partial_i\bar{\beta}(\gamma,\sigma,x)\Big), ∂ i b ˉ γ ( x ) = σ : σ = γ ∑ ( δ iσ β ˉ ( σ , γ , x ) + Σ σ ∂ i β ˉ ( σ , γ , x ) − δ iγ β ˉ ( γ , σ , x ) − Σ γ ∂ i β ˉ ( γ , σ , x ) ) ,
∂ j ∂ i b ˉ γ ( x ) = ∑ σ : σ ≠ γ ( δ i σ ∂ j β ˉ ( σ , γ , x ) + δ j σ ∂ i β ˉ ( σ , γ , x ) + Σ σ ∂ j ∂ i β ˉ ( σ , γ , x ) − δ i γ ∂ j β ˉ ( γ , σ , x ) − δ j γ ∂ i β ˉ ( γ , σ , x ) − Σ γ ∂ j ∂ i β ˉ ( γ , σ , x ) ) . \partial_j\partial_i\bar{b}^\gamma(x)=\sum_{\sigma:\sigma\neq\gamma}\Big(\delta_{i\sigma}\,\partial_j\bar{\beta}(\sigma,\gamma,x)+\delta_{j\sigma}\,\partial_i\bar{\beta}(\sigma,\gamma,x)+\Sigma^\sigma\,\partial_j\partial_i\bar{\beta}(\sigma,\gamma,x)-\delta_{i\gamma}\,\partial_j\bar{\beta}(\gamma,\sigma,x)-\delta_{j\gamma}\,\partial_i\bar{\beta}(\gamma,\sigma,x)-\Sigma^\gamma\,\partial_j\partial_i\bar{\beta}(\gamma,\sigma,x)\Big). ∂ j ∂ i b ˉ γ ( x ) = σ : σ = γ ∑ ( δ iσ ∂ j β ˉ ( σ , γ , x ) + δ jσ ∂ i β ˉ ( σ , γ , x ) + Σ σ ∂ j ∂ i β ˉ ( σ , γ , x ) − δ iγ ∂ j β ˉ ( γ , σ , x ) − δ jγ ∂ i β ˉ ( γ , σ , x ) − Σ γ ∂ j ∂ i β ˉ ( γ , σ , x ) ) .
Moreover b ˉ \bar{b} b ˉ agrees on Δ l × A \Delta^l\times\mathcal{A} Δ l × A with the aggregate state drift of β \beta β .
(ii) (First-order bounds.) ∣ ∂ i b ˉ γ ( x ) ∣ ≤ l ( B + K ) |\partial_i\bar{b}^\gamma(x)|\le l\,(B+K) ∣ ∂ i b ˉ γ ( x ) ∣ ≤ l ( B + K ) for all i ∈ { 1 , … , l + m } i\in\{1,\dots,l+m\} i ∈ { 1 , … , l + m } and all x ∈ Δ l × A x\in\Delta^l\times\mathcal{A} x ∈ Δ l × A ; and if in addition A \mathcal{A} A is convex , then
∣ b ˉ γ ( x ) − b ˉ γ ( y ) ∣ ≤ l + m l ( B + K ) d ( x , y ) for all x , y ∈ Δ l × A . |\bar{b}^\gamma(x)-\bar{b}^\gamma(y)|\le \sqrt{l+m}\;l\,(B+K)\;d(x,y)\qquad\text{for all }x,y\in\Delta^l\times\mathcal{A}. ∣ b ˉ γ ( x ) − b ˉ γ ( y ) ∣ ≤ l + m l ( B + K ) d ( x , y ) for all x , y ∈ Δ l × A .
(iii) (Second-order bounds and uniform continuity.) ∣ ∂ j ∂ i b ˉ γ ( x ) ∣ ≤ 3 l K |\partial_j\partial_i\bar{b}^\gamma(x)|\le 3\,l\,K ∣ ∂ j ∂ i b ˉ γ ( x ) ∣ ≤ 3 l K for all i , j ∈ { 1 , … , l + m } i,j\in\{1,\dots,l+m\} i , j ∈ { 1 , … , l + m } and x ∈ Δ l × V x\in\Delta^l\times V x ∈ Δ l × V ; and for every real ε > 0 \varepsilon>0 ε > 0 there is a real δ > 0 \delta>0 δ > 0 such that ∣ ∂ j ∂ i b ˉ γ ( x ) − ∂ j ∂ i b ˉ γ ( y ) ∣ ≤ ε |\partial_j\partial_i\bar{b}^\gamma(x)-\partial_j\partial_i\bar{b}^\gamma(y)|\le\varepsilon ∣ ∂ j ∂ i b ˉ γ ( x ) − ∂ j ∂ i b ˉ γ ( y ) ∣ ≤ ε for all i , j ∈ { 1 , … , l + m } i,j\in\{1,\dots,l+m\} i , j ∈ { 1 , … , l + m } , all γ ∈ { 1 , … , l } \gamma\in\{1,\dots,l\} γ ∈ { 1 , … , l } , and all x , y ∈ Δ l × V x,y\in\Delta^l\times V x , y ∈ Δ l × V with d ( x , y ) ≤ δ d(x,y)\le\delta d ( x , y ) ≤ δ .