Researcher profile

Jin-Wu Jiang

Jin-Wu Jiang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
66works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

66 published item(s)

preprint2021arXiv

Misfit strain effect on the thermal expansion coefficient of graphene/MoS$_2$ van der Waals heterostructures

Van der Waals heterostructures such as graphene/MoS$_2$ are promising candidates for plenty of optical or electronic applications, owing to advanced properties inherited from the constitutional atomic layers. Thermal expansion is an important phenomenon to be considered for the thermal stability of the van der Waals heterstructure as temperature commonly rises during the operation of nano devices. In the present work, the thermal expansion coefficient for the graphene/MoS$_2$ heterostructure is investigated by molecular dynamics simulations, and the effect from the unavoidable misfit strain on the thermal expansion coefficient is revealed. The misfit strain can tune the thermal expansion coefficient by a factor of six, and this effect is quite robust in sense that it is not sensitive to the size or direction of the heterostructure. An analytic formula is derived to directly relate the thermal expansion coefficient to the misfit strain of the heterostructure, which qualitatively agrees with the numerical results although the analytic formula underestimates the misfit strain effect. Further analysis discloses that the misfit strain can efficiently engineer the thermal induced ripples, which serves as the key mechanism for the misfit strain effect on the thermal expansion coefficient. These findings provide valuable information for the thermal stability of van der Waals heterostructures and shall be benefit for practical applications of van der Waals heterostructure based nano devices.

preprint2020arXiv

A physical description for the monoclinic phase of zirconia based on the Tersoff potential

Zirconia is well-known for plenty of important morphologys with Zr coordination varying from sixfold in the octagonal phase to eightfold in the cubic or tetragonal phase. The development of empirical potentials to describe these zirconia morphologys is an important issue but a long-standing challenge, which becomes a bottleneck for the theoretical investigation of large zirconia structures. In contrast to the standard core-shell model, we develop a new potential for zirconia through the combination of long-range Coulomb interaction and bond order Tersoff model. The bond order characteristic of the Tersoff potential enables it to be well suited for the description of these zirconia morphologys with different coordination numbers. In particular, the complex monoclinic phase with two inequivalent oxygens, that is difficult to be described by most existing empirical potentials, can be well captured by this newly developed potential. It is shown that this potential can provide reasonable predictions for most static and dynamic properties of various zirconia morphologys. Besides its clear physical essence, this potential is at least one order faster than core-shell based potentials in the molecule dynamics simulation, as it discards the concept of the ultralight shell that demands for an extremely small time step. We also provide potential scripts for the widely used packages GULP and LAMMPS.

preprint2017arXiv

Molecular Dynamics Simulations for Anisotropic Thermal Conductivity of Borophene

The present work carries out molecular dynamics simulations to compute the thermal conductivity of the borophene nanoribbon and the borophene nanotube using the Muller-Plathe approach. We investigate the thermal conductivity of the armchair and zigzag borophenes, and show the strong anisotropic thermal conductivity property of borophene. We compare the results of the borophene nanoribbon and the borophene nanotube, and find the thermal conductivity of the borophene is structure dependent.

preprint2016arXiv

Auxetic Nanomaterials: Recent Progress and Future Development

Auxetic materials (materials with negative Poisson's ratio) and nanomaterials have independently been for many years two of the most active research fields in material science. Recently, these formerly independent fields have begun to intersect in new and interesting ways due to the recent discovery of auxeticity in nanomaterials like graphene, metal nanoplates, black phosphorus, and others. Here we review the research emerging at the intersection of auxeticity and nanomaterials. We first survey the atomistic mechanisms, both intrinsic and extrinsic, that have been found, primarily through atomistic simulations, to cause auxeticity in nanomaterials. We then outline the available experimental evidence for auxetic nanomaterials. In order to lay the groundwork for future work in this exciting area, we close by discussing several future prospects as well as the current challenges in this field.

preprint2016arXiv

Buckled Graphene for Efficient Energy Harvest, Storage, and Conversion

Buckling is one of the most common phenomena in atomic-thick layered structures like graphene. While the buckling phenomenon usually causes disaster for most nano-devices, we illustrate one positive application of the buckled graphene for energy harvest, storage, and conversion. More specifically, we perform molecular dynamical simulations to show that the buckled graphene can be used to collect the wasted mechanical energy and store the energy in the form of internal knotting potential. Through strain engineering, the knotting potential can be converted into useful kinetic (thermal) energy that is highly concentrated at the free edges of the buckled graphene. The present study demonstrates potential applications of the buckled graphene for converting the dispersed wasted mechanical energy into the concentrated useful kinetic (thermal) energy.

preprint2016arXiv

Intrinsic Negative Poisson's Ratio for Single-Layer Graphene

Negative Poisson's ratio (NPR) materials have drawn significant interest because the enhanced toughness, shear resistance and vibration absorption that typically are seen in auxetic materials may enable a range of novel applications. In this work, we report that single-layer graphene exhibits an intrinsic NPR, which is robust and independent of its size and temperature. The NPR arises due to the interplay between two intrinsic deformation pathways (one with positive Poisson's ratio, the other with NPR), which correspond to the bond stretching and angle bending interactions in graphene. We propose an energy-based deformation pathway criteria, which predicts that the pathway with NPR has lower energy and thus becomes the dominant deformation mode when graphene is stretched by a strain above 6%, resulting in the NPR phenomenon.

preprint2016arXiv

Negative Poisson's Ratio in Single-Layer Graphene Ribbons

The Poisson's ratio characterizes the resultant strain in the lateral direction for a material under longitudinal deformation. Though negative Poisson's ratios (NPR) are theoretically possible within continuum elasticity, they are most frequently observed in engineered materials and structures, as they are not intrinsic to many materials. In this work, we report NPR in single-layer graphene ribbons, which results from the compressive edge stress induced warping of the edges. The effect is robust, as the NPR is observed for graphene ribbons with widths smaller than about 10 nm, and for tensile strains smaller than about 0.5%, with NPR values reaching as large as -1.51. The NPR is explained analytically using an inclined plate model, which is able to predict the Poisson's ratio for graphene sheets of arbitrary size. The inclined plate model demonstrates that the NPR is governed by the interplay between the width (a bulk property), and the warping amplitude of the edge (an edge property), which eventually yields a phase diagram determining the sign of the Poisson's ratio as a function of the graphene geometry.

preprint2016arXiv

Origin of Hinge-Like Mechanism in Single-Layer Black Phosphorus: the Angle-Angle Cross Interaction

The single-layer black phosphorus is characteristic for its puckered configuration that possesses the hinge-like mechanism, which leads to the highly anisotropic in-plane Poisson's ratios and the negative out-of-plane Poisson's ratio. We reveal that the hinge-like mechanism can be attributed to the angle-angle cross interaction, which, combined with the bond stretching and angle bending interactions, is able to provide a good description of the mechanical properties in the single-layer black phosphorus. We also propose a nonlinear angle-angle cross interaction, which follows the form of the Stillinger-Weber potential and is advantageous for molecular dynamics simulations of single-layer black phosphorous under large deformations.

preprint2016arXiv

Tunable negative Poisson's ratio in hydrogenated graphene

We perform molecular dynamics simulations to investigate the effect of hydrogenation on the Poisson's ratio of graphene. It is found that the value for Poisson's ratio in graphene can be effectively tuned from positive to negative by varying the percentage of hydrogenation. Specifically, the Poisson's ratio decreases with the increase of the percentage of hydrogenation, and reaches a minimum value of -0.04 at the percentage of hydrogenation about 50%. The Poisson's ratio starts to increase with a further increase of the percentage of hydrogenation. The appearance of minimum negative Poisson's ratio in the hydrogenated graphene is attributed to the suppression of the hydrogenation-induced ripples during stretching of graphene. Our results demonstrate that hydrogenation is a valuable approach for the tuning of the Poisson's ratio from positive to negative in graphene.

preprint2015arXiv

A Gaussian Treatment for the Friction Issue of Lennard-Jones Potential in Layered Materials: Application to Friction between Graphene, MoS2 and Black Phosphorus

The Lennard-Jones potential is widely used to describe the interlayer interactions within layered materials like graphene. However, it is also widely known that this potential strongly underestimates the frictional properties for layered materials. Here we propose to supplement the Lennard-Jones potential by a Gaussian-type potential, which enables more accurate calculations of the frictional properties of two-dimensional layered materials. Furthermore, the Gaussian potential is computationally simple as it introduces only one additional potential parameter that is determined by the interlayer shear mode in the layered structure. The resulting Lennard-Jones-Gaussian potential is applied to compute the interlayer cohesive energy and frictional energy for graphene, MoS2, black phosphorus, and their heterostructures.

preprint2015arXiv

A Review on Flexural Mode of Graphene: Lattice Dynamics, Thermal Conduction, Thermal Expansion, Elasticity, and Nanomechanical Resonance

Single-layer graphene is so flexible that its flexural mode (also called the ZA mode, bending mode, or out-of-plane transverse acoustic mode) is important for its thermal and mechanical properties. Accordingly, this review focuses on exploring the relationship between the flexural mode and thermal and mechanical properties of graphene. We first survey the lattice dynamic properties of the flexural mode, where the rigid translational and rotational invariances play a crucial role. After that, we outline contributions from the flexural mode in four different physical properties or phenomena of graphene -- its thermal conductivity, thermal expansion, Young's modulus, and nanomechanical resonance. We explain how graphene's superior thermal conductivity is mainly due to its three acoustic phonon modes at room temperature, including the flexural mode. Its coefficient of thermal expansion is negative in a wide temperature range resulting from the particular vibration morphology of the flexural mode. We then describe how the Young's modulus of graphene can be extracted from its thermal fluctuations, which are dominated by the flexural mode. Finally, we discuss the effects of the flexural mode on graphene nanomechanical resonators, while also discussing how the essential properties of the resonators, including mass sensitivity and quality factor, can be enhanced.

preprint2015arXiv

A Stillinger-Weber Potential for Single-Layer Black Phosphorus, and the Importance of Cross-Pucker Interactions for Negative Poisson's Ratio and Edge Stress-Induced Bending

The distinguishing structural feature of single-layer black phosphorus is its puckered structure, which leads to many novel physical properties. In this work, we first present a new parameterization of the Stillinger-Weber potential for single-layer black phosphorus. In doing so, we reveal the importance of a cross-pucker interaction term in capturing its unique mechanical properties, such as a negative Poisson's ratio. In particular, we show that the cross-pucker interaction enables the pucker to act as a re-entrant hinge, which expands in the lateral direction when it is stretched in the longitudinal direction. As a consequence, single-layer black phosphorus has a negative Poisson's ratio in the direction perpendicular to the atomic plane. As an additional demonstration of the impact of the cross-pucker interaction, we show that it is also the key factor that enables capturing of the edge stress-induced bending of single-layer black phosphorus that has been reported in ab initio calculations.

preprint2015arXiv

An Analytic Study of Strain Engineering the Electronic Bandgap in Single-Layer Black Phosphorus

We present an analytic study, based on the tight-binding approximation, of strain effects on the electronic bandgap in single-layer black phosphorus. We obtain an expression for the variation of the bandgap induced by a general strain type that includes both tension in and out of the plane and shear, and use this to determine the most efficient strain direction for different strain types, along which the strongest bandgap manipulation can be achieved. We find that the strain direction that enables the maximum manipulation of the bandgap is not necessarily in the armchair or zigzag direction. Instead, to achieve the strongest bandgap modulation, the direction of the applied mechanical strain is dependent on the type of applied strain.

preprint2015arXiv

Anisotropic Young's Modulus for Single-Layer Black Phosphorus: The Third Principle Direction Besides Armchair and Zigzag

We derive an analytic formula for the Young's modulus in single-layer black phosphorus using the valence force field model. By analyzing the directional dependence for the Young's modulus, we explore the third principle direction with direction angle phi_tp = 0.268pi besides armchair and zigzag directions. The maximum Young's modulus value is in the third principle direction. More specifically, the Young's modulus is 52.2 N/m, 85.4 N/m, and 111.4 N/m in the armchair direction, zigzag direction, and the third principle direction, respectively. This new principle direction is of significance for future discussions of other anisotropic properties in the single-layer black phosphorus.

preprint2015arXiv

Effects of electron-phonon interaction on thermal and electrical transport through molecular nano-conductors

The topic of this review is the effects of electron-phonon interaction (EPI) on the transport properties of molecular nano-conductors. A nano-conductor connects to two electron leads and two phonon leads, possibly at different temperatures or chemical potentials. The EPI appears only in the nano-conductor. We focus on its effects on charge and energy transport. We introduce three approaches. For weak EPI, we use the nonequilibrium Green's function method to treat it perturbatively. We derive the expressions for the charge and heat currents. For weak system-lead couplings, we use the quantum master equation approach. In both cases, we use a simple single level model to study the effects of EPI on the system's thermoelectric transport properties. It is also interesting to look at the effect of currents on the dynamics of the phonon system. For this, we derive a semi-classical generalized Langevin equation to describe the nano-conductor's atomic dynamics, taking the nonequilibrium electron system, as well as the rest of the atomic degrees of freedom as effective baths. We show simple applications of this approach to the problem of energy transfer between electrons and phonons.

preprint2015arXiv

Elastic Bending Modulus for Single-Layer Black Phosphorus

We derive an analytic formula for the elastic bending modulus of single-layer black phosphorus (SLBP) based on the valence force field model. The obtained elastic bending modulus is 4.8028 eV and 7.9905 eV along the armchair and zigzag directions in the SLBP, respectively. These values are obviously larger than the bending modulus of 1.4 eV in graphene due to the intrinsic finite thickness for SLBP. Our derivation analytically illustrates that the elastic bending modulus of the SLBP is proportional to the square of the intrinsic thickness of the SLBP.

preprint2015arXiv

Mechanical Strain Effects on Black Phosphorus Nanoresonators

We perform classical molecular dynamics to investigate the effects of mechanical strain on single-layer black phosphorus nanoresonators at different temperatures. We find that the resonant frequency is highly anisotropic in black phosphorus due to its intrinsic puckered configuration, and that the quality factor in the armchair direction is higher than in the zigzag direction at room temperature. The quality factors are also found to be intrinsically larger than graphene and MoS2 nanoresonators. The quality factors can be increased by more than a factor of two by applying tensile strain, with uniaxial strain in the armchair direction being most effective. However, there is an upper bound for the quality factor increase due to nonlinear effects at large strains, after which the quality factor decreases. The tension induced nonlinear effect is stronger along the zigzag direction, resulting in a smaller maximum strain for quality factor enhancement.

preprint2015arXiv

Parametrization of Stillinger-Weber Potential Based on Valence Force Field Model: Application to Single-Layer MoS2 and Black Phosphorus

We propose to parametrize the Stillinger-Weber potential for covalent materials starting from the valence force field model. All geometrical parameters in the Stillinger-Weber potential are determined analytically according to the equilibrium condition for each individual potential term, while the energy parameters are derived from the valence force field model. This parametrization approach transfers the accuracy of the valence force field model to the Stillinger-Weber potential. Furthermore, the resulting Stilliinger-Weber potential supports for stable molecular dynamics simulations, as each potential term is at energy minimum state separately at the equilibrium configuration. We employ this procedure to parametrize Stillinger-Weber potentials for the single-layer MoS2 and black phosphorous. The obtained Stillinger-Weber potentials predict accurate phonon spectrum and mechanical behaviors. We also provide input scripts of these Stillinger-Weber potentials used by publicly available simulation packages including GULP and LAMMPS.

preprint2015arXiv

Registry Effect on the Thermal Conductivity of Few-Layer Graphene

We perform molecular dynamics simulations to study the registry effect on the thermal conductivity of few-layer graphene. The interlayer interaction is described by either the Lennard-Jones potential or the registry-dependent potential. Our calculations show that the thermal conductivity in few-layer graphene from both potentials are close to each other, i.e the registry effect is essentially not important. It is because the thermal transport in few-layer graphene is mainly limited by the interlayer breathing mode, which is insensitive to the registry.

preprint2015arXiv

Strain Engineering for Thermal Conductivity of SWCNT Forests

We perform classical molecular dynamics simulations to investigate the mechanical compression effect on the thermal conductivity of the single-walled carbon nanotube (SWCNT) forest, in which SWCNTs are closely aligned and parallel with each other. We find that the thermal conductivity can be linearly enhanced by increasing compression before the buckling of SWCNT forests, but the thermal conductivity decreases quickly with further increasing compression after the forest is buckled. Our phonon mode analysis reveals that, before buckling, the smoothness of the inter-tube interface is maintained during compression, and the inter-tube van der Waals interaction is strengthened by the compression. Consequently, the twisting-like mode (good heat carrier) is well preserved and its group velocity is increased by increasing compression, resulting in the enhancement of the thermal conductivity. The buckling phenomenon changes the circular cross section of the SWCNT into ellipse, which causes effective roughness at the inter-tube interface for the twisting motion. As a result, in ellipse SWCNTs, the radial breathing mode (poor heat carrier) becomes the most favorable motion instead of the twisting-like mode, leading to the quick decrease of the thermal conductivity with further increasing compression after buckling.

preprint2015arXiv

The Effects of Free Edge Interaction-Induced Knotting on the Buckling of Monolayer Graphene

Edge effects play an important role for many properties of graphene. While most works have focused on the effects from isolated free edges, we present a novel knotting phenomenon induced by the interactions between a pair of free edges in graphene, and investigate its effect on the buckling of monolayer graphene. Upon compression, the buckling of graphene starts gradually in the form of two buckling waves from the warped edges. The collision of these two buckling waves results in the creation of a knot structure in graphene. The knot structure enables the buckled graphene to exhibit two unique post-buckling characteristics. First, it induces a five-fold increase in graphene's mechanical stiffness during the buckling process. Second, the knotted structure enables graphene to exhibit a mechanically stable post-buckling regime over a large (3%) compressive strain regime, which is significantly larger than the critical buckling strain of about 0.5%. The combination of these two effects enables graphene to exhibit an unexpected post-buckling stability that has previously not been reported. We predict that numerical simulations or experiments should observe two distinct stress strain relations for the buckling of identical graphene samples, due to the characteristic randomness in the formation process of the knot structure.

preprint2015arXiv

Thermal Conduction in Single-Layer Black Phosphorus: Highly Anisotropic?

The single-layer black phosphorus is characteristic for its puckered structure, which has leaded to highly anisotropy in its optical, electronic, and mechanical properties. We report, using the non-equilibrium Green's function approach and the first-principles method, that the anisotropy in the thermal conduction is very weak in the single-layer black phosphorus -- the difference between two in-plane directions is less than 4%. Our phonon calculations disclose that the two in-plane acoustic phonon branches have lower group velocities in the direction perpendicular to the pucker, as the black phosphorus is softer in this direction, leading to a weakening effect for the thermal conductance in the perpendicular direction. However, the out-of-plane acoustic phonon branch behaviors abnormally; i.e., the group velocity of this phonon branch is higher in the perpendicular direction, although the single-layer black phosphorus is softer in this direction. The abnormal behavior of the out-of-plane phonon branch is closely related to the highly anisotropic Poisson's ratio in the single-layer black phosphorus. As a result of the counteraction between the in-plane phonon modes and the out-of-plane phonon mode, the thermal conductance in the perpendicular direction is weaker than the parallel direction, but the anisotropy is pretty small.

preprint2014arXiv

A Fast Uniaxial Compression of the Single-Layer MoS2

The Euler buckling theorem states that the buckling critical strain is an inverse square function of the length for a thin plate in the static compression process. However, the suitability of this theorem in the dynamical process is unclear, so we perform molecular dynamics simulations to examine the applicability of the Euler buckling theorem in case of a fast compression of the single-layer MoS2. We find that the Euler buckling theorem is not applicable in such dynamical process, as the buckling critical strain becomes a length-independent constant in the buckled system with many ripples. However, the Euler buckling theorem can be resumed in this dynamical process after restricting the theorem to an individual ripple in the buckled structure.

preprint2014arXiv

Interlayer Breathing and Shear Modes in Few-Layer Black Phosphorus

The interlayer breathing and shear modes in few-layer black phosphorus are investigated for their symmetry and lattice dynamical properties. The symmetry groups for the even-layer and odd-layer few-layer black phosphorus are utilized to determine the irreducible representation and the infrared and Raman activity for the interlayer modes. The valence force field model is applied to calculate the eigenvectorw and frequencies for the interlayer breathing and shear modes, which are explained using the atomic chain model. The anisotropic puckered configuration for black phosphorus leads to a highly anisotropic frequency for the two interlayer shear modes. More specifically, the frequency for the shear mode in the direction perpendicular to the pucker is less than half of the shear mode in the direction parallel with the pucker. We also report a set of interlayer modes having the same frequency for all few-layer black phosphorus with layer number N=3i with integer i, because of their collective vibrational displacements. The optical activity of the collective modes supports possible experimental identification for these modes.

preprint2014arXiv

Mechanical Properties of MoS2/Graphene Heterostructures

We perform classic molecular dynamics simulations to comparatively investigate the mechanical properties of single-layer MoS2 and a graphene/MoS2/graphene heterostructure under uniaxial tension. We show that the lattice mismatch between MoS2 and graphene will lead to an spontaneous strain energy in the interface. The Young's modulus of MoS2 can be enhanced by a factor of five by sandwiching it between two graphene layers. While the stiffness is enhanced, the yield strain of the MoS2 is reduced due to lateral buckling of the outer graphene layers due to the applied mechanical tension.

preprint2014arXiv

Mechanical Properties of Single-Layer Black Phosphorus

The mechanical properties of single-layer black phosphrous under uniaxial deformation are investigated using first-principles calculations. Both Young's modulus and the ultimate strain are found to be highly anisotropic and nonlinear as a result of its quasi-two-dimensional puckered structure. Specifically, the in-plane Young's modulus is 44.0 GPa in the direction perpendicular to the pucker, and 92.7 GPa in the parallel direction. The ultimate strain is 0.48 and 0.20 in the perpendicular and parallel directions, respectively.

preprint2014arXiv

MoS2 Nanoresonators: Intrinsically Better Than Graphene?

We perform classical molecular dynamics simulations to examine the intrinsic energy dissipation in single-layer MoS$_{2}$ nanoresonators, where a point of emphasis is to compare its dissipation characteristics with those of single-layer graphene. Our key finding is that MoS$_{2}$ nanoresonators exhibit significantly lower energy dissipation, and thus higher quality (Q)-factors by at least a factor of four below room temperature, than graphene. Furthermore, this high Q-factor endows MoS$_{2}$ nanoresonators with a higher figure of merit, defined as frequency times Q-factor, despite a resonant frequency that is $50%$ smaller than graphene for the same size. By utilizing arguments from phonon-phonon scattering theory, we show that this reduced energy dissipation is enabled by the large energy gap in the phonon dispersion of MoS$_{2}$, which separates the acoustic phonon branches from the optical phonon branches, leading to a preserving mechanism for the resonant oscillation of MoS$_{2}$ nanoresonators. We further investigate the effects of tensile mechanical strain and nonlinear actuation on the Q-factors, where the tensile strain is found to counteract the reductions in Q-factor that occur with higher actuation amplitudes. Overall, our simulations illustrate the potential utility of MoS$_{2}$ for high frequency sensing and actuation applications.

preprint2014arXiv

Negative Poisson's Ratio in Single-Layer Black Phosphorus

The Poisson's ratio is a fundamental mechanical property that relates the resulting lateral strain to applied axial strain. While this value can theoretically be negative, it is positive for nearly all materials, though negative values have been observed in so-called auxetic structures. However, nearly all auxetic materials are bulk materials whose microstructure has been specifically engineered to generate a negative Poisson's ratio. In the present work, we report using first principles calculations the existence of a negative Poisson's ratio in a single-layer, two-dimensional material, black phosphorus. In contrast to engineered bulk auxetics, this behavior is intrinsic for single layer black phosphorus, and originates from its unique, puckered structure, where the pucker can be regarded as a re-entrant structure that is comprised of two coupled orthogonal hinges. As a result of this atomic structure, a negative Poisson's ratio is observed in the out-of-plane direction under uniaxial deformation in the direction parallel to the pucker, with the Poisson's ratio becoming increasingly negative with both increased tension and compression. The puckered structure also results in highly anisotropic in-plane Poisson's ratios, which are found to be 0.4 in the direction perpendicular and 1.28 in the direction parallel to the pucker.

preprint2014arXiv

Phonon Bandgap Engineering of Strained Monolayer MoS2

The phonon band structure of monolayer MoS2 is characteristic for a large energy gap between acoustic and optical branches, which protects the vibration of acoustic modes from being scattered by optical phonon modes. Therefore, the phonon bandgap engineering is of practical significance for the manipulation of phonon-related mechanical or thermal properties in monolayer MoS2. We perform both phonon analysis and molecular dynamics simulations to investigate the tension effect on the phonon bandgap and the compression induced instability of the monolayer MoS2. Our key finding is that the phonon bandgap can be narrowed by the uniaxial tension, and is completely closed at epsilon=0.145; while the biaxial tension only has limited effect on the phonon bandgap. We also demonstrate the compression induced buckling for the monolayer MoS2. The critical strain for buckling is extracted from the band structure analysis of the flexure mode in the monolayer MoS2 and is further verified by molecular dynamics simulations and the Euler buckling theory. Our study illustrates the uniaxial tension as an efficient method for manipulating the phonon bandgap of the monolayer MoS2, while the biaxial compression as a powerful tool to intrigue buckling in the monolayer MoS2.

preprint2014arXiv

Phonon Modes in Single-Walled Molybdenum Disulphide (MoS2) Nanotubes: Lattice Dynamics Calculation and Molecular Dynamics Simulation

We study the phonon modes in single-walled MoS$_{2}$ nanotubes via the lattice dynamics calculation and molecular dynamics simulation. The phonon spectra for tubes of arbitrary chiralities are calculated from the dynamical matrix constructed by the combination of an empirical potential with the conserved helical quantum numbers $(κ, n)$. In particular, we show that the frequency ($ω$) of the radial breathing mode is inversely proportional to the tube diameter ($d$) as $ω=665.3/d$ {cm$^{-1}$}. The eigen vectors of the first twenty lowest-frequency phonon modes are illustrated. Based on these eigen vectors, we demonstrate that the radial breathing oscillation is disturbed by phonon modes of three-fold symmetry initially, and the tube is squashed by the modes of two-fold symmetry eventually. Our study provides fundamental knowledge for further investigations of the thermal and mechanical properties of the MoS$_{2}$ nanotubes.

preprint2014arXiv

Polar Surface Effects on the Thermal Conductivity in ZnO Nanowires: a Shell-Like Surface Reconstruction-Induced Preserving Mechanism

We perform molecular dynamics (MD) simulations to investigate the effect of polar surfaces on the thermal transport in zinc oxide (ZnO) nanowires. We find that the thermal conductivity in nanowires with free polar (0001) surfaces is much higher than in nanowires that have been stabilized with reduced charges on the polar (0001) surfaces, and also hexagonal nanowires without any transverse polar surfaces. From normal mode analysis, we show that the higher thermal conductivity is due to a shell-like reconstruction that occurs for the free polar surfaces. This shell-like reconstruction suppresses twisting motion in the nanowires such that the bending phonon modes are not scattered by the other phonon modes, and leads to substantially higher thermal conductivity in the ZnO nanowire with free polar surfaces. Furthermore, the auto-correlation function of the normal mode coordinate is utilized to extract the phonon lifetime, which leads to a concise explanation for the higher thermal conductivity in ZnO nanowires with free polar surfaces. Our work demonstrates that ZnO nanowires without polar surfaces, which exhibit low thermal conductivity, are more promising candidates for thermoelectric applications than nanowires with polar surfaces (and also high thermal conductivity).

preprint2014arXiv

Preserving the Q-Factors of ZnO Nanoresonators via Polar Surface Reconstruction

We perform molecular dynamics simulations to investigate the effect of polar surfaces on the quality (Q)-factors of zinc oxide (ZnO) nanowire-based nanoresonators. We find that the Q-factors in ZnO nanoresonators with free polar (0001) surfaces is about one order of magnitude higher than in nanoresonators that have been stabilized with reduced charges on the polar (0001) surfaces. From normal mode analysis, we show that the higher Q-factor is due to a shell-like reconstruction that occurs for the free polar surfaces. This shell-like reconstruction suppresses twisting motion in the nanowires such that the mixing of other modes with the resonant mode of oscillation is minimized, and leads to substantially higher Q-factors in the ZnO nanoresonators with free polar surfaces.

preprint2014arXiv

The Buckling of Single-Layer MoS2 Under Uniaxial Compression

Molecular dynamics simulations are performed to investigate the buckling of single-layer MoS2 under uniaxial compression. The strain rate is found to play an important role on the critical buckling strain, where higher strain rate leads to larger critical strain. The critical strain is almost temperature-independent for T<50 K, and it increases with increasing temperature for T>50 K owning to the thermal vibration assisted healing mechanism on the buckling deformation. The length-dependence of the critical strain from our simulations is in good agreement with the prediction of the Euler buckling theory.

preprint2013arXiv

A Surface Stacking Fault Energy Approach to Predicting Defect Nucleation in Surface-Dominated Nanostructures

We present a surface stacking fault (SSF) energy approach to predicting defect nucleation from the surfaces of surface-dominated nanostructure such as FCC metal nanowires. The approach leads to a criteria that predicts the initial yield mechanism via either slip or twinning depending on whether the unstable twinning energy or unstable slip energy is smaller as determined from the resulting SSF energy curve. The approach is validated through a comparison between the SSF energy calculation and low-temperature classical molecular dynamics simulations of copper nanowires with different axial and transverse surface orientations, and cross sectional geometries. We focus on the effects of the geometric cross section by studying the transition from slip to twinning previously predicted in moving from a square to rectangular cross section for $\ <100\ > /\{100\}$ nanowires, and also for moving from a rhombic to truncated rhombic cross sectional geometry for $\ <110\ >$ nanowires. \hsp{We also provide the important demonstration that the criteria is able to predict the correct deformation mechanism when full dislocation slip is considered concurrently with partial dislocation slip and twinning. This is done in the context of rhombic aluminum nanowires which do not show a tensile reorientation due to full dislocation slip.} We show that the SSF energy criteria successfully predicts the initial mode of surface-nucleated plasticity at low temperature, while also discussing the effects of strain and temperature on the applicability of the criterion.

preprint2013arXiv

Adsorbate Migration Effects on Continuous and Discontinuous Temperature-Dependent Transitions in the Quality Factors of Graphene Nanoresonators

We perform classical molecular dynamics simulation to investigate the mechanisms underpinning the unresolved, experimentally-observed temperature-dependent scaling transition in the quality factors of graphene nanomechanical resonators (GNMR). Our simulations reveal that the mechanism underlying this temperature scaling phenomenon is the out-of-plane migration of adsorbates on GNMRs. Specifically, the migrating adsorbate undergoes frequent collisions with the GNMR, which strongly influences the resulting mechanical oscillation, and thus the quality factors. We also predict a discontinuous transition in the quality factor at a lower critical temperature, which results from the in-plane migration of the adsorbate. Overall, our work clearly demonstrates the strong effect of adsorbate migration on the quality factors of GNMRs.

preprint2013arXiv

Elastic Bending Modulus of Single-Layer Molybdenum Disulphide (MoS2): Finite Thickness Effect

We derive, from an empirical interaction potential, an analytic formula for the elastic bending modulus of single-layer MoS2 (SLMoS2). By using this approach, we do not need to define or estimate a thickness value for SLMoS2, which is important due to the substantial controversy in defining this value for two-dimensional or ultrathin nanostructures such as graphene and nanotubes. The obtained elastic bending modulus of 9.61 eV in SLMoS2 is significantly higher than the bending modulus of 1.4 eV in graphene, and is found to be within the range of values that are obtained using thin shell theory with experimentally obtained values for the elastic constants of SLMoS2. This increase in bending modulus as compared to monolayer graphene is attributed, through our analytic expression, to the finite thickness of SLMoS2. Specifically, while each monolayer of S atoms contributes 1.75 eV to the bending modulus, which is similar to the 1.4 eV bending modulus of monolayer graphene, the additional pairwise and angular interactions between out of plane Mo and S atoms contribute 5.84 eV to the bending modulus of SLMoS2.

preprint2013arXiv

Mechanical Oscillation of Kinked Silicon Nanowires: a Natural Nanoscale Spring

We perform classical molecular dynamics simulations to demonstrate the application of kinked silicon nanowires (KSiNWs) as nanoscale springs. The spring-like oscillation in gigahertz frequency range is successfully actuated using a similar procedure as the actuation of a classical mass spring oscillator. We detect the spring-like mechanical oscillation and some other low-frequency oscillations by the energy spectrum analysis, where a dimensional crossover phenomenon is observed for the transverse mode in KSiNWs with decreasing aspect ratio. Our findings shed light on the elastic properties of the KSiNW and open a way for its application in nanomechanical devices.

preprint2013arXiv

Modulation of Thermal Conductivity in Kinked Silicon Nanowires: Phonon Interchanging and Pinching Effects

We perform molecular dynamics simulations to investigate the reduction of the thermal conductivity by kinks in silicon nanowires. The reduction percentage can be as high as 70% at room temperature. The temperature dependence of the reduction is also calculated. By calculating phonon polarization vectors, two mechanisms are found to be responsible for the reduced thermal conductivity: (1) the interchanging effect between the longitudinal and transverse phonon modes and (2) the pinching effect, i.e a new type of localization, for the twisting and transverse phonon modes in the kinked silicon nanowires. Our work demonstrates that the phonon interchanging and pinching effects, induced by kinking, are brand new and effective ways in modulating heat transfer in nanowires, which enables the kinked silicon nanowires to be a promising candidate for thermoelectric materials.

preprint2013arXiv

Molecular Dynamics Simulations of Single-Layer Molybdenum Disulphide (MoS2): Stillinger-Weber Parametrization, Mechanical Properties, and Thermal Conductivity

We present a parameterization of the Stillinger-Weber potential to describe the interatomic interactions within single-layer MoS2 (SLMoS2). The potential parameters are fitted to an experimentally-obtained phonon spectrum, and the resulting empirical potential provides a good description for the energy gap and the crossover in the phonon spectrum. Using this potential, we perform classical molecular dynamics simulations to study chirality, size, and strain effects on the Young&#39;s modulus and the thermal conductivity of SLMoS2. We demonstrate the importance of the free edges on the mechanical and thermal properties of SLMoS2 nanoribbons. Specifically, while edge effects are found to reduce the Young&#39;s modulus of SLMoS2 nanoribbons, the free edges also reduce the thermal stability of SLMoS2 nanoribbons, which may induce melting well below the bulk melt temperature. Finally, uniaxial strain is found to efficiently manipulate the thermal conductivity of infinite, periodic SLMoS2.

preprint2013arXiv

Orientation Dependent Thermal Conductance in Single-Layer MoS2

We investigate the thermal conductivity in the armchair and zigzag MoS$_{2}$ nanoribbons, by combining the non-equilibrium Green&#39;s function approach and the first-principles method. A strong orientation dependence is observed in the thermal conductivity. Particularly, the thermal conductivity for the armchair MoS$_{2}$ nanoribbon is about 673.6 Wm$^{-1}$K$^{-1}$ in the armchair nanoribbon, and 841.1 Wm$^{-1}$K$^{-1}$ in the zigzag nanoribbon at room temperature. By calculating the Caroli transmission, we disclose the underlying mechanism for this strong orientation dependence to be the fewer phonon transport channels in the armchair MoS$_{2}$ nanoribbon in the frequency range of [150, 200] {cm$^{-1}$}. Through the scaling of the phonon dispersion, we further illustrate that the thermal conductivity calculated for the MoS$_{2}$ nanoribbon is esentially in consistent with the superior thermal conductivity found for graphene.

preprint2013arXiv

Size-Sensitive Young&#39;s modulus of Kinked Silicon Nanowires

We perform both classical molecular dynamics simulations and beam model calculations to investigate the Young&#39;s modulus of kinked silicon nanowires (KSiNWs). The Young&#39;s modulus is found to be highly sensitive to the arm length of the kink and is essentially inversely proportional to the arm length. The mechanism underlying the size dependence is found to be the interplay between the kink angle potential and the arm length potential, where we obtain an analytic relationship between the Young&#39;s modulus and the arm length of the KSiNW. Our results provide insight into the application of this novel building block in nanomechanical devices.

preprint2013arXiv

Superior thermal conductivity and extremely high mechanical strength in polyethylene chains from {\it ab initio} calculation

The upper limit of the thermal conductivity and the mechanical strength are predicted for the polyethylene chain, by performing the {\it ab initio} calculation and applying the quantum mechanical non-equilibrium Green&#39;s function approach. Specially, there are two main findings from our calculation: (1). the thermal conductivity can reach a high value of 310 W/K/m in a 100 nm polyethylene chain at room temperature; (2). the Young&#39;s modulus in the polyethylene chain is as high as 374.5 GPa, and the polyethylene chain can sustain $32.85% \pm 0.05%$ (ultimate) strain before undergoing structural phase transition into gaseous ethylene.

preprint2013arXiv

Why twisting angles are diverse in graphene Moire patterns?

The interlayer energy of the twisting bilayer graphene is investigated by the molecular mechanics method using both the registry-dependent potential and the Lennard-Jones potential. Both potentials show that the interlayer energy is independent of the twisting angle $θ$, except in the two boundary regions $θ\approx 0$ or $60^{\circ}$, where the interlayer energy is proportional to the square of the twisting arc length. The calculation results are successfully interpreted by a single atom model. An important information from our findings is that, from the energy point of view, there is no preference for the twisting angle in the experimental bilayer graphene samples, which actually explains the diverse twisting angles in the experiment.

preprint2012arXiv

A comparative study of two molecular mechanics models based on harmonic potentials

We show that the two molecular mechanics models, the stick-spiral and the beam models, predict considerably different mechanical properties of materials based on energy equivalence. The difference between the two models is independent of the materials since all parameters of the beam model are obtained from the harmonic potentials. We demonstrate this difference for finite width graphene nanoribbons and a single polyethylene chain comparing results of the molecular dynamics (MD) simulations with harmonic potentials and the finite element method with the beam model. We also find that the difference strongly depends on the loading modes, chirality and width of the graphene nanoribbons, and it increases with decreasing width of the nanoribbons under pure bending condition. The maximum difference of the predicted mechanical properties using the two models can exceed 300% in different loading modes. Comparing the two models with the MD results of AIREBO potential, we find that the stick-spiral model overestimates and the beam model underestimates the mechanical properties in narrow armchair graphene nanoribbons under pure bending condition.

preprint2012arXiv

Enhancing the Mass Sensitivity of Graphene Nanoresonators Via Nonlinear Oscillations: The Effective Strain Mechanism

We perform classical molecular dynamics simulations to investigate the enhancement of the mass sensitivity and resonant frequency of graphene nanomechanical resonators that is achieved by driving them into the nonlinear oscillation regime. The mass sensitivity as measured by the resonant frequency shift is found to triple if the actuation energy is about 2.5 times the initial kinetic energy of the nanoresonator. The mechanism underlying the enhanced mass sensitivity is found to be the effective strain that is induced in the nanoresonator due to the nonlinear oscillations, where we obtain an analytic relationship between the induced effective strain and the actuation energy that is applied to the graphene nanoresonator. An important implication of this work is that there is no need for experimentalists to apply tensile strain to the resonators before actuation in order to enhance the mass sensitivity. Instead, enhanced mass sensitivity can be obtained by the far simpler technique of actuating nonlinear oscillations of an existing graphene nanoresonator.

preprint2012arXiv

How Does Folding Modulate Thermal Conductivity of Graphene

We study thermal transport in folded graphene nanoribbons using molecular dynamics simulations and the non-equilibrium Green&#39;s function method. It is found that the thermal conductivity of flat graphene nanoribbons can be modulated by folding and changing interlayer couplings. The analysis of transmission reveals that the reduction of thermal conductivity is due to scattering of low frequency phonons by the folds. Our results suggest that folding can be utilized in the modulation of thermal transport properties in graphene and other two dimensional materials.

preprint2012arXiv

Why edge effects are important on the intrinsic loss mechanisms of graphene nanoresonators?

Molecular dynamics simulations are performed to investigate edge effects on the quality factor of graphene nanoresonators with different edge configurations and of various sizes. If the periodic boundary condition is applied, very high quality factors ($3\times10^{5}$) are obtained for all kinds of graphene nanoresonators. However, if the free boundary condition is applied, quality factors will be greatly reduced by two effects resulting from free edges: the imaginary edge vibration effect and the artificial effect. Imaginary edge vibrations will flip between a pair of doubly degenerate warping states during the mechanical oscillation of nanoresonators. The flipping process breaks the coherence of the mechanical oscillation of the nanoresonator, which is the dominant mechanism for extremely low quality factors. There is an artificial effect if the mechanical oscillation of the graphene nanoresonator is actuated according to an artificial vibration (non-natural vibration of the system), which slightly reduce the quality factor. The artificial effect can be eliminated by actuating the mechanical oscillation according to a natural vibration of the nanoresonator. Our simulations provide an explanation for the recent experiment, where the measured quality factor is low and varies between identical samples with free edges.

preprint2011arXiv

A nonequilibrium Green&#39;s function study of thermoelectric properties in single-walled carbon nanotubes

The phonon and electron transport in single-walled carbon nanotubes (SWCNT) are investigated using the nonequilibrium Green&#39;s function approach. In zigzag SWCNT ($n$, 0) with $mod(n,3)\not=0$, the thermal conductance is mainly attributed to the phonon transport, while the electron only has few percentage contribution. The maximum value of the figure of merit ($ZT$) is about 0.2 in this type of SWCNT. The $ZT$ is considerably larger in narrower SWCNT because of enhanced Seebeck coefficient. $ZT$ is smaller in the armchair SWCNT, where Seebeck coefficient is small due to zero band gap. It is found that the cluster isotopic doping can reduce the phonon thermal conductance obviously and enhance the value of $ZT$. The uniaxial elongation and compress strain depresses phonons in whole frequency region, leading to the reduction of the phonon thermal conductance in whole temperature range. Interestingly, the elongation strain can affect the phonon transport more seriously than the compress strain, because the high frequency $G$ mode is completely filtered out under elongation strain $ε>0.05$. The strain also has important effect on the subband edges of the electron band structure by smoothing the steps in the electron transmission function. The $ZT$ is decreased by strain as the reduction in the electronic conductance overcomes the reduction in the thermal conductance.

preprint2011arXiv

A theoretical study of thermal conductivity in single-walled boron nitride nanotubes

We perform a theoretical investigation on the thermal conductivity of single-walled boron nitride nanotubes (SWBNT) using the kinetic theory. By fitting to the phonon spectrum of boron nitride sheet, we develop an efficient and stable Tersoff-derived interatomic potential which is suitable for the study of heat transport in sp2 structures. We work out the selection rules for the three-phonon process with the help of the helical quantum numbers $(κ, n)$ attributed to the symmetry group (line group) of the SWBNT. Our calculation shows that the thermal conductivity $κ_{\rm ph}$ diverges with length as $κ_{\rm ph}\propto L^β$ with exponentially decaying $β(T)\propto e^{-T/T_{c}}$, which results from the competition between boundary scattering and three-phonon scattering for flexure modes. We find that the two flexure modes of the SWBNT make dominant contribution to the thermal conductivity, because their zero frequency locates at $κ=\pmα$ where $α$ is the rotational angle of the screw symmetry in SWBNT.

preprint2011arXiv

A universal exponential factor in the dimensional crossover from graphene to graphite

A universal exponential factor, $γ_{c}=π/2$, is disclosed for the dimensional crossover of few-layer graphene (FLG) from two-dimensional graphene to three-dimensional graphite. $γ_{c}$ is found by analyzing available experimental data on different properties of FLG with varying thickness. A theoretical study on the phonon spectrum of the vertical acoustic mode in FLG is carried out to further check this exponential factor $γ_{c}$. Interestingly, the same exponential factor appears in the dimensional crossover of the phonon mode. It turns out that the exponential factor $γ_{c}$ is related to the homogeneous Helmholtz-like molal equation in the mass transfer with a first order chemical reaction. The finding should provide valuable information for experimentalists and theorists in the future investigation on thickness dependent properties of FLG.

preprint2011arXiv

Bright and dark modes induced by graphene bubbles

Through a lattice dynamics analysis, it is revealed that the bubble plays a role of energy shield in the graphene, which helps to split the normal modes into two categories of distinct topological nature, namely the bright and dark modes. The topological invariants, Euler characteristic, of the bright and dark modes are 1 and 0, respectively. For bright modes, the energy is confined inside the bubble, so this type of modes are sensitive to the shape of the bubble; while opposite phenomenon is observed for the dark modes. The different behavior from these two types of normal modes is examined and verified in the process of phonon thermal transport. The bright and dark modes are expected to be distinguished in experiment with existing scanning force microscope techniques, and they should play significant roles in many other physical processes.

preprint2011arXiv

First principle study of the thermal conductance in graphene nanoribbon with vacancy and substitutional silicon defect

The thermal conductance in graphene nanoribbon with a vacancy or silicon point defect (substitution of C by Si atom) is investigated by non-equilibrium Green&#39;s function (NEGF) formalism combined with first-principle calculations density-functional theory with local density approximation. An efficient correction to the force constant matrix is presented to solve the conflict between the long-range character of the {\it ab initio} approach and the first-nearest-neighboring character of the NEGF scheme. In nanoribbon with a vacancy defect, the thermal conductance is very sensitive to the position of the vacancy defect. A vacancy defect situated at the center of the nanoribbon generates a saddle-like surface, which greatly reduces the thermal conductance by strong scattering to all phonon modes; while an edge vacancy defect only results in a further reconstruction of the edge and slightly reduces the thermal conductance. For the Si defect, the position of the defect plays no role for the value of the thermal conductance, since the defective region is limited within a narrow area around the defect center.

preprint2011arXiv

Graphene-based tortional resonator from molecular dynamics simulation

Molecular dynamics simulations are performed to study graphene-based torsional mechanical resonators. The quality factor is calculated by $Q_{F}=ωτ/2π$, where the frequency $ω$ and life time $τ$ are obtained from the correlation function of the normal mode coordinate. Our simulations reveal the radius-dependence of the quality factor as $Q_{F}=2628/(22R^{-1}+0.004R^{2})$, which yields a maximum value at some proper radius $R$. This maximum point is due to the strong boundary effect in the torsional resonator, as disclosed by the temperature distribution in the resonator. Resulting from the same boundary effect, the quality factor shows a power law temperature-dependence with power factors bellow 1.0. The theoretical results supply some valuable information for the manipulation of the quality factor in future experimental devices based on the torsional mechanical resonator.

preprint2011arXiv

Joule heating and thermoelectric properties in short single-walled carbon nanotubes: electron-phonon interaction effect

The electron-phonon interaction (EPI) effect in single-walled carbon nanotube is investigated by the nonequilibrium Green&#39;s function approach within the Born approximation. Special attention is paid to the EPI induced Joule heating phenomenon and the thermoelectric properties in both metallic armchair (10, 10) tube and semiconductor zigzag (10, 0) tube. For Joule heat in the metallic (10, 10) tube, the theoretical results for the breakdown bias voltage is quite comparable with the experimental value. It is found that the Joule heat can be greatly enhanced by increasing the chemical potential, while the role of the temperature is not so important for Joule heat. In zigzag (10, 0) tube, the Joule heat is smaller than the armchair tube, resulting from nonzero band gap in the electron band structure. For the electronic conductance $G_{e}$ and electron thermal conductance $σ_{el}$, the EPI has important effect at higher temperature or higher chemical potential. Compared with ballistic transport, there is an opposite tendency for $G_{e}$ to decrease with increasing temperature after EPI is considered. This is due to the dominant effect of the electron phonon scattering mechanism in the electron transport in this situation. There is an interesting `electron-drag&#39; phenomenon for the phonon thermal conductance in case of low temperature and high chemical potential, where phonons are dragged by electrons from low temperature region into high temperature region through EPI effect.

preprint2011arXiv

Manipulation of heat current by the interface between graphene and white graphene

We investigate the heat current flowing across the interface between graphene and hexagonal boron nitride (so-called white graphene) using both molecular dynamics simulation and nonequilibrium Green&#39;s function approaches. These two distinct methods discover the same phenomena that the heat current is reduced linearly with increasing interface length, and the zigzag interface causes stronger reduction of heat current than the armchair interface. These phenomena are interpreted by both the lattice dynamics analysis and the transmission function explanation, which both reveal that the localized phonon modes at interfaces are responsible for the heat management. The room temperature interface thermal resistance is about $7\times10^{-10}$m$^{2}$K/W in zigzag interface and $3.5\times10^{-10}$m$^{2}$K/W in armchair interface, which directly results in stronger heat reduction in zigzag interface. Our theoretical results provide a specific route for experimentalists to control the heat transport in the graphene and hexagonal boron nitride compound through shaping the interface between these two materials.

preprint2011arXiv

Minimum thermal conductance in graphene and boron nitride superlattice

The minimum thermal conductance versus supercell size ($d_{s}$) is revealed in graphene and boron nitride superlattice with $d_{s}$ far below the phonon mean free path. The minimum value is reached at a constant ratio of $d_{s}/L\approx 5%$, where $L$ is the total length of the superlattice; thus the minimum point of $d_{s}$ depends on $L$. The phenomenon is attributed to the localization property and the number of confined modes in the superlattice. With the increase of $d_{s}$, the localization of the confined mode is enhanced while the number of confined modes decreases, which directly results in the minimum thermal conductance.

preprint2011arXiv

Molecular dynamics simulation for heat transport in thin diamond nanowires

The phonon thermal conductivity in diamond nanowires (DNW) is studied by molecular dynamics simulation. It is found that the thermal conductivity in narrower DNW is lower and does not show obvious temperature dependence; a very small value (about 2.0 W/m/K) of thermal conductivity is observed in ultra-narrow DNW, which may be of potential applications in thermoelectric devices. These two phenomena are probably due to the dominant surface effect and phonon confinement effect in narrow DNW. Our simulation reveals a high anisotropy in the heat transport of DNW. Specifically, the thermal conductivity in DNW along [110] growth direction is about five times larger than that of [100] and [111] growth directions. The anisotropy is believed to root in the anisotropic group velocity for acoustic phonon modes in DNW along three different growth directions.

preprint2011arXiv

Phonon modes at finite temperature in graphene and single-walled carbon nanotubes

The phonon modes at finite temperature in graphene and single-walled carbon nanotubes (SWCNT) are investigated by the mode projection technique combined with molecular dynamics simulation. It is found that the quadratic phonon spectrum of the flexural mode in graphene is renormalized into linear at finite temperature by the phonon-phonon scattering. Possible influence of this renormalization on the electron-phonon interaction and phonon thermal transport is analyzed. For the SWCNT, however the quadratic property of the phonon spectrum for the flexural mode can survive at finite temperature under protection of the two-fold degeneracy of the flexural mode in SWCNT due to its quasi-one-dimensional structure. The frequency and life time of optical phonon modes are also studied at different temperatures. In particular, the life time of the in-plane optical phonons in graphene and the axial optical phonon in SWCNT can be described by an uniform formula $τ(T)=560/T$, which coincides with the experimental results.

preprint2011arXiv

Thermal expansion in multiple layers of graphene

In this file, we apply the nonequilibrium Green&#39;s function method to calculate the coefficient of thermal expansion (CTE) in multiple layers of graphene. We focus on the effect from different layer number $N$. Our main prediction is that: the increase of $N$ can either enhance or weaken CTE, depending on the strength of the substrate interaction $γ$. If $γ< ε$, where $ε$ is the inter-layer interaction, the CTE will increase with increasing $N$. Otherwise, if $γ> ε$, CTE will decrease with increasing $N$.

preprint2011arXiv

Topological effect on thermal conductivity in graphene

The topological effect on thermal conductivity is investigated through the comparison among graphene nanoribbons, carbon nanotubes and the Mobius-like graphene strips (MGS), by molecular dynamics simulation. It is found that the thermal conductivity of MGS is less than one half of that of graphene nanoribbons. The underlying mechanism whereby MGS acquire such low thermal conductivity may be attributable to the enhanced phonon-phonon scattering and localization property, which are induced by the nontrivial topology of Mobius strip. Moreover by counting in the dimensions of MGS, a lower length/width ratio reduces its thermal conductivity, as the phonon-phonon scattering and localization within might be further elevated.

preprint2010arXiv

Isotopic effects on the thermal conductivity of graphene nanoribbons: localization mechanism

Thermal conductivity of graphene nanoribbons (GNR) with length 106~Å and width 4.92~Å after isotopic doping is investigated by molecular dynamics with quantum correction. Two interesting phenomena are found: (1) isotopic doping reduces thermal conductivity effectively in low doping region, and the reduction slows down in high doping region; (2) thermal conductivity increases with increasing temperature in both pure and doped GNR; but the increasing behavior is much more slowly in the doped GNR than that in pure ones. Further studies reveal that the physics of these two phenomena is related to the localized phonon modes, whose number increases quickly (slowly) with increasing isotopic doping in low (high) isotopic doping region.

preprint2010arXiv

Mechanism of phonon localized edge modes

The phonon localized edge modes are systematically studied, and two conditions are proposed for the existence of the localized edge modes: (I) coupling between different directions ($x$, $y$ or $z$) in the interaction; (II) different boundary conditions in three directions. The generality of these two conditions is illustrated by different lattice structures: one-dimensional (1D) chain, 2D square lattice, 2D graphene, 3D simple cubic lattice, 3D diamond structure, etc; and with different potentials: valence force field model, Brenner potential, etc.

preprint2010arXiv

Self-repairing in single-walled carbon nanotubes by heat treatment

Structure transformation by heat treatment in single-walled carbon nanotubes (SWCNT) is investigated using molecular dynamics simulation. The critical temperature for the collapse of pure SWCNT is as high as 4655 K due to strong covalent carbon-carbon bonding. Above 2000 K, the cross section of SWCNT changes from circle to ellipse. The self-repairing capability is then investigated and two efficient processes are observed for the SWCNT to repair themselves. (1) In the first mechanism, vacancy defects aggregate to form a bigger hole, and a bottleneck junction is constructed nearby. (2) In the second mechanism, a local curvature is generated around the isolate vacancy to smooth the SWCNT. Benefit from the powerful self-repairing capability, defective SWCNT can seek a stable configuration at high temperatures; thus the critical temperature for collapse is insensitive to the vacancy defect density.

preprint2010arXiv

Thermal contraction in silicon nanowires at low temperatures

The thermal expansion effect of silicon nanowires (SiNW) in [100], [110] and [111] directions with different sizes is theoretically investigated. At low temperatures, all SiNW studied exhibit thermal contraction effect due to the lowest energy of the bending vibration mode which has negative effect on the coefficient of thermal expansion (CTE). The CTE in [110] direction is distinctly larger than the other two growth directions because of the anisotropy of the bending mode in SiNW. Our study reveals that CTE decreases with an increase of the structure ratio $γ=length/diameter$, and is negative in whole temperature range with $γ=1.3$.