Structural, stratigraphic and geochemical studies of the Horwood Peninsula - Gander Bay Area, Northeast Newfoundland /

The Horwood Peninsula - Gander Bay area is located at NE Newfoundland in the Botwood Zone (Williams et a1., 1974) or in the Dunnage Zone (Williams, 1979) of the Central Mobile Belt of the Newfoundland Appalachians. The area is underlain by Middle Ordovician to possible Lower Silurian rocks of the Davidsville and Indian Islands Groups, respectively. Three conformable formations named informally : the Mafic Volcanic Formation, the Greywacke and Siltstone Formation and the Black Slate Formation, have been recognized in the Davidsville Group. The Greywacke and the Black Slate Formations pass locally into a Melange Formation. From consideration of regional structure and abundant locally-derived mafic volcanic olisto- 1iths in the melange, it is considered to have originated by gravity sliding rather than thrusting. Four formations have been recognized in the Indian Islands Group. They mainly contain silty slate and phyllite, grey cherty siltstone, green to red micaceous siltstone and limestone horizons. Repetition of lithological units by F1 folding are well-demonstrated in one of formations in this Group. The major structure in this Group on the Horwood Peninsula is interpreted to be a synclinal complex. The lithology of this Group is different from the Botwood Group to the west and is probably Late Ordovician and/or Early Silurian in age. The effects of soft-sediment deformation can be seen from the lower part of the Davidsville Group to the middle part of the Indian Islands Group indicating continuous and/or episodic slumping and sliding activities throughout the whole area. However, no siginificant depOSitional and tectonic break that could be assigned to the Taconian Orogeny has been recognized in this study. Three periods of tectonic deformation were produced by the Acadian Orogeny. Double boudinage in thin dikes indicates a southeast-northwest sub-horizontal compression and main northeast-southwest sub-horizontal extension during the D1 deformation. A penetrative, axial planar slaty cleavage (Sl) and tight to isocJ.ina1 F1 folds are products of this deformation. The D2 and D3 deformations formed S2 and S3 fabrics associated with crenulations and kink bands which are well-shown in the slates and phyllites of the Indian Islands Group. The D2 and D3 deformations are the products of vertical and northeast-southwest horizontal shortening respectively. The inferred fault between the Ordovician slates (Davidsville Group) and the siltstones (Indian Islands Group) suggested by Williams (1963, 1964b, 1972, 1978) is absent. Formations can be followed without displacement across this inferred fault. Chemically, the pillow lavas, mafic agglomerates, tuff beds and diabase dikes are subdivided into three rock suites : (a) basaltic komatiite (Beaver Cove Assemblage), (b) tholeiitic basalt (diabase dikes), (c) alkaline basalt (Shoal Bay Assemblage). The high Ti02 , MgO, Ni contents and bimodal characteristic of the basaltic komatiite in the area are comparable to the Svartenhuk Peninsula at Baffin Bay and are interpreted to be the result of an abortive volcano-tectonic rift-zone in a rear-arc basin. Modal and chemical analyses of greywackes and siltstones show the trend of maturity of these rocks increasing from poorly sorted Ordovician greywackes to fairly well-sorted Silurian siltstones. Rock fragments in greywackes indicate source areas consisting of plagiogranite, low grade metamorphic rocks and ultramafic rocks. Rare sedimentary structures in both Groups indicate a southeasterly provenance. Trace element analyses of greywackes also reveal a possible island-arc affinity.

A treatment of atomic mean square displacement in higher order perturbation theory

A general derivation of the anharmonic coefficients for a periodic lattice invoking the special case of the central force interaction is presented. All of the contributions to mean square displacement (MSD) to order 14 perturbation theory are enumerated. A direct correspondance is found between the high temperature limit MSD and high temperature limit free energy contributions up to and including 0(14). This correspondance follows from the detailed derivation of some of the contributions to MSD. Numerical results are obtained for all the MSD contributions to 0(14) using the Lennard-Jones potential for the lattice constants and temperatures for which the Monte Carlo results were calculated by Heiser, Shukla and Cowley. The Peierls approximation is also employed in order to simplify the numerical evaluation of the MSD contributions. The numerical results indicate the convergence of the perturbation expansion up to 75% of the melting temperature of the solid (TM) for the exact calculation; however, a better agreement with the Monte Carlo results is not obtained when the total of all 14 contributions is added to the 12 perturbation theory results. Using Peierls approximation the expansion converges up to 45% of TM• The MSD contributions arising in the Green's function method of Shukla and Hubschle are derived and enumerated up to and including 0(18). The total MSD from these selected contributions is in excellent agreement with their results at all temperatures. Theoretical values of the recoilless fraction for krypton are calculated from the MSD contributions for both the Lennard-Jones and Aziz potentials. The agreement with experimental values is quite good.

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.

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.

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.

The effects of plants(Typha Latifolia)and root-bed medium on the treatment of domestic sewage within a vertical flow constructed wetland /

The effect that plants {Typha latifolia) as well as root-bed medium physical and chemical characteristics have on the treatment of primary treated domestic wastewater within a vertical flow constructed wetland system was investigated. Five sets of cells, with two cells in each set, were used. Each cell was made of concrete and measured 1 .0 m X 1 .0 m and was 1.3 m deep. Four different root-bed media were tested : Queenston Shale, Fonthill Sand, Niagara Shale and a Michigan Sand. Four of the sets contained plants and a single type of root-bed medium. The influence of plants was tested by operating a Queenston Shale set without plants. Due to budget constraints no replicates were constructed. All of the sets were operated independently and identically for twenty-eight months. Twelve months of data are presented here, collected after 16 months of continuous operation. Root-bed medium type did not influence BOD5 removal. All of the sets consistently met Ontario Ministry of Environment (MOE) requirements (<25 mg/L) for BOD5 throughout the year. The 12 month average BOD5 concentration from all sets with plants was below 2.36 mg/L. All of the sets were within MOE discharge requirements (< 25 mg/L) for suspended solids with set effluent concentrations ranging from 1.53 to 14.80 mg/L. The Queenston Shale and Fonthill Sand media removed the most suspended solids while the Niagara Shale set produced suspended solids. The set containing Fonthill Sand was the only series to meet MOE discharge requirements (< Img/L) for total phosphorus year-round with a twelve month mean effluent concentration of 0.23 mg/L. Year-round all of the root-bed media were well below MOE discharge requirements (< 20mg/L in winter and < 10 mg/L in sumnner) for ammonium. The Queenston Shale and Fonthill Sand sets removed the most total nitrogen. Plants had no effect on total nitrogen removal, but did influence how nitrogen was cycled within the system. Plants increased the removal of suspended solids by 14%, BOD5 by 10% and total phosphorus by 22%. Plants also increased the amount of dissolved oxygen that entered the system. During the plant growing season removal of total phosphorus was better in all sets with plants regardless of media type. The sets containing Queenston Shale and Fonthill Sand media achieved the best results and plants in the Queenston Shale set increased treatment efficiency for every parameter except nitrogen. Vertical flow wetland sewage treatment systems can be designed and built to consistently meet MOE discharge requirements year-round for BOD5, suspended solids, total phosphorus and ammonium. This system Is generally superior to the free water systems and sub-surface horizontal flow systems in cold climate situations.