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. 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 contributemore »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.« less
  2. Magnetic excitations in van der Waals (vdW) materials, especially in the two-dimensional (2D) limit, are an exciting research topic from both the fundamental and applied perspectives. Using temperature-dependent, magneto-Raman spectroscopy, we identify the hybridization of two-magnon excitations with two phonons in manganese phosphorus triselenide (MnPSe 3 ), a magnetic vdW material that hosts in-plane antiferromagnetism. Results from first-principles calculations of the phonon and magnon spectra further support our identification. The Raman spectra’s rich temperature dependence through the magnetic transition displays an avoided crossing behavior in the phonons’ frequency and a concurrent decrease in their lifetimes. We construct a model basedmore »on the interaction between a discrete level and a continuum that reproduces these observations. Our results imply a strong hybridization between each phonon and a two-magnon continuum. This work demonstrates that the magnon-phonon interactions can be observed directly in Raman scattering and provides deep insight into these interactions in 2D magnetic materials.« less
  3. 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 (algebraicmore »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 of 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
  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 themore »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.« less
  5. 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 themore »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.« less