The knowledge of diffusion parameters provides important understanding of many physical and mechanical properties of materials. In most of the applications silicides are grown by a diffusion controlled process mainly in thin film condition. Because of this reason, most of the studies till date are available in thin film condition. Although more than one phase is present in all these systems, mainly disilicides were found at the interface. In this thesis bulk interdiffusion studies are conducted by coupling pure refractory metals (group IVB, VB and VIB elements) with single crystal Si. Several phase layers grow between binary refractory metal and Si systems. The layer thicknesses of the phases are measured from the microstructures. Composition profiles were measured in electron probe microanalyzer. Different diffusion parameters are estimated such as parabolic growth constants, integrated diffusion coefficients, activation energy for diffusion and ratio of tracer diffusivities of the components are estimated. Growth mechanisms of the phases are discussed with the help of diffusion parameters. The atomic mechanism of the diffusion is discussed considering crystal structure of the phases along with possible defects present. Solid diffusion couple experiments are conducted to analyse the growth mechanism of the phases and the diffusion mechanism of the components in the Ti-Si system. The calculation of the parabolic growth constant and of the integrated diffusion coefficients substantiates that the analysis is intrinsically prone to erroneous conclusions if it is based just on the parabolic growth constants determined for a multiphase interdiffusion zone. The location of the marker plane is detected based on the uniform grain morphology in the TiSi2 phase, which indicates that this phase grows mainly because of Si diffusion. The growth mechanism of the phases and morphological evolution in the interdiffusion zone are explained with the help of imaginary diffusion couples. The activation enthalpies for the integrated diffusion coefficient of TiSi2 and the Si tracer diffusion are calculated as 190±9 and 170±12 kJ/mol, respectively. The crystal structure, details on the nearest neighbours of the elements and the relative mobilities of the components indicate that the vacancies are mainly present on the Si sublattice. Diffusion controlled growth of the phases in the Hf-Si and Zr-Si are studied by bulk diffusion couple technique. Only two phases grow in the interdiffusion zone, although several phases are present in both the systems. The location of the Kirkendall marker plane detected based on the grain morphology indicates that the disilicides grow by the diffusion of Si. Diffusion of the metal species in these phases is negligible. This indicates that vacancies are present mainly on the Si sublattice. The activation energies for integrated diffusion coefficients in the HfSi2 and ZrSi2 are estimated as 394 ± 37 and 346 ± 34 kJ/mol, respectively. The same is calculated for the HfSi phase as 485±42 kJ/mol. The activation energies for Si tracer diffusion in the HfSi2 and ZrSi2 phases are estimated as 430 ± 36 and 348 ± 34 kJ/mol, respectively. We conducted interdiffusion studies to understand the atomic mechanism of the diffusing species and the growth mechanism of the phases. Integrated diffusion coefficients and the ratio of tracer diffusion coefficients were estimated for these analyses. The activation energies for the integrated diffusion coefficients were calculated as 550 ± 70 and 410 ± 39 kJ/mol in the TaSi2 and the Ta5Si3 phases, respectively. In the TaSi2 phase, Ta has a slightly lower but comparable diffusion rate with respect to Si, although no TaTa bonds are present in the crystal. In the Ta5Si3 phase, Si has higher diffusion rate, which is rather unusual, if we consider the atoms in the nearest-neighbor positions for both the elements. The ratio of Si to Ta tracer diffusion coefficients is found to be lower in the Si-rich phase, TaSi2, compared to the Si-lean phase, Ta5Si3, which is also unusual. This indicates the type of structural defects present. An analysis on the growth mechanism of the phases indicates that duplex morphology and the Kirkendall marker plane should only be present in the TaSi2 phase. This is not present in the Ta5Si3 phase because of the very high growth rate of the TaSi2 phase, which consumes most of the Ta5Si3 phase layer. The problems in the calculation method used previously by others in this system are also explained. Experiments are conducted in the W-Si system to understand the diffusion mechanism of the species. The activation energies for integrated diffusion are found to be 152±7 and 301±40 kJ/mol in the WSi2 and W5Si3 phases, respectively. In both the phases, Si has a much higher diffusion rate compared to W. The result found in the WSi2 phase is not surprising, if we consider the nearest neighbors in the crystal. However, it is rather unusual to find that Si has higher diffusion rate in the W5Si3 phase, indicating the presence of high concentration of Si antisites in this phase. In the group IVB, VB and VIB M-Si systems are considered to show an interesting pattern in diffusion of components with the change in atomic number in a particular group. MSi2 and M5Si3 are considered for this discussion. Except in the Ta-Si system, activation energy for integrated diffusion of MSi2 is always lower than M5Si3. Interestingly, in both the phases, the relative mobilities measured by the ratio of tracer D* diffusion coefficients, S i decreases with the increase in atomic number in both the DM* groups. Both the phases have similar crystal structures in a particular group in which these parameters are calculated. In both the phases Si has higher diffusion rate compared to M. Absence of any M-M bonds in MSi2 and increase in the diffusivities of M with the increase in atomic number substantiates the increasing concentration of M anti-sites and higher interactions of M with vacancies. Only one or two Si-Si bonds are present in M5Si3, however, the higher diffusion rate of Si indicates the presence of vacancies mainly D* on its sublattice. On the other hand, increase in S i with increasing atomic number in DM* Both the groups substantiates increasing interactions of M and vacancies.