Boltzmann–Matano analysis

The Boltzmann-Matano method is used to convert the partial differential equation resulting from Fick's law of diffusion into a more easily solved ordinary differential equation, which can then be applied to calculate the diffusion coefficient as a function of concentration.

Ludwig Boltzmann worked on Fick's second law to convert it into an ordinary differential equation, whereas Chujiro Matano performed experiments with diffusion couples and calculated the diffusion coefficients as a function of concentration in metal alloys. Specifically, Matano proved that the diffusion rate of A atoms into a B atom crystal lattice is a function of the amount of A atoms already in the B lattice.

The importance of the classic Boltzmann-Matano method consists in the ability to extract diffusivities from concentration-distance data. These methods, also known as inverse methods, have both proven to be reliable, convenient and accurate with the assistance of modern computational techniques.

Boltzmann’s Transformation

Boltzmann’s Transformation converts Fick's second law into an easily solvable ordinary differential equation. Assuming a diffusion coefficient D that is in general a function of concentration c, Fick's second law is:

${\frac {\partial c}{\partial t}}={\frac {\partial }{\partial x}}\overbrace {\left[D(c){\frac {\partial c}{\partial x}}\right]} ^{\text{Flux}}$

where t is time and x is distance.

Boltzmann's transformation consists in introducing a variable ξ, defined as a combination of t and x:

$\xi ={\frac {x}{2{\sqrt {t}}}}$

The partial derivatives of ξ are:

${\frac {\partial \xi }{\partial t}}=-{\frac {x}{4t^{3/2}}}=-{\frac {\xi }{2t}}$

${\frac {\partial \xi }{\partial x}}={\frac {1}{2{\sqrt {t}}}}$

To introduce ξ into Fick's law, we express its partial derivatives in terms of ξ, using the chain rule:

${\frac {\partial c}{\partial t}}={\frac {\partial c}{\partial \xi }}{\frac {\partial \xi }{\partial t}}=-{\frac {\xi }{2t}}{\frac {\partial c}{\partial \xi }}$

${\frac {\partial c}{\partial x}}={\frac {\partial c}{\partial \xi }}{\frac {\partial \xi }{\partial x}}={\frac {1}{2{\sqrt {t}}}}{\frac {\partial c}{\partial \xi }}$

Inserting these expressions into Fick's law produces the following modified form:

$-{\frac {\xi }{2t}}{\frac {\partial c}{\partial \xi }}={\frac {1}{2{\sqrt {t}}}}{\frac {\partial }{\partial x}}\left[D(c){\frac {\partial c}{\partial \xi }}\right]$

Note how the time variable in the right-hand side could be taken outside of the partial derivative, since the latter regards only variable x.

It is now possible to remove the last reference to x by using again the same chain rule used above to obtain ∂ξ/∂x:

$-{\frac {\xi }{2t}}{\frac {\partial c}{\partial \xi }}={\frac {1}{4t}}{\frac {\partial }{\partial \xi }}\left[D(c){\frac {\partial c}{\partial \xi }}\right]$

Because of the appropriate choice in the definition of ξ, the time variable t can now also be eliminated, leaving ξ as the only variable in the equation, which is now an ordinary differential equation:

$-2\xi {\frac {\mathrm {d} c}{\mathrm {d} \xi }}={\frac {\mathrm {d} }{\mathrm {d} \xi }}\left[D(c){\frac {\mathrm {d} c}{\mathrm {d} \xi }}\right]$

This form is significantly easier to solve numerically, and one only needs to perform a back-substitution of t or x into the definition of ξ to find the value of the other variable.

The Parabolic Law

Observing the previous equation, a trivial solution is found for the case dc/dξ=0, that is when concentration is constant over ξ. This can be interpreted as the rate of advancement of a concentration front being proportional to the square root of time ($x\propto {\sqrt {t}}$), or, equivalently, to the time necessary for a concentration front to arrive at a certain position being proportional to the square of the distance ($t\propto x^{2}$); the square term gives the name parabolic law.

Matano’s method

Chuijiro Matano applied Boltzmann's transformation to obtain a method to calculate diffusion coefficients as a function of concentration in metal alloys. Two alloys with different concentration would be put into contact, and annealed at a given temperature for a given time t, typically several hours; the sample is then cooled to ambient temperature, and the concentration profile is virtually "frozen". The concentration profile c at time t can then be extracted as a function of the x coordinate.

In Matano's notation, the two concentrations are indicated as c_L and c_R (L and R for left and right, as shown in most diagrams), with the implicit assumption that c_L > c_R; this is however not strictly necessary as the formulas hold also if c_R is the larger one. The initial conditions are:

$c=c_{L}\qquad \forall x<0$

$c=c_{R}\qquad \forall x>0$

Also, the alloys on both sides are assumed to stretch to infinity, which means in practice that they are large enough that the concentration at their other ends is unaffected by the transient for the entire duration of the experiment.

To extract D from Boltzmann's formulation above, we integrate it from ξ=+∞, where c=c_R at all times, to a generic ξ^*; we can immediately simplify dξ, and with a change of variables we get:

$-2\int _{c_{R}}^{c^{*}}\xi \mathrm {d} c=\int _{c=c_{R}}^{c=c^{*}}\mathrm {d} \left[D(c)\left({\frac {\mathrm {d} c}{\mathrm {d} \xi }}\right)\right]$

We can translate ξ back into its definition and bring the t terms out of the integrals, as t is constant and given as the time of annealing in the Matano method; on the right-hand side, extraction from the integral is trivial and follows from definition.

$-{\frac {1}{2t}}\int _{c_{R}}^{c^{*}}x\mathrm {d} c=\left[D(c)\left({\frac {\mathrm {d} c}{\mathrm {d} x}}\right)\right]_{c=c_{R}}^{c=c^{*}}$

We know that dc/dx → 0 as c → c_R, that is the concentration curve "flattens out" when approaching the limit concentration value. We can then rearrange:

$D(c^{*})=-{\frac {1}{2t}}{\frac {\int _{c_{R}}^{c^{*}}x\mathrm {d} c}{(\mathrm {d} c/\mathrm {d} x)_{c=c^{*}}}}$

Knowing the concentration profile c(x) at annealing time t, and assuming it is invertible as x(c), we can then calculate the diffusion coefficient for all concentrations between c_R and c_L.

The Matano interface

The last formula has one significant shortcoming: no information is given about the reference according to which x should be measured. It was not necessary to introduce one as Boltzmann's transformation worked fine without a specific reference for x; it is easy to verify that the Boltzmann transformation holds also when using x-X_M instead of plain x.

X_M is often indicated as the Matano interface, and is in general not coincident with x=0: since D is in general variable with concentration c, the concentration profile is not necessarily symmetric. Introducing X_M in the expression for D(c^*) above, however, introduces a bias that appears to make the value of D completely arbitrary, depending on which X_M we choose.

The values for X_M, however, are limited to a single one by physical meaningfulness. Since the denominator term dc/dx goes to zero for c → c_L (as the concentration profile flattens out), the integral in the numerator must also tend to zero in the same conditions, otherwise D(c_L) would tend to infinity, which is not physically meaningful (note that, strictly speaking, this does not guarantee that D does not tend to infinity, but it is only a necessary condition). The condition is then:

$0=\int _{c_{R}}^{c_{L}}(x-X_{M})\mathrm {d} c$

$X_{M}={\frac {1}{c_{L}-c_{R}}}\int _{c_{R}}^{c_{L}}x\mathrm {d} c$

In other words, X_M is the average position weighed on concentrations, and can be easily found from the concentration profile providing it is invertible to the form x(c).

Application of Boltzmann-Matano Method

Take for example two pairs of different alloys (A-B alloys) that are joined to form two classical diffusion couples. On the left, the alloy has a mole fraction of 0.25 of component B. And on the right, the alloy has a mole fraction of 0.75 of component B. We are given the following concentration-distance data.

We are assuming that one pair of alloys, exhibits a constant diffusivity, and upon diffusion yields a concentration field described by the usual linear Grube-Jedele diffusion couple solution.

$C(x,t)=0.75-0.25erfc\left({\frac {x}{4Dt}}\right)$

And then for the second pair of alloys, we assume it has a concentration-dependent diffusivity. For this example, we choose t=0.5.

$C\left(x,0.5\right)=0.25\left[tanh(x)+2\right]$

These equations are both ideal solutions for each of the pairs of alloys. Identical numerical procedures were applied to the two concentration fields, to implement the Boltzmann-Matano procedure for this example. Steps were taken for this particular example to determine the Boltzmann-Matano solution as well as the mass balance for the Matano’s interface. The following plot is the result for the calculation for the diffusivity for each set of alloys.

This example shows that when implementing the Boltzmann-Matano method, one can inversely solve for the diffusivity from given concentration-distance data. However, large errors occur at the end point data and could propose problems in any complex concentration profiles. Note as a reminder, this approach is for the classic diffusion couple. And when implemented in complex data, the errors seen in this ideal example may be significantly large and unacceptable in actual experimental data.

Sources

• M. E. Glicksman, Diffusion in Solids: Field Theory, Solid-State Principles, and Applications, Wiley, New York, 2000.
• Matano, Chujiro. "On the Relation between the Diffusion-Coefficients and Concentrations of Solid Metals (The Nickel-Copper System)". Japanese Journal of Physics. Jan. 16, 1933.