skip to main content

Title: Phonon spectrum and phonon focusing in coarse-grained atomistic simulations
This work investigates the accuracy, efficiency, and applicability of coarse-grained (CG) atomistic methods in simulation of phonon dynamics. First, we compute and compare phonon dispersion relations in CG models with those in atomically resolved models, using the concurrent atomistic-continuum (CAC). The CG atomistic models using the CAC method are shown to reproduce long-wavelength phonons with great accuracy, while capturing the dynamics of some short-wavelength phonons that are usually inaccessible to CG methods. We then present CG simulation results of the propagation of heat pulses in Si with the interaction between atoms being modelled with the Stillinger-Weber potential; the experimentally observed phonon-focusing patterns in the (1 0 0) and (1 1 1) planes of Si crystals are reproduced. The accuracy and efficiency of the CAC method in CG simulation of acoustic and optical phonon branches are quantified with respect to atomically-resolved molecular dynamics simulations. The applicability and limitations of concurrent multiscale methods in the simulation of phonon transport across atomistic-continuum interface are investigated. Possible ways to overcome the limitations are discussed.
Authors:
; ; ; ; ;
Award ID(s):
1761512
Publication Date:
NSF-PAR ID:
10093449
Journal Name:
Computational materials science
Volume:
162
ISSN:
0927-0256
Sponsoring Org:
National Science Foundation
More Like this
  1. Recent developments in generalized continuum modeling methods ranging from coarse-grained atomistics to micromorphic theory offer potential to make more intimate physical contact with dislocation field problems framed at length scales on the order of microns. We explore a range of discrete dynamical and continuum mechanics approaches to crystal plasticity that are relevant to modeling behavior of populations of dislocations. Predictive atomistic and coarse-grained atomistic models are limited in terms of length and time scales that can be accessed; examples of the latter are discussed in terms of interactions of multiple dislocations in heterogeneous systems. Generalized continuum models alleviate restrictions to a significant extent in modeling larger scales of dislocation configurations and reactions, and are useful to consider effects of dislocation configuration on strength at characteristic length scales of sub-micron and above; these models require a combination of bottomup models and top-down experimental information to inform parameters and model form. The concurrent atomistic-continuum (CAC) method is extended to model complex multicomponent alloy systems using an average atom approach. Examples of CAC are presented, along with potential to assist in informing parameters of a recently developed micropolar crystal plasticity model based on a set of sub-micron dislocation field problems. Prospects for furthermore »developments are discussed.« less
  2. We developed a new coarse-grained (CG) molecular dynamics force field for polyacrylamide (PAM) polymer based on fitting to the quantum mechanics (QM) equation of state (EOS). In this method, all nonbond interactions between representative beads are parameterized using a series of QM-EOS, which significantly improves the accuracy in comparison to common CG methods derived from atomistic molecular dynamics. This CG force-field has both higher accuracy and improved computational efficiency with respect to the OPLS atomistic force field. The nonbond components of the EOS were obtained from cold-compression curves on PAM crystals with rigid chains, while the covalent terms that contribute to the EOS were obtained using relaxed chains. For describing PAM gels we developed water–PAM interaction parameters using the same method. We demonstrate that the new CG-PAM force field reproduces the EOS of PAM crystals, isolated PAM chains, and water–PAM systems, while successfully predicting such experimental quantities as density, specific heat capacity, thermal conductivity and melting point.
  3. In this work, we investigate misfit dislocations in PbTe/PbSe heteroepitaxial systems using the concurrent atomistic–continuum (CAC) method. A potential model containing the long-range Coulombic interaction and short-range Buckingham potential is developed for the system. By considering the minimum potential energy of relaxed interface structures for various initial conditions and PbTe layer thicknesses, the equilibrium structure of misfit dislocations and the dislocation spacings in PbTe/PbSe(001) heteroepitaxial thin films are obtained as a function of the PbTe layer thicknesses grown on a PbSe substrate. The critical layer thickness above which misfit dislocations inevitably form, the structure of the misfit dislocations at the interfaces, and the dependence of average dislocation spacing on PbTe layer thickness are obtained and discussed. The simulation results provide an explanation for the narrowing of the spread of the distribution of misfit dislocation spacing as layer thickness increases in PbTe/PbSe(001) heteroepitaxy.
  4. In this work, we investigate misfit dislocations in PbTe/PbSe heteroepitaxial systems using the concurrent atomistic–continuum (CAC) method. A potential model containing the long-range Coulombic interaction and short-range Buckingham potential is developed for the system. By considering the minimum potential energy of relaxed interface structures for various initial conditions and PbTe layer thicknesses, the equilibrium structure of misfit dislocations and the dislocation spacings in PbTe/PbSe(001) heteroepitaxial thin films are obtained as a function of the PbTe layer thicknesses grown on a PbSe substrate. The critical layer thickness above which misfit dislocations inevitably form, the structure of the misfit dislocations at the interfaces, and the dependence of average dislocation spacing on PbTe layer thickness are obtained and discussed. The simulation results provide an explanation for the narrowing of the spread of the distribution of misfit dislocation spacing as layer thickness increases in PbTe/PbSe(001) heteroepitaxy.
  5. We present data-driven coarse-grained (CG) modeling for polymers in solution, which conserves the dynamic as well as structural properties of the underlying atomistic system. The CG modeling is built upon the framework of the generalized Langevin equation (GLE). The key is to determine each term in the GLE by directly linking it to atomistic data. In particular, we propose a two-stage Gaussian process-based Bayesian optimization method to infer the non-Markovian memory kernel from the data of the velocity autocorrelation function (VACF). Considering that the long-time behaviors of the VACF and memory kernel for polymer solutions can exhibit hydrodynamic scaling (algebraic decay with time), we further develop an active learning method to determine the emergence of hydrodynamic scaling, which can accelerate the inference process of the memory kernel. The proposed methods do not rely on how the mean force or CG potential in the GLE is constructed. Thus, we also compare two methods for constructing the CG potential: a deep learning method and the iterative Boltzmann inversion method. With the memory kernel and CG potential determined, the GLE is mapped onto an extended Markovian process to circumvent the expensive cost of directly solving the GLE. The accuracy and computational efficiency ofmore »the proposed CG modeling are assessed in a model star-polymer solution system at three representative concentrations. By comparing with the reference atomistic simulation results, we demonstrate that the proposed CG modeling can robustly and accurately reproduce the dynamic and structural properties of polymers in solution.« less