On the equation of state and atomic mean square displacement of crystals

We have presented a Green's function method for the calculation of the atomic mean square displacement (MSD) for an anharmonic Hamil toni an . This method effectively sums a whole class of anharmonic contributions to MSD in the perturbation expansion in the high temperature limit. Using this formalism we have calculated the MSD for a nearest neighbour fcc Lennard Jones solid. The results show an improvement over the lowest order perturbation theory results, the difference with Monte Carlo calculations at temperatures close to melting is reduced from 11% to 3%. We also calculated the MSD for the Alkali metals Nat K/ Cs where a sixth neighbour interaction potential derived from the pseudopotential theory was employed in the calculations. The MSD by this method increases by 2.5% to 3.5% over the respective perturbation theory results. The MSD was calculated for Aluminum where different pseudopotential functions and a phenomenological Morse potential were used. The results show that the pseudopotentials provide better agreement with experimental data than the Morse potential. An excellent agreement with experiment over the whole temperature range is achieved with the Harrison modified point-ion pseudopotential with Hubbard-Sham screening function. We have calculated the thermodynamic properties of solid Kr by minimizing the total energy consisting of static and vibrational components, employing different schemes: The quasiharmonic theory (QH), ).2 and).4 perturbation theory, all terms up to 0 ().4) of the improved self consistent phonon theory (ISC), the ring diagrams up to o ().4) (RING), the iteration scheme (ITER) derived from the Greens's function method and a scheme consisting of ITER plus the remaining contributions of 0 ().4) which are not included in ITER which we call E(FULL). We have calculated the lattice constant, the volume expansion, the isothermal and adiabatic bulk modulus, the specific heat at constant volume and at constant pressure, and the Gruneisen parameter from two different potential functions: Lennard-Jones and Aziz. The Aziz potential gives generally a better agreement with experimental data than the LJ potential for the QH, ).2, ).4 and E(FULL) schemes. When only a partial sum of the).4 diagrams is used in the calculations (e.g. RING and ISC) the LJ results are in better agreement with experiment. The iteration scheme brings a definitive improvement over the).2 PT for both potentials.

On the path integral formulation and the evaluation of quantum statistical averages

Four problems of physical interest have been solved in this thesis using the path integral formalism. Using the trigonometric expansion method of Burton and de Borde (1955), we found the kernel for two interacting one dimensional oscillators• The result is the same as one would obtain using a normal coordinate transformation, We next introduced the method of Papadopolous (1969), which is a systematic perturbation type method specifically geared to finding the partition function Z, or equivalently, the Helmholtz free energy F, of a system of interacting oscillators. We applied this method to the next three problems considered• First, by summing the perturbation expansion, we found F for a system of N interacting Einstein oscillators^ The result obtained is the same as the usual result obtained by Shukla and Muller (1972) • Next, we found F to 0(Xi)f where A is the usual Tan Hove ordering parameter* The results obtained are the same as those of Shukla and Oowley (1971), who have used a diagrammatic procedure, and did the necessary sums in Fourier space* We performed the work in temperature space• Finally, slightly modifying the method of Papadopolous, we found the finite temperature expressions for the Debyecaller factor in Bravais lattices, to 0(AZ) and u(/K/ j,where K is the scattering vector* The high temperature limit of the expressions obtained here, are in complete agreement with the classical results of Maradudin and Flinn (1963) .

Sedimentology and micromorphology of drumlin sediments at Chimney Bluffs State Park, New York, United States of America

The drumlin sediments at Chimney Bluffs, New York appear to represent a block-inmatrix style glacial melange. This melange comprises sand stringers, lenses and intraclasts juxtaposed in an apparently massive diamicton. Thin section examination of these glacigenic deposits has revealed microstructures indicative of autokinetic subglacial defonnation which are consistent with a deformable bed origin for the diamicton. These features include banding and. necking of matrix grains, oriented plasma fabrics and the formation of pressure shadows at the long axis ends of elongate clasts. Preservation of primary stratification within the sand intraclasts appears to suggest that these features were pre-existing up-ice deposits that were frozen, entrained, then deposited as part of a defonning till layer beneath an advancing ice sheet. Multi-directional micro-shearing within the sand blocks is thought to reflect the frozen nature of the sand units in such a high strain environment. It is also contended that dewatering of the sediment pile leading to the eventual immobilisation of the defonning till layer was responsible for opening sub-horizontal fissures within the diamicton. These features were subsequently infilled with mass flow poorly sorted sands and silts which were subjected to ductile defonnation during the waning stages of an actively deforming till layer. Microstructures indicative of the dewatering processes in the sand units include patches of fine-grained particles within a coarser-grained matrix and the presence of concentrated zones of translocated clays. However, these units were probably confined within an impermeable diamicton casing that prevented massive pore water influxes from the deforming till layer~ Hence, these microstructures probably reflect localised dewatering of the sand intraclasts. A layered subglacial shear zone model is proposed for the various features exhibited by the drumlin sediments. The complexity of these structures is explained in terms of ii superposing deformation styles in response to changing pore water pressures. Constructional glaciotectonics, as implied by the occurrence of sub-horizontal fissuring, is suggested as the mechanism for the stacking of the sand intraclast units within the diamicton. The usefulness of micromorphology in complimenting the traditional sedimentology of glacigenic deposits is emphasised by the current study. An otherwise massive diamicton was shown to contain microstructures indicative of the very high strain rates expected in a complexly deforming till layer. . It is quite obvious from this investigation that the classification of diamictons needs to be re-examined for evidence of microstructures that could lead to the re-interpretation of diamicton forming processes. RESUME Le pacquet de sediments drumlinaire de Chimney Bluffs, New York, represent un "bloc-en-matrice" genre de melange glaciale. Des structures microscopique comprennent l'evidence pour la defonnation intrinseque attribuee a l'origine lit non resistant du drumlin. PreselVation des structures primaires au coeur des blocs arenaces suggere que ceux sont des depots preexistant qui furent geles, entraines et par la suite sedimentes au milieu d'une couche de debris sous-glaciaires en voie de deformation. Des failles microscopiques a l'interieur des blocs arenaces appuient aussi l'idee d'un bloc cohesif (c'est-a-dire gele) au centre d'un till non resistant. Des implications significatives s'emergent de cette etude pour les conditions sous-glaciaire et les processus de la formation des drumlin.

The measured variation of the Debye-Waller factor of aluminum from 295K to 815K by using the energy dispersive x-ray diffraction technique

The Energy Dispersive X-ray Diffraction System at Brock University has been used to measure the intensities of the diffraction lines of aluminum powder sample as a function of temperature. At first, intensity measurements at high temperature were not reproducible. After some modifications have been made, we were able to measure the intensities of the diffraction lines to 815K, with good accuracy and reproducibility. Therefore the changes of the Debye-Waller factor from room temperature up to 815K for aluminum were determined with precision. Our results are in good agreement with those previously published.

Thermodynamic properties and Debye-Waller factor of fee materials

We have calculated the equation of state and the various thermodynamic properties of monatomic fcc crystals by minimizing the Helmholtz free energy derived in the high temperature limit for the quasiharmonic theory, QH, and the lowest-order (cubic and quartic), 'A2, anharmonic terms of the perturbation theory, PT. The total energy in each case is obtained by adding the static energy. The calculation of the thermal properties was carried out for a nearest-neighbour central-force model of the fcc lattice by means of the appropriate thermodynamic relations. We have calculated the lattice constant, the thermal expansion, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the isothermal and adiabatic bulk moduli, and the Griineisen parameter, for the rare-gas solids Kr and Xe, and gold. Morse potential and modified Morse potential were each used to represent the atomic interaction for the three fcc materials. For most of the calculated thermodynamic properties from the QH theory, the results for Kr and Xe with the modified Morse potential show an improvement over the results for the Morse potential when compared with the experimental data. However, the results of the 'A 2 equation of state with the modified Morse potential are in good agreement with experiment only in the case of the specific heat at constant volume and at constant pressure. For Au we have calculated the lattice contribution from the QH and 'A 2 PT and the electronic contribution to the thermal properties. The electronic contribution was taken into account by using the free electron model. The results of the thermodynamic properties calculated with the modified Morse potential were similar to those obtained with the Morse potential. U sing the minimized equation of state we also calculated the Mossbauer recoilless fraction for Kr and Xe and the Debye-Waller factor (DWF) for Pb, AI, eu, Ag, and Au. The Mossbauer recoilless fraction was obtained for the above two potentials and Lennard-Jones potential. The L-J potential gives the best agreement with experiment for Kr. No experimental data exists for Xe. At low temperature the calculated DWF results for Pb, AI, and eu show a good agreement with experimental values, but at high temperature the experimental DWF results increase very rapidly. For Ag the computed values were below the expected results at all temperatures. The DWF results of the modified Morse potential for Pb, AI, eu and Ag were slightly better than those of the Morse potential. In the case of Au the calculated values were in poor agreement with experimental results. We have calculated the quasiharmonic phonon dispersion curves for Kr, Xe, eu, Ag, and Au. The calculated and experimental results of the frequencies agree quite well for all the materials except for Au where the longitudinal modes show serious discrepancies with the experimental results. In addition, the two lowest-order anharmonic contributions to the phonon frequency were derived using the Green's function method. The A 2 phonon dispersion curves have been calculated only for eu, and the results were similar to those of the QH dispersion curves. Finally, an expression for the Griineisen parameter "( has been derived from the anharmonic frequencies, and calculated for these materials. The "( results are comparable with those obtained from the thermodynamic definition.

An anharmonic contribution to the Helmholtz free energy O(lambda 6)

The anharmonic contributions of order A6 to the Helmholtz free energy for a crystal in which every atom is on a site of inversion symmetry, have been evaluated The cor~esponding diagrams in the various orders of the perturbation theory have been presented The validity of the expressions given is for high temperatures. Numerical calculations for the diagrams which contribute to the free energy have been worked out for a nearest-n~ighbour central-force model of a facecentered cubic lattice in the high-temperature limit and in the leading term and the Ludwig approximations. The accuracy of the Ludwig approximation in evaluating the Brillouin-zone sums has been investigated. Expansion for all diagrams in the high-temperature limit has been carried out The contribution to the specific heat involves a linear as well as cubic term~ We have applied Lennard-Jones, Morse and Exponential 6 types of potentials. A comparison between the contribution to the free energy of order A6 to that of order A4 has been made.