Adopt the setting and notation of Linearisation of a Perturbed Controlled Aggregate Flow Along a Comparison Pair: Exact Variation-of-Constants Identity and Residual Bound : natural numbers l ≥ 2 l\ge2 l ≥ 2 , m ≥ 1 m\ge1 m ≥ 1 , the control set A ⊆ R m \mathcal{A}\subseteq\mathbb{R}^m A ⊆ R m , real numbers B ≥ 0 B\ge0 B ≥ 0 and T > 0 T>0 T > 0 , the transition-rate family β \beta β with its twice continuously differentiable extension of derivative bound K K K , the labels c = ( σ , γ ) ∈ L c=(\sigma,\gamma)\in\mathcal{L} c = ( σ , γ ) ∈ L with vectors v c v_c v c , the label rates ψ c \psi_c ψ c , state gradients g c g^{c} g c , drift Jacobian E \mathcal{E} E , drift b b b , the constants Λ 2 \Lambda_2 Λ 2 , Λ 3 \Lambda_3 Λ 3 , Λ E \Lambda_{\mathcal{E}} Λ E , the notions measurable and bounded measurable for maps on [ 0 , T ] [0,T] [ 0 , T ] , the comparison pair ( S , A ) (S,\mathsf{A}) ( S , A ) with E u ⋆ = E ( S u , A u ) \mathcal{E}^{\star}_u=\mathcal{E}(S_u,\mathsf{A}_u) E u ⋆ = E ( S u , A u ) continuous in u u u , the two-parameter fundamental solution Φ E ( t , u ) \Phi^{\mathcal{E}}(t,u) Φ E ( t , u ) and its bound Φ ˉ \bar{\Phi} Φ ˉ , and, for a control path a a a and paths x , y x,y x , y , the objects e e e , ρ \rho ρ , d \mathsf{d} d and L t ( ⋅ ) L_t(\cdot) L t ( ⋅ ) defined there. Adopt also from Open-Loop Aggregate Solution Driven by Aggregate Transition Clocks a natural number N ≥ 1 N\ge1 N ≥ 1 , the aggregate lattice G N ⊆ Δ l \mathbb{G}_N\subseteq\Delta^l G N ⊆ Δ l , clock families, control paths, and open-loop aggregate solutions with their consumed clock times C c \mathsf{C}^{c} C c ; and from Cell-Count Form of a Compensated Counting-Path Functional Along a Time Change: Window-Discrepancy Bounds for the Time-Change and Partial-Cell Errors , for a counting path p p p and a real number R > 0 R>0 R > 0 , the compensated path M ˉ ( x ) = p ( x ) − x \bar{M}(x)=p(x)-x M ˉ ( x ) = p ( x ) − x , the window discrepancy D i s c w ( p ) \mathrm{Disc}_w(p) Disc w ( p ) (formed with this R R R ), the compensated functional Λ ( ⋅ ) \Lambda(\cdot) Λ ( ⋅ ) , and, for a bounded measurable H : [ 0 , T ] → R l H:[0,T]\to\mathbb{R}^l H : [ 0 , T ] → R l , the norm ∥ H ∥ 1 = ∫ [ 0 , T ] ∣ H u ∣ d u \lVert H\rVert_1=\int_{[0,T]}|H_u|\,du ∥ H ∥ 1 = ∫ [ 0 , T ] ∣ H u ∣ d u . Write ∣ ⋅ ∣ |\cdot| ∣ ⋅ ∣ for the Euclidean norm , x ⋅ y x\cdot y x ⋅ y for the dot product , and N \sqrt{N} N for the nonnegative square root .
Cells and clocks. Let R > 0 R>0 R > 0 be a real number with R ≥ N B T R\ge NBT R ≥ NBT ; for every label c c c let J c ≥ 1 J_c\ge1 J c ≥ 1 be a natural number and 0 = b 0 c < b 1 c < ⋯ < b J c c = R 0=\mathsf{b}^{c}_0<\mathsf{b}^{c}_1<\dots<\mathsf{b}^{c}_{J_c}=R 0 = b 0 c < b 1 c < ⋯ < b J c c = R real numbers, with cell lengths μ c , j = b j c − b j − 1 c \mu_{c,j}=\mathsf{b}^{c}_j-\mathsf{b}^{c}_{j-1} μ c , j = b j c − b j − 1 c ; let L \mathsf{L} L be the set of pairs q = ( c , j ) q=(c,j) q = ( c , j ) with 1 ≤ j ≤ J c 1\le j\le J_c 1 ≤ j ≤ J c , write μ q = μ c , j \mu_q=\mu_{c,j} μ q = μ c , j and μ max = max q ∈ L μ q \mu_{\max}=\max_{q\in\mathsf{L}}\mu_q μ m a x = max q ∈ L μ q (the cells of Injection Weights from a Bounded Profile Along Mean-Field Clocks: Label-Rate Form of the Fluctuation Covariance, Weight Bounds, the Prior Quadratic Form, and the Step-Function Injection ). Let p = ( p c ) c ∈ L p=(p^{c})_{c\in\mathcal{L}} p = ( p c ) c ∈ L be a clock family, a a a a control path, x 0 ∈ G N x_0\in\mathbb{G}_N x 0 ∈ G N , and let Σ : [ 0 , T ] → R l \Sigma:[0,T]\to\mathbb{R}^l Σ : [ 0 , T ] → R l be an open-loop aggregate solution for the data ( p , a , x 0 ) (p,a,x_0) ( p , a , x 0 ) , with consumed clock times C c \mathsf{C}^{c} C c ; define the cell counts K c , j = p c ( b j c ) − p c ( b j − 1 c ) \mathsf{K}_{c,j}=p^{c}(\mathsf{b}^{c}_j)-p^{c}(\mathsf{b}^{c}_{j-1}) K c , j = p c ( b j c ) − p c ( b j − 1 c ) and write K q = K c , j \mathsf{K}_q=\mathsf{K}_{c,j} K q = K c , j . Define the mean-field label rates and mean-field clocks
ϕ c ( t ) = ψ c ( S t , A t ) , C ˉ t c = N ∫ [ 0 , t ] ϕ c ( u ) d u ( c ∈ L , t ∈ [ 0 , T ] ) , \phi_c(t)=\psi_c(S_t,\mathsf{A}_t),\qquad \bar{\mathsf{C}}^{c}_t=N\int_{[0,t]}\phi_c(u)\,du\qquad(c\in\mathcal{L},\ t\in[0,T]), ϕ c ( t ) = ψ c ( S t , A t ) , C ˉ t c = N ∫ [ 0 , t ] ϕ c ( u ) d u ( c ∈ L , t ∈ [ 0 , T ]) ,
as in Injection Weights from a Bounded Profile Along Mean-Field Clocks: Label-Rate Form of the Fluctuation Covariance, Weight Bounds, the Prior Quadratic Form, and the Step-Function Injection , whose cell boundaries b j c b^{c}_j b j c are written b j c \mathsf{b}^{c}_j b j c here (each ϕ c \phi_c ϕ c is measurable with values in [ 0 , B ] [0,B] [ 0 , B ] by claim 1). Let y : [ 0 , T ] → Δ l y:[0,T]\to\Delta^l y : [ 0 , T ] → Δ l be a measurable map with y t = y 0 + ∫ [ 0 , t ] b ( y u , a u ) d u y_t=y_0+\int_{[0,t]}b(y_u,a_u)\,du y t = y 0 + ∫ [ 0 , t ] b ( y u , a u ) d u for every t t t (the flow ), and let c ∈ R l \mathbf{c}\in\mathbb{R}^l c ∈ R l (the estimand direction ; the bold letter is unrelated to the labels c c c ). Put e t = Σ t − y t e_t=\Sigma_t-y_t e t = Σ t − y t , and let ρ u \rho_u ρ u and d u \mathsf{d}_u d u be the residual and the control-gradient discrepancy of Linearisation of a Perturbed Controlled Aggregate Flow Along a Comparison Pair: Exact Variation-of-Constants Identity and Residual Bound for x = Σ x=\Sigma x = Σ and this y y y . Assume the clock-discrepancy hypothesis : a real number w 1 ≥ 0 w_1\ge0 w 1 ≥ 0 with
(CD) ∣ C u c − C ˉ u c ∣ ≤ w 1 for every c ∈ L and u ∈ [ 0 , T ] . \textbf{(CD)}\qquad|\mathsf{C}^{c}_u-\bar{\mathsf{C}}^{c}_u|\le w_1\qquad\text{for every }c\in\mathcal{L}\text{ and }u\in[0,T]. (CD) ∣ C u c − C ˉ u c ∣ ≤ w 1 for every c ∈ L and u ∈ [ 0 , T ] .
For c ∈ L c\in\mathcal{L} c ∈ L and u ∈ [ 0 , T ] u\in[0,T] u ∈ [ 0 , T ] put H u c = Φ E ( T , u ) E u ⋆ v c ∈ R l H^{c}_u=\Phi^{\mathcal{E}}(T,u)\,\mathcal{E}^{\star}_u\,v_c\in\mathbb{R}^l H u c = Φ E ( T , u ) E u ⋆ v c ∈ R l . Finally define the entry times and the cell coefficients : for q = ( c , j ) ∈ L q=(c,j)\in\mathsf{L} q = ( c , j ) ∈ L with C ˉ T c ≥ b j c \bar{\mathsf{C}}^{c}_T\ge \mathsf{b}^{c}_j C ˉ T c ≥ b j c let τ ˉ q \bar{\tau}_q τ ˉ q be the least u ∈ [ 0 , T ] u\in[0,T] u ∈ [ 0 , T ] with C ˉ u c ≥ b j c \bar{\mathsf{C}}^{c}_u\ge \mathsf{b}^{c}_j C ˉ u c ≥ b j c (it exists by claim 1) and put α q = c ⋅ ( Φ E ( T , τ ˉ q ) v c ) \alpha_q=\mathbf{c}\cdot\bigl(\Phi^{\mathcal{E}}(T,\bar{\tau}_q)\,v_c\bigr) α q = c ⋅ ( Φ E ( T , τ ˉ q ) v c ) ; for q = ( c , j ) q=(c,j) q = ( c , j ) with C ˉ T c < b j c \bar{\mathsf{C}}^{c}_T<\mathsf{b}^{c}_j C ˉ T c < b j c put α q = 0 \alpha_q=0 α q = 0 .
1. (Perturbation form of the solution.) Each ϕ c \phi_c ϕ c is measurable with values in [ 0 , B ] [0,B] [ 0 , B ] ; Σ \Sigma Σ is bounded measurable with values in G N ⊆ Δ l \mathbb{G}_N\subseteq\Delta^l G N ⊆ Δ l , and e e e is bounded measurable with ∣ e t ∣ ≤ 2 |e_t|\le2 ∣ e t ∣ ≤ 2 ; each C c \mathsf{C}^{c} C c and each C ˉ c \bar{\mathsf{C}}^{c} C ˉ c is continuous and nondecreasing on [ 0 , T ] [0,T] [ 0 , T ] with values in [ 0 , N B T ] ⊆ [ 0 , R ] [0,NBT]\subseteq[0,R] [ 0 , NBT ] ⊆ [ 0 , R ] ; for every q = ( c , j ) q=(c,j) q = ( c , j ) with C ˉ T c ≥ b j c \bar{\mathsf{C}}^{c}_T\ge \mathsf{b}^{c}_j C ˉ T c ≥ b j c the set { u ∈ [ 0 , T ] : C ˉ u c ≥ b j c } \{u\in[0,T]:\bar{\mathsf{C}}^{c}_u\ge \mathsf{b}^{c}_j\} { u ∈ [ 0 , T ] : C ˉ u c ≥ b j c } has a least element τ ˉ q \bar{\tau}_q τ ˉ q and equals [ τ ˉ q , T ] [\bar{\tau}_q,T] [ τ ˉ q , T ] ; and, with the perturbation
m t = 1 N ∑ c ∈ L v c M ˉ c ( C t c ) , M ˉ c ( x ) = p c ( x ) − x ( x ∈ [ 0 , R ] ) , \mathfrak{m}_t=\frac{1}{N}\sum_{c\in\mathcal{L}}v_c\,\bar{M}^{c}(\mathsf{C}^{c}_t),\qquad \bar{M}^{c}(x)=p^{c}(x)-x\quad(x\in[0,R]), m t = N 1 c ∈ L ∑ v c M ˉ c ( C t c ) , M ˉ c ( x ) = p c ( x ) − x ( x ∈ [ 0 , R ]) ,
which is bounded measurable with m 0 = 0 \mathfrak{m}_0=0 m 0 = 0 , one has Σ t = x 0 + ∫ [ 0 , t ] b ( Σ u , a u ) d u + m t \Sigma_t=x_0+\int_{[0,t]}b(\Sigma_u,a_u)\,du+\mathfrak{m}_t Σ t = x 0 + ∫ [ 0 , t ] b ( Σ u , a u ) d u + m t for every t ∈ [ 0 , T ] t\in[0,T] t ∈ [ 0 , T ] . Consequently Linearisation of a Perturbed Controlled Aggregate Flow Along a Comparison Pair: Exact Variation-of-Constants Identity and Residual Bound applies to x = Σ x=\Sigma x = Σ , y y y , a a a and m \mathfrak{m} m .
2. (Linear response in cell-count form.) N c ⋅ L T ( m ) = 1 N ∑ c ∈ L c ⋅ Λ c ( C c ) \sqrt{N}\,\mathbf{c}\cdot L_T(\mathfrak{m})=\frac{1}{\sqrt{N}}\sum_{c\in\mathcal{L}}\mathbf{c}\cdot\Lambda^{c}(\mathsf{C}^{c}) N c ⋅ L T ( m ) = N 1 ∑ c ∈ L c ⋅ Λ c ( C c ) , where Λ c \Lambda^{c} Λ c is the compensated functional of Cell-Count Form of a Compensated Counting-Path Functional Along a Time Change: Window-Discrepancy Bounds for the Time-Change and Partial-Cell Errors formed from p c p^{c} p c , R R R , s = T s=T s = T , v = v c v=v_c v = v c and H = H c H=H^{c} H = H c (bounded measurable with ∥ H c ∥ 1 ≤ 2 Φ ˉ 2 Λ E T \lVert H^{c}\rVert_1\le\sqrt{2}\,\bar{\Phi}^{2}\Lambda_{\mathcal{E}}T ∥ H c ∥ 1 ≤ 2 Φ ˉ 2 Λ E T ); and
∣ N c ⋅ L T ( m ) − 1 N ∑ q ∈ L α q ( K q − μ q ) ∣ ≤ ∣ c ∣ N ∑ c ∈ L ( 2 + ∥ H c ∥ 1 ) ( D i s c w 1 ( p c ) + D i s c μ max ( p c ) ) = : E c e l l . \Bigl|\sqrt{N}\,\mathbf{c}\cdot L_T(\mathfrak{m})-\frac{1}{\sqrt{N}}\sum_{q\in\mathsf{L}}\alpha_q\,(\mathsf{K}_q-\mu_q)\Bigr|\ \le\ \frac{|\mathbf{c}|}{\sqrt{N}}\sum_{c\in\mathcal{L}}\bigl(\sqrt{2}+\lVert H^{c}\rVert_1\bigr)\bigl(\mathrm{Disc}_{w_1}(p^{c})+\mathrm{Disc}_{\mu_{\max}}(p^{c})\bigr)=:\mathsf{E}_{\mathrm{cell}} . N c ⋅ L T ( m ) − N 1 q ∈ L ∑ α q ( K q − μ q ) ≤ N ∣ c ∣ c ∈ L ∑ ( 2 + ∥ H c ∥ 1 ) ( Disc w 1 ( p c ) + Disc μ m a x ( p c ) ) =: E cell .
3. (Estimand linearisation.) Writing e ˉ = sup u ∈ [ 0 , T ] N ∣ e u ∣ \bar{e}=\sup_{u\in[0,T]}\sqrt{N}\,|e_u| e ˉ = sup u ∈ [ 0 , T ] N ∣ e u ∣ ,
∣ N c ⋅ ( Σ T − y T ) − 1 N ∑ q ∈ L α q ( K q − μ q ) ∣ ≤ ∣ c ∣ Φ ˉ 2 N ∣ x 0 − y 0 ∣ + ∣ c ∣ Φ ˉ 2 E r e s + E c e l l , \Bigl|\sqrt{N}\,\mathbf{c}\cdot(\Sigma_T-y_T)-\frac{1}{\sqrt{N}}\sum_{q\in\mathsf{L}}\alpha_q\,(\mathsf{K}_q-\mu_q)\Bigr|\ \le\ |\mathbf{c}|\,\bar{\Phi}^{2}\,\sqrt{N}\,|x_0-y_0|+|\mathbf{c}|\,\bar{\Phi}^{2}\,\mathsf{E}_{\mathrm{res}}+\mathsf{E}_{\mathrm{cell}}, N c ⋅ ( Σ T − y T ) − N 1 q ∈ L ∑ α q ( K q − μ q ) ≤ ∣ c ∣ Φ ˉ 2 N ∣ x 0 − y 0 ∣ + ∣ c ∣ Φ ˉ 2 E res + E cell ,
where
E r e s = N ∫ [ 0 , T ] ∣ ρ u ∣ d u ≤ 2 l ( l − 1 ) Λ 2 T e ˉ 2 N + 2 l ( l − 1 ) Λ 3 e ˉ ∫ [ 0 , T ] ∣ y u − S u ∣ d u + 2 e ˉ ∫ [ 0 , T ] d u d u , \mathsf{E}_{\mathrm{res}}=\sqrt{N}\int_{[0,T]}|\rho_u|\,du\ \le\ \sqrt{2}\,l(l-1)\Lambda_2T\,\frac{\bar{e}^{2}}{\sqrt{N}}+\sqrt{2}\,l(l-1)\Lambda_3\,\bar{e}\int_{[0,T]}|y_u-S_u|\,du+\sqrt{2}\,\bar{e}\int_{[0,T]}\mathsf{d}_u\,du , E res = N ∫ [ 0 , T ] ∣ ρ u ∣ d u ≤ 2 l ( l − 1 ) Λ 2 T N e ˉ 2 + 2 l ( l − 1 ) Λ 3 e ˉ ∫ [ 0 , T ] ∣ y u − S u ∣ d u + 2 e ˉ ∫ [ 0 , T ] d u d u ,
the two last integrands being bounded measurable.