Capturing material toughness by molecular simulation: accounting for large yielding effects and limits

The inherent computational cost of molecular simulations limits their use to the study of nanometric systems with potentially strong size effects. In the case of fracture mechanics, size effects due to yielding at the crack tip can affect strongly the mechanical response of small systems. In this paper we consider two examples: a silica crystal for which yielding is limited to a few atoms at the crack tip, and a nanoporous polymer for which the process zone is about one order of magnitude larger. We perform molecular simulations of fracture of those materials and investigate in particular the system and crack size effects. The simulated systems are periodic with an initial crack. Quasi-static loading is achieved by increasing the system size in the direction orthogonal to the crack while maintaining a constant temperature. As expected, the behaviors of the two materials are significantly different. We show that the behavior of the silica crystal is reasonably well described by the classical framework of linear elastic fracture mechanics (LEFM). Therefore, one can easily upscale engineering fracture properties from molecular simulation results. In contrast, LEFM fails capturing the behavior of the polymer and we propose an alternative analysis based on cohesive crack zone models. We show that with a linear decreasing cohesive law, this alternative approach captures well the behavior of the polymer. Using this cohesive law, one can anticipate the mechanical behavior at larger scale and assess engineering fracture properties. Thus, despite the large yielding of the polymer at the scale of the molecular simulation, the cohesive zone analysis offers a proper upscaling methodology.MIT Energy InitiativeShell Oil CompanySchlumberger Limite

Surface Chemistry and Atomic-Scale Reconstruction of Kerogen–Silica Composites

Interest in gas shale, a novel source rock of natural gas, has increased tremendously in recent years. Better understanding of the kerogen-rock interaction is of crucial importance for efficient gas extraction and, hence, asset management. In this study, we explore the possible chemical bonds between kerogen and silica, one of the most predominant mineral constituents of gas shale, by means of quantum chemistry. Energetically favorable bond formation reactions are found between alcoholic hydroxyl, carboxylate, and aldehyde groups, as well as aliphatic double bonds of kerogen and the silica surface. The performance of a reactive force field was also assessed in a representative set of chemical reactions and found to be satisfactory. The potential impact of bond formation reactions between the two phases on the actual kerogen-silica interface is discussed as a function of the kerogen type, maturity, and density. Finally, a methodology aiming to reconstruct realistic kerogen-silica interfaces is presented