Etude numérique des corrections d'échelle au comportement dominant à l'équilibre et hors de l'équilibre

Jean‐Charles Walter

Published Oct 15, 2009• 0 citations

Abstract

La première partie a pour sujet le comportement à l'équilibre du modèle d'Ising pour d>4. Dans un premier temps, nous étudions le comportement thermique dans le cadre du comportement d'échelle étendu. Par interpolation de données numériques en dimensions cinq à huit, nous obtenons un développement décrivant la susceptibilité dans toute la phase haute température. Dans un second temps, nous étudions les effets de taille finie. Les résultats numériques obtenus pour le modèle d'Ising 5d sont compatibles avec une croissance anormale de la longueur de corrélation pour des conditions de bords libres. La seconde partie a pour sujet le vieillissement dans les systèmes de spins 2d complètement frustrés. Dans un premier temps, nous étudions le vieillissement du modèle d'Ising complètement frustré 2d lors d'une trempe depuis la phase haute température jusqu'à la température critique. La présence de défauts topologiques, comme pour le modèle XY, se manifeste par des corrections logarithmiques lors de la croissance de la longueur caractéristique. Dans un second temps, nous étudions le vieillissement du modèle XY complètement frustré 2d. Lors d'une trempe depuis l'état fondamental jusque dans la ligne critique, le vieillissement des spins est bien décrit par les ondes de spins. Lors d'une trempe depuis la phase haute température jusqu'à la température BKT des spins et jusqu'à la température de brisure de symétrie de la chiralité, nous estimons les grandeurs universelles des deux variables. Les résultats pour la chiralité sont incompatibles avec la classe d'universalité du modèle d'Ising 2d. Des corrections logarithmiques sont également présentes.

DOI: 10.70675/8d6559bfz3321z467fzb36az47db36f23b23

Related papers

Comportement critique de modèles bidimensionnels en présence de perturbations inhomogènes

Farkas Adam Bagaméry

Cette thèse est consacrée à l'étude de l'influence de différentes perturbations sur le comportement critique de systèmes bidimensionnels. Dans une première partie nous considérons le modèle d'Ising 2D placé dans un champ de surface aléatoire qui constitue une perturbation marginalement non-pertinente au point critique. Nous étudions le problème à l'aide de simulations Monte Carlo pour différentes amplitudes de désordre. Nous déterminouns des exposants critiques effectifs dont la dépendance en température ou en taille est en accord avec les résultats de théorie des champs et interprétable en termes de corrections logarithmiques au comportement du système pur. Dans la seconde partie, nous étudions le modèle d'Ising sur le réseau carré avec un désordre corrélé biaxial qui constitue une perturbation pertinente d'après le critère de Harris étendu. Le problème est abordé par simulations Monte Carlo extensives. L'exposant n prend la valeur attendue dans un système avec désordre corrélé isotrope tandis que la densité d'aimantation a une dimension supérieure à celle du système pur. Les propriétés conformes des densités d'énergie et d'aimantation sont aussi examinées. Finalement, nous considérons le comportement critique à l'interface entre deux sous-systèmes appartenant à des classes d'universalité différentes mais ayant la même température critique. Le problème est étudié dans le cadre d'une théorie de champ moyen et en utilisant des arguments d'échelle. La grande variété de comportements d'interface obtenue est confrontée à une étude numériques sur des interfaces de Potts 2D à q états (2< q<4).

Jan 1, 20060 citations

Critical behaviour of two-dimensional models under the influence of inhomogeneous perturbations

Farkas Adam Bagaméry, Loïc Turban, Ferenc Iglói

Jan 1, 20060 citations

Out-of-equilibrium dynamics in classical field theories and Ising spin models

Hugo Ricateau

Dynamique hors équilibre des théories classiques des champs et des modèles de spin d’Ising Cette thèse est constituée de deux parties indépendantes. Dans le premier chapitre, nous introduisons une méthode numérique permettant d'intégrer des équations aux dérivées partielles représentant la dynamique Hamiltonienne de théories des champs. Cette méthode est un intégrateur multi-symplectique qui préserve localement le tenseur énergie-impulsion sur de très longues périodes de temps et avec précision. Son principal avantage est d'être extrêmement simple tout en restant bien définie localement. Nous la mettons à l'épreuve sur le cas particulier du modèle phi^4 en 1+1 dimensions; nous expliquons également comment l'implémenter en dimensions supérieures. De plus, nous faisons une présentation géométrique de la structure multi-symplectique et nous introduisons une construction permettant de résoudre le problème de dégénérescence pouvant l'affecter.Le second chapitre traite d'aspects hors équilibre dans les systèmes statistiques: nous nous intéressons en particulier à la question de l'impact d'un taux de refroidissement fini lors d'une trempe à travers une transition de phase du second ordre. Pour décrire plus fidèlement le régime hors équilibre qui se produit avant la transition de phase, nous étendons le mécanisme dit de Kibble-Zurek. Nous décrivons comment la taille caractéristique des objets géométriques présents dans le système dépend du temps et du taux de refroidissement; ceci, avant et une fois le point critique atteint. Ces prédictions théoriques sont mises à l'épreuve sur l'exemple du modèle d'Ising ferromagnétique. Nous décrivons également les propriétés géométriques des domaines qui apparaissent dans le système au cours de la dynamique de refroidissement.

Sep 29, 20170 citations

TRANSITIONS DE PHASE EN DIMENSIONS FRACTALES

Pascal Monceau

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".

Dec 16, 20040 citations

Dynamique et exposants critiques dans les verres de spin ising

L. Bernardi

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

Jan 1, 19960 citations

Etude par simulations numériques des transitions de phase dans les modèles de spins XY désordonnés

J. Maucourt

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.

Jan 1, 19970 citations