#
APS March Meeting 2013

## Volume 58, Number 1

##
Monday–Friday, March 18–22, 2013;
Baltimore, Maryland

### Session G11: Invited Session: Concurrent Multiple Length-Scale Modelling

11:15 AM–2:15 PM,
Tuesday, March 19, 2013

Room: 310

Sponsoring
Unit:
DCOMP

Chair: Qing Peng, Rensselear Polytechnical Institute

Abstract ID: BAPS.2013.MAR.G11.4

### Abstract: G11.00004 : Coarse-graining molecular dynamics models using an extended Galerkin method*

1:03 PM–1:39 PM

Preview Abstract
View Presentation
Abstract

####
Author:

Xiantao Li

(The Pennsylvania State University)

I will present a systematic approach to coarse-grain molecular dynamics models for solids. The coarse-grained models are derived by Galerkin projection to a sequence of Krylov subspaces. On the coarsest space, the model corresponds to a finite element discretization of the continuum elasto-dynamics model. On the other hand, the projection to the finest space yields the full molecular dynamics description. The models in between serve as a smooth transition between the two scales.
We start with a molecular dynamics (MD) model, $m_i\ddot{\mathbf x}_i= -\frac{\partial V}{\partial \mathbf x_i}$.
First, let $Y_0$ be the approximation space for the continuum model. By projecting the MD model onto the subspace, we obtain a coarse-grained model, $ M \ddot{\mathbf q} = F(\mathbf q)$.
Using the Cauchy-Born approximation, this model can be shown to coincide with
the finite element representation of the continuum elastodynamics model.
This model has limited accuracy near lattice defects. One natural idea is to switch to the MD model in regions surround local defect. As a result, one creates an interface between
the continuum and atomistic description, where coupling conditions are needed. Direct coupling methods may involve enforcing constraints or mixing the energy or forces. Such an approach
may suffer from large phonon reflections at the interface, and introduce large modeling error.
In order to seamlessly couple this model to MD, we successively expand the approximation space to the Krylov spaces,
$ K_\ell = Y_0 + A Y_0 + \cdots + A^\ell Y_0$.
Here $A$ is the force constant matrix, computed from the atomistic model. Due to the translational invariance, only a smaller number of such matrices need to be computed. By projecting the MD model onto this new subspace, we obtain an extended system,
$M \ddot{\mathbf q} = F_0(\mathbf q, \xi_1, \cdots, \xi_{\ell}),
\ddot{\xi}_1= F_1(\mathbf q, \xi_1, \cdots, \xi_{\ell}),
\cdots \cdots,
\ddot{\xi}_{\ell}= F_\ell(\mathbf q, \xi_1, \cdots, \xi_{\ell}).$
The additional variables $\xi_j$ represent the coefficients in the extended approximation space.
Using this systematic approach, one can build a hierarchy of models with increasing accuracy, each of which is a well-posed model. At the top of the hierarchy is the continuum model, represented on a finite element mesh. Then, on the same mesh, we obtain higher order approximations of MD. In the limit $\ell \to N$, the full MD description is recovered.

*Supported by National Natural Science Foundation grant DMS1016582.

To cite this abstract, use the following reference: http://meetings.aps.org/link/BAPS.2013.MAR.G11.4