Data type: Software (2026)
Availability: Public
Licence: Creative Commons Attribution 4.0 International
Molecular dynamics (MD) simulations are widely used to provide insights into fracture mechanisms while maintaining chemical specificity. However, particle-based techniques such as MD are limited in terms of accessible length scales and applicable boundary conditions, which restricts the investigation of fracture phenomena in typical engineering settings. In an attempt to overcome these limitations, we apply the partitioned-domain Capriccio method to couple atomistic MD samples representing silica glass with the finite element (FE) method. With this approach, we perform mode I (rectangular panel under tension, three-, and four-point bending), mode II as well as mode III (rectangular panel under in-plane or out-of-plane shear) simulations. Thereby, we investigate multiple criteria to identify the onset of crack propagation based on the virial stress, the number of pair interactions, the kinetic energy/temperature, the crack velocity, and the crack opening displacement. The approach presented provides quantitatively plausible results for the critical stress intensity factors KIc, KIIc, and KIIIc. This contribution shows that the Capriccio method is a flexible means of performing fracture simulations that take into account boundary conditions typical of experimental test setups with atomistic specificity near the crack tip. While also pointing out the current limitations of the Capriccio method, we demonstrate its potential to integrate atomistic insights into FE models with significantly larger overall dimensions.