Thermal properties of molybdenum disulfide $\left({\mathrm{MoS}}_2\right)$ have recently attracted attention related to fundamentals of heat propagation in strongly anisotropic materials, and in the context of potential applications to optoelectronics and thermoelectrics. Multiple empirical potentials have been developed for classical molecular dynamics (MD) simulations of this material, but it has been unclear which provides the most realistic results. Here, we calculate lattice thermal conductivity of single- and multilayer pristine ${\mathrm{MoS}}_2$ by employing three different thermal transport MD methods: equilibrium, nonequilibrium, and homogeneous nonequilibrium ones. We mainly use the Graphics Processing Units Molecular Dynamics code for numerical calculations, and the Large-scale Atomic/Molecular Massively Parallel Simulator code for crosschecks. Using different methods and computer codes allows us to verify the consistency of our results and facilitate comparisons with previous studies, where different schemes have been adopted. Our results using variants of the Stillinger-Weber potential are at odds with some previous ones and we analyze the possible origins of the discrepancies in detail. We show that, among the potentials considered here, the reactive empirical bond order (REBO) potential gives the most reasonable predictions of thermal transport properties as compared to experimental data. With the REBO potential, we further find that isotope scattering has only a small effect on thermal conduction in ${\mathrm{MoS}}_2$ and the in-plane thermal conductivity decreases with increasing layer number and saturates beyond about three layers. We identify the REBO potential as a transferable empirical potential for MD simulations of ${\mathrm{MoS}}_2$ which can be used to study thermal transport properties in more complicated situations such as in systems containing defects or engineered nanoscale features. This work establishes a firm foundation for understanding heat transport properties of ${\mathrm{MoS}}_2$ using MD simulations.

The thermal transfer across Pd-H₂O interface is dependent on water structure, where the partial ordered water structure rises the interfacial thermal conductance.

Grain boundaries (GBs) are inevitable defects in large-area MoS₂ samples, and they greatly affect the thermal transport of MoS₂.

Understanding the structural transition of ionic liquids (ILs) confined in a nanospace is imperative for the application of ILs in energy storage, gas separation, and other chemical engineering techniques. In this work, the quantitative relations between the properties and height of the nanochannel (H) for the ([Emim]⁺[TF₂N]⁻) IL are explored through molecular dynamics simulations. Interestingly, the entropy of the confined IL exhibits a nonmonotonic behavior as H increases: initially increasing for H < 1.0 nm and then decreasing for 1.0 < H < 1.1 nm, followed by increasing again for H > 1.1 nm; it finally approaches that of liquid bulk ILs. The vibrational spectrum of the confined IL is analyzed to investigate the nature of nonmonotonic entropy, showing that the liquidity and partial solidity will be respectively attenuated and enhanced as H decreases from 5.0 to 0.75 nm. Moreover, the hydrogen bond (HB) network and external force are also calculated, showing similar nonmonotonic behaviors when compared with the thermodynamic properties. The entropy gain of the confined IL originates from the reduced HB interactions, weaker external force, and partial solid nature, where more phase space sampling for ILs inside a bilayer graphene nanochannel (BLGC) can be achieved. All the above relations demonstrate that there exists a critical height of the nanochannel (H_CR = 1.0 nm) at which the confined IL possesses weaker HB interaction, higher entropy, and better stability. The critical height of the nanochannel is also identified in the analysis of the local structures of cation head groups and anions, indicating that the confined IL could have a faster in-plane diffusive ability. These factors can serve as key indicators in quantitatively characterizing the mechanism for the structural transition of ILs inside a nanochannel and facilitate the rational design of nanopores and nanochannels to regulate the properties and structures of ILs in practical application scenarios.

In this paper, a new nanostructure is proposed, namely, the knitted graphene-nanoribbon sheet (KGS), which consists of zigzag and/or armchair graphene nanoribbons. The knitting technology is introduced to graphene nanotechnology to produce large area graphene sheets. Compared with pristine graphene, the chirality of a knitted graphene-nanoribbon sheet is much more flexible and can be designed on demand. The mechanical properties of KGSs are investigated by molecular dynamics simulations, including the effect of vacancies. With hydrogen atoms saturating the ribbon edges, the structure (KGS + H) is found to be of significant mechanical robustness, whose fracture does not rely on the critical bonds. The fracture strain of KGS + H remains nearly unchanged as long as there remains a single defect-free graphene nanoribbon in the tensile direction. This graphene nano knitting technique is experimentally feasible, inspired by a recent demonstration by Fournier et al. [Phys. Rev. B, 2011, 84, 035435] of lifting a single molecular wire using a combined frequency-modulated atomic force and tunnelling microscope.