L'objet de cette these est l'etude des transitions de phase dans les modeles de spins xy desordonnes. Les outils numeriques que nous avons utilises sont la methode monte carlo et la methode de groupe de renormalisation de paroi de domaine. Les resultats que nous avons obtenus sont le fruit de simulations realisees a partir de codes originaux mis au point sur machines paralleles. Le premier type de modele que nous avons etudie est un modele a desordre de liaisons : le modele bimodal ou regne une competition des interactions de signe aleatoire gele entre premiers voisins. Une etude monte carlo suivie d'une analyse d'echelle de taille finie nous a permis de determiner le diagramme de phase du modele en dimension 2 ainsi que l'ensemble de ses exposants critiques. Dans la region de desordre fort, les transitions de verre de spin et de verre chiral qui se produisent a temperature nulle, sont decouplees. Dans le cas du modele tridimensionnel a desordre maximum, nous avons grace a l'apport d'un algorithme tres efficace de recherche de l'etat fondamental, obtenu des resultats d'energies de defaut de paroi de domaine qui montrent que la transition de verre chiral se produit a temperature finie et que la transition de verre de spin se produit probablement a temperature tres faible. Nos resultats montrent que la dimension critique inferieure du modele est proche de 3. Le deuxieme type de modele aborde est un modele ou le desordre est introduit par l'existence de dephasages aleatoires geles entre les spins. Par une utilisation poussee de l'hypothese d'universalite, nous avons confirme une recente theorie qui exclu l'existence d'une phase reentrante dans le cas bidimensionnel. Ainsi, la phase kosterlitz-thouless est stable jusqu'a temperature nulle en dessous du desordre critique. En dimension 3, nous avons montre par une etude de groupe de renormalisation de paroi de domaine que la transition de verre de spin a lieu a temperature finie.
L'objet de cette these est l'etude des transitions de phase dans les modeles de spins xy desordonnes. Les outils numeriques que nous avons utilises sont la methode monte carlo et la methode de groupe de renormalisation de paroi de domaine. Les resultats que nous avons obtenus sont le fruit de simulations realisees a partir de codes originaux mis au point sur machines paralleles. Le premier type de modele que nous avons etudie est un modele a desordre de liaisons : le modele bimodal ou regne une competition des interactions de signe aleatoire gele entre premiers voisins. Une etude monte carlo suivie d'une analyse d'echelle de taille finie nous a permis de determiner le diagramme de phase du modele en dimension 2 ainsi que l'ensemble de ses exposants critiques. Dans la region de desordre fort, les transitions de verre de spin et de verre chiral qui se produisent a temperature nulle, sont decouplees. Dans le cas du modele tridimensionnel a desordre maximum, nous avons grace a l'apport d'un algorithme tres efficace de recherche de l'etat fondamental, obtenu des resultats d'energies de defaut de paroi de domaine qui montrent que la transition de verre chiral se produit a temperature finie et que la transition de verre de spin se produit probablement a temperature tres faible. Nos resultats montrent que la dimension critique inferieure du modele est proche de 3. Le deuxieme type de modele aborde est un modele ou le desordre est introduit par l'existence de dephasages aleatoires geles entre les spins. Par une utilisation poussee de l'hypothese d'universalite, nous avons confirme une recente theorie qui exclu l'existence d'une phase reentrante dans le cas bidimensionnel. Ainsi, la phase kosterlitz-thouless est stable jusqu'a temperature nulle en dessous du desordre critique. En dimension 3, nous avons montre par une etude de groupe de renormalisation de paroi de domaine que la transition de verre de spin a lieu a temperature finie.
Parmi les méthodes de la théorie du groupe de renormalisation, les développements en epsilon sont basés sur des calculs dans l'espace réciproque et permettent de calculer les exposants critiques associés aux transitions magnétiques du second ordre pour des valeurs non entières de la dimension d'espace. Une interprétation physique naturelle consiste à se demander comment se comporte un système de spins en interaction dans un espace de dimension fractale. Or les structures fractales sont construites par itération d'une cellule génératrice dont les détails sont donc présents à plusieurs échelles ; la question qui se pose alors est de savoir ce qui se passe lorsque l'invariance par translation est remplacée par une invariance d'échelle géométrique. La convergence vers la limite thermodynamique se produit en même temps que le processus d'itération construit la structure. De ce fait, des simulations Monte Carlo fiables de ces transitions de phase n'ont pu être menées à bien que récemment, puisqu'elles nécessitent la simulation de très grandes tailles, lesquelles varient comme des séries géométriques avec l'étape d'itération. C'est en utilisant des algorithmes non locaux dits “d'amas” (Wolff, Swendsen-Wang), capables de réduire le ralentissement critique de manière significative, et des méthodes d'histogrammes pour traiter les données des simulations Monte-Carlo que j'ai tout d'abord réalisé ces études. Il s'avère que le calcul précis des exposants critiques est rendu encore plus difficile par le fait que l'analyse en tailles finies du modèle d'Ising souffre de corrections d'échelle qui peuvent affecter fortement le comportement de certaines grandeurs thermodynamiques, en particulier lorsque la dimension fractale tend vers 1. J'ai montré que ces corrections d'échelle sont en partie liées à la très forte inhomogénéité du réseau sous-jacent (due à l'existence de trous sur plusieurs ordres de grandeurs) et à la concomitance de la construction du fractal avec la convergence vers la limite thermodynamique. Les résultats que j'ai obtenus pour les exposants critiques, ou leurs bornes, sont toujours compatibles avec la relation d'hyperscaling dans laquelle on substitue la dimension de Hausdorff à la dimension d'espace. Le comportement critique en dimension non entière se décrit dans le cadre de l'universalité faible. Cela se manifeste par un désaccord net entre les exposants que j'ai obtenus par les méthodes Monte Carlo et les développements en epsilon. Les exposants critiques ne dépendent pas seulement de la dimension d'espace, des propriétés de symétrie du paramètre d'ordre et de la portée des interactions, mais aussi des propriétés géométriques de la structure fractale : Très récemment des calculs précis d'exposants critiques m'ont permis de montrer que des classes d'universalité différentes sont en général nécessaires pour décrire le comportement du modèle d'Ising sur des fractals de même dimension et de lacunarités différentes. Un tel résultat généralise le concept d'universalité faible proposé par Masuo Suzuki. L'hypothèse d'homogénéité qui sous-tend les lois d'échelle permettant de décrire un comportement critique se dérive par renormalisation. La procédure de renormalisation dans l'espace direct est naturelle dans les fractals, puisqu'elle suit exactement le processus inverse de construction de la structure. Avec mon étudiant Pai-Yi Hsiao, nous avons mené à bien l'étude du modèle d'Ising par une méthode de renormalisation Monte-Carlo sur une structure fractale de dimension voisine de 1, 89 ; il s'avère que l'exposant associé à l'une des directions propres peut être calculé avec une très bonne précision et est en accord avec les résultats de l'analyse en tailles finies. En revanche, la convergence est très lente dans l'autre direction, ce qui est lié aux corrections d'échelle mises en évidence lors de cette analyse. La cinétique stochastique associée à la formation des amas construits par l'algorithme de Wolff sous tend la compréhension du phénomène de ralentissement critique. J'ai montré que les distributions des tailles des amas de Wolff ont une propriété d'homogénéité qui fait intervenir l'un des exposants associé à une des directions propres du processus de renormalisation. Par ailleurs, les distributions des tensions de surface des amas vérifient une propriété analogue dans laquelle intervient un nouvel exposant critique. L'étude des fonctions d'autocorrélation m'a permis de calculer précisément les exposants dynamiques de Wolff lorsque la température critique est connue, et d'éclaircir l'évolution du ralentissement critique avec la dimension et la connectivité. Dans le cas de systèmes invariants par translation, l'ordre de la transition ferromagnétique du modèle de Potts est lié au nombre d'états de spin ; le passage du premier au second ordre est attendu pour des dimensions non entières. Par ailleurs, la présence de désordre peut, dans certaines conditions, induire une transition du second ordre à partir d'un système qui en présente une du premier. L'étude du comportement critique du modèle de Potts sur des fractals est donc particulièrement intéressante, puisque deux des paramètres qui le déterminent (dimensionnalité et désordre structurel) sont liés. Avec mon étudiant Pai-Yi Hsiao, nous avons montré que la transition associée au modèle de Potts à trois états sur une structure fractale de dimension voisine de 1, 89 est du second ordre. Les difficultés attendues lorsqu'on augmente le nombre d'états de spins se font déjà nettement sentir : Les corrections d'échelle empêchent de calculer la température critique avec une très bonne précision. Nous n'avons donc pu donner que des bornes pour certains exposants ; nous avons cependant clairement mis en évidence la différence entre les classes d'universalité de Potts à 2 et 3 états. L'étude de la percolation en dimension non entière est liée à la fois à celle du modèle de Potts et aux algorithmes d'amas. Elle est basée sur l'étude des moments de la distribution de taille des amas, ce qui nécessite la localisation de pics en fonction de la probabilité d'occupation. J'ai pu montrer que les corrections d'échelle n'affectent pratiquement pas le comportement des pics avec la taille des structures, et proposé de les interpréter en termes de "seuil effectif".
Le sujet de cette these porte sur les verres de spin ising, qui sont les archetypes des systemes complexes en physique des solides. Bien que la phase verre de spin presente des aspects remarquables tels que le vieillissement nous nous sommes plutot interesse au comportement des verres de spin a la transition. Dans la premiere partie de cette these les verres de spin sont introduit d'un point de vue experimental, afin de degager les concepts clefs, a savoir la frustration et le desordre. Muni de ces concepts il est alors possible d'aborder les modeles theoriques. Enfin pour conclure la partie introductive une presentation des bases de la methode monte carlo est donnee. La seconde partie est la presentation de nos resultats. Dans un premier temps nous avons etudie l'influence de la distribution des interactions sur les exposants critiques. Dans le cadre de la theorie des phenomenes critiques les exposants critiques ne doivent varier qu'avec la dimension d'espace et le nombre de composantes du parametre d'ordre ; c'est la notion de classe d'universalite. Nous apportons ici une preuve numerique de la non-universalite pour les verres de spin ising. Dans un second temps nous nous sommes interesse aux proprietes des verres de spin pour des temperatures superieures a la temperature de transition. En effet l'etude de la relaxation dans les verres de spin indique l'existence d'un regime non-exponentiel compris dans une gamme de temperatures allant de la temperature de transition a une temperature environ quatre fois superieure. Afin de determiner cette temperature nous avons utilise plusieurs techniques numeriques. Enfin nous avons etudie l'influence du desordre sur un systeme completement frustre. Ces deux dernieres etudes nous ont alors permis de donner une explication simple, basee sur l'espace des phases, de la transition verre de spin
Dans cette these, nous etudions des modeles de spin avec des interactions a longue portee de type ferromagnetique, decroissantes avec la distance entre les particules, de la forme 1/r d + , en utilisant les techniques de simulation de monte carlo en general, combinees aux analyses des lois d'echelle en taille finie. Deux regimes sont consideres : le regime classique ou champ moyen (0 < < d/2) et le regime non classique (d/2 < < d). Les differentes predictions des analyses du groupe de renormalisation sont exposees dans le premier chapitre pour ce type d'interactions. Le chapitre i est consacre aux rappels sur les phenomenes critiques et le chapitre ii a la presentation des differentes techniques utilisees dans cette these. A savoir la methode de simulation monte carlo, la technique des histogrammes de monte carlo et la methode des fonctions de green appliquee au ferromagnetisme. Dans le chapitre iii, nous etudions le modele d'ising avec spin infini a une et deux dimensions. On trouve une transition de phase de type continu ; les exposants critiques sont estimes pour toutes les valeurs etudiees de. En general, nos resultats sont en accord avec la tendance predite par la theorie du groupe de renormalisation dans la region dite non classique. Ce n'est pas le cas dans la region champ moyen. Le modele de potts est traite dans le chapitre iv, a une, deux et trois dimensions. Non seulement le systeme unidimensionel presente une transition de phase mais, pour une valeur de donnee, la nature de cette transition depend du nombre d'etats q, comme dans le cas des interactions a courte portee en dimension deux et trois. Les resultats obtenus a 2d et 3d suggerent egalement l'existence d'une frontiere entre une zone de transition du premier ordre et une zone de transition critique dans le plan des parametres (,q) qualitativement semblable au cas 1d. Finalement, nous considerons les modeles de spins vectoriels, les modeles xy et d'heisenberg unidimensionnels dans le chapitre v. Les systemes presentent un ordre a longue portee a temperature finie pour inferieur a 1. Pour le cas particulier = 1, les resultats de la simulation numerique tendent a confirmer l'existence d'une transition de phase avec une aimantation nulle. Quant a la theorie des fonctions de green, il y aurait une transition de phase dans le cas du modele heisenberg a temperature finie mais pas pour le modele x-y. Certaines predictions du groupe de renormalisation ne sont pas verifiees. Par exemple, nous trouvons que la temperature critique n'a pas un comportement lineaire en mais plutot en loi de puissance.
Nous avons aborde au cours de ce travail de these differents aspects du desordre et de la frustration dans les systemes de spins. La premiere partie est dediee aux effets de la frustration sur le comportement critique du modele de potts bidimensionnel a q = 3 etats sur le reseau dit de villain. Utilisant des variables grassmaniennes et une analogie avec le modele des dimeres, nous analysons l'etat fondamental du systeme en fonction de x = j a n t i f e r r o/j f e r r o. Nous trouvons que malgre la degenerescence infinie du fondamental et le fait que l'entropie residuelle soit finie, l'ordre a longue portee emerge dans le systeme completement frustre (x = 1) a temperature critique finie, contrairement au cas q = 2. La technique des matrices de transfert est utilisee afin d'acceder au diagramme de phase du systeme. L'evolution des exposants critiques en fonction de x dans l'intervalle 1, 1 conduit a la conclusion que le systeme evolue dans la meme classe d'universalite que son homologue non-frustre. Le diagramme de phase est verifie utilisant la simulation monte-carlo. Les exposants critiques au point x = 1 sont obtenus utilisant la methode et l'analyse en taille finie. Le deuxieme aspect aborde dans cette these est l'influence du faible desordre sur le point tricritique du modele blume-capel (nomme bc par la suite) en trois dimensions. Par l'intermediaire des simulations monte-carlo intensives, nous suivons la trajectoire du point tricritique du modele de bc en fonction du desordre gele de liens. Nos simulations verifient que la transition de premier ordre est presente dans ce systeme jusqu'a un seuil de desordre des liens. Les techniques des histogrammes et l'analyse en taille finie donnent acces aux exposants critiques du modele desordonne. Nous verifions ainsi que la classe d'universalite change suite a la presence du desordre gele dans le systeme. Les memes conclusions sont retrouvees dans le cas d'un potentiel d'anisotropie variant selon une distribution bimodale.
D ans cette présente thèse de doctorat et dans le but de découvrir le comportement magnétique des matériaux en fonction de di érents paramètres extérieurs, nous avons proposé quelques modèles de spin sur réseaux pour l’étude. Pour ce fait, nous avons donné en détails, dans le premier chapitre, l’origine du magnétisme dans les di érents matériaux et leur comportement en fonction du changement des paramètres. Dans le second chapitre, nous avons mis l’accent sur un phénomène extraordinaire qui se manifeste au niveau des matériaux qui est la transition de phase. Cette dernière n’est qu’un changement d’état subi à cause des modifications au niveau des paramètres. Nous avons cité les di érentes classifications des transitions de phase notamment la théorie d’Ehrenfest et celle de Landau. Nous avons donné un peu d’intérêt à la deuxième classification avec plus de détails. Nous avons présenté aussi les phénomènes tricritiques qui mènent au calcul des di érents exposants critiques. Nous avons expliqué aussi la notion des classes d’universalité auxquelles appartiennent les di érents systèmes étudiés. Dans le troisième chapitre, nous avons présenté notre travail de fond. Nous avons défini le modèle à étudier à savoir le modèle de Blume-Emery-Gri ths. nous avons donné un aperçu sur les méthodes qu’on a utilisé à citer la méthode de Monte-Carlo et la méthode (finite size scaling) TMFSS basée sur la matrice de transfert. Ces méthodes ont montré leurs preuves dans plusieurs études antérieures. Le diagramme de phase obtenu, riche avec ses caractéristiques, exhibe une ligne de transition de second ordre. Nous avons justifié la présence des di érentes phases par étude numérique de l’évolution thermique de l’aimantation du système, du moment quadrupolaire et du moment octupolaire. Les résultats sont en accord avec ceux des états fondamentaux. En utilisant la méthode de la matrice de transfert, nous avons étudié le comportement critique du système. Nous avons calculé l’exposant critique pour quelques valeurs des constantes d’interaction. Le quatrième chapitre a été consacré à l’étude des états fondamentaux de deux modèles importants. Le premier est le modèle de Blume-Emery-Gri ths à spin-1. Les interactions sont antiferromagnétiques avec le premier plus proche voisin. Nous avons considéré un système triangulaire avec sa particularité de frustration. Nous avons fait une étude analytique qui nous a donné les énergies des di érentes configurations. Une étude numérique nous a montré les phases les plus stables qu’on a représenté sur un diagramme de phase. Le deuxième modèle était le modele de Blume-Capel à spin-3/2 sur un réseau carré. Nous avons suivi les mêmes étapes d’étude. Comme perspective, il serait intéressant de compléter les deux derniers travaux par des études numériques. Nous prévoyons utiliser les deux méthodes exactes à savoir la méthode de Monte-Carlo et celle de la matrice de transfert. Nous essayerons de calculer les di érents exposants critiques et de déterminer leur classe d’universalité.
L'objet de cette these est l'etude des proprietes statiques et dynamiques du modele de spins xxo bimodal sur reseau carre. A partir de simulations numeriques intensives, nous obtenons son diagramme de phase (l,x) classique a temperature nulle ou x est la concentration des defauts magnetiques et l leur force. Nous mettons en evidence, par le calcul de differents parametres d'ordre et des fonctions de correlation spin-spin, une transition qui separe une phase de decroissance algebrique des correlations et un verre de spin isotrope caracterise par une decroissance exponentielle des correlations. Nos resultats suggerent egalement une aimantation nulle pour toute concentration x finie des defauts a la limite thermodynamique. Nous etudions ensuite la dynamique du systeme. Le calcul des quantites hydrodynamiques qui la caracterisent montre l'existence d'ondes de spin de dispersion lineaire dans la limite des basses frequences. L'obtention du facteur de structure dynamique du systeme dans l'approximation harmonique des ondes de spin permet d'acceder aux composantes dans le plan xy et transverse des fluctuations magnetiques. Nous en deduisons differentes quantites qui permettent de discuter du comportement dynamique du systeme. Nous procedons egalement a une modelisation de l'ensemble de ces proprietes dynamiques dans le cadre de l'approximation cpa (coherent potential approximation). Les resultats obtenus par cette approche, qui permet d'alleger notablement la charge numerique et d'acceder a de plus grandes tailles de systemes, sont en tres bon accord avec nos resultats des calculs exacts. Enfin nous etudions l'effet de la frustration sur les fluctuations quantiques du systeme en calculant la premiere correction en 1/s aux grandeurs hydrodynamiques. Nous montrons l'existence, dans le cas ou s=1/2, d'une phase quantique desordonnee a temperature nulle caracterisee par une constante de raideur nulle. Nous proposons egalement une approche originale pour investiguer les proprietes du systeme dans cette phase