New Analytical Method of Interdiffusion Problems

The nonlinear diffusion equation for a binary system interdiffusion was analytically solved in the previous work. The theoretical relation of Kirkendall effect was also derived in the previous work. These new results have not yet been concretely applied to actual diffusion problems. In the present work, it is revealed that the previous results reproduce the experimental concentration profile by taking account of the movement of diffusion region space. It is thus actually confirmed that any problems of binary system interdiffusion can be solved by the new analytical method if even diffusivities of self-diffusion and impurity diffusion in the materials concerned are given. The method for solving interdiffusion problems of many elements system, which is extremely important for the development of new useful materials, is also reasonably discussed. Further, it is revealed that the concept of intrinsic diffusion is unsuitable for the diffusion theory. The fundamental theory of diffusion discussed here will be useful for analyzing actual diffusion problems in future.


Introduction
The diffusion problem is one of the most fundamental and important research subjects in the material science field.The diffusion research has been thus widely and actively performed in accordance with the industry requirements for the development of new useful materials [1] [2] [3] [4] [5].However, the progress of fundamental theory of diffusion had been hardly seen until recently.A lot of research papers accepting the Darken equation [6], which is mathematically wrong in the derivation process [7], have been still reported [8] [9] [10] [11] [12].
It was thus required to establish a new method for solving interdiffusion problems.In that situation, the analytical method of the nonlinear diffusion equation, which is applicable to analyzing interdiffusion problems, was reported [13] [14] [15].Nevertheless, it does not seem that the method is really applied to analyzing results of diffusion experiments.In the present study, therefore, the analytical theory of interdiffusion problems is concretely expanded and is applied to experimental results.
First of all, we mention here fundamental concept of the diffusion theory in mathematics.In general, the usual experiments of interdiffusion between elements I and II are performed within such a temperature region that the normalized concentrations I C and II C satisfy the relation of The diffusion equations for I C and II C in the interdiffusion field are and where I D and II D are diffusivities for the elements I and II.Substituting Equation (1) into Equations ((2) and (3)) yields The above equations mean using the so-called interdiffusion coefficient of Here, note that Equation ( 5) is valid only in the differential equation for I C or II C [15] [16].The relation between diffusion fluxes I II I There is thus no doubt that the Kirkendall effect (K effect) relevant to the essence of interdiffusion mechanism is caused by I II 0 J J + ≠ [15]  [16] [17] [18].Nevertheless, the mathematical basics mentioned here had not been noticed in the long history of diffusion.
In order to understand the K effect, therefore, the intrinsic diffusion coeffi- [12], the equation itself is not only mathematically wrong in the derivation process but also entirely meaningless because of using the intrinsic diffusion coefficients nonexistent in the recent diffusion theory [7] [15].In the present work, the physical meaning of K effect is reasonably discussed compared with experimental results, regardless of the Darken equation.As discussed at the conclusion section, the Gauss divergence theory indicates that the K effect may occur in the diffusion problems as universal phenomena.
Applying the analytical solutions of Equation ( 4) to experimental results of interdiffusion problems [17] [19] [20] [21], we confirmed that the concentration profiles are reasonably reproduced in the present work.Hereafter, as mentioned above, the Darken equation as well as the concept of intrinsic diffusion will disappear from the theory of interdiffusion problems.On the other hand, the new method discussed here will become one of the most useful methods for analyzing interdiffusion problems.

Mathematical Theory of Analytical Method
The diffusion Equations ((2) and ( 3)) in the time and space ( ) where the suffixes I and II are removed.When the diffusivity depends on the concentration, it had been believed until recently that the mathematical solutions of the nonlinear diffusion Equation ( 7) are impossible.However, the mathematical method for solving Equation ( 7) was established in the previous work [13] [14].This means that the diffusion problems including a many elements system interdiffusion are essentially solved [18].The matters necessary for the present study are briefly summarized in the following.
Boltzmann transformed Equation ( 7) into the ordinary differential equation of using the parabolic law x t ζ = [22].In the previous work [13], Equation (8) was rewritten as a formula of diffusion flux in the parabolic space yielding ( ) ( ) ( ) where ( ) ( ) In mathematics, the dependence of diffusivity on the concentration means

1) Impurity diffusion
In case of the impurity diffusion, since the diffusivity of Equation ( 9) corresponds approximately to a constant value where ( ) The solution of Equation ( 11) is obtained as for the initial condition of

2) Binary system interdiffusion
By solving simultaneously Equations (( 9) and ( 10)), mathematical solutions of Equation ((8) or ( 9)) are possible (See Ref. [13]).As shown in Figure 1 ,  t x   against the origin of ( ) , t x .The notation • means an inert marker.The partial figures (a), (b)   and (c) represent the initial state at a room temperature, a state at a high temperature during diffusion treatment and a state at a room temperature after diffusion treatment, respectively. where We confirmed that Equations (( 13) and ( 14)) agree well with results of the empirical Boltzmann Matano method [13] [22] [23].
In order to specify the general solutions, we must determine the initial and/or boundary values of diffusivities in the interdiffusion problem.The general method for determining them had not been experimentally and theoretically known until recently.However, they can be reasonably obtained by using the self-diffusion coefficient and impurity diffusion coefficient, since the analytical solutions are obtained.Here, note that Equations (( 13) and ( 14)) yield the dependence of diffusivity on the concentration given by ( ) ( ) As can be easily seen, problems of binary system interdiffusion in the material composed of an arbitrary rate between elements I and II are solved by investigating the interdiffusion problem between a pure material I and a pure material II.In that case, the above initial and/or boundary values for Equations (( 13) and ( 14)) are physically accepted as of the material j itself and impurity diffusion coefficient imp j D of the material j in the other material.The solutions of interdiffusion problem between pure materials I and II are thus obtained as ( ) and ( ) where If we can experimentally obtain even diffusivities of self-diffusion and impurity diffusion in a material concerned, initial and/or boundary values necessary for solving an interdiffusion problem are thus obtained from using Equation (16).
For the interdiffusion problem of an arbitrary diffusion couple between the material A composed of and the material B composed of (16) shows that the initial and/or boundary values of diffusivities are The initial and/or boundary values of diffusivities corresponding to each element in a material composed of an arbitrary concentration rate for elements I and II were thus obtained as Equations ((17a) and (17b)).Therefore, the diffusivity profile of Equation ( 13) is obtained by using these diffusivity values of Equations ((17a) and (17b)).At the same time, the concentration profile of Equation ( 14) is also obtained by using these diffusivity values.As a matter of course, the diffusivity values of Equations ((17a) and (17b)) are also applied to int j D , j α and j β in Equations (( 13) and ( 14)).
In addition, the K effect shows that the diffusion region space, which is composed of vacancies and/or interstices among micro particles in a material, moves through the migration of their micro particles [17].It is also physically considered that the diffusion region space interacts with the free space near the surface of specimen considered to be the sink and source [18].This means that the coordinate transformation of the diffusion equation is necessary for analyzing interdiffusion problems [15] [16].Using the solutions obtained here, the relation of K effect eff x ∆ was obtained as ( ) ( )

3) Determination of initial and/or boundary values of diffusivities in a ternary system interdiffusion
It is revealed that the so-called interdiffusion coefficient in a many elements system is meaningful only in the differential equation of diffusion [15] [18].In the following, we first investigate the interdiffusion problem between the material A composed of elements I and II and the pure material B composed of an element III.Their initial concentrations are expressed as in the diffusion region.
In that case, the initial and/or boundary values of diffusivities in the ternary system, The present solutions of a ternary system interdiffusion problem are thus easily obtained as follows [18].
)   Further, using the following relations of ( ) { } ( ) ( ) the dependences of diffusivity on the concentration are obtained as ( are reasonably valid in the diffusion region of the interdiffusion problem mentioned above.
The initial and/or boundary values of diffusivities in the ternary system interdiffusion, which are applicable to the material A composed of the concentration rate In the above theory, the K effect is obtained as ( ) where As can be seen from the above analytical method, the successive calculations yield the solutions of a ( ) elements system interdiffusion problem by solving an interdiffusion problem between a material composed of N elements and a pure material.
H. Cho et al.

Analysis of Experimental Results
Kirkendall found that Zn atoms diffuse faster than Cu atoms in the Zn-Cu alloy [17] [19] [20].The molybdenum wire set at a point on the initial interface between the pure copper and the brass (Cu-70%, Zn-30%) was used as an inert marker then.The diffusion treatment was performed at the temperature 880˚C during 119 days in the experiment.In the present work, we apply the above theory to the historical experimental results first known as the K effect.The initial values of concentration in Figure 1 are then.In that case, however, it is difficult to apply the analytical method resulting in Equations ((17a) and (17b)), since the pure Zn melts at a temperature of 420˚C.We cannot thus understand the selfdiffusivity Zn B,self D in the above theory.On the other hand, the diffusivity of a metal near the temperature of melting point becomes about ) is thus valid in the present case [15] [18].For the diffusion length Dt µ , 2 µ = is used as a matter of convenience for Equation (18).Here, the impurity diffusivity of Zn in the pure Cu given by using the defined Zn int D and Zn β .
In Figure 2, the black curve denotes the concentration profile of Zn given by Equation (29) of the fixed coordinate system at the temperature 880˚C and diffusion time  Considering the inert characteristic of marker, the movement of an inert marker shows that the diffusion region space moves to the inverse orientation against the movement of diffusion particles.As discussed in the previous works, the coordinate transformation of the diffusion equation is thus indispensable for understanding the K effect [15] [16] [18].In the present case, therefore, the origin of fixed coordinate system ( ) , t x is set at the point of mass center on the in- itial interface and the origin of moving coordinate system ( )  is set at a point of space on the initial interface under the condition of t t =  and where v is a velocity of the origin of ( ) , t x   against the origin of ( ) = ∫  into Equation (29), the concentra- tion profile for 0 x > at a room temperature is expressed as ( )

Zn
Zn sft Zn int 3 1 erf 20 2 temperature fall from 880˚C to a room temperature after diffusion treatment, the supersaturated vacancies for 0 x > flow to the negative direction of x axis because of the unsaturated state of vacancies in the region of 0 x < , and at the same time they also flow to the specimen surface considered to be the sink and source of vacancies then.The inert marker does not move to the direction of x axis even if vacancies flow to the specimen surface perpendicular to the x axis.
Therefore, the inert marker on the

Discussion and Conclusions
Using the analytical solutions expressed by Equations ((29) and ( 30)), the interdiffusion problem between the pure copper and brass alloy was reasonably solved in the present study, regardless of the Darken equation.It was concretely confirmed that the Darken equation, which has been widely used for analyzing problems of binary system interdiffusion, is not only actually unnecessary but also theoretically unsuitable for analyzing the interdiffusion problems.By applying the diffusion theory obtained previously to the actual interdiffusion problem, the necessity of coordinate transformation for the diffusion equation was concretely confirmed in the present work.At the same time, it was also confirmed that the concept of intrinsic diffusion is nonexistent from the beginning.
The quantity of excess vacancies corresponding to ( ) x x − ∆ returns to the negative direction of x axis after the diffusion treatment.In other words, the quantity of excess vacancies corresponding to eff x ∆ flows toward the specimen surface.The K effect thus corresponds to the quantity of vacancies absorbed by the specimen surface considered to be the sink and source of vacancies.
In the region of pure copper, it is considered that diffusion behavior is approximately acceptable as an impurity diffusion mechanism if we neglect the effect of vacancies for interdiffusion problems.In fact, the result based on the concept was reasonably obtained.This gives evidence that the analytical method for solving interdiffusion problems of many elements system is valid, since we can determine initial and/or boundary values of diffusivities as discussed in the Section 2.
For the developments of new useful materials, solving interdiffusion problems between materials is fundamentally one of the most important research subjects.
The new analytical method discussed here will be widely applicable to analyzing the interdiffusion problems of many elements system.Behavior of vacancies plays extremely important role in the interdiffusion problems.In order to understand further detailed diffusion behavior, solving the diffusion equation of vacancies in the diffusion region will be necessary [15] [18].
Here, the conclusions obtained from the present work are summarized as follows.
1) It was confirmed that the concept of intrinsic diffusion is not only unnecessary but also wrong from a viewpoint of mathematical physics.The Darken equation should not be used for interdiffusion problems.
2) The Gauss divergence theorem indicates that the diffusion flux ( ) 3) The equation of is also valid only in the differential equation of diffusion.On the other hand, when we solve Equation (4) in accordance with the given initial and/or boundary values different from each other, the diffusivi-.When we use solutions of Equations (4) for the diffusion flux J used for analyzing interdiffusion problems[6] [8][9] [10][11]

Figure 1 .
Figure 1.Schematic figure of interdiffusion phenomena.The coordinate systems of ( ) , t x and ( ) , t x  are set at a point of the mass center and at a point of space on the initial interface between materials A and B, respectively.The relation between their coordinate systems resulting from the movement of diffusion region space is conceived as t t =  and

D
17a) and (17b)) obtained by analyzing the problems of binary system interdiffusion are rewritten as I of III in the material A. In the same manner, those values diffusion coefficients of elements I, II in the material B and the self-diffusion coefficient of the material B.

) 4 )
Ternary system interdiffusionFor a ternary system interdiffusion problem between materials composed of an arbitrary concentration rate, for example, a material A composed of I that case, we must determine the initial and/or boundary values of diffusivities for each element in the materials A and B. In Equation(20), we first determine the point of A

D
manner, we can determine the initial and/or boundary values of diffusivities which are applicable to the ternary system in- terdiffusion in the material B composed of the concentration rate I initial and/or boundary values, the analytical solutions of the ternary system interdiffusion are reasonably obtained as ( ) kind[8].The  value of Zn B,self D is determined as follows in the present work.The relation of K effect generalized from Equation (18) is used for determining the value of Zn B,self D compared with the empirical relation of eff

D
value in the present study[21].Substituting the initial values, as a physically reasonable value in the present study.Here, the solution of basic diffusion Equation(7) for the present problem is expressed as

5 4 .
284 10 s t = × .In the same figure, the red curve denotes the concentration profile of Zn at a room temperature after the diffusion treatment and

Figure 2 .
Figure 2. Concentration profile of Zn in the diffusion region of Zn-Cu alloy.The notation ○ indicates an experimental value at a given point in the diffusion region.The black curve shows Equation (29) at the temperature 880˚C and the time

,
t x .Equation (27) indicates that the velocity v of the origin of ( ) since the excess vacancies disappear during a temperature fall after the diffusion treatment.In other words, the 0 x =  axis and the inert marker on the sft after the diffusion treatment[15] [16].The red curve approximately agrees with the experimental results for 0 dotted red curve is over shifted in the region of 0 x < compared the experimental results.Resulting from the diffusion behavior of Zn atoms, the material B side of diffusion couple is in a supersaturated state of vacancies, while the material A side is in an unsaturated state of vacancies.The red curve of Equation (30) should be thus accepted only in the region for 0 x > .In the present case, the material A is the pure Cu and it seems that Zn atoms on the interface nearA x x = in the diffusion region of A B x x x ≤ ≤diffuse into the pure Cu as impurities, if we neglect the vacancy behavior in the diffusion field.This indicates that we can thus adopt 0 x < as a first approximation.The dotted blue curve denotes the concentration of Zn for 0 x < at a room temperature.Further, taking account of the effect of vacancies for 0 x < , the blue curve shifts to the positive direction of the x axis under the condition of the shift resulting from the unsaturated vacancies for 0 x < .The shifted blue curve agrees approximately with the experimental results for 0 x < at the room temperature after the diffusion treatment.As a result, this gives an evidence for the validity of present method for determining an initial and/or boundary value of diffusivity discussed in the section 2.Using Equations ((29), (30)), the concentration profile of experimental results was physically reproduced as shown by the red curve for 0 x > and the blue curve for 0 x < in Figure2.However, the correlation between eff x ∆ and sft x is considerably complicated as discussed in the following.As mentioned above, the diffusion region space moves to the positive direction of x axis in the present case.The inert maker on the initial interface moves to sft the temperature fall because of the vacancy flow along the x axis.H. Cho et al.
]. Here, the flux ( ) J t means a movement of space within the diffusion region and the K effect gives evidence for its validity.The flux eq J independent of the time and space plays an impor- tant role in the self-diffusion mechanism[15].It was also confirmed that the as Equation (4)-(10) at the page 2124 in Ref.[8] should be revised as