A Caltech Library Service

Finite-temperature extension of the quasicontinuum method using Langevin dynamics: entropy losses and analysis of errors

Marian, J. and Venturini, G. and Hansen, B. L. and Knapp, J. and Ortiz, M. and Campbell, G. H. (2010) Finite-temperature extension of the quasicontinuum method using Langevin dynamics: entropy losses and analysis of errors. Modelling and Simulation in Materials Science and Engineering, 18 (1). Art. No. 015003. ISSN 0965-0393.

Full text is not posted in this repository. Consult Related URLs below.

Use this Persistent URL to link to this item:


The concurrent bridging of molecular dynamics and continuum thermodynamics presents a number of challenges, mostly associated with energy transmission and changes in the constitutive description of a material across domain boundaries. In this paper, we propose a framework for simulating coarse dynamic systems in the canonical ensemble using the quasicontinuum method (QC). The equations of motion are expressed in reduced QC coordinates and are strictly derived from dissipative Lagrangian mechanics. The derivation naturally leads to a classical Langevin implementation where the timescale is governed by vibrations emanating from the finest length scale occurring in the computational cell. The equations of motion are integrated explicitly via Newmark's (β = 0; γ = 1/2) method, which is parametrized to ensure overdamped dynamics. In this fashion, spurious heating due to reflected vibrations is suppressed, leading to stable canonical trajectories. To estimate the errors introduced by the QC reduction in the resulting dynamics, we have quantified the vibrational entropy losses in Al uniform meshes by calculating the thermal expansion coefficient for a number of conditions. We find that the entropic depletion introduced by coarsening varies linearly with the element size and is independent of the nodal cluster diameter. We rationalize the results in terms of the system, mesh and cluster sizes within the framework of the quasiharmonic approximation. The limitations of the method and alternatives to mitigate the errors introduced by coarsening are discussed. This work represents the first of a series of studies aimed at developing a fully non-equilibrium finite-temperature extension of QC.

Item Type:Article
Related URLs:
URLURL TypeDescription
Ortiz, M.0000-0001-5877-4824
Additional Information:© 2010 Institute of Physics and IOP Publishing Limited. Print publication: Issue 1 (January 2010); received 29 April 2009; in final form 15 September 2009; published 10 December 2009. This work performed under LDRD project Project 06-SI-005 under the auspices of the US Department of Energy by Lawrence Livermore National Laboratory under contract DE-AC52- 07NA27344. GV and MO gratefully acknowledge the support of the Department of Energy through Caltech’s PSAAP Center for the Predictive Simulation of the Dynamic Response of Materials.
Funding AgencyGrant Number
Department of Energy (DOE)DE-AC52-07NA27344
Issue or Number:1
Record Number:CaltechAUTHORS:20100121-141104452
Persistent URL:
Official Citation:Finite-temperature extension of the quasicontinuum method using Langevin dynamics: entropy losses and analysis of errors J Marian, G Venturini, B L Hansen, J Knap, M Ortiz and G H Campbell 2010 Modelling Simul. Mater. Sci. Eng. 18 015003 (31pp) doi: 10.1088/0965-0393/18/1/015003
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:17272
Deposited By: Jason Perez
Deposited On:22 Jan 2010 21:54
Last Modified:24 Nov 2020 00:47

Repository Staff Only: item control page