Prevalence of overweight and obesity in children from urban and rural Niagara

Purpose: The influence of environment in the development of overweight and obesity is an ongoing concern. This investigation examined the influence of urbanization on the rates of childhood overweight and obesity. Method: 2167 (1090M, 1077F) grade four children from 75 schools in Ontario's Niagara Region were sampled. A sophisticated algorithm overlaying electoral boundaries, population densities, and the knowledge of community members was used to classify schools into one of three location categories: urban {N= 1588), urban fringe {N= 379), and rural (A^= 234). Each subject was measured for: height, weight, and aerobic performance (Leger). Physical activity was evaluated with the self-report Participation Questionnaire (free-time and organized sport activities), and teacher's evaluations of student activity. Overweight (overweight and obesity combined) was measured both as a continuous (BMI) and categorical variable (BMI category), to evaluate the prevalence by location. A multivariate analysis was used to test for a suppression effect. Results: BMI and BMI category did not differ significantly by location or gender, and no evidence of a gender interaction existed. According to both a linear and logistic regression, physical activity or fitness levels did not suppress the influence of location on BMI and BMI category. Age, gender, free-time activity, organized sports, fitness level, and number of siblings, were all found to significantly influence overweight. Conclusions: It is plausible that the prevalence of overweight does not differ in urban and rural children from the Niagara Region. Further investigation is recommended, examining subjects by individual location of residence, in multiple regions throughout Ontario.

An experimental investigation of dechanneling in copper single crystals

This investigation comprises a comparison of experimental and theoretical dechanneling of MeV protons in copper single crystals. Dechanneling results when an ion's transverse energy increases to the value where the ion can undergo small impact parameter collisions with individual atoms. Depth dependent dechanneling rates were determined as functions of lattice temperature, ion beam energy and crystal axis orientation. Ion beam energies were IMeV and 2MeV,temperatures ranged from 35 K to 280 K and the experiment was carried out along both the (lOa) and <110) axes. Experimental data took the form of aligned and random Rutherford backscattered energy spectra. Dechanneling rates were extracted from these spectra using a single scattering theory that took explicit account of the different stopping powers experienced by channeled and dechanneled ions and also included a correction factor to take into account multiple scattering effects along the ion's trajectory. The assumption of statistical equilibrium and small angle scattering of the channeled ions allows a description of dechanneling in terms of the solution of a diffusion like equation which contains a so called diffusion function. The diffusion function is shown to be related to the increase in average transverse energy. Theoretical treatments of increase in average transverse energy due to collisions of projectiles with channel electrons and thermal perturbations in the lattice potential are reviewed. Using the diffusion equation and the electron density in the channel centre as a fitting parameter dechanneling rates are extracted. Excellent agreement between theory and experiment has been demonstrated. Electron densities determined in the fitting procedure appear to be realistic. The surface parameters show themselves to be good indicators of the quality of the crystal.

Behavioural, mensural, and ecological correlates of polygyny in the eastern meadowlark (Sturnella magna)

During 1982 and 1983 I studied male attributes and attributes of the territory of male Eastern Meadowlarks (Sturnella magna) in order to determine whether there was a correlation between any of the attributes investigated and the number of females attracted by a male. Seventeen males, nine of which were polygynous and eight monogamous, were studied in 1982 and sixteen males.of which .. seven were polygynous and nine ~onogamous, were studied in 1983. The study was conducted in Short Hills Park, 10 km southwest of St. Catharines, Ontario and was designed to compare two hypotheses: the "sexy son" hypothesis (Weatherhead and Robertson,1977) and the polygyny threshold model (Verner and Willson,1966, Orians, 1969). Male attributes investigated were male size and song behaviour. Six measures of male size were taken: weight, flattened and natural wing chord length, culmen length, bill depth and length of the tarsometatarsus. In 1983 song repertoire size and song versatility measures were investigated. Attributes of the territory studied were: territory size, density of plant stems, percentage plant cover and measures of vegetation structure. In 1983 Arthropods were collected from each territory and sorted according to taxonomic group and size. During 1983, territory attributes were sampled twice, once early and once later in the nesting season. Analysis of data involved univariate comparisons between monogamous and polygynous males using T-tests and multivariate comparisons were made using discriminant function analysis (DFA) and principle components analysis (PCA).No correlations were found between the number of females attracted with, .ny measure of male size or with me, .sures of song versatili or size of song repertoire. Also no correlation was found between terri size and the number of females nesting on a terri . Some attributes of the male's terri id distinguish between monog,mous and po s males of thistudy. Analysis of Arthropod numbers showed that e~ .eran counts were significantly great~r on polygynous territories, a1 the total numb~rs of Arthropods collected showed no s fico .nt differences between territories of monogamous and po males. DFA chose ear teran and Hymenopteran counts as multivariate discriminators; both variables we' e more vegetation revealed that there were no univariate differences between the two groups of males fOT 1982 stem densities, but ~ spp. and Solidago spp. were chosen DFA as multivariate discriminators. The total number of plant stems and of Vicia spp. stems were s ficantly the early 1983 ing on monogamous territories for however DFA found no multivariate discriminators" Variables concerned with the overall aspects of vegetation structure showed significant differences between territories of monogamous and polygynous males. DFA of the 1982 sampling of vegetation structure showed significantly greater mat depth and vegetation height on polygynous territories, a finding which was not supported, however, by peA. For the early 1983 sampling period, plant height was greater on polygynous territories. Multivariate analysis identified greater green cover on polygynous territories, greater ground cover on monogamous territories, and greater depth of mat material on monogamous territories as discriminators between territories of monogamous and polygynous males. A DFA on the major variables of the study showed no significant difference between the territories of monogamous and polygynous male Meadowlarks. Of the correlations found, some were for non-prey Arthr~ods, for cover plants with very small samples sizes, or for variables which were greater for monogamous males during one sampling period and polygynous males during the next. While multivariate discriminators were found, peA showed no grouping of monogamous or polygynous males according to any of the variables investigated. On the basis of the univariate and multivariate analysis of major variables, I concluded that there were no correlations between the number of females attracted with male attributes and no unambiguous correlation with attributes of the territory. My study does not unequivocally support either the "sexy son" or the polygyny threshold hypothesis.

Ecological interactions between zooplankton and photosynthetic bacteria in Crawford Lake, Ontario

The present study was carried out to test the hypothesis that photosynthetic bacteria contribute a large portion of the food of filter feeding zooplankton populations in Crawford Lake, Ontario. The temporal and spatial variations of both groups of organisms are strongly dependent on one another. 14 By using C-Iabelled photosynthetic bacteria. the ingestion and clearance rates of Daphnia pulex, ~. rosea, and Keratella spp were estimated during summer and fall of 1982. These quantitative estimations of zooplankton ingestion and clearence rates on photosynthetic bacteria comprised an original addition to the literature. Photosynthetic bacteria comprised a substantial portion of the diet of all four dominant zooplankton species. The evidence for this is based on the ingestion and clearance rates of the dominant zooplankton species. Ingestion rates of D. pulex and D. rosea ranged 5 5 -1 -1 - -- 5 - -- 5 from 8.3X10 -1 to 14.6XlO -1 cells.ind. hr and 8.1X10 to 13.9X10 cells.ind. hr • Their clearance rates ranged from 0.400 to 1.000 -1 -1 -1 -1 ml.ind. hr. and 0.380 to 0.930 ml.ind. hr • The ingestion and clearance -1 -1 -1 -1 rates of Keratella spp were 600 cell.ind. hr and 0.40 ul.ind. hr respectively. Clearance rates were inversely proportional to the concentration of food cells and directly proportional to the body size of the animals. It is believed that despite the very short reg~neration times of photosynthetic bacteria (3-8 hours) their population densities were controlled in part by the feeding rates of the dominant zooplankton in Crawford Lake. By considering the regeneration times of photosynthetic bacteria and the population clearance rates of zooplankton, it was estimated that between 16 to 52% and 11 to 35% of the PHotosynthetic bacteria were' consumed· by Daphnia· pulex. and Q.. rosea per day. The temporal and spatial distribution of Daphnia pulex, !.. rosea, Keratella quadrata, K. coChlearis and photosynthetic bacteria in Crawford Lake were also investigated during the period of October, 1981 to December, 1982. The photosynthetic bacteria in the lake, constituted a major food source for only those zooplankton Which tolerate anaerobic conditions. Changes in temperature and food appeared to correlate with the seasonal changes in zooplankton density. All four dominant species of zooplankton were abundant at the lake's surface (O-4m) during winter and spring and moved downwards with the thermocline as summer stratification proceeded. Photosynthetic bacteria formed a 2 m thick layer at the chemocline. The position of this photosynthetic bacterial J-ayer changed seasonally. In the summer, the bacterial plate moved upwards and following fall mixing it moved downwards. A vertical shift of O.8m (14.5 to 15.3m) was recorded during the period of June to December. The upper limit of the photosynthetic bacteria in the water column was controlled by dissolved oxygen, and sulfide concentrations While their lower limit was controlled by light intensity. A maximum bacterio- 1 chlorophyll concentration of 81 mg Bchl.l was recorded on August 9, 1981. The seasonal distribution of photosynthetic bacteria was controlledinpart' by ·theg.-"z1ai'_.Q;~.zoopl. ank:tCm;-.Qther -ciactors associated with zooplankton grazing were oxygen and sulfide concentrations.

Effects of intracerebroventricular bombesin administration on local cerebral glucose utilization in the restrained and unrestrained rat

Intracerebroventricular (ICV) administration of bombesin (BN) induces a syndrome characterized by stereotypic locomotion and grooming, hyperactivity and sleep elimination, hyperglycemia and hypothermia, hyperhemodynamics, feeding inhibition, and gastrointestinal function changes. Mammalian BN-like peptides (MBNs), e.g. gastrin-releasing peptide (GRP), Neuromedin C (NMC), and Neuromedin B (NMB), have been detected in the central nervous system. Radio-labeled BN binds to specific sites in discrete cerebral regions. Two specific BN receptor subtypes (GRP receptor and NMB receptor) have been identified in numerous brain regions. The quantitative 2-[14C]deoxyglucose ([14C]20G) autoradiographic method was used to map local cerebral glucose utilization (LCGU) in the rat brain following ICV injection of BN (vehicle, BN O.1Jlg, O.5Jlg). At each dose, experiments were conducted in freely moving or restrained conditions to determine whether alterations in cerebral function were the result of BN central administration, or were the result of BN-induced motor stereotypy. The anteroventral thalamic nucleus (AV) (p=O.029), especially its ventrolateral portion (AVVL) (pdensities of BN receptors, the internal granular layer of the olfactory bulb (IGr) (p=O.028), and the suprachiasmatic nucleus (SCh) (p=O.003) exhibited BN treatment effects. BN effects on LCGU were also observed in the median eminence (ME) (p=O.011). Restraint, however, decreased LCGU in the lateral dorsal thalamic nucleus, ventrolateral and dorsomedial parts (LOVL and LOOM) (p=O.044, p=O.009), and the lateral geniculate (LG) (p=O.027). In sum, BN induced a marked and highly localized alteration in cerebral metabolism within parts of the anterior thalamus, which is the principle relay in the limbic circuitry. BN effects were also observed in IGr, Mi, SCh, and ME. Effects of restraint were found in LOVL, LOOM, and LG. It is suggested that increased LCGU in AV and AVVL may be the result of functional change in the limbic circuitry and the hypothalamus caused by BN receptor functional modification. In IGr, increased LCGU following BN administration is considered to be mainly the result the activation of NMB receptor, a subtype of BN receptors. In SCh, increased LCGU is believed to be caused both by BN effects on the thalamic, the hypothalamic, and the limbic functions and by activation of GRP receptor, another BN receptors subtype found in SCh. In ME, increased LCGU is suggested to be caused by BN effects on the hypothalamic functions, especially those related to the neuroendocrine functions. None of the alterations seen in these regions reflects the emission of stereotyped motor behaviors. Rather, they reflect a direct influence of BN central administration upon functioning of the cerebral regions influenced by BN administration. The restraint effects seen in LO, including LOOM and LOVL, are suggested to be the result of altered behavioral expression. The restraint effects seen in LG is suggested to be the result of reduced locomotion.

Fimbriae of coprinus cinereus, schizophyllum commune and ustilago violacea = structural aspects and a role in conjugation

Surface fibrils (fimbriae) have been observed on fungi from every major group. Fimbriae are thought to be involved in the following cell to cell interactions: conjugation, flocculation and adhesion. Several higher fungi exibit two other types of interactions: hyphal fusion (anastomosis) and clamp connection formation. As a prelude to examining the role of fimbriae in these processes, the fimbriae of two fungi that undergo these fusion events were examined. Electron microscopy studies revealed that Coprinus cinereus and Schizophyllum commune are fimbriated. C. cinereus fimbriae were 5 nm in diameter and 0.5 to 20 11m in length. Fimbriae of C. cinereus oidia were more numerous and longer than those of the hyphal stage. S. commune fimbriae were also 5 nm in diameter, but were only 0.5 to 2 11m in length. There was an unequal distribution of fimbriae on the hyphal surfaces of S. commune . Fimbriae were sparsely distributed over the entire hyphal surface, with higher densities of fibrils present on the side growths of the hyphae found in the older sections of the mycelium. Antiserum raised against Ustilago violacea fimbrial protein (AU) crossreacted strongly with 37 and 39 kd C. cinereus mycelial proteins. In contrast, AU bound very weakly to 89 and 92 kd S. commune mycelial proteins. Since AU cross-reacted poorly with S. commune fimbrial proteins, it was impossible to further characterize the fimbriae of this specIes. The 37 and 39 kd C. cinereus proteins, were isolated by electroelution and were shown to be able to form fibrils the same diameter as oidial fimbriae. The 37 kd protein was shown to be composed of several proteins with isoelectric points ranging from pH 6.1 to 7.63. Furthermore, the 37 kd protein was found to be multimeric, while the 39 kd protein was not. These results strongly suggested that the 37 kd protein is the structural fimbrial protein of C. cine reus . Finally, a series of experiments were designed to determine whether fimbriae are required for conjugation in U. violacea Conjugation was inhibited significantly with AU, but not with pre-immune serum or AU preincubated with purified fimbrial protein. Thus, it was concluded that fimbriae play a central role in mating in this organism.

Gypsy moths (Lymantria dispar) in the Niagara Region : population density variation, introduction of an entomopathogenic fungus (Entomophaga maimaiga) and occurrence of nuclear polyhedrosis virus

The gypsy moth, Lymantria dispar, a major defoliator of broad leaf trees, was accidentally introduced into North America in 1869. Much interest has been generated regarding the potential of using natural pathogens for biological control of this insect. One of these pathogens, a highly specific fungus, Entomophaga maimaiga, was accredited with causing major epizootics in populations of gypsy moth across the north-eastern United States in 1989 and 1990 and is thought to be spreading northwards into Canada. This study examined gypsy moth population densities in the Niagara Region. The fungus, .E.. maimaiga, was artificially introduced into one site and the resulting mortality in host populations was noted over two years. The relationship between fungal mortality, host population density and occurrence of another pathogen, the nuclear polyhedrosis virus (NPV), was assessed. Gypsy moth population density was assessed by counting egg masses in 0.01 hectare (ha) study plots in six areas, namely Louth, Queenston, Niagara-on-the-Lake, Shorthills Provincial Park, Chippawa Creek and Willoughby Marsh. High variability in density was seen among sites. Willoughby Marsh and Chippawa Creek, the sites with the greatest variability, were selected for more intensive study. The pathogenicity of E. maimaiga was established in laboratory trials. Fungal-infected gypsy moth larvae were then released into experimental plots of varying host density in Willoughby Marsh in 1992. These larvae served as the inoculum to infect field larvae. Other larvae were injected with culture medium only and released into control plots also of varying host density. Later, field larvae were collected and assessed for the presence of .E.. maimaiga and NPV. A greater proportion of larvae were infected from experimental plots than from control plots indicating that the experimental augmentation had been successful. There was no relationship between host density and the proportion of infected larvae in either experimental or control plots. In 1992, 86% of larvae were positive for NPV. Presence and intensity of NPV infection was independent of fungal presence, plot type or interaction of these two factors. Sampling was carried out in the summer of 1993, the year after the introduction, to evaluate the persistence of the pathogen in the environment. Almost 50% of all larvae were infected with the fungus. There was no difference between control and experimental plots. Data collected from Willoughby Marsh indicated that there was no correlation between the proportion of larvae infected with the fungus and host population density in either experimental or control plots. About 10% of larvae collected from a nearby site, Chippawa Creek, were also positive for .E.. maimaiga suggesting that low levels of .E.. maimaiga probably occurred naturally in the area. In 1993, 9.6% of larvae were positive for NPV. Again, presence or absence of NPV infection was independent of fungal presence plot type or interaction of these two factors. In conclusion, gypsy moth population densities were highly variable between and within sites in the Niagara Region. The introduction of the pathogenic fungus, .E.. maimaiga, into Willoughby Marsh in 1992 was successful and the fungus was again evident in 1993. There was no evidence for existence of a relationship between fungal mortality and gypsy moth density or occurrence of NPV. The results from this study are discussed with respect to the use of .E.. maimaiga in gypsy moth management programs.

Individual differences and factors affecting male behaviour in field crickets

Individual differences in male sexual behav~our and the factors influencing calling behaviour were studied in the field crickets Gryllus 2 integer and Q. veletis. In a large (13m) outdoor arena individually numbered adult male ~~ integer started calling at three to five days of age but thereafter the age of individual G. integer males did not affect nightly calling duration. Calling also did not correlate with individual weight. In this study individual male calling was continuously distributed from 0 hrs. per night to 3.5 hrs. per night, on average. A temporal effect on the number of G. integer males calling was observed. The number of males calling through the night was uniform, but a sharp increase in the number calling was observed in the early morning. No difference in calling times was observed between the night and dawn callers. AlsC)' males calling at dawn usually didnotc'all during the preceeding night. Calling and reproductive success in 1979 demonstrated a negative logarithmic relationship while in the 1980(initial) population a negative linear relationship was observed. No relationship was seen in the 1980 high density population. The ratio of non-callers to callers also affected the mating of individuals in the 1979 and1980(initial) densities:-non~callers (males calling .5 hrs. per night, on average, or less) obtained more females when the population contained a high number of callers, this being a negative logarithmic relationship to, No such relationship was observed in the 1980 high density population. Individual displacement varied nightly and was not correlated to amount of calling or reproductive success of individual G. integer males. G. integer males were displa~ed more when in a higher density in the outdoor arena Male G. integer and G. veletis behaviours were also observed in an indoor arena at different densities and, in G. veletis, with respect to female presence. When females were present in the arena, in G. veletis, male calling was reduced. Males of both species called less, on average, when in ~ higher density, than when they were in a lower density. Male displacement of both species increased on average when in a higher density as compared to displacement in a lower density. Aggression was measured by aggressive call-ing and fighting and was studied in regards to density.G. integer demonstrated less aggression in all but one comparison at higher density. No difference was observed in the ratio of aggressive calling to f.ighting comparison in G. integer. G. veletis demonstrated mixed results. No difference in aggression between densities was observed in comparisons. Less.aggression did occur in higher densities when comparisons invol.ved fighting behaviour. Male behaviour represents a competitive strategy against ot~er males, strategy being defined as a genetic (in part) alternative to other strategies. In this sense, the factors of time, density, male-male aggression, and female presence are conditions demonstrated to affect male behaviour in G. integer and G. veletis. Individual male differences and other considerations suggest that alternative male behaviours are represented by at least two conditional strategies. This possibility, and the transient 'or stable nature of genetic polymorphisms in field cricket behaviour are considered.

Male reproductive competition in the field crickets : gryllus veletis and G. pennsylvanicus

Sexual behavior in the field crickets, Gryllus veletis and G. pennsylvanicus , was studied in outdoor arenas (12 m2) at high and low levels of population density in 1983 and 1984. Crickets were weighed, individually marked, and observed from 2200 until 0800 hrs for at least 9 continuous nights. Calling was measured at 5 min intervals, and movement and matings were recorded hourly. Continuous 24 hr observations were also conducted,·and occurrences of aggressive and courtship songs were noted. The timing of males searching, calling, courting, and fighting for females should coincide with female movement and mating patterns. For most samples female movement and matings occurred at night in the 24 hr observations and were randomly distributed with time for both species in the 10 hr observations. Male movement for G. veletis high density only was enhanced at night in the 24 hr observations, however, males called more at night in both species at high and low densities. Male movement was randomly distributed with time in the 10 hr observations, and calling increased at dawn for the G. pennsylvanicus 1984 high density sample, but was randomly distributed in other samples. Most courtship and aggression songs in the 24 hr observations were too infrequent for statistical testing and generally did not coincide with matings. Assuming residual reproductive value, and costs attached to a male trait in terms of future reproductive success decline with age, males should behave in more costly ways with age; by calling and moving more with age. Consequently, mating rates should increase with age. Female behavior may not change with age. G. veletis , females moved more with age at both low density samples, however, crickets moved less with age at high density. G. pennsylvanicus females moved more with age in the 1984 low density sample, whereas crickets moved less with age in the 1983 high density sample. For both species males in the 1984 high density samples called less with age. For G. pennsylvanicus in 1983 calling and mating rates increased with age. Mating rates decreased with age for G. veletis males in the high density sample. Aging may not affect cricket behavior. As population density increases fewer calling sites become available, costs of territoriality increase, and matings resulting from non-calling behavior should increase. For both species the amount of calling and in G. veletis the distance travelled per night was not different between densities. G. pennsylvanicus males and females moved more at low density. At the same deneity levels there were no differences in calling, mating, and, movement rates in G. veletis , however, G. pennsylvanicus males moved more at high density in 1983 than 1984. There was a positive relationship between calling and mating for the G. pennsylvanicus low density sample only, and selection was acting directly to increase calling. For both species no relationships between movement and mating success was found, however, the selection gradient on movement in the G. veletis high density population was significant. The intensity of selection was not significant and was probably due to the inverse relationship between displacement and weight. Larger males should call more, mate more, and move less than smaller males. There were no correlations between calling and individual weight, and an inverse correlation between movement and size in the G. veletis high density population only. In G. pennsylvanicus , there was a positive correlation between individual weight and mating, but, some correlate of weight was under counter selection pressure and-prevented significance of the intensity of selection. In contrast, there was an inverse correlation in the G.·veletis low density B sample. Both measures of selection intensities were significant and showed that weight only was under selection pressures. An inverse correlation between calling and movement was found for G. veletis at low density only. Because males are territorial, females are predicted to move more than males, however, if movement is a mode of male-male reproductive competition then males may move more than females. G. pennsylvanicus males moved more than females in all samples, however, G. veletis males and females moved similar distances at all densities. The variation in relative mating success explained by calling scores, movement, and weight for both species and all samples were not significant In addition, for both species and all samples the intensity of selection never equalled the opportunity for selection.

Mating behaviours of the field cricket (Gryllus integer) at different levels of male density

The reproductive behaviour of the field cricket, Gryllus integer, was systematically observed in indoor arenas to determine the extent of female Choice and male-male competition at different sex ratios representing two male densities (12:6 and 6:6). The costs and benefits to males and females in those two densities were analyzed according to the theory of the evolution o£ leks. Observations were conducted during the dark hours when most calling occurred since hourly rates of courtship song and mating did not fluctuate significantly over a 24 h period. Female mating rates were not significantly different between densities, therefore males at high densities were not advantaged because of increased female tendencies to mate when social stimulation was increased. Mean rates of acoustical signalling (calling and courtin"g) did not differ significantly between densities. Mean rates of fighting by males at the high density were significantly greater than those of males at the low density. Mating benefits associated with callin~courting and fighting were measured. Mating rates did not vary with rates of calling at either density. Calling was not a prerequisite to mating. Courtship song preceded all matings. There was a significant power fit between male mating and courting rates, and male mating and fighting rates at the low, but not at the high, density. Density differences in the benefits associated with increased courting and fighting may relate, in part, to greater economic defensibility and monopoly of females due to reduced male competition at the low density. Dominant males may be preferentially chosen by females or better able to monopolize mating opportunities than subordinate males. Three criteria were used to determine whether dominant males were preferentially chosen by females. The number of matings by males who won fights (within 30 min of mating) was significantly greater than the number of matings by males who were defeated in such fights. Mating rates did not vary significantly with rates of winning at either density. There was a significant power fit between male mating rates and the percentage of fights a male won (irrespective of his fighting-frequency) at the low density. The mean duration a male guarded the female after mating did not vary significantly between densities. There was a significant linear relationship between the duration a spermatophore was retained and the duration a male guarded the female after mating. Courtship song apparently stimulated spermatophore removal. Male guarding involved inter-male aggression and reduced courtship attempts by other males. Males at the high density received no apparent reproductive benefits associated with increased social stimulation. Conclusive evidence for preferential choice of males by females, using the criteria examined here, is lacking. Males at the lower density had fewer competitors and could monopolize females more effectively.

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.

Application of Reptation Quantum Monte Carlo and Related Methods to LiH

This work investigates mathematical details and computational aspects of Metropolis-Hastings reptation quantum Monte Carlo and its variants, in addition to the Bounce method and its variants. The issues that concern us include the sensitivity of these algorithms' target densities to the position of the trial electron density along the reptile, time-reversal symmetry of the propagators, and the length of the reptile. We calculate the ground-state energy and one-electron properties of LiH at its equilibrium geometry for all these algorithms. The importance sampling is performed with a single-determinant large Slater-type orbitals (STO) basis set. The computer codes were written to exploit the efficiencies engineered into modern, high-performance computing software. Using the Bounce method in the calculation of non-energy-related properties, those represented by operators that do not commute with the Hamiltonian, is a novel work. We found that the unmodified Bounce gives good ground state energy and very good one-electron properties. We attribute this to its favourable time-reversal symmetry in its target density's Green's functions. Breaking this symmetry gives poorer results. Use of a short reptile in the Bounce method does not alter the quality of the results. This suggests that in future applications one can use a shorter reptile to cut down the computational time dramatically.

Generalization of the BTK Theory to the Case of Finite Quasiparticle Lifetimes

A generalization to the BTK theory is developed based on the fact that the quasiparticle lifetime is finite as a result of the damping caused by the interactions. For this purpose, appropriate self-energy expressions and wave functions are inserted into the strong coupling version of the Bogoliubov equations and subsequently, the coherence factors are computed. By applying the suitable boundary conditions to the case of a normal-superconducting interface, the probability current densities for the Andreev reflection, the normal reflection, the transmission without branch crossing and the transmission with branch crossing are determined. Accordingly the electric current and the differential conductance curves are calculated numerically for Nb, Pb, and Pb0.9Bi0.1 alloy. The generalization of the BTK theory by including the phenomenological damping parameter is critically examined. The observed differences between our approach and the phenomenological approach are investigated by the numerical analysis.