Depuis la nuit des temps, la matière est connue sous trois états, soient l'état solide, l'état liquide et l'état gazeux. Depuis, de nouveaux états ont été découverts. Les études au sein des laboratoires des professeurs Soldera et Lebel se concentrent sur ces états particuliers. Ils travaillent surtout sur les transitions de phases menant à ces états. Lorsqu'on parle de solide, on fait généralement référence à l'état cristallin, qui est un réseau bien structuré où les interactions à longue portée sont importantes et bien définies. Il existe toutefois une autre forme de solide où le désordre semble régner. Les interactions à longue portée ne sont pas prépondérantes. Il s'agit de l'état vitreux. Le passage de l'état liquide à l'état vitreux est appelé la transition vitreuse, et est caractérisé par une température, T[indice inférieur g] . Les principaux composés menant à l'état vitreux sont les polymères et les verres inorganique [i.e. inorganiques] tels que l'oxyde de silice, très répendus [i.e. répandus] et grandement étudiés. Les verres organiques sont une autre famille de composés qui présentent une phase amorphe stable. Le développement des technologies de pointe a suscité un intérêt pour les verres moléculaires dû à leurs propriétés intéressantes telles que leur structure bien définie et une bonne stabilité thermique. La T[indice inférieur g] détermine le domaine de températures où il peut être utilisé selon l'application recherchée. Les verres moléculaires ont beaucoup de potentiel en industrie, mais ils constituent une famille de matériaux fonctionnels méconnus. S'il reste beaucoup à apprendre sur les polymères, tout reste à apprendre sur les verres moléculaires. Il y a moins d'une vingtaine d'années que la triazine a commencé à être étudiée comme précurseur de verres moléculaires. La synthèse de verres efficaces peut être fastidieuse et ne procure pas nécessairement les propriétés recherchées. C'est pour cette raison que cette étude a été entreprise. L'idée de base de ce travail vise à trouver des moyens de prédire la T[indice inférieur g] , la première propriété recherchée, et ce, en se basant seulement sur la structure moléculaire. Pour relever ce défi, la voie considérée est la simulation atomistique faisant appel à un modèle tout-atome. Cette méthode a déjà fait ses preuves dans l'étude de la T[indice inférieur g] de polymères vinyliques, par les travaux du groupe du professeur Soldera. Bien que le principe de base soit resté le même, il [a] fallu apporter de légères modifications quant à son application aux systèmes étudiés. Suite à la paramétrisation du champ de forces utilisé, pcff, la simulation d'une quinzaine de molécules dont la T[indice inférieur g] expérimentale est connue a pu être conduite. Les résultats obtenus ont permis de révéler une linéarité entre les valeurs obtenues par simulation et les données expérimentales. Il est important de noter que la méthode a permis de distinguer les différents verres moléculaires étudiés bien que la seule différence entre les composés soit le groupement fonctionnel. Cet accord a permis d'effectuer diverses analyses au niveau moléculaire afin de mieux comprendre les facteurs qui influencent les valeurs de T[indice inférieur g] , et donc in fine d'apporter des indices pour une meilleure compréhension de cette transition fort complexe. Les fonctions de distribution radiale ainsi que la rotation des différentes fonctionnalités par rapport au noyau commun triazine, ont ainsi été regardées.

Dans les alliages 3d, les propriétés magnétiques des solutés peuvent être extrêmement sensibles aux environnements chimiques locaux, et avoir un impact crucial sur diverses propriétés thermodynamiques et cinétiques. Afin de comprendre les propriétés fondamentales de ces alliages, la première partie de ce travail est dédiée à l’étude ab-initio des effets de l’environnement chimique local sur l’état magnétique des solutés de Mn dans le Fe-Mn. Diverses configurations contenant du Mn, isolé ou sous forme d’amas, en présence de lacunes ou d’impuretés interstitielles sont étudiées et leur configuration magnétique de plus basse énergie est déterminée. Un modèle effectif d’interactions est paramétré à partir des données ab-initio afin d’étudier les propriétés des alliages Fe-Mn à température finie. Les propriétés clés sont identifiées et le modèle est validé à basse température en reproduisant les résultats ab-initio. L’utilisation de ce modèle couplé à des simulations Monte Carlo permet de simuler l’évolution chimique des alliages Fe-Mn en fonction de la température et de la concentration en Mn, tout en relaxant la structure magnétique en temps réel. Afin d’illustrer les possibles applications de ce modèle, diverses propriétés sont étudiées telles que la dépendance en concentration de la température de Curie ou encore l’évolution en température de l’énergie de mélange et de l’ordre atomique à courte distance. Puisque dans ces alliages la diffusion est en général régie par mécanisme lacunaire, nous proposons aussi un formalisme prenant en compte explicitement les effets de l’ordre magnétique local sur les propriétés des lacunes. Par simulations Monte Carlo de traceurs, cette approche prédit la dépendance en température de l’auto-diffusion dans le Fe en excellent accord avec les études expérimentales. La déviation de la loi d’Arrhénius proche de la température de Curie est directement prédite, ainsi que le changement de pente entre les régimes ferromagnétique et paramagnétique. La précision du modèle de Ruch, couramment utilisé dans la littérature, est discutée au vu des résultats obtenus. Enfin, cette approche est appliquée à la diffusion d’un soluté de Mn dans le Fe pur et comparée aux résultats expérimentaux. / In 3d alloys, magnetic properties of solutes can be extremely sensitive to local chemical environments and can have a crucial impact on various thermodynamic and kinetic properties. In order to properly understand the fundamental properties of these alloys, the first part of this work is dedicated to the study of the effects of local chemical environment on the magnetic state of Mn solutes in bcc Fe-Mn by means of Density Functional Theory. Namely, configurations containing Mn, being isolated or forming a cluster, and in the presence of vacancies or interstitial impurities are investigated and their lowest-energy magnetic configuration is determined. The ab-initio data produced are then used to parameterize an effective interaction model in order to study the properties of Fe-Mn alloys at finite temperature. The key features of Fe-Mn alloys are identified, and the model is validated at low temperature by reproducing ab-initio predictions. Using this model coupled to Monte Carlo simulations, we simulate the chemical evolution of Fe-Mn properties depending on temperature and Mn concentration while relaxing the magnetic structure on-the-fly. In order to illustrate the validity and the applicability of the model, we examine certain finite temperature properties of bcc Fe-Mn alloys such as the concentration dependence of the Curie temperature or the temperature evolution of the mixing energy and the atomic short-range order. Since diffusion in Fe and Fe-Mn alloys is generally ruled by vacancy-mechanism, we also propose a formalism to take explicitly into account the properties of vacancies in the interaction model and the effect of local magnetic state on these properties. Using tracer diffusion Monte Carlo simulations, this approach predicts the temperature dependence of self-diffusion in bcc Fe in excellent agreement with experimental results, including the deviation from Arrhenius law around the Curie temperature and the change of slope between the ferromagnetic and paramagnetic regimes. The accuracy of the widely used Ruch model is discussed in the light of the present results. Finally, we apply this approach to the diffusion of a Mn solute in bcc Fe and compare with experimental results.

Les cristaux liquides constituent un état de la matière intermédiaire entre les solides et les liquides. De ce fait, ils allient à la fois la fluidité des liquides et les propriétés anisotropes des solides. Ce qui les rend à la fois excitants et compliqués à étudier est leur grande diversité. En effet, il existe plus d'une dizaine de phases différentes, chacune caractérisée par une structure et des propriétés particulières. De plus, de petits changements dans la structure des molécules peuvent mener à un polymorphisme liquide cristallin modifié ou des intervalles de stabilité thermique bien différents. Au sein du laboratoire, pour essayer d'y voir plus clair et de comprendre quels sont les facteurs microscopiques qui influencent le polymorphisme liquide cristallin, une technique relativement récente a été utilisée. Il s'agit de la dynamique moléculaire au niveau atomistique. Cet outil permet d'étudier des molécules dans des conditions qui tendent à représenter au mieux des conditions expérimentales, avec une température, une pression, une densité, etc. En simulant des systèmes qui ont été caractérisés expérimentalement, des corrélations entre le comportement expérimental et simulé peuvent être extraites. L'hypothèse que nous souhaitons vérifier dans ces travaux est la suivante: les variations de l'énergie non liante de Coulomb, un terme énergétique qui peut être extrait des simulations, peuvent être corréler à des observations expérimentales. Cette hypothèse se base sur les travaux précédents effectués au sein du laboratoire sur la phase smectique C (SmC). Dans ces travaux, la plage de stabilité thermique de la phase SmC a été corrélée aux valeurs des énergies non liantes de Coulomb. Dans le cadre de cette étude, le même protocole de simulation est utilisé pour étudier la phase smectique A (SmA), dans le but de confirmer les observations faites sur la phase SmC. Pour ce faire, deux familles de mésogènes (molécules possédant une ou plusieurs phases liquides cristallines dans leurs polymorphisme) sont étudiées. Ces deux familles diffèrent par de petits détails structuraux qui ont pourtant une grande influence sur le polymorphisme expérimental. Il s'agit, dans le cas de ces deux familles, de la longueur d'une chaîne alkyle ou alkoxy sur le mésogène. Pour la première famille de mésogènes, le type de phase SmA observé est différent selon la longueur de la chaîne. Pour la deuxième famille, des mésogènes formés par liaison halogène, il a été observé expérimentalement que la longueur de la chaîne a un impact sur la plage de stabilité thermique de la phase SmA. Ces deux familles sont donc tout à fait appropriées pour une étude par simulation au niveau atomistique. En ce qui concerne les résultats de la première famille, le mésogène présentant une phase SmA différente des autres mésogènes est associé à une énergie de Coulomb également différente. Pour relier ces énergies à l'organisation des molécules au sein des phases, un raisonnement basé sur la dépendance en distance de cette interaction est utilisé. Si une attraction existe entre les mésogènes, une énergie relativement plus basse correspond à des mésogènes plus proches les uns des autres. Au contraire, si les molécules ont tendance à s'éloigner les unes des autres, une énergie relativement plus élevée correspond à une distance intermoléculaire plus faible. Nous avons pu tirer des conclusions quant aux positions relatives des mésogènes et montrer qu'elles pouvaient être reliées à la phase observée expérimentalement. Pour la deuxième famille, nous observons que le mésogène avec la plage de stabilité thermique la plus grande est celui pour lequel les énergies de Coulomb sont les plus faibles, donc celui pour lequel l'arrangement initial imposé (arrangement SmA) est le plus stable énergétiquement. Ces deux projets d'étude de la phase SmA montrent que l'hypothèse mise de l'avant dans les travaux précédents portant sur la phase SmC est vérifiée. La combinaison simulation/expérience fournit des données sur les deux échelles nécessaires, à savoir la molécule et la phase macroscopique, pour essayer d'identifier les liens existants entre ces deux mondes.

Les phyllosilicates sont des minéraux d'importance dans l'étude des roches métamorphiques et dans les sites de confinement des déchets liés aux activités anthropiques. Or, les modèles thermodynamiques actuels des argiles ne permettent pas de reproduire les variations de volume observées lors de leur déshydratation, ni leur changement de composition lors du métamorphisme. Les modèles thermodynamiques des micas sont en revanche bien contraints à haute température mais ne permettent pas des estimations thermobarométriques précises à moins de 350° C environ.<br />Nous proposons deux modèles thermodynamiques pour les smectites, illites, interstratifiés illites / smectites et micas. Ces modèles permettent respectivement de calculer la composition des phases stables à basse température, incluant les argiles, et l'estimation thermobarométrique de leurs conditions de cristallisation. Ils prennent tous deux en compte l'hydratation des phyllosilicates, variable en fonction de la composition, la pression, la température et l'activité d'eau. L'évaluation des propriétés thermodynamiques des équivalents hydratés des pôles purs des micas et des paramètres de solution solide nécessaires a été réalisée pour reproduire au mieux les contraintes expérimentales de déshydratation, nature des phases stables, calorimétrie ainsi que les estimations des conditions de cristallisation d'une compilation d'analyses, dont la gamme s'étend de la diagenèse à la haute pression - haute température. Ces modèles sont appliqués au calcul de diagramme de phase dans des systèmes chimiques simples et à l'estimation des conditions de cristallisation de phyllosilicates de nombreux échantillons naturels, y compris à partir de cartographies élémentaires.<br />Nous avons testé plusieurs méthodes d'estimation des propriétés thermodynamiques des phases pour lesquelles les contraintes expérimentales et / ou du milieu naturel sont insuffisantes. Aucune d'entre elles ne permet une estimation de l'enthalpie standard de formation directement utilisable à des fins thermobarométriques. Nous proposons cependant une approche pour améliorer cette situation. De plus, nous avons utilisé une méthode de simulation atomistique pour évaluer l'enthalpie de mélange le long de deux solutions solides d'intérêt pour la pétrologie de basse température. Les résultats sont compatibles avec les observations du milieu naturel, et le calcul du solvus entre les pôles muscovite et pyrophyllite confirme l'importance de l'hydratation pour la stabilité des argiles.

Du fait de l'augmentation de la demande de produits pour les applications grand public, industrielles et automobiles, des mémoires embarquées fiables et à faible coût de fabrication sont de plus en plus demandées. Dans ce contexte, les mémoires split-gate à piégeage discret sont proposées pour des microcontrôleurs. Elles combinent l'avantage d'une couche de stockage discrète et de la con guration split-gate. Durant ce travail de recherche, des mémoires split-gate à couche de piégeage discret ayant des longueurs de grille de 20nm sont présentées pour la première fois. Celles-ci on été réalisées avec des nanocristaux de silicium (Si-nc), du nitrure de silicium (SiN) ou un hybride Si-nc/SiN avec diélectrique de control de type SiO2 ou AlO et sont comparées en termes de performances lors des procédures d'eff acement et de rétention. Ensuite, la miniaturisation des mémoires split-gate à piégeage de charge est étudié, en particulier au travers de l'impact de la réduction de la longueur de grille sur la fenêtre de mémorisation, la rétention et la consommation. Le rôle des défauts dans le diélectrique de contrôle (alumine) utilisé dans les mémoires de type TANOS a été étudié. Des travaux ont été menés pour déterminer l'origine des pièges dans ce matériau, par le biais de la simulation atomistique ainsi que d'analyses physico-chimiques précises. Nous avons montré que la concentration de pièges dans AlO pouvait être réduite par ajustement des conditions de procédé de fabrication, débouchant ainsi sur l'amélioration de la rétention dans les mémoires à piégeage de charge. Ce résultat est convenable pour les applications de type embarqué

Du fait de l'augmentation de la demande de produits pour les applications grand public, industrielles et automobiles, des mémoires embarquées fiables et à faible coût de fabrication sont de plus en plus demandées. Dans ce contexte, les mémoires split-gate à piégeage discret sont proposées pour des microcontrôleurs. Elles combinent l'avantage d'une couche de stockage discrète et de la con guration split-gate. Durant ce travail de recherche, des mémoires split-gate à couche de piégeage discret ayant des longueurs de grille de 20nm sont présentées pour la première fois. Celles-ci on été réalisées avec des nanocristaux de silicium (Si-nc), du nitrure de silicium (SiN) ou un hybride Si-nc/SiN avec diélectrique de control de type SiO2 ou AlO et sont comparées en termes de performances lors des procédures d'eff acement et de rétention. Ensuite, la miniaturisation des mémoires split-gate à piégeage de charge est étudié, en particulier au travers de l'impact de la réduction de la longueur de grille sur la fenêtre de mémorisation, la rétention et la consommation. Le rôle des défauts dans le diélectrique de contrôle (alumine) utilisé dans les mémoires de type TANOS a été étudié. Des travaux ont été menés pour déterminer l'origine des pièges dans ce matériau, par le biais de la simulation atomistique ainsi que d'analyses physico-chimiques précises. Nous avons montré que la concentration de pièges dans AlO pouvait être réduite par ajustement des conditions de procédé de fabrication, débouchant ainsi sur l'amélioration de la rétention dans les mémoires à piégeage de charge. Ce résultat est convenable pour les applications de type embarqué / Due to the increasing demand for consumer, industrial and automotive products, highly reliable, and low integration cost embedded memories are more and more required. In this context, split-gate charge trap memories were proposed for microcontroller products, combining the advantage of a discrete storage layer and of the split-gate con guration. In this thesis, split-gate charge trap memories with electrical gate length down to 20nm are presented for the 1st time. Silicon nanocristals (Si-nc), or silicon nitride (SiN) and hybrid Si-nc/SiN based split-gate memories, with SiO2 or AlO control dielectrics, are compared in terms of program erase and retention. Then, the scalability of split-gate charge trap memories is studied, investigating the impact of gate length reduction on the memory window, retention and consumption. We thus studied the role of defects on alumina control dielectric employed in TANOS-like memory. We used atomistic simulation, consolidated by a detailed alumina physico-chemical material analysis, to investigate the origin of traps in alumina. We showed that the trap concentration in AlO can be decreased by adjusting the process conditions leading to improved retention behaviour in charge trap memory, suitable for embedded applications.

Dans l'industrie nucléaire, un verre borosilicaté complexe est utilisé pour confiner les produits de fission et les actinides mineurs. Sous irradiation, la structure et les propriétés mécaniques du verre évoluent. Dans cette étude, les simulations atomistiques et multi-échelles des trois verres borosilicatés simplifiés ont été lancées pour comprendre l'origine de l'évolution de leur comportement à la fracturation sous irradiation. Sous l'irradiation, l'élasticité diminue tandis que la plasticité augmente. La fracturation se produit suite à la formation et à la coalescence de cavités. Les modifications de structure sous l'irradiation conduisent à un retard de la coalescence et de la décohésion du verre irradié. Plusieurs phénomènes complexes se superposent pour expliquer ce comportement, notamment les modifications de la répartition des cavités, la mobilité des sodiums, et l'organisation des réseaux boratés et silicatés. Selon la nature du phénomène prépondérant, la ténacité peut augmenter ou diminuer sous irradiation.

La croissance de couches minces d'oxydes sur les alliages métalliques et les métaux purs est un phénomène ayant fait l'objet d'un grand nombre d'études expérimentales. Les structures des couches d'oxydes sont bien connues, ainsi que les modèles mathématiques servant à modéliser les aspects macroscopiques de la croissance. Cependant, les détails des mécanismes de germination de la couche d'oxyde dans les premiers stades de la corrosion ainsi que sa croissance ultérieure demeurent peu ou mal connus. La simulation atomistique apparaît comme une alternative pour évaluer les différents mécanismes proposés et appréhender l'influence des différents paramètres physico-chimiques. Le développement d'un tel outil de simulation a démarré au LPCS avec la thèse de M. Legrand. En se basant sur l'exemple de l'alliage FeCr, un modèle informatique tridimensionnel dit "modèle de Legrand", permettant de simuler la dissolution sélective et la passivation des alliages binaires a été réalisé. L'évolution dynamique est basée sur une technique de type Monte Carlo classique. Le logiciel permet de simuler l'évolution d'un alliage quelconque, d'une composition et d'une structure cristallographique donnée. Il prend en compte la diffusion des atomes sur la surface, leur dissolution et le blocage de la dissolution par formation d'une couche de passivation. Cet outil était adapté pour la simulation des premiers stades de la corrosion. L'objectif de ce travail est d'améliorer ce modèle existant, afin de simuler l'évolution de la couche passive sur une échelle de temps plus longue. A l'issue de ce travail, de nombreux apports ont été effectués. Ainsi, l'introduction d'un champ de force MEAM (Modified Embedded Atom Method) pour le calcul des barrières de diffusion et de dissolution, a permis de remplacer les probabilités de diffusion empiriques par des probabilités calculées, et de mettre en évidence la diffusion préférentielle des Cr vers leurs semblables. L'introduction d'une dynamique de simulation Monte Carlo Cinétique (KMC) a permis une prise en compte réaliste de l'évolution cinétique du modèle avec des temps de simulation reliés au temps réel. Enfin Un second réseau cristallographique RVO (réseau virtuel d'oxyde) tridimensionnel, correspondant à celui de la couche passive (Cr2O3) a été implémenté, ainsi qu'une interface graphique pour un meilleur suivi de la simulation. Les résultats obtenus lors des simulations sont en accord avec les observations expérimentales: passivation totale à partir du Fe-16Cr, enrichissement en Cr de la couche passive, allures des courbes cinétiques, influence du champ électrique, mise en évidence l'apparition de cavités sous la couche passive.

Le présent manuscrit présente le développement méthodologique du calcul de la tension de surface d’interfaces solide-liquide via des simulations de dynamique moléculaire. Après une courte présentation des avancées dans le domaine du calcul de la tension de surface pour les interfaces fluide-fluide et solide-fluide, les principales méthodes de calcul de la tension de surface d’un point de vue théorique sont montrées et généralisées pour le cas des interfaces solide-liquide, puis mises en oeuvre dans le cas de simulations de dynamique moléculaire. Un système école, constitué d’une feuille de graphène pour la phase solide et d’un bain de méthane pour la phase liquide, est ensuite étudié pour observer l’influence des artefacts de simulation sur le calcul de la tension de surface, montrant en particulier des effets de taille bien plus importants que pour le cas des interfaces liquide-liquide. Un autre système constitué d’une tranche de cuivre pour la phase solide, et d’un bain de méthane pour la phase liquide, a permis d’étudier l’effet inédit aux systèmes solide-liquide appelé anisotropie, montrant en particulier l’importance du caractère tensoriel de la tension de surface pour ce type de système. L’influence des paramètres du potentiel croisé entre les atomes de cuivre et de méthane est ensuite étudié. Finalement, deux systèmes applicatifs sont abordés, d’une part le système graphène-eau permettant d’étudier les effets de l’interaction électrostatique, et d’autre part un système constitué d’un solide explosif, le 1,3,5-triamino-2,4,6-trinitrobenzène (TATB) en contact avec un bain polymère pour la phase liquide, représentatif d’un cas réel d’intérêt. / This manuscript presents the methodological development of surface tension calculation of solid-liquidinterfaces via molecular dynamics simulations. After a short presentation of the advances in the field ofsurface tension calculation for fluid-fluid and solid-fluid interfaces, the main methods of surface tensioncalculation from a theoretical point of view are shown and generalized for solid-liquid interfaces, thenimplemented in the case of molecular dynamics simulations. A school system, consisting of a graphenesheet for the solid phase and a methane bath for the liquid phase, is then studied to observe the influenceof simulation artifacts on the surface tension calculation, showing in particular much larger size effectsthan in the case of liquid-liquid interfaces. Another system consisting of a copper slice for the solid phaseand a methane bath for the liquid phase made it possible to study the novel effect of solid-liquid systemscalled anisotropy, showing in particular the importance of the tensor character of the surface tension forthis type of system. The influence of the parameters of the cross potential between copper and methaneatoms is then studied. Finally, two application systems are discussed, on the one hand the graphene-watersystem for studying the effects of electrostatic interaction, and on the other hand a system consisting ofan explosive solid, 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) in contact with a polymer bath for theliquid phase, representing a real case of interest.

Dans le cadre de la vitrification de déchets nucléaires, la connaissance et la maîtrise des propriétés structurales et dynamiques des verres incluant les radioéléments sont importantes (dans le liquide et le solide). Elles influencent notamment la qualité du colis de verre, la durée de vie du procédé de vitrification et la quantité de produits de fission qu’il est possible d’introduire. Des simulations de dynamique moléculaire ont été réalisées afin d’analyser l’influence de la composition des matrices vitreuses sur les propriétés structurales et dynamiques. Le verre nucléaire industriel R7T7 étant composé d’un trop grand nombre d’oxydes (une trentaine) pour être simulé correctement, un verre simplifié, composé des oxydes SiO2, B2O3 et Na2O (majoritaires du verre R7T7) a été préféré. L’ajout de La2O3 permet de simuler l’impact des produits de fission et les actinides mineurs sur les propriétés de la matrice vitreuse. Les deux systèmes SiO2-B2O3-Na2O et SiO2-B2O3-Na2O-La2O3 ont permis d’étudier l’effet du sodium et du lanthane sur différentes propriétés. Au cours de ces travaux, un champ de force polarisable a été développé pour réaliser les simulations. Nos calculs à température ambiante ont permis de reproduire les résultats expérimentaux de la structure, de la répartition BIII/B IV, de la densité, avec un bon accord. Une étude a été menée dans le liquide pour l’analyse de la viscosité et de la conductivité électrique. La spéciation B IV/B III et l’influence des changements structuraux sur la densité avec l’augmentation de la température ont également été observées au cours de trempes thermiques. Les limites actuelles de cette approche sont également décrites. / As result of the nuclear waste vitrification, the knowledge and understanding of the dynamic and structural properties of glasses, including the behavior of radionuclides, is important (in liquid and solid phases). It can influence the glass waste properties, the lifetime of the vitrification process and the amount of radionuclides introduced in the glass matrix. Molecular dynamic simulations have been done to study the influence of the glass matrix composition into the structural and dynamic properties of the glass. A simplified glass, with 3 major oxides of the R7T7 glass such as SiO2, B2O3 and Na2O, have been used to simulate the R7T7 industrial nuclear glass (a 30 oxides glass). The inclusion of La2O3 allows us to simulate the impact of fission products and minor actinides into the properties of the glass matrix. Both systems, the SiO2-B2O3-Na2O and SiO2-B2O3-Na2O-La2O3, allow us to study the sodium and lanthanum effect on the properties of the glass. During this work, a polarizable force field has been developed to do these simulations. The results obtained at room temperature let us reproduce the experimental results of the structure, the distribution of BIII/BIV and the density. A study has been done on the viscosity and electrical conductivity of the liquid. The distribution BIV/BIII and the influence of the structural changes on the density along with the temperature have also been observed with thermal quenchings. The current limits of this approach are also described.