## Abstract

As for the measurement of diffusion coefficient in transparent liquids by digital holographic interferometry, there are many kinds depending on the mathematical model and experimental setup. The method of using the distance of the peaks in concentration difference profile to determine diffusion coefficient by Mach-Zehnder interferometric experimental setup, is easy to handle. In order to improve the accuracy and convenience of measurement in this method, we combine procedures of hologram analysis used by Mialdun *et al* (2011) and those by He *et al* (2009). By using polynomial to fit the continuous phase difference (interference phase) of object beam, it can be more convenient and accurate to determine the distance between the two peaks. Besides, the shift of initial time has been discussed as a separated topic at the first time and two functions for minimization have been given for determination of the initial time in this paper. With the development, diffusion coefficient of KCl in water at 0.33mol/L and 25 °C has been measured. The diffusion coefficient is 1.844 × 10^{−9} m^{2}/s with the uncertainty of ± 0.012 × 10^{−9} m^{2}/s. Our measurement has more similar result with other methods than holographic interferometry. The relative uncertainty of diffusion coefficient in our experiment is less than 1% and total uncertainty of temperature is within ± 0.04 K.

© 2015 Optical Society of America

## 1. Introduction

The study of diffusivity is important in many fields such as chemical engineering, pollution control, biology and separation of isotopes. Knowledge of diffusivity is necessary for the design of chemical equipment and for the mass transfer studies [1,2]. Diffusion coefficients can be obtained by experimental methods or empirical correlations. However, experimental methods can give more accurate results than empirical correlations. Diffusion coefficients in transparent liquid mixtures can be measured by optical methods including conventional interferometry, holographic interferometry and electronic speckle pattern interferometry. Holographic interferometry is one of the most widely used techniques for diffusivity [3,4] because of requiring simpler elements and less effort than conventional interferometry [5] and being a whole-field technique [1,6,7]. Compared with holographic interferometry, the digital one is more popular because it doesn’t need wet processing, has no limits on taken frames and has the fast and simple nature.

The measurement of diffusion coefficients by digital holographic interferometry can be classified depending on mathematical models. There are four methods which are widely researched. The first one is using the distance of the peaks in concentration difference profile between two times to determine diffusion coefficient, which was put forward by Bochner and Pipman in double-exposure holographic interferometry [5]. The second one is using the location of the fringes in interference fringes, which was put forward by Gabelmann-Gray and Fenichel [8,9], and developed by Ruiz-Bevia *et al* [10]. The third one is using relationship between the deflection angle of object beam and the concentration distribution in diffusion cell, which was put forward by Chhaniwal, Anand, *et al* [11,12]. The last one is by fitting equations for the experimental concentration profiles at different times, which was put forward by Mialdun *et al* [13].

The second method used by Gabelmann-Gray and Fenichel has two disadvantages. One is that it requires the central bright fringe in the middle of the interference fringes [9], which is not always true because of a small wedge angle in the diffusion direction between the two windows in diffusion cell [5], or something else. The other is that it requires the time interval between two exposures to be short [10], in which situation it’s hard to detect the change of concentration accurately. Then Ruiz-Bevia *et al* had solved the problems by using integral formula and pairs of fringes of the same interference order [10]. But the method by Ruiz-Bevia *et al* is complex because it needs modifying the interface position between two solutions and calculating diffusion coefficients over and over. In the third method there’s no the reference beam [11,12], which makes the experimental setup simpler and possible to be an instrument. However, there is one backward. It’s hard to handle because it requires the deflection large enough to be detected accurately and small enough to keep the object beam collimated. The forth method can combine the measurement of diffusion coefficients and that of Soret coefficients. However, it’s more complex. The first method is relatively simpler to handle.

The method using the distance of the peaks in concentration difference (phase difference) profile, can be classified as several kinds depending on the experimental setup, such as multiple beam interferometric method [14], Michelson interferometric method [14], fringe projection method [14], lensless Fourier transform method [1,15], Mach-Zehnder interferometric method [16] and so on. As Chhaniwal *et al* [14] described, the former three methods all use the shift of the interference fringes between two times to determine the distance. Lensless Fourier transform method uses spherical wave as the reference beam, which makes it easier to get phase information of object beam by Fresnel-Kirchoff diffraction integral. In this method the phase difference between two times is obtained and a numerical reference wave-front tilted at a slight angle is used to convert the infinite fringes to finite ones, which has inflection points corresponding to the two peaks in concentration difference profile. Then the finite fringes are skeletalized for determination of the distance. If the one-dimensional phase difference profile can be obtained numerically, it can be more convenient and accurate for measurement of the distance. In the Mach-Zehnder interferometric method, a plane wave is used as the reference beam, which leads to that it’s convenient to pick up the phase information of object beam through filtering in frequency-domain, such as Mialdun *et al* [13] and He *et al* [16] did. In this method, it’s very easy for hologram analysis and to obtain the one-dimensional phase difference profile numerically.

Based on the above, the method of using the distance of two peaks in phase difference profile by digital holographic interferometry and Mach-Zehnder interferometric experimental setup is easy to handle, and the obtained one-dimensional phase difference profile makes it convenient and accurate to measure the distance. In this method there are two quantities to be measured for determination of diffusion coefficients, the diffusion time and distance between the peaks. In order to improve the accuracy of measurement by this method, it’s necessary to have appropriate procedures of hologram analysis and accurate determination of the diffusion time.

In the method there are two things to think about in the procedures of hologram analysis. One is the determination of the phase difference between two times, and the other is the unwrapping of phase difference. He *et al* [16] and Mialdun *et al* [13] have researched on these. They have different ideas on the procedures. The procedure for determination of phase difference used by He *et al*, is more convenient without a spectrum moving in Fourier domain to eliminate the carrier frequency *f*_{0}, compared with the one used by Mialdun *et al*. However, the unwrapping procedure used by Mialdun *et al* is simpler and has more accurate results than the one by He *et al*, which can be recognized through the unwrapping results in their papers. In a word, a combination of the procedures used by He *et al* and those used by Mialdun *et al* is made in this work.

In order to determine diffusion time accurately, the initial time for diffusion needs to be researched. The shift of the initial time in the method used by Bochner and Pipman was thought to be about 10 seconds [5]. In the method they chose freezing the solutions before putting them to contact, which cannot be used widely. In the experiment by Mialdun *et al* [17], the shift of the initial time was large to about 6 min in some condition and it was varying with concentration interval between solutions for diffusion. These mean that the shift of the initial time will be different depending on diffusion cell and the concentration interval between solutions. For example, a shift of the initial time with one minute will make a difference of 2.5% for diffusion coefficients, two minutes 5%, when the two exposure times are chosen at *t*_{1} = 30min, *t*_{2} = 50min. When the diffusion time is larger, the influence of the initial time is smaller. However, there is limitation for *t*_{2} (*t*_{2} > *t*_{1}) in order to satisfy the condition of one-dimensional infinite diffusion in the finite diffusion cell. At the same time, the time interval between two exposure times should be larger in order to detect the change of concentration accurately when diffusion time is larger, because the change of concentration between certain time interval will be smaller as time is increasing. Therefore, the first exposure time *t*_{1} can’t be chosen very large, which maybe makes obvious influence on measurement results with the shift of initial time. In order to eliminate this uncertain error, the initial time needs to be determined. In the work by Mialdun *et al* [13], the initial time is determined by fitting equation for many experimental concentration profiles at different times through varying diffusion coefficient and the initial time. However, this fitting way is not suitable for the method used in this paper. And the problem hasn’t been discussed in detail as a separate topic, even in the work by Mialdun *et al* [13,17]. In this paper, the reason for the shift of initial time is explained, and two functions for minimization are given to determine the initial time.

With the two kinds of development for the method, the diffusion coefficient of KCl in water at the concentration of 0.33mol/L and 25 °C is measured.

## 2. Theory

In a one-dimensional infinite medium the solution of the diffusion equation at time *t* for two liquids initially separated at *y* = 0, with concentrations *c*_{1} and *c*_{2} is

*D*is the diffusion constant with the difference of concentration

*c*

_{1}-

*c*

_{2}being small. Then the concentration difference between two times

*t*

_{1}and

*t*

_{2}can be expressed as:

With a plane object wave going through the diffusing solutions perpendicular to the diffusion direction, the information of concentration difference will be introduced into the phase of the object beam. Because refractive index of solutions is varying linearly with the concentration for small interval of concentration [5], the phase difference of object beam can be expressed as

*φ*is the phase of object beam,

*L*is the thickness of diffusing liquid, Δ

*n*is the difference of refractive index,

*λ*is the wavelength of light in vacuum,

*m*is the mean value of the derivative for the applied concentration range, and Δ

*φ*

_{0}is the phase difference without change of concentration. Because Δ

*φ*

_{0}is the same when

*y*is changed, the phase difference has the same profile with the concentration difference. The maximum and minimum of the concentration difference (phase difference) between two times

*t*

_{1}and

*t*

_{2}in the diffusion direction can be obtained by the formula

*w*is the distance between the maximum and minimum of the concentration difference (two peaks in the concentration difference profile). Therefore, we get the formula to calculate diffusion coefficient

With the help of digital holographic interferometry, we can get the phase difference profile of object beam to obtain the distance *w*. Finally, we can calculate diffusion coefficient with the three parameters *w*, *t*_{1} and *t*_{2}.

With one plane wave split into a reference beam and an object beam, we make the object beam go through the diffusion cell and interfere with the reference beam at some angle. Then we use a CCD camera with lens to record the hologram. Through analysis of the hologram, we can obtain the object beam.

The reference beam at the CCD plane can be expressed in complex amplitude, with *z*-direction perpendicular to the CCD plane

*R*

_{0}is the real amplitude,

*k*is the wave vector,

_{R}*ϕ*is the projection at the CCD plane from the phase distribution in the plane perpendicular to wave vector

_{R}*k*,

_{R}*θ*is the angle between wave vector and

_{Rx}*x*-direction, and

*θ*is the angle between wave vector and

_{Ry}*y*-direction. The object beam can be expressed as

The interference of the reference beam and the object beam will result in a spatially varying intensity distribution at the CCD plane. This hologram can be expressed as

*y*-direction parallel with the diffusion direction, making both

*θ*and

_{Ox}*θ*equal to 90 degrees. The hologram we get at time

_{Oy}*t*, can be converted into the Fourier domain

*F*means the Fourier transform,

*f*and

_{x}*f*are the frequency, and ⊗ means convolution. Because the reference beam and the object beam are plane waves,

_{y}*O*

_{0}and

*R*

_{0}don’t change with coordinate.

*G*

_{0}(

*f*,

_{x}*f*) is nearly located at origin of coordinate. If there is not concentration difference in the diffusion cell,

_{y}*ϕ*-

_{O}*ϕ*almost doesn’t change with coordinate and

_{R}*G*(

*f*-cos

_{x}*θ*/

_{Rx}*λ*,

*f*-cos

_{y}*θ*/

_{Ry}*λ*) presents a point at (cos

*θ*/

_{Rx}*λ*, cos

*θ*/

_{Ry}*λ*) in Fourier domain. Thinking about the concentration difference in the diffusion cell,

*ϕ*-

_{O}*ϕ*is the function of

_{R}*y*. So

*G*(

*f*-cos

_{x}*θ*/

_{Rx}*λ*,

*f*

_{y}-cos

*θ*/

_{Ry}*λ*) is actually a line in Fourier domain, and

*G**(

*f*+ cos

_{x}*θ*/

_{Rx}*λ*,

*f*+ cos

_{y}*θ*/

_{Ry}*λ*) is also a line. With appropriate angle

*θ*and

_{Rx}*θ*, the three parts

_{Ry}*G*

_{0},

*G*and

*G** are separated in Fourier domain. In a word, we can convert the hologram of object and reference beam into Fourier domain, obtain

*G*or

*G** by appropriate filter, and get

*O*

_{0}

*R*

_{0}exp[

*j*(2

*πy*cos

*θ*/

_{Ry}*λ*+ 2

*πx*cos

*θ*/

_{Rx}*λ*+

*ϕ*-

_{O}*ϕ*)] or

_{R}*O*

_{0}

*R*

_{0}exp[-

*j*(2

*πy*cos

*θ*/

_{Ry}*λ*+ 2

*πx*cos

*θ*/

_{Rx}*λ*+

*ϕ*-

_{O}*ϕ*)] through inverse Fourier transform.

_{R}Because both *θ _{Ox}* and

*θ*are equal to 90 degrees, which means the object beam is perpendicular to the CCD plane, the phase distribution of object beam

_{Oy}*φ*is equal to

*ϕ*. Then we can get the phase difference of object beam between two times

_{O}*t*

_{1}and

*t*

_{2}using the method by He

*et al*[16].

*θ*/

_{Rx}*λ*, cos

*θ*/

_{Ry}*λ*) or (-cos

*θ*/

_{Rx}*λ*, -cos

*θ*/

_{Ry}*λ*) to origin in Fourier domain.

The phase difference Δ*φ* calculated by above function is wrapped, which means that it belongs to the range (-*π*, *π*). We need convert it into the continuous phase difference, which makes it convenient and accurate for determination of distance between the peaks in phase difference numerically. The continuous phase difference obtained by He *et al* [16], has larger deviation with the theoretical phase difference profile, which can be recognized in the paper. However, the unwrapping procedure by Mialdun *et al* [13], is simple and has good results. Therefore, we unwrap the phase difference to construct the continuous natural phase as Mialdun *et al* did. If the difference is less than *π*, the phase remains unchanged. If the phase of the latter pixel is more than *π* bigger than the former in *y*-direction, the phases of the latter pixel and all pixels behind are minus 2*π*. If the phase of the latter pixel is more than *π* smaller than the former in *y*-direction, the phases of the latter pixel and all pixels behind are plus 2*π*.

## 3. The initial time

Figure 1 shows the sketch of the diffusion cell used in this work. Firstly, the lighter solution at some concentration is filled into the diffusion cell and all capillary tubes for draining air. Then the heavier solution is filled into the tubes through liquid inlet with exhaust port (exhaust port and liquid outlet) closed and liquid outlet open until the condition shown as Fig. 1(a). Subsequently, the heavier solution is admitted slowly from below lifting the lighter solution through liquid inlet with exhaust port open and liquid out closed, which is shown in Fig. 1(b) under the ideal condition. However, there will be mixing process with bulk-flow when the heavier solution is filled into the diffusion cell through tubes, because the section of the diffusion cell is much larger than that of tubes. This mixing process will make the interface between two solutions into a mixing area, which is shown in Fig. 1(c).

Although the time of mixing between two solutions is short, the speed of mass transfer by mixing with bulk-flow is much faster than that by diffusion. When the shift of initial time due to mixing with bulk-flow is transferred to that due to diffusion, it would be large. Therefore, the shift of initial time needs to be considered in the measurement. If it counts to be zero (time recording starts) when the interface between two solutions is established, the initial time *t*_{0} would be negative (*t* = *t*^{count}-*t*_{0} with *t*^{count} the recorded time and *t* the actual diffusion time). In the diffusion cell used by this work, there is some lighter solution in tubes when it’s started to lift the lighter solution, which is shown in Fig. 1(a). Therefore, it takes time until the interface is established from starting to lift the lighter solution. If it counts to be zero when the heavier solution is admitted to lift the lighter solution, the initial time *t*_{0} would be negative or positive depending on design of the diffusion cell and tubes.

Considering the influence of the initial time *t*_{0}, diffusion time *t* can be expressed as *t* = *t*^{count} – *t*_{0}, with *t*^{count} being the recorded time. Then the Eq. (4) can be expressed as

*t*is equal to

*t*

_{2}

^{count}-

*t*

_{1}

^{count}. As ∆

*t*is fixed, the distance

*w*is varying with time

*t*

_{2}

^{count}and there are two parameters

*D*and

*t*

_{0}. In the curve of function

*w*varying with time

*t*

_{2}

^{count}, the parameter

*D*determines the slope of the curve and the parameter

*t*

_{0}determines the shift of the curve, which can be recognized through Eq. (12) and Fig. 2. On the condition of ∆

*t*= 50 min and

*D*= 1.845 × 10

^{−9}m

^{2}/s, there are several curves with different initial time

*t*

_{0}, which is shown in Fig. 2(a). On the condition of ∆

*t*= 50 min and

*t*

_{0}= 0, there are several curves with different

*D*, which is shown in Fig. 2(b).

Based on the above, diffusion coefficient *D* and the initial time *t*_{0} are two separated parameters. Therefore, it can be used for determination of *t*_{0} to minimize the deviation between the experimental curve of *w* and the calculated one by varying *D* and *t*_{0}, when ∆*t* is fixed. The function for minimization is expressed as

*w*

^{cal}is calculated by Eq. (12) and

*w*

^{exp}is measured in the experimental phase difference profile. When the initial time is determined in this way, the diffusion coefficient is obtained at the same time.

Considering the influence of the initial time *t*_{0}, the Eq. (5) can be expressed as

*w*is varying with

*t*

_{2}

^{count}to keep the diffusion coefficient constant when ∆

*t*is fixed. Supposing that the initial time

*t*

_{0}is equal to 5 min and diffusion coefficient is equal to 1.845 × 10

^{−9}m

^{2}/s, a series of distance

*w*at different

*t*

_{2}

^{count}can be obtained by Eq. (12) with ∆

*t*being 50 min. Then the diffusion coefficients can be calculated with the series of distance

*w*. If shift of the initial time isn’t taken into consideration, the calculation results of diffusion coefficient by Eq. (5) will be different as

*t*

_{2}

^{count}is increasing, which is shown in Fig. 3. The same is for the situation

*t*

_{0}= −5 min, which is also displayed in Fig. 3. If the initial time is taken into consideration, the results of diffusion coefficient will be constant as Fig. 3 shows.

In Fig. 3, it can be found that calculation results of diffusion coefficient is varying with time without considering shift of initial time, when the initial time is not zero. Taking shift of initial time into consideration, the calculation results will be constant with the right initial time. Therefore, it can be used for determination of *t*_{0} to minimize difference among calculation results of diffusion coefficient by varying *t*_{0}. The function for minimization is expressed as

*D*

^{cal}is the calculation result by Eq. (15) and ${\overline{D}}^{\text{cal}}$ is the average value of different calculation results. Through minimization of the function expressed as (16), the average diffusion coefficient and initial time can be obtained at the same time.

## 4. Experiment

Figure 4 shows the schematic of the digital holographic interferometry setup employed for our measurement of diffusion coefficient. A spatially filtered and expanded He-Ne laser beam is split into a reference and an object beam. After passing through the diffusion cell perpendicular to the diffusion direction, the object beam interferes with the reference beam and makes a hologram, which is recorded by a CCD camera with lens. The CCD camera (FA-21-1M120), has 1024 × 1024 pixels with the pixel size of 7.4 × 7.4 μm^{2}. One pixel for the CCD camera corresponds to 43.79 um in the diffusion cell.

The experimental cell is shown like Fig. 5. The diffusion cell has three disks with a rectangular cell to hold the diffusing liquids in the middle disk. The rectangular cell measures 1.2 cm × 1.2 cm × 6 cm with glass plates on two sides to allow the laser light to pass through. The diffusion cell is put in water bath to be kept at the constant temperature.

There are some capillary tubes cycling in water bath for preheating liquid and these tubes are connected to the liquid inlet and have about half volume of the rectangular cell for diffusion. At first, the lighter solution at some concentration is filled into diffusion cell and all capillary tubes. Then the heavier solution with another concentration replaces the lighter liquid in the tubes through liquid inlet 3 and liquid outlet 4. After 20 minutes, the heavier solution is admitted slowly from below lifting the lighter liquid until the interface reaching the middle of the diffusion cell through liquid outlet 1 and liquid inlet 3. It counts to be zero when heavier solution starts to lift the lighter.

Water bath is provided and controlled by F33-EHJULABO. The fluctuation of temperature in the water bath is within 0.02 K. The temperature measured by the Pt100 (Fluke) has an uncertainty of ± 0.02 K. Considering both the fluctuation of temperature and the uncertainty of the temperature sensor Pt100, the total uncertainty for the temperature in the diffusion cell is within 0.04 K.

With this experimental setup, we measure the diffusion coefficient of KCl in water at the concentration of 0.33mol/L and 25 °C. The two kinds of concentration are 0.28 and 0.38 mol/L, as Ruiz-Bevia *et al* [10] chose.

## 5. Procedures of hologram analysis

The hologram is recorded as a picture every 5 minutes, as Fig. 6(a) shows. Each hologram is converted into Fourier domain, which is shown in Fig. 6(b). In the Fourier domain there are three parts, which include an ellipse and two lines. The ellipse corresponds to *G*_{0} and noise, whereas the two lines correspond to *G* and *G**. A rectangular filter is used to obtain *G* or *G**. The size and location of the rectangular filter window is chosen manually to cover one line, as Fig. 6(c) displays. After filtering, an inverse Fourier transform is used to get the phase information of the object beam, *O*_{0}*R*_{0}exp[*j*(2*πy*cos*θ _{Ry}*/

*λ*+ 2

*πx*cos

*θ*/

_{Rx}*λ*+

*ϕ*-

_{O}*ϕ*)] or

_{R}*O*

_{0}

*R*

_{0}exp[-

*j*(2

*πy*cos

*θ*/

_{Ry}*λ*+ 2

*πx*cos

*θ*/

_{Rx}*λ*+

*ϕ*-

_{O}*ϕ*)].

_{R}After getting the phase information of the object beam at times *t*_{1} and *t*_{2}, the phase difference (interference phase) can be obtained with Eq. (10), which is shown in Fig. 7. Through the procedure expressed above, we can unwrap the phase difference into continuous phase like Fig. 8.

After getting the continuous phase difference, we can use polynomial to fit the phase difference in *y*-direction (diffusion direction). The result for one profile in *y*-direction is shown in Fig. 9, in which we can find that the fitting curve and the original phase difference coincide very well. The fitting value in each pixel is obtained for determination of the peaks in the phase difference profile.

After getting the location of the peaks at each *x* coordinate, we can get the distribution of the peaks of concentration difference (phase difference), which is shown in Fig. 10(a). Through this figure we can make the conclusion that the assumption of one-dimensional diffusion is satisfied in the experiment. Besides, we can get the distance *w* of the peaks in concentration difference profile at each *x* coordinate like Fig. 10(b). With the average in *x*-direction, we can get the average *w* for calculating the diffusion coefficient.

Because the size and location of the rectangular filter window is chosen manually depending on the figure in Fourier domain, the measurement result of the distance *w* is different between two times of analysis. In order to overcome the error, the analysis of each pair of holograms is treated twice to get average value of distance *w*. The average value is *w*, whereas the stand deviation is *σ* and the relative deviation is *δw* = *σ*/*w*. When the relative deviation is large, it means that noise has obvious influence on the hologram analysis. This phenomenon is always happened when the time interval between two exposure times is relatively short. When the relative deviation is small, it means that the hologram analysis has good accuracy.

## 6. Experimental results

We measure the diffusion coefficient of KCl in water at 0.33mol/L and 25 °C. The hologram is recorded as a picture every 5 minutes. Through any two pictures, we can get a distance *w* between the two peaks of phase difference. Each pair of holograms has been treated twice to get the average. The average value is *w*, whereas the stand deviation is *σ* and the relative deviation is *δw* = *σ*/*w*. The measurement results are shown in Table 1. The left is for the condition ∆*t* = 40 min, whereas the right is for the condition ∆*t* = 50 min. The relative deviation in Table 1 is less than 0.5%, which means that the results have good accuracy.

With a series of distance *w* at different *t*_{2}^{count}, the diffusion coefficient can be determined by minimizing function *F*_{1} when ∆*t* is fixed. This way to determine diffusion coefficient is called the first way. And the way by minimizing function *F*_{2} is called the second way. With some pairs of *w* and *t*_{2}^{count} at fixed ∆*t*, it can lead to a pair of diffusion coefficient and initial time with the first way. The same is for the second way. In order to certify the stability of results, numbers of pairs of *w* and *t*_{2}^{count} are changed from 20 to 9, as Table 2 and Table 3 show. The measurement results of diffusion coefficient by both first and second way are shown in Table 2 and Table 3. The distribution of diffusion coefficients calculated by two ways is shown in Fig. 11.

With the average of all results for diffusion coefficients, we obtain the final diffusion coefficient of KCl in water at the concentration of 0.33mol/L and 25 °C, which is 1.844 × 10^{−9} m^{2}/s with the uncertainty of ± 0.012 × 10^{−9} m^{2}/s. The relative uncertainty of diffusion coefficient is less than 1%. The contrast of our measurement result with literature data is shown in Table 4. The deviation of our result with Harned and Gosting is positive and less than 0.3%, whereas the deviation with Szydlowska and Ruiz-Bevia, who used holographic interferometry, is negative and about 1%.

## 7. Conclusion

As for the measurement of diffusion coefficient in transparent liquids by digital holographic interferometry, there are many kinds depending on mathematical model and experimental setup. The method of using the distance of the peaks in concentration difference profile to determine diffusion coefficient by Mach-Zehnder interferometric experimental setup, is easy to handle. In order to improve the accuracy and convenience of hologram analysis in this method, a combination of procedures for hologram analysis used by Mialdun *et al* and those used by He *et al*, has been made. After getting continuous phase difference of object beam between two times, the experimental phase difference profile is fitted by polynomial. It turns out that the fitting curve has good coincidence with the experimental. This makes it convenient and accurate to determine the distance of the peaks in concentration difference profile.

In order to determine diffusion time accurately, the shift of the initial time needs to be considered. When the heavier solution is injected into the diffusion cell to lift the lighter solution, there will be mixing process with bulk-flow because the section of the diffusion cell is much larger than the tubes for filling. Although the time of mixing between two solutions is short, shift of the initial time will be large when the influence of mixing with bulk-flow on the initial time is transferred to that of diffusion. In the method used by this work, two functions for minimization are given to determine the initial time.

With the development, diffusion coefficient of KCl in water at 0.33mol/L and 25 °C has been measured. The measurement result of diffusion coefficient is 1.844 × 10^{−9} m^{2}/s with the uncertainty of ± 0.012 × 10^{−9} m^{2}/s. Our measurement has more similar result with other methods than holographic interferometry. The relative uncertainty of diffusion coefficient in our experiment is less than 1%. The total uncertainty of temperature in our experiment is within ± 0.04 K.

## Acknowledgments

This work was supported by the Major State Basic Research Development Program of China (973 Program No. 2015CB251502).

## References and links

**1. **A. Anand, V. K. Chhaniwal, and C. S. Narayanamurthy, “Diffusivity studies of transparent liquid solutions by use of digital holographic interferometry,” Appl. Opt. **45**(5), 904–909 (2006). [CrossRef] [PubMed]

**2. **P. P. Chikode, S. J. Pawar, V. J. Fulari, and M. B. Dongare, “Study of diffusion process in sucrose solution by using double exposure holographic interferometry,” J. Opt. **36**(4), 157–167 (2007).

**3. **R. Riquelme, I. Lira, C. Perez-Lopez, J. A. Rayas, and R. Rodriguez-Vera, “Interferometric measurement of a diffusion coefficient: comparison of two methods and uncertainty analysis,” J. Phys. D Appl. Phys. **40**(9), 2769–2776 (2007). [CrossRef]

**4. **J. Colombani and J. Bert, “Holographic interferometry for the study of liquids,” J. Mol. Liq. **134**(1-3), 8–14 (2007). [CrossRef]

**5. **N. Bochner and J. Pipman, “A simple method of determining diffusion constants by holographic interferometry,” J. Phys. D Appl. Phys. **9**(13), 1825–1830 (1976). [CrossRef]

**6. **Y. Zhang, J. Zhao, J. Di, H. Jiang, Q. Wang, J. Wang, Y. Guo, and D. Yin, “Real-time monitoring of the solution concentration variation during the crystallization process of protein-lysozyme by using digital holographic interferometry,” Opt. Express **20**(16), 18415–18421 (2012). [CrossRef] [PubMed]

**7. **J. M. Desse, P. Picart, and P. Tankam, “Digital three-color holographic interferometry for flow analysis,” Opt. Express **16**(8), 5471–5480 (2008). [CrossRef] [PubMed]

**8. **L. Gabelmann-Gray and H. Fenichel, “Holographic interferometric study of liquid diffusion,” Appl. Opt. **18**(3), 343–345 (1979). [CrossRef] [PubMed]

**9. **H. Fenichel, H. Frankena, and F. Groen, “Experiments on diffusion in liquids using holographic interferometry,” Am. J. Phys. **52**(8), 735–738 (1984). [CrossRef]

**10. **F. Ruiz-Bevia, A. Celdran-Mallol, C. Santos-Garcia, and J. Fernandez-Sempere, “Holographic interferometric study of free diffusion: a new mathematical treatment,” Appl. Opt. **24**(10), 1481–1484 (1985). [CrossRef] [PubMed]

**11. **V. K. Chhaniwal, A. Anand, and B. S. Chakrabarty, “Diffusion studies in transparent liquid mediums utilizing polarization imaging,” Opt. Lasers Eng. **46**(12), 888–892 (2008). [CrossRef]

**12. **V. Chhaniwal, C. S. Narayanamurthy, and A. Anand, “Imaging of mass transfer process using artificial fringe deflection,” Opt. Eng. **53**(7), 074106 (2014). [CrossRef]

**13. **A. Mialdun and V. Shevtsova, “Measurement of the Soret and diffusion coefficients for benchmark binary mixtures by means of digital interferometry,” J. Chem. Phys. **134**(4), 044524 (2011). [CrossRef] [PubMed]

**14. **V. K. Chhaniwal, A. Anand, S. Girhe, D. Patil, N. Subrahmanyam, and C. S. Narayanamurthy, “New optical techniques for diffusion studies in transparent liquid solutions,” J. Opt. A, Pure Appl. Opt. **5**(5), S329–S337 (2003). [CrossRef]

**15. **G. Sheoran, A. Anand, and C. Shakher, “Lensless Fourier transform digital holographic interferometer for diffusivity measurement of miscible transparent liquids,” Rev. Sci. Instrum. **80**(5), 053106 (2009). [CrossRef] [PubMed]

**16. **M. He, Y. Guo, Q. Zhong, and Y. Zhang, “A new method of processing Mach-Zehnder interference fringe data in determination of diffusion coefficient,” Int. J. Thermophys. **30**(6), 1823–1837 (2009). [CrossRef]

**17. **A. Mialdun, V. Sechenyh, J. C. Legros, J. M. Ortiz de Zárate, and V. Shevtsova, “Investigation of Fickian diffusion in the ternary mixture of 1,2,3,4-tetrahydronaphthalene, isobutylbenzene, and dodecane,” J. Chem. Phys. **139**(10), 104903 (2013). [CrossRef] [PubMed]

**18. **H. S. Harned and R. L. Nuttall, “The differential diffusion coefficient of Potassium Chloride in aqueous solutions,” J. Am. Chem. Soc. **71**(4), 1460–1463 (1949). [CrossRef]

**19. **J. Gosting, “A study of the diffusion of Potassium Chloride in water at 25° with the Gouy interference method,” J. Am. Chem. Soc. **72**(10), 4418–4422 (1950). [CrossRef]

**20. **J. Szydlowska and B. Janowska, “Holographic measurement of diffusion coefficients,” J. Phys. D Appl. Phys. **15**(8), 1385–1393 (1982). [CrossRef]

**21. **J. A. Rard and D. G. Miller, “Mutual diffusion coefficients of NaCl_{2}-H_{2}O and CaCl- H_{2}O at 25 °C from Rayleigh interferometry,” J. Chem. Eng. Data **25**, 211–215 (1980). [CrossRef]

**22. **V. M. M. Lobo, A. C. F. Ribeiro, and L. M. P. Verissimo, “Diffusion coefficients in aqueous solutions of Potassium Chloride at high and low concentrations,” J. Mol. Liq. **78**(1-2), 139–149 (1998). [CrossRef]