---
title: 'Caputo Derivative: Definition & Applications'
url: https://www.emergentmind.com/topics/caputo-derivative
type: topic
---

# Caputo Derivative: Definition & Applications

The Caputo derivative is a fractional calculus operator that extends the concept of integer-order differentiation to non-integer (fractional) orders, while preserving compatibility with classical initial and boundary value problems. It plays a fundamental role in the analysis and modeling of systems exhibiting memory and hereditary properties, particularly within the frameworks of anomalous transport, viscoelasticity, fractional diffusion, and fractional variational calculus.

## 1. Definition and Fundamental Formulation

For a function \( f \) defined on the interval \([a, b]\) and sufficiently smooth (e.g., \( f \in C^n([a,b]) \)), the left Caputo fractional derivative of order \(\alpha\) (\(0 < \alpha < 1\)) is defined as
\[
{}^C_{a}D_{t}^\alpha f(t) = \frac{1}{\Gamma(1-\alpha)} \int_a^t (t-\tau)^{-\alpha} f'(\tau)\,d\tau,
\]
where \(\Gamma\) denotes the Gamma function. The right Caputo derivative of order \(\beta\) is defined analogously:
\[
{}^C_{t}D_{b}^\beta f(t) = -\frac{1}{\Gamma(1-\beta)} \int_t^b (\tau-t)^{-\beta} f'(\tau)\,d\tau.
\]
These definitions generalize to arbitrary order \( \alpha > 0 \), with
\[
{}^C_{a}D_{t}^\alpha f(t) = \frac{1}{\Gamma(n-\alpha)} \int_a^t (t-\tau)^{n-\alpha-1} f^{(n)}(\tau)\,d\tau, \qquad n-1 < \alpha < n.
\]
The Caputo derivative differs structurally from the Riemann–Liouville derivative in that the integer derivative appears inside the convolution kernel, yielding compatibility with standard (classical) initial conditions and ensuring that the fractional derivative of a constant vanishes [1109.0658, 1007.2937].

## 2. Key Analytical Properties

- **Linearity:** For scalars \( c_1, c_2 \) and appropriately smooth functions \( f, g \),
  \[
  {}^C_{a}D_{t}^\alpha(c_1 f + c_2 g)(t) = c_1\,{}^C_{a}D_{t}^\alpha f(t) + c_2\,{}^C_{a}D_{t}^\alpha g(t).
  \]
- **Vanishing on constants:** \( {}^C_{a}D_{t}^\alpha c = 0 \) for constant \(c\).
- **Classical limit and compatibility:** As \( \alpha \to 1^- \),
  \[
  \lim_{\alpha \to 1^-} {}^C_{a}D_{t}^\alpha f(t) = f'(t).
  \]
  For integer order \(n\), the Caputo derivative reduces to the ordinary \(n\)th derivative.
- **Relation to Riemann–Liouville derivative:**
  \[
  {}^C_{a}D_{t}^\alpha f(t) = {}_{a}D_{t}^\alpha f(t) - \frac{f(a)}{\Gamma(1-\alpha)}(t-a)^{-\alpha}.
  \]
- **Fractional integration by parts:** For \(0 < \alpha < 1\), the Caputo operator satisfies
  \[
  \int_a^b g(t)\,{}^C_{a}D_{t}^\alpha f(t)\,dt = \left[ f(t) I_b^{1-\alpha}g(t) \right]_{t=a}^{t=b} + \int_a^b f(t) {}_{t}D_{b}^\alpha g(t)\,dt,
  \]
  where \(I_b^{1-\alpha}\) is the right-sided Riemann–Liouville fractional integral [1007.2937, 1109.0658, 1002.3790].

## 3. Combined Caputo Derivative and Extensions

Malinowska & Torres introduced the combined Caputo derivative as a convex combination of the left and right Caputo operators:
\[
{}^C D^{\alpha,\beta}_\gamma f(t) = \gamma\,{}^C_{a}D_{t}^\alpha f(t) + (1-\gamma)\,{}^C_{t}D_{b}^\beta f(t),
\]
for \(\gamma \in [0,1]\), \(0 < \alpha, \beta < 1\). Special cases recover standard left (\(\gamma=1\)) and right (\(\gamma=0\)) Caputo derivatives, and the symmetric (Klimek) derivative (\(\alpha=\beta\), \(\gamma=1/2\)). This operator is linear in \(f\), admits a suitable integration-by-parts formula, and generalizes variational principles built on the Caputo model [1109.4664, 1007.0743].

## 4. Fractional Variational Calculus and Euler–Lagrange Equations

In the Caputo calculus, the fractional variational principle is formulated for functionals of the form
\[
J[y] = \int_a^b L(t, y(t), {}^C_{a}D_{t}^\alpha y(t), {}^C_{t}D_{b}^\beta y(t))\,dt,
\]
where \(L\) may depend on left and right Caputo derivatives of \(y\). The associated Euler–Lagrange equation is
\[
\frac{\partial L}{\partial y} + {}_{t}D_{b}^\alpha\left( \frac{\partial L}{\partial ({}^C_{a}D_{t}^\alpha y)} \right) + {}_{a}D_{t}^\beta\left( \frac{\partial L}{\partial ({}^C_{t}D_{b}^\beta y)} \right) = 0,
\]
with analogous forms in the combined Caputo framework [1109.4664, 1007.0743, 1109.0658, 1210.0705].

The boundary conditions—so-called transversality conditions—require careful handling, as the presence of fractional derivatives naturally introduces nonlocal contributions, realized through fractional integrals of the Lagrangian derivatives at the interval endpoints [1002.3790, 1109.4664].

For isoperimetric problems, additional constraints are incorporated via augmented Lagrangians, with Lagrange multipliers and complementary slackness conditions as in the classical theory [1109.4664, 1007.0743, 1007.2937].

## 5. Numerical Discretization and Computational Algorithms

The Caputo derivative's nonlocal structure leads to quasi-memory effects in numerical implementations. The L1 scheme, the shifted Grünwald–Letnikov method, and sum-of-exponentials approximations are standard approaches to discretization. For a uniform mesh \(t_n = n h\), the L1 approximation reads
\[
D_{C,h}^\alpha f_n = \frac{1}{\Gamma(2-\alpha)\,h^\alpha} \sum_{k=0}^{n-1} \left[ (n-k)^{1-\alpha} - (n-k-1)^{1-\alpha} \right] (f_{k+1}-f_k).
\]
Recent results quantify truncation errors and stability for both fixed and small fractional orders \(\alpha\), and establish unconditional convergence for implicit schemes applied to fractional PDEs [2107.10489, 2108.02999, 1806.03421, 1901.03159].

In time-fractional PDE solvers, efficient history evaluation, such as the Fast Implicit Double Recurrence (FIDR) algorithm, is essential to mitigate computational cost, especially in the small \(\alpha\) regime. Discrete convolution coefficients must be carefully analyzed for positivity and monotonicity to guarantee stability and equipartition of error [2108.02999, 1901.03159].

## 6. Generalizations: Variable and Distributed Order Caputo Derivatives

Extensions to variable-order and distributed-order settings have been developed to address systems where the fractional order itself may vary in space, time, or be governed by a probability distribution. The variable order Caputo derivatives take the form
\[
{}^{C}_aD_t^{\alpha(t)} f(t) = \frac{1}{\Gamma(1-\alpha(t))} \int_a^t (t-\tau)^{-\alpha(t)} f'(\tau) d\tau,
\]
with specific expansions and error control for numerical schemes [1511.02017].

The distributed-order Caputo derivative is formalized as
\[
D_t^{(\mu)} u(t) = \int_0^1 \mu(\alpha) D_t^\alpha u(t) d\alpha,
\]
where \(\mu(\alpha)\) is a nonnegative integrable weight, with existence and regularity results for corresponding fractional diffusion equations derived under minimal regularity of data and coefficients [1706.05591].

## 7. Analytical and Applied Implications

The Caputo derivative provides an operator-theoretic and functional-analytic framework for nonlocal models, with spectral characterizations linked to fractional powers of Volterra and elliptic operators in appropriate Sobolev spaces [1411.7289]. It is suitable both for modeling physical processes exhibiting nonlocality or memory and for rigorous mathematical analysis, due to its compatibility with classical data and well-posedness structures.

Monotonicity results for Caputo derivatives extend classical criteria: monotonicity of a function on an interval is characterized by the sign of all its Caputo derivatives of order \(\alpha \in (0,1)\) on that interval [1510.03285]. The Caputo calculus also enables generalizations of the fundamental lemma of the calculus of variations and derivation of Euler–Lagrange equations involving only Caputo derivatives, thus eliminating the need for Riemann–Liouville terms in boundary conditions [1210.0705].

In computational contexts, the Caputo fractional derivative informs the design of stable, efficient iterative optimization algorithms, such as Caputo-fractional gradient descent, which exhibit smoothing properties and close links to Tikhonov regularization in quadratic cases [2104.02259].

The Caputo derivative's analytical tractability, compatibility with physical initial and boundary data, and adaptability to numerical schemes ensure its centrality in the theory and applications of fractional differential equations and variational problems across mathematical physics, engineering, and applied analysis [1109.4664, 1007.2937, 1002.3790, 1109.0658, 2107.10489, 1411.7289, 2108.02999, 1510.03285, 1706.05591, 1511.02017, 1210.0705, 1806.03421, 1901.03159, 1711.07126].

Source: https://www.emergentmind.com/topics/caputo-derivative