Application of [Lambda] to the fourth perturbation theory in calculating the equation of state of rare gas solids and fcc metals

We have calculated the thermodynamic properties of monatomic fcc crystals from the high temperature limit of the Helmholtz free energy. This equation of state included the static and vibrational energy components. The latter contribution was calculated to order A4 of perturbation theory, for a range of crystal volumes, in which a nearest neighbour central force model was used. We have calculated the lattice constant, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the adiabatic and the isothermal bulk modulus, and the Gruneisen parameter, for two of the rare gas solids, Xe and Kr, and for the fcc metals Cu, Ag, Au, Al, and Pb. The LennardJones and the Morse potential were each used to represent the atomic interactions for the rare gas solids, and only the Morse potential was used for the fcc metals. The thermodynamic properties obtained from the A4 equation of state with the Lennard-Jones potential, seem to be in reasonable agreement with experiment for temperatures up to about threequarters of the melting temperature. However, for the higher temperatures, the results are less than satisfactory. For Xe and Kr, the thermodynamic properties calculated from the A2 equation of state with the Morse potential, are qualitatively similar to the A 2 results obtained with the Lennard-Jones potential, however, the properties obtained from the A4 equation of state are in good agreement with experiment, since the contribution from the A4 terms seem to be small. The lattice contribution to the thermal properties of the fcc metals was calculated from the A4 equation of state, and these results produced a slight improvement over the properties calculated from the A2 equation of state. In order to compare the calculated specific heats and bulk moduli results with experiment~ the electronic contribution to thermal properties was taken into account~ by using the free electron model. We found that the results varied significantly with the value chosen for the number of free electrons per atom.

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.

Isolation of DNA sequences with homology to a degenerate FaRP oligonucleotide from the crayfish Procambarus clarkii

The neuropeptide Th1RFamide with the sequence Phe-Met-Arg-Phe-amide was originally isolated in the clam Macrocallista nimbosa (price and Greenberg, 1977). Since its discovery, a large family ofFl\1RFamide-related peptides termed FaRPs have been found to be present in all major animal phyla with functions ranging from modulation of neuronal activity to alteration of muscular contractions. However, little is known about the genetics encoding these peptides, especially in invertebrates. As FaRP-encoding genes have yet to be investigated in the invertebrate Malacostracean subphylum, the isolation and characterization ofFaRP-encoding DNA and mRNA was pursued in this project. The immediate aims of this thesis were: (1) to amplify mRNA sequences of Procambarus clarkii using a degenerate oligonucleotide primer deduced from the common amino acid sequence ofisolated Procambarus FaRPS, (2) to determine if these amplification products encode FaRP gene sequences, and (3) to create a selective cDNA library of sequences recognized by the degenerate oligonucleotide primer. The polymerase chain reaction - rapid amplification of cDNA ends (PCR-RACE) is a procedure in which a single gene-specific primer is used in conjunction with a generalized 3' or 5' primer to amplify copies ofthe region between a single point in the transcript and the 3' or 5' end of cDNA of interest (Frohman et aI., 1988). PCRRACE reactions were optimized with respect to primers used, buffer composition, cycle number, nature ofgenetic substrate to be amplified, annealing, extension and denaturation temperatures and times, and use of reamplification procedures. Amplification products were cloned into plasmid vectors and recombinant products were isolated, as were the recombinant plaques formed in the selective cDNA library. Labeled amplification products were hybridized to recombinant bacteriophage to determine ligated amplification product presence. When sequenced, the five isolated PCR-RACE amplification products were determined not to possess FaRP-encoding sequences. The 200bp, 450bp, and 1500bp sequences showed homology to the Caenorhabditis elegans cosmid K09A11, which encodes for cytochrome P450; transfer-RNA; transposase; and tRNA-Tyr, while the 500bp and 750bp sequences showed homology with the complete genome of the Vaccinia virus. Under the employed amplification conditions the degenerate oligonucleotide primer was observed to bind to and to amplify sequences with either 9 or 10bp of 17bp identity. The selective cDNA library was obselVed to be of extremely low titre. When library titre was increased, white. plaques were isolated. Amplification analysis of eight isolated Agt11 sequences from these plaques indicated an absence of an insertion sequence. The degenerate 17 base oligonucleotide primer synthesized from the common amino acid sequence ofisolated Procambarus FaRPs was thus determined to be non-specific in its binding under the conditions required for its use, and to be insufficient for the isolation and identification ofFaRP-encoding sequences. A more specific primer oflonger sequence, lower degeneracy, and higher melting temperature (TJ is recommended for further investigation into the FaRP-encoding genes of Procambarlls clarkii.

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.

On the formulation and the calculation of the harmonic contributions to the Debye-Waller factor in metals (sodium)

The algebraic expressions for the anharmonic contributions to the Debye-Waller factor up to 0(A ) and 0 L% ) £ where ^ is the scattering wave-vector] have been derived in a form suitable for cubic metals with small ion cores where the interatomic potential extends to many neighbours. This has been achieved in terms of various wave-vector dependent tensors, following the work of Shukla and Taylor (1974) on the cubic anharmonic Helmholtz free energy. The contribution to the various wave-vector dependent tensors from the coulomb and the electron-ion terms in the interatomic metallic potential has been obtained by the Ewald procedure. All the restricted multiple whole B r i l l o u i n zone (B.Z.) sums are reduced to single whole B.Z. sums by using the plane wave representation of the delta function. These single whole B.Z. sums are further reduced to the •%?? portion of the B.Z. following Shukla and Wilk (1974) and Shukla and Taylor (1974). Numerical calculations have been performed for sodium where the Born-Mayer term in the interatomic potential has been neglected because i t is small £ Vosko (1964)3 • *n o^er to compare our calculated results with the experimental results of Dawton (1937), we have also calculated the r a t io of the intensities at different temperatures for the lowest five reflections (110), (200), (220), (310) and (400) . Our calculated quasi-harmonic results agree reasonably well with the experimental results at temperatures (T) of the order of the Debye temperature ( 0 ). For T » © ^ 9 our calculated anharmonic results are found to be in good agreement with the experimental results.The anomalous terms in the Debye-Waller factor are found not to be negligible for certain reflections even for T ^ ©^ . At temperature T yy Op 9 where the temperature is of the order of the melting temperature (Xm) » "the anomalous terms are found to be important almost for all the f i ve reflections.

Petrography, geochemistry and structure of the Timmins area, Ontario

Regional structural analysis of the Timmins area indicates four major periods of tectonic deformation. The DI deformation is characterized by a series of isoclinal FI folds which are outlined in the study area by bedding, pillow tops and variolitic flows. The D2 deformation developed the Porcupine Syncline and refolded the Fl folds about a NE. axis. A pervasive S2 foliation developed during low grade (greenschist) regional metamorphism associated with the D2 deformation. The S2 foliation developed south of the Destor-Porcupine Break. The third phase of tectonic D3 deformation is recognized by the development of a S3 sub-horizontal crenulation cleavage which developed on the plane of the S2 foliation. No meso scopic folds are associated with this deformation. The 8 3 crenulation cleavage is observed south of the Destor-Porcupine Break. The D4 tectonic deformation is recorded as a subvertical S4 crenulation cleavage which developed on the plane of the S2 foliation and also offsets the S3 crenulation cleavage. Macroscopic F4 folds have refolded the F2 axial plane. No metamorphic recrystallization is associated with this deformation. The S4 crenulation cleavage is observed south of the Destor-Porcupine Break. Petrographic evidence indicates that the Timmins area has been subjected to pervasive regional low grade (greenschist) metamorphism which has recrystallized the original mineralogy. South of the study are~ the Donut Lake ultramafic lavas have been subjected to contact medium grade (amphibolite facies) metamorphism associated with the intrusion of the Peterlong Lake Complex. The Archean volcanic rocks of the Timmins area have been subdivided into komatiitic, tholeiitic and calcalkaline suites based on Zr, Ti0 2 and Ni. The three elements were used because of their r e lative immobility during subsequent metamorphic events. Geochemical observations in the Timmins area indicates that the composition of the Goose Lake and Donut Lake Formations are a series of peridotitic, pyroxenitic and basaltic komatiites. The Lower Schumacher Formation is a sequence of basaltic komatiites while the upper part of the Lower Schumacher Formation is an intercalated sequence of basaltic komatiites and low Ti0 2 tholeiites. The variolitic flows are felsic tholeiites in composition and geochemical evidenc e sugg ests that they developed as a n immiscible splitting of a tholeiitic magma. The Upper Schumacher Formation is a sequence of tholeiitic rocks dis p laying a mild iron enrichment. The Krist and Boomerang Formations are the felsic calc-alkaline rocks of the study area which are characteristically pyroclastic. The Redstone Fo rmation is dominantly a calc-alkali ne sequence of volcani c rocks whose minor mafic end me mbers exposed in 1t.he study hav e basaltic komatiitic compositions. Geochemical evidence sugges ts that the Keewatin-type se dimentary rocks have a composition similar to a quartz diorite or a granodiorite. Fi e l d obs ervations and petrographic evidence suggests that they were derived fr om a distal source and now repr esent i n part a turbidite sequence. The Timiskaming-type sedimentary rocks approach the c omp osi t ion of the felsic calc-alkaline rocks of the study area . The basal conglomerate in the study are a sugge s ts that th e uni t was derived fr om a proximal source. Petrographic and ge ochemical evidence suggests that the peridotitic and pyroxenitic komatiites originated as a 35-55% partial melt within the mantle, in excess of 100 Km. depth. The melt ros e as a diapir with the subsequent effusion of the ultramafic lavas, The basaltic komatiites and tholeiitic rocks originated in the mantle from lesser degrees of partial melting and fractionated in low pressure chambers. Geochemical evidence suggests a "genetic link" between the basaltic komatiites and tholeiites, The calc-alkaline rocks developed as a result of the increa.se In PO in the magma chamber. The felsic calcalkaline rocks are a late stage effusion possibly the last major volcanic eruptions in the area.

Petrography, Petrochemistry and petrogenesis of Huronian volcanic rocks of the Elliot Lake region, Ontario

In the Elliot lake region of northern Ontario, Yolcanlc lava piles represent the lowermost units of the Huronian SUpergroup. These rocks general1y trend east-west and belong to the Elliot lake Group. They are s1tuated on the north and south limbs or the QuIrke lake Syncline. The volcanIc rocks of this study contain a secondary minerai assemblage consisting of actinolite, biotite, chlorIte, eptdote/cllnozoislte tttanomagnettte and calcite characteristic of greenschist metamorphism. Compilation of data suggests that metamorphism of the volcanic rocks proceeded between 325- and 425-C and between 2.4 and 4.7 kb. Geochemtcally these lavas represent tholeiitic and calc-alkaline assemblages. The tholeiites are character1sttcally enriched tn Fe and Tt and consist mainly of basalts, basaltic andesites and andesites. These rocks are believed to have formed by the partIal melting of a peridottte source at low P-T. In contrast, the calc-alkaline rocks are depleted in Fe and TI, but show a signIficant enrichment In 51 and Zr; andesIte Is the major rock type for thIs assemblage. I·t Is postUlated that the calc-alkalIne sU1te of rocks was the result of eIther the partIal meltIng of abasaltic·magma at shallow depth, or the melttng of s1al1c crustal materIal due to the added we1ght of tholeiitIc material on an unstable crust and to downwarplng processes Inttlated by convection cells.

Sedimentology and stratigraphy of the Lower Silurian Grimsby formation, in subsurface, Lake Erie and Southwestern Ontario

Since the first offshore Lake Erie well was drilled in 1941, the Grimsby and Thorold formations of the Cataract Group have been economically important to the oil and gas industry of Ontario. The Cataract Group provides a significant amount of Ontario's gas production primarily from wells located on Lake Erie. The Grimsby - Thorold formations are the result of nearshore estuarine processes influenced by tides on a prograding shelf and are composed of subtidal channel complexes, discrete tidal channels, mud flats and non-marine deposits. Deposition was related to a regressive - transgressive cycle associated with eustatic sea level changes caused by the melting and resurgence of continental glaciation centred in Africa in the Late Ordovician/Early Silurian. Grimsby deposition began during a regression with the deposition of subtidal channel complexes incised into the marine deposits of the Cabot Head Formation. The presence of mud drapes and mud couplets suggest that these deposits were influenced by tides. These deposits dominate the lower half of the Grimsby. Deposition continued with a change from these subtidal channel complexes to laterally migrating, discrete, shallow tidal channels and mud flats. These were in turn overlain by the non-marine deposits of the Thorold Formation. Grimsby - Thorold deposition ended with a major transgression replacing siliciclastic deposition with primarily carbonate deposition. Sediment was sourced from the east and southeast and associated with a continuation of the Taconic Orogeny into the Early Silurian. The fluvial head of the estuary prograded from a shoreline that was located in western New York and western Pennsylvania running NNE-SSW and then turning NW-SE and paralleling the present day Lake Erie shoreline. iii The facies attributed to the Grimsby - Thorold formations can be ascribed to the three zones within the tripartite zonation suggested by Dalrymple et ale (1992) for estuaries, that is, a marine-dominated facies, a mixed energy facies, and a facies that is dominated by fluvial processes. Also, sediments within the Grimsby - Thorold are commonly fining upwards sequences which are common in estuarine settings whereas deltaic deposits are normally composed of coarsening upwards sequences in a vertical wedge shape with coarser material near the head. The only coarsening observed was in the Thorold Formation and attributed to non-marine deposition by palynological evidence. The presence of a lag deposit at the base of the sediments of the Grimsby Thorold formations suggests that they were incised into the Cabot Head Formation. Further, the thickness of Early Silurian sediments located between the top of the Queenston Formation, where Early Silurian sedimentation began, to the top of the Reynales - Irondequoit formation are constant whether the Grimsby - Thorold formations are present or not. Also, cross-sections using a sand body located in the Cabot Head Formation for correlation further imply that the Grimsby Formation has been incised into the previous deposits of the Cabot Head.

Geochemistry and petrogenesis of the McElroy and Larder Lake assemblages, Abitibi Greenstone Belt, Northeastern Ontario

The McElroy and Larder Lake assemblages, located in the southern Abitibi Greenstone Belt are two late Archean metavolcanic sequences having markedly contrasting physical characteristics arid are separated from one another by a regional fault. An assemblage is an informal term which describes stratified volcanic and/or sedimentary rock units built during a specific time period in a similar depositional or volcanic setting and are commonly bounded by faults, unconformities or intrusions. The petrology and petrogenesis of these assemblages have been investigated to determine if a genetic link exists between the two adjacent assemblages. The McElroy assemblage is homoclinal sequence of evolved massive and pillowed fl.ows, which except for the basal unit represents a progressively fractionated volcanic pile. From the base to the top of the assemblage the lithologies include Fe-tholeiitic, dendritic flows; komatiite basaltic, ultramafic flows; Mg-tholeiitic, leucogabbro; Mg-tholeiitic, massive flows and Fe-tholeiitic, pillowed flows. Massive flows range from coarse grained to aphanitic and are commonly plagioclase glomerophyric. The Larder Lake assemblage consists of komatiitic, Mg-rich and Fe-rich tholeiitic basalts, structurally disrupted by folds and faults. Tholeiitic rocks in the Larder Lake assemblage range from aphanitic to coarse grained massive and pillowed flows. Komatiitic flows contain both spinifex and massive textures. Geochemical variability within both assemblages is attributed to different petrogenetic histories. The lithologies of the McElroy assemblage were derived by partial melting of a primitive mantle source followed by various degrees of crystal fractionation. Partial melting of a primitive mantle source generated the ultramafic flows and possibly other flows in the assemblage. Fractionation of ultramafic flows may have also produced the more evolved McElroy lithologies. The highly evolved, basal, dendritic flow may represent the upper unit 3 of a missing volcanic pile in which continued magmatism generated the remaining McElroy lithologies. Alternatively, the dendritic flows may represent a primary lava derived from a low degree (10-15%) partial melt of a primitive mantle source which was followed by continued partial melting to generate the ultramafic flows. The Larder Lake lithologies were derived by partial melting of a komatiitic source followed by gabbroic fractionation. The tectonic environment for both assemblages is interpreted to be an oceanic arc setting. The McElroy assemblage lavas were generated in a mature back arc setting whereas the Larder Lake lithologies were produced during the early stages of komatiitc crust subduction. This setting is consistent with previous models involving plate tectonic processes for the generation of other metavolcanic assemblages in the Abitibi Greenstone Belt.

Physical volcanology, sedimentology, stratigraphy and petrochemistry of the Berry Creek metavolcanics: an Archean calc-alkaline complex, Lake of the Woods, Ontario

The steeply dipping, isoclinally folded early Precambrian (Archean) Berry Creek Metavolcanic Complex comprises primary to resedimented pyroclastic, epiclastic and autoclastic deposits. Tephra erupted from central volcanic edifices was dumped by mass flow mechanisms into peripheral volcanosedimentary depressions. Sedimentation has been essentially contemporaneous with eruption and transport of tephra. The monolithic to heterolithic tuffaceous horizons are interpreted as subaerial to subaqueous pumice and ash flows, secondary debris flows, lahars, slump deposits and turbidites. Monolithic debris flows, derived from crumble breccia and dcme talus, formed during downslope collapse and subsequent gravity flowage. Heterolithic tuff, lahars and lava flow morphologies suggest at least temporary emergence of the edifice. Local collapse may have accompanied pyroclastic volcanism. The tephra, produced by hydromagmatic to magmatic eruptions, were rapidly transported, by primary and secondary mechanisms, to a shallow littoral to deep water subaqueous fan developed upon the subjacent mafic metavolcanic platform. Deposition resulted from traction, traction carpet, and suspension sedimentation from laminar to turbulent flows. Facies mapping revealed proximal (channel to overbank) to distal facies epiclastics (greywackes, argillite) intercalated with proximal vent to medial fan facies crystal rich ash flows, debris flows, bedded tuff and shallow water to deep water lava flows. Framework and matrix support debris flows exhibit a variety of subaqueous sedimentary structures, e.g., coarse tail grading, double grading, inverse to normal grading, graded stratified pebbly horizons, erosional channels. Pelitic to psammitic AE turbidites also contain primary stru~tures, e.g., flames, load casts, dewatering pipes. Despite low to intermediate pressure greenschist to amphibolite grade metamorphism and variably penetrative deformation, relicts of pumice fragments and shards were recognized as recrystallized quartzofeldspathic pseudomorphs. The mafic to felsic metavolcanics and metasediments contain blasts of hornblende, actinolite, garnet, pistacitic epidote, staurolite, albitic plagioclase, and rarely andalusite and cordierite. The mafic metavolcanics (Adams River Bay, Black River, Kenu Lake, Lobstick Bay, Snake Bay) display _holeiitic trends with komatiitic affinities. Chemical variations are consistent with high level fractionation of olivine, plagioclase, amphibole, and later magnetite from a parental komatiite. The intermediate to felsic (64-74% Si02) metavolcanics generally exhibit calc-alkaline trends. The compositional discontinuity, defined by major and trace element diversity, can be explained by a mechanism involving two different magma sources. Application of fractionation series models are inconsistent with the observed data. The tholeiitic basalts and basaltic andesites are probably derived by low pressure fractionation of a depleted (high degree of partial melting) mantle source. The depleted (low Y, Zr) calc-alkaline metavolcanics may be produced by partial melting of a geochemically evolved source, e.g., tonalitetrondhjemite, garnet amphibolite or hydrous basalt.

Phonon spectra and temperature variation of thermodynamic properties of fcc metals via Finnis-Sinclair type many body potentials: Sutton-Chen and improved Sutton-Chen models

Volume(density)-independent pair-potentials cannot describe metallic cohesion adequately as the presence of the free electron gas renders the total energy strongly dependent on the electron density. The embedded atom method (EAM) addresses this issue by replacing part of the total energy with an explicitly density-dependent term called the embedding function. Finnis and Sinclair proposed a model where the embedding function is taken to be proportional to the square root of the electron density. Models of this type are known as Finnis-Sinclair many body potentials. In this work we study a particular parametrization of the Finnis-Sinclair type potential, called the "Sutton-Chen" model, and a later version, called the "Quantum Sutton-Chen" model, to study the phonon spectra and the temperature variation thermodynamic properties of fcc metals. Both models give poor results for thermal expansion, which can be traced to rapid softening of transverse phonon frequencies with increasing lattice parameter. We identify the power law decay of the electron density with distance assumed by the model as the main cause of this behaviour and show that an exponentially decaying form of charge density improves the results significantly. Results for Sutton-Chen and our improved version of Sutton-Chen models are compared for four fcc metals: Cu, Ag, Au and Pt. The calculated properties are the phonon spectra, thermal expansion coefficient, isobaric heat capacity, adiabatic and isothermal bulk moduli, atomic root-mean-square displacement and Gr\"{u}neisen parameter. For the sake of comparison we have also considered two other models where the distance-dependence of the charge density is an exponential multiplied by polynomials. None of these models exhibits the instability against thermal expansion (premature melting) as shown by the Sutton-Chen model. We also present results obtained via pure pair potential models, in order to identify advantages and disadvantages of methods used to obtain the parameters of these potentials.

(A) Photoregulation of DNA Functions by Cyclic Azobenzene-tethered Oligonucleotides (B) Site-specific Fluorescent Labeling of DNA using Inverse Electron Demand Diels-Alder Reaction between trans-Cyclooctene Derivatives and BODIPY-Tetrazine Adducts

(A) Most azobenzene-based photoswitches require UV light for photoisomerization, which limit their applications in biological systems due to possible photodamage. Cyclic azobenzene derivatives, on the other hand, can undergo cis-trans isomerization when exposed to visible light. A shortened synthetic scheme was developed for the preparation of a building block containing cyclic azobenzene and D-threoninol (cAB-Thr). trans-Cyclic azobenzene was found to thermally isomerize back to the cis-form in a temperature-dependent manner. cAB-Thr was transformed into the corresponding phosphoramidite and subsequently incorporated into oligonucleotides by solid phase synthesis. Melting temperature measurement suggested that incorporation of cis-cAB into oligonucleotides destabilizes DNA duplexes, these findings corroborate with circular dichroism measurement. Finally, Fluorescent Energy Resonance Transfer experiments indicated that trans-cAB can be accommodated in DNA duplexes. (B) Inverse Electron Demand Diels-Alder reactions (IEDDA) between trans-olefins and tetrazines provide a powerful alternative to existing ligation chemistries due to its fast reaction rate, bioorthogonality and mutual orthogonality with other click reactions. In this project, an attempt was pursued to synthesize trans-cyclooctene building blocks for oligonucleotide labeling by reacting with BODIPY-tetrazine. Rel-(1R-4E-pR)-cyclooct-4-enol and rel-(1R,8S,9S,4E)-Bicyclo[6.1.0]non-4-ene-9-ylmethanol were synthesized and then transformed into the corresponding propargyl ether. Subsequent Sonogashira reactions between these propargylated compounds with DMT-protected 5-iododeoxyuridine failed to give the desired products. Finally a methodology was pursued for the synthesis of BODIPY-tetrazine conjugates that will be used in future IEDDA reactions with trans-cyclooctene modified oligonucleotides.

The Ligand and Membrane-binding Behaviour of the Phosphatidylinositol Transfer Proteins (PITPa & PITPb)

Human Class I phosphatidylinositol transfer proteins (PITPs) exists in two forms: PITPα and PITPβ. PITPs are believed to be lipid transfer proteins based on their capacity to transfer either phosphatidylinositol (PI) or phosphatidylcholine (PC) between membrane compartments in vitro. In Drosophila, the PITP domain is found to be part of a multi-domain protein named retinal degeneration B (RdgBα). The PITP domain of RdgBα shares 40 % sequence identity with PITPα and has been shown to possess PI and PC binding and transfer activity. The detailed molecular mechanism of ligand transfer by the human PITPs and the Drosophila PITP domain remains to be fully established. Here, we investigated the membrane interactions of these proteins using dual polarization interferometry (DPI). DPI is a technique that measures protein binding affinity to a flat immobilized lipid bilayer. In addition, we also measured how quickly these proteins transfer their ligands to lipid vesicles using a fluorescence resonance energy transfer (FRET)-based assay. DPI investigations suggest that PITPβ had a two-fold higher affinity for membranes compared to PITPα. This was reflected by a four-fold faster ligand transfer rate for PITPβ in comparison to PITPα as determined by the FRET assay. Interestingly, DPI analysis also demonstrated that PI-bound human PITPs have lower membrane affinity compared to PC-bound PITPs. In addition, the FRET studies demonstrated the significance of membrane curvature in the ligand transfer rate of PITPs. The ligand transfer rate was higher when the accepting vesicles were highly curved. Furthermore, when the accepting vesicles contained phosphatidic acid (PA) which have smaller head groups, the transfer rate increased. In contrast, when the accepting vesicles contained phosphoinositides which have larger head groups, the transfer rate was diminished. However, PI, the favorite ligand of PITPs, or the presence of anionic lipids did not appear to influence the ligand transfer rate of PITPs. Both DPI and FRET examinations revealed that the PITP domain of RdgBα was able to bind to membranes. However, the RdgBα PITP domain appears to be a poor binder and transporter of PC.

Les gènes et protéines BSP chez la souris et l’humain : clonage, caractérisation, expression sous forme recombinante et étude des fonctions biologiques

L’infertilité affecte environ 15% des couples en âge de se reproduire. Dans près de la moitié des cas, des facteurs masculins sont à la base de l’infertilité, quoique les causes exactes demeurent souvent inconnues. Les spermatozoïdes de mammifères subissent une série d’étapes de maturation avant d’acquérir la capacité de féconder un ovocyte. Les premiers changements ont lieu à l’intérieur de l’épididyme, où les spermatozoïdes gagnent la capacité de se mouvoir ainsi que de reconnaître et d’interagir avec l’ovocyte. Suite à l’éjaculation, ils doivent subir une seconde série de modifications à l’intérieur du tractus génital femelle, nommée capacitation. Nous avons préalablement démontré que chez le bovin, la famille de protéines BSP (Binder of SPerm) est essentielle à la capacitation. Des homologues des BSP ont aussi été isolés du fluide séminal de porc, de bouc, de bélier, de bison et d’étalon. Malgré la détection d’antigènes apparentés aux BSP dans le fluide séminal de souris et d’humain, les homologues des BSP n’ont jamais été caractérisés chez ces espèces. Nous avons émis l’hypothèse que des homologues des BSP seraient exprimés chez la souris et l’humain et joueraient un rôle dans la maturation des spermatozoïdes. Nous avons démontré que des séquences homologues aux BSP sont présentes dans les génomes murin et humain. Le génome murin contient trois séquences; Bsph1, Bsph2a et Bsph2b, tandis qu’une seule séquence (BSPH1) a été identifée chez l’humain. Les séquences d’ADNc de Bsph1, Bsph2a et BSPH1 ont été clonées, tandis que Bsph2b serait probablement un pseudogène. Les trois gènes sont exprimés uniquement dans l’épididyme et font partie d’une sous-famille distincte à l’intérieur de la famille des BSP. Chez les ongulés, les BSP sont exprimées par les vésicules séminales, sont ajoutées aux spermatozoïdes lors de l’éjaculation et représentent une proportion significative des protéines du plasma séminal. Au contraire, les BSP épididymaires ne sont retrouvées qu’en faibles quantités dans le fluide séminal. L’étude de leur rôle dans les fonctions spermatiques était donc plus difficile que chez les ongulés, où l’isolement des protéines natives du plasma séminal à l’aide de techniques de chromatographie était possible. Afin d’étudier sa fonction, nous avons exprimé BSPH1 recombinante dans E. coli. Les ponts disulfure des domaines de type-II caractéristiques de ces protéines ont fait en sorte que l’expression de BSPH1 fusionnée à une étiquette hexahistidine ou glutathion-S-transférase a donné lieu à des protéines insolubles dans les corps d’inclusion. La production de BSPH1 soluble a été possible grâce à l’ajout d’une étiquette thiorédoxine et l’expression dans une souche au cytoplasme oxidatif. BSPH1 a été purifiée par affinité et sa liaison aux partenaires connus des BSP, la phosphatidylcholine, les lipoprotéines de faible densité et la membrane des spermatozoïdes, suggérait que la protéine recombinante possédait sa conformation native et pouvait être utilisée pour des essais fonctionnels. La forme native de BSPH1 a été détectée dans le plasma séminal humain suite au fractionnement par gel filtration. La liaison de BSPH1 native à une colonne d’affinité à l’héparine a indiqué qu’elle partage aussi cette propriété de liaison avec la famille des BSP, et pourrait lier les GAGs semblables à l’héparine du tractus génital féminin. Une colonne d’immunoaffinité anti-BSPH1 a été préparée à l’aide d’anticorps générés contre des protéines recombinantes, et a permis d’isoler BSPH1 native à partir d’extraits de spermatozoïdes humains. Nos résultats montrent que BSPH1 native serait localisée dans les microdomaines « rafts » de la membrane. Sa masse moléculaire apparente était de 32 kDa, ce qui est supérieur à la masse prédite selon sa séquence en acides aminés, indiquant la présence probable de modifications post-traductionnelles, ou d’une migration anormale. L’effet de BSPH1 recombinante et des anticorps anti-BSPH1 sur la motilité, la viabilité et la capacitation a aussi été étudié. Les deux dernières variables ont été mesurées par un essai de cytométrie en flux, optimisé dans cette étude. Aucun effet des protéines recombinantes ou des anticorps sur la motilité et la viabilité des spermatozoïdes n’a été noté. Quoiqu’une stimulation modeste, quoique significative, de la capacitation ait été observée à la plus faible concentration de BSPH1, les concentrations plus élevées n’ont pas montré d’effet. De la même manière, les anticorps anti-BSPH1 n’ont pas eu d’effet significatif sur la capacitation. Ces résultats suggèrent que BSPH1 produite dans E. coli n’affecte pas la capacitation de façon marquée. Cependant, puisque BSPH1 native possède probablement des modifications post-traductionnelles, une protéine recombinante produite dans des cellules de mammifères pourrait affecter les fonctions spermatiques. De manière alternative, les BSP épididymaires remplissent peut-être un rôle différent dans les fonctions spermatiques que celles sécrétées par les vésicules séminales des ongulés. Les résultats décrits dans cette thèse pourraient contribuer à améliorer le diagnostic de l’infertilité masculine, ainsi que les techniques de reproduction assistée et éventuellement, pourraient mener au développement de contraceptifs masculins.