A methodology for determining the local mechanical properties of model atomistic glassy polymeric nanostructured...

We propose a methodology for calculating the distribution of the mechanical properties in model atomistic polymer-based nanostructured systems. The use of atomistic simulations is key in unravelling the fundamental mechanical behavior of composite materials. Most simulations involving the mechanical properties of polymer nanocomposites (PNCs) concern their global (average) properties, which are typically extracted by applying macroscopic strain on the boundaries of the simulation box and calculating the total (global) stress by invoking the Virial formalism over all atoms within the simulation box; thus, extracting the pertinent mechanical properties from the corresponding stress-strain relation. However, in order to probe the distribution of mechanical properties within heterogeneous multi-component polymer-based systems, a detailed computation of stress and strain fields within specific sub-domains is necessary. For example, it is well known for multi-component nanostructured systems, such as PNCs, that the mechanical behavior of the polymer/nanofiller interphases, or interfaces, is crucial for determining the global mechanical properties of the composite materials. Here we propose a new methodology to probe the distribution of mechanical properties by directly computing the (local) stress and strain at the atomic level, and averaging over user-defined subdomains. The workflow of our computational method possesses the following features:•Calculating the stress and strain per atom (or per particle) for nanostructured microscopic (here atomistic) model configurations, under an imposed applied deformation.•Averaging the local, per-atom defined, stress and strain on user-defined subdomains within the nanostructured model system.•Predicting the mechanical properties within the specific subdomains, focusing on polymer/solid interphases.