Thermal conductivity of defective graphene: an efficient molecular dynamics study based on graphics processing units

Thermal conductivity of wrinkled graphene ring with defects

Graphene rings have great prospects in the fields of biological modulators, electrochemical biosensors, and resonators, but are prone to wrinkling which can affect their physical properties. This work establishes a theoretical model predicting the torsional wrinkling behavior of defective monolayer graphene rings, which provides direct understanding and reliable accuracy of the wrinkle levels. Then the thermal conductivity of wrinkled graphene rings is studied considering different wrinkle levels, defect concentrations and radii. It is found that with increased radius, defect concentration and torsional angle, the ratio of wrinkle amplitude to wavelength increases gradually. Vacancy defects and radii have more significant influences on the thermal conductivity than torsional wrinkles. The main influence mechanism of wrinkles and defects on thermal conductivity is revealed by phonon density of state. This work provides theoretical guidance for thermal manipulation based on the wrinkle-tuning approach.

Thermal transport in a defective pillared graphene network: insights from equilibrium molecular dynamics simulation

Graphene-based hybrid nanostructures have great potential to be ideal candidates for developing tailored thermal transport materials. In this study, we perform equilibrium molecular dynamics simulations employing the Green-Kubo method to investigate the influence of topological defects in three-dimensional pillared graphene networks. Similar to single-layer graphene and carbon nanotubes, the thermal conductivity (k) of pillared graphene systems exhibits a strong correlation with the system size (L), following a power-law relation k ∼ Lα, where α ranges from 0.12 to 0.15. Our results indicate that the vacancy defects significantly reduce thermal conductivity in pillared graphene systems compared to Stone-Wales defects. We observe that, beyond defect concentration, the location of the defects also plays a crucial role in determining thermal conductivity. We further analyze the phonon vibrational spectrum and the phonon participation ratio to obtain more insight into the thermal transport in the defective pillared graphene network. In most scenarios, longitudinal and flexural acoustic phonons experience significant localization within the 15-45 THz frequency range in the defective pillared graphene system.

Molecular Dynamics Simulation and Experimental Study of Mechanical Properties of Graphene-Cement Composites

To investigate the mechanical properties of graphene (G) and calcium silicate hydrate (C-S-H) composites in different directions, molecular dynamics (MD) simulations and experiments were used, and the effects of temperature, loading rate, and graphene defects were also investigated. The experimental results show that the addition of graphene can improve the flexural, compressive, and tensile strength of the composite. The results of molecular dynamics simulation show that the addition of graphene in x and z directions can enhance the tensile strength of G/C-S-H in three directions, while the addition of graphene in y direction can reduce the tensile strength of G/C-S-H. At the same time, the tensile strength of G/C-S-H decreases with the increase in temperature and increases with the increase in loading rate. Meanwhile, the mechanical properties of G/C-S-H can be improved using a certain concentration of monatomic vacancy defects, diatomic vacancy defects, and S-W defects.

Molecular Dynamics Simulation on In-Plane Thermal Conductivity of Graphene/Hexagonal Boron Nitride van der Waals Heterostructures

Graphene, hexagonal boron nitride (h-BN), and their heterostructures are promising thermal interface materials due to the outstanding thermal properties of graphene and h-BN. For the heterostructures, extensive work has mainly focused on the thermal transport of two-dimensional (2D) graphene/h-BN (GBN) in-plane heterostructures in which graphene and h-BN are bonded at the interface. In this study, we investigate the thermal conductivity of three-dimensional (3D) GBN van der Waals (vdW) heterostructures by means of nonequilibrium molecular dynamics (NEMD) simulations. Unlike the 2D GBN in-plane heterostructure, the 3D GBN vdW heterostructure consists of three layers where graphene is sandwiched by two h-BN sheets via vdW forces. Various techniques, including hydrogen-functionalization, vacancy defects, tensile strain, interlayer coupling strength, layer numbers of h-BN, size effect, and temperature, are extensively explored to find an effective route for the modulation of the thermal conductivity. It is found that the thermal conductivity of the triple-layer GBN vdW heterostructure is very sensitive to these extrinsic factors. Of these, hydrogen-functionalization is the most effective method. A low hydrogen coverage of 1% in the sandwiched graphene can lead to 55% reduction in the thermal conductivity of the vdW heterostructure. Vacancy defects on graphene exert a more significant effect on the thermal conductivity reduction for the vdW heterostructure than B or N vacancies in the outer h-BN layers. This work reveals the physical mechanism for manipulating the thermal transport along the GBN vdW heterostructures via structural modification and provides a useful guideline for designing novel thermal management devices based on the GBN vdW heterostructures.

GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations

We present our latest advancements of machine-learned potentials (MLPs) based on the neuroevolution potential (NEP) framework introduced in [Fan et al., Phys. Rev. B 104, 104309 (2021)] and their implementation in the open-source package GPUMD.We increase the accuracy of NEP models both by improving the radial functions in the atomic-environment descriptor using a linear combination of Chebyshev basis functions and by extending the angular descriptor with some four-body and five-body contributions as in the atomic cluster expansion approach.We also detail our efficient implementation of the NEP approach in graphics processing units as well as our workflow for the construction of NEP models, and we demonstrate their application in large-scale atomistic simulations.By comparing to state-of-the-art MLPs, we show that the NEP approach not only achieves above-average accuracy but also is far more computationally efficient.These results demonstrate that the GPUMD package is a promising tool for solving challenging problems requiring highly accurate, large-scale atomistic simulations.To enable the construction of MLPs using a minimal training set, we propose an active-learning scheme based on the latent space of a pre-trained NEP model.Finally, we introduce three separate Python packages, GPYUMD, CALORINE, and PYNEP, which enable the integration of GPUMD into Python workflows.

Strong reduction of thermal conductivity of WSe2with introduction of atomic defects

The thermal conductivities of pristine and defective single-layer tungsten diselenide (WSe₂) are investigated by using equilibrium molecular dynamics method. The thermal conductivity of WSe₂increases dramatically with size below a characteristic of ~5 nm and levels off for broader samples and reaches a constant value of ~2 W/mK. By introducing atomic vacancies, we discovered that the thermal conductivity of WSe₂is significantly reduced. In particular, the W vacancy has a greater impact on thermal conductivity reduction than Se vacancies: the thermal conductivity of pristine WSe₂is reduced by ~60% and ~70% with the adding of ~1% of Se and W vacancies, respectively. The reduction of thermal conductivity is found to be related to the decrease of mean free path (MFP) of phonons in the defective WSe₂. The MFP of WSe₂decreases from ~4.2 nm for perfect WSe₂to ~2.2 nm with the addition of 0.9% Se vacancies. More sophisticated types of point defects, such as vacancy clusters and anti-site defects, are explored in addition to single vacancies and are found to dramatically renormalize the phonons. The reconstruction of the bonds leads to localized phonons in the forbidden gap in the phonon density of states which leads to a drop in thermal conduction. This work demonstrates the influence of different defects on the thermal conductivity of single-layer WSe₂, providing insight into the process of defect-induced phonon transport as well as ways to improve heat dissipation in WSe₂-based electronic devices.

Anomalous thermal conductivity enhancement in low dimensional resonant nanostructures due to imperfections

Nanophononic metamaterials have broad applications in fields such as heat management, thermoelectric energy conversion, and nanoelectronics. Phonon resonance in pillared low-dimensional structures has been suggested to be a feasible approach to reduce thermal conductivity (TC). In this work, we study the effects of imperfections in pillared nanostructures based on graphene nanoribbons (GNR), using classical molecular dynamics simulations and harmonic lattice dynamics. The TC of perfect pillared GNR is only about 13% of that of pristine GNR due to the strong phonon resonant hybridization in pillared GNR. However, introducing imperfections such as vacancy defects and mass mismatch between the pillars and the base material, and alloy disorder in the pillars, can weaken the resonant hybridization and abnormally increase the TC. We show that both vacancy defects and mass mismatch can reduce the penetration of the resonant modes from the pillars into the base material, while the alloy disorder in the pillars can scatter the phonons inside them, which turns regular resonance into a random one with weaker hybridization. Our work provides useful insight into the phonon resonance mechanisms in experimentally relevant low dimensional nanostructures containing various imperfections.

Anisotropic lattice thermal conductivity in topological semimetal ZrGeX(XS, Se, Te): a first-principles study

Topological semimetals have attracted significant attentions owing to their potential applications in numerous fields such as low-power electron devices and quantum computation, which are closely related to their thermal transport properties. In this work, the phonon transport properties of topological Dirac nodal-line semimetals ZrGeX(X= S, Se, Te) with the PbClF-type structures are systematically studied using the first-principles calculations combined with the Boltzmann transport theory. The obtained lattice thermal conductivities show an obvious anisotropy, which is caused by the layer structures of ZrGeX(X= S, Se, Te). The room-temperature lattice conductivity of ZrGeTe alongcdirection is found to be as low as 0.24 W m^-1 K^-1, indicating that it could be of great significance in the fields of thermal coating materials and solar cell absorber. In addition, we extract each phonon branch from group velocities, phonon scattering rates, Grüneisen parameters, and phase space volumes to investigate the mechanism underlying the low thermal conductivity. It is concluded that the difference of thermal conductivities of three materials may be caused by the number of scattering channels and the effect of anharmonic. Furthermore, the phonon mean free path alongadirection is relatively longer. Nanostructures or polycrystalline structures may be effective to reduce the thermal conductivity and improve the thermoelectric properties.

Why thermal conductivity of CaO is lower than that of CaS: a study from the perspective of phonon splitting of optical mode

Generally speaking, for materials with the same structure, the thermal conductivity is higher for lighter atomic masses. However, we found that the thermal conductivity of CaO is lower than that of CaS, despite the lighter atomic mass of O than S. To uncover the underlying physical mechanisms, the thermal conductivity of CaM (M = O, S, Se, Te) and the corresponding response to strain is investigated by performing first-principles calculations along with the phonon Boltzmann transport equation. For unstrained system, the order of thermal conductivity is CaS > CaO > CaSe > CaTe. This order remains unchanged in the strain range of -2% to 5%. When the compressive strain is larger than 2%, the thermal conductivity of CaO surpasses that of CaS and becomes the highest thermal conductivity material among the four compounds. By analyzing the mode-dependent phonon properties, the phonon lifetime is found to be dominant over other influential factors and leads to the disparate response of thermal conductivity under strain. Moreover, the changing trend of three-phonon scattering phase space is consistent with that of phonon lifetime, which is directly correlated to the phonon frequency gap induced by the LO-TO splitting. The variation of Born effective charge is found to be opposite for CaM. The Born effective charge of CaO decreases with tensile strain increasing, demonstrating stronger charge delocalization and lower ionicity, while the Born effective charges of CaS, CaSe, and CaTe show a dramatic increase. Such variation indicates that the bonding nature can be effectively tuned by external strain, thus affecting the phonon anharmonic properties and thermal conductivity. The difference of bonding nature is further confirmed by the band structure. Our results show that the bonding nature of CaM can be modulated by external strain and leads to disparate strain dependent thermal conductivity.