Animated Solution for Physics - Atoms and Nuclei: At time t=0, a material is composed of two radioactive atoms A and B, where NA(0)=2NB(0). The decay constant of both kind of radioactive atoms is λ. However, A disintegrates to B and B disintegrates to C. Which of the following figures represents the evolution of NB(t)/NB(0) with respect to time t?
[NA(0)=Number of A atoms at t=0NB(0)=Number of B atoms at t=0]
Select Answer:
Visualized Solution
Initial Setup
\text{Decay Chain: } A \xrightarrow{\lambda} B \xrightarrow{\lambda} C
\text{Initial Conditions: } N_A(0) = 2N_B(0)
Rate Equation
\text{Rate of change of } B:
\frac{dN_B(t)}{dt} = \text{Rate of formation} - \text{Rate of decay}
The problem of successive radioactive decays is a beautiful intersection of nuclear physics and calculus. When a parent nucleus decays into a radioactive daughter nucleus, the population of the daughter nucleus is governed by a delicate balance between its rate of formation and its rate of decay. Let's dive into the mathematics of this relay race!
Analyzing the Setup
We are given a decay chain where atom A disintegrates into B, and B further disintegrates into C. Both A and B share the exact same decay constant, λ.
Initially, at t=0, the sample contains twice as many A atoms as B atoms. Mathematically, this is written as:
NA(0)=2NB(0)
Our goal is to find how the ratio NB(0)NB(t) evolves over time. To do this, we must construct the rate equation for NB(t).
The Master Equation
The number of B atoms changes due to two simultaneous processes:
1. Formation:B is created when A decays. The rate of formation is λNA(t).
2. Decay:B is destroyed when it decays into C. The rate of decay is λNB(t).
Combining these, the net rate of change of B is:
dtdNB(t)=λNA(t)−λNB(t)
We already know that A undergoes simple exponential decay, so NA(t)=NA(0)e−λt. Substituting this into our rate equation gives:
dtdNB(t)=λNA(0)e−λt−λNB(t)
The Integrating Factor Magic
This is a classic first-order linear differential equation. Let's rearrange it into the standard form:
dtdNB(t)+λNB(t)=λNA(0)e−λt
To solve this, we multiply the entire equation by an integrating factor, which is e∫λdt=eλt. This brilliant mathematical trick collapses the left side into the exact derivative of a product:
dtd(NB(t)eλt)=λNA(0)
Now, we substitute our initial condition NA(0)=2NB(0):
dtd(NB(t)eλt)=2λNB(0)
Integrating both sides with respect to time t:
NB(t)eλt=2λNB(0)t+C
To find the constant of integration C, we apply the initial condition at t=0:
NB(0)⋅e0=0+C⟹C=NB(0)
Final Calculation and Graphing
Substituting C back into our equation and dividing by eλt, we get the explicit function for NB(t):
NB(t)=NB(0)(1+2λt)e−λt
Dividing by NB(0) gives the exact ratio we need to plot:
NB(0)NB(t)=(1+2λt)e−λt
To determine the shape of the graph, we need to find if and where this ratio reaches a maximum. We differentiate the expression with respect to t and set it to zero:
dtd[(1+2λt)e−λt]=2λe−λt−λ(1+2λt)e−λt=0
Factoring out λe−λt:
λe−λt(2−1−2λt)=0
1−2λt=0⟹t=2λ1
At this peak time, the value of the ratio is:
NB(0)NB=(1+2λ⋅2λ1)e−1/2=e2≈1.21
Since 1.21>1, the graph starts at 1, rises to a peak at t=2λ1, and then decays exponentially towards zero. This perfectly matches the visual representation in Option (c).