A Sensitivity-based Framework for Calibrating Coupled Ordinary Differential Equations under Model Discrepancy
Abstract
Coupled systems of ordinary differential equations (ODEs) are widely used to model complex physical and biological processes. In these applications, ODE model parameters must be calibrated to field data to enable prediction and parameter inference. In practice, the governing ODEs are an imperfect approximation to the true system, leading to non-negligible model discrepancy that must be incorporated to avoid biased parameter estimates. However, it is well known that highly flexible discrepancy models can confound discrepancy and simulator parameters, leading to poor identifiability. As a result, practitioners often must impose problem-specific prior constraints to adequately regularize the discrepancy, which can be challenging and cumbersome. In this work, we propose a novel calibration framework for coupled, multi-output ODE systems that enforces automatic, model-structural constraints on the discrepancy, improving identifiability without requiring application-specific discrepancy priors. Forward-model gradients are computed via sensitivity equations, solved jointly with the ODE system as a coupled initial value problem. Posterior inference is performed using an adaptive Metropolis-within-Gibbs sampler tailored to the resulting constrained posterior. We demonstrate the approach on three model systems: a mass-spring oscillator, a model of infectious disease spread, and a neuron firing model. We show that our method leads to improved parameter inference and predictive accuracy compared to standard calibration approaches.
Related
- Derivative-Informed Dimension Reduction of Stochastic Differential Equations: A Path-Space Sensitivity Approach
- A Unified Space-Time Variational Framework for Wave Equations on Differential Forms
- QGPINNs: A Physics-Informed Neural Network Framework for Nonlocal Differential Equations on Quantum Graphs
- An intrinsic finite element framework for scalar- and vector-valued partial differential equations on general manifolds
- Synthetic geometry of differential equations: I. Jets and comonad structure
BibTeX
@misc{colebank2026sensitivitybased,
title = {{A Sensitivity-based Framework for Calibrating Coupled Ordinary Differential Equations under Model Discrepancy}},
author = {Mitchel J. Colebank and William Consagra},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.09261},
eprint = {2610.09261},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.09261},
note = {CC BY-NC-SA 4.0},
}
RIS
TY - UNPB
TI - A Sensitivity-based Framework for Calibrating Coupled Ordinary Differential Equations under Model Discrepancy
AU - Mitchel J. Colebank
AU - William Consagra
PY - 2026
DA - 2026/10/07
PB - arXiv
DO - 10.48550/arxiv.2610.09261
UR - https://arxiv.org/abs/2610.09261
AB - Coupled systems of ordinary differential equations (ODEs) are widely used to model complex physical and biological processes. In these applications, ODE model parameters must be calibrated to field data to enable prediction and parameter inference. In practice, the governing ODEs are an imperfect approximation to the true system, leading to non-negligible model discrepancy that must be incorporated to avoid biased parameter estimates. However, it is well known that highly flexible discrepancy models can confound discrepancy and simulator parameters, leading to poor identifiability. As a result, practitioners often must impose problem-specific prior constraints to adequately regularize the discrepancy, which can be challenging and cumbersome. In this work, we propose a novel calibration framework for coupled, multi-output ODE systems that enforces automatic, model-structural constraints on the discrepancy, improving identifiability without requiring application-specific discrepancy priors. Forward-model gradients are computed via sensitivity equations, solved jointly with the ODE system as a coupled initial value problem. Posterior inference is performed using an adaptive Metropolis-within-Gibbs sampler tailored to the resulting constrained posterior. We demonstrate the approach on three model systems: a mass-spring oscillator, a model of infectious disease spread, and a neuron firing model. We show that our method leads to improved parameter inference and predictive accuracy compared to standard calibration approaches.
ER -