Mes travaux de recherche s'inscrivent dans le cadre de la physique théorique des liquides et recourent  aux simulations numériques  (Monte-Carlo et dynamique moléculaire classique) pour étudier les propriétés de différents systèmes  :
 

• Propriétés physico-chimiques, thermodynamiques, structurales et dynamiques des liquides silicatés et carbonatés d'intérêt géologique

• Comportements génériques et universalité en fragmentation : fragmentation des noyaux et agrégats atomiques

• Gélification et transition vitreuse dans les suspensions colloïdales

• Structure en terme d'amas (percolation) des fluides simples supercritiques

Propriétés physico-chimiques, thermodynamiques, structurales et dynamiques des liquides silicatés et carbonatés d'intérêt géologique

en collaboration avec Bertrand Guillot (DR), Nicolas Folliet (post-doc), Boris Mantisi (ATER), Elsa Desmaele (Doctorante) et Thomas Dufils (Doctorant).
 

Les roches silicatées qui constituent le manteau terrestre sont soumises à des conditions thermodynamiques extrêmes. Le gradient thermique entre la base du manteau et la  croûte crée des mouvements convectifs (convection solide) à l'origine de la tectonique des plaques et dont une des manifestations la plus spectaculaire est le volcanisme. En effet, la pression diminuant sous l'effet des mouvements ascendants, une fusion partielle du manteau peut avoir lieu localement produisant des liquides magmatiques silicatés. Des gaz dissous dans ces magmas (H2O, CO2, gaz rares...) sont alors transportés du manteau à la surface et participent ainsi aux cycles biogéochimiques (comme celui du carbone). Pour mieux appréhender la géodynamique de notre planète et son évolution, il est nécessaire de connaitre les propriétés physico-chimiques des liquides magmatiques aux conditions régnant dans le manteau (1400-3000K et 0-1350 kbar).  C'est dans ce but que nous développons des champs de force permettant de calculer par dynamique moléculaire les propriétés des magmas naturels. Nous avons ainsi étudié la solubilité du carbone et des gaz rares dans ces magmas en collaboration avec des collègues géochimistes.

Comportements génériques et universalité en fragmentation

en collaboration avec Heli Hietala, Serguei Mechkov et François Sausset, (étudiants en thèse ou en stage de Master M2 au LPTMC).

Des noyaux atomiques aux supernovae, en passant par les assiettes et les roches, tout système physique peut être fragmenté, à condition de lui apporter une énergie suffisante par rapport à son énergie de cohésion. La fragmentation joue un rôle fondamental dans des domaines aussi variés que la détermination du diagramme de phase de la matière nucléaire, les éruptions volcaniques, ou la formation des planètes. Par ailleurs, environ 30% de la consommation totale d'énergie dans le monde est utilisée lors de processus industriels de fragmentation. Malgré cette diversité d'échelles de longueurs et de matériaux, la fragmentation présente des comportements génériques, peut-être universels. Ainsi, la distribution en taille des fragments obéit fréquemment à une loi de puissance caractérisée par un exposant noté τ. Cette loi d'échelle, mise en évidence par des expériences en laboratoire ou des observations de phénomènes naturels, suggère l'existence de classes d'universalité en fragmentation. Mais à ce jour, les théories et modèles ne parviennent pas à rendre compte de la dispersion des valeurs mesurées de l'exposant τ (entre 1 et 2.5 environ). La question de l'universalité en fragmentation reste donc ouverte.

La dynamique moléculaire est un outil idéal pour étudier ce phénomène irréversible et fortement hors d'équilibre qui pourrait dépendre de nombreux paramètres comme la dimension de l'espace, la forme, la taille ou le matériau de l'objet.

Le système modèle que nous proposons est constitué de particules bidimensionnelles interagissant par un potentiel central de type Lennard-Jones, avec une attraction à courte portée pour décrire les interactions à une échelle mésoscopique. Dans un premier temps, nous avons choisi d'étudier la fragmentation d'un disque envoyé contre un mur avec une certaine vitesse d'impact, le paramètre de contrôle de cette étude. Les calculs de dynamique moléculaire permettent alors de suivre l'évolution du système et de calculer la distribution en taille des éventuels fragments. Notons que nos résultats sont peu sensibles au choix du potentiel d'interaction, ce qui indique une certaine universalité, confirmée par des expériences.

A faible vitesse d'impact, le disque rebondit sur le mur après avoir subit des déformations élastiques, voire plastiques. Lorsque la vitesse d'impact augmente, le disque se fragmente.

Gélification et transition vitreuse dans les suspensions colloïdales

en collaboration avec Antonio Coniglio, Antonio de Candia, Emanuela del Gado, Analisa Fierro et Marco Tarzia de l'université Federico II de Naples, Italie.

Une suspension colloïdale est un ensemble de particules mésoscopiques dispersées dans un fluide porteur. Ces particules peuvent être solides ou liquides (on parle alors d'émulsions) et leur taille est typiquement comprise entre quelques nanomètres et quelques micromètres. Cette définition très générale regroupe une large variété de systèmes physico-chimiques présents dans la nature (comme le sang ou le lait) et dans le monde industriel (par exemple la peinture et les produits cosmétiques et agro-alimentaires).

En modifiant les propriétés physico-chimiques d'une suspension colloïdale, les expérimentateurs peuvent contrôler les interactions entre les particules et en particulier choisir la portée et l'intensité de l'interaction attractive, qui est typique des fluides simples atomiques. Ainsi, depuis une dizaine d'années, des expériences sur les systèmes colloïdaux ont mis en évidence des comportements thermodynamiques (transitions de phase) et dynamiques (transition vitreuse) rappelant ceux des systèmes atomiques. Mais à la différence de ces derniers, la grande taille des particules - comparée à celle des atomes du fluide porteur - non seulement facilite une observation directe et l'étude de leurs propriétés dynamiques, mais donne également lieu à des phénomènes mal compris comme la gélification colloïdale : à basse fraction volumique, et lorsque l'attraction entre les particules est élevée (ou à basse température), les particules ont tendance à former des amas dispersés dans une phase désordonnée (« cluster phase ») observée par microscopie confocale (voir figure 1). Puis, en augmentant la fraction volumique, le système présente les propriétés viscoélastiques et la dynamique lente d'un gel.

Pour mieux comprendre les mécanismes sous-jacents à la gélification colloïdale, nous avons proposé un modèle simple, basé sur un potentiel d'interaction de type DLVO présentant une attraction à très courte portée et une barrière répulsive à longue portée, due aux charges résiduelles observées à la surface des particules. L'étude de ce modèle par des calculs de dynamique moléculaire a porté sur :

• les propriétés thermodynamiques : détermination du diagramme de phase, révélant la présence d'une phase désordonnée et de différentes phases modulées : une phase colonnaire constituée de colonnes parallèles dont les axes sont placés sur un réseau triangulaire et une phase lamellaire (voir figure 2).

• les propriétés structurales : calculs du facteur de structure, de la distribution en taille des amas et de leurs rayons de giration (estimation de la dimension fractale).

• les propriétés dynamiques : calculs des fonctions de corrélation à deux temps et d'autocorrélation des liens entre particules permettant d'estimer le temps de relaxation structurel et la durée de vie des liens.

Nos résultats suggèrent le scénario suivant de la gélification colloïdale, résultant de la compétition entre attraction à courte portée et répulsion à longue portée : en augmentant la fraction volumique, des amas stables et compacts (« cluster phase », phase désordonnée) s'agrègent malgré la répulsion et forment des structures allongées qui à leur tour s'enchevêtrent pour former un réseau percolant (le gel) décrit par les exposants critiques de la percolation aléatoire. La gélification se produit lorsque le temps de relaxation structurel, qui augmente avec la fraction volumique, devient du même ordre de grandeur que la durée de vie des liens entre particules. L'étude du diagramme de phase montre que le gel est métastable par rapport aux phases modulées (colonnaire et lamellaire). Dans une suspension colloïdale, les inhomogénéités locales, la polydispersité des particules ou la viscosité du fluide porteur empêcheraient donc le système d'atteindre l'équilibre thermodynamique et favoriseraient la formation du gel. Ce scénario a été confirmé par les récentes expériences d' Andrew I. Campbell et de ses collaborateurs à l'Université de Bristol.

Structure en terme d'amas (percolation) des fluides simples supercritiques

Le travail réalisé au cours de ma thèse était d'établir la correspondance entre le diagramme de phase d'un fluide simple et sa structure microscopique en termes d'amas. Les calculs de dynamiques moléculaire et Monte-Carlo que nous avons réalisés sur des systèmes de grandes tailles suggèrent que des amas énergétiquement stables existent dans la phase supercritique, c'est à dire pour des températures supérieures à la température critique. Ces amas qui sont plus froids que le fluide, engendrent une ligne de percolation, caractérisée par les exposants de la percolation aléatoire à trois dimensions. Cette ligne, qui part du point critique thermodynamique, est pratiquement iso-énergétique. Par ailleurs, un amas macroscopique non percolant se forme au passage de la courbe de coexistence liquide-gaz pour des densités inférieures à la densité critique. Cette interprétation géométrique de la condensation renforce l'idée que ces amas sont des objets physiques. Des expériences réalisables dans des systèmes macroscopiques sont proposées pour mettre en évidence la présence de ces amas stables et en particulier l'existence de cette ligne de percolation. Le cas des systèmes de petites tailles est également considéré afin de discuter des expériences de fragmentation de noyaux et agrégats atomiques. En effet, par la simulation de l'expansion d'une gouttelette de fluide, nous avons montré que les amas stables étaient les précurseurs des fragments observés lors d'un processus de fragmentation. Bien que les les comportements critiques soient atténués, les propriétés des amas stables restent valables dans ces petits systèmes. Enfin, l'analyse d'une expérience de fragmentation de noyaux atomiques montre un très bon accord entre les données expérimentales et les prédictions théoriques. Outre son intéret sur le plan fondamental, ce travail pourrait permettre une meilleure compréhension des propriétés de solubilité des fluides supercritiques, utilisés dans de nombreuses applications industrielles.

Physique statistique de N. Sator et N. Pavloff (Vuibert, 2016)

                              Table des matières

Revues internationales avec comité de lecture   

L3 :

Introduction aux probabilités et aux statistiques

Introduction au principe variationnel et à la mécanique analytique

Introduction à la physique statistique

M1:

Mécanique statistique (parcours PAD, physique à distance) :  Vidéos de Physique statistique

M2 :

Du mouvement brownien à la modélisation financière

Travaux dirigés en simulation numérique en physique statistique

Travaux dirigés en physique statistique des fluides classiques

Introduction aux probabilités et aux statistiques   (L3)                  Licence de sciences et technologies        Mention Physique    Parcours L3-Phytem Ce cours présente des notions de base sur le calcul de probabilités et les variables aléatoires, quelques distributions de probabilité (gaussienne, binomiale, Poisson) et le théorème de la limite centrale.                                                             Cours  (30 pages en format pdf)                           Plan du cours                  TD (4 pages en format pdf)
	Introduction au principe variationnel et à la mécanique analytique(L3)
	Licence de sciences et technologies        Mention Physique    Parcours L3-Phytem Ces notes de cours présentent une introduction au principe variationnel et à  la mécanique analytique illustrée par quelques problèmes.                     Cours (36 pages en format pdf)            TD (2 pages en format pdf)
	Introduction à la physique statistique (L3)   Licence de sciences et technologies        Mention Physique    Parcours L3-Phytem L'objectif de ce cours est de présenter aux étudiants de L3 les bases de la mécanique statistique classique illustrée par des applications variées en thermodynamique et en physique.           Plan du cours et bibliographie           TD1   Grands nombres et espace des phases                                          TD2   Microétats quantiques et classiques           TD3   L'ensemble microcanonique           TD4   L'ensemble canonique I           TD5   L'ensemble canonique II           TD6 (DM)  L'ensemble grand-canonique                TD7   Particules en interaction I                    TD8   Particules en interaction II         (Energie moyenne d'Onsager)                       Un site web présentant des simulations en physique statistique    et un applet de percolation de sites sur réseau carré.       Quelques articles importants :              Modern Ergodic Theory, Joël Lebowitz and Oliver Penrose,            Physics Today, 26 (February 1973), 155-175 pdf            Boltzmann's entropy and time's arrow, Joël Lebowitz,            Physics Today, 46 (September 1993), 32-38 pdf            Boltzmann's Approach to Statistical Mechanics, Sheldon Goldstein,            Lecture Notes in Physics, 574 (2001), 39-54 pdf            Science of Chaos or Chaos in Science, Jean Bricmont,            Physicalia Magazine, 3-4 (1995), 159-208 pdf            Ludwig Boltzmann, les théories physiques et les atomes, Anouk Barberousse,            Images de la physique, (2007)  pdf
	Du mouvement brownien à la modélisation financière (M2)
	Master de sciences et technologies Spécialité "Systèmes dynamique et statistiques de la matière complexe" Parcours "Modélisation, statistique et algorithmique"   Depuis plus d'un siècle, le cours des actifs financiers est étudié sous l'angle des marches aléatoires. Le mouvement brownien et le phénomène de diffusion sont donc intimement liés à la modélisation financière. Cette UE propose une initiation à la finance de marché basée sur des concepts physiques et mathématiques utilisés pour modéliser l'évolution des actifs financiers. Dans ce cours, nous nous intéresserons au mouvement brownien, à travers l'équation de Langevin, et au phénomène de diffusion décrit en terme d'une marche aléatoire à la limite continue. Nous aborderons ensuite le processus de Wiener et le lemme d'Îto,  qui nous permettront d'introduire le modèle de Black-Scholes, pierre angulaire de la modélisation financière.     Cours : Plan du cours et bibliographie     Travaux dirigés : TD1     Le mouvement brownien : Einstein et Langevin TD2     Marches aléatoires discrètes TD3     Equations de la diffusion et de Fokker-Planck TD4     Modèle de Black-Scholes   TD5     L'approche risque-neutre et le modèle de Black       Travaux pratiques sur ordinateur Les séances de TP auront lieu de 14h à 18h les mardi 15 novembre 2016 22 novembre 2016 29 novembre 2016   6 décembre 2016 dans la salle informatique 111 du couloir 22-23  au 1er étage   TP1     Marches aléatoires TP2     Cours d'un actif financier TP3     Le modèle de Black-Scholes Data 1 Data 2 Data 3 Data 4
	Travaux dirigés en simulation numérique en physique statistique (M2) Master 2 de Sciences et Technologies      Mention Physique  et Applications Spécialité "Systèmes dynamiques et statistiques de la matière complexe" Parcours "Modélisation, statistique et algorithmique" Cours de Pascal Viot Travaux dirigés (en format pdf)           TD1      Systèmes modèles : Le fluide simple et le modèle d'Ising           TD2      Simulations Monte-Carlo 1 : Probabilité d'acceptation           TD3      Simulations Monte-Carlo 2 : Algorithme de Wang-Landau et méthode d'insetion de Widom           TD4      Dynamique  Moléculaire à température constante           TD5      Analyse en taille finie pour les transitions de phases continues           TD6      Système hors équilibre : Modèle de diffusion-coagulation           TD7      Système hors équilibre : Modèle ASEP                        Travaux dirigés en physique statistique des fluides classiques (M2) Master 2 de Sciences et Technologies      Mention Physique  et Applications Spécialité "Concepts fondamentaux de la physique" Parcours "Physique des liquides et matière molle" Cours de Claude Aslangul Travaux dirigés (en format pdf)           TD1      Fluctuations           TD2      Fluides de sphères dures et particules colloïdales           TD3      Développement du viriel           TD4      Transition de phase liquide-gaz et amas de particules           TD5      Marche aléatoire sur réseau