Organic–inorganic interfaces exist in many natural or synthetic materials, such as mineral–protein interfaces found in bone and epoxy–silica interfaces found in concrete construction. Here, we report a model to predict the intrinsic strength between organic and inorganic materials, based on a molecular dynamics simulation approach combined with the metadynamics method, used to reconstruct the free energy surface between attached and detached states of the bonded system and scaled up to incorporate it into a continuum model. We apply this technique to model an epoxy–silica system that primarily features nonbonded and nondirectional van der Waals and Coulombic chemical interactions. The intrinsic strength between epoxy and silica derived from the molecular level is used to predict the structural behavior of epoxy–silica interface at the macroscopic length scale by invoking a finite element approach using a cohesive zone model which shows a good agreement with existing experimental results.
