Solid state NMR chemical shifts as an alternative to diffraction data in the determination of SIC polytypic structures

Silicon carbide, which has many polytypic modifications of a very simple and very symmetric structure, is an excellent model system for exploring, the relationship between chemical shift, long-range dipolar shielding, and crystal structure in network solids. A simple McConnell equation treatment of bond anisotropy effects in a poly type predicts chemical shifts for silicon and carbon sites which agree well with the experiment, provided that contributions from bonds up to 100 A are included in the calculation. The calculated chemical shifts depend on three factors: the layer stacking sequence, electrical centre of gravity, and the spacings between silicon and carbon layers. The assignment of peaks to lattice sites is proved possible for three polytypes (6H, 15R, and 3C). The fact that the calculated chemical shifts are very sensitive to layer spacings provides us a potential way to detennine and refine a crystal structure. In this work, the layer spacings of 6H SiC have been calculated and are within X-ray standard deviations. Under this premise, the layer spacings of 15R have been detennined. 29Si and 13C single crystal nmr studies of 6H SiC polytype indicate that all silicons and carbons are magnetically anisotropic. The relationship between a magnetic shielding tensor component and layer spacings has been derived. The comparisons between experimental and semi-empirical chemical shielding tensor components indicate that the paramagnetic shielding of silicon should be included in the single crystal chemical shift calculation.

An application of cell-cluster theory to a rare gas crystal /|n[by] K. Westera. - 260 St. Catharines [Ont.] : Dept. of Physics, Brock University,

The Lennard-Jones Devonshire 1 (LJD) single particle theory for liquids is extended and applied to the anharmonic solid in a high temperature limit. The exact free energy for the crystal is expressed as a convergent series of terms involving larger and larger sets of contiguous particles called cell-clusters. The motions of all the particles within cell-clusters are correlated to each other and lead to non-trivial integrals of orders 3, 6, 9, ... 3N. For the first time the six dimensional integral has been calculated to high accuracy using a Lennard-Jones (6-12) pair interaction between nearest neighbours only for the f.c.c. lattice. The thermodynamic properties predicted by this model agree well with experimental results for solid Xenon.

Molecular dynamics calculation of mean square displacement in alkali metals and rare gas solids and comparison with lattice dynamics

Molec ul ar dynamics calculations of the mean sq ua re displacement have been carried out for the alkali metals Na, K and Cs and for an fcc nearest neighbour Lennard-Jones model applicable to rare gas solids. The computations for the alkalis were done for several temperatures for temperature vol ume a swell as for the the ze r 0 pressure ze ro zero pressure volume corresponding to each temperature. In the fcc case, results were obtained for a wide range of both the temperature and density. Lattice dynamics calculations of the harmonic and the lowe s t order anharmonic (cubic and quartic) contributions to the mean square displacement were performed for the same potential models as in the molecular dynamics calculations. The Brillouin zone sums arising in the harmonic and the quartic terms were computed for very large numbers of points in q-space, and were extrapolated to obtain results ful converged with respect to the number of points in the Brillouin zone.An excellent agreement between the lattice dynamics results was observed molecular dynamics and in the case of all the alkali metals, e~ept for the zero pressure case of CSt where the difference is about 15 % near the melting temperature. It was concluded that for the alkalis, the lowest order perturbation theory works well even at temperat ures close to the melting temperat ure. For the fcc nearest neighbour model it was found that the number of particles (256) used for the molecular dynamics calculations, produces a result which is somewhere between 10 and 20 % smaller than the value converged with respect to the number of particles. However, the general temperature dependence of the mean square displacement is the same in molecular dynamics and lattice dynamics for all temperatures at the highest densities examined, while at higher volumes and high temperatures the results diverge. This indicates the importance of the higher order (eg. ~* ) perturbation theory contributions in these cases.

Monte Carlo study of the XY-model on quasi-periodic lattices

Monte Carlo Simulations were carried out using a nearest neighbour ferromagnetic XYmodel, on both 2-D and 3-D quasi-periodic lattices. In the case of 2-D, both the unfrustrated and frustrated XV-model were studied. For the unfrustrated 2-D XV-model, we have examined the magnetization, specific heat, linear susceptibility, helicity modulus and the derivative of the helicity modulus with respect to inverse temperature. The behaviour of all these quatities point to a Kosterlitz-Thouless transition occuring in temperature range Te == (1.0 -1.05) JlkB and with critical exponents that are consistent with previous results (obtained for crystalline lattices) . However, in the frustrated case, analysis of the spin glass susceptibility and EdwardsAnderson order parameter, in addition to the magnetization, specific heat and linear susceptibility, support a spin glass transition. In the case where the 'thin' rhombus is fully frustrated, a freezing transition occurs at Tf == 0.137 JlkB , which contradicts previous work suggesting the critical dimension of spin glasses to be de > 2 . In the 3-D systems, examination of the magnetization, specific heat and linear susceptibility reveal a conventional second order phase transition. Through a cumulant analysis and finite size scaling, a critical temperature of Te == (2.292 ± 0.003) JI kB and critical exponents of 0:' == 0.03 ± 0.03, f3 == 0.30 ± 0.01 and I == 1.31 ± 0.02 have been obtained.

Phonon dispersion curves and atomic mean square displacement for several fcc and bcc materials

The atomic mean square displacement (MSD) and the phonon dispersion curves (PDC's) of a number of face-centred cubic (fcc) and body-centred cubic (bcc) materials have been calclllated from the quasiharmonic (QH) theory, the lowest order (A2 ) perturbation theory (PT) and a recently proposed Green's function (GF) method by Shukla and Hiibschle. The latter method includes certain anharmonic effects to all orders of anharmonicity. In order to determine the effect of the range of the interatomic interaction upon the anharmonic contributions to the MSD we have carried out our calculations for a Lennard-Jones (L-J) solid in the nearest-neighbour (NN) and next-nearest neighbour (NNN) approximations. These results can be presented in dimensionless units but if the NN and NNN results are to be compared with each other they must be converted to that of a real solid. When this is done for Xe, the QH MSD for the NN and NNN approximations are found to differ from each other by about 2%. For the A2 and GF results this difference amounts to 8% and 7% respectively. For the NN case we have also compared our PT results, which have been calculated exactly, with PT results calculated using a frequency-shift approximation. We conclude that this frequency-shift approximation is a poor approximation. We have calculated the MSD of five alkali metals, five bcc transition metals and seven fcc transition metals. The model potentials we have used include the Morse, modified Morse, and Rydberg potentials. In general the results obtained from the Green's function method are in the best agreement with experiment. However, this improvement is mostly qualitative and the values of MSD calculated from the Green's function method are not in much better agreement with the experimental data than those calculated from the QH theory. We have calculated the phonon dispersion curves (PDC's) of Na and Cu, using the 4 parameter modified Morse potential. In the case of Na, our results for the PDC's are in poor agreement with experiment. In the case of eu, the agreement between the tlleory and experiment is much better and in addition the results for the PDC's calclliated from the GF method are in better agreement with experiment that those obtained from the QH theory.

Polytypism and Silicon carbide : a solid state nuclear magnetic resonance study

A survey of predominantly industrial silicon carbide has been carried out using Magic Angle Spinning nuclear magnetic resonance (MAS nmr); a solid state technique. Three silicon carbide polytypes were studied; 3C, 6H, and 15R. The 13C and 29 Si MAS nmr spectra of the bulk SiC sample was identified on the basis of silicon (carbon) site type in the d iff ere n t pol Y t Y pes • Out to 5.00 A fro mac en t r a lsi 1 i con (0 r carbon) atom four types of sites were characterized using symmetry based calculations. This method of polytype analysis was also considered, in the prelminary stages, for applications with other polytypic material; CdBr 2 , CdI 2 , and PbI 2 " In an attempt to understand the minor components of silicon carbide, such as its surface, some samples were hydrofluoric acid washed and heated to extreme temperatures. Basically, an HF removable species which absorbs at -110 ppm (Si0 2 ) in the 29 Si MAS nmr spectrum is found in silicon carbide after heating. Other unidentified peaks observed at short recycle delays in some 29 Si MAS nmr spectra are considered to be impurities that may be within the lattice. These components comprise less than 5% of the observable silicon. A Tl study was carried out for 29 Si nuclei in a 3C ii polytype sample, using the Driven Equilibrium Single-Pulse Observation of T1 (DESPOT) technique. It appears as though there are a number of nuclei that have the same chemical shift but different T1 relaxation times. The T1 values range from 30 seconds to 11 minutes. Caution has to be kept when interpreting these results because this is the first time that DESPOT has been used for solid samples and it is not likely in full working order. MAS nmr indicates that the 13C and 29 Si ~sotropic chemical shifts of silicon carbide appear to have a reciprocal type of relationship_ Single crystal nmr analysis of a 6H sample is accordance with this finding when only the resultant isotropic shift is considered. However, single crystal nmr also shows that the actual response of the silicon and carbon nuclear environment to the applied magnetic field at various angles is not at all reciprocal. Such results show that much more single crystal nmr work is required to determine the actual behavior of the local magnetic environment of the SiC nuclei.