A generalized Irving-Kirkwood formula for the calculation of stress in molecular dynamics models.

In non-equilibrium molecular dynamics simulations, continuum mechanics quantities can be computed from the position and momentum of the particles based on the classical Irving-Kirkwood formalism. For practical purposes, the implementations of Irving-Kirkwood formulas often involve a spatial averaging using a smooth kernel function. The resulting formula for the stress has been known as Hardy stress. Usually results obtained this way still need to be further processed to reduce the fluctuation, e.g., by ensemble or time averaging. In this paper we extend Hardy's formulas by systematically incorporating both spatial and temporal averaging into the expression of continuum quantities. The derivation follows the Irving-Kirkwood formalism, and the average quantities still satisfy conservation laws in continuum mechanics. We will discuss the selection of kernel functions and present several numerical tests.

[1]  Xiaobo Nie,et al.  Resolving singular forces in cavity flow: multiscale modeling from atomic to millimeter scales. , 2005, Physical review letters.

[2]  T. J. Delph,et al.  Stress calculation in atomistic simulations of perfect and imperfect solids , 2001 .

[3]  Harold W. Hatch,et al.  Molecular modeling of mechanical stresses on proteins in glassy matrices: formalism. , 2012, The Journal of chemical physics.

[4]  T. Delph.  Near-surface stresses in silicon(001) , 2008 .

[5]  Seth Root,et al.  Continuum predictions from molecular dynamics simulations: Shock waves , 2003 .

[6]  J. Kirkwood,et al.  The Statistical Mechanical Theory of Transport Processes. IV. The Equations of Hydrodynamics , 1950 .

[7]  C. Dafermos.  Hyberbolic Conservation Laws in Continuum Physics , 2000 .

[8]  White,et al.  Molecular-dynamics simulations of void collapse in shocked model-molecular solids. , 1994, Physical review. B, Condensed matter.

[9]  D. H. Tsai.  The virial theorem and stress calculation in molecular dynamics , 1979 .

[10]  Xiaobo Nie,et al.  A continuum and molecular dynamics hybrid method for micro- and nano-fluid flow , 2004, Journal of Fluid Mechanics.

[11]  J. Molinari,et al.  Multiscale modeling of two-dimensional contacts. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  Gary P. Morriss,et al.  Statistical Mechanics of Nonequilibrium Liquids , 2008 .

[13]  Jonathan A. Zimmerman,et al.  Calculation of stress in atomistic simulation , 2004 .

[14]  Youping Chen,et al.  Coarse-grained simulations of single-crystal silicon , 2009 .

[15]  B. Holian,et al.  Test of a new heat-flow equation for dense-fluid shock waves. , 2010, The Journal of chemical physics.

[16]  Youping Chen,et al.  Connecting molecular dynamics to micromorphic theory. (I). Instantaneous and averaged mechanical variables , 2003 .

[17]  David J. Srolovitz,et al.  Local stress calculation in simulations of multicomponent systems , 2009, J. Comput. Phys..

[18]  Mark F. Horstemeyer,et al.  Atomistic Finite Deformation Simulations: A Discussion on Length Scale Effects in Relation to Mechanical Stresses , 1999 .

[19]  M. Baskes,et al.  Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals , 1984 .

[20]  Berend Smit,et al.  Molecular Dynamics Simulations , 2002 .

[21]  H. Johnson,et al.  Stress evolution to steady state in ion bombardment of silicon , 2008 .

[22]  Dick Bedeaux,et al.  ON THE PHYSICAL INTERPRETATION OF FIELDS IN CONTINUUM MECHANICS , 1993 .

[23]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[24]  Dick Bedeaux,et al.  Continuum equations of balance via weighted averages of microscopic quantities , 1994, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[25]  Youping Chen,et al.  Connecting molecular dynamics to micromorphic theory. (II). Balance laws , 2003 .

[26]  Björn Engquist,et al.  Heterogeneous multiscale methods for stiff ordinary differential equations , 2005, Math. Comput..

[27]  Min Zhou,et al.  A new look at the atomic level virial stress: on continuum-molecular system equivalence , 2003, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[28]  O. Kum.  Orientation effects in shocked nickel single crystals via molecular dynamics , 2001 .

[29]  T.J Delph.  Local stresses and elastic constants at the atomic scale , 2005, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[30]  E. Weinan,et al.  Heterogeneous multiscale method for the modeling of complex fluids and micro-fluidics , 2005 .

[31]  R. Hardy,et al.  Formulas for determining local properties in molecular‐dynamics simulations: Shock waves , 1982 .

[32]  Youping Chen.  Local stress and heat flux in atomistic systems involving three-body forces. , 2006, The Journal of chemical physics.

[33]  Ping Sheng,et al.  Molecular scale contact line hydrodynamics of immiscible flows. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[34]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[35]  Matthew Tirrell,et al.  Molecular dynamics of flow in micropores , 1987 .

[36]  Ellad B. Tadmor,et al.  A Unified Interpretation of Stress in Molecular Systems , 2010, 1008.4819.

[37]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[38]  S. Varadhan,et al.  Hydrodynamical limit for a Hamiltonian system with weak noise , 1993 .

[39]  Richard B. Lehoucq,et al.  Translation of Walter Noll’s “Derivation of the Fundamental Equations of Continuum Thermodynamics from Statistical Mechanics” , 2010 .

[40]  J. Barrat,et al.  On the study of local-stress rearrangements during quasi-static plastic shear of a model glass: Do local-stress components contain enough information? , 2007, The European physical journal. E, Soft matter.

[41]  H. Spohn.  Large Scale Dynamics of Interacting Particles , 1991 .

[42]  David L. McDowell,et al.  Coarse-grained atomistic simulation of dislocations , 2011 .

[43]  Mark F. Horstemeyer,et al.  A large deformation atomistic study examining crystal orientation effects on the stress–strain relationship , 2002 .

[44]  Phil Attard,et al.  Molecular dynamics study of the effect of atomic roughness on the slip length at the fluid-solid boundary during shear flow. , 2004, Langmuir : the ACS journal of surfaces and colloids.

[45]  E Weinan,et al.  A multiscale coupling method for the modeling of dynamics of solids with application to brittle cracks , 2010, J. Comput. Phys..

[46]  Seth Root,et al.  Continuum Properties from Molecular Simulations , 2002 .

[47]  C. T. White,et al.  Orientation dependence of shock-induced chemistry in diamond , 2002 .

[48]  James B. Adams,et al.  Interatomic Potentials from First-Principles Calculations: The Force-Matching Method , 1993, cond-mat/9306054.

[49]  Reese E. Jones,et al.  A material frame approach for evaluating continuum variables in atomistic simulations , 2008, J. Comput. Phys..

[50]  R. Clausius,et al.  XVI. On a mechanical theorem applicable to heat , 1870 .

[51]  E Weinan,et al.  Multiscale modeling of the dynamics of solids at finite temperature , 2005 .