Structural solution using molecular dynamics: Fundamentals and a case study of epoxy-silica interface

In this paper, the molecular dynamics (MD) simulation technique is described in the context of structural mechanics applications, providing a fundamental understanding of the atomistic approach, and demonstrating its applicability. Atomistic models provide a bottom-up description of material properties and processes, and MD simulation is capable of solving the dynamic evolution of equilibrium and non-equilibrium processes. The applicability of the technique to structural engineering problems is demonstrated through an interface debonding problem in a multi-layered material system usually encountered in composite structures. Interface debonding may lead to a possible premature failure of fiber reinforced polymer (FRP) bonded reinforced concrete (RC) structural elements subjected to moisture. Existing knowledge on meso-scale fracture mechanics may not fully explain the weakening of the interface between concrete and epoxy, when the interface is under moisture; there is a need to study the moisture affected debonding of the interface using a more fundamental approach that incorporates chemistry in the description of materials. The results of the atomistic modeling presented in this paper show that the adhesive strength (in terms of energy) between epoxy and silica is weakened in the presence of water through its interaction with epoxy. This is correlated with the existing meso-scale experimental data. This example demonstrates that MD simulation can be effectively used in studying the durability of the system through an understanding of how materials interact with the environment at the molecular level. In view of the limitation of MD simulation on both length- and time-scales, future research may focus on the development of a bridging technique between MD and finite element modeling (FEM) to be able to correlate the results from the nano- to the macro-scale.

[1]  Philippe H. Geubelle,et al.  The elastic modulus of single-wall carbon nanotubes: a continuum analysis incorporating interatomic potentials , 2002 .

[2]  Oral Büyüköztürk,et al.  Fracture characterization of concrete/epoxy interface affected by moisture , 2010 .

[3]  Tomotsugu Shimokawa,et al.  Matching conditions in the quasicontinuum method: Removal of the error introduced at the interface between the coarse-grained and fully atomistic region , 2004 .

[4]  E. B. Tadmor,et al.  Quasicontinuum models of interfacial structure and deformation , 1998 .

[5]  Markus J. Buehler,et al.  Fracture mechanics of protein materials , 2007 .

[6]  O. Büyüköztürk,et al.  Debonding of FRP plated concrete: A tri-layer fracture treatment , 2006 .

[7]  C. Brooks.  Computer simulation of liquids , 1989 .

[8]  D. Osguthorpe,et al.  Structure and energetics of ligand binding to proteins: Escherichia coli dihydrofolate reductase‐trimethoprim, a drug‐receptor system , 1988, Proteins.

[9]  Ronald E. Miller,et al.  Atomistic/continuum coupling in computational materials science , 2003 .

[10]  Oral Büyüköztürk,et al.  Peel and Shear Fracture Characterization of Debonding in FRP Plated Concrete Affected by Moisture , 2006 .

[11]  M. Ortiz,et al.  An adaptive finite element approach to atomic-scale mechanics—the quasicontinuum method , 1997, cond-mat/9710027.

[12]  Huajian Gao,et al.  Crack nucleation and growth as strain localization in a virtual-bond continuum , 1998 .

[13]  Harold S. Park,et al.  A surface Cauchy–Born model for nanoscale materials , 2006 .

[14]  G. I. Bell.  Models for the specific adhesion of cells to cells. , 1978, Science.

[15]  J. Brédas,et al.  Molecular Dynamics Simulations of Intercalated Poly(ε-Caprolactone)-Montmorillonite Clay Nanocomposites , 2004 .

[16]  Ellad B. Tadmor,et al.  Hybrid Continuum Mechanics and Atomistic Methods for Simulating Materials Deformation and Failure , 2007 .

[17]  R.K. Kalia,et al.  Multiscale simulation of nanosystems , 2001, Comput. Sci. Eng..

[18]  M. Buehler,et al.  Energy landscape, structure and rate effects on strength properties of alpha-helical proteins , 2010, Journal of physics. Condensed matter : an Institute of Physics journal.

[19]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[20]  Ronald E. Miller,et al.  The Quasicontinuum Method: Overview, applications and current directions , 2002 .

[21]  M. Ortiz,et al.  An analysis of the quasicontinuum method , 2001, cond-mat/0103455.

[22]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[23]  M.G.B. Drew,et al.  The art of molecular dynamics simulation , 1996 .

[24]  C. Catlow,et al.  Potential models for ionic oxides , 1985 .

[25]  W F van Gunsteren,et al.  Molecular mechanics in biology: from structure to function, taking account of solvation. , 1994, Annual review of biophysics and biomolecular structure.

[26]  F. Momany,et al.  Validation of the general purpose QUANTA ®3.2/CHARMm® force field , 1992 .

[27]  A. Hagler,et al.  Derivation of force fields for molecular mechanics and dynamics from ab initio energy surfaces. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[28]  Arnold T. Hagler,et al.  On the use of quantum energy surfaces in the derivation of molecular force fields , 1994 .

[29]  Ming-Jing Hwang,et al.  Derivation of Class II Force Fields. 2. Derivation and Characterization of a Class II Force Field, CFF93, for the Alkyl Functional Group and Alkane Molecules , 1994 .

[30]  B. Alder,et al.  Phase Transition for a Hard Sphere System , 1957 .

[31]  Norman A. Fleck,et al.  Brittle fracture of adhesive joints , 1992 .

[32]  Harold S. Park,et al.  Surface stress effects on the resonant properties of metal nanowires: The importance of finite deformation kinematics and the impact of the residual surface stress , 2008 .

[33]  Markus J. Buehler,et al.  Hierarchies, multiple energy barriers, and robustness govern the fracture mechanics of α-helical and β-sheet protein domains , 2007, Proceedings of the National Academy of Sciences.

[34]  A. Voter.  Hyperdynamics: Accelerated Molecular Dynamics of Infrequent Events , 1997 .

[35]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[36]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[37]  Markus J Buehler,et al.  A realistic molecular model of cement hydrates , 2009, Proceedings of the National Academy of Sciences.

[38]  Maurizio Fermeglia,et al.  Computer simulation of nylon-6/organoclay nanocomposites: prediction of the binding energy , 2003 .

[39]  Harold S. Park,et al.  Bridging Scale Methods for Nanomechanics and Materials , 2006 .

[40]  Noam Bernstein,et al.  Mixed finite element and atomistic formulation for complex crystals , 1999 .

[41]  Yuuki Ryoji,et al.  Stress based criterion for an interface crack kinking out of the interface in dissimilar materials , 1992 .

[42]  Stephanie L. Walkup,et al.  Guide for the Design and Construction of Externally Bonded FRP Systems for Strengthening Concrete Structures (ACI 440.2R-02) , 2005 .

[43]  M A Ryan,et al.  Molecular modeling of polymer composite-analyte interactions in electronic nose sensors. , 2003, Sensors and actuators. B, Chemical.

[44]  D. Beveridge,et al.  Molecular dynamics studies of axis bending in d(G5-(GA4T4C)2-C5) and d(G5-(GT4A4C)2-C5): effects of sequence polarity on DNA curvature. , 1999, Journal of molecular biology.

[45]  Markus J. Buehler,et al.  Atomistically Informed Mesoscale Model of Alpha-Helical Protein Domains , 2009 .

[46]  M. Ortiz,et al.  Quasicontinuum analysis of defects in solids , 1996 .

[47]  T. Nishioka,et al.  Concepts of Separated J-Integrals, Separated Energy Release Rates, and the Component Separation Method of the J-Integral for Interfacial Fracture Mechanics , 2003 .

[48]  Terry R. Stouch,et al.  Effects of Switching Functions on the Behavior of Liquid Water in Molecular Dynamics Simulations , 1994 .

[49]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[50]  Jacob Fish,et al.  Bridging the scales in nano engineering and science , 2006 .

[51]  Maurizio Fermeglia,et al.  Estimation of the Binding Energy in Random Poly(Butylene terephtalate-co-thiodiethylene terephtalate) Copolyesters/Clay Nanocomposites via Molecular Simulation , 2004 .

[52]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[53]  J. R. Grigera,et al.  The effect of stereochemistry upon carbohydrate hydration. A molecular dynamics simulation of β-d-galactopyranose and (α,β)-d-talopyranose , 1994 .

[54]  Ming-Jing Hwang,et al.  Derivation of class II force fields. I. Methodology and quantum force field for the alkyl functional group and alkane molecules , 1994, J. Comput. Chem..

[55]  Markus J. Buehler,et al.  Atomistic Modeling of Materials Failure , 2008 .

[56]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[57]  O. Büyüköztürk,et al.  Progress on understanding debonding problems in reinforced concrete and steel members strengthened using FRP composites , 2004 .

[58]  S. L. Mayo,et al.  DREIDING: A generic force field for molecular simulations , 1990 .

[59]  Z. Suo,et al.  Mixed mode cracking in layered materials , 1991 .

[60]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[61]  Klaus Fiedler,et al.  An Extension of the Consistent Valence Force Field (CVFF) with the Aim to Simulate the Structures of Vanadium Phosphorus Oxides and the Adsorption of n‐Butane and of 1‐Butene on their Crystal Planes , 2002 .

[62]  F. Tsobnang,et al.  Molecular mechanics studies of specific interactions in organomodified clay nanocomposite , 2005 .