Abstract
Heterostructures based on two-dimensional (2D) materials have attracted intense attention in recent decades due to their unusual and tunable physics/chemical properties, which can be converted into promising engineering applications ranging from electronics, photonics, and phononics to energy recovery. A fundamental understanding of thermal transport in 2D heterostructures is crucial importance for developing micro-nano devices based on them. In this review, we summarized the recent advances of thermal transport in 2D heterostructures. Firstly, we introduced diverse theoretical approaches and experimental techniques for thermal transport in low-dimensional materials. Then we briefly reviewed the thermal properties of various 2D single-phase materials beyond graphene such as hexagonal boron nitride (h-BN), phosphorene, transition metal dichalcogenides (TMDs) and borophene, and emphatically discussed various influencing factors including structural defects, mechanical strain, and substrate interactions. Moreover, we highlighted thermal conduction control in tailored nanosystems—2D heterostructures and presented the associated underlying physical mechanisms, especially interface-modulated phonon dynamics. Finally, we outline their significant applications in advanced thermal management and thermoelectrics conversion, and discuss a number of open problems on thermal transport in 2D heterostructures.
Introduction
The excellent physicochemical properties of graphene (superior thermal stability, high conductivity, favorable mechanical strength, broadband absorption) have inspired great interest not only in its fundamental and applied aspects, but also in exploring other two-dimensional (2D) materials. Subsequently, numerous post-graphene 2D materials ranging from hexagonal boron nitride (h-BN) (Jo et al., 2013), silicene (Vogt et al., 2012), transition-metal dichalcogenide (TMDC) like MoS2 (Smithe et al., 2018) and MoSe2 (Zhang et al., 2013), to phospherene (Chen et al., 2018c) and borophene (Ranjan et al., 2019) have been synthesized, and present diversified electronic, thermal and optical properties, which could be complementary to those of graphene. For example, compared to the zero band-gap of graphene, MoS2 exhibit direct band gap, making it an attractive material for nanoelectronic applications like field-effect transistors with a high on–off ratio and low power consumption (Chang et al., 2013). Another example, stanene, a single-layer buckled honeycomb structure of Tin atom, is demonstrated to be near-room-temperature quantum anomalous Hall effect (Wu et al., 2014) and ultra-low thermal conductivity (Cherukara et al., 2016), which is particularly suited for thermoelectric devices. Furthermore, 2D organic carbon nitride materials as graphene derivatives exhibited high photocatalytic activity owing to their interface interaction, and can be used as the next-generation highly effective photocatalyst (Zhao et al., 2018a; Zhao et al., 2021).
Recently, research on heterostructures composed of distinct 2D materials has captured much attention because their emergence could bring novel exciting physical properties beyond either of the individual component monolayer structures and provide additional degrees of freedom for device engineering. Normally, heterostructures are divided into two classes: in-plane heterostructures where two 2D materials are stitched seamlessly into a shared plane with an atomically sharp 1D interface; van der Waals (vdW) heterostructures where distinct 2D materials are stacked vertically with non-bonded vdW forces between adjacent layers. Benefited from the progresses in the microtechnology, high-quality graphene/h-BN in-plane heterostructures can be synthesized via chemical vapor deposition (CVD) method, and their compositional and structural diversities could translate into greater freedom for tuning the physical properties such as the presence of a metal-insulator transition (Ci et al., 2010; Liu et al., 2014d). Based on the same preparation method, the integration of MoS2, MoSe2 or WSe2 monolayers into a heterostructure has provided the possibility to produce strong localized photoluminescence and in-plane p-n junctions, showing perspectives in phototransistors and field-effect transistor (Huang et al., 2014; Li et al., 2015). The advantages of in-plane heterostructures mainly contain simpler band alignment, more distinct phase separation and strong interfacial action (Li et al., 2015). Compared to in-plane heterostructures which demands the matched lattice structures of constituent materials, vdW heterostructures intrinsically create an ultra-clean interface regardless of the typical interfacial lattice-matching constraints, which is crucial for tunneling devices that suffer from interfacial defects and dislocations. Additionally, vdW heterostructures offer new mechanisms for manipulating the properties of emergent device via the stacking sequence, the relative rotation between adjacent layers as well as interlayer spacing. A few pioneering works have reported prototypical field-effect-transistor with graphene-based (Bertolazzi et al., 2013; Georgiou et al., 2013) or phosphorus vdW heterostructures (Kang et al., 2016), photodiodes based on MoS2/WSe2 stacks (Roy et al., 2015), and high-efficiency photovoltaic cells based on graphene/WS2 vertical junctions (Iannaccone et al., 2018). On the other hand, the non-bonded vdW interactions of adjacent layers can preserve the electronic properties of each layer, and help resolve a critical challenge that the poor electronic quality of 2D layered materials deposited on SiO2 substrates (Yankowitz et al., 2019). For instance, the quality of graphene encapsulated in h-BN is dramatically improved, and possesses high mobility and ballistic transport over tens of micrometers, which approaches those of suspended graphene (Wang et al., 2013a). It is reported h-BN/phosphorus/h-BN sandwiched heterostructure exists high field-effect mobility and high on-off ratios at room temperature, and the embed phosphorus remains intact and preserves high quality (Chen et al., 2015b). In a word, the physicochemical properties of different heterostructures are closely linked to material composition and structural configuration, representing the significance of structure-property relations.
Thermal transport in 2D heterostructures is not only of a fundamental issue in physics but also crucial to engineering applications. On the one hand, heat dissipation arising from charge carrier scattering processes has to be considered in the development of such nanostructures for electronic and photoelectric devices because the generation of waste Joule heat often leads to the formation of high-temperature hot spots in devices (Wang et al., 2017d), which severely affect the stability and performance of devices. Hence, the choice of 2D building blocks with high thermal conductivity helps to alleviate Joule heat dissipation, and a deeper understanding of the underlying mechanisms for heat conduction across the heterogeneous interface would aid in the design of functional devices. On the other hand, low thermal conductivity materials are necessary for the application in thermoelectric devices converting waste heat into electrical energy. The introduction of 2D heterostructures probably causes increasing Seebeck coefficient and lowering thermal conductivity based on carrier-energy filtering or quantum confinement effects. Compared to dissipating heat and thermoelectric generator, the design of functional thermal devices such as rectifier, transistor, and logical gate to handle heat signals requires controllable thermal transport. In fact, the structural diversities combined with structural defect and external stress field can realize the rich diversity of thermal properties in 2D heterostructures, and enables them to be used in different thermal fields. From another perspective, the interfaces of heterostructures make the thermal transport mechanisms deviate from that in single 2D materials and conventional bulk materials, either through inducing a contact thermal resistance or new phonon modes. For example, it is reported that the thermal conductivity of graphene and graphene nanoribbons (GNRs) monotonously decreases with the tensile strain due to phonon softening (Wei et al., 2011; Shen et al., 2014; Guo et al., 2020), while the load of tensile strain of graphene/h-BN in-plane heterostructures significantly enhance the interfacial thermal conductance owing to improving the alignment of out-of-plane phonon modes near the interface (Ong et al., 2016). Further research demonstrates that the thermal conductivity of graphene/h-BN superlattices first increases and then decreases with increasing the tensile strain, which stems from the competition between in-plane phonons softening and flexural phonons stiffening (Zhu and Ertekin, 2015). Furthermore, some abnormal thermal phenomena such as thermal rectification (TR) (Yuan et al., 2015; Zhang et al., 2018), negative differential thermal resistance (NDTR) (Chen et al., 2017b), thermal bandgap (Zhu and Ertekin, 2014b) and thermophoresis (Mandelli et al., 2017) have been reported in 2D heterostructures, and sparked heated debate.
As we know, phonons are the main heat carriers in semiconductors like TMDCs and silicene, and insulators like h-BN. Specifically, the ultra-high thermal conductivity of graphene dominated by phonons has raised the exciting prospect for heat removal in electronic devices and peculiar phonon transport properties, as reviewed by (Nika and Balandin, 2012; Xu et al., 2014c; Al and Farías, 2016; Xu et al., 2016b; Nika and Balandin, 2017) which has also laid the solid basis for other emerging 2D materials. Phonon transport in novel 2D materials and their heterostuctures probably differs from that in graphene because of crystal structures and atom species, suggesting that the existing knowledge of graphene couldn’t be transplanted to them directly. Fortunately, sustaining advances of theory and experiment offers microscopic details from the atomistic level to the device scale. To date, the thermal properties of 2D materials beyond graphene have been discussed in many previous review articles. Regarding to the progress of theoretical approaches and experimental techniques in 2D materials, please refer to (Cahill et al., 2014; Wang et al., 2017d; Zeng et al., 2018b; Fu et al., 2019) For the phonon transport in phosphorene, borophenes, TMDs and other 2D materials, please refer to (Gu et al., 2018; Qin and Hu, 2018; Chen and Chen, 2020a; Li et al., 2020a; Zhao et al., 2020b) For a detailed survey of thermoelectricity applications of 2D materials, please refer to (Dollfus et al., 2015; Zhang and Zhang, 2017a; Li et al., 2020b). In this review, we not only summarizes the existing results of phonon transport in pristine 2D materials and their heterostuctures, but also gives a physical picture of various phonon scattering mechanisms induced by interior imperfections or external field, to pave the way for better thermal management in practical devices. In Theoretical Approaches and Experimental Measurements, we introduced diverse theoretical approaches and experimental techniques for investigating heat conduction in low-dimensional nanostructures. In Thermal Transport in Some Typical 2D Materials, we briefly reviewed the unique thermal properties of several typical 2D materials. In Thermal Transport in 2D Heterostructures, we highlighted the manipulation of heat conduction in tailored nanosystems—2D heterostructures and present the associated underlying transport mechanisms, especially interface-modulated phonon dynamics. Finally, we show the potential applicability of 2D heterostructures in advanced thermal and thermoelectrics devices, and give the conclusions and a brief outlook.
Theoretical Approaches and Experimental Measurements
Theoretical Methods
Currently, different theoretical method have been developed to study the thermal transport properties of 2D materials, principally including the Boltzmann transport equation (BTE), atomistic Green’s functions (AGF) and molecular dynamics (MD) simulations. Note that the latter two are especially suitable for describing phonon transmission/scattering across an interface in 2D heterostructures.
Molecular Dynamics Simulations
MD simulations with classical potentials are a powerful technique to handle many-body problems at the atomistic level. Here, all-order anharmonic effect is inherently contained in atomistic interactions, differing from that in BTE and AGF. At present, the MD method is widely adopted to study phonon transport problems such as structural effect and various influences like doping, defect, strain and substate. Note that the accuracy of classical MD calculations depends on the quality of the inter-atomic potentials. Recently, intensive efforts have been conducted to construct reliable potential function from the calculations of first principles (Rowe et al., 2018). In practice, equilibrium MD (EMD) and nonequilibrium MD (NEMD) simulations are adopted to extract the thermal properties of 2D materials. In EMD, the calculation of thermal conductivity is based on the Green–Kubo formula or Einstein Relation ration (Sellan et al., 2010; Sevik et al., 2011). For NEMD, both sides of sample are connected to a hot and cold reservoir, respectively. After sufficient run-time, the stationary-state heat flux and temperature are obtained to calculate thermal conductivity in terms of Fourier’s law of heat conduction,
Obviously, the obtained thermal conductivity using the two methods does not provide microscopic perspective into thermal conduction in crystals. For this, a spectral energy density (SED) technique on the base of harmonic lattice dynamics and EMD calculations has been developed to extract mode–resolved thermal conductivity (Thomas et al., 2010; Qiu and Ruan, 2012a). However, the contribution of some long wavelengths phonons can’t be contained in the simulation results because of size effect (Xie et al., 2016; Gu et al., 2018). As a result, the calculated thermal conductivity is sensitive to the size of simulation domain. Note that the NEMD method has also been applied to study interface-modulated phonon dynamics such as graphene/MoS2 (Liu et al., 2017a), graphene/phosphorene in-plane heterostructures (Liu et al., 2018a) and graphene/h-BN (Xu et al., 2014a) MoS2/WSe2 (Zhang et al., 2017b), phosphorene/graphene vdW heterostructures (Liang et al., 2018).
The presence of phonon-interface scattering will produce thermal resistance, also called as Kapitza thermal resistance. The phenomenon could be partly explained by the acoustic mismatch model (AMM) and the diffuse mismatch model (DMM), which is good for some simple material (Kakodkar and Feser, 2017). While for 2D heterostructrues, both kinds of theoretical models are likely to fail because of the dimensionality mismatch and limited thickness of actual system. In addition to this, the models couldn’t give a detailed picture for interfacial thermal transport as AGF calculation. For this, a spectral decomposition of interfacial thermal conductance method based on NEMD simulations is developed by (Sääskilahti et al., 2015) and extended by (Zhou and Hu, 2017c), offering a detailed description of inelastic effects. Under nonequilibrium steady state, the spectrally resolved heat flux between two atoms i and j is given bywhere is the simulation time, is the discrete Fourier transforms of atomic velocity and is the two-order force constant matrix calculated by the finite displacement. “L” and “R” denotes the left and right sides along the interface. The phonon transmission function can be obained in terms of the spectral heat fluxwhere is the Boltzmann parameter, and is the applied temperature bias. At low temperatures, inelastic effects caused by the anharmonicity were negligible and the spectral conductance matched the results obtained from the harmonic Green’s function method. Under this case, this transmission function is simply equal to , the number of propagating modes. The anharmonic phonon effect incorporated in NEMD render dependent on the sample length , which can be phenomenologically taken into account through the relation (Sääskilahti et al., 2015).where is the effective phonon mean free path (MFP).
Additionally, the phonon wave-packet (PWP) method based on wavelet transforms and MD simulation has also developed to calculate the phonon transmission coefficient. A detailed derivation process is referred to the articles by (Schelling and Phillpot, 2003). The advantage of the PWP method is that the phonon scattering process could be observed directly from the simulation, not just the transmission result as the AGF. The basic idea of PWP is to form a phonon wave packet so that it is localized in both real space and reciprocal space. To produce a wave packet at wave vector , the initial displacement from equilibrium of the ith atom in the nth unit cell along the a direction iswhere is the wave amplitude, is the mass of the ith atom, is the polarization vector, and are the positions of the packet center. It is worth noting that because of the very small amplitude of the incoming wave packets, only elastic scattering is captured, thus the frequencies of incoming and outgoing phonons are always the same. Inelastic scattering is very difficult to pick out the signal of wave-packet against the background of the thermal vibrations.
Atomistic Green’s Functions
As we know, phonons have wave nature. The AGF method can rigorously include wave effects on a discrete atomic lattice. Actually, this method was initially developed to handle quantum electron transport in nanostructures (Fan and Chen, 2010; Chen et al., 2018a; Deng et al., 2018; Zeng et al., 2018a; Chen et al., 2019a; Zhou et al., 2019; Fan et al., 2020). By making a few careful substitutions, the approach has been extended to study phonon transport in a wide variety of nanostructures (Xu et al., 2009; Peng and Chen, 2014; Zhou et al., 2017b; He et al., 2019) especially suitable for low-dimensional hetrostuctures such as Si/Ge (Tian et al., 2012b), graphene/h-BN (Peng et al., 2017), MoS2/metal (Yan et al., 2016b) interfaces and others (Ma et al., 2016). Here, we simply introduce this method from the numerical simulation point of view. The phonon Hamiltonian of entire system is consisting of three parts, namely Left (L), Center (C) and Right (R) regions, the matrix form can be written aswhere is a column vector consisting of all the atomic displacement in region, is nonlinear part of the interaction, is the coupling matrix between the central region and right lead. denotes coupled harmonic oscillators. Here is the force constant matrix obtained by empirical force field or the first principle calculations. Within the harmonic approximation (), the retarded/advanced Green’s function for central region is given byWhere is an identity matrix, and the is the self-energy of left and right leads. The retarded/advanced surface Green’s functions for the leads is defined as
The surface Green’s function usually solved recursively or using the decimation technique. According to Caroli formula, the phonon transmission coefficient as phonon frequency is represented aswhere , In terms of the Landauer formula under ballistic limit, the phonon thermal conductance is obtained from:
Here is the Bose–Einstein distribution, is the maximum phonon frequency. Note that the AGF method couldn’t directly provide the mode–resolved thermal condcutance. Even so, it is superior to the elastic wave continuum model, which is normally used at low temperatures (Peng et al., 2010; Xie et al., 2013; Xie et al., 2014b; Zhang et al., 2014b). Recently, a convenience technique developed by (Ong et al., 2016) has been widely used for the investigation of interface thermal transport, where the phonon reflection, absorption and transmission could be obtained, as described next. The first step is defined a non-Hermitian matrix that acts on the eigenstates of left or right lead. Here, take the left lead for example:
Diagonalizing one can obtain the eigenvector matrix and the eigenvalue matrix that consist of the normalized eigenstates and eigenvalue. These matrices are needs for calculation of individual phonon transmission. Next, the transmission matrix is computed by:where and . Once these matrices are known, the transmission coefficient of individual phonon mode on the left or right leads is given by:
By summing over the subscript m and n, the total phonon transmission in Eq. 9 can be obtained. If there is significant nonlinear effect in the central scattering region, the third-order force constant can be consierded as a perturbation of harmonic system. The three-phonon scattering self-energy is computed by using the following relations:where represents the Hilbert transform, is the third-order force constant. Finally, the effective transmission function is computed by:
Boltzmann Transport Equation
From the particle-based picture perspective, phonon transport is governed by the BTE, where phonons obey the Bose-Einstein distribution function at thermal equilibrium but deviating from under a temperature gradient. Here indicates the phonon mode of wave vector and phonon polarization . The deviating distribution function can be obtained from the BTE:where is the group velocity of phonon mode. Assuming a small enough in most practical situations, can be expanded to in first order so that , where depends on linearly. If only the scattering mechanism of two- and three-phonon processes are considered, Eq. 1 is written as:where and is the relaxation time of mode , which corresponds to the relaxation time approximation (RTA).where denotes absorption processes that the combined energy of two incident phonons leads to one new phonon (), and describes emission processes in which the energy of one incident phonon is split among two phonons (). A detailed descriptions of and is available in (Li et al., 2014b). Eq. 2 is numerically solved by iterating from a zero-th order approximation start of . The iterative procedure has a great influence on diamond, graphene and carbon nantubes (CNTs) in which the normal processes play a significant role in the phonon-phonon scattering (Lindsay et al., 2009; Lindsay et al., 2010; Lee et al., 2015a; Lee and Lindsay, 2017). However, for Si, Ge and amorphous materials with strong Umklapp scattering (Ward and Broido, 2010; Zhou et al., 2020a) iterating to convergence typically has little effect on the obtained thermal conductivity. Apart from phonon-phonon scattering, there may be various other phonon scattering processes, including phonon-defect (), phonon-impurity (), and phonon-boundary (). Based on Matthiessen’s rule where various scattering mechanisms are presumed to be independent of each other, the total relaxation time is given by . By the above solution, the lattice thermal conductivity can be obtainedhere is the specific heat of the mode . At present, the BTE in combination with first-principles calculations becomes popular to predict the thermal conductivity of various 2D materials, showing reasonable agreement with experimental measurements (Lindsay and Broido, 2011; Fugallo et al., 2014; Gu and Yang, 2014; Lindsay et al., 2014). It is writhing pointed out that the over-simplified of phonon scattering processes and a finite number of points causes obvious ambiguity in the calculation of thermal conductivity. Another point worth noting is that the coupling effect of various scattering mechanisms play important roles in phonon transport.
Experimental Measurement
The measurement of thermal conductivity or interfacial thermal conductance (ITC) in 2D nanostructures is helpful to grasp the law of phonon transport. At present, the optothermal Raman and microbridge techniques are the most common method for measuring the thermal properties of 2D materials. In addition the 3ω and time-domain thermoreflectance (TDTR) techniques have been adopted to obtain ITC. Because of space limitation, we briefly introduce some measurement results for 2D materials and discuss the measurement errors.
With the aid of optothermal Raman technique, the thermal conductivity measurement of single-layer graphene was reported by (Balandin et al., 2008), and the obtained value ranges from 2,000 to 5,000 Wm−1 K−1 at room temperature. Analogously, the thermal conductivity of bilayer graphene was also identified (Li et al., 2014a), coinciding with the theoretical predictions (Nika et al., 2014). Furthermore (Malekpour et al., 2016), reported that the room-temperature thermal conductivity of defective graphene reduces from 1800 to 400 Wm−1K−1 when the defect density induced by electron beam irradiation slightly increases from 0.0005 to 0.005%. (Chen et al., 2012) measured the thermal conductivity of graphene with different 13C isotope concentrations, and revealed that 50% concentration results in a 50% reduction in thermal conductivity. Except graphene, this technique has also applied to black phosphorous (Luo et al., 2015) MoS2 and MoSe2 (Zhang et al., 2015c), In a word, although the uncertainty of the optical absorption coefficient, weak temperature-dependence of Raman peak shift and uncertainty of contact thermal resistance will affect the measurement accuracy, the optothermal Raman measurement is still a simple and effective method to measure the thermal conductivity of 2D structures.
The microbridge technique was initially used to measure the thermal conductivity of 1D nanostructures (Li et al., 2003), and then was applied in 2D materials (Pettes et al., 2011). For example, (Xu et al., 2014b) experimentally found the logarithmic dependence of thermal conductivity on the sample length L () in suspended graphene. Moreover, (Jang et al., 2013) reported that there is positive correlation between thermal conductivity and the number of layers in few-layer graphene, and the value of ten-layer graphene is over 1,000 Wm−1 K−1. On the other hand, the optothermal Raman measurement indicated that the thermal conductivity of 4-layer MoS2 is 44–50 W m−1 K−1 at 300 K. (Jo et al., 2014). Analogously, (Lee et al., 2015b) verified that the thermal conductivity anisotropic ratio approaches two at 300 K. Note that the presence of contact resistance, polymer residue and inconsistency of sample qualities limit the precision of measurement.
The ITC measurement of in-plane heterostructures is quite challenging compared to that of vdW interfaces. There exist far fewer measurements on the ITC. Usually, the measured ITC of graphene/substrate interface is in the range from 20 to 100 MWm−2 K−1 at 300 K. Taking this value as reference, the measured ITC MoS2/substrate interface is considerably smaller. Using the TDTR-based approach, (Lee et al., 2015b) reported that the ITC of MoS2/metal thin-films interface is ∼26 MW m−2 K−1. Moreover, the optothermal Raman measurement indicated that the ITC of MoS2/SiO2 interface is about 1.94 MW m−2 K−1 (Taube et al., 2015). Analogously, (Behranginia et al., 2018) reported that the ITC between WSe2/SiO2 interface ranges from 10 to 32 MW m−2 K−1 via analyzing the temperature-dependence of the low-frequency shear mode. (Chen et al., 2014) measured the ITC of vertically stacked graphene/h-BN interface and found that the value reaches 7.4 MW/m2K at room temperature, which stems from the remaining organic residues at the interface.
Thermal Transport in Some Typical 2D Materials
Graphene
It is well known that single-layer graphene possesses ultra-high room-temperature thermal conductivity. The results of thermal conductivity measurements for graphene have been discussed in Molecular Dynamics Simulations. Note that the measured values are closely related to experimental approaches, external condition (strain, substrate), sample size (length, layer-thickness) and quality (like point-defect, isotopic abundance, grain size). Hence, theoretical simulations could provide deeper insight into concerning the thermal property of graphene. Using the PBT method combined in conjunction with the force-constant model, (Lindsay et al., 2010) claimed that the flexural acoustic (ZA) phonons contribute 70–80% of heat conduction in graphene, which is ascribed to the existence of selection rules in phonon scattering processes. The phenomenon has also been confirmed by others research groups based on similar theoretical approach (Fugallo et al., 2014; Lindsay et al., 2014). Using the BTE method and Raman technique, (Seol et al., 2010) reported that the thermal transport contribution from ZA modes for graphene is as large as 77%, and this ratio is decreased by the substrate or organic residues (Pettes et al., 2011). Recently, using the iterative BTE method, (Lee et al., 2015a) found that there is an obvious hydrodynamic phonon transport in graphene over a wide temperature range due to the presence of strong N scattering processes. Subsequently, (Machida et al., 2020) experimentally showed an extremely high thermal conductivity in the very thin graphite samples (8.5-micrometer-thick) at higher temperature, which extends the hydrodynamic regime from cryogenic to room temperatures. That is, there are three phonon transport behaviors, namely ballistic, diffusive and hydrodynamic transport. Furthermore, (Al and Farías, 2016) outlined the phonon dynamics in supported graphene, and revealed the physical picture of phonon-substrate scattering. (Nika and Balandin, 2017) minutely discussed the mode-resolved thermal conductivity in graphene. In addition, our previous work highlighted the manipulation of thermal transport in graphene from the viewpoint of the wave/particle nature of phonons (Chen and Chen, 2020a). (Dollfus et al., 2015) summarized experimental and theoretical results about the thermoelectric properties of graphene and its derivatives, and described how the poor thermoelectric performance of graphene achieves high thermoelectric figures by tuning on a wide range both the thermal and electronic transport based on different strategies of nanostructuring and bandgap engineering. Those studies on thermal transport of graphene have stimulated great interest in exploring other emerging 2D materials, as shown in Figure 1.
FIGURE 1
Hexagonal Boron Nitride
Monolayer h-BN is also called “white graphene” because of their similar honeycomb structure. Unlike zero band-gap graphene, h-BN with 5.8 eV band-gap has been viewed as the top ranked candidate for the dielectric material in electronic devices or the new nonmetal catalyst in the field of photocatalysis (Zhao et al., 2019). It was reported that the electron mobility of graphene electronic devices supported on few-layer h-BN far outweighs that supported on silicon dioxide (SiO2) because of the atomically smooth surface of h-BN (Serrano et al., 2007). In addition, h-BN also holds high chemical stability, strong mechanical strength and good thermal property. The basal-plane thermal conductivity of high-quality bulk h-BN sample could reach up to 390 W m−1 K−1 at 300 K (Geick et al., 1966), which shows great potential in current-generation dielectric materials. For sing-layer h-BN, the BTE calculations indicated that the system’s thermal conductivity is about 600 Wm−1 K−1 (Lindsay and Broido, 2011), but the EMD results showed that the value is 400 Wm−1 K−1 (Sevik et al., 2011).
Experimentally, because the intensity of Raman peaks in h-BN is very weak, the optothermal Raman method may not be applicable for h-BN. (Zhou et al., 2014) found only 5% optical absorption for nine-layer h-BN in optothermal Raman measurment, and the obtained thermal conductivity was 227–280 m−1 K−1. Using the microbridge method, (Jo et al., 2013) reported that the thermal conductivity of an 11-layer sample reaches ∼360 W m−1 K−1 at 300 K, approaching the value of bulk h-BN (Lindsay and Broido, 2011). Actually, the accurate measurement of thermal properties of 2D heterostructures still has a long way to go.
Silicene
Silicene possesses a buckled honeycomb atomic arrangement of Si atoms. There is a spin-orbit coupling of 1.55 meV owing to the breaking of mirror symmetry in this bucked structure. In addition, silicene can be easily integrated into electronic devices (Liu et al., 2013a). From the view of the heat dissipation of electronic devices, the influence of buckled structure on phonon transport is well worth exploring. Based on the BTE method, (Gu and Yang, 2015) demonstrated the thermal conductivity of silicene is only a fifth of that of bulk silicon (140 W m−1 K−1). Moreover, it was also found the contribution of ZA modes to the thermal conductivity of silicene is about 7.5%. Similarly, (Xie et al., 2014a) found the relaxation times of acoustic phonon modes with long-wavelength are now close to zero at , accounting for the low thermal conductivity for silicene.
Based on different MD simulations, the calculated thermal conductivity for silicene ranged from 5 to 50 W m−1 K−1, which could be ascribed to the adopted empirical potentials including the Tersoff potential, embedded-atom method (MEAM) potential, optimized Stillinger-Weber potential (Zhang et al., 2014a; Liu et al., 2014b). Compared with planar 2D materials, the structure buckling has important effect on the phonon transport in silicene via providing additional scattering or breaking the symmetry selection rule. As a result, there is an unusual thermal response to the applied tensile strain in silicene: a small strain can enhance the thermal conductivity while further increase will cause the decrease of thermal conductivity. That is differentiated from the monotone decreasing with tensile strain for graphene (Hu et al., 2013). This can be explained as follows: for a small strain, the bucked configuration becomes flat because of bond rotation, rendering the enhancement of in-plane stiffness and facilitating phonon transport; under large strain, the Si-Si bonds are significantly elongated and the in-plane stiffness weakens, lowering the thermal conductivity. However, the experimental study of thermal transport in silicene hasn’t been reported by now.
Transition Metal Dichalcogenides
Unlike graphene and h-BN with single atomic layer, a transition metal layer is sandwiched between two chalcogenide atomic layers for TMDs. In general, these TMDs possess a finite bandgap, and are suitable for many electronic and photonic applications (Zhao et al., 2020a). Owing to the difference of lattice structure between TMDs and graphene, TMDs demonstrates different thermal transport properties (Zhao et al., 2020b). Clearly, the strong covalent bonds between the carbon atoms result in a high thermal conductivity of graphene. By comparison, the transition metal-chalcogen bond is weak because of its special sandwich structure, raising the particularity of phonon properties. Besides the crystal structure, the atomic masses are another factor affecting the thermal conductivity of TMDs. A large mass difference between transition metal and chalcogen atoms often induces substantial phonon-gap (Jiang, 2014). For instance, the phonon gap in WS2 achieves 3.3 THz, and some three-phonon processes (like acoustic + acoustic optical) are strongly suppressed owing to the energy conservation requirement; while the phonon gap is only 1.4 THz in MoS2, far less than the frequency range (0∼6.9 THz) of acoustic branches (Gu and Yang, 2014). As a result, long phonon relaxation time is observed in WS2, and the thermal conductivity of WS2 (142 Wm−1 K−1) is significantly greater than that of MoS2 (103 Wm−1 K−1).
MoS2 and MoSe2 as the most typical TMDs have been experimentally fabricated (Keyshar et al., 2017), and their thermal transport properties have also been widely explored. (Wei et al., 2014) claimed that the thermal conductivity of MoS2 is about 26 Wm−1 K−1 at room temperature using BTE method under the RTA assumption. Analogously, (Gu and Yang, 2014) reported the thermal conductivity of MoSe2 is 54 Wm−1 K−1 under the iterative solution, which is over 20% than that under the RTA assumption. Furthermore, the relative contribution of different phonon polarization branches to thermal conductivity resembles the case in silicene rather than graphene because MoS2 is not an individual atomic plane. (Cai et al., 2014) predicted a room-temperature thermal conductivity of 23.2W Wm−1K−1 for MoS2 by solving the AGF. Based on NEMD simulations, (Hong et al., 2016b) showed that the thermal conductivities of MoS2 and MoSe2 are 110.43 and 43.88 Wm−1 K−1, respectively. Lately, 2D MXenes, namely hexagonal transition metal carbides/nitrides, are highly attractive because of their promising technological applications such as energy storage, thermoelectric conversion and flexible devices (Ma et al., 2014; Sarikurt et al., 2018). Clearly, the performance of aforementioned devices strongly depends on the thermal transport properties of constituent material. As we know, low thermal conductivity and excellent conductor are desired in thermoelectric devices, while micro-nano electronic devices require the material with high thermal conductivity. Some previous studies demonstrated that the heat dissipation is a major challenge for MXene-based devices owing to their large aspect ratios and heterogeneous structures. For example, (Hemmat et al., 2019) found the ITC of Ti3C2Tz/SiO2 interface is the range of 10–27 MW m−2 K−1 under different annealing temperatures, and a similar range of ITC was obtained in Ti3C2/SiO2 interface with or without AlOx encapsulation (Yasaei et al., 2018). On the other hand, the in-plane thermal conductivity of MXenes has also been studied. For instance, (Guo et al., 2018) demonstrated that the thermal transport performance in luorine-terminated Ti2CTz exceeds that in oxygen-terminated one, but is comparable to that of bare Ti2C using BTE method. Analogously, (Gholivand et al., 2019) reported that compared to oxygen-terminated Ti3C2Tz, the fluorine-terminated structure has remarkable advantage in heat dissipation. In another theoretical study, (Sarikurt et al., 2018) investigated the influences of transition metals and surface functionalization on the thermal conductivity of M2CO2, and the obtained results showed that their thermal conductivities range from 15-60 Wm−1 K−1.
Experimentally, the optothermal Raman method was used to measure the room-temperature thermal conductivities of 1-layer and 11-layer MoS2, and the obtained values were ∼34.5 Wm−1 K−1 (Liu et al., 2018a) and ∼52 Wm−1 K−1 (Sahoo et al., 2013), respectively, which exceeds the results calculated from the NEMD method (Wang et al., 2016). Meanwhile, the measured value for 4-layer MoS2 was in the range of 44∼50 Wm−1 K−1 by thermal bridge method (Jo et al., 2014). In addition, the temperature-dependent thermal conductivity of substrate-supported MoS2 was extracted from the optothermal Raman method, and the measured room-temperature thermal conductivity was 62.2 Wm−1 K−1, but reduced to 7.45 W Wm−1 K−1 at 450 K (Taube et al., 2015). Apart from MoS2, (Peimyoo et al., 2015) used the optothermal Raman method to determine the thermal conductivity of monolayer and bilayer CVD-grown WS2. The measured values are 32 and 53 Wm−1 K−1 for monolayer and bilayer WS2 at 300 K, respectively. Analogously, (Zhang et al., 2015c) reported that the room-temperature thermal conductivity of monolayer MoSe2 is 59 ± 18 Wm−1 K−1. The similar phonon and Raman scattering in different kinds of TMDs have been discussed earlier. It is believed that other than phonon frequency due to the different mass between Mo and W atom, the common features of 1 L, 2 L and bulk structures in both MoS2 and WS2 should be similar. (Yan et al., 2013) demonstrated a low thermal conductivity of 9 W m−1 K−1 for TaSe2 with a thickness of 45 nm. The experimental discrepancies were primarily due to the sample quality like surface disorders, and the measurement error. Most recently, (Liu and Li, 2018c) experimentally observed that the thermal conductivity of Ti3C2Tz to be 55.88 Wm−1 K−1, Note that neither explicit configuration nor surface termination atoms are clarified in this study.
Phosphorene
Phosphorene as a 2D puckered honeycomb structure possesses high carrier mobility and large direct band-gap (∼1.5 eV) (Liu et al., 2014c), indicating abundant nanoelectronic applications. In addition, phosphorene has shown promising application in thermoelectrics devices (Fei et al., 2014). Actually, either electrical or thermoelectrical devices strong depend on the thermal transport properties of phosphorene. Hence, the investigation of the thermal conductivity of such 2D material is necessary.
Theoretically, the thermal conductivity of phosphorene has been explored using different methods including EMD (Xu et al., 2015), BTE under RAT (Zhu et al., 2014a) or iterative rule (Qin et al., 2016). For example, (Qin et al., 2016) showed that the simulated thermal conductivity along zigzag (ZZ) and armchair (AC) directions are respectively 15.33 Wm−1 K−1 and 4.59 Wm−1 K−1 using the iterative BTE method. The low thermal conductivity could be understood from the softening transverse optical (TO) phonons resulted from the long-ranged interaction in phosphorene structure. Combined the spectral energy density (SED) technique with EMD simulations, the ZZ and AC -direction thermal conductivity were predicted to be 152.7 and 33 Wm−1 K−1 in phosphorene, respectively (Xu et al., 2015). (Zhu et al., 2014a) used the RTA to predict the thermal conductivity of black phosphorene of 84 Wm−1 K−1 (24 Wm−1 K−1) in the ZZ (AC) direction (armchair) at 300 K. Although the results reported are discrepant from different research groups, the thermal transport anisotropy behaviors have been clearly demonstrated for all techniques. To understand the phenomenon, the thermal conductivity contribution from ZA modes has also been calculated, which corresponds to 16% (8%) for the case of the ZZ (AC) direction (Qin et al., 2016). That is analogous to silicene (7.5%) but different from graphene (75%), because the hinge-like structure of phosphorene broken the selection rule in three phonon processes. Recently, (Zhao et al., 2018b) thought that the large anisotropic thermal transport stems from phonon group velocity rather than relaxation time. To be specific, the speed of sound velocity is positively correlated with the Young’s modulus E, and the group velocity of low-frequency phonons could be treated as the speed of sound. Meanwhile, the three-point bending method was employed to measure the Young’s modulus of black phosphorene along crystallographic directions. It was found that the thermal conductivity anisotropy ratio is closely related to the Young’s modulus one. Moreover, further enhancement of the anisotropy ratio can be realized by introducing phononic crystal (Xu and Zhang, 2016a) or loading tensile strain (Zhu et al., 2018). The latter regulation has also been experimentally proved, where Raman sensitive peak , and modes have opposite responses when the tensile strain was respectively applied to ZZ- and AC-directions.
Considering the limitations of synthesis technique, only few experiments were performed to extract the thermal conductivity of different thicknesses phosphorene films. Based on the measurement of micro-Raman spectroscopy, (Luo et al., 2015) demonstrated that the thermal conductivity of a 15 nm thick sample is ∼40 Wm−1 K−1 (20 Wm−1 K−1) along AC (ZZ) direction and decreases with the decrease of film thick. That is, the thermal conductivity anisotropic level has negative relation with film thickness, as shown in Figure 1. Another example, (Jang et al., 2015) adopted the TDTR method to investigate the thermal transport properties of silicon-supported thin films, and claimed that the thermal conductivity anisotropic ratio is about 2.5. Using the thermal bridge method, (Lee et al., 2015b) showed that when the thickness of phosphorous flakes reduces from 300 to 50 nm, the ZZ (AC) direction thermal conductivity drops to 12 from 27 Wm−1 K−1 (5 from 15 Wm−1 K−1). The first-principle calculation revealed the anisotropic thermal transport stems from orientation-dependent phonon dispersion and phonon–phonon scattering based on first-principle calculation. By comparison, the out-of-plane thermal transport is independent of the thickness of phosphorene film. Interestingly, (Su and Zhang, 2015) performed the optothermal Raman method to evaluate the substrate effect on the thermal conductivity phosphorene film, and the measured results indicated that the introduction of SiO2/Si substrate can enhance the thermal conductivity by more than one fold. The novel phenomenon seems to contradict the generally held notion held notion that the substrates have negative effect on the thermal conductivity of 2D materials.
Borophenes
Borophene is quite unique because there are rich hexagonal holes distributions in the triangular lattice such as the structure of borophene is mixed by triangle lattice and hexagonal lattice. At present borophene sheets with different phases and borophene nanoribbons have been experimentally synthesized on various substrates using the molecular beam epitaxy (MBE) approaches. For example, (Mannix et al., 2015) fabricated borophene sheets on the Ag (111) surface and found the planar islands with single atoms thickness occur at substrate temperature 570 and 650 K. Combined with theoretical calculations, the two phases of borophene can be identified as and sheets. Moreover, (Li et al., 2018d) chose Al (111) substrate to grow boron because more charge can be transferred from Al (111) to boron atoms, which help enhance the stability of honeycomb borophene. The experimental results are very exciting, in which the perfect honeycomb-shape islands with single atom thickness form at substrate temperature 500 K and the lattice constant (0.29 nm) agrees well with the calculated value of freestanding honeycomb borophene. In others words, the type and temperature of substrate and is important for the formation of borophene morphology like 3D clusters or planar borophene. The successful preparation of several borophene structures provides opportunities for further exploring their electronic properties and device applications.
It is well known that light atomic mass and strong interatomic interaction are the crucial factors for the high thermal conductance of a nanostructure. Compared to carbon atom, the atomic mass of boron is light, and the Young’s module of borophene is close to that of graphene (Li et al., 2018a), indicating that borophene probably possesses superior thermal conductivity as high as that of graphene. Using AGF method, (Zhou et al., 2017a) found the phonon thermal conductance of borophene along the AC direction is as high as 7.8 nWK−1 at 300 K n the ballistic transport region, which is almost twice that of graphene. Moreover, using the ReaxFF potential, (Mortazavi et al., 2017) predicted the thermal conductivity of borophene to be about 147 and 75.9 Wm−1 K−1 along AC and ZZ directions, while it is about 233.3 and 102.5 Wm−1 K−1 using the SW potential. Here, the low thermal conductivity for borophene was explained by the following two reasons. On the one hand, small extension of acoustic dispersions and out-of-plane symmetry breaking leaded to stronger anharmonic scattering (Sun et al., 2016). On the other hand, the N-processes have very little contribution to thermal conductivity. All of that causes short phonon lifetime, and lowers thermal conductivity.
Thermal Transport in 2D Heterostructures
In-Plane Heterostructures
2D heterostructures provide possibilities for high-performance optoelectronic devices. Modern advanced fabrication techniques like molecular beam epitaxy and epitaxial growth techniques enable the in-plane heterostructures with atomically sharp interfaces (Liu et al., 2013b; Gao et al., 2013; Li et al., 2015). Nevertheless, because of lattice matching constraints, only certain combinations of materials can enable continued heteroepitaxial growth. On the other hand, the presence of interfaces in heterostructures also induces some novel properties that are distinct from the parent materials. In particular, the lattice mismatch-induced strain would change the force environment of atoms, which cause the deformation of 2D crystalline structures and the variation of phonon spectra. At present, the ITC measurement of in-plane heterostructures is quite challenging compared to that of vdW interfaces. Hence, theoretical simulations provide a convenient way to understand interfacial thermal transport of in-plane heterostructures.
Benefitting from the mall lattice mismatch between grpahene and h-BN, their in-plane heterostructures as a typical case have been successfully synthesized (Gao et al., 2013). The simulated results of EMD method showed that the thermal conductivity along the vertical-interface direction is limited by the less conductive component (h-BN) and is insensitive to the details of the interface (Song and Medhekar, 2013), while that along the parallel-interface direction stabilises a constant which approaches the average of the parent materials, and gradually declines with increasing BN composition. The vibrational eigen-mode analysis reveals that increasing the BN composition would cause intermediate-frequency phonons to be localised in graphene domain (Chen et al., 2016b). (Hong et al., 2016c) studied the effect of sample length and ambient temperature on the ITC in graphen/h-BN in-plane heterostructures via NEMD method, and showed that the ITC varies from 1.9 to 4.5 GW/m2K as the sample length in the range from 20∼100 nm and rising temperature could obviously enhance the ITC. Based on the same method and analogous models, (Liu et al., 2019a) found that increasing the number of atomic row can elevate (lower) the ITC in the zigzag (armchair) interface topography. To understand the mechanism of phonon transmission/scattering at interface, (Feng et al., 2017) calculated the temperatures of phonons under the framework of NEMD simulations. Taking graphene/h-BN planar interfaces as example, the temperature jump of the ZA mode is only 10∼20% of the overall temperature jump because there is considerable overlap between the acoustic branches of graphene and h-BN, while that of the ZO mode is as high as 220% of the overall one, stemming from the negligible coupling at the interface, as shown in Figure 2. Moreover, when interface mixing is introduced, the contribution ratio of ZA mode to the overall temperature jump increases from 19 to 41% and that of the in-plane modes decrease, indicating that a certain number of heat transfer channels shifts from the ZA mode to the TA/LA modes. Using AGF method, (Ong et al., 2016) demonstrated that alternating the bond orientation and type at the interface significantly affect the scattering and transmission of higher frequency modes. In the case of the zigzag-bond interface, the bonding type is able to dramatically alter the preferred orientation of the transverse optical (TO) phonons.
FIGURE 2
During the crystal preparation process, a series of defects like single-vacancy (SV), stone−wales (SW) ineluctably exist in heterostructures. (Sutter et al., 2011) experimentally prepared graphene/h-BN in-plane heterostructures via a sequential CVD growth method, and pointed out that the interface of prepared heterostrucutres will be intermixing when the growth temperature becomes high. Based on NEMD simulations, (Liu et al., 2016) showed that the topological defect near the graphene/h-BN interface could promote the interfacial thermal transport, which may be attributable to the localization of the stress fields. Analogously, (Li et al., 2018b) systematically studied the influences of SV and SW interfacial defects on the ITC in the in-plane heterostructures. The results indicated that the ITC quickly decreases with SV concentration but slowly with SW defects, and ultimately remain unchanged for large concentration. The phonon spectra and stress distribution analysis revealed that the regression of in-plane modes induced by SV defects results in the reduced thermal conductivity. Moreover, (Gao and Xu, 2017) reported that the existence of SW interfacial defects can realize the regulation of phonon transport direction, being equivalent to applying mechanical tensile strain. More recently, (Ni et al., 2019) reported the dependence of thermal conductivity on interface composition diffusion in graphene/h-BN in-plane heterostructures, and demonstrated that the system’s thermal conductivity is closely related to the interface composition diffusion length because of the Anderson localization of phonons. (Liang et al., 2020) found that the ITC decreases with the increase of layer number n and saturates at n = 3 in multilayer in-plane graphen/h-BN heterostructures. Even so, the asymptotic value still surpasses that in monolayer heterostructures.
Graphene/h-BN superlattices as special hetrostructures were forecasted to own numerous excellent properties like good thermoelectric performance (Yokomizo and Nakamura, 2013) and opto-electronics properties (Liu et al., 2013b). Using AGF method, (Jiang et al., 2011b) reported that there is a minimum in the period length dependence of the thermal conductance of graphene/h-BN superlattice. In general, the thermal transport characteristic of a superlattice indicates that thermal conductivity decreases first, and then increases with the increase of period length. In principle, to the left of the minimum value, heat conduction mostly relies on the coherent phonon transport. In turn, it is dominated by coherent phonons induced by interfacial modulation, which has been verified experimentally in 3D superlattices (Ravichandran et al., 2014) and theoretically in other 2D superlattices (Liu et al., 2014b). It is worth noting that the thermal conductivity of superlattices often decreases remarkably compared to that of parent materials, suggesting promising potential thermoelectric applications. Furthermore, (Zhu and Ertekin, 2014b) proposed that the initial reduction of thermal conductivity in graphene/h-BN superlattices is principally accounted for the effect of harmonic wave like phonon band gap and flat phonon branches. Subsequently, weakening interface scattering is conducive to thermal transport. Moreover, they also found that the critical period length is insensitive to the sample length. In addition, (Da Silva et al., 2016) carried out normal mode decomposition (NMD) approach to compute the thermal conductivity of graphene/h-BN superlattices with short , and showed that the phonons with frequency range 0∼25 THz accounts for over 90% of the total thermal conductivity regardless of the period length. Moreover, the reducing thermal conductivity with increasing and was attributed to the decrease of phonon group velocities rather than relaxation times. Recently, (Chen et al., 2016c) found that the critical period length decreases at higher temperatures in the superlattices, suggesting the suppression of phonon wave effect. Interestingly, for certain period lengths, the system’s thermal conductivity at 200 K is very close to that at 300 K, contrary to previous studies that high temperature normally leads to the decrease of thermal conductivity, as shown in Figure 3. The reason for this unusual phenomenon was identified as the emergence of strong phonon wave interference under lower temperature. Thereafter, (Wang et al., 2017c) observed that the temperature profile behaves wave-like fluctuation in pristine/nitrogenated holey graphene superlattice at 100 K, which is also believed to be the result of the phonon wave interference. (Mu et al., 2015) reported that the existence of smooth interface will disrupts the coherence of phonons in supperlattice structure, and then lowers its thermal conductivity. Furthermore, the phonon transport properties of other 2D superlattices including silicene/germanene superlattice (Wang et al., 2017b) and MoS2/MoSe2 superlattice (Ding et al., 2018a) have also been studied using different simulation methods. These previous studies provide valuable insights on the phonon transport in in-plane heterostructures.
FIGURE 3
Except for graphene/h-BN in-plane heterostructures, some previous research also focused on the interfacial thermal transport across other 2D heterojunctions. Combined the first-principles with NEMD methods, (Liu et al., 2017a) indicated that the ITC of graphene/MoS2 in-plane heterostructure is higher than covalently bonded graphene/metal interfaces. Here, every Mo-C bond could act as a heat transfer channel. (Liu et al., 2014a) explored the dependence of ITC on temperature and strain in graphene/silicene in-plane heterostructure, and concluded that high temperature and small strain is beneficial to improve the ITC. In this work, it is worth noting that when the heat flux imposed into the system is lower than 42 GW m−2, the value ITC almost varies with increasing the heat flux, corresponding to the Fourier heat conduction; while it exceeds the critical value, a low-frequency kinetic wave is stimulated, constructing an extra channel for the non-Fourier heat transfer. (Qin et al., 2019) found that the ITC of MoS2/WSe2 in-plane heterostructure is about 10% of that of graphene/h-BN heterojunction because of strong phonon scatterings at the interface, which could be useful for thermoelectric applications. (Liu et al., 2018a) claimed that the ITC of phosphorene-based coplanar heterostructures doesn’t present obvious anisotropy as the thermal conductivity of phosphorene because the interfacial anharmonicity results in a “mixing” effect for phonon transport. Moreover, the dependence of ITC on strain is very sensitive to the interface morphology, which can be understood by the analysis of stress distribution and phonon spectra overlap.
van der Waals Heterostructures
Benefitting from the progress of ingenious experimental methods, the synthesis of vdW heterostructures via self-assembly technology has achieved an unprecedented degree of control and accuracy, and provides new insights to explore structure-function relationship (Frisenda et al., 2018) and to assemble complex devices (Abidi et al., 2018; Ding et al., 2018b) which are able to be tuned by stacking sequence, orientation, interlayer coupling, and external field. For example, graphene/hBN vdW heterostructures were prepared via an all-dry poly (methyl methacrylate) transfer process using polyvinylalcohol as the water-soluble sacrificial layer, which can avoid the polymer residual problem and prevent contamination from chemical solutions (Tien et al., 2016). Another example is that monolayer WS2 and MoS2 can be directly deposited on h-BN by CVD, which maintain their intrinsic bandgap of 2.01 and 1.85 eV (Fu et al., 2016), indicating a high quality of heterostructures with uniform distribution and a clean interface being developed. Thermal properties of vdW heterostructures is of great interest because it is highly relevant to the functionality, performance and reliability of such heterostructures-enabled devices.
For bilayer graphene/h-BN heterostructures, there are periodic moire patterns on graphene layer. Three typical stacking configurations including AA, AB, AB′ stacking have been verified in experimental results. Moreover, different moire patterns may also affect the phonon transport in such heterostructures. In terms of the analysis of phonon dispersion, (Jung et al., 2012) found that the AB stacking is superior to the AA stacking for graphene/h-BN vdW heterostructures. (Zhang et al., 2017c) carried out the NEMD simulations to calculate the thermal conductivity of h-BN-supported graphene, and demonstrated that the system’s thermal conductivity is reduced only 23% compared to the suspended case. That is because the force environment of graphene is barely affected by h-BN substrate. Additionally, they also discovered that the thermal transport properties strongly depend on the stacking ways of vdW heterostructures. Based on AGF method, (Yan et al., 2016a) showed that when the interlayer spacing between grapheme and h-BN layers is small, the perpendicular heat conduction is determined by in-plane acoustic modes. Another important finding in this work is that the ITC could be enhanced by more than 50% from C-N matched to C-B matched interface. Using the same method, (Wang et al., 2013b) found that the thermal conductance of graphene coupled with h-BN greatly reduces at low temperatures but can be still preserved as high as 96% that of graphene at room temperature. This could be explained by the fact that the h-BN coupling affects the out-of-plane vibrations, linearizing the flexural phonon modes and suppressing the transmittance. However, it has a minor effect on the in-plane vibrations, leaving the TA and LA phonons unchanged. Furthermore, the effect of different layers of h-BN (N) on the thermal conductance of graphene has also been studied, and the results showed that the presence of multi-layers h-BN () hardly affect the heat conduction performance of graphene, indicating that the design of the graphene/h-BN vdW heterostuctures is robust from the standpoint of heat dissipation. Analogously, (Pak and Hwang, 2016) predicted that the thermal conductivity of graphene changes very little after introducing h-BN substrate. The reason was that the strong interlayer coupling could minimize the internal stress fluctuations of graphene, as shown in Figure 4. Moreover, the ITC at gaphene/h-BN interface was computed to be 4.5 M Wm−2 K−1. Acoording to lattice dynamics theory, (Zou and Cao, 2017) thought that the phonon group velocity of graphene supported by h-BN have minor changes in comparison with suspend graphene because of weak interlayer coupling. Meanwhile, the lifetimes of ZA phonons had significant reduction because of the invalidity of the selection rule. Using a fast transient technique, (Zhang et al., 2015a) observed that hydrogenating graphene can improve the ITC of grapheme/h-BN vdW heterostructures by more than 70%. That is because the hydrogenation could enhance the coupling degree between in-plane and out-of-plane phonons, and widens phonon channels. Furthermore, the reported ITC monotonically decreased with temperature and interatomic bond strengths. (Chen et al., 2014) used the Raman spectroscopy technique to measure the ITC of graphene/h-BN vdW heterostructures and the obtained ITC was 7.4 M WK−1 m−2 at room temperature, which exceeds the value of graphene/SiC interface.
FIGURE 4
Apart from graphene/h-BN vdW heterostructure, the ITC of other vdW heterostructures have also been investigated. Based on MD simulations, (Pei et al., 2017) showed that constructing the graphene/phosphorene/graphene sandwiched structures can elevate the thermal stability and thermal conductivity of phosphorene. (Liang et al., 2018) found that the in-plane thermal transport in graphene/phosphorene vdW heterostructures exists obvious anisotropy and increasing the interlayer coupling could enhance the thermal conductivity of phosphorene layer because of phonons stiffening. With the aid of NEMD simulations, (Zhang et al., 2016) found defect, hydrogenation and compressive strain can increase the ITC of graphene/phosphorene/grapheme multilayer structures. By BTE method, the room-temperature thermal conductivity of MoS2/MoSe2 bilayer heterostructure was predicted to be 25.4 Wm−1 K−1 (Ma et al., 2019), which is intermediate between MoSe2 and MoS2. Moreover, it was also found that increasing temperature will decrease the in-plane thermal conductivity, but increase the ITC. (Gao et al., 2016) showed that the interlayer coupling can reduce the in-plane thermal conductivity of graphene layer in graphene/MoS2 bilayer heterostructure, but barely affects the thermal transport in MoS2 layer. Similar phenomenon has been reported in the h-BN/silicene bilayer heterostructure by (Cai, et al., 2016). More importantly, they also found that the electronic and phonon transport are dominated by the silicone and h-BN layers, repectively, which realizes the decupling of phonon and electron transport. (Hong et al., 2018) summarized systematically the ITC of different grapheme-based heterostructurs, and the obtained results range from 2.8 to 25 M WK−1 m−2 at temperature range from 100 to 700 K. In a word, the ITC can be modulated by temperature, contact pressure, hydrogenation and vacancy defects. Experimentally, the ITC measurement of vdW heterostructures has been reported in previous studies. Using the Raman optothermal technique, (Liu et al., 2017c) observed that the ITC of MoS2/h-BN vdW heterostructure (∼17 MW m−2 K−1) is a third of that of graphene/h-BN one (∼52 MW m−2 K−1). That is because the lighter mass of carbon atoms leads to larger frequency range that can be available for phonon transmission from the graphene to h-BN region. Another example, (Kim et al., 2018) used the space-resolved Raman spectroscopy technique to show that the doping level regulated by high electric field is beneficial to the spatial heat dissipation of graphene/h-BN devices, particularly close to the temperature of graphene breakdown.
Manipulation of Thermal Transport in 2D Materials and Their Heterostructres
The manipulation of thermal transport in 2D heterostructres is crucial for the design of micro-nano devices based on them. However, since thermal transport is carried by a large population of phonons spanning a frequency spectrum and existing diverse interactions, controlling phonon transport is actually a complex task. Basically, the phonon scattering processes could be enhanced by introducing structural defects like isotopes, vacancy or dislocation, substrates, chemisorption, and strain. Note that using the mentioned-above controlling method, only the phonons with short wavelengths gets scattered efficiently by above tuning methods, while some long-wavelength phonons are nearly not influenced. Recently, a new thinking for manipulating phonon transport on the base of the wave nature of phonons like interference or local resonance has drawn considerable attention. Besides, the pseudospins and topological effects of phonons as new quantum degrees of freedom might be utilized to manipulate phonons (Gao et al., 2018; Liu et al., 2019b; Xiong et al., 2018).
From the Perspective of Phonon-Particle Picture
During the preparation of 2D materials, point-defect, isotopes and grain boundaries, inevitably resides in samples. The presence of imperfections normally introduces phonon scattering and thus causes decreasing thermal conductivity. The effect of randomly distributed isotopes was studied in h-BN via BTE calculations (Lindsay and Broido, 2011), demonstrating that the reduced thermal conductivity is the result of high-frequency phonon modes scattering induced by mass difference. Using AGF method, (Mingo et al., 2010) found that the presence of isotope clustering could strengthen the cross section scattering. Analogously, (Morooka et al., 2008) demonstrated two interesting phonon transport behavior in GNRs with SW-defects: the first is that there is edge heat flux and the second is the appearance of circulating heat flux near SW defect. Subsequently, (Jiang et al., 2011a) showed that the probability of phonon transmission across the Si-doped region in graphene is very low, causing a significant reduction of thermal conductivity. (Huang et al., 2013) reported the influence of extended line-defects on the thermal conductance of GNRs, and showed that the modulation of the relative orientation between line-defect and transport direction can substantially change the system’s thermal conductance. With aid of the NEMD simulations, (Ding et al., 2015) indicated the thermal conductivity of MoS2 is reduced by over 60% for a 0.5% mono-Mo vacancy concentration. The analysis of SDE revealed that the reduction of thermal conductivity in defective MoS2 is mainly determined by the phonon relaxation time. Analogously, (Zhou and Chen, 2015) found that the presence of nanoholes in graphyne nanoribbons will result in the intensive localization of phonons. In particular, a thought of “nanodomains in 2D alloy” has been proposed (Gu and Yang, 2016), because alloying can strongly scatter the high-frequency phonons, as shown in Figure 5. Moreover, (Ding et al., 2016) found that an unexpected enhancement of the ITC is observed when vacancy defects are introduced at the interface layers. The defect-induced enhancement in thermal transport is attributed to the increased friction and thus enhanced phonon coupling at the interfaces. Using the suspended thermal bridge method, (Aiyiti et al., 2018) experimentally verified that the introduction of defect is an effective approach to lower the thermal conductivity of few layer MoS2.
FIGURE 5
Mechanical strain as a regular method has been widely applied to modulate the physical properties of 2D materials. In general, the thermal conductivities of 2D planar nanostructures like graphene or graphyne monotonically decrease with increasing the tensile strain because of phonon softening, which lowers the phonon group velocity (Wei et al., 2011). By contrast, in the buckled 2D materials like silicene or phosphorene, the dependence of thermal conductivity on the tensile strain shows a transition from the initial increase to subsequent reduction variation. For instance, based on MD simulations, (Fan et al., 2017) demonstrated that the loading of tensile strain can increase the relaxation times of ZA phonon modes in graphene. Using BTE method under RTA, (Kuang et al., 2016) reported that the applied tensile strain can either reduce or increase the thermal conductivity of graphene, which depends on the sample length. Apart from graphene, (Jiang et al., 2013) demonstrated that the thermal conductivity of MoS2 monotonically decreases with the increase of tensile strain. (Xie et al., 2016) performed the BTE simulations to study the controllable thermal conductivity in monolayer silicene by mechanical strain. Here, there was an up-and-down relation between the obtained thermal conductivity and the tensile strain, because a small tensile strain is able to partially recover the selection rule that only exists in 2D planar nanostructure. While the tensile strain becomes large, the appearance of phonon softening lowers the phonon group velocity. Moreover, using NEMD method, (Zhu and Ertekin, 2015) reported that the thermal conductivity of graphene/h-BN superlattice initially increases to a peak value, after which it then decreases with further strain. The non-monotonic curve was considered as a competition between in-plane softening and flexural stiffening of phonons. (Chen et al., 2017c) showed that as a compressive strain of 0.18 is applied in the graphene/black phosphorus vdW heterostructures, the ITC increase more than 15-fold compared to the strain-free value. The significant enhancement of ITC induced by compressive strain could be understood from the following two aspects. The first is that the applied normal pressure on the multilayer structures clearly strengths the interlayer coupling. The second is that the compressive strain also enhances the shear interaction between the interlayers and thus improves the thermal transport contribution from in-plane phonons. Experimentally, (Guo et al., 2020) reported the tensile strain dependence of the thermal conductivity of monolayer graphene using optothermal Raman method. The measured results demonstrated that under a 0.12% biaxial tensile strain, the system’s thermal conductivity drops approximately by 20% at 350 K, which verifies the previous theoretical literatures prediction (Varshney et al., 2018). Similar phenomenon has also been reported in polycrystalline multilayer graphene based on the microbridge measurement (Zeng et al., 2020c). Additionly, using TDTR method, (Meng et al., 2019) showed that a 9% cross-plane compressive strain can cause the 7-fold increase of the cross-plane thermal conductivity of multilayer MoS2, which is available for improving MoS2-based device performance and stability, as shown in Figure 6. The analysis of lattice dynamics and phonon spectroscopy revealed that the drastic increase of cross-plane thermal conductivity with pressure is mainly the result of strengthened interlayer forces and enhanced group velocity of longitudinal acoustic phonons, while the saturation of thermal conductivity above 15 GPa is associated with the combined effects from increasing group velocity and reduced phonon lifetimes. It is conceivable that the observed phenomena should occur in most vdW heterosturctues.
FIGURE 6
The third common approach to modulate the thermal conductivity of 2D heterostructres is chemical functionalization. Using SED technology, (Pei et al., 2011; Kim et al., 2012) found that the hydrogenation coverage dependence of thermal conductivity in functionalized graphene shows a U-shaped variation, similar to the composition dependence of thermal conductivity in Si/Ge alloy (Xie et al., 2020). (Mu et al., 2014) reported that a 5% oxidation coverage leads to a 90% reduction of thermal conductivity in graphene samples. It can be explained from the following two aspects: one is that the existence of the oxygen atoms adds phonon-impurity scattering; the other is that the reflection symmetry of pristine graphene plane is destructed owing to oxygen coverage, which enhances the scattering of flexural modes. In some cases, hydrogenation is beneficial for the thermal transport in 2D materials. According to the BTE calculated results, (Wu et al., 2016) demonstrated that the hydrogenations can enhance the phonon transport in penta-graphene (up to 76% increase), which ascribes to the wakening anharmonicity of C-C bond. Using EMD method, (Yang et al., 2018) found that when the coverage of the adsorbed groups exceeds 0.1, the thermal conductivity of functioned h-BN remains constant (about 100 Wm−1 K−1), which are still sufficiently thermally conductive as fillers of composites. (Hong et al., 2016a) showed that the ITR of graphene/phosphorene bilayer heterostructure increases monotonically with the hydrogen coverage in graphene layer. The maximum addition of 84.5% is observed at 12% hydrogen coverage. The enhanced thermal transport can be understood that an extra H–P heat dissipation channel is created in addition to that between C–P atoms and the absorbed H atoms on graphene behave as scattering centers, thereby enhancing the graphene’s lateral to flexural direction phonon coupling which indirectly strengthens the thermal transmission from graphene to phosphorene. (Liang et al., 2018) found that the in-plane thermal conductivity of phosphorene layer in graphene/phosphorene vdW heterostrucutres is significantly improved by hydrogenating the graphene layer.
When 2D materials or related heterostructres are integrated for device applications, their thermal transport properties have been drastically altered on account of the phonon-substrate scattering. The heat generated in device geometry may flow in the sample, and is dissipated directly into the substrate through their interfaces. (Sadeghi et al., 2013) experimentally observed that when the number of layers in SiO2–supported few-layer graphene achieves 34, the system’s thermal conductivity almost equated with the natural graphite, because there are long-wavelength out-of-plane phonons and partially diffusive interface-phonon scattering in thick sample. In general, the substrate will cause the enhancement of optical phonon scattering in 2D materials (Berciaud et al., 2009). In terms of SED analysis, (Qiu et al., 2012b) argued that the substrate effect is responsible for the reduction of the phonon relaxation time of ZA modes. Interestingly, (Su and Zhang, 2015) experimentally demonstrated that the SiO2 substrate could double the thermal conductivity of phosphorene film, which differs from the case in SiO2-supported CNTs (Ong et al., 2011). Based on EMD simulations, (Chen et al., 2016a) revealed that the degree of thermal transport anisotropy can be significantly enhanced in supported phosphorene. That was because the substrate effect significantly suppresses the ZA modes that dominate the thermal transport along armchair direction. (Li et al., 2018c) found that the ITC of graphene/h-BN in-plane heterostructrue is enhanced as it couples with a SiO2 substrate, which ascribes to the increase of thermal transport channels. Similar simulated results have been reported in asymmetric graphene/h-BN in-plane heterosturctures (Sandonas et al., 2017).
From the Perspective of Phonon-Wave Picture
Considering that phonons are waves of the atomic lattice, it is natural to think phonons interference, which could be used to regulate thermal transport (Chen et al., 2005; Xie et al., 2018a; Xie et al., 2018b; Zeng et al., 2020a). Actually, in periodic superlattice, phonons would be reflected by the periodic interfere as that the photons in photonic crystals. To produce the interference patterns, it should satisfy the following two conditions. First, the interface roughness is much less than the phonon wavelength, in which the phase of phonons could be preserved after undergoing interface scattering. However, the wavelengths of most phonons at room temperature are rather short in superlattice. Hence, atomically smooth interface are required for distinct phonon interference. Second, the spacing of two adjacent interfaces (namely period length) should also be relatively short so that propagating phonons will not be scattered before being reflected at the interfaces. For this, the coherent phonons consisting of the constructive interference of incident and reflected phonons could preserve their phases across multiple interfaces. That is, when the period length is sufficiently long, phonons are difficult to produce interference patterns, and so the phonon transport is diffusive and corresponding thermal conductivity monotonically increases with the number of scattering interfaces. From the kinetic theory perspective, the phonon interference might generate a phonon-gap, and decline the system’s thermal conductivity to some extent. More importantly, the interference could flatten the phonon dispersion in a large frequency range, which affects the group velocity, as discussed in In-Plane Heterostructures. Note that the strict condition of phonon interference actually limits the applications of one-dimensional superlattice.
Lately, a metamaterial constructed by arrays of pillars on the nanofilms/nanowires has attracted more and more attention. The presence of nanopillars causes the formation of various resonant modes in the whole frequency range. If the resonant modes have the same polarization and frequency with the propagating modes from the base material, both kinds of phonons will hybridize due to the band anticrossing effect, and the phonon group velocity is lowered. The advantage of this design is that the long-wavelength phonons can be effectively controlled. For instance, using BTE method, (Davis and Hussein, 2014) calculated the thermal conductivity of nanopillared silicon thin films, and found that the obtained value is less than half of that in pristine nanowires. This is attributed to the hybridizations taking place at both subwavelength and superwavelength frequencies. Moreover, based on the NEMD simulations, (Xiong et al., 2016) reported that introducing resonant structures to alloying structures could achieve ultra-low thermal conductivity (0.9 Wm−1 K−1). It can be understood that the local resonance effect impedes the transport of low-frequency phonons and alloying can effectively scatter high-frequency phonons. Subsequent studies also confirmed this conclusion (Liu et al., 2018b; Nomura et al., 2018). Apart from silicon-based nanostrucutres, (Ma et al., 2018) reported a reduction of half thermal conductivity in nanopillared GNRs compared to pristine nanoribbons. It was also found that isotopic doping in nanopillars is conducive to phonon transport, because doping could weaken the hybridization degree. By designing the pristine/branched GNR junctions, (Chen et al., 2018b) proposed a high-performance thermal rectifier where its TR ratio is as high as 400% under a small temperature bias. The underlying mechanism was found to be that when the longitudinal phonons exists obvious local resonance for the negative heat flux. Experimentally, (Anufriev et al., 2017) observed a reduced thermal conductivity in silicon films by attaching aluminum pillars. However, they argued that the reduction of thermal conductivity originates from interface roughness induced by the metal deposition rather than phonon local resonance. In other words, the validity of phonon local resonance lowering thermal conductivity is still not clear. Furthermore, the design 2D hole-based phononic crystal is another feasible method to manipulate phonon transport. (Feng and Ruan, 2016) claimed that the thermal conductivity of graphene phononic crystal (GPC) is three orders of magnitude lower than the pristine one and decreases exponentially with increasing porosity. The vibrational eigen-mode analysis revealed that the phonon localization and phonon back scattering around the nanopores are the major factor responsible for the ultra-low thermal conductivity of GPC. Similar results have also been reported in some other studies (Hu et al., 2019). Analogously, (Cui et al., 2018a) reported that the thermal conductivity of silicene phononic crystal is not sensitive to temperature but strongly depends on the porosity ratio, as that in 3D phononic crystal (Hopkins et al., 2011). There were three reasons for that, namely increasing phonon scatterings, inducing phonon localization and phonon coherent effects.
Phonon-Driven Emerging Applications of 2D Heterostructres
Thermal Management
From the point of view of informatics, phonons as heat carriers can be also used for data storage and processing, namely called phononics. Various thermal devices (corresponding to the electronic counterparts) like rectifier, transistor, and logic gate have been conceptualized (Li et al., 2012; Jia et al., 2020), making the design of thermal logic circuits possible. Then, some of them are gradually realized via designing different sophisticated nanojunctions. To be specific, thermal rectifier as basic components that enables heat flow preferentially in one direction has been explored both theoretically and experimentally. For example, (Chang et al., 2006) experimentally observed the TR phenomenon in gradual mass-loading BN nanotubes opened the active regulation of heat flux in low-dimensional materials (Chang et al., 2006). After that, numerous architecture models like graded core-shell nanowires (Liu et al., 2014e), asymmetric nanobibbons (Zhang and Zhang, 2011; Sandonas et al., 2017; Li et al., 2018c) and phase change materials (Kang et al., 2018) have been proposed to use as nanoscale thermal rectifiers. For example, the TR ratio of the single-layer Y-shaped GNRs was predicted to be about 45% under a normalized temperature bias () of 0.5 and increases with the number of layer (Zhang and Zhang, 2011). With the aid of AGF method, (Ouyang et al., 2010) investigated the ballistic phonon transport in three-terminal GNRs and showed that the TR effect rises quickly with increasing the structural asymmetry. Moreover, the corner shape of nanojunctions was found to play an important role in the rectification effect, and the TR ratio could approach 200% via the parameters optimization. (Wang et al., 2014) confirmed that the TR effect in asymmetric GNRs gradually weakens with increasing the sample width. The mechanism origin was that phonon lateral confinement results in the mismatch of phonon spectra, phonon edge localization and the temperature and space dependence of thermal conductivity, all of which colludes to produce asymmetric phonon transport. (Tian et al., 2012a) experimentally designed a novel solid-state thermal rectifier based on the triangular-shaped graphene oxide. The measured TR ratio was up to 121% under a temperature bias of 12 K and increased with the increasing of the angle.
Apart from geometric asymmetry and mass asymmetry, the introduction of functionalized or defected engineering is also able to realize efficient TR effect. (Pal and Puri, 2014) proposed asymmetrically polymer-functionalized carbon nanotubes (CNTs) to obtain high TR ratio. When the applied temperature bias was 40 K, the TR ratio could be as high as 204% in this complex structure. From this work, it can be summarized that the structural asymmetry and the strong temperature dependence of vibration density of states are two necessary conductions to obtain high TR ratio. (Melis et al., 2015) reported that the TR ratio is about 54% in pristine/hydrogenated graphene heterojuction with triangular morphology. (Hu et al., 2017) proposed the thought of improving TR performance by tendering multi thermal rectifiers. (Nobakht et al., 2018) studied TR in defect-engineered graphene with asymmetric holing arrangements. It was found that the TR ratio increases as porosity increases, and can reach as high as 78% at room temperature with triangular hole arrangement with porosity of 75%, enabling effective and precise TR. Using microbridge method, (Wang et al., 2017a) experimentally demonstrated the measured TR ratio is about 26% in asymmetric defect graphene, which is ascribed to temperature and space dependence of thermal conductivity. It is worth pointing out that the aforementioned good TR effect commonly requires large temperature bias and delicate nanostructures namely complicated synthetic procedures, which limits the practical applications. Recently, phase change materials show promising potential in high-performance thermal rectifier. Utilizing the crystalline polymer nanofiber/cross-linked polymer nanofiber heterojunctions, (Zhang and Luo, 2015b) constructed a new thermal rectifier and its workflow was that the polymer nanofiber portion undergoes a phase transition and its thermal conductivity dramatically alters under temperature gradient. As a result, a 20 K temperature bias stimulated ∼100% TR ratio. Detailed spectral analysis revealed that the large TR ratio was ascribed to by the vibrational mismatch of sulfur and carbon atoms. (Kang et al., 2018) developed an analytical framework upon a material architecture for actualizing one type of nonlinear thermal transport, and applied the analytical framework to a junction comprising two phase-change materials (polyethylene and vanadium dioxide). The results showed that such heterojunction has a maximal theoretical rectification factor of approximately 140%. However, associated problems include slow phase change and narrow temperature range, thus hindering practical applications.
Actually, constructing nanojunctions is an effective method to implement efficient thermal rectifier. Base on NEMD simulations, (Yang et al., 2017) demonstrated that the TR ratio is as high as 1,200% in graphene/CNTs junction under 300 K temperature bias, approaching the actual application, which is considered as the formation of standing wave along negative heat flux direction. (Chen et al., 2019b) reported a high TR ratio of 334% in carbon/h-BN heteronanotubes can achieve 334%. Similar phenomenon has also been found in gaphene/h-BN in-plane heterostructures (Chen et al., 2016d). Note that a moderate C/BN ratio and parallel arrangement help improve the TR ratio in CBN hybrids (Chen et al., 2016b). (Dong et al., 2019) reported that the TR ratio in Au/CNTs heterojunctions decreases with the number of the molecular bridges but increases with the length of CNTs, reaching 375% under 60 K temperature bias. (Sandonas et al., 2017) indicated that the TR ratio in supported graphene/h-BN in-plane heterojunctions increases from 8.8 to 79% by lowering the substrate temperature. Recently, (Wei et al., 2019) proposed a new thermal rectifier base on multilayer graphene-based with interlayer gradient functionalization. The TR mechanism was that the vocational of each graphene layer exists obvious difference. (Chen et al., 2020b) reported a distinct TR effect in asymmetric graphene/h-BN vdW heterostructure. Moreover, the TR ratio increased rapidly with the interlayer coupling level, as shown in Figure 7A. Detailed spectral analyses revealed that for strong interlayer coupling, the low-frequency phonons presents the characteristic of localized modes along negative heat flux direction, while , medium-frequency phonons become high-efficient heat conduction channels, as illustrated in Figures 7B,C.
FIGURE 7
Apart from TR effect, another crucial thermal transport phenomenon is negative differential thermal resistance (NDTR) in which heat flux is inversely correlated with thermal bias. Actually, NDTR is the foundation not only for thermal switching, but also for thermal memory. By now, the phenomenon has been widely investigated in non-linear lattice models (Chan et al., 2014). Some obtained results showed that the existence of interface effect is a requirement for generating NDTR. Subsequently, (Ai et al., 2012) found that the phenomenon of NDTR only exists in single-layer GNRs with appropriate width. For multi-layer GNRs, NDTR regime becomes gradually smaller and even disappears with increasing the number of layers. By designing sandwiched nanofluids, (Li et al., 2019) demonstrated that the heat flux is suppressed and even prohibited when temperature bias becomes sufficiently large, because nanofluids are adsorbed on the solid surface under negative temperature bias, and the free fluid molecules is decreased. Furthermore, (Chen et al., 2017b) studied the nonlinear thermal transport in graphene/h-BN in-plane heterostructure and showed that when surpasses a threshold, NDTR phenomenon appears. It could be understood that out-of-plane kinetic wave is stimulated in graphene domain under large temperature bias, which results in the out-of-plane vibration mismatch between graphene and h-BN regions, impeding interfacial thermal transfer, as illustrated in Figure 8. It was also found that the regime of NDTR becomes small as the sample length increases, but barely changes with the sample width. That is, the wave-dominated thermal transport mechanism becomes dominant when the sample length is small or the applied temperature bias is large. To realize logic calculations uisng phonons, (Pal and Puri, 2015) designed a thermal AND gate consisting of asymmetric GRN and the test of a truth table demonstrated that with either inputs at “0” state, the output value is in the range 0∼0.27, which is comparable to silicon-based logic gate (0∼0.12). Note that the running speed of the thermal device reaches 100 ps owing to ultrahigh phonon group velocity in graphene. (Liu et al., 2015) experimentally demonstrated a thermal modulator consisting of a partially clamped single-layer graphene, and found the high/low heat flux ratio can remarkably is as high as 150% by regulating the external pressure. The reason for the decrease of heat flux was that the phonon spectra mismatch of the unclamped and clamped graphene sections leads to strong phonon-interface scattering. Another experiment, (Sood et al., 2018) prepared a thermal transistor based on multi-layer MoS2 driven by reversible electro-chemical intercalation of Li ions. The thermal device is in the following: during the charging the cell, Li ions continue to enter the multilayer MoS2 and some 2H phase MoS2 layers are converted into the 1T phase, namely 2H-1T phase interface, which can lowers phonon group velocity and enhances phonon scattering rates. Moreover, lithiation may induce disordered stacking among multilayer MoS2 and hinders the interlayer thermal transport, causing the decrease of ITC from 15 to 1.6 MW m−2 K−1. For the delithiation step, the cross-plane thermal conductance increases back to the pre-lithiation value. The obtained on/off ratio is about 10, opening the possibility of exactly controlling heat flux, as shown in Figure 9.
FIGURE 8
FIGURE 9
Thermoelectric Conversion
Thermoelectric materials enable heat energy to directly convert into electric power. The efficiency of thermoelectric conversion is determined by a dimensionless parameter, namely figure-of-merit (), where is Seebeck coefficient, is electrical conductance, is absolute temperature, and is electron and phonon thermal conductance, respectively. It is not difficult to find that enhancing power factor () and reducing thermal conductivity are both effective strategies to promote thermoelectric performance. However, according to the Wiedeman-Franz law, , and depend on each other. In contrast, inhibiting phonon transport is more viable realizing high conversion efficiency. In the last few decades, a huge amount of effort has been conducted to reduce the thermal conductivity of thermoelectric materials via nanostructure engineering (Zhou et al., 2020b).
Reducing dimensionality as an important strategy for enhancing thermoelectric performance has been proposed by (Heremans et al., 2013), and its working mechanisms are as follows: 1) the size-quantization in low-dimensional nanostrucutres facilitates increasing the Seebeck coefficient; 2) the phonon-boundary scatterings can reduces the system’s thermal conductivity. In fact, it has reported that the electronic transmission exists obvious fluctuation within the first conductance band for asymmetric GNRs, and consequently the Seebeck coefficient is enhanced (Xie et al., 2012; Cao et al., 2019). (Yokomizo and Nakamura, 2013) demonstrated that the value of in graphene/h-BN superlattices is twenty times larger than that in graphene. (Mazzamuto et al., 2011) reported that the electron transport exits resonant tunneling behavior in mixed armchair/zigzag GNRs, which corresponds to a ZT value of 1 at room temperature. On the other hand, (Tran et al., 2015) found the phonon conductance is significantly inhibited in GNRs with BNNRs branches, but the electron transmission is equivalent to that in pristine GNRs, similar to that in vdW graphene junction (Hung Nguyen et al., 2014). Moreover, (Cui et al., 2018b) reported that the ZT value of graphyne nanoribbons with edge disorder can achieve 2.5 because of the strong phonon-boundary scatterings. Along the way, the MoS2/WSe2 (Jia et al., 2019) and graphene/MoS2 vdW heterostructures (Sadeghi et al., 2016) have been proposed to construct thermoelectronic devices with excellent performance. Lately, the electronic devices based on organic single-molecules have been shown excellent electronic properties such as quantum tunneling and magnetic change (Chen et al., 2017a; Deng et al., 2017; Xie et al., 2017; Zhang et al., 2019) and some of them can be also used as the competitive thermoelectric materials (Wu et al., 2019; Wu et al., 2020; Zeng et al., 2020b). Actually, sometimes it should balance the power output and thermoelectric conversion efficiency. That is, blindly reducing the system’s thermal conductivity is not always better (Narducci, 2011). As shown in Figure 10C, (Barrios-Vargas et al., 2017) reported an 80% reduction for the thermal conductivity of polycrystalline graphene/h-BN in-plane heterostructures. Meanwhile, the electrical conductivity also declined rapidly. As a result, the calculated ZT value was only 0.01. Experimentally, (Mahapatra et al., 2017) measured the electron/phonon transport in the twisted bilayer graphene, and demonstrated that the cross-plane thermoelectric transport is driven by the scattering of electrons and breathing phonon modes that represent a unique “phonon drag” effect across atomic distances, as illustrated in Figure 10A. Moreover, it was also found that the value of in gold/graphene/WSe2 vdW heterostructrues approaches 100 (Rosul et al., 2019).
FIGURE 10
FIGURE 11
To sum up, 2D heterostructures have great potential for the applications of thermoelectric conversion. However, theoretical calculations still have large limitations due to the limitation of computation cost and the handing question of van der Waals force in the DFT calculation. In particular, the electron–phonon coupling was not considered in the calculation process. In fact, some previous studies have indicated that the electron–phonon coupling can severely affect the electron transport properties such as spin filter and magnetoresistances (Liu et al., 2017b; Li and Chen, 2017). Again, (Zhu et al., 2019) found that the presence of bipolar effect in SnP3 could lower the system thermal conductivity and help enhance the ZT value (∼3.4). Hence, new theoretical method including multiparticle coupling is necessary to predict the thermoelectric performance of nanostructurs.
Conclusion and Outlook
In this review, we have introduced various experimental and theoretical methods for the thermal transport in low-dimensional materials. Then, the thermal transport properties of 2D materials and their in-plane/vdW heterostructures were discussed. Furthermore, we have summarized different thoughts to manipulate thermal conductivity. In practical terms, the introduction of isotopes, dislocation, substrates, chemisorption, and strain can strength diffuse phonon scattering (particle-like), which is valid with a limited range of high-frequency phonons. Meanwhile, the wave nature of phonons like interference or local resonance can regulate low-frequency phonons. The diversity of thermal properties in 2D heterostructures indicates important potential applications in thermoelectric, thermal dissipation and thermal devices fields, which needs further investigations.
Advanced fabrication techniques have made multistage circuit integration of 2D heterostructures possible, and resolving the thermal dissipation across the interfaces also looms ahead when designing nanodevices. Considering that heat transfer in semiconductor device is carried by a large population of phonons spanning a frequency spectrum and existing diverse interactions, controlling phonon transport is actually a complex task in these nanostructures and some important questions are still not clear. First, the phonon transport in nanostructures under coupling multi-field effect like electron-phonon or magneton-phonon couplings is in the primary stage of development. Some already obtained results are based on either theoretical models or based on multiple hypotheses, and haven’t been verified by experiments. Second, because of the multi-scale character of phonon transport, the accurate control of phonon transport in specific frequency range remains surprisingly difficult. Moreover, the anomalous non-Fourier transport caused by the excitation of the kinetic wave needs further confirmation in experiment. Third, on the experiment side, the accurate measurement of thermal properties of 2D heterostructures still has a long way to go.
Funding
This work was supported by the National Key Research and Development Program of China (No. 2017YFB0701602), by the National Natural Science Foundation of China (Nos. 11674092, 11904161), by Hunan Provincial Natural Science Foundation of China (No. 2018JJ3421), by the Scientific Research Fund of Hunan Provincial Education Department (No. 18C0454) and by the Open Fund of Hunan Provincial Key Laboratory of Advanced Materials for New Energy Storage and Conversion (No. 2018TP1037-202004).
Statements
Author contributions
X-KC, Y-JZ and K-QC conceived the manuscript. All authors reviewed and discussed the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AbidiI. H.WengL.-T.WongC. P. J.TyagiA.GanL.DingY.et al (2018). New approach to unveiling individual atomic layers of 2D materials and their heterostructures. Chem. Mater.30, 1718–1728. 10.1021/acs.chemmater.7b05371
2
AiB. Q.ZhongW. R.HuB. (2012). Dimension dependence of negative differential thermal resistance in graphene nanoribbons. J. Phys. Chem. C116, 13810–13815. 10.1021/jp303431k
3
AiyitiA.HuS.WangC.XiQ.ChengZ.XiaM.et al (2018). Thermal conductivity of suspended few-layer MoS2. Nanoscale10, 2727–2734. 10.1039/c7nr07522g
4
Al TalebA.FaríasD. (2016). Phonon dynamics of graphene on metals. J. Phys. Condens. Matter28, 103005. 10.1088/0953-8984/28/10/103005
5
AnufrievR.YanagisawaR.NomuraM. (2017). Aluminium nanopillars reduce thermal conductivity of silicon nanobeams. Nanoscale9, 15083–15088. 10.1039/c7nr05114j
6
BalandinA. A.GhoshS.BaoW.CalizoI.TeweldebrhanD.MiaoF.et al (2008). Superior thermal conductivity of single-layer graphene. Nano Lett.8, 902–907. 10.1021/nl0731872
7
Barrios-VargasJ. E.MortazaviB.CummingsA. W.Martinez-GordilloR.PrunedaM.ColomboL.et al (2017). Electrical and thermal transport in coplanar polycrystalline graphene–hBN heterostructures. Nano Lett.17, 1660–1664. 10.1021/acs.nanolett.6b04936
8
BehranginiaA.HemmatZ.MajeeA. K.FossC. J.YasaeiP.AksamijaZ.et al (2018). Power dissipation of WSe2 field-effect transistors probed by low-frequency Raman thermometry. ACS Appl. Mater. Interfaces10, 24892–24898. 10.1021/acsami.8b04724
9
BerciaudS.RyuS.BrusL. E.HeinzT. F. (2009). Probing the intrinsic properties of exfoliated graphene: Raman spectroscopy of free-standing monolayers. Nano Lett.9, 346–352. 10.1021/nl8031444
10
BertolazziS.KrasnozhonD.KisA. (2013). Nonvolatile memory cells based on MoS2/graphene heterostructures. ACS Nano7, 3246–3252. 10.1021/nn3059136
11
CahillD. G.BraunP. V.ChenG.ClarkeD. R.FanS.GoodsonK. E.et al (2014). Nanoscale thermal transport. II. 2003–2012. App. Phys. Rev.1, 011305 . 10.2172/1160069
12
CaiY.LanJ.ZhangG.ZhangY. W. (2014). Lattice vibrational modes and phonon thermal conductivity of monolayer MoS2. Phy. Rev. B.89 (3), 035438. 10.1103/physrevb.89.035438
13
CaiY.PeiQ. X.ZhangG.ZhangY. W. (2016). Decoupled electron and phonon transports in hexagonal boron nitride-silicene bilayer heterostructure. J. Appl. Phys.119 (6), 065102. 10.1063/1.4941534
14
CaoX. H.WuD.FengY.-X.ZhouW.-X.TangL.-M.ChenK.-Q. (2019). Effect of electrophilic substitution and destructive quantum interference on the thermoelectric performance in molecular devices. J. Phys. Condens. Matter31, 345303. 10.1088/1361-648x/ab2299
15
ChanH. K.HeD.HuB. (2014). Scaling analysis of negative differential thermal resistance. Phys. Rev. E89 (5), 052126. 10.1103/physreve.89.052126
16
ChangC. W.OkawaD.MajumdarA.ZettlA. (2006). Solid-state thermal rectifier. Science314, 1121–1124. 10.1126/science.1132898
17
ChangH. Y.YangS.LeeJ.TaoL.HwangW.-S.JenaD.et al (2013). High-performance, highly bendable MoS2 transistors with high-k dielectrics for flexible low-power systems. ACS Nano, 7, 5446–5452. 10.1021/nn401429w
18
ChenX.-K.ChenK.-Q. (2020a). Thermal transport of carbon nanomaterials. J. Phys. Condens. Matter32, 153002. 10.1088/1361-648x/ab5e57
19
ChenK. Q.LiW. X.DuanW.ShuaiZ.GuB. L. (2005). Effect of defects on the thermal conductivity in a nanowire. Phys. Rev. B72 (4), 045422. 10.1103/physrevb.72.045422
20
ChenS.WuQ.MishraC.KangJ.ZhangH.ChoK.et al (2012). Thermal conductivity of isotopically modified graphene. Nat. Mater.11, 203–207. 10.1038/nmat3207
21
ChenC. C.LiZ.ShiL.CroninS. B. (2014). Thermal interface conductance across a graphene/hexagonal boron nitride heterojunction. App. Phy. Lett.104 (8), 081908. 10.1063/1.4866335
22
ChenX.WuY.WuZ.HanY.XuS.WangL.et al (2015b). High-quality sandwiched black phosphorus heterostructure and its quantum oscillations. Nat. Commun.6, 7315. 10.1038/ncomms8315
23
ChenJ.ChenS.GaoY. (2016a). Anisotropy enhancement of thermal energy transport in supported black phosphorene. J. Phys. Chem. Lett.7, 2518–2523. 10.1021/acs.jpclett.6b00858
24
ChenX.-K.XieZ.-X.ZhouW.-X.ChenK.-Q. (2016b). The thermal conductivity in hybridised graphene and boron nitride nanoribbons modulated with strain. J. Phys. D Appl. Phys.49, 115301. 10.1088/0022-3727/49/11/115301
25
ChenX. K.XieZ. X.ZhouW. X.TangL. M.ChenK. Q. (2016c). Phonon wave interference in graphene and boron nitride superlattice. Appl. Phys. Lett.109 (2), 023101. 10.1063/1.4958688
26
ChenX.-K.XieZ.-X.ZhouW.-X.TangL.-M.ChenK.-Q. (2016d). Thermal rectification and negative differential thermal resistance behaviors in graphene/hexagonal boron nitride heterojunction. Carbon100, 492–500. 10.1016/j.carbon.2016.01.045
27
ChenT.LiH.ZhangY.LiuD.ChaoY.WangL. (2017a). Parity effects induced by the resonant electronic states coupling in polyacetylene-based devices. J. Electron. Mater.46, 5121–5126. 10.1007/s11664-017-5501-1
28
ChenX. K.LiuJ.PengZ. H.DuD.ChenK. Q. (2017b). A wave-dominated heat transport mechanism for negative differential thermal resistance in graphene/hexagonal boron nitride heterostructures. Appl. Phys. Lett.110 (9), 091907. 10.1063/1.4977776
29
ChenY.ZhangY.CaiK.JiangJ.ZhengJ.-C.ZhaoJ.et al (2017c). Interfacial thermal conductance in graphene/black phosphorus heterogeneous structures. Carbon117, 399–410. 10.1016/j.carbon.2017.03.011
30
ChenT.GuoC.XuL.LiQ.LuoK.LiuD.et al (2018a). Modulating the properties of multi-functional molecular devices consisting of zigzag gallium nitride nanoribbons by different magnetic orderings: a first-principles study. Phys. Chem. Chem. Phys.20, 5726. 10.1039/c7cp07467k
31
ChenX.-K.LiuJ.XieZ.-X.ZhangY.DengY.-X.ChenK.-Q. (2018b). A local resonance mechanism for thermal rectification in pristine/branched graphene nanoribbon junctions. Appl. Phys. Lett.113, 121906. 10.1063/1.5053233
32
ChenX.ChenC.LeviA.HoubenL.DengB.YuanS.et al (2018c). Large-velocity saturation in thin-film black phosphorus transistors. ACS Nano, 12, 5003–5010. 10.1021/acsnano.8b02295
33
ChenT.XuL.LiQ.LiX.LongM. (2019a). Direction and strain controlled anisotropic transport behaviors of 2D GeSe-phosphorene vdW heterojunctions. Nanotechnology30, 445703. 10.1088/1361-6528/ab375b
34
ChenX.-K.XieZ.-X.ZhangY.DengY.-X.ZouT.-H.LiuJ.et al (2019b). Highly efficient thermal rectification in carbon/boron nitride heteronanotubes. Carbon148, 532–539. 10.1016/j.carbon.2019.03.073
35
ChenX.-K.PangM.ChenT.DuD.ChenK.-Q. (2020b). Thermal rectification in asymmetric graphene/hexagonal boron nitride van der Waals heterostructures. ACS Appl. Mater. Interfaces12, 15517–15526. 10.1021/acsami.9b22498
36
CherukaraM. J.NarayananB.KinaciA.SasikumarK.GrayS. K.ChanM. K. Y.et al (2016). Ab initio-based bond order potential to investigate low thermal conductivity of stanene nanostructures. J. Phys. Chem. Lett.7, 3752–3759. 10.1021/acs.jpclett.6b01562
37
CiL.SongL.JinC.JariwalaD.WuD.LiY.et al (2010). Atomic layers of hybridized boron nitride and graphene domains. Nat. Mater.9, 430–435. 10.1038/nmat2711
38
CuiL.ShiS.LiZ.WeiG.DuX. (2018a). Reduction of thermal conductivity in silicene nanomesh: insights from coherent and incoherent phonon transport. Phys. Chem. Chem. Phys.20, 27169–27175. 10.1039/c8cp03993c
39
CuiX.OuyangT.LiJ.HeC.TangC.ZhongJ. (2018b). Enhancing the thermoelectric performance of gamma-graphyne nanoribbons by introducing edge disorder. Phys. Chem. Chem. Phys.20, 7173–7179. 10.1039/c7cp08154e
40
da SilvaC.SaizF.RomeroD. A.AmonC. H. (2016). Coherent phonon transport in short-period two-dimensional superlattices of graphene and boron nitride. Phy. Rev. B93, 125427. 10.1103/physrevb.93.125427
41
DavisB. L.HusseinM. I. (2014). Thermal conductivity reduction by local resonance. Phys. Rev. Lett.112 (5), 055505. 10.1103/physrevlett.112.055505
42
DengY.-X.ChenS.-Z.ZengY.ZhouW.-X.ChenK.-Q. (2017). Large spin rectifying and high-efficiency spin-filtering in superior molecular junction. Org. Electron.50, 184–190. 10.1016/j.orgel.2017.07.046
43
DengY.-X.ChenS.-Z.ZengY.FengY.ZhouW.-X.TangL.-M.et al (2018). Spin gapless semiconductor and half-metal properties in magnetic penta-hexa-graphene nanotubes. Org. Electron.63, 310–317. 10.1016/j.orgel.2018.09.046
44
DingZ.PeiQ.-X.JiangJ.-W.ZhangY.-W. (2015). Manipulating the thermal conductivity of monolayer MoS2 via lattice defect and strain engineering. J. Phys. Chem. C119, 16358–16365. 10.1021/acs.jpcc.5b03607
45
DingZ.PeiQ.-X.JiangJ.-W.HuangW.ZhangY.-W. (2016). Interfacial thermal conductance in graphene/MoS2 heterostructures. Carbon96, 888–896. 10.1016/j.carbon.2015.10.046
46
DingG.HeJ.GaoG. Y.YaoK. (2018a). Two-dimensional MoS2-MoSe2 lateral superlattice with minimized lattice thermal conductivity. J. Appl. Phys.124, 165101. 10.1063/1.5051067
47
DingY.ZhouN.GanL.YanX.WuR.AbidiI. H.et al (2018b). Stacking-mode confined growth of 2H-MoTe2/MoS2 bilayer heterostructures for UV-vis-IR photodetectors. Nanomater. Energy49, 200–208. 10.1016/j.nanoen.2018.04.055
48
DollfusP.NguyenV. H.Saint-MartinJ. (2015). Thermoelectric effects in graphene nanostructures. J. Phys. Condens. Matter27, 133204. 10.1088/0953-8984/27/13/133204
49
DongY.DiaoC.SongY.ChiH.SinghD. J.LinJ. (2019). Molecular bridge thermal diode enabled by vibrational mismatch. Phys. Rev. Appl.11, 024043. 10.1115/mnhmt2019-3905
50
FanZ. Q.ChenK. Q. (2010). Negative differential resistance and rectifying behaviors in phenalenyl molecular device with different contact geometries. Appl. Phys. Lett.96 (5), 053509. 10.1063/1.3309708
51
FanZ.PereiraL. F. C.HirvonenP.ErvastiM. M.ElderK. R.DonadioD.et al (2017). Thermal conductivity decomposition in two-dimensional materials: application to graphene. Phys. Rev. B95, 144309. 10.1103/physrevb.95.144309
52
FanZ.-Q.ZhangZ.-H.YangS.-Y. (2020). High-performance 5.1 nm in-plane Janus WSeTe Schottky barrier field effect transistors. Nanoscale12, 21750. 10.1039/d0nr05269h
53
FeiR.FaghaniniaA.SoklaskiR.YanJ.-A.LoC.YangL. (2014). Enhanced thermoelectric efficiency via orthogonal electrical and thermal conductances in phosphorene. Nano Lett.14, 6393–6399. 10.1021/nl502865s
54
FengT.RuanX. (2016). Ultra-low thermal conductivity in graphene nanomesh. Carbon101, 107–113. 10.1016/j.carbon.2016.01.082
55
FengT.YaoW.WangZ.ShiJ.LiC.CaoB.et al (2017). Spectral analysis of nonequilibrium molecular dynamics: spectral phonon temperature and local nonequilibrium in thin films and across interfaces. Phys. Rev. B95, 195202. 10.1103/physrevb.95.195202
56
FrisendaR.Navarro-MoratallaE.GantP.Pérez De LaraD.Jarillo-HerreroP.GorbachevR. V.et al (2018). Recent progress in the assembly of nanodevices and van der Waals heterostructures by deterministic placement of 2D materials. Chem. Soc. Rev.47, 53–68. 10.1039/c7cs00556c
57
FuL.SunY.WuN.MendesR. G.ChenL.XuZ.et al (2016). Direct growth of MoS2/h-BN heterostructures via a sulfide-resistant alloy. ACS Nano10, 2063–2070. 10.1021/acsnano.5b06254
58
FuY.HanssonJ.LiuY.ChenS.ZehriA.SamaniM. K.et al (2019). Graphene related materials for thermal management. 2D Mater7 (1), 012001. 10.1088/2053-1583/ab48d9
59
FugalloG.CepellottiA.PaulattoL.LazzeriM.MarzariN.MauriF. (2014). Thermal conductivity of graphene and graphite: collective excitations and mean free paths. Nano Lett.14, 6109–6114. 10.1021/nl502059f
60
GaoY.XuB. (2017). Controllable interface junction, in-plane heterostructures capable of mechanically mediating on-demand asymmetry of thermal transports. ACS Appl. Mater. Interfaces9, 34506–34517. 10.1021/acsami.7b11508
61
GaoY.ZhangY.ChenP.LiY.LiuM.GaoT.et al (2013). Toward single-layer uniform hexagonal boron nitride-graphene patchworks with zigzag linking edges. Nano Lett.13, 3439–3443. 10.1021/nl4021123
62
GaoY.LiuQ.XuB. (2016). Lattice mismatch dominant yet mechanically tunable thermal conductivity in bilayer heterostructures. ACS Nano10, 5431–5439. 10.1021/acsnano.6b01674
63
GaoM.ZhangW.ZhangL. (2018). Nondegenerate chiral phonons in graphene/hexagonal boron nitride heterostructure from first-principles calculations. Nano Lett.18, 4424–4430. 10.1021/acs.nanolett.8b01487
64
MaR.PerryC. H.RupprechtG. (1966). Normal modes in hexagonal boron nitride. Phys. Rev.146, 543. 10.1103/physrev.146.543
65
GeorgiouT.JalilR.BelleB. D.BritnellL.GorbachevR. V.MorozovS. V.et al (2013). Vertical field-effect transistor based on graphene–WS2 heterostructures for flexible and transparent electronics. Nat. Nanotechnol.8, 100–103. 10.1038/nnano.2012.224
66
GholivandH.FuladiS.HemmatZ.Salehi-KhojinA.Khalili-AraghiF. (2019). Effect of surface termination on the lattice thermal conductivity of monolayer Ti3C2Tz MXenes. J. Appl. Phys.126, 065101. 10.1063/1.5094294
67
GuX.YangR. (2014). Phonon transport in single-layer transition metal dichalcogenides: a first-principles study. Appl. Phys. Lett.105, 131903. 10.1063/1.4896685
68
GuX.YangR. (2015). First-principles prediction of phononic thermal conductivity of silicene: a comparison with graphene. J. Appl. Phys.117 (2), 025102. 10.1063/1.4905540
69
GuX.YangR. (2016). Phonon transport in single-layer MO1−x Wx S2 alloy embedded with WS2 nanodomains. Phys. Rev. B94, 075308. 10.1103/physrevb.94.075308
70
GuX.WeiY.YinX.LiB.YangR. (2018). Colloquium: phononic thermal properties of two-dimensional materials. Rev. Mod. Rhys90, 041002. 10.1103/revmodphys.90.041002
71
GuoZ.MiaoN.ZhouJ.PanY.SunZ. (2018). Coincident modulation of lattice and electron thermal transport performance in MXenes via surface functionalization. Phys. Chem. Chem. Phys.20, 19689–19697. 10.1039/c8cp02564a
72
GuoM.QianY.QiH.BiK.ChenY. (2020). Experimental measurements on the thermal conductivity of strained monolayer graphene. Carbon157, 185–190. 10.1016/j.carbon.2019.10.027
73
HeJ.LiD.YingY.FengC.HeJ.ZhongC.et al (2019). Orbitally driven giant thermal conductance associated with abnormal strain dependence in hydrogenated graphene-like borophene. Npj Comput. Mater5, 47. 10.1038/s41524-019-0183-2
74
HemmatZ.YasaeiP.SchultzJ. F.HongL.MajidiL.BehranginiaA.et al (2019). Tuning thermal transport through atomically thin Ti3C2TzMXene by current annealing in vacuum. Adv. Funct. Mater.29, 1805693. 10.1002/adfm.201805693
75
HeremansJ. P.DresselhausM. S.BellL. E.MorelliD. T. (2013). When thermoelectrics reached the nanoscale. Nat. Nanotechnol.8, 471. 10.1038/nnano.2013.129
76
HongY.ZhangJ.ZengX. C. (2016a). Interlayer thermal conductance within a phosphorene and graphene bilayer. Nanoscale8, 19211–19218. 10.1039/c6nr07977f
77
HongY.ZhangJ.ZengX. C. (2016b). Thermal conductivity of monolayer MoSe2 and MoS2. J. Phys. Chem. C120, 26067–26075. 10.1021/acs.jpcc.6b07262
78
HongY.ZhangJ.ZengX. C. (2016c). Thermal contact resistance across a linear heterojunction within a hybrid graphene/hexagonal boron nitride sheet. Phys. Chem. Chem. Phys.18, 24164–24170. 10.1039/c6cp03933b
79
HongY.JuM. G.ZhangJ.ZengX. C. (2018). Phonon thermal transport in a graphene/MoSe2 van der Waals heterobilayer. Phys. Chem. Chem. Phys.20, 2637–2645. 10.1039/c7cp06874c
80
HopkinsP. E.ReinkeC. M.SuM. F.OlssonR. H.IIIShanerE. A.LesemanZ. C.et al (2011). Reduction in the thermal conductivity of single crystalline silicon by phononic crystal patterning. Nano Lett.11, 107–112. 10.1021/nl102918q
81
HuM.ZhangX.PoulikakosD. (2013). Anomalous thermal response of silicene to uniaxial stretching. Phy. Rev. B87, 195417. 10.1103/physrevb.87.195417
82
HuS.AnM.YangN.LiB. (2017). A series circuit of thermal rectifiers: an effective way to enhance rectification ratio. Small13, 1602726. 10.1002/smll.201602726
83
HuS.ZhangZ.JiangP.RenW.YuC.ShiomiJ.et al (2019). Disorder limits the coherent phonon transport in two-dimensional phononic crystal structures. Nanoscale11, 11839–11846. 10.1039/c9nr02548k
84
HuangH.XuY.ZouX.WuJ.DuanW. (2013). Tuning thermal conduction via extended defects in graphene. Phys. Rev. B87, 205415. 10.1103/physrevb.87.205415
85
HuangC.WuS.SanchezA. M.PetersJ. J. P.BeanlandR.RossJ. S.et al (2014). Lateral heterojunctions within monolayer MoSe2-WSe2 semiconductors. Nat. Mater.13, 1096–1101. 10.1038/nmat4064
86
Hung NguyenV.NguyenM. C.NguyenH.-V.Saint-MartinJ.DollfusP. (2014). Enhanced thermoelectric figure of merit in vertical graphene junctions. Appl. Phys. Lett.105, 133105. 10.1063/1.4896915
87
IannacconeG.BonaccorsoF.ColomboL.FioriG. (2018). Quantum engineering of transistors based on 2D materials heterostructures. Nat. Nanotechnol.13, 183–191. 10.1038/s41565-018-0082-6
88
JangW.BaoW.JingL.LauC. N.DamesC. (2013). Thermal conductivity of suspended few-layer graphene by a modified T-bridge method. Appl. Phys. Lett.103, 133102. 10.1063/1.4821941
89
JangH.WoodJ. D.RyderC. R.HersamM. C.CahillD. G. (2015). Anisotropic thermal conductivity of exfoliated black phosphorus. Adv. Mater.27, 8017–8022. 10.1002/adma.201503466
90
JiaP. Z.ZengY. J.WuD.PanH.CaoX. H.ZhouW. X.et al (2019). Excellent thermoelectric performance induced by interface effect in MoS2/MoSe2 van der Waals heterostructure. J. Phys.: condens. Matter32 (5), 055302. 10.1088/1361-648x/ab4cab
91
JiaP.-Z.WuD.ZhangQ.-Q.ZhouW.-X.FanZ.-Q.FengY.-X.et al (2020). Design of thermal metamaterials with excellent thermal control functions by using functional nanoporous graphene. Phys. Status Solidi RRL14, 2000333. 10.1002/pssr.202000333
92
JiangJ.-W.WangB.-S.WangJ.-S. (2011a). First principle study of the thermal conductance in graphene nanoribbon with vacancy and substitutional silicon defects. Appl. Phys. Lett.98, 113114. 10.1063/1.3567768
93
JiangJ. W.WangJ. S.WangB. S. (2011b). Minimum thermal conductance in graphene and boron nitride superlattice. Appl. Phys. Lett.99(4), 043109. 10.1063/1.3619832
94
JiangJ.-W.ParkH. S.RabczukT. (2013). Molecular dynamics simulations of single-layer molybdenum disulphide (MoS2): Stillinger-Weber parametrization, mechanical properties, and thermal conductivity. J. Appl. Phys.114, 064307. 10.1063/1.4818414
95
JiangJ. W. (2014). Phonon bandgap engineering of strained monolayer MoS2. Nanoscale6, 8326–8333. 10.1039/c4nr00279b
96
JoI.PettesM. T.KimJ.WatanabeK.TaniguchiT.YaoZ.et al (2013). Thermal conductivity and phonon transport in suspended few-layer hexagonal boron nitride. Nano Lett.13, 550–554. 10.1021/nl304060g
97
JoI.PettesM. T.OuE.WuW.ShiL. (2014). Basal-plane thermal conductivity of few-layer molybdenum disulfide. Appl. Phys. Lett.104, 201902. 10.1063/1.4876965
98
JungJ.QiaoZ.NiuQ.MacDonaldA. H. (2012). Transport properties of graphene nanoroads in boron nitride sheets. Nano Lett.12, 2936–2940. 10.1021/nl300610w
99
KakodkarR. R.FeserJ. P. (2017). Probing the validity of the diffuse mismatch model for phonons using atomistic simulations. Phys. Rev. B95, 125434. 10.1103/physrevb.95.125434
100
KangJ.JariwalaD.RyderC. R.WellsS. A.ChoiY.HwangE.et al (2016). Probing out-of-plane charge transport in black phosphorus with graphene-contacted vertical field-effect transistors. Nano Lett., 16, 2580–2585. 10.1021/acs.nanolett.6b00144
101
KangH.YangF.UrbanJ. J. (2018). Thermal rectification via heterojunctions of solid-state phase-change materials. Phys. Rev. Appl.10 (2), 024034. 10.1103/physrevapplied.10.024034
102
KeysharK.BergM.ZhangX.VajtaiR.GuptaG.ChanC. K.et al (2017). Experimental determination of the ionization energies of MoSe2, WS2, and MoS2 on SiO2 using photoemission electron microscopy. ACS Nano11, 8223–8230. 10.1021/acsnano.7b03242
103
KimJ. Y.LeeJ.-H.GrossmanJ. C. (2012). Thermal transport in functionalized graphene. ACS Nano6, 9050–9057. 10.1021/nn3031595
104
KimD.KimH.YunW. S.WatanabeK.TaniguchiT.RhoH.et al (2018). Energy dissipation mechanism revealed by spatially resolved Raman thermometry of graphene/hexagonal boron nitride heterostructure devices. 2D Mater5 (2), 025009. 10.1088/2053-1583/aaab14
105
KuangY.LindsayL.ShiS.WangX.HuangB. (2016). Thermal conductivity of graphene mediated by strain and size. Heat Mass Tran.101, 772–778. 10.1016/j.ijheatmasstransfer.2016.05.072
106
LeeS.LindsayL. (2017). Hydrodynamic phonon drift and second sound in a (20, 20) single-wall carbon nanotube. Phys. Rev. B95, 184304. 10.1103/physrevb.95.184304
107
LeeS.BroidoD.EsfarjaniK.ChenG. (2015a). Hydrodynamic phonon transport in suspended graphene. Nat. Commun.6, 6290. 10.1038/ncomms7290
108
LeeS.YangF.SuhJ.YangS.LeeY.LiG.et al (2015b). Anisotropic in-plane thermal conductivity of black phosphorus nanoribbons at temperatures higher than 100 K. Nat. Commun.6, 9573. 10.1038/ncomms9573
109
LiB.-L.ChenK.-Q. (2017). Effects of electron-phonon interactions on the spin-dependent Seebeck effect in graphene nanoribbons. Carbon119, 548–554. 10.1016/j.carbon.2017.04.069
110
LiD.WuY.KimP.ShiL.YangP.MajumdarA. (2003). Thermal conductivity of individual silicon nanowires. Appl. Phys. Lett.83, 2934–2936. 10.1063/1.1616981
111
LiN.RenJ.WangL.ZhangG.HänggiP.LiB. (2012). Colloquium: phononics: manipulating heat flow with electronic analogs and beyond. Rev. Mod. Phys.84, 1045. 10.1103/revmodphys.84.1045
112
LiH.YingH.ChenX.NikaD. L.CocemasovA. I.CaiW.et al (2014a). Thermal conductivity of twisted bilayer graphene. Nanoscale6, 13402–13408. 10.1039/c4nr04455j
113
LiW.CarreteJ.KatchoN. A.MingoN. (2014b). ShengBTE: a solver of the Boltzmann transport equation for phonons. Comput. Phys. Commun.185, 1747–1758. 10.1016/j.cpc.2014.02.015
114
LiM.-Y.ShiY.ChengC.-C.LuL.-S.LinY.-C.TangH.-L.et al (2015). Epitaxial growth of a monolayer WSe2-MoS2 lateral p-n junction with an atomically sharp interface. Science349, 524–528. 10.1126/science.aab4097
115
LiD.HeJ.DingG.TangQ.YingY.HeJ.et al (2018a). Stretch-driven increase in ultrahigh thermal conductance of hydrogenated borophene and dimensionality crossover in phonon transmission. Adv. Funct. Mater.28, 1801685. 10.1002/adfm.201801685
116
LiM.ZhengB.DuanK.ZhangY.HuangZ.ZhouH. (2018b). Effect of defects on the thermal transport across the graphene/hexagonal boron nitride interface. J. Phys. Chem. C122, 14945–14953. 10.1021/acs.jpcc.8b02750
117
LiT.TangZ.HuangZ.YuJ. (2018c). Substrate effects on the thermal performance of in-plane graphene/hexagonal boron nitride heterostructures. Carbon130, 396–400. 10.1016/j.carbon.2018.01.017
118
LiW.KongL.ChenC.GouJ.ShengS.ZhangW.et al (2018d). Experimental realization of honeycomb borophene. Sci. Bull.63, 282–286. 10.1016/j.scib.2018.02.006
119
LiF.WangJ.XiaG.LiZ. (2019). Negative differential thermal resistance through nanoscale solid-fluid-solid sandwiched structures. Nanoscale11, 13051–13057. 10.1039/c9nr01606f
120
LiD.GaoJ.ChengP.HeJ.YinY.HuY.et al (2020a). 2D boron sheets: structure, growth, and electronic and thermal transport properties. Adv. Funct. Mater.30, 1904349. 10.1002/adfm.201904349
121
LiD.GongY.ChenY.LinJ.KhanQ.ZhangY.et al (2020b). Recent progress of two-dimensional thermoelectric materials. Nano-Micro Lett.12, 36. 10.1007/s40820-020-0374-x
122
LiangT.ZhangP.YuanP.ZhaiS. (2018). In-plane thermal transport in black phosphorene/graphene layered heterostructures: a molecular dynamics study. Phys. Chem. Chem. Phys.20, 21151–21162. 10.1039/c8cp02831a
123
LiangT.ZhouM.ZhangP.YuanP.YangD. (2020). Multilayer in-plane graphene/hexagonal boron nitride heterostructures: insights into the interfacial thermal transport properties. Int. J. Heat Mass Tran.151, 119395. 10.1016/j.ijheatmasstransfer.2020.119395
124
LindsayL.BroidoD. A. (2011). Enhanced thermal conductivity and isotope effect in single-layer hexagonal boron nitride. Phys. Rev. B84, 155421. 10.1103/physrevb.84.155421
125
LindsayL.BroidoD. A.MingoN. (2009). Lattice thermal conductivity of single-walled carbon nanotubes: beyond the relaxation time approximation and phonon-phonon scattering selection rules. Phys. Rev. B80, 125407. 10.1103/physrevb.80.125407
126
LindsayL.BroidoD. A.MingoN. (2010). Flexural phonons and thermal transport in graphene. Phys. Rev. B82, 115427. 10.1103/physrevb.82.115427
127
LindsayL.LiW.CarreteJ.MingoN.BroidoD. A.ReineckeT. L. (2014). Phonon thermal transport in strained and unstrained graphene from first principles. Phys. Rev. B89, 155426. 10.1103/physrevb.89.155426
128
LiuR.LiW. (2018c). High-thermal-stability and high-thermal-conductivity Ti3C2Tx MXene/poly (vinyl alcohol) (PVA) composites. ACS Omega3, 2609–2617. 10.1021/acsomega.7b02001
129
LiuH.GaoJ.ZhaoJ. (2013a). Silicene on substrates: a way to preserve or tune its electronic properties. J. Phys. Chem. C117, 10353–10359. 10.1021/jp311836m
130
LiuZ.MaL.ShiG.ZhouW.GongY.LeiS.et al (2013b). In-plane heterostructures of graphene and hexagonal boron nitride with controlled domain sizes. Nat. Nanotechnol.8, 119. 10.1038/nnano.2012.256
131
LiuB.BaimovaJ. A.ReddyC. D.DmitrievS. V.LawW. K.FengX. Q.et al (2014a). Interface thermal conductance and rectification in hybrid graphene/silicene monolayer. Carbon79, 236–244. 10.1016/j.carbon.2014.07.064
132
LiuB.ReddyC. D.JiangJ.ZhuH.BaimovaJ. A.DmitrievS. V. (2014b). Thermal conductivity of silicene nanosheets and the effect of isotopic doping. J. Phys. D Appl. Phys.47, 165301. 10.1088/0022-3727/47/16/165301
133
LiuH.NealA. T.ZhuZ.LuoZ.XuX.TománekD.et al (2014c). Phosphorene: an unexplored 2D semiconductor with a high hole mobility. ACS Nano8, 4033–4041. 10.1021/nn501226z
134
LiuL.ParkJ.SiegelD. A.McCartyK. F.ClarkK. W.DengW.et al (2014d). Heteroepitaxial growth of two-dimensional hexagonal boron nitride templated by graphene edges. Science343, 163–167. 10.1126/science.1246137
135
LiuY.-Y.ZhouW.-X.TangL.-M.ChenK.-Q. (2014e). An important mechanism for thermal rectification in graded nanowires. Appl. Phys. Lett.105, 203111. 10.1063/1.4902427
136
LiuX.ZhangG.ZhangY.-W. (2015). Graphene-based thermal modulators. Nano Res.8, 2755–2762. 10.1007/s12274-015-0782-2
137
LiuX.ZhangG.ZhangY.-W. (2016). Topological defects at the graphene/h-BN interface abnormally enhance its thermal conductance. Nano Lett.16, 4954–4959. 10.1021/acs.nanolett.6b01565
138
LiuX.GaoJ.ZhangG.ZhangY.-W. (2017a). MoS2-graphene in-plane contact for high interfacial thermal conduction. Nano Res.10, 2944–2953. 10.1007/s12274-017-1504-8
139
LiuY.-Y.LiB.-L.ChenS.-Z.JiangX.ChenK.-Q. (2017b). Effect of room temperature lattice vibration on the electron transport in graphene nanoribbons. Appl. Phys. Lett.111, 133107. 10.1063/1.4999127
140
LiuY.OngZ. Y.WuJ.ZhaoY.WatanabeK.TaniguchiT.et al (2017c). Thermal conductance of the 2D MoS2/h-BN and graphene/h-BN interfaces. Sci. Rep.7, 43886. 10.1038/srep43886
141
LiuX.GaoJ.ZhangG.ZhangY.-W. (2018a). Design of phosphorene/graphene heterojunctions for high and tunable interfacial thermal conductance. Nanoscale, 10, 19854–19862. 10.1039/c8nr06110f
142
LiuY.-Y.ZengY.-J.JiaP.-Z.CaoX.-H.JiangX.ChenK.-Q. (2018b). An efficient mechanism for enhancing the thermoelectricity of nanoribbons by blocking phonon transport in 2D materials. J. Phys. Condens. Matter30, 275701. 10.1088/1361-648x/aac7f5
143
LiuF.ZouR.HuN.NingH.YanC.LiuY.et al (2019a). Enhancement of thermal energy transport across the graphene/h-BN heterostructure interface. Nanoscale11, 4067–4072. 10.1039/c8nr10468a
144
LiuY.XuY.DuanW. (2019b). Three-dimensional topological states of phonons with tunable pseudospin physics. Research2019, 5173580. 10.34133/2019/5173580
145
LuoZ.MaassenJ.DengY.DuY.GarreltsR. P.LundstromM. S.et al (2015). Anisotropic in-plane thermal conductivity observed in few-layer black phosphorus. Nat. Commun.6, 8572. 10.1038/ncomms9572
146
MaZ.HuZ.ZhaoX.TangQ.WuD.ZhouZ.et al (2014). Tunable band structures of heterostructured bilayers with transition-metal dichalcogenide and MXene monolayer. J. Phys. Chem. C118, 5593–5599. 10.1021/jp500861n
147
MaD.DingH.MengH.FengL.WuY.ShiomiJ.et al (2016). Nano-cross-junction effect on phonon transport in silicon nanowire cages. Phys. Rev. B94, 165434. 10.1103/physrevb.94.165434
148
MaD.WanX.YangN. (2018). Unexpected thermal conductivity enhancement in pillared graphene nanoribbon with isotopic resonance. Phys. Rev. B98, 245420. 10.1103/physrevb.98.245420
149
MaJ.-J.ZhengJ.-J.ZhuX.-L.LiuP.-F.LiW.-D.WangB.-T. (2019). First-principles calculations of thermal transport properties in MoS2/MoSe2 bilayer heterostructure. Phys. Chem. Chem. Phys.21, 10442–10448. 10.1039/c9cp01702j
150
MachidaY.MatsumotoN.IsonoT.BehniaK. (2020). Phonon hydrodynamics and ultrahigh-room-temperature thermal conductivity in thin graphite. Science367, 309–312. 10.1126/science.aaz8043
151
MahapatraP. S.SarkarK.KrishnamurthyH. R.MukerjeeS.GhoshA. (2017). Seebeck coefficient of a single van der Waals junction in twisted bilayer graphene. Nano Lett.17, 6822–6827. 10.1021/acs.nanolett.7b03097
152
MalekpourH.RamnaniP.SrinivasanS.BalasubramanianG.NikaD. L.MulchandaniA.et al (2016). Thermal conductivity of graphene with defects induced by electron beam irradiation. Nanoscale8, 14608–14616. 10.1039/c6nr03470e
153
MandelliD.LevenI.HodO.UrbakhM. (2017). Sliding friction of graphene/hexagonal–boron nitride heterojunctions: a route to robust superlubricity. Sci. Rep.7, 10851. 10.1038/s41598-017-10522-8
154
MannixA. J.ZhouX.-F.KiralyB.WoodJ. D.AlducinD.MyersB. D.et al (2015). Synthesis of borophenes: anisotropic, two-dimensional boron polymorphs. Science350, 1513–1516. 10.1126/science.aad1080
155
MazzamutoF.NguyenV. H.ApertetY.CaërC.ChassatC.Saint-MartinJ.et al (2011). Enhanced thermoelectric properties in graphene nanoribbons by resonant tunneling of electrons. Phys. Rev. B83, 235426. 10.1103/physrevb.83.235426
156
MelisC.BarbarinoG.ColomboL. (2015). Exploiting hydrogenation for thermal rectification in graphene nanoribbons. Phys. Rev. B92, 245408. 10.1103/physrevb.92.245408
157
MengX.PandeyT.JeongJ.FuS.YangJ.ChenK.et al (2019). Thermal conductivity enhancement in MoS2 under extreme strain. Phys. Rev. Lett.122, 155901. 10.1103/physrevlett.122.155901
158
MingoN.EsfarjaniK.BroidoD. A.StewartD. A. (2010). Cluster scattering effects on phonon conduction in graphene. Phys. Rev. B81 (4), 045408. 10.1103/physrevb.81.045408
159
MorookaM.YamamotoT.WatanabeK. (2008). Defect-induced circulating thermal current in graphene with nanosized width. Phys. Rev. B77, 033412.
160
MortazaviB.LeM.-Q.RabczukT.PereiraL. F. C. (2017). Anomalous strain effect on the thermal conductivity of borophene: a reactive molecular dynamics study. Phys. E93, 202–207. 10.1016/j.physe.2017.06.012
161
MuX.WuX.ZhangT.GoD. B.LuoT. (2014). Thermal transport in graphene oxide–from ballistic extreme to amorphous limit. Sci. Rep.4, 3909. 10.1038/srep03909
162
MuX.ZhangT.GoD. B.LuoT. (2015). Coherent and incoherent phonon thermal transport in isotopically modified graphene superlattices. Carbon83, 208–216. 10.1016/j.carbon.2014.11.028
163
NarducciD. (2011). Do we really need high thermoelectric figures of merit? A critical appraisal to the power conversion efficiency of thermoelectric materials. Appl. Phys. Lett.99, 102104. 10.1063/1.3634018
164
NiY.ZhangH.HuS.WangH.VolzS.XiongS. (2019). Interface diffusion-induced phonon localization in two-dimensional lateral heterostructures. Int. J. Heat Mass Tran.144, 118608. 10.1016/j.ijheatmasstransfer.2019.118608
165
NikaD. L.BalandinA. A. (2012). Two-dimensional phonon transport in graphene. J. Phys.: condens. Matter24 (23) 233203. 10.1088/0953-8984/24/23/233203
166
NikaD. L.BalandinA. A. (2017). Phonons and thermal transport in graphene and graphene-based materials. Rep.Prog. Phys.80 (3), 036502. 10.1088/1361-6633/80/3/036502
167
NikaD. L.CocemasovA. I.BalandinA. A. (2014). Specific heat of twisted bilayer graphene: engineering phonons by atomic plane rotations. App. Phy. Lett.105, 031904. 10.1063/1.4890622
168
NobakhtA. Y.GandomiY. A.WangJ.BowmanM. H.MarableD. C.GarrisonB. E.et al (2018). Thermal rectification via asymmetric structural defects in graphene. Carbon132, 565–572. 10.1016/j.carbon.2018.02.087
169
NomuraM.ShiomiJ.ShigaT.AnufrievR. (2018). Thermal phonon engineering by tailored nanostructures. J. Appl. Phys.57, 080101. 10.1364/opj.2018.31acj4
170
OngZ. Y.PopE.ShiomiJ. (2011). Reduction of phonon lifetimes and thermal conductivity of a carbon nanotube on amorphous silica. Phys. Rev. B84, 165418. 10.1103/physrevb.84.165418
171
OngZ. Y.ZhangG.ZhangY. W. (2016). Controlling the thermal conductance of graphene/h−BN lateral interface with strain and structure engineering. Phys. Rev. B93 (7), 075406. 10.1103/physrevb.93.075406
172
OuyangT.ChenY.XieY.WeiX. L.YangK.YangP.et al (2010). Ballistic thermal rectification in asymmetric three-terminal graphene nanojunctions. Phys. Rev. B82, 245403. 10.1103/physrevb.82.245403
173
PakA. J.HwangG. S. (2016). Theoretical analysis of thermal transport in graphene supported on hexagonal boron nitride: the importance of strong adhesion due to π-bond polarization. Phys. Rev. Appl.6 (3), 034015. 10.1103/physrevapplied.6.034015
174
PalS.PuriI. K. (2014). Thermal rectification in a polymer-functionalized single-wall carbon nanotube. Nanotechnology25, 345401. 10.1088/0957-4484/25/34/345401
175
PalS.PuriI. K. (2015). Thermal and gate using a monolayer graphene nanoribbon. Small11, 2910–2917. 10.1002/smll.201303888
176
PeiQ.-X.ShaZ.-D.ZhangY.-W. (2011). A theoretical analysis of the thermal conductivity of hydrogenated graphene. Carbon49, 4752–4759. 10.1016/j.carbon.2011.06.083
177
PeiQ.-X.ZhangX.DingZ.ZhangY.-Y.ZhangY.-W. (2017). Thermal stability and thermal conductivity of phosphorene in phosphorene/graphene van der Waals heterostructures. Phys. Chem. Chem. Phys.19, 17180–17186. 10.1039/c7cp02553j
178
PeimyooN.ShangJ.YangW.WangY.CongC.YuT. (2015). Thermal conductivity determination of suspended mono-and bilayer WS2 by Raman spectroscopy. Nano Res.8, 1210–1221. 10.1007/s12274-014-0602-0
179
PengX.-F.ChenK.-Q. (2014). Thermal transport for flexural and in-plane phonons in graphene nanoribbons. Carbon77, 360–365. 10.1016/j.carbon.2014.05.039
180
PengX. F.ChenK. Q.WanQ.ZouB. S.DuanW. (2010). Quantized thermal conductance at low temperatures in quantum wire with catenoidal contacts. Phys. Rev. B81, 195317. 10.1103/physrevb.81.195317
181
PengX.-F.ZhouX.TanS.-H.WangX.-J.ChenL.-Q.ChenK.-Q. (2017). Thermal conductance in graphene nanoribbons modulated by defects and alternating boron-nitride structures. Carbon113, 334–339. 10.1016/j.carbon.2016.11.066
182
PettesM. T.JoI.YaoZ.ShiL. (2011). Influence of polymeric residue on the thermal conductivity of suspended bilayer graphene. Nano Lett.11, 1195–1200. 10.1021/nl104156y
183
QinG.HuM. (2018). Thermal transport in phosphorene. Small14, 1702465. 10.1002/smll.201702465
184
QinG.ZhangX.YueS. Y.QinZ.WangH.HanY.et al (2016). Resonant bonding driven giant phonon anharmonicity and low thermal conductivity of phosphorene. Phy. Rev. B94, 165445. 10.1103/physrevb.94.165445
185
QinH.PeiQ.-X.LiuY.ZhangY.-W. (2019). The mechanical and thermal properties of MoS2-WSe2 lateral heterostructures. Phys. Chem. Chem. Phys.21, 15845–15853. 10.1039/c9cp02499a
186
QiuB.RuanX. (2012a). Reduction of spectral phonon relaxation times from suspended to supported graphene. Appl. Phys. Lett.100, 193101. 10.1063/1.4712041
187
QiuB.WangY.ZhaoQ.RuanX. (2012b). The effects of diameter and chirality on the thermal transport in free-standing and supported carbon-nanotubes. Appl. Phys. Lett.100, 233105. 10.1063/1.4725194
188
RanjanP.SahuT. K.BhushanR.YamijalaS. S.LateD. J.KumarP.et al (2019). Freestanding borophene and its hybrids. Adv. Mater.31, 1900353. 10.1002/adma.201900353
189
RavichandranJ.YadavA. K.CheaitoR.RossenP. B.SoukiassianA.et al (2014). Crossover from incoherent to coherent phonon scattering in epitaxial oxide superlattices. Nat. Mater.13, 168–172. 10.1038/nmat3826
190
RosulM. G.LeeD.OlsonD. H.LiuN.WangX.HopkinsP. E.et al (2019). Thermionic transport across gold-graphene-WSe2 van der Waals heterostructures. Sci. Adv.5, eaax7827. 10.1126/sciadv.aax7827
191
RoweP.CsányiG.AlfèD.MichaelidesA. (2018). Development of a machine learning potential for graphene. Phys. Rev. B97 (5), 054303. 10.1103/PhysRevB.97.054303
192
RoyT.TosunM.CaoX.FangH.LienD.-H.ZhaoP.et al (2015). Dual-gated MoS2/WSe2 van der Waals tunnel diodes and transistors. ACS Nano9, 2071–2079. 10.1021/nn507278b
193
SääskilahtiK.OksanenJ.VolzS.TulkkiJ. (2015). Frequency-dependent phonon mean free path in carbon nanotubes from nonequilibrium molecular dynamics. Phys. Rev. B91, 115426. 10.1103/physrevb.91.115426
194
SadeghiM. M.JoI.ShiL. (2013). Phonon-interface scattering in multilayer graphene on an amorphous support. Proc. Natl. Acad. Sci. U. S. A110, 16321–16326. 10.1073/pnas.1306175110
195
SadeghiH.SangtarashS.LambertC. J. (2016). Cross-plane enhanced thermoelectricity and phonon suppression in graphene/MoS2 van der Waals heterostructures. 2D Mater4 (1), 015012. 10.1088/2053-1583/4/1/015012
196
SahooS.GaurA. P. S.AhmadiM.GuinelM. J.-F.KatiyarR. S. (2013). Temperature-dependent Raman studies and thermal conductivity of few-layer MoS2. J. Phys. Chem. C117, 9042–9047. 10.1021/jp402509w
197
SandonasL. M.Cuba-SupantaG.GutierrezR.DianatA.LandauroC. V.CunibertiG. (2017). Enhancement of thermal transport properties of asymmetric Graphene/h-BN nanoribbon heterojunctions by substrate engineering. Carbon124, 642–650. 10.1016/j.carbon.2017.09.025
198
SarikurtS.ÇakırD.KeçeliM.SevikC. (2018). The influence of surface functionalization on thermal transport and thermoelectric properties of MXene monolayers. Nanoscale10, 8859–8868. 10.1039/c7nr09144c
199
SchellingP. K.PhillpotS. R. (2003). Multiscale simulation of phonon transport in superlattices. J. Appl. Phys.93, 5377–5387. 10.1063/1.1561601
200
SellanD. P.LandryE. S.TurneyJ. E.McGaugheyA. J.AmonC. H. (2010). Size effects in molecular dynamics thermal conductivity predictions. Phys. Rev. B81, 214305. 10.1103/physrevb.81.214305
201
SeolJ. H.JoI.MooreA. L.LindsayL.AitkenZ. H.PettesM. T.et al (2010). Two-dimensional phonon transport in supported graphene. Science, 328, 213–216. 10.1126/science.1184014
202
SerranoJ.BosakA.ArenalR.KrischM.WatanabeK.TaniguchiT.et al (2007). Vibrational properties of hexagonal boron nitride: inelastic X-ray scattering and ab initio calculations. Phy. Rev. Lett.98 (9), 095503. 10.1103/PHYSREVLETT.98.095503
203
SevikC.KinaciA.HaskinsJ. B.ÇağınT. (2011). Characterization of thermal transport in low-dimensional boron nitride nanostructures. Phys. Rev. B, 84 (8), 085409. 10.1103/physrevb.84.085409
204
ShenX.LinX.JiaJ.WangZ.LiZ.KimJ.-K. (2014). Tunable thermal conductivities of graphene oxide by functionalization and tensile loading. Carbon80, 235–245. 10.1016/j.carbon.2014.08.062
205
SmitheK. K. H.EnglishC. D.SuryavanshiS. V.PopE. (2018). High-field transport and velocity saturation in synthetic monolayer MoS2. Nano Lett.18, 4516–4522. 10.1021/acs.nanolett.8b01692
206
SongJ.MedhekarN. V. (2013). Thermal transport in lattice-constrained 2D hybrid graphene heterostructures. J. Phys. Condens. Matter25, 445007. 10.1088/0953-8984/25/44/445007
207
SoodA.XiongF.ChenS.WangH.SelliD.ZhangJ.et al (2018). An electrochemical thermal transistor. Nat. Commun.9, 4510. 10.1038/s41467-018-06760-7
208
SuL.ZhangY. (2015). Temperature coefficients of phonon frequencies and thermal conductivity in thin black phosphorus layers. App. Phys. Lett.107 (7), 071905. 10.1063/1.4928931
209
SunH.LiQ.WanX. G. (2016). First-principles study of thermal properties of borophene. Phys. Chem. Chem. Phys.18, 14927–14932. 10.1039/c6cp02029a
210
SutterP.LahiriJ.AlbrechtP.SutterE. (2011). Chemical vapor deposition and etching of high-quality monolayer hexagonal boron nitride films. ACS Nano5, 7303–7309. 10.1021/nn202141k
211
TaubeA.JudekJ.ŁapińskaA.ZdrojekM. (2015). Temperature-dependent thermal properties of supported MoS2 monolayers. ACS Appl. Mater. Interfaces7, 5061–5065. 10.1021/acsami.5b00690
212
ThomasJ. A.TurneyJ. E.IutziR. M.AmonC. H.McGaugheyA. J. (2010). Predicting phonon dispersion relations and lifetimes from the spectral energy density. Phys. Rev. B81, 081411. 10.1115/ihtc14-22262
213
TianH.XieD.YangY.RenT. L.ZhangG.WangY. F.et al (2012a). A novel solid-state thermal rectifier based on reduced graphene oxide. Sci. Rep.2, 523. 10.1038/srep00523
214
TianZ.EsfarjaniK.ChenG. (2012b). Enhancing phonon transmission across a Si/Ge interface by atomic roughness: first-principles study with the Green's function method. Phys. Rev. B86, 235304. 10.1103/physrevb.86.235304
215
TienD. H.ParkJ.-Y.KimK. B.LeeN.ChoiT.KimP.et al (2016). Study of graphene-based 2D-heterostructure device fabricated by all-dry transfer process. ACS Appl. Mater. Interfaces8, 3072–3078. 10.1021/acsami.5b10370
216
TranV.-T.Saint-MartinJ.DollfusP. (2015). High thermoelectric performance in graphene nanoribbons by graphene/BN interface engineering. Nanotechnology26 (49), 495202. 10.1088/0957-4484/26/49/495202
217
VarshneyV.LeeJ.BrownJ. S.FarmerB. L.VoevodinA. A.RoyA. K. (2018). Effect of length, Diameter, chirality, Deformation, and strain on contact Thermal conductance between single-Wall carbon nanotubes. Front. Mater5, 17. 10.3389/fmats.2018.00017
218
VogtP.De PadovaP.QuaresimaC.AvilaJ.FrantzeskakisE.AsensioM. C.et al (2012). Silicene: compelling experimental evidence for graphenelike two-dimensional silicon. Phys. Rev. Lett.108, 155501. 10.1103/physrevlett.108.155501
219
WangL.MericI.HuangP. Y.GaoQ.GaoY.TranH.et al (2013a). One-dimensional electrical contact to a two-dimensional material. Science342, 614–617. 10.1126/science.1244358
220
WangX.HuangT.LuS. (2013b). High performance of the thermal transport in graphene supported on hexagonal boron nitride. Appl. Phys. Express6 (7), 075202. 10.7567/apex.6.075202
221
WangY.VallabhaneniA.HuJ.QiuB.ChenY. P.RuanX. (2014). Phonon lateral confinement enables thermal rectification in asymmetric single-material nanostructures. Nano Lett.14, 592–596. 10.1021/nl403773f
222
WangY.ZhangK.XieG. (2016). Remarkable suppression of thermal conductivity by point defects in MoS2 nanoribbons. Appl. Surf. Sci.360, 107–112. 10.1016/j.apsusc.2015.10.235
223
WangH.HuS.TakahashiK.ZhangX.TakamatsuH.ChenJ. (2017a). Experimental study of thermal rectification in suspended monolayer graphene. Nat. Commun.8, 15843. 10.1038/ncomms15843
224
WangX.HongY.ChanP. K. L.ZhangJ. (2017b). Phonon thermal transport in silicene-germanene superlattice: a molecular dynamics study. Nanotechnology28, 255403. 10.1088/1361-6528/aa71fa
225
WangX.WangM.HongY.WangZ.ZhangJ. (2017c). Coherent and incoherent phonon transport in a graphene and nitrogenated holey graphene superlattice. Phys. Chem. Chem. Phys.19, 24240–24248. 10.1039/c7cp04219a
226
WangY.XuN.LiD.ZhuJ. (2017d). Thermal properties of two dimensional layered materials. Adv. Funct. Mater.27, 1604134. 10.1002/adfm.201604134
227
WardA.BroidoD. A. (2010). Intrinsic phonon relaxation times from first-principles studies of the thermal conductivities of Si and Ge. Phys. Rev. B81 (8) 085205. 10.1103/physrevb.81.085205
228
WeiN.XuL.WangH.-Q.ZhengJ.-C. (2011). Strain engineering of thermal conductivity in graphene sheets and nanoribbons: a demonstration of magic flexibility. Nanotechnology22, 105705. 10.1088/0957-4484/22/10/105705
229
WeiX.WangY.ShenY.XieG.XiaoH.ZhongJ.et al (2014). Phonon thermal conductivity of monolayer MoS2: a comparison with single layer graphene. Appl. Phys. Lett.105, 103902. 10.1063/1.4895344
230
WeiA.LahkarS.LiX.LiS.YeH. (2019). Multilayer graphene-based thermal rectifier with interlayer gradient functionalization. ACS Appl. Mater. Interfaces11, 45180–45188. 10.1021/acsami.9b11762
231
WuS. C.ShanG.YanB. (2014). Prediction of near-room-temperature quantum anomalous Hall effect on honeycomb materials. Phys. Rev. Lett.113, 256401. 10.1103/physrevlett.113.256401
232
WuX.VarshneyV.LeeJ.ZhangT.WohlwendJ. L.RoyA. K.et al (2016). Hydrogenation of penta-graphene leads to unexpected large improvement in thermal conductivity. Nano Lett.16, 3925–3935. 10.1021/acs.nanolett.6b01536
233
WuD.CaoX.-H.ChenS.-Z.TangL.-M.FengY.-X.ChenK.-Q.et al (2019). Pure spin current generated in thermally driven molecular magnetic junctions: a promising mechanism for thermoelectric conversion. J. Mater. Chem. A7, 19037–19044. 10.1039/c9ta04642a
234
WuD.CaoX. H.JiaP. Z.ZengY. J.FengY. X.TangL. M.et al (2020). Excellent thermoelectric performance in weak-coupling molecular junctions with electrode doping and electrochemical gating. Sci. China Phys. Mech63, 276811. 10.1007/s11433-019-1528-y
235
XieZ. X.TangL. M.PanC. N.LiK. M.ChenK. Q.DuanW. (2012). Enhancement of thermoelectric properties in graphene nanoribbons modulated with stub structures. Appl. Phys. Lett.100, 073105. 10.1103/physrevb.93.075404
236
XieG.GuoY.LiB.YangL.ZhangK.TangM.et al (2013). Phonon surface scattering controlled length dependence of thermal conductivity of silicon nanowires. Phys. Chem. Chem. Phys.15, 14647–14652. 10.1039/c3cp50969a
237
XieH.HuM.BaoH. (2014a). Thermal conductivity of silicene from first-principles. Appl. Phys. Lett.104, 131906. 10.1063/1.4870586
238
XieZ.-X.ZhangY.YuX.LiK.-M.ChenQ. (2014b). Ballistic thermal conductance by phonons through superlattice quantum-waveguides. J. Appl. Phys.115, 104309. 10.1063/1.4868595
239
XieH.OuyangT.GermaneauÉ.QinG.HuM.BaoH. (2016). Large tunability of lattice thermal conductivity of monolayer silicene via mechanical strain. Phys. Rev. B93, 075404. 10.1103/physrevb.93.075404
240
XieF.FanZ.-Q.ChenK.-Q.ZhangX.-J.LongM.-Q. (2017). Influence of anchoring groups on single-molecular junction conductance: theoretical comparative study of thiol and amine. Org. Eelctron50, 198–203. 10.1016/j.orgel.2017.07.052
241
XieG.DingD.ZhangG. (2018a). Phonon coherence and its effect on thermal conductivity of nanostructures. Adv. Phys. X3, 1480417. 10.1080/23746149.2018.1480417
242
XieG.JuZ.ZhouK.WeiX.GuoZ.CaiY.et al (2018b). Ultra-low thermal conductivity of two-dimensional phononic crystals in the incoherent regime. Npj Comput. Mater4, 21. 10.1038/s41524-018-0076-9
243
XieZ.-X.YuX.ChenX.-K.ZhouW.-X.ShiY.-M.ZhangL.-F. (2020). Modulation of thermal transport in AlxGa1−xAs alloy nanowires with varying compositions. Appl. Phys. Lett.116 (14), 143102. 10.1063/5.0003961
244
XiongS.SääskilahtiK.KosevichY. A.HanH.DonadioD.VolzS. (2016). Blocking phonon transport by structural resonances in alloy-based nanophononic metamaterials leads to ultralow thermal conductivity. Phys. Rev. Lett.117 (2), 025503. 10.1103/physrevlett.117.025503
245
XiongG.XingY.ZhangL. (2018). Interfacial thermal transport via one-dimensional atomic junction model. Front. Energy Res.6, 6. 10.3389/fenrg.2018.00006
246
XuW.ZhangG. (2016a). Remarkable reduction of thermal conductivity in phosphorene phononic crystal. J. Phys. Condens. Matter28, 175401. 10.1088/0953-8984/28/17/175401
247
XuY.ChenX.GuB.-L.DuanW. (2009). Intrinsic anisotropy of thermal conductance in graphene nanoribbons. Appl. Phys. Lett.95, 233116. 10.1063/1.3272678
248
XuW.ZhangG.LiB. (2014a). Interfacial thermal resistance and thermal rectification between suspended and encased single layer graphene. J. Appl. Phys.116, 134303. 10.1063/1.4896733
249
XuX.PereiraL. F.WangY.WuJ.ZhangK.ZhaoX.et al (2014b). Length-dependent thermal conductivity in suspended single-layer graphene. Nat. Commun.5, 3689. 10.1038/ncomms4689
250
XuY.LiZ.DuanW. (2014c). Thermal and thermoelectric properties of graphene. Small10, 2182–2199. 10.1002/smll.201303701
251
XuW.ZhuL.CaiY.ZhangG.LiB. (2015). Direction dependent thermal conductivity of monolayer phosphorene: parameterization of Stillinger-Weber potential and molecular dynamics study. J. Appl. Phys.117, 214308. 10.1063/1.4922118
252
XuX.ChenJ.LiB. (2016b). Phonon thermal conduction in novel 2D materials. J. Phys. Condens. Matter28, 483001. 10.1088/0953-8984/28/48/483001
253
YanZ.JiangC.PopeT. R.TsangC. F.StickneyJ. L.GoliP.et al (2013). Phonon and thermal properties of exfoliated TaSe2 thin films. J. Appl. Phys.114, 204301. 10.1063/1.4833250
254
YanZ.ChenL.YoonM.KumarS. (2016a). Phonon transport at the interfaces of vertically stacked graphene and hexagonal boron nitride heterostructures. Nanoscale8, 4037–4046. 10.1039/c5nr06818e
255
YanZ.ChenL.YoonM.KumarS. (2016b). The role of interfacial electronic properties on phonon transport in two-dimensional MoS2 on metal substrates. ACS Appl. Mater. Interfaces8, 33299–33306. 10.1021/acsami.6b10608
256
YangX.YuD.CaoB. (2017). Giant thermal rectification from single-carbon nanotube-graphene junction. ACS Appl. Mater. Interfaces9, 24078–24084. 10.1021/acsami.7b04464
257
YangN.ZengX.LuJ.SunR.WongC.-P. (2018). Effect of chemical functionalization on the thermal conductivity of 2D hexagonal boron nitride. Appl. Phys. Lett.113, 171904. 10.1063/1.5050293
258
YankowitzM.MaQ.Jarillo-HerreroP.LeRoyB. J. (2019). van der Waals heterostructures combining graphene and hexagonal boron nitride. Nat. Rev. Phys.1, 112–125. 10.1038/s42254-018-0016-0
259
YasaeiP.HemmatZ.FossC. J.LiS. J.HongL.BehranginiaA.et al (2018). Enhanced thermal boundary conductance in few-layer Ti3C2 MXene with encapsulation. Adv. Mater.30, 1801629. 10.1002/adma.201801629
260
YokomizoY.NakamuraJ. (2013). Giant Seebeck coefficient of the graphene/h-BN superlattices. Appl. Phys. Lett.103, 113901. 10.1063/1.4820820
261
YuanK.SunM.WangZ.TangD. (2015). Tunable thermal rectification in silicon-functionalized graphene nanoribbons by molecular dynamics simulation. Int. J. Therm. Sci.98, 24–31. 10.1016/j.ijthermalsci.2015.07.004
262
ZengJ.ChenK.-Q.TongY.-X. (2018a). Covalent coupling of porphines to graphene edges: quantum transport properties and their applications in electronics. Carbon127, 611–617. 10.1016/j.carbon.2017.11.047
263
ZengY. J.LiuY. Y.ZhouW. X.ChenK. Q. (2018b). Nanoscale thermal transport: theoretical method and application. Chin. Phys. B27 (3), 036304. 10.1088/1674-1056/27/3/036304
264
ZengY.-J.WuD.CaoX.-H.FengY.-X.TangL.-M.ChenK.-Q. (2020a). Significantly enhanced thermoelectric performance of molecular junction by twist angle dependent phonon interference effect. J. Mater. Chem. A8, 11884–11891. 10.1039/d0ta02423f
265
ZengY. J.WuD.CaoX. H.ZhouW. X.TangL. M.ChenK. Q. (2020b). Nanoscale organic thermoelectric materials: measurement, theoretical models, and optimization strategies. Adv. Funct. Mater.30, 1903873. 10.1002/adfm.201903873
266
ZengY.LoC.-L.ZhangS.ChenZ.MarconnetA. (2020c). Dynamically tunable thermal transport in polycrystalline graphene by strain engineering. Carbon158, 63–68. 10.1016/j.carbon.2019.11.060
267
ZhangT.LuoT. (2015b). Thermal diodes: giant thermal rectification from polyethylene nanofiber thermal diodes. Small11, 4656. 10.1002/smll.201570222
268
ZhangG.ZhangH. (2011). Thermal conduction and rectification in few-layer graphene Y Junctions. Nanoscale3, 4604–4607. 10.1039/c1nr10945f
269
ZhangG.ZhangY.-W. (2017a). Thermoelectric properties of two-dimensional transition metal dichalcogenides. J. Mater. Chem. C5, 7684–7698. 10.1039/c7tc01088e
270
ZhangY.ChangT.-R.ZhouB.CuiY.-T.YanH.LiuZ.et al (2013). Direct observation of the transition from indirect to direct bandgap in atomically thin epitaxial MoSe2. Nat. Nanotechnol.9, 111–115. 10.1038/nnano.2013.277
271
ZhangX.XieH.HuM.BaoH.YueS.QinG.et al (2014a). Thermal conductivity of silicene calculated using an optimized Stillinger-Weber potential. Phy. Rev. B89, 054310. 10.1103/physrevb.89.054310
272
ZhangY.XieZ.-X.YuX.WangH.-B.LiK.-M. (2014b). Ballistic thermal transport in a cylindrical semiconductor nanowire modulated with bridge contacts. Journal of Applied Physics116, 144304. 10.1063/1.4897548
273
ZhangJ.HongY.YueY. (2015a). Thermal transport across graphene and single layer hexagonal boron nitride. J. Appl. Phys.117, 134307. 10.1063/1.4916985
274
ZhangX.SunD.LiY.LeeG.-H.CuiX.ChenetD.et al (2015c). Measurement of lateral and interfacial thermal conductivity of single-and bilayer MoS2 and MoSe2 using refined optothermal Raman technique. ACS Appl. Mater. Interfaces7, 25923–25929. 10.1021/acsami.5b08580
275
ZhangY.-Y.PeiQ.-X.MaiY.-W.LaiS.-K. (2016). Interfacial thermal conductance in multilayer graphene/phosphorene heterostructure. J. Phys. D Appl. Phys.49, 465301. 10.1088/0022-3727/49/46/465301
276
ZhangJ.HongY.WangX.YueY.XieD.JiangJ.et al (2017b). Phonon thermal properties of transition-metal dichalcogenides MoS2 and MoSe2 heterostructure. J. Phys. Chem. C121, 10336–10344. 10.1021/acs.jpcc.7b02547
277
ZhangZ.HuS.ChenJ.LiB. (2017c). Hexagonal boron nitride: a promising substrate for graphene with high heat dissipation. Nanotechnology28, 225704. 10.1088/1361-6528/aa6e49
278
ZhangY.PeiQ.-X.WangC.-M.YangC.ZhangY.-W. (2018). Interfacial thermal conductance and thermal rectification of hexagonal BCnN/graphene in-plane heterojunctions. J. Phys. Chem. C122, 22783–22789. 10.1021/acs.jpcc.8b08015
279
ZhangY.ChenS.-Z.XieZ.-X.YuX.DengY.-X.ChenX.-H.et al (2019). Spin-resolved transport properties of DNA base multi-functional electronic devices. Phys. Lett.383, 2069–2075. 10.1016/j.physleta.2019.03.034
280
ZhaoG.ChengY.WuY.XuX.HaoX. (2018a). New 2D carbon nitride organic materials synthesis with huge-application prospects in CN photocatalyst. Small14, 1704138. 10.1002/smll.201704138
281
ZhaoY.ZhangG.NaiM. H.DingG.LiD.LiuY.et al (2018b). Probing the physical origin of anisotropic thermal transport in black phosphorus nanoribbons. Adv. Mater.30, 1804928. 10.1002/adma.201804928
282
ZhaoG.WangA.HeW.XingY.XuX. (2019). 2D new nonmetal photocatalyst of sulfur‐doped h‐BN nanosheeets with high photocatalytic activity. Adv. Mater. Interfaces6, 1900062. 10.1002/admi.201900062
283
ZhaoG.ChengY.SunP.MaW.HaoS.WangX.et al (2020a). Biocarbon based template synthesis of uniform lamellar MoS2 nanoflowers with excellent energy storage performance in lithium-ion battery and supercapacitors. Electrochim. Acta331, 135262. 10.1016/j.electacta.2019.135262
284
ZhaoY.CaiY.ZhangL.LiB.ZhangG.ThongJ. T. L. (2020b). Thermal transport in 2D semiconductors-considerations for device applications. Adv. Funct. Mater.30, 1903929. 10.1002/adfm.201903929
285
ZhaoG.HaoS.GuoJ.XingY.ZhangL.XuX. (2021). Design of p-n homojunctions in metal-free carbon nitride photocatalyst for overall water splitting. Chin. J. Catal.42, 501–509. 10.1016/s1872-2067(20)63670-1
286
ZhouW.-X.ChenK.-Q. (2015). Enhancement of thermoelectric performance in β-graphyne nanoribbons by suppressing phononic thermal conductance. Carbon85, 24–27. 10.1016/j.carbon.2014.12.059
287
ZhouY.HuM. (2017c). Full quantification of frequency-dependent interfacial thermal conductance contributed by two- and three-phonon scattering processes from nonequilibrium molecular dynamics simulations. Phys. Rev. B95, 115313. 10.1103/physrevb.95.115313
288
ZhouH.ZhuJ.LiuZ.YanZ.FanX.LinJ.et al (2014). High thermal conductivity of suspended few-layer hexagonal boron nitride sheets. Nano Res.7, 1232–1240. 10.1007/s12274-014-0486-z
289
ZhouH.CaiY.ZhangG.ZhangY. W. (2017a). Superior lattice thermal conductance of single-layer borophene. Npj 2D Mater. Appl.1, 14. 10.1038/s41699-017-0018-2
290
ZhouY.-H.ZhangX.JiC.LiuZ.-M.ChenK.-Q. (2017b). The length and hydrogenation effects on electronic transport properties of carbon-based molecular wires. Org. Electron.51, 332–340. 10.1016/j.orgel.2017.09.031
291
ZhouY.ZhengX.ChengZ.-Q.ChenK.-Q. (2019). Current superposition law realized in molecular devices connected in parallel. J. Phys. Chem. C123, 10462–10468. 10.1021/acs.jpcc.9b01812
292
ZhouW. X.ChengY.ChenK. Q.XieG.WangT.ZhangG. (2020a). Thermal conductivity of amorphous materials. Adv. Funct. Mater.30, 1903829. 10.1002/adfm.201903829
293
ZhouW.-X.WuD.XieG.ChenK.-Q.ZhangG. (2020b). α-Ag2S: a ductile thermoelectric material with high ZT. ACS Omega5, 5796–5804. 10.1021/acsomega.9b03929
294
ZhuT.ErtekinE. (2014b). Phonon transport on two-dimensional graphene/boron nitride superlattices. Phys. Rev. B90, 195209. 10.1103/physrevb.90.195209
295
ZhuT.ErtekinE. (2015). Resolving anomalous strain effects on two-dimensional phonon flows: the cases of graphene, boron nitride, and planar superlattices. Phys. Rev. B91, 205429. 10.1103/physrevb.91.205429
296
ZhuL.ZhangG.LiB. (2014a). Coexistence of size-dependent and size-independent thermal conductivities in phosphorene. Phy. Rev. B90, 214302. 10.1103/physrevb.90.214302
297
ZhuW.LiangL.RobertsR. H.LinJ.-F.AkinwandeD. (2018). Anisotropic electron–phonon interactions in angle-resolved Raman study of strained black phosphorus. ACS Nano12, 12512–12522. 10.1021/acsnano.8b06940
298
ZhuX.-L.LiuP.-F.ZhangJ.ZhangP.ZhouW.-X.XieG.et al (2019). Monolayer SnP3: an excellent p-type thermoelectric material. Nanoscale11, 19923–19932. 10.1039/c9nr04726c
299
ZouJ. H.CaoB. Y. (2017). Phonon thermal properties of graphene on h-BN from molecular dynamics simulations. Appl. Phys. Lett.110 (10), 103106. 10.1063/1.4978434
Summary
Keywords
two-dimensional materials, heterostructures, thermal conductivity, phonon property, thermoelectric
Citation
Chen X-K, Zeng Y-J and Chen K-Q (2020) Thermal Transport in Two-Dimensional Heterostructures. Front. Mater. 7:578791. doi: 10.3389/fmats.2020.578791
Received
01 July 2020
Accepted
09 November 2020
Published
07 December 2020
Volume
7 - 2020
Edited by
Weifeng Li, Shandong University, China
Reviewed by
Gang Zhao, University of Jinan, China
Yang Zhao, Western University, Canada
Updates
Copyright
© 2020 Chen, Chen and Zeng.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ke-Qiu Chen, keqiuchen@hnu.edu.cn
This article was submitted to Energy Materials, a section of the journal Frontiers in Materials
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.