The effect of counterpoise correction and relaxation energy term to the internal rotation barriers: application to the BF3 ...NH3 and C2H4 ...SO2 dimers

The relevance of the fragment relaxation energy term and the effect of the basis set superposition error on the geometry of the BF3⋯NH3 and C2H4⋯SO2 van der Waals dimers have been analyzed. Second-order Møller-Plesset perturbation theory calculations with the d95(d,p) basis set have been used to calculate the counterpoise-corrected barrier height for the internal rotations. These barriers have been obtained by relocating the stationary points on the counterpoise-corrected potential energy surface of the processes involved. The fragment relaxation energy can have a large influence on both the intermolecular parameters and barrier height. The counterpoise correction has proved to be important for these systems

Counterpoise-corrected geometries and harmonic frequencies of N-body clusters: application to (HF)n (n=3,4)

A study was conducted on the methods of basis set superposition error (BSSE)-free geometry optimization and frequency calculations in clusters larger than a dimer. In particular, three different counterpoise schemes were critically examined. It was shown that the counterpoise-corrected supermolecule energy can be easily obtained in all the cases by using the many-body partitioning of energy

A quantum molecular similarity analysis of changes in molecular electron density caused by basis set flotation and electric field application

Quantum molecular similarity (QMS) techniques are used to assess the response of the electron density of various small molecules to application of a static, uniform electric field. Likewise, QMS is used to analyze the changes in electron density generated by the process of floating a basis set. The results obtained show an interrelation between the floating process, the optimum geometry, and the presence of an external field. Cases involving the Le Chatelier principle are discussed, and an insight on the changes of bond critical point properties, self-similarity values and density differences is performed

Relations among several nuclear and electronic density functional reactivity indexes

A set of connections among several nuclear and electronic indexes of reactivity in the framework of the conceptual Density Functional Theory by using an expansion ofthe energy functional in terms of the total number of electrons and the normal coordinates within a canonical ensemble was derived. The relations obtained provided explicit links between important quantities related to the chemical reactivity of a system. This paper particularly demonstrates that the derivative of the electronic energy with respect to the external potential of a system in its equilibrium geometry was equal to the negative of the nuclear repulsion derivative with respect to the external potential

Experimental and computational study on the synthesis and characterization of self-assembling meso/macroporous materials in highly concentrated emulsions

Report for the scientific sojourn carried out at Massachusetts General Hospital Cancer Center-Harvard Medical School, Estats Units, from 2010 to 2011. The project aims to study the aggregation behavior of amphiphilic molecules in the continuous phase of highly concentrated emulsions, which can be used as templates for the synthesis of meso/macroporous materials. At this stage of the project, we have investigated the self-assembly of diblock and triblock surfactants under the effect of a confined geometry being surrounded by the droplets of the dispersed phase. These droplets limit the growth of the aggregates, deeply modify their orientation and hence alter their spatial arrangement as compared to the self-assembly taking place far enough from any boundary surface, that is in the bulk. By performing Monte Carlo simulations, we have showed that the interface between the dispersed and continuous phases as well as its shape has a significant impact on the structural order of the resulting aggregates and hence on the potential applications of highly concentrated emulsions as reaction media, drug delivery systems, or templates for meso/macroporous materials. Due to the combined effect of symmetry breaking and morphological frustration, very intriguing structures, such as square columnar liquid crystals, twisted X-shaped aggregates, and helical phases of cylindrical aggregates, never observed in the bulk for the same model surfactant, have been found. The presence of other more conventional structures, such as micelles and cubic and hexagonal liquid crystals, formed at low and high amphiphilic concentrations, respectively, further enhance the interest on this already rich aggregation behavior.

Dose reconstruction in the context of tomotherapy

In vivo dosimetry is a way to verify the radiation dose delivered to the patient in measuring the dose generally during the first fraction of the treatment. It is the only dose delivery control based on a measurement performed during the treatment. In today's radiotherapy practice, the dose delivered to the patient is planned using 3D dose calculation algorithms and volumetric images representing the patient. Due to the high accuracy and precision necessary in radiation treatments, national and international organisations like ICRU and AAPM recommend the use of in vivo dosimetry. It is also mandatory in some countries like France. Various in vivo dosimetry methods have been developed during the past years. These methods are point-, line-, plane- or 3D dose controls. A 3D in vivo dosimetry provides the most information about the dose delivered to the patient, with respect to ID and 2D methods. However, to our knowledge, it is generally not routinely applied to patient treatments yet. The aim of this PhD thesis was to determine whether it is possible to reconstruct the 3D delivered dose using transmitted beam measurements in the context of narrow beams. An iterative dose reconstruction method has been described and implemented. The iterative algorithm includes a simple 3D dose calculation algorithm based on the convolution/superposition principle. The methodology was applied to narrow beams produced by a conventional 6 MV linac. The transmitted dose was measured using an array of ion chambers, as to simulate the linear nature of a tomotherapy detector. We showed that the iterative algorithm converges quickly and reconstructs the dose within a good agreement (at least 3% / 3 mm locally), which is inside the 5% recommended by the ICRU. Moreover it was demonstrated on phantom measurements that the proposed method allows us detecting some set-up errors and interfraction geometry modifications. We also have discussed the limitations of the 3D dose reconstruction for dose delivery error detection. Afterwards, stability tests of the tomotherapy MVCT built-in onboard detector was performed in order to evaluate if such a detector is suitable for 3D in-vivo dosimetry. The detector showed stability on short and long terms comparable to other imaging devices as the EPIDs, also used for in vivo dosimetry. Subsequently, a methodology for the dose reconstruction using the tomotherapy MVCT detector is proposed in the context of static irradiations. This manuscript is composed of two articles and a script providing further information related to this work. In the latter, the first chapter introduces the state-of-the-art of in vivo dosimetry and adaptive radiotherapy, and explains why we are interested in performing 3D dose reconstructions. In chapter 2 a dose calculation algorithm implemented for this work is reviewed with a detailed description of the physical parameters needed for calculating 3D absorbed dose distributions. The tomotherapy MVCT detector used for transit measurements and its characteristics are described in chapter 3. Chapter 4 contains a first article entitled '3D dose reconstruction for narrow beams using ion chamber array measurements', which describes the dose reconstruction method and presents tests of the methodology on phantoms irradiated with 6 MV narrow photon beams. Chapter 5 contains a second article 'Stability of the Helical TomoTherapy HiArt II detector for treatment beam irradiations. A dose reconstruction process specific to the use of the tomotherapy MVCT detector is presented in chapter 6. A discussion and perspectives of the PhD thesis are presented in chapter 7, followed by a conclusion in chapter 8. The tomotherapy treatment device is described in appendix 1 and an overview of 3D conformai- and intensity modulated radiotherapy is presented in appendix 2. - La dosimétrie in vivo est une technique utilisée pour vérifier la dose délivrée au patient en faisant une mesure, généralement pendant la première séance du traitement. Il s'agit de la seule technique de contrôle de la dose délivrée basée sur une mesure réalisée durant l'irradiation du patient. La dose au patient est calculée au moyen d'algorithmes 3D utilisant des images volumétriques du patient. En raison de la haute précision nécessaire lors des traitements de radiothérapie, des organismes nationaux et internationaux tels que l'ICRU et l'AAPM recommandent l'utilisation de la dosimétrie in vivo, qui est devenue obligatoire dans certains pays dont la France. Diverses méthodes de dosimétrie in vivo existent. Elles peuvent être classées en dosimétrie ponctuelle, planaire ou tridimensionnelle. La dosimétrie 3D est celle qui fournit le plus d'information sur la dose délivrée. Cependant, à notre connaissance, elle n'est généralement pas appliquée dans la routine clinique. Le but de cette recherche était de déterminer s'il est possible de reconstruire la dose 3D délivrée en se basant sur des mesures de la dose transmise, dans le contexte des faisceaux étroits. Une méthode itérative de reconstruction de la dose a été décrite et implémentée. L'algorithme itératif contient un algorithme simple basé sur le principe de convolution/superposition pour le calcul de la dose. La dose transmise a été mesurée à l'aide d'une série de chambres à ionisations alignées afin de simuler la nature linéaire du détecteur de la tomothérapie. Nous avons montré que l'algorithme itératif converge rapidement et qu'il permet de reconstruire la dose délivrée avec une bonne précision (au moins 3 % localement / 3 mm). De plus, nous avons démontré que cette méthode permet de détecter certaines erreurs de positionnement du patient, ainsi que des modifications géométriques qui peuvent subvenir entre les séances de traitement. Nous avons discuté les limites de cette méthode pour la détection de certaines erreurs d'irradiation. Par la suite, des tests de stabilité du détecteur MVCT intégré à la tomothérapie ont été effectués, dans le but de déterminer si ce dernier peut être utilisé pour la dosimétrie in vivo. Ce détecteur a démontré une stabilité à court et à long terme comparable à d'autres détecteurs tels que les EPIDs également utilisés pour l'imagerie et la dosimétrie in vivo. Pour finir, une adaptation de la méthode de reconstruction de la dose a été proposée afin de pouvoir l'implémenter sur une installation de tomothérapie. Ce manuscrit est composé de deux articles et d'un script contenant des informations supplémentaires sur ce travail. Dans ce dernier, le premier chapitre introduit l'état de l'art de la dosimétrie in vivo et de la radiothérapie adaptative, et explique pourquoi nous nous intéressons à la reconstruction 3D de la dose délivrée. Dans le chapitre 2, l'algorithme 3D de calcul de dose implémenté pour ce travail est décrit, ainsi que les paramètres physiques principaux nécessaires pour le calcul de dose. Les caractéristiques du détecteur MVCT de la tomothérapie utilisé pour les mesures de transit sont décrites dans le chapitre 3. Le chapitre 4 contient un premier article intitulé '3D dose reconstruction for narrow beams using ion chamber array measurements', qui décrit la méthode de reconstruction et présente des tests de la méthodologie sur des fantômes irradiés avec des faisceaux étroits. Le chapitre 5 contient un second article intitulé 'Stability of the Helical TomoTherapy HiArt II detector for treatment beam irradiations'. Un procédé de reconstruction de la dose spécifique pour l'utilisation du détecteur MVCT de la tomothérapie est présenté au chapitre 6. Une discussion et les perspectives de la thèse de doctorat sont présentées au chapitre 7, suivies par une conclusion au chapitre 8. Le concept de la tomothérapie est exposé dans l'annexe 1. Pour finir, la radiothérapie «informationnelle 3D et la radiothérapie par modulation d'intensité sont présentées dans l'annexe 2.

Petrology and chemistry of the Chalten Plutonic Complex and implications on the magmatic and tectonic evolution of the Southernmost Andes (Patagonia) during the Miocene

Résumé Scientific:Pétrologie et Géochimie du Complexe Plutonique de Chaltén et les conséquences pour l'évolution magmatique et tectonique du Andes du Sud (Patagonia) pendant le MiocèneLe sujet de cette thèse est le Complexe Plutonique de Chaltén (CHPC), situé à la frontière entre le Chili et l'Argentine, en Patagonie (49°15'S). Ce complexe s'est mis en place au début du Miocène, dans un contexte de changements tectoniques importants. La géométrie et la vitesse de migration des plaques en Patagonie a été modifiée suite l'ouverture de la plaque Farallon il y a 25Ma (Pardo-Casas and Molnar 1987) et la subduction de la ride active du Chili sous la plaque sud-américaine il y a 14Ma (Cande and Leslie 1986). Les effets de cette reconfiguration tectonique sur la morphologie et le magmatisme de la plaque supérieure sont encore sujets à discussion. Dans ce contexte, un groupe d'intrusions miocènes - telle que le CHPC - est particulièrement intriguant, car en position transitionnelle entre le batholithe patagonien et l'arc volcanique cénozoïque et récent à l'ouest, et les laves de plateau de Patagonie à l'est (Fig. 1). A cause de leur position tectonique transitoire, ces plutons isolés hors du batholithe représentent un endroit clé pour comprendre les interactions entre la tectonique à large échelle et le magmatisme en Patagonie. Ici, je présente de nouvelles données de terrain, petrologiques, géochimiques et géochronologiques dans le but de caractériser la nature du CHPC, qui était largement inconnu avant cette étude, dans le but de tester l'hypothèse de migration de l'arc et erosion par subduction.Les résultats de l'investigation géochimique (chapitre 2) montrent que le CHPC n'est qu'un exemple parmi les plutons isolés d'arrière arc ave une composition calco-alcaline caractéristique, c-à-d une signature d'arc. La plupart de ces plutons isolés ont une composition alcaline. Le CHPC, contrairement, a une signature calco-alcaline avec Κ intermédiaire, tel que le batholithe patagonien et la plupart des roches volcaniques quaternaires liées à l'arc le long des Andes.De nouvelles données géochronologiques U-Pb de haute précision sur des zircons, acquis par TIMS, sur le CHPC donnent des âges entre 17.0 et 16.4Ma. Les âges absolus sont en accord avec la séquence intrusive déduite des relations de terrain (chapitre 1). Ces données sont les premières contraintes d'âge U-Pb sur le CHPC. Elles montrent clairement que l'histoire magmatique du CHPC n'a pas de lien direct avec la subduction de la ride à cette latitude (Cande and Leslie 1986), car le complexe est au moins 6Ma plus ancien.Une comparaison en profondeur avec les autres intrusions d'âge Miocène en Patagonie révèlent - pour la première fois - une évolution temporelle intéressante. Il y a une tendance E-W distincte au magmatisme calco-alcalin entre 20-16Ma avec une diminution de l'âge vers l'est - le CHPC est l'expression la plus orientale de cette tendance. Je suggère que la relation espace-temps reflète une migration vers l'est (vers le continent) de l'arc magmatique. Je propose que le facteur principal contrôlant cette migration est la subduction rapide suite à la reconfiguration de la vitesse des plaques tectoniques après l'ouverture la plaque Farallon (à ~26Ma) qui résulterait en une déformation importante ainsi qu'à des taux élevés d'érosion dans la fosse de subduction.Les rapports d'isotopes radiogéniques (Pb, Sr, Nd) élevés, une signature 6018 basse et un rapport Th/La élevé sont des paramètres distinctifs pour les roches mafiques du CHPC. Le modèle isotopique présenté (chapitre 2) suggère que cette signature reflète une contamination de la source, dans le coin de manteau, plutôt qu'une contamination crustale. La signature des éléments en trace du CHPC indiquent que le coin de manteau a été contaminé par des composés terrigènes, le plus vraisemblablement par des sédiments paléozoïques.Les travaux de terrain, la pétrographie et la géothermobarométrie ont été utilisés dans le but de comprendre l'histoire interne du CHPC (chapitre 3). Ces données suggèrent deux niveaux distincts de cristallisation : l'un dans la croûte moyenne (6 à 4.5kbar) et l'autre à un niveau peu profond (3.5 à 2kbar). La modélisation isotopique AFC de la contamination crustale indique des taux variables d'assimilation, qui ne sont pas corrélés avec le degré de différenciation. Cela suggère que différents volumes de magma se sont différenciés en profondeur, de façon indépendante. Cela implique que le CHPC se serait formés en plusieurs puises de magmas provenant d'au moins trois sources différentes. Les textures des granodiorites et des granites indiquent des teneurs élevées en cristaux avant la mise en place et, par conséquent, des températures d'emplacement faibles. Les observations de terrain montrent que les roches mafiques sont déformées, alors que ce n'est pas le cas pour les granodiorites et granites (plus jeunes). La déformation des roches mafiques est encore sujet de recherche, afin de savoir si elle est liée à la déformation régionale en régime compressif ou à l'emplacement lui-même. Cependant, la mise en place de grand volume de magma felsique riche en cristaux suggère un régime d'extension.Scientific Abstract:Petrology and chemistry of the Chaltén Plutonic Complex and implications on the magmatic and tectonic evolution of the Southernmost Andes (Patagonia) during the MioceneThe subject of this thesis is the Chaltén Plutonic Complex (CHPC) located at the frontier between Chile and Argentina in Patagonia (at 49° 15 'Southern latitude). This complex intruded during early Miocene in a context of major tectonics changes. The plate geometry of Patagonia has been modified by changes in the plate motions after the break up of the Farallôn plate at 25Ma (Pardo-Casas and Molnar 1987) and by the subduction of the Chile spreading Ridge beneath South-America at 14 Ma (Cande and Leslie 1986). The effects of this tectonic setting on the morphology and the magmatism of the overriding plate are a matter of on-going discussion. Particularly intriguing in this context is a group of isolated Miocene intrusions - like the CHPC - which are located in a transitional position between the Patagonian Batholith and the Cenozoic and Recent volcanic arc in the West, and the Patagonian plateau lavas in the East (Fig. 1). Due to their transient tectonic position these isolated plutons outside the batholith represent a key to understanding the interaction between global-scale tectonics and magmatism in Patagonia. Here, I present new field, penological, geochemical and geochronological data to characterize the nature of the CHPC, which was largely unknown before this study, in order to test the hypothesis of time- transgressive magmatism.The results of the geochemical investigation (Chapter 2) show that the CHPC is only one among these isolated back-arc plutons with a characteristic calc-alkaline composition, i.e. arc signature. Most of these isolated intrusives have an alkaline character. The CHPC, in contrast, has a medium Κ calc-alkaline signature, like the Patagonian batholith and most of the Quaternary arc-related volcanic rocks along the Andes.New high precision TIMS U-Pb zircon dating of the CHPC yield ages between 17.0 to 16.4 Ma. The absolute ages support the sequence of intrusion relations established in the field (Chapter 1). These data are the first U-Pb age constraints on the CHPC, and clearly show that the magmatic history of CHPC has no direct link to the subduction of the ridge, since this complex is at least 6 Ma older than the time of collision of the Chile ridge at this latitude (Cande and Leslie 1986).An in-depth comparison with other intrusion of Miocene age in Patagonia reveals - for the first time - an interesting temporal pattern. There is a distinct E-W trend of calc-alkaline magmatism between 20-16 Ma with the younging of ages in the East - the CHPC is the easternmost expression of this trend. I suggest that this time-space relation reflects an eastward (landward) migration of the magmatic arc. I propose that main factor controlling this migration is the fast rates of subduction after the major reconfigurations of the plate tectonic motions after the break up of the Farallôn Plate (at -26 ) resulting in strong deformation and high rates of subduction erosion.High radiogenic isotope ratios (Pb, Sr, Nd) ratios, low 5018 signature and high Th/La ratios in mafic rocks are distinctive features of the CHPC. The presented isotopic models (Chapter 2) suggest that this signature reflects source contamination of the mantle wedge rather than crustal contamination. The trace element signature of the CHPC indicates that the mantle wedge was contaminated with a terrigenous component, most likely from Paleozoic sediments.Fieldwork, petrography and geothermobarometry were used to further unravel the internal history of the CHPC (Chapter 3). These data suggest two main levels of crystallization: one a mid crustal levels (6 to 4.5 kbar) and other a shallow level (3.5 to 2 kbar). Isotopic AFC modeling of crustal contamination indicate variable rates of assimilation, which are not correlated with the degree of differentiation. This suggests that different batches of magma differentiate independently at depths. This implies that the CHPC would have formed by several pulses of magmas from at least 3 different sources. Textures of granodiorites and granites indicate a high content of crystals previous to the emplacement and consequently low emplacement temperatures. Field observations show that the mafic rocks are deformed, whereas the (younger) granodiorites and granites are not. It is still subject of investigation whether the deformation of the mafic rocks is related to regional deformation during a compressional regime or to the emplacement it self. However, the emplacement of huge amount of crystal rich felsic magmas suggests an extensional regime.Résumé Grand PublicPétrologie et Géochimie du Complexe Plutonique de Chaltén et les conséquences pour l'évolution magmatique et tectonique du Andes du Sud (Patagonia) pendant le MiocèneLe Complexe Plutonique de Chaltén (CHPC) est un massif montagneux situé à 49°S à la frontière entre le Chili et l'Argentine, en Patagonie (région la plus au sud de l'Amérique du Sud). Il est composé de montagnes qui peuvent atteindre plus de 3000 mètres d'altitude, telles que le Cerro Fitz Roy (3400m) et le Cerro Torre (3100m). Ces montagnes sont composées de roches plutoniques, c.-à-d. des magmas qui se sont refroidis et ont cristallisés sous la surface terrestre.La composition chimique de ces roches montre que les magmas, qui ont formé ce complexe plutonique, font partie d'un volcanisme d'arc. Celui-ci se forme lorsqu'une plaque océanique plonge sous une plaque continentale. Les géologues appellent ce processus « subduction ». Dans un tel scénario, le manteau terrestre, qui se fait prendre entre ces deux plaques, fond pour former ainsi du magma. Ce magma remonte à travers la plaque continentale vers la surface. Si celui-ci atteint la surface, il forme les roches volcaniques, comme par exemple des laves. S'il n'atteint pas la surface, le magma se refroidit pour former finalement les roches plutoniques.Le long de la marge ouest d'Amérique du Sud, la plaque Nazca - qui se situe au sud-est de la plaque océanique pacifique - passe en dessous de la plaque d'Amérique du Sud. La bordure ouest du sud de la plaque sud-américaine a également été affectée par d'autres processus tectoniques, tels que des changements dramatiques dans les déplacements de plaques (il y a 25Ma) et la collision de la ride du Chili (depuis 15 Ma jusqu'à aujourd'hui). Ces caractéristiques tectoniques et magmatiques font de cette région un haut lieu pour les géologues. La plaque Nazca, s'est formée suite à l'ouverture d'une plaque océanique plus ancienne, il y a 25Ma. Cette ouverture est liée aux vitesses de subduction les plus rapides jamais connues. La ride du Chili est l'endroit où le sol de l'Océan Pacifique s'ouvre, formant deux plaques océaniques : les plaques Nazca et Antarctique. La ride du Chili subducte sous la plaque sud-américaine depuis 15Ma, en association avec la formation de grands volumes de magma ainsi que des changements morphologiques importants. La question de savoir lequel de ces changements tectoniques globaux affecte la géologie et la géographie de Patagonie a été, et est encore, discutée pendant de nombreuses années. De nombreux chercheurs suggèrent que la plupart des caractéristiques morphologiques et magmatiques en Patagonie sont liés à la subduction de la ride du Chili, mais cette suggestion est encore débattue comme le montre notre étude.Le batholithe de Patagonie du sud (SPB) est un énorme massif composé de roches plutoniques et il s'étend tout au long de la côte ouest de Patagonie (au sud de 47°S). Ces roches correspondent certainement aux racines d'un ancien arc volcanique, qui a été soulevé et érodé. Le CHPC, ainsi que d'autres petites intrusions dans la région, se situe dans une position exotique, à 100km à l'est du SPB. Certains chercheurs suggèrent que ces intrusions pourraient être liées à la subduction de la ride du Chili.Afin de débattre de cette problématique, nous avons utilisé différentes méthodes géochronologiques pour déterminer l'âge du CHPC et le comparer (a) à l'âge des roches intrusives similaires du SPB et (b) à l'âge de la collision de la ride du Chili. Dans ce travail, nous prouvons que le CHPC s'est formé au moins 7Ma avant la collision avec la ride du Chili. Sur la base des âges du CHPC et de la composition chimique de ses roches et minéraux, nous proposons que le CHPC fait partie d'un arc volcanique ancien. La migration de l'arc volcanique plus profondément dans le continent résulte de la grande vitesse de subduction entre 25 et lOMa. Des caractéristiques évidentes pour un tel processus - telles qu'une déformation importante et une vitesse d'érosion élevée - peuvent être rencontrées tout au long de la bordure ouest de l'Amérique du sud.

Morphodynamic analysis of cerebral aneurysm pulsation from time-resolved rotational angiography

This paper presents a technique to estimate and model patient-specific pulsatility of cerebral aneurysms over onecardiac cycle, using 3D rotational X-ray angiography (3DRA) acquisitions. Aneurysm pulsation is modeled as a time varying-spline tensor field representing the deformation applied to a reference volume image, thus producing the instantaneousmorphology at each time point in the cardiac cycle. The estimated deformation is obtained by matching multiple simulated projections of the deforming volume to their corresponding original projections. A weighting scheme is introduced to account for the relevance of each original projection for the selected time point. The wide coverage of the projections, together with the weighting scheme, ensures motion consistency in all directions. The technique has been tested on digital and physical phantoms that are realistic and clinically relevant in terms of geometry, pulsation and imaging conditions. Results from digital phantomexperiments demonstrate that the proposed technique is able to recover subvoxel pulsation with an error lower than 10% of the maximum pulsation in most cases. The experiments with the physical phantom allowed demonstrating the feasibility of pulsation estimation as well as identifying different pulsation regions under clinical conditions.

Méthodes sismiques actives adaptées à l'étude de milieux hétérogènes : application aux glissements de terrain

Résumé Les glissements de terrain représentent un des principaux risques naturels dans les régions montagneuses. En Suisse, chaque année les glissements de terrains causent des dégâts qui affectent les infrastructures et ont des coûts financiers importants. Une bonne compréhension des mécanismes des glissements peut permettre d'atténuer leur impact. Celle-ci passe notamment par la connaissance de la structure interne du glissement, la détermination de son volume et de son ou ses plans de glissement. Dans un glissement de terrain, la désorganisation et la présence de fractures dans le matériel déplacé engendre un changement des paramètres physiques et en particulier une diminution des vitesses de propagation des ondes sismiques ainsi que de la densité du matériel. Les méthodes sismiques sont de ce fait bien adaptées à l'étude des glissements de terrain. Parmi les méthodes sismiques, l'analyse de la dispersion des ondes de surface est une méthode simple à mettre en oeuvre. Elle présente l'avantage d'estimer les variations des vitesses de cisaillement avec la profondeur sans avoir spécifiquement recours à l'utilisation d'une source d'onde S et de géophones horizontaux. Sa mise en oeuvre en trois étapes implique la mesure de la dispersion des ondes de surface sur des réseaux étendus, la détermination des courbes de dispersion pour finir par l'inversion de ces courbes. Les modèles de vitesse obtenus à partir de cette procédure ne sont valides que lorsque les milieux explorés ne présentent pas de variations latérales. En pratique cette hypothèse est rarement vérifiée, notamment pour un glissement de terrain dans lequel les couches remaniées sont susceptibles de présenter de fortes hétérogénéités latérales. Pour évaluer la possibilité de déterminer des courbes de dispersion à partir de réseaux de faible extension des mesures testes ont été effectuées sur un site (Arnex, VD) équipé d'un forage. Un profil sismique de 190 m de long a été implanté dans une vallée creusée dans du calcaire et remplie par des dépôts glacio-lacustres d'une trentaine de mètres d'épaisseur. Les données acquises le long de ce profil ont confirmé que la présence de variations latérales sous le réseau de géophones affecte l'allure des courbes de dispersion jusqu'à parfois empêcher leur détermination. Pour utiliser l'analyse de la dispersion des ondes de surface sur des sites présentant des variations latérales, notre approche consiste à déterminer les courbes de dispersions pour une série de réseaux de faible extension, à inverser chacune des courbes et à interpoler les différents modèles de vitesse obtenus. Le choix de la position ainsi que de l'extension des différents réseaux de géophones est important. Il tient compte de la localisation des hétérogénéités détectées à partir de l'analyse de sismique réfraction, mais également d'anomalies d'amplitudes observées sur des cartes qui représentent dans le domaine position de tir - position du récepteur, l'amplitude mesurée pour différentes fréquences. La procédure proposée par Lin et Lin (2007) s'est avérée être une méthode efficace permettant de déterminer des courbes de dispersion à partir de réseaux de faible extension. Elle consiste à construire à partir d'un réseau de géophones et de plusieurs positions de tir un enregistrement temps-déports qui tient compte d'une large gamme de distances source-récepteur. Au moment d'assembler les différentes données une correction de phase est appliquée pour tenir compte des hétérogénéités situées entre les différents points de tir. Pour évaluer cette correction nous suggérons de calculer pour deux tir successif la densité spectrale croisée des traces de même offset: Sur le site d'Arnex, 22 courbes de dispersions ont été déterminées pour de réseaux de géophones de 10 m d'extension. Nous avons également profité du forage pour acquérir un profil de sismique verticale en ondes S. Le modèle de vitesse S déduit de l'interprétation du profil de sismique verticale est utilisé comme information à priori lors l'inversion des différentes courbes de dispersion. Finalement, le modèle en deux dimension qui a été établi grâce à l'analyse de la dispersion des ondes de surface met en évidence une structure tabulaire à trois couches dont les limites coïncident bien avec les limites lithologiques observées dans le forage. Dans celui-ci des argiles limoneuses associées à une vitesse de propagation des ondes S de l'ordre de 175 m/s surmontent vers 9 m de profondeur des dépôts de moraine argilo-sableuse caractérisés par des vitesses de propagation des ondes S de l'ordre de 300 m/s jusqu'à 14 m de profondeur et supérieur ou égal à 400 m/s entre 14 et 20 m de profondeur. Le glissement de la Grande Combe (Ballaigues, VD) se produit à l'intérieur du remplissage quaternaire d'une combe creusée dans des calcaires Portlandien. Comme dans le cas du site d'Arnex les dépôts quaternaires correspondent à des dépôts glacio-lacustres. Dans la partie supérieure la surface de glissement a été localisée à une vingtaine de mètres de profondeur au niveau de l'interface qui sépare des dépôts de moraine jurassienne et des dépôts glacio-lacustres. Au pied du glissement 14 courbes de dispersions ont été déterminées sur des réseaux de 10 m d'extension le long d'un profil de 144 m. Les courbes obtenues sont discontinues et définies pour un domaine de fréquence de 7 à 35 Hz. Grâce à l'utilisation de distances source-récepteur entre 8 et 72 m, 2 à 4 modes de propagation ont été identifiés pour chacune des courbes. Lors de l'inversion des courbes de dispersion la prise en compte des différents modes de propagation a permis d'étendre la profondeur d'investigation jusqu'à une vingtaine de mètres de profondeur. Le modèle en deux dimensions permet de distinguer 4 couches (Vs1 < 175 m/s, 175 m/s < Vs2 < 225 m/s, 225 m/s < Vs3 < 400 m/s et Vs4 >.400 m/s) qui présentent des variations d'épaisseur. Des profils de sismiques réflexion en ondes S acquis avec une source construite dans le cadre de ce travail, complètent et corroborent le modèle établi à partir de l'analyse de la dispersion des ondes de surface. Un réflecteur localisé entre 5 et 10 m de profondeur et associé à une vitesse de sommation de 180 m/s souligne notamment la géométrie de l'interface qui sépare la deuxième de la troisième couche du modèle établi à partir de l'analyse de la dispersion des ondes de surface. Abstract Landslides are one of the main natural hazards in mountainous regions. In Switzerland, landslides cause damages every year that impact infrastructures and have important financial costs. In depth understanding of sliding mechanisms may help limiting their impact. In particular, this can be achieved through a better knowledge of the internal structure of the landslide, the determination of its volume and its sliding surface or surfaces In a landslide, the disorganization and the presence of fractures in the displaced material generate a change of the physical parameters and in particular a decrease of the seismic velocities and of the material density. Therefoe, seismic methods are well adapted to the study of landslides. Among seismic methods, surface-wave dispersion analysis is a easy to implement. Through it, shearwave velocity variations with depth can be estimated without having to resort to an S-wave source and to horizontal geophones. Its 3-step implementation implies measurement of surface-wave dispersion with long arrays, determination of the dispersion curves and finally inversion of these curves. Velocity models obtained through this approach are only valid when the investigated medium does not include lateral variations. In practice, this assumption is seldom correct, in particular for landslides in which reshaped layers likely include strong lateral heterogeneities. To assess the possibility of determining dispersion curves from short array lengths we carried out tests measurements on a site (Arnex, VD) that includes a borehole. A 190 m long seismic profile was acquired in a valley carved into limestone and filled with 30 m of glacio-lacustrine sediments. The data acquired along this profile confirmed that the presence of lateral variations under the geophone array influences the dispersion-curve shape so much that it sometimes preventes the dispersion curves determination. Our approach to use the analysis of surface-wave dispersion on sites that include lateral variations consists in obtaining dispersion curves for a series of short length arrays; inverting each so obtained curve and interpolating the different obtained velocity model. The choice of the location as well as the geophone array length is important. It takes into account the location of the heterogeneities that are revealed by the seismic refraction interpretation of the data but also, the location of signal amplitude anomalies observed on maps that represent, for a given frequency, the measured amplitude in the shot position - receiver position domain. The procedure proposed by Lin and Lin (2007) turned out to be an efficient one to determine dispersion curves using short extension arrays. It consists in building a time-offset from an array of geophones with a wide offset range by gathering seismograms acquired with different source-to-receiver offsets. When assembling the different data, a phase correction is applied in order to reduce static phase error induced by lateral variation. To evaluate this correction, we suggest to calculate, for two successive shots, the cross power spectral density of common offset traces. On the Arnex site, 22 curves were determined with 10m in length geophone-arrays. We also took advantage of the borehole to acquire a S-wave vertical seismic profile. The S-wave velocity depth model derived from the vertical seismic profile interpretation is used as prior information in the inversion of the dispersion-curves. Finally a 2D velocity model was established from the analysis of the different dispersion curves. It reveals a 3-layer structure in good agreement with the observed lithologies in the borehole. In it a clay layer with a shear-wave of 175 m/s shear-wave velocity overlies a clayey-sandy till layer at 9 m depth that is characterized down to 14 m by a 300 m/s S-wave velocity; these deposits have a S-wave velocity of 400 m/s between depths of 14 to 20 m. The La Grand Combe landslide (Ballaigues, VD) occurs inside the Quaternary filling of a valley carved into Portlandien limestone. As at the Arnex site, the Quaternary deposits correspond to glaciolacustrine sediments. In the upper part of the landslide, the sliding surface is located at a depth of about 20 m that coincides with the discontinuity between Jurassian till and glacio-lacustrine deposits. At the toe of the landslide, we defined 14 dispersion curves along a 144 m long profile using 10 m long geophone arrays. The obtained curves are discontinuous and defined within a frequency range of 7 to 35 Hz. The use of a wide range of offsets (from 8 to 72 m) enabled us to determine 2 to 4 mode of propagation for each dispersion curve. Taking these higher modes into consideration for dispersion curve inversion allowed us to reach an investigation depth of about 20 m. A four layer 2D model was derived (Vs1< 175 m/s, 175 m/s <Vs2< 225 m/s, 225 m/s < Vs3 < 400 m/s, Vs4> 400 m/s) with variable layer thicknesses. S-wave seismic reflection profiles acquired with a source built as part of this work complete and the velocity model revealed by surface-wave analysis. In particular, reflector at a depth of 5 to 10 m associated with a 180 m/s stacking velocity image the geometry of the discontinuity between the second and third layer of the model derived from the surface-wave dispersion analysis.

TECTONICS OF THE HELVETIC NAPPE SYSTEM (W SWITZERLAND): CONTROL OF STRAIN LOCALIZATION AND BASEMENT-COVER DEFORMATION. Insights from models based on continuum mechanics

The Helvetic nappe system in Western Switzerland is a stack of fold nappes and thrust sheets em-placed at low grade metamorphism. Fold nappes and thrust sheets are also some of the most common features in orogens. Fold nappes are kilometer scaled recumbent folds which feature a weakly deformed normal limb and an intensely deformed overturned limb. Thrust sheets on the other hand are characterized by the absence of overturned limb and can be defined as almost rigid blocks of crust that are displaced sub-horizontally over up to several tens of kilometers. The Morcles and Doldenhom nappe are classic examples of fold nappes and constitute the so-called infra-Helvetic complex in Western and Central Switzerland, respectively. This complex is overridden by thrust sheets such as the Diablerets and Wildhörn nappes in Western Switzerland. One of the most famous example of thrust sheets worldwide is the Glariis thrust sheet in Central Switzerland which features over 35 kilometers of thrusting which are accommodated by a ~1 m thick shear zone. Since the works of the early Alpine geologist such as Heim and Lugeon, the knowledge of these nappes has been steadily refined and today the geometry and kinematics of the Helvetic nappe system is generally agreed upon. However, despite the extensive knowledge we have today of the kinematics of fold nappes and thrust sheets, the mechanical process leading to the emplacement of these nappe is still poorly understood. For a long time geologist were facing the so-called 'mechanical paradox' which arises from the fact that a block of rock several kilometers high and tens of kilometers long (i.e. nappe) would break internally rather than start moving on a low angle plane. Several solutions were proposed to solve this apparent paradox. Certainly the most successful is the theory of critical wedges (e.g. Chappie 1978; Dahlen, 1984). In this theory the orogen is considered as a whole and this change of scale allows thrust sheet like structures to form while being consistent with mechanics. However this theoiy is intricately linked to brittle rheology and fold nappes, which are inherently ductile structures, cannot be created in these models. When considering the problem of nappe emplacement from the perspective of ductile rheology the problem of strain localization arises. The aim of this thesis was to develop and apply models based on continuum mechanics and integrating heat transfer to understand the emplacement of nappes. Models were solved either analytically or numerically. In the first two papers of this thesis we derived a simple model which describes channel flow in a homogeneous material with temperature dependent viscosity. We applied this model to the Morcles fold nappe and to several kilometer-scale shear zones worldwide. In the last paper we zoomed out and studied the tectonics of (i) ductile and (ii) visco-elasto-plastic and temperature dependent wedges. In this last paper we focused on the relationship between basement and cover deformation. We demonstrated that during the compression of a ductile passive margin both fold nappes and thrust sheets can develop and that these apparently different structures constitute two end-members of a single structure (i.e. nappe). The transition from fold nappe to thrust sheet is to first order controlled by the deformation of the basement. -- Le système des nappes helvétiques en Suisse occidentale est un empilement de nappes de plis et de nappes de charriage qui se sont mis en place à faible grade métamorphique. Les nappes de plis et les nappes de charriage sont parmi les objets géologiques les plus communs dans les orogènes. Les nappes de plis sont des plis couchés d'échelle kilométrique caractérisés par un flanc normal faiblement défor-mé, au contraire de leur flanc inverse, intensément déformé. Les nappes de charriage, à l'inverse se caractérisent par l'absence d'un flanc inverse bien défini. Elles peuvent être définies comme des blocs de croûte terrestre qui se déplacent de manière presque rigide qui sont déplacés sub-horizontalement jusqu'à plusieurs dizaines de kilomètres. La nappe de Mordes et la nappe du Doldenhorn sont des exemples classiques de nappes de plis et constitue le complexe infra-helvétique en Suisse occidentale et centrale, respectivement. Ce complexe repose sous des nappes de charriages telles les nappes des Diablerets et du Widlhörn en Suisse occidentale. La nappe du Glariis en Suisse centrale se distingue par un déplacement de plus de 35 kilomètres qui s'est effectué à la faveur d'une zone de cisaillement basale épaisse de seulement 1 mètre. Aujourd'hui la géométrie et la cinématique des nappes alpines fait l'objet d'un consensus général. Malgré cela, les processus mécaniques par lesquels ces nappes se sont mises en place restent mal compris. Pendant toute la première moitié du vingtième siècle les géologues les géologues ont été confrontés au «paradoxe mécanique». Celui-ci survient du fait qu'un bloc de roche haut de plusieurs kilomètres et long de plusieurs dizaines de kilomètres (i.e., une nappe) se fracturera de l'intérieur plutôt que de se déplacer sur une surface frictionnelle. Plusieurs solutions ont été proposées pour contourner cet apparent paradoxe. La solution la plus populaire est la théorie des prismes d'accrétion critiques (par exemple Chappie, 1978 ; Dahlen, 1984). Dans le cadre de cette théorie l'orogène est considéré dans son ensemble et ce simple changement d'échelle solutionne le paradoxe mécanique (la fracturation interne de l'orogène correspond aux nappes). Cette théorie est étroitement lié à la rhéologie cassante et par conséquent des nappes de plis ne peuvent pas créer au sein d'un prisme critique. Le but de cette thèse était de développer et d'appliquer des modèles basés sur la théorie de la méca-nique des milieux continus et sur les transferts de chaleur pour comprendre l'emplacement des nappes. Ces modèles ont été solutionnés de manière analytique ou numérique. Dans les deux premiers articles présentés dans ce mémoire nous avons dérivé un modèle d'écoulement dans un chenal d'un matériel homogène dont la viscosité dépend de la température. Nous avons appliqué ce modèle à la nappe de Mordes et à plusieurs zone de cisaillement d'échelle kilométrique provenant de différents orogènes a travers le monde. Dans le dernier article nous avons considéré le problème à l'échelle de l'orogène et avons étudié la tectonique de prismes (i) ductiles, et (ii) visco-élasto-plastiques en considérant les transferts de chaleur. Nous avons démontré que durant la compression d'une marge passive ductile, a la fois des nappes de plis et des nappes de charriages peuvent se développer. Nous avons aussi démontré que nappes de plis et de charriages sont deux cas extrêmes d'une même structure (i.e. nappe) La transition entre le développement d'une nappe de pli ou d'une nappe de charriage est contrôlé au premier ordre par la déformation du socle. -- Le système des nappes helvétiques en Suisse occidentale est un emblement de nappes de plis et de nappes de chaînage qui se sont mis en place à faible grade métamoiphique. Les nappes de plis et les nappes de charriage sont parmi les objets géologiques les plus communs dans les orogènes. Les nappes de plis sont des plis couchés d'échelle kilométrique caractérisés par un flanc normal faiblement déformé, au contraire de leur flanc inverse, intensément déformé. Les nappes de charriage, à l'inverse se caractérisent par l'absence d'un flanc inverse bien défini. Elles peuvent être définies comme des blocs de croûte terrestre qui se déplacent de manière presque rigide qui sont déplacés sub-horizontalement jusqu'à plusieurs dizaines de kilomètres. La nappe de Morcles and la nappe du Doldenhorn sont des exemples classiques de nappes de plis et constitue le complexe infra-helvétique en Suisse occidentale et centrale, respectivement. Ce complexe repose sous des nappes de charriages telles les nappes des Diablerets et du Widlhörn en Suisse occidentale. La nappe du Glarüs en Suisse centrale est certainement l'exemple de nappe de charriage le plus célèbre au monde. Elle se distingue par un déplacement de plus de 35 kilomètres qui s'est effectué à la faveur d'une zone de cisaillement basale épaisse de seulement 1 mètre. La géométrie et la cinématique des nappes alpines fait l'objet d'un consensus général parmi les géologues. Au contraire les processus physiques par lesquels ces nappes sont mises en place reste mal compris. Les sédiments qui forment les nappes alpines se sont déposés à l'ère secondaire et à l'ère tertiaire sur le socle de la marge européenne qui a été étiré durant l'ouverture de l'océan Téthys. Lors de la fermeture de la Téthys, qui donnera naissance aux Alpes, le socle et les sédiments de la marge européenne ont été déformés pour former les nappes alpines. Le but de cette thèse était de développer et d'appliquer des modèles basés sur la théorie de la mécanique des milieux continus et sur les transferts de chaleur pour comprendre l'emplacement des nappes. Ces modèles ont été solutionnés de manière analytique ou numérique. Dans les deux premiers articles présentés dans ce mémoire nous nous sommes intéressés à la localisation de la déformation à l'échelle d'une nappe. Nous avons appliqué le modèle développé à la nappe de Morcles et à plusieurs zones de cisaillement provenant de différents orogènes à travers le monde. Dans le dernier article nous avons étudié la relation entre la déformation du socle et la défonnation des sédiments. Nous avons démontré que nappe de plis et nappes de charriages constituent les cas extrêmes d'un continuum. La transition entre nappe de pli et nappe de charriage est intrinsèquement lié à la déformation du socle sur lequel les sédiments reposent.

ParA2, a Vibrio cholerae chromosome partitioning protein, forms left-handed helical filaments on DNA.

Most bacterial chromosomes contain homologs of plasmid partitioning (par) loci. These loci encode ATPases called ParA that are thought to contribute to the mechanical force required for chromosome and plasmid segregation. In Vibrio cholerae, the chromosome II (chrII) par locus is essential for chrII segregation. Here, we found that purified ParA2 had ATPase activities comparable to other ParA homologs, but, unlike many other ParA homologs, did not form high molecular weight complexes in the presence of ATP alone. Instead, formation of high molecular weight ParA2 polymers required DNA. Electron microscopy and three-dimensional reconstruction revealed that ParA2 formed bipolar helical filaments on double-stranded DNA in a sequence-independent manner. These filaments had a distinct change in pitch when ParA2 was polymerized in the presence of ATP versus in the absence of a nucleotide cofactor. Fitting a crystal structure of a ParA protein into our filament reconstruction showed how a dimer of ParA2 binds the DNA. The filaments formed with ATP are left-handed, but surprisingly these filaments exert no topological changes on the right-handed B-DNA to which they are bound. The stoichiometry of binding is one dimer for every eight base pairs, and this determines the geometry of the ParA2 filaments with 4.4 dimers per 120 A pitch left-handed turn. Our findings will be critical for understanding how ParA proteins function in plasmid and chromosome segregation.

Large-scale evolutionary trends of Acrochordiceratidae Arthaber, 1911 (Ammonoidea, Middle Triassic) and Cope's rule

Directed evolution of life through millions of years, such as increasing adult body size, is one of the most intriguing patterns displayed by fossil lineages. Processes and causes of such evolutionary trends are still poorly understood. Ammonoids (externally shelled marine cephalopods) are well known to have experienced repetitive morphological evolutionary trends of their adult size, shell geometry and ornamentation. This study analyses the evolutionary trends of the family Acrochordiceratidae Arthaber, 1911 from the Early to Middle Triassic (251228 Ma). Exceptionally large and bed-rock-controlled collections of this ammonoid family were obtained from strata of Anisian age (Middle Triassic) in north-west Nevada and north-east British Columbia. They enable quantitative and statistical analyses of its morphological evolutionary trends. This study demonstrates that the monophyletic clade Acrochordiceratidae underwent the classical evolute to involute evolutionary trend (i.e. increasing coiling of the shell), an increase in its shell adult size (conch diameter) and an increase in the indentation of its shell suture shape. These evolutionary trends are statistically robust and seem more or less gradual. Furthermore, they are nonrandom with the sustained shift in the mean, the minimum and the maximum of studied shell characters. These results can be classically interpreted as being constrained by the persistence and common selection pressure on this mostly anagenetic lineage characterized by relatively moderate evolutionary rates. Increasing involution of ammonites is traditionally interpreted by increasing adaptation mostly in terms of improved hydrodynamics. However, this trend in ammonoid geometry can also be explained as a case of Copes rule (increasing adult body size) instead of functional explanation of coiling, because both shell diameter and shell involution are two possible paths for ammonoids to accommodate size increase.

On prediction of individual sequences

Sequential randomized prediction of an arbitrary binary sequence isinvestigated. No assumption is made on the mechanism of generating the bit sequence. The goal of the predictor is to minimize its relative loss, i.e., to make (almost) as few mistakes as the best ``expert'' in a fixed, possibly infinite, set of experts. We point out a surprising connection between this prediction problem and empirical process theory. First, in the special case of static (memoryless) experts, we completely characterize the minimax relative loss in terms of the maximum of an associated Rademacher process. Then we show general upper and lower bounds on the minimaxrelative loss in terms of the geometry of the class of experts. As main examples, we determine the exact order of magnitude of the minimax relative loss for the class of autoregressive linear predictors and for the class of Markov experts.

An attempt to refine rockfall hazard zoning based on the kinetic energy, frequency and fragmentation degree

Rockfall hazard zoning is usually achieved using a qualitative estimate of hazard, and not an absolute scale. In Switzerland, danger maps, which correspond to a hazard zoning depending on the intensity of the considered phenomenon (e.g. kinetic energy for rockfalls), are replacing hazard maps. Basically, the danger grows with the mean frequency and with the intensity of the rockfall. This principle based on intensity thresholds may also be applied to other intensity threshold values than those used in Switzerland for rockfall hazard zoning method, i.e. danger mapping. In this paper, we explore the effect of slope geometry and rockfall frequency on the rockfall hazard zoning. First, the transition from 2D zoning to 3D zoning based on rockfall trajectory simulation is examined; then, its dependency on slope geometry is emphasized. The spatial extent of hazard zones is examined, showing that limits may vary widely depending on the rockfall frequency. This approach is especially dedicated to highly populated regions, because the hazard zoning has to be very fine in order to delineate the greatest possible territory containing acceptable risks.