1
|
Candela R, Gelin S, Mousseau N, Veiga RGA, Domain C, Perez M, Becquart CS. Investigating the kinetics of the formation of a C Cottrell atmosphere around a screw dislocation in bcc iron: a mixed-lattice atomistic kinetic Monte-Carlo analysis. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:065704. [PMID: 33137791 DOI: 10.1088/1361-648x/abc6c3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We present a mixed-lattice atomistic kinetic Monte-Carlo algorithm (MLKMC) that integrates a rigid-lattice AKMC approach with the kinetic activation-relaxation technique (k-ART), an off-lattice/self-learning AKMC. This approach opens the door to study large and complex systems adapting the cost of identification and evaluation of transition states to the local environment. To demonstrate its capacity, MLKMC is applied to the problem of the formation of a C Cottrell atmosphere decorating a screw dislocation in α-Fe. For this system, transitions that occur near the dislocation core are searched by k-ART, while transitions occurring far from the dislocation are computed before the simulation starts using the rigid-lattice AKMC. This combination of the precision of k-ART and the speed of the rigid-lattice makes it possible to follow the onset of the C Cottrell atmosphere and to identify interesting mechanisms associated with its formation.
Collapse
Affiliation(s)
- R Candela
- Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207-UMET-Unité Matériaux et Transformations, F-59000 Lille, France
- Laboratoire commun EDF-CNRS Etude et Modélisation des Microstructures pour le Vieillissement des, Matériaux (EM2VM), France
| | - S Gelin
- Département de Mathématiques et de Génie Industriel and Regroupement Québécois sur les Matériaux de Pointe, Polytechnique Montréal, Case Postale 6079, Succursale Centre-Ville, Montreal (QC) H3C 3A7, Canada
| | - N Mousseau
- Département de Physique and Regroupement Québécois sur les Matériaux de Pointe, Université de Montréal, Case Postale 6128, Succursale Centre-Ville, Montreal (QC) H3C 3J7, Canada
| | - R G A Veiga
- Universidade Federal do ABC, Center for Engineering, Modeling, and Social Applied Sciences (CECS), Av dos Estados, 5001, Santa Terezinha, CEP 09210-580, Santo André/SP, Brazil
| | - C Domain
- Laboratoire commun EDF-CNRS Etude et Modélisation des Microstructures pour le Vieillissement des, Matériaux (EM2VM), France
- EDF-R&D, Département Matériaux et Mécanique des Composants (MMC), Les Renardières, F-77818, Moret sur Loing Cedex, France
| | - M Perez
- Université de Lyon, INSA-Lyon, MATEIS, UMR CNRS 5510, 25 Avenue Jean Capelle, F-69621, Villeurbanne, France
| | - C S Becquart
- Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207-UMET-Unité Matériaux et Transformations, F-59000 Lille, France
- Laboratoire commun EDF-CNRS Etude et Modélisation des Microstructures pour le Vieillissement des, Matériaux (EM2VM), France
| |
Collapse
|
2
|
Jay A, Huet C, Salles N, Gunde M, Martin-Samos L, Richard N, Landa G, Goiffon V, De Gironcoli S, Hémeryck A, Mousseau N. Finding Reaction Pathways and Transition States: r-ARTn and d-ARTn as an Efficient and Versatile Alternative to String Approaches. J Chem Theory Comput 2020; 16:6726-6734. [DOI: 10.1021/acs.jctc.0c00541] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Antoine Jay
- LAAS-CNRS, Université de Toulouse, CNRS, F-31555 Toulouse, France
| | - Christophe Huet
- ISAE-SUPAERO, Université de Toulouse, F-31555 Toulouse, France
| | - Nicolas Salles
- CNR-IOM, Democritos National Simulation Center, Istituto Officina dei Materiali, c/o SISSA, via Bonomea 265, IT-34136 Trieste, Italy
| | - Miha Gunde
- LAAS-CNRS, Université de Toulouse, CNRS, F-31555 Toulouse, France
| | - Layla Martin-Samos
- CNR-IOM, Democritos National Simulation Center, Istituto Officina dei Materiali, c/o SISSA, via Bonomea 265, IT-34136 Trieste, Italy
| | | | - Georges Landa
- LAAS-CNRS, Université de Toulouse, CNRS, F-31555 Toulouse, France
| | - Vincent Goiffon
- ISAE-SUPAERO, Université de Toulouse, F-31555 Toulouse, France
| | - Stefano De Gironcoli
- CNR-IOM, Democritos National Simulation Center, Istituto Officina dei Materiali, c/o SISSA, via Bonomea 265, IT-34136 Trieste, Italy
- SISSA, via Bonomea 265, IT-34136 Trieste, Italy
| | - Anne Hémeryck
- LAAS-CNRS, Université de Toulouse, CNRS, F-31555 Toulouse, France
| | - Normand Mousseau
- Département de Physique and Regroupement québécois sur les matriaux de pointe, Département de Physique, Université de Montréal, C.P. 6128, succursale Centre-ville H3C 3J7 Montréal Canada Montréal, Canada
| |
Collapse
|
3
|
Dufresne A, Arayro J, Zhou T, Ioannidou K, Ulm FJ, Pellenq R, Béland LK. Atomistic and mesoscale simulation of sodium and potassium adsorption in cement paste. J Chem Phys 2018; 149:074705. [PMID: 30134717 DOI: 10.1063/1.5042755] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
An atomistic and mesoscopic assessment of the effect of alkali uptake in cement paste is performed. Semi-grand canonical Monte Carlo simulations indicate that Na and K not only adsorb at the pore surface of calcium silicate hydrates (C-S-H) but also adsorb in the C-S-H hydrated interlayer up to concentrations of the order of 0.05 and 0.1 mol/kg, respectively. Sorption of alkali is favored as the Ca/Si ratio of C-S-H is reduced. Long timescale simulations using the Activation Relaxation Technique indicate that characteristic diffusion times of Na and K in the C-S-H interlayer are of the order of a few hours. At the level of individual grains, Na and K adsorption leads to a reduction of roughly 5% of the elastic moduli and to volume expansion of about 0.25%. Simulations using the so-called primitive model indicate that adsorption of alkali ions at the pore surface can reduce the binding between C-S-H grains by up to 6%. Using a mesoscopic model of cement paste, the combination of individual grain swelling and changes in inter-granular cohesion was estimated to lead to overall expansive pressures of up to 4 MPa-and typically of less than 1 MPa-for typical alkali concentrations observed at the proximity of gel veins caused by the alkali-silica reaction.
Collapse
Affiliation(s)
- Alice Dufresne
- Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Jack Arayro
- MultiScale Materials Science for Energy and Environment (MSE2), The Joint CNRS/MIT/Aix-Marseille University Laboratory, UMI CNRS 3466, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Tingtao Zhou
- Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Katerina Ioannidou
- Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Franz-Josef Ulm
- Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Roland Pellenq
- Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Laurent Karim Béland
- MultiScale Materials Science for Energy and Environment (MSE2), The Joint CNRS/MIT/Aix-Marseille University Laboratory, UMI CNRS 3466, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
4
|
Trochet M, Sauvé-Lacoursière A, Mousseau N. Algorithmic developments of the kinetic activation-relaxation technique: Accessing long-time kinetics of larger and more complex systems. J Chem Phys 2017; 147:152712. [DOI: 10.1063/1.4995426] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mickaël Trochet
- Département de Physique and Regroupement Québécois sur les Matériaux de Pointe, Université de Montréal, C.P. 6128, Succursale Centre-Ville, Montréal H3C 3J7, Québec, Canada
| | - Alecsandre Sauvé-Lacoursière
- Département de Physique and Regroupement Québécois sur les Matériaux de Pointe, Université de Montréal, C.P. 6128, Succursale Centre-Ville, Montréal H3C 3J7, Québec, Canada
| | - Normand Mousseau
- Département de Physique and Regroupement Québécois sur les Matériaux de Pointe, Université de Montréal, C.P. 6128, Succursale Centre-Ville, Montréal H3C 3J7, Québec, Canada
| |
Collapse
|
5
|
Direct Observation of Defect Range and Evolution in Ion-Irradiated Single Crystalline Ni and Ni Binary Alloys. Sci Rep 2016; 6:19994. [PMID: 26829570 PMCID: PMC4734288 DOI: 10.1038/srep19994] [Citation(s) in RCA: 110] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2015] [Accepted: 12/22/2015] [Indexed: 12/25/2022] Open
Abstract
Energetic ions have been widely used to evaluate the irradiation tolerance of structural materials for nuclear power applications and to modify material properties. It is important to understand the defect production, annihilation and migration mechanisms during and after collision cascades. In this study, single crystalline pure nickel metal and single-phase concentrated solid solution alloys of 50%Ni50%Co (NiCo) and 50%Ni50%Fe (NiFe) without apparent preexisting defect sinks were employed to study defect dynamics under ion irradiation. Both cross-sectional transmission electron microscopy characterization (TEM) and Rutherford backscattering spectrometry channeling (RBS-C) spectra show that the range of radiation-induced defect clusters far exceed the theoretically predicted depth in all materials after high-dose irradiation. Defects in nickel migrate faster than in NiCo and NiFe. Both vacancy-type stacking fault tetrahedra (SFT) and interstitial loops coexist in the same region, which is consistent with molecular dynamics simulations. Kinetic activation relaxation technique (k-ART) simulations for nickel showed that small vacancy clusters, such as di-vacancies and tri-vacancies, created by collision cascades are highly mobile, even at room temperature. The slower migration of defects in the alloy along with more localized energy dissipation of the displacement cascade may lead to enhanced radiation tolerance.
Collapse
|
6
|
N'Tsouaglo GK, Béland LK, Joly JF, Brommer P, Mousseau N, Pochet P. Probing Potential Energy Surface Exploration Strategies for Complex Systems. J Chem Theory Comput 2015; 11:1970-7. [PMID: 26574398 DOI: 10.1021/ct501032v] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The efficiency of minimum-energy configuration searching algorithms is closely linked to the energy landscape structure of complex systems, yet these algorithms often include a number of steps of which the effect is not always clear. Decoupling these steps and their impacts can allow us to better understand both their role and the nature of complex energy landscape. Here, we consider a family of minimum-energy algorithms based, directly or indirectly, on the well-known Bell-Evans-Polanyi (BEP) principle. Comparing trajectories generated with BEP-based algorithms to kinetically correct off-lattice kinetic Monte Carlo schemes allow us to confirm that the BEP principle does not hold for complex systems since forward and reverse energy barriers are completely uncorrelated. As would be expected, following the lowest available energy barrier leads to rapid trapping. This is why BEP-based methods require also a direct handling of visited basins or barriers. Comparing the efficiency of these methods with a thermodynamical handling of low-energy barriers, we show that most of the efficiency of the BEP-like methods lie first and foremost in the basin management rather than in the BEP-like step.
Collapse
Affiliation(s)
- Gawonou Kokou N'Tsouaglo
- Département de physique and Regroupement québécois sur les matériaux de pointe, Université de Montréal , C.P. 6128, Succursale Centre-Ville, Montréal, H3C 3J7 Québec Canada
| | - Laurent Karim Béland
- Département de physique and Regroupement québécois sur les matériaux de pointe, Université de Montréal , C.P. 6128, Succursale Centre-Ville, Montréal, H3C 3J7 Québec Canada
| | - Jean-François Joly
- Département de physique and Regroupement québécois sur les matériaux de pointe, Université de Montréal , C.P. 6128, Succursale Centre-Ville, Montréal, H3C 3J7 Québec Canada
| | - Peter Brommer
- Centre for Predictive Modelling, School of Engineering, University of Warwick , Coventry CV4 7AL, United Kingdom
| | - Normand Mousseau
- Département de physique and Regroupement québécois sur les matériaux de pointe, Université de Montréal , C.P. 6128, Succursale Centre-Ville, Montréal, H3C 3J7 Québec Canada.,Laboratoire de Physique Théorique de la Matière Condensée, Université Pierre et Marie Curie , Boite 121, 4, Place Jussieu, 75252 Paris Cedex 05, France
| | - Pascal Pochet
- Université de Grenoble Alpes , INAC-SP2M, L_Sim, F-38000 Grenoble, France.,CEA, INAC-SP2M, Atomistic Simulation Laboratory, F-38000 Grenoble, France
| |
Collapse
|