DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity

Abstract

Here, we compare time averaging and ensemble averaging as two different methods for phase space sampling in molecular dynamics (MD) calculations of thermal conductivity. For the comparison, we calculate thermal conductivities of solid argon and silicon structures, using equilibrium MD. Moreover, we introduce two different schemes for the ensemble averaging approach and show that both can reduce the total simulation time as compared to time averaging. It is also found that velocity rescaling is an efficient mechanism for phase space exploration. Although our methodology is tested using classical MD, the approaches used for generating independent trajectories may find their greatest utility in computationally expensive simulations such as first principles MD. For such simulations, where each time step is costly, time averaging can require long simulation times because each time step must be evaluated sequentially and therefore phase space averaging is achieved through sequential operations. Conversely, with ensemble averaging, phase space sampling can be achieved through parallel operations, since each trajectory is independent. For this reason, particularly when using massively parallel architectures, ensemble averaging can result in much shorter simulation times (similar to 100-200X), but exhibits similar overall computational effort.

Authors:
 [1];  [2];  [1]
  1. Georgia Inst. of Technology, Atlanta, GA (United States)
  2. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States); Energy Frontier Research Centers (EFRC) (United States). Solid-State Solar-Thermal Energy Conversion Center (S3TEC)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1185896
Alternate Identifier(s):
OSTI ID: 1286997
Grant/Contract Number:  
AC05-00OR22725
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Applied Physics
Additional Journal Information:
Journal Volume: 117; Journal Issue: 04; Journal ID: ISSN 0021-8979
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
74 ATOMIC AND MOLECULAR PHYSICS; 75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY

Citation Formats

Gordiz, Kiarash, Singh, David J., and Henry, Asegun. Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity. United States: N. p., 2015. Web. doi:10.1063/1.4906957.
Gordiz, Kiarash, Singh, David J., & Henry, Asegun. Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity. United States. https://doi.org/10.1063/1.4906957
Gordiz, Kiarash, Singh, David J., and Henry, Asegun. Thu . "Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity". United States. https://doi.org/10.1063/1.4906957. https://www.osti.gov/servlets/purl/1185896.
@article{osti_1185896,
title = {Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity},
author = {Gordiz, Kiarash and Singh, David J. and Henry, Asegun},
abstractNote = {Here, we compare time averaging and ensemble averaging as two different methods for phase space sampling in molecular dynamics (MD) calculations of thermal conductivity. For the comparison, we calculate thermal conductivities of solid argon and silicon structures, using equilibrium MD. Moreover, we introduce two different schemes for the ensemble averaging approach and show that both can reduce the total simulation time as compared to time averaging. It is also found that velocity rescaling is an efficient mechanism for phase space exploration. Although our methodology is tested using classical MD, the approaches used for generating independent trajectories may find their greatest utility in computationally expensive simulations such as first principles MD. For such simulations, where each time step is costly, time averaging can require long simulation times because each time step must be evaluated sequentially and therefore phase space averaging is achieved through sequential operations. Conversely, with ensemble averaging, phase space sampling can be achieved through parallel operations, since each trajectory is independent. For this reason, particularly when using massively parallel architectures, ensemble averaging can result in much shorter simulation times (similar to 100-200X), but exhibits similar overall computational effort.},
doi = {10.1063/1.4906957},
journal = {Journal of Applied Physics},
number = 04,
volume = 117,
place = {United States},
year = {Thu Jan 29 00:00:00 EST 2015},
month = {Thu Jan 29 00:00:00 EST 2015}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 17 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Thoroughly sampling sequence space: Large-scale protein design of structural ensembles
journal, December 2002

  • Larson, Stefan M.; England, Jeremy L.; Desjarlais, John R.
  • Protein Science, Vol. 11, Issue 12
  • DOI: 10.1110/ps.0203902

Investigation of tensile response and thermal conductivity of boron-nitride nanosheets using molecular dynamics simulations
journal, June 2012

  • Mortazavi, Bohayra; Rémond, Yves
  • Physica E: Low-dimensional Systems and Nanostructures, Vol. 44, Issue 9
  • DOI: 10.1016/j.physe.2012.05.007

Climate model dependence and the replicate Earth paradigm
journal, December 2012


Thermal conductivity of diamond and related materials from molecular dynamics simulations
journal, October 2000

  • Che, Jianwei; Çağın, Tahir; Deng, Weiqiao
  • The Journal of Chemical Physics, Vol. 113, Issue 16
  • DOI: 10.1063/1.1310223

Modeling solid-state chemistry: Interatomic potentials for multicomponent systems
journal, March 1989


Thermal conductivity of solid argon from molecular dynamics simulations
journal, February 2004

  • Tretiakov, Konstantin V.; Scandolo, Sandro
  • The Journal of Chemical Physics, Vol. 120, Issue 8
  • DOI: 10.1063/1.1642611

Numerical comparisons of three recently proposed algorithms in the protein folding problem
journal, May 1997


Molecular Dynamics Study of Solid Thin-Film Thermal Conductivity
journal, March 2000

  • Lukes, J. R.; Li, D. Y.; Liang, X. -G.
  • Journal of Heat Transfer, Vol. 122, Issue 3
  • DOI: 10.1115/1.1288405

How to improve the accuracy of equilibrium molecular dynamics for computation of thermal conductivity?
journal, May 2010


Thermal rectification in pristine-hydrogenated carbon nanotube junction: A molecular dynamics study
journal, April 2014

  • Gordiz, Kiarash; Mehdi Vaez Allaei, S.
  • Journal of Applied Physics, Vol. 115, Issue 16
  • DOI: 10.1063/1.4873124

Molecular dynamics simulation of bulk silicon under strain
journal, June 2008


Fast Parallel Algorithms for Short-Range Molecular Dynamics
journal, March 1995


Thermal Conductivity of Solid Argon by Classical Molecular Dynamics
journal, January 1998


First-principles molecular dynamics simulation of water dissociation on TiO2 (110)
journal, October 1996


Thermal Vibrations of Atoms in Cubic Crystals II: The Amplitude of Atomic Vibrations
journal, November 1955


Thermal Conductivity of MgO Periclase from Equilibrium First Principles Molecular Dynamics
journal, September 2009


A B I NITIO M OLECULAR D YNAMICS WITH D ENSITY F UNCTIONAL T HEORY
journal, October 2002


Molecular dynamics study of the lattice thermal conductivity of Kr/Ar superlattice nanowires
journal, June 2004


Molecular Dynamics Simulation of Thermal Conduction in Nanoporous thin Films
journal, January 2004


Thermal interface conductance in Si/Ge superlattices by equilibrium molecular dynamics
journal, May 2012


Molecular dynamics simulations for the prediction of thermal conductivity of bulk silicon and silicon nanowires: Influence of interatomic potentials and boundary conditions
journal, August 2011

  • Abs da Cruz, Carolina; Termentzidis, Konstantinos; Chantrenne, Patrice
  • Journal of Applied Physics, Vol. 110, Issue 3
  • DOI: 10.1063/1.3615826

Ab Initio Molecular Dynamics Simulations for Collision between C 60 and Alkali-Metal Ions: A Possibility of Li@ C 60
journal, May 1996


Spectral Phonon Transport Properties of Silicon Based on Molecular Dynamics Simulations and Lattice Dynamics
journal, February 2008

  • Henry, Asegun S.; Chen, Gang
  • Journal of Computational and Theoretical Nanoscience, Vol. 5, Issue 2
  • DOI: 10.1166/jctn.2008.2454

Thermal Conductivity of Individual Single-Wall Carbon Nanotubes
journal, September 2006

  • Lukes, Jennifer R.; Zhong, Hongliang
  • Journal of Heat Transfer, Vol. 129, Issue 6
  • DOI: 10.1115/1.2717242

Accelerated molecular dynamics force evaluation on graphics processing units for thermal conductivity calculations
journal, May 2013


Molecular dynamics simulation of thermal conductivity of silicon nanowires
journal, October 1999

  • Volz, Sebastian G.; Chen, Gang
  • Applied Physics Letters, Vol. 75, Issue 14
  • DOI: 10.1063/1.124914

Thermal conductivity decomposition and analysis using molecular dynamics simulations. Part I. Lennard-Jones argon
journal, April 2004


Molecular-dynamics simulation of thermal conductivity of silicon crystals
journal, January 2000


GPU accelerated molecular dynamics simulation of thermal conductivities
journal, February 2007

  • Yang, Juekuan; Wang, Yujuan; Chen, Yunfei
  • Journal of Computational Physics, Vol. 221, Issue 2
  • DOI: 10.1016/j.jcp.2006.06.039

A sampling problem in molecular dynamics simulations of macromolecules.
journal, April 1995

  • Clarage, J. B.; Romo, T.; Andrews, B. K.
  • Proceedings of the National Academy of Sciences, Vol. 92, Issue 8
  • DOI: 10.1073/pnas.92.8.3288

Exploring phase space
journal, May 2000


Thermal conductivity decomposition and analysis using molecular dynamics simulations
journal, April 2004


High Thermal Conductivity of Single Polyethylene Chains Using Molecular Dynamics Simulations
journal, December 2008


Charge-polarized interfacial superlattices in marginally twisted hexagonal boron nitride
journal, January 2021


Works referencing / citing this record:

Phonon Transport at Crystalline Si/Ge Interfaces: The Role of Interfacial Modes of Vibration
journal, March 2016

  • Gordiz, Kiarash; Henry, Asegun
  • Scientific Reports, Vol. 6, Issue 1
  • DOI: 10.1038/srep23139

Phonon transport at interfaces: Determining the correct modes of vibration
journal, January 2016

  • Gordiz, Kiarash; Henry, Asegun
  • Journal of Applied Physics, Vol. 119, Issue 1
  • DOI: 10.1063/1.4939207

Interface conductance modal analysis of lattice matched InGaAs/InP
journal, May 2016

  • Gordiz, Kiarash; Henry, Asegun
  • Applied Physics Letters, Vol. 108, Issue 18
  • DOI: 10.1063/1.4948520

Phonon transport at interfaces between different phases of silicon and germanium
journal, January 2017

  • Gordiz, Kiarash; Henry, Asegun
  • Journal of Applied Physics, Vol. 121, Issue 2
  • DOI: 10.1063/1.4973573

Using Green-Kubo modal analysis (GKMA) and interface conductance modal analysis (ICMA) to study phonon transport with molecular dynamics
journal, February 2019

  • Seyf, Hamid Reza; Gordiz, Kiarash; DeAngelis, Freddy
  • Journal of Applied Physics, Vol. 125, Issue 8
  • DOI: 10.1063/1.5081722

Interface conductance modal analysis of a crystalline Si-amorphous SiO 2 interface
journal, April 2019

  • Gordiz, Kiarash; Muraleedharan, Murali Gopal; Henry, Asegun
  • Journal of Applied Physics, Vol. 125, Issue 13
  • DOI: 10.1063/1.5085328

A formalism for calculating the modal contributions to thermal interface conductance
journal, September 2015