Tutorials › Population Pharmacokinetics › Infusion Start, Stop, And Washout

Choosing the Structural Model · Tutorial 97 of 599

Infusion Start, Stop, And Washout

Learn to construct piecewise infusion profiles, calculate concentrations after stopping an infusion, and determine how long washout takes to a specified threshold.

Intermediate 10 min read

What You'll Learn

  • State the one-compartment assumptions used for infusion start and stop calculations.
  • Derive the concentration response during infusion from the amount-balance equation.
  • Construct a continuous piecewise concentration profile at an infusion stop.
  • Calculate post-infusion concentrations using elapsed time since stopping.
  • Find time to a specified washout threshold and interpret the result.
  • Recognize common timing, continuity, and model-assumption errors.

Why infusion start and stop times matter

A constant-rate infusion does not make concentration jump to its eventual plateau. Concentration rises over time as drug enters the body while elimination removes drug. When the infusion stops, the input disappears, but drug already in the body remains and continues to be eliminated. The concentration therefore changes its pattern at the stop time: it rises toward a plateau during infusion and declines afterward.

The calculation is naturally piecewise: use one equation while drug is being infused and another after the input ends. The two equations must agree at the stop time. This tutorial derives those expressions for the one-compartment model and uses them to calculate post-infusion decline and washout.

Definition: Let \(C(t)\) be the drug concentration at time \(t\), \(V\) the apparent volume of distribution, \(CL\) the clearance, and \(R_0\) a constant infusion rate. In the one-compartment model with first-order elimination, \(k=CL/V\) is the elimination rate constant. During infusion, \(dA/dt=R_0-CL\,C(t)\), where \(A(t)=V C(t)\). After the infusion stops, the input rate is zero.

Concentration while an infusion is running

Assume the infusion begins at time \(t_s\), continues at the constant rate \(R_0\), and the concentration at its start is \(C_s\). The model assumes instantaneous mixing in one compartment, constant \(V\) and \(CL\), linear first-order elimination, and a known infusion rate. There are no additional doses or changes in these parameters during the interval.

Since \(A=VC\), the amount-balance equation during infusion becomes \(V\,dC/dt=R_0-CL\,C\). Dividing by \(V\) gives \(dC/dt+kC=R_0/V\). The constant-input equilibrium concentration is \(R_0/CL\), because at equilibrium input equals elimination: \(R_0=CL\,C\). Solving the differential equation from the starting concentration \(C_s\) gives:

$$ C(t)=C_s e^{-k(t-t_s)}+\frac{R_0}{CL}\left(1-e^{-k(t-t_s)}\right), \qquad t_s\leq t\leq t_e. $$

Here \(t_e\) is the stop time. The elapsed infusion time is \(t-t_s\), not necessarily the clock time \(t\). If the compartment is drug-free at the start, \(C_s=0\), and the equation reduces to \(C(t)=(R_0/CL)(1-e^{-k(t-t_s)})\). The concentration approaches \(R_0/CL\) asymptotically; it does not reach that value at a finite time under this model.

Worked Example: Concentration After Starting an Infusion

Suppose a hypothetical patient has \(V=20\) L and \(CL=2.0\) L/h. An infusion begins at time zero at \(R_0=8.0\) mg/h, and the initial concentration is zero. Find the concentration after 3.0 h.

First, \(k=CL/V=(2.0\ \text{L/h})/(20\ \text{L})=0.1000\ \text{h}^{-1}\). The constant-input equilibrium concentration is \(R_0/CL=(8.0\ \text{mg/h})/(2.0\ \text{L/h})=4.0\) mg/L. Substitution into the start-of-infusion equation gives:

$$ C(3.0)=4.0\left(1-e^{-(0.1000)(3.0)}\right) =4.0(1-0.7408) =1.037\ \text{mg/L}. $$

Thus the predicted concentration is 1.037 mg/L after 3.0 h. The units are consistent: \(k\) has units of inverse time, so its product with elapsed time is dimensionless, and \(R_0/CL\) has units of concentration.

Stopping the infusion and joining the pieces

At \(t_e\), the infusion rate changes from \(R_0\) to zero. The model contains no instantaneous removal of drug at that moment, so concentration is continuous: the value just after stopping equals the value just before stopping. Denote this shared concentration by \(C_e=C(t_e)\). After the stop, the amount-balance equation is \(dA/dt=-CL\,C\), or \(dC/dt=-kC\). Its solution, initialized at \(C_e\), is:

$$ C(t)=C_e e^{-k(t-t_e)},\qquad t\geq t_e. $$

A complete profile for a single infusion is therefore obtained by calculating \(C_e\) from the infusion equation and using it as the initial condition for the post-infusion equation:

$$ C(t)= \begin{cases} C_s e^{-k(t-t_s)}+\dfrac{R_0}{CL}\left(1-e^{-k(t-t_s)}\right), & t_s\leq t\leq t_e,\\[6pt] C_e e^{-k(t-t_e)}, & t\geq t_e, \end{cases} \qquad C_e=C_s e^{-k(t_e-t_s)}+\dfrac{R_0}{CL}\left(1-e^{-k(t_e-t_s)}\right). $$

The equations agree at the boundary: setting \(t=t_e\) in the second piece gives \(C_e e^0=C_e\), exactly the value calculated from the first piece. The slope can change at the stop because the input changes, even though concentration itself does not jump.

Worked Example: Concentration During and After an Infusion

For a hypothetical one-compartment model, let \(V=20\) L, \(CL=3.0\) L/h, and \(R_0=12\) mg/h. The infusion starts at zero concentration and stops after 6.0 h. Find the concentration at stopping and 4.0 h later.

The elimination rate constant is \(k=3.0/20=0.1500\ \text{h}^{-1}\), and the equilibrium concentration during infusion is \(R_0/CL=12/3.0=4.0\) mg/L. At the stop time:

$$ C_e=4.0\left(1-e^{-(0.1500)(6.0)}\right) =4.0(1-0.4066) =2.374\ \text{mg/L}. $$

Four hours after stopping, the elapsed time in the decline equation is \(4.0\) h:

$$ C(t_e+4.0)=2.3737e^{-(0.1500)(4.0)} =2.3737(0.5488) =1.303\ \text{mg/L}. $$

The concentration is 2.374 mg/L when the infusion stops and 1.303 mg/L four hours later. The post-stop calculation uses the concentration at the stop as its starting value; it does not continue the infusion equation after the input has ended.

Washout to a specified concentration

Under first-order elimination, post-infusion concentration declines exponentially. For any positive threshold \(C_{\text{th}}\) below the concentration at stopping, solve \(C_{\text{th}}=C_e e^{-k\Delta t}\) for the time \(\Delta t\) after stopping:

$$ \Delta t=\frac{\ln(C_e/C_{\text{th}})}{k}. $$

This result requires \(C_e>C_{\text{th}}>0\). If the threshold is defined as a fraction \(f\) of the stop concentration, where \(0<f<1\), then \(C_{\text{th}}=fC_e\), so \(\Delta t=-\ln(f)/k\). For example, the time for concentration to fall to one half of its stop value is \(\ln(2)/k\), the elimination half-life. A positive exponential approaches zero indefinitely; in this mathematical model, it never becomes exactly zero at a finite time. “Washout” therefore needs a stated threshold or fraction to have a definite time.

Worked Example: Time to a Washout Threshold

Use the stop concentration from the preceding example, \(C_e=2.3737\) mg/L, with \(k=0.1500\ \text{h}^{-1}\). How long after stopping does concentration take to fall below 0.100 mg/L?

The threshold is positive and below \(C_e\), so the threshold-time equation applies. The time at which concentration reaches the threshold is:

$$ \Delta t=\frac{\ln(2.3737/0.100)}{0.1500} =\frac{\ln(23.737)}{0.1500} \approx\frac{3.167}{0.1500} =21.11\ \text{h}. $$

The predicted concentration reaches 0.100 mg/L about 21.11 h after the infusion stops; it is below that threshold at later times. This is a model-based washout time, not a claim that drug is completely absent from the body.

Residual concentration and interrupted inputs

An infusion need not begin from a drug-free state. If \(C_s\) is already present at the start, the starting concentration contributes \(C_s e^{-k(t-t_s)}\), while the infusion contributes the remaining term. This follows directly from the solution above and prevents the common error of discarding residual drug when a new infusion begins.

Worked Example: Starting an Infusion With Residual Drug

Suppose a hypothetical patient has \(V=25\) L and \(CL=2.5\) L/h. A new infusion of 10 mg/h begins when the concentration is 1.20 mg/L. Find the concentration after 5.0 h, then 3.0 h after the infusion stops.

Here \(k=2.5/25=0.1000\ \text{h}^{-1}\), and \(R_0/CL=10/2.5=4.0\) mg/L. At 5.0 h, while the infusion is running:

$$ C(5.0)=1.20e^{-0.500}+4.0(1-e^{-0.500}) =1.20(0.6065)+4.0(0.3935) =0.7278+1.574 =2.302\ \text{mg/L}. $$

This value is the concentration at stopping. Three hours later:

$$ C(t_e+3.0)=2.3017e^{-(0.1000)(3.0)} =2.3017(0.7408) =1.705\ \text{mg/L}. $$

The residual concentration contributes during the infusion, so the stop concentration is greater than it would be if the infusion had started from zero. After stopping, the entire concentration declines from 2.302 mg/L according to the post-infusion equation.

Checks, comparisons, and limitations

A useful way to check a piecewise calculation is to verify the time origin, boundary value, direction of change, and units. During infusion from a concentration below \(R_0/CL\), concentration rises toward that equilibrium; after stopping, a positive concentration declines. At the stop, both pieces must give the same concentration. The exponent must be dimensionless, and all concentrations compared or added must use the same units.

The infusion profile differs from an IV bolus profile because infusion supplies drug over an interval rather than instantaneously. The constant-rate infusion equation describes the rising phase; after the input stops, the same first-order elimination process gives a simple exponential decline. A bolus-style decline equation is appropriate only after the relevant input has ended, and its starting concentration must be the concentration at that time. More complicated input schedules can be handled by applying the same logic over each interval with a defined input rate and carrying the boundary concentration forward.

These expressions are only as appropriate as their assumptions. A changing infusion rate requires a different input function or separate intervals; nonlinear elimination, changing clearance or volume, incomplete mixing, or an incorrect record of start and stop times can make the predicted curve misleading. Evidence of systematic disagreement between predictions and observed concentrations around the start or stop may indicate a timing problem, a misspecified input history, or inadequacy of the one-compartment assumptions. The equations alone cannot distinguish among those explanations.

Common mistake: Using the total time since the infusion began as the elapsed time in the post-infusion equation, or continuing to use the infusion equation after stopping. The correction is to calculate the stop concentration first, then use time since the stop: \(C(t)=C_e e^{-k(t-t_e)}\). Also, washout is not a finite time to exactly zero; report the threshold or fraction used.

Key takeaway: Calculate concentration during the infusion using elapsed time since its start. At the stop, carry the resulting concentration forward unchanged, then calculate decline using elapsed time since stopping. For a washout time, specify a positive threshold or remaining fraction.

Check Your Understanding

Try each problem, then compare your work with the solutions below. The data are hypothetical.

  1. An infusion begins at zero concentration with \(V=32\) L, \(CL=4.0\) L/h, and \(R_0=16\) mg/h. What is the predicted concentration after 2.0 h?
  2. An infusion with \(V=10\) L, \(CL=2.0\) L/h, and \(R_0=6.0\) mg/h starts from zero and stops after 4.0 h. Find the stop concentration and the concentration 2.0 h later.
  3. If \(k=0.1000\ \text{h}^{-1}\), how long after stopping does concentration take to fall to 10% of its stop value?
  4. Why is it incorrect to say that the model predicts complete washout at a finite time?

Solutions. (1) \(k=4.0/32=0.1250\ \text{h}^{-1}\) and \(R_0/CL=4.0\) mg/L. Thus \(C(2.0)=4.0(1-e^{-0.250})=4.0(1-0.7788)=0.8848\) mg/L. (2) \(k=2.0/10=0.2000\ \text{h}^{-1}\) and \(R_0/CL=3.0\) mg/L. At stopping, \(C_e=3.0(1-e^{-0.800})=3.0(1-0.4493)=1.652\) mg/L. Two hours later, \(C=1.652e^{-0.400}=1.107\) mg/L. (3) For \(f=0.10\), \(\Delta t=-\ln(0.10)/0.1000=23.03\) h. (4) After stopping, \(C=C_e e^{-k\Delta t}\), which remains positive for every finite \(\Delta t\) when \(C_e>0\); it approaches zero asymptotically. A finite washout time therefore requires a defined positive threshold.

Clinical Biostats