Social Icons

PRACTICAL 4 : Determination of Diffusion Coefficient






TITLE:

Determination of Diffusion Coefficient

DATE OF EXPERIMENT:

3th  May 2013

OBJECTIVES:

1.      To study the diffusion of molecules in an agar medium.
2.      To determine the diffusion coefficient of crystal violet and bromothymol blue at different temperature.

INTRODUCTION:
           
            Diffusion is a process of the movement of solute molecules in the response to applied force, occurring in a system which is not in equilibrium. The examples of applied forces are temperature, concentration, pressure and electric current. According to Fick’s law, the diffusion flux is proportional to the minus gradient of concentrations. The solute molecules move from region of higher concentration to region of lower concentration. Fick’s law of diffusion can be used to calculate the diffusion coefficient. Fick’s law states that the flux of material (amount dm in time dt) across a given plane (area A) is proportional to the concentration gradient dc/dx.

(1)     dm = -DA(dc/dx)dt


where D is the diffusion coefficient or diffusivity for the solute, m2s-1.

            If a solution containing neutral particles with the concentration M0, is placed within cylindrical tube next to a water column, diffusion can be stated as

                                        (2)     M = M0 exp (–x2/4Dt)



where M is the concentration at distance x from the intersection between water and solution that is measure at time t.

By changing equation (2) to its logarithmic form, we get

ln M = ln M0 –x2/4Dt

OR
                                       (3) 2.303 x 4D (log10 M0 – log10 M) t = x2



Thus a plot of x2 against t can produce a straight line that passes through the origin with the slope 2.303 x 4D (log10 Ma – log10 M). From here D can be calculated.

            If the particles in the solution are assumed to be spherical, their size and molecular weight can be calculated by the Stokes-Einstein equation.

(4) D = kT/6pha

where
·         k is the Boltzmann constant 1.38 x 1023 Jk-1,
·         T is the temperature in Kelvin,
·         p is the viscosity of the solvent in Nm-2s, and
·         a is the radius of particle in M.

The volume of a spherical particle is 4/3 pa3, thus its weight M is equivalent to 4/3 pa3Nr (r = density)

It is known that molecular weight M=mN (N is Avogadro’s number 6.023 x 1023 mol-1).

  (5) M = 4/3 pa3Nr


            If charged particles are involved in the diffusion, equation (3) needs to be modified to include potential gradient effect that exists between the solution and solvent. However, this can be overcome by adding a little sodium chloride into the solvent to prevent the formation of this potential gradient.

            Agar gels contain a partially strong network of molecules that are penetrated by water. The water molecules form a continuous phase around the gel. Thus, the solute molecules can diffuse freely in the water if chemical interactions and adsorption effects do not exist entirely. Therefore, the gel forms a suitable support system to be used in diffusion studies for molecules in a medium of water.




APPARATUS:

10ml, 50ml, 100ml and 250ml measuring cylinders, 1000ml beaker, 14 test tubes with stoppers, weighing boat, spatula, glass rod, electronic balance, hot plate and stirrer, test tube rack, dropper, filter funnel,marker and water bath.

MATERIALS:

Jelly Powder, Ringlet’s solution, 1 : 200 crystal violet indicator, 1 : 400 crystal violet indicator, 1 : 600 crystal violet indicator, 1 : 500 000 crystal violet indicator, 1 : 200 bromothymol blue solution, 1 : 400 bromothymol blue solution, 1 : 600 bromothymol blue solution, 1 : 500 000 bromothymol blue solution, and distilled water.






PROCEDURES:


1.      Agar in Ringer’s solution was prepared as four tablets were crushed by using pestle and mortar.

2.      Then, the crushed tablets were transfered into the 1000ml beaker and 500ml of distilled water was added into the beaker to dissolve the tablets.

3.      7g of agar powders were weighed on electronic balance and 425ml Ringer’s solution was measured in a suitable measuring cylinders.

4.       Agar powders and Ringer’s solution were mixed in the beaker. The mixed solution was heated until the solution completely dissolved and boiled.

5.       The solution was continuously stirred during the heating process to prevent the jelly powder from embed in the mixture.

6.       During the heating process, 14 test tubes were labelled according to the concentrations and types of the solutes, which are bromotyhmol blue and crystal violet. The heating process was stopped when the mixture become clear dissolved solution.

7.       20ml of agar solution was prepared in 14 labelled test tubes.

8.      5ml of 1:500 000 crystal violet solution was added into a test tube containing hot agar solution labelled with 1:500 000 crystal violet at 28°C whereas 5ml of 1:500 000 bromothymol blue solution was added into test tube containing the hot gel solution labelled with 1:500 00 bromothymol blue at 37°C. These two test tubes will be used as the standard to measure the colour distance resulting from the solute diffusion.

9.      The agar solution in the rest of test tubes were allowed to cool at room temperature until they become solid agar.

10.  5ml of each crystal violet solution was added into the gels that were prepared in 6 test tubes according to their concentration.

11.  This step was repeated by using bromothymol blue as an indicator. All test tubes were closed to prevent evaporation and stored at temperature 28°C and 37°C.

12.  The distance of crystal violet and bromothymol blue solution diffused in the agar were measured and recorded for seven days.


RESULTS:

a)  Crystal Violet

System
Time (seconds)
x, cm
x2, cm2
Slope of graph
D, cm2s-1
Temp., o C
Average Diffusion Coefficient, D cm2s-1
1:200








0
0
0
3.057 X 10-5
 9.7662 x 10-7
28







 2.8652 x 10-7






86400
1.3
1.69
172800
2.1
4.41
259200
2.8
7.84
345600
3.3
10.89
432000
3.7
13.69
518400
4.1
16.81
604800
4.3
18.49



1:400








0
0
0
1.911 X 10-5
6.6986 x
10-7
28







86400
1.2
1.44
172800
1.8
3.24
259200
2.4
5.76
345600
2.7
7.29
432000
3.0
9.00
518400
3.3
10.89
604800
3.4
11.56



1:600





0
0
0
1.488X 10-5
5.5281 x
10-7
28







86400
0.75
0.56
172800
1.6
2.56
259200
2
4.0
345600
2.3
5.29
432000
2.5
6.25
518400
2.8
7.84
604800
3.0
9.0


1:200








0
0
0
 3.652 X 10-5
 1.16669 x10-6
37







9.2097 x 10-7






86400
1.6
2.56
172800
2.4
5.76
259200
3.0
9.00
345600
3.6
12.96
432000
4.0
16.00
518400
4.4
19.36
604800
4.7
22.09



1:400








0
0
0
2.646 X 10-5
 9.27497 x 10-7
37







86400
1.4
1.96
172800
2.0
4.00
259200
2.6
6.76
345600
3.0
9.00
432000
3.4
11.56
518400
3.8
14.44
604800
4.0
16.00



1:600







0
0
0
1.800 X 10-5
6.6871 x 10-7
37






86400
1.0
1.00
172800
1.7
2.89
259200
2.2
4.84
345600
2.5
6.25
432000
2.8
7.84
518400
3.1
9.61
604800
3.3
10.89












 








b)  Bromothymol Blue


System
Time (seconds)
x, cm
x2, cm2
Slope of graph
D, cm2S-1
Temp, o C
Average Diffusion Coefficient, D cm2S-1
A (1:200)








0
0
0
2.917 X 10-5








9.31886 x 10-7








28
























 1.25302 x
10-6







86400
1.5
2.25
172800
2.3
5.29
259200
2.7
7.29
345600
3.2
10.24
432000
3.6
12.96
518400
3.9
15.21
604800
4.2
17.64



B (1:400)








0
0
0
2.388 X 10-5








8.3703 x 10-7
28







86400
1.3
1.69
172800
1.9
3.61
259200
2.5
6.25
345600
2.9
8.41
432000
3.2
10.24
518400
3.6
12.96
604800
3.8
14.44



C (1:600)








0
0
0
 5.357 X 10-6








1.99017 x 10-6
28







86400
0.4
0.16
172800
0.8
0.64
259200
1.1
1.21
345600
1.3
1.69
432000
1.4
1.96
518400
1.6
2.56
604800
1.8
3.24



D (1:200)








0
0
0
3.652 X 10-5








1.16669 x 10-6 
37





















1.3313 x 10-6











86400
1.7
2.89
172800
2.5
6.25
259200
3.1
9.61
345600
3.6
12.96
432000
4
16
518400
4.4
19.36
604800
4.7
22.09



E (1:400)








0
0
0
2.388 X 10-5








8.3703 x
10-7







37







86400
1.3
1.69
172800
2.0
4.00
259200
2.5
6.25
345600
2.9
8.41
432000
3.2
10.24
518400
3.5
12.25
604800
3.8
14.44



F (1:600)








0
0
0
8.003 X 10-6








1.9901 x 10-6
37







86400
0.7
0.49
172800
1.1
1.21
259200
1.4
1.96
345600
1.6
2.56
432000
1.9
3.61
518400
2.1
4.41
604800
2.2
4.84

















CALCULATIONS:

From equation:
2.303 x 4D (log 10 Mo – log 10 M) t = X²
Hence the slope of the graph = 2.303 x 4D (log 10 Mo - log 10 M)

1.      Crystal violet system with dilution 1:200 (28 0C)
Slope= 3.057 x 10-5cm2/sec
            M = 1:500000                       M0 = 1:200
    = 1 / 500000                          = 1 / 200
    = 2 x 10-6                                         = 5 x 10-3
           
2.303x4D [log 10 (5x10-3)-log 10 (2x10-6)] = 3.057×10-5 cm2/sec
D = 9.7662×10-7 cm2/sec

2.      Crystal violet system with dilution 1:400 (28ºC)
Slope = 1.911×10-5 cm2/sec
            M = 1:500000                      Mo = 1:400
     = 1 / 500000                         = 1 / 400
     = 2 x 10-6                              = 2.5 x 10-3

2.303x4D [log 10 (2.5x10-3)-log 10 (2x10-6)] = 1.911×10-5 cm2/sec
D = 6.6986×10-7 cm2/sec

3.      Crystal violet system with dilution 1:600 (28 0C)
     Slope=1.488×10-5 cm2/sec

            M = 1:500000                      Mo = 1:600
     = 1 / 500000                         = 1 / 600
     = 2 x 10-6                                        = 1.67 x 10-3
    

2.303x4D [log 10 (1.67x10-3)-log 10 (2x10-6 )] = 1.488×10-5 cm2/sec
D = 5.5281×10-7 cm2/sec

Average of Diffusion Coefficient, m²/sec for Crystal violet system at 28ºC
= (9.7662×10-7 cm2/sec+6.6986×10-7 cm2/sec+5.5281×10-7 cm2/sec) / 3
= 2.8652×10-7 cm2/sec

4.      Crystal violet system with dilution 1:200 (37ºC)
           Slope =3.652×10-5 cm2/sec

            M = 1:500000                       M0 = 1:200
    = 1 / 500000                          = 1 / 200
    = 2 x 10-6                                         = 5 x 10-3
           
            2.303x4D [log 10 (5x10-3)-log 10 (2x10-6 )]  = 3.652×10-5 cm2/sec
 D=1.16669×10-6 cm2/sec

5.      Crystal violet system with dilution 1:400 (37ºC)
Slope =2.646×10-5 cm2/sec

            M = 1:500000                      Mo = 1:400
     = 1 / 500000                         = 1 / 400
     = 2 x 10-6                              = 2.5 x 10-3
           
2.303x4D [log 10 (2.5x10-3)-log 10 (2x10-6)] = 2.646×10-5 cm2/sec
D = 9.267497×10-7 cm2/sec

6. Crystal violet system with dilution 1:600 (37ºC)
Slope = 1.8×10-5 cm2/sec

            M = 1:500000                      Mo = 1:600
     = 1 / 500000                         = 1 / 600
     = 2 x 10-6                                        = 1.67 x 10-3 

2.303x4D [log 10 (1.67x10-3)-log 10 (2x10-6 )] = 1.8×10-5 cm2/sec

D = 6.871×10-7 cm2/sec

     Average of Diffusion Coefficient, m²/sec for Crystal violet system at 37ºC
      = (1.16669×10-6 cm2/sec + 9.267497×10-7 cm2/sec +6.871×10-7 cm2/sec) / 3
      = 9.2097×10-7 cm2/sec

# Value of D37°C for crystal violet using the equation
D28°C/D37°C = T28°C/T37°C
2.8652×10-7/D37°C = (28°C + 273.15K)/(37°C + 273.15K)
                    D 37°C = 2.9508×10-7 cm2/sec

7.      Bromotymol blue system with dilution 1:200 (28ºC)

Slope  = 2.917×10-5 cm2/sec

M = 1:500000                      Mo = 1:200
= 1 / 500000                          = 1 / 200
= 2 x 10-6                                          = 5 x 10-3
           

          2.303x4D [log 10 (5x10-3)-log 10 (2x10-6)] = 2.917×10-5 cm2/sec
D = 9.1886×10-7 cm2/sec

8.      Bromotymol blue system with dilution 1:400 (28ºC)
Slope = 2.388×10-5 cm2/sec

            M = 1:500000                      Mo = 1:400
    = 1 / 500000                          = 1 / 400
    = 2 x 10-6                                         = 2.5 x 10-3
2.303x4D [log 10 (2.5x10-3)-log 10 (2x10-6 )] = 2.388×10-5 cm2/sec
             D= 8.3703×10-7 cm2/sec
     9. Bromotymol blue system with dilution 1:600 (28ºC)
Slope =5.357×10-6 cm2/sec
             M = 1:500000                      Mo = 1:600
     = 1 / 500000                         = 1 / 600
     = 2 x 10-6                                        = 1.67 x 10-3

          2.303x4D [log 10 (1.67x10-3)-log 10 (2x10-6 )] = 5.357×10-6 cm2/sec
D = 1.9907×10-7 cm2/sec

Average of Diffusion Coefficient, m²/sec for Bromothymol blue system at 28ºC
      = (9.1886×10-7 cm2/sec + 8.3703×10-7  cm2/sec +11.9907×10-7 cm2/sec) / 3
      = 1.253×10-6 cm2/sec

10.  Bromotymol blue system with dilution 1:200 (37ºC)
Slope = 3.625×10-5 cm2/sec


            M = 1:500000                       Mo = 1:200
    = 1 / 500000                          = 1 / 200
    = 2 x 10-6                                         = 5 x 10-3

2.303x4D [log 10 (5x10-3)-log 10 (2x10-6 )] = 3.625×10-5 cm2/sec
D = 1.1669×10-6 cm2/sec

11. Bromotymol blue system with dilution 1:400 (37ºC)
      Slope  =2.388×10-5 cm2/sec
            M = 1:500000                      Mo = 1:400
     = 1 / 500000                         = 1 / 400
     = 2 x 10-6                              = 2.5 x 10-3
          
          2.303x4D [log 10 (2.5x10-3)-log 10 (2x10-6)] = 2.3888×10-5 cm2/sec
D = 8.3703×10-7 cm2/sec

12. Bromotymol blue system with dilution 1:600 (37ºC)
Slope  =8.003×10-6 cm2/sec

            M = 1:500000                      Ma = 1:600
     = 1 / 500000                         = 1 / 600
     = 2 x 10-6                                        = 1.67 x 10-3

2.303x4D [log 10 (1.67x10-3)-log 10 (2x10-6 )] = 8.003×10-6 cm2/sec
            D = 1.9901×10-6 cm2/sec
  
 Average of Diffusion Coefficient, m²/sec for Bromotymol blue system at 37ºC
      = (1.1669×10-6 cm2/sec +8.3703×10-7cm2/sec +1.9901×10-6 cm2/sec) / 3
      = 1.3313×10-6 cm2/sec


# Value of D37°C for bromothymol blue using the equation
D28°C/D37°C = T28°C/T37°C
1.253×10-6 /D37°C = (28°C + 273.15K)/(37°C + 273.15K)
                    D 37°C = 1.2904×10-6 cm2/sec



DISCUSSION:

                Diffusion is the movement of solute particles in a solid from a high concentration area to a low concentration area, resulting in the uniform distribution of the substance. Diffusion is a process which is not due to the force action but due to the random movements of atoms. In each diffusion reaction the flux of matter is equal to the conductivity multiplied by a driving force. Conductivity is the mobility of the diffusing species or known as diffusivity. The presence of concentration gradient acts as the driving force for the solute particles to move.

            Fick’s First Law states that the flux, dm, across a membrane of unit area, A, is proportional to the concentration gradient, dc/dx, and is expressed by

dm = -DA(dc/dx)dt

D is the diffusion coefficient or diffusivity for the solute. The value of diffusivity describes the rate of diffusion and therefore its unit is m2s-1. Diffusion coefficient depends on several factors such as temperature, pressure, composition, physical state, structure of the phase, and oxygen fugacity. The higher the value of diffusion coefficient, the easier for that solute to penetrate through the continuous phase.

The negative sign indicates that the direction of the diffusive flux is down the concentration gradient. Diffusive flux goes from high concentration region to low concentration region with the gradient going the opposite direction, from low to high concentration. Fick’s First Law only applies to steady state flux where the concentration gradient is uniform.

            In this practical, a solution containing neutral molecules with the concentration M0 is placed within a cylindrical tube. The diffusion can be expressed as

M = Mo exp (x²/4Dt) 

The equation is derived into logarithmic form and produce new equation,

ln M = ln Mo – (x²/4Dt)

                                          or    2.303 x 4D (log 10 Mo- log 10 M) t = x²

Therefore, a graph of x2 against time produces a straight line that passes through origin with the slope of 2.303 x 4D (log 10 Mo- log 10 M). D can be calculated from the gradient obtained.

            The results show that the diffusion coefficient, D, for both crystal violet and bromothymol blue is increasing from 1:600<1:400<1:200 which mean that 1:200 concentration is the easiest to diffuse as 1:200 concentration is the most concentrated solution. Therefore, the solute with high concentration is easier to diffuse through the solid medium compared to the low concentration. This is because the concentration gradient is high for concentrated solute and forces the solute particles to diffuse further through the medium.

            If we compare the diffusion coefficient in the term of temperature between 28°C and 37°C, the D for 37°C has the larger value than that of 28°C.  The same phenomena occurred for both crystal violet and bromothymol blue solution. Temperature is one of the factors that affect the rate of diffusion. As the temperature increases, the amount of energy available for diffusion is increased. There would be increase in molecules' mobility (kinetic energy). The molecules move faster and there will be more spontaneous spreading of the material which means that diffusion occurs quicker. Thus, the diffusivity increases as the temperature increases.

The value of D37°C can be determined if we already know the value for D28°C by using the equation below

D28°C/D37°C = T28°C/T37°C

where T is the temperature in Kelvin.
However, there are the differences between the value of D37°C that we get from the experiment and from the theoretical value. For crystal violet, 9.2099×10-7 cm2s-1 is the experimental value of D37°C whereas 7.5502×10-7 cm2s-1 is the theoretical value that obtained from the equation above. Meanwhile, for bromothymol blue the D37°C from theoretical and experimental is 6.7559×10-7 cm2s-1 and 7.6703×10-7 cm2s-1 respectively. Generally, the experimental value for D37°C is higher than theoretical value. This is maybe due to the some errors that occurred during the experiment is carried out. The temperature may be set higher than 37°C that cause the rate of diffusion to increase. Besides, the errors also due to the room temperature at 28°C that is not constant throughout the experiment and the viscosity of agar in the test tube that is not uniform. Misreading or parallax error might occur.    

According to the equation:
M   = 4/3πa3 N Þ
a3   = 3M / 4πN Þ
a     = 3√ 3M / 4πN Þ    --------------------------- (1)
Where M= molecular weight, a = radius of particle, N=Avogadro's number (6.02x1023),
Þ =density.
Another equation:                    
D = kT / 6 πηa --------------------------- (2)
Substituting (1) into the equation (2)
D = kT / 6 πηa
    = kT / 6 πη 3√ 3M/4πNÞ
    = kT 3√4πNÞ / 6 πη 3√ 3M ------------------------- (3)
Where D=diffusion coefficient, η= viscosity, T=temperature, k=Boltzmann constant                (1.38 x1023 Jk-1)

According to equation (3), D is inversely proportional to M. So, D increases when M decreases. Crystal violet with molecular formula C25N3H30Cl has molecular weight of 407.979 g mol-1 while bromothymol blue solution with molecular formula C27H28Br2O5S has molecular weight of 624.38 g mol−1. Molecular weight is how much mass each particle has or how heavy it is. The heavier the particle, the slower it diffuses into solidified agar solution, assuming energy of the system remains constant. Crystal violet solution diffuses faster than bromothymol blue since the value of diffusion coefficient for crystal violet is higher than bromothymol blue at both temperatures. Crystal violet solution diffuses easily because its molecular weight is much smaller compared to that of bromothymol blue and thus easy for them to penetrate through gel medium. 

CONCLUSION:

From the experiment, the diffusion coefficient of crystal violet at 28°C is 7.3311×10-7 cm2s-1 and at 37°C is 9.2099×10-7 cm2s-1. Meanwhile, for bromothymol blue solution the diffusion coefficient at 28°C is 6.5599×10-7 cm2s-1 and at 37°C is 7.6703×10-7 cm2s-1. The diffusion coefficient of crystal violet is generally higher than bromothymol blue because of the difference inmolecular weight. In term of temperature, diffusivity at 37°C is much higher than diffusivity at 28°C. This is due to the rapid movement of the solute particles as the temperature increases.

REFERENCES:

3.      physics.nyu.edu/grierlab/methods/node11.htm

1 comment:

  1. hello..i want to ask you something..where did you get the value of log m ?

    ReplyDelete