# Bulletin of the American Physical Society

# 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

## Follow Us |
## Engage
Become an APS Member |
## My APS
Renew Membership |
## Information for |
## About APSThe American Physical Society (APS) is a non-profit membership organization working to advance the knowledge of physics. |

© 2018 American Physical Society
| All rights reserved | Terms of Use
| Contact Us

**Headquarters**
1 Physics Ellipse, College Park, MD 20740-3844
(301) 209-3200

**Editorial Office**
1 Research Road, Ridge, NY 11961-2701
(631) 591-4000

**Office of Public Affairs**
529 14th St NW, Suite 1050, Washington, D.C. 20045-2001
(202) 662-8700