Search

Plastic crystals, owing to their exceptional properties, are gradually finding applications in solid-state refrigeration and ferroelectric fields. However, their inherently low thermal conductivity restricts their utilization in electronic devices. This study demonstrates that applying megapascal pressure of hot pressing can enhance the thermal conductivity of plastic crystal films. Most importantly, it induces significant anisotropy in thermal conductivity. Such anisotropy in thermal conductivity is beneficial for specialized thermal management applications, such as directing heat flow paths in electronic devices. In this study, [(CH3)4N][FeCl4] PCs films were prepared by hot pressing. At a pressure of 16 MPa, the ratio of in-plane to cross-plane thermal conductivity in the film reaches a remarkable 5.5. This is attributed to the preferential orientation along the (002) crystal plane induced by uniaxial pressure, leading to the formation of a layered structure and the creation of a flat and dense film. Furthermore, according to molecular dynamics simulations, the thermal conductivity along the [100] and [010] directions (parallel to the (002) crystal plane) is higher than in other directions. Therefore, significant modulation of anisotropy in thermal conductivity is achieved in [(CH3)4N][FeCl4] films by applying uniaxial hot pressing pressure. This phenomenon has the potential to greatly broaden the application of plastic crystals in the field of flexible electronic devices.

The influence of hydrated cation-{\pi} interaction forces on the adsorption and filtration capabilities of graphene-based membrane materials is significant. However, the lack of interaction potential between hydrated Cs+ and graphene limits the scope of adsorption studies. Here, it is developed that a deep neural network potential function model to predict the interaction force between hydrated Cs+ and graphene. The deep potential has DFT-level accuracy, enabling accurate property prediction. This deep potential is employed to investigate the properties of the graphene surface solution, including the density distribution, mean square displacement, and vibrational power spectrum of water. Furthermore, calculations of the molecular orbital electron distributions indicate the presence of electron migration in the molecular orbitals of graphene and hydrated Cs+, resulting in a strong electrostatic interaction force. The method provides a powerful tool to study the adsorption behavior of hydrated cations on graphene surfaces and offers a new solution for handling radionuclides.

Modulation of thermal conductivity has become a hotspot in the field of heat conduction. A novel strategy based on targeted phonon excitation has been recently proposed for efficient and reversible modulation of thermal conductivity. In this article, the effectiveness of that strategy is further evaluated on hexagonal boron nitride through ab initio methods. Results indicate that thermal conductivity can be increased from 885 W m-1 K-1 to 1151 W m-1 K-1 or decreased to 356 W m-1 K-1, thereby broadening the scope of applicability of this strategy., Comment: 12 pages, 3 figures

Machine learning's application in solar-thermal desalination is limited by data shortage and inconsistent analysis. This study develops an optimized dataset collection and analysis process for the representative solar still. By ultra-hydrophilic treatment on the condensation cover, the dataset collection process reduces the collection time by 83.3%. Over 1,000 datasets are collected, which is nearly one order of magnitude larger than up-to-date works. Then, a new interdisciplinary process flow is proposed. Some meaningful results are obtained that were not addressed by previous studies. It is found that Radom Forest might be a better choice for datasets larger than 1,000 due to both high accuracy and fast speed. Besides, the dataset range affects the quantified importance (weighted value) of factors significantly, with up to a 115% increment. Moreover, the results show that machine learning has a high accuracy on the extrapolation prediction of productivity, where the minimum mean relative prediction error is just around 4%. The results of this work not only show the necessity of the dataset characteristics' effect but also provide a standard process for studying solar-thermal desalination by machine learning, which would pave the way for interdisciplinary study.

The thermal transport across inorganic/organic interfaces attracts interest for both academic and industry due to its widely applications in flexible electronics etc. Here, the interfacial thermal conductance of inorganic/organic interfaces consisting of silicon and polyvinylidene fluoride is systematically investigated by molecular dynamics simulations. Interestingly, it is demonstrated that a modified silicon surface with hydroxyl groups can drastically enhance the conductance by 698%. These results are elucidated based on interfacial couplings and lattice dynamics insights. This study not only provides feasible strategies to effectively modulate the interfacial thermal conductance of inorganic/organic interfaces but also deepens the understanding of the fundamental physics underlying phonon transport across interfaces.

The thermal management have been widely focused due to broad applications. Generally, the deformation can largely tune the thermal transport. The main challenge of flexible electronics/ materials is to maintain thermal conductance under large deformation. This work investigates the thermal conductance of a nano-designed Si metamaterial constructed with curved nanobeams by molecular dynamics simulation. Interestingly, it shows that the thermal conductance of the nano-designed Si metamaterial is insensitive under a large deformation (strain~-41%). The new feature comes from the designed curved nanobeams which makes a quasi-zero stiffness. Further calculations show that, when under a large deformation, the average stress in nanobeam is ultra-small (<151 MPa) and its phonon density of states are little changed. This work provides valuable insights on multifunction, such as both stable thermal and mechanical properties, of nano-designed metamaterials.

Thermal conductivity is a critical material property in numerous applications, such as those related to thermoelectric devices and heat dissipation. Effectively modulating thermal conductivity has become a great concern in the field of heat conduction. In this study, a quantum strategy is proposed to modulate thermal conductivity by exciting targeted phonons. The results show that the thermal conductivity of graphene can be tailored in the range of 1559 W/m-K (49%) to 4093 W/m-K (128%), compared with the intrinsic value of 3189 W/m-K. A similar trend is also observed for graphene nanoribbons. The results are obtained through both ab initio calculations and molecular dynamics simulations. This brand-new quantum strategy to modulate thermal conductivity paves a way for quantum heat conduction.

In a current-carrying single-molecular junction (SMJ), a hierarchy of hybrid energy transport processes takes place under a highly nonequilibrium situation, including energy transfer from electrons to molecular vibrations via electron-vibration interaction, energy redistribution within different vibrational modes via anharmonic coupling, and eventual energy transport to surrounding electrodes. A comprehensive understanding of such processes is a prerequisite for their potential applications as single-molecular devices. $Ab$ $initio$ current-induced molecular dynamics (MD) is an ideal approach to address this complicated problem. But the computational cost hinders its usage in systematic study of realistic SMJs. Here, we achieve orders of magnitude improvement in the speed of MD simulation by employing machine-learning potential with accuracy comparable to density functional theory. Using this approach, we show that SMJs with graphene electrodes generate order of magnitude less heating than those with gold electrodes. Our work illustrates the superior heat transport property of graphene as electrodes for SMJs, thanks to its better phonon spectral overlap with molecular vibrations., Comment: 8 pages with supplemental materials

We develop a unified theory for the second sound in two dimensional materials. Previously studied drifting and driftless second sound are two limiting cases of the theory, corresponding to the drift and diffusive part of the energy flux, respectively. We find that due to the presence of quadratic flexural phonons the drifting second sound does not exist in the thermodynamic limit, while the driftless mode is less affected. This is understood as a result of infinite effective inertia of flexual phonons, due to their constant density states and divergent Bose-Einstein distribution in the long wave length limit. Consequently, the group velocity of the drifting mode is smaller than that of the driftless mode. However, upon tensile strain, the velocity of drifting mode becomes larger. Both of them increase with tensile strain due to the linearization of the flexural phonon dispersion. Our results clarify several puzzles encountered previously and pave the way for exploring wave-like heat transport beyond hydrodynamic regime.

Thermal conductance of single molecular junctions at room temperature has been measured recently using picowatt-resolution scanning probes. However, fully understanding thermal transport in a much wider temperature range is needed for the exploration of energy transfer at single-molecular limit and the development of single-molecular devices. Here, employing a semiclassical Langevin molecular dynamics method, a comparative study is performed on the thermal transport of an alkane chain between Au and graphene electrodes, respectively. We illustrate the different roles of quantum statistics and anharmonic interaction in the two types of junctions. For a graphene junction, quantum statistics is essential at room temperature, while the anharmonic interaction is negligible. For a Au junction, it is the other way. Our study paves the way for theoretically understanding thermal transport of realistic single-molecular junctions in the full temperature range by including both quantum statistics and anharmonic interaction within one theoretical framework., Comment: 7 pages, 7 figures

Solar still is an eco-friendly and convenient desalination system that can provide fresh water for remote areas and emergencies. The energy efficiency and productivity of conventional solar still are unsatisfying and need improvement, which requires a deep understanding of the heat and mass transfer process in solar still. In this work, the effect of the inner air circulation on the system's heat and mass transfer performance and energy efficiency are studied theoretically and experimentally. The theoretical results reveal that a weak acceleration of the air circulation inside the SS will significantly increase its performance, due to the improved mass transfer process. By enhancing the inner air circulation, the evaporation and condensation in the solar still can reach up to the limit, and the theoretical energy efficiency reaches up to 87%, 91.5%, and 94.5%, for the input power density at 300 W/m2, 500 W/m2, and 700 W/m2, respectively. Besides, lower ambient temperature and higher ambient convective heat transfer coefficient will decrease the energy efficiency. Given the heat loss, the experimental energy efficiencies are only 3% to 6% lower than the theoretical results, which indicates that the great performance predicted by the theory can be realized in practical application. This work provides a new understanding and strategy for improving the performance of the solar still.

Epoxy resins are widely used polymer matrices for numerous applications. Despite substantial advances, the molecular-level knowledge-base required to exploit these materials to their full potential remains limited. A deeper comprehension of structure/property relationships in epoxy resins at the molecular level is critical to progressing these efforts. It can be laborious, if not impractical, to elucidate these relationships based on experiments alone. Here, molecular dynamics simulations are used to calculate and compare thermal conductivities and mechanical properties of an exemplar epoxy resin, Bisphenol F cross-linked with Diethyl Toluene Diamine, revealing these inter-relationships. Both elastic modulus and thermal transport of the epoxy resin show an increase with greater cross-linking. Specifically, decomposition of the thermal conductivity into different force contributions suggests that the bonded term contributes to an increase in the heat flux. These outcomes provide a foundation for designing and fabricating customized epoxy resins with desirable thermal and mechanical attributes.

As a critical way to modulate thermal transport in nanostructures, phonon resonance hybridization has become an issue of great concern in the field of phonon engineering. In this work, we optimized phonon transport across graphene nanoribbon and obtained minimized thermal conductance by means of designing pillared nanostructures based on resonance hybridization. Specifically, the optimization of thermal conductance was performed by the combination of atomic Green` s function and Bayesian optimization. Interestingly, it is found that thermal conductance decreases non-monotonically with the increasing of number for pillared structure, which is severed as resonator and blocks phonon transport. Further mode-analysis and atomic Green` s function calculations revealed that the anomalous tendency originates from decreased phonon transmission in a wide frequency range. Additionally, nonequilibrium molecular dynamics simulations are performed to verify the results with the consideration of high-order phonon scattering. This finding provides novel insights into the control of phonon transport in nanostructures.

Investigating thermal transport in van der Waals heterostructure is of scientific interest and practical importance for their applications in a broad range. In this work, thermal conductivity of one-dimensional heterostructure consisting of carbon and boron nitride nanotubes is systematically investigated via molecular dynamics simulations. Thermal conductivity is found to have strong dependences on temperature, length and diameter. In addition, the axial strain and intensity of van der Waals interaction are demonstrated to be able to modulate thermal conductivity up to about 43% and 37%, respectively. Moreover, the dependence of thermal conductivity on the chirality of componential nanotubes is studied. These results are explained based on lattice dynamics insights. This work not only provides feasible strategies to modulate thermal conductivity, but also enhances the understanding of the fundamental physics of phonon transport in one-dimensional heterostructure., Comment: 19 pages, 6 figures

Tuning interfacial thermal conductance has been a key task for the thermal management of nanoelectronic devices. Here, it is studied how the interfacial thermal conductance is great influenced by modulating the mass distribution of the interlayer of one-dimensional atomic chain. By nonequilibrium Green's function and machine learning algorithm, the maximum/minimum value of thermal conductance and its corresponding mass distribution are calculated. Interestingly, the mass distribution corresponding to the maximum thermal conductance is not a simple function, such as linear and exponential distribution predicted in previous works, it is similar to a sinusoidal curve around linear distribution for larger thickness interlayer. Further, the mechanism of the abnormal results is explained by analyzing the phonon transmission spectra and density of states. The work provides deep insight into optimizing and designing interfacial thermal conductance by modulating mass distribution of interlayer atoms.

Exploring the mechanism of interfacial thermal transport and reducing the interfacial thermal resistance is of great importance for thermal management and modulation. Herein, the interfacial thermal resistance between overlapped graphene nanoribbons is largely reduced by adding bonded carbon chains by performing molecular dynamics simulations. And the analytical model (cross-interface model, CIM) is utilized to analyze and explain the two-dimensional thermal transport mechanism at cross-interface. An order of magnitude reduction in interfacial thermal resistance is found as the graphene nanoribbons are bonded by just one carbon chain. Interestingly, the decreasing rate of interfacial thermal resistance slows down gradually with the increasing of the number of carbon chains, which can be explained by the proposed theoretical relationship based on CIM. Moreover, by the comparison of CIM and traditional simplified model, the accuracy of CIM is verified and demonstrated in overlapped graphene nanoribbons. This work provides a new way to improve the interfacial thermal transport and reveal the essential mechanism for low-dimensional materials applied in thermal management.

Solar still is a convenient off-grid device for desalination, which can provide fresh water for families, ships, islands and so on. The conventional inclined solar still (ISS) suffers from low efficiency and low productivity. To improve the performance of solar still, a flat solar still (FSS) is proposed, which has a working principle similar to the solar cell. The condensate water in FSS is collected by the capillary grid attached under the ultra-hydrophilic glass cover, instead of by gravity. Therefore, FSS avoids the inclined structure and is much more compact than ISS. The daily productivity of FSS reaches up to 4.3 kg/m2. Theoretical analysis shows that the enhanced mass transfer in FSS by the compact structure is an important factor for high performance. More interestingly, FSS can also be easily extended to more stage for latent heat recovery. The results show that the daily productivity of a double-stage FSS reaches up to 7 kg/m2, which is much higher than the conventional solar still. FSS paves a new way in designing and optimizing of solar still., Comment: 34

Phonon engineering focuses on heat transport modulation on atomic-scale. Different from reported methods, it is shown that electric field can also modulate heat transport in ferroelectric polymers, poly(vinylidene fluoride), by both simulation and measurement. Interestingly, thermal conductivities of poly(vinylidene fluoride) array can be enhanced by a factor of 3.25 along the polarization direction by simulation. The semi-crystalline poly(vinylidene fluoride) film can be also enhanced by a factor of 1.5 which is found by both simulation and measurement. The morphology and phonon property analysis reveal that the enhancement arises from the higher inter-chain lattice order, stronger inter-chain interaction, higher phonon group velocity and suppressed phonon scattering. This study offers a new modulation strategy with quick response and without fillers.

Presented is a new generation prediction model of a tubular solar still (TSS) productivity utilizing two machine learning (ML) techniques, namely:Random forest (RF) and Artificial neural network (ANN). Prediction models were conducted based on experimental data recorded under Egyptian climate. Meteorological and operational thermal parameters were utilized as input layers. Moreover, Bayesian optimization algorithm (BOA) was used to obtain the optimal performance of RF and ANN models. In addition, these models results were compared to those of a multilinear regression (MLR) model. As resulted, experimentally, the average value accumulated productivity was 4.3 L/(m2day). For models results, RF was less sensitive to hyper parameters than ANN as ANN performance could be significantly improved by BOA more than RF. In addition, RF achieved better prediction performance of TSS on the current dataset. The determination coefficients (R2) of RF and ANN were 0.9964 and 0.9977, respectively, which were much higher than MLR models, 0.9431. Based on the robustness performance and high accuracy, RF is recommended as a stable method for predicting the productivity of TSS.

Tailoring thermal properties with nanostructured materials can be of vital importance for many applications. Generally classical phonon size effects are employed to reduce the thermal conductivity, where strong phonon scattering by nanostructured interfaces or boundaries can dramatically supress the heat conduction. When these boundaries or interfaces are arranged in a periodic pattern, coherent phonons may have interference and modify the phonon dispersion, leading to dramatically reduced thermal conductivity. Such coherent phonon transport has been widely studied for superlattice films and recently emphasized for periodic nanoporous patterns. Although the wave effects have been proposed for reducing the thermal conductivity, more recent experimental evidence shows that such effects can only be critical at an ultralow temperature, i.e., around 10 K or below. At room temperature, the impacted phonons are mostly restricted to hypersonic modes that contribute little to the thermal conductivity. In this review, the theoretical and experimental studies of periodic porous structures are summarized and compared. The general applications of periodic nanostructured materials are further discussed.

The application of low-dimensional materials for heat dissipation requires a comprehensive understanding of the thermal transport at the cross interface, which widely exists in various composite materials and electronic devices. In this work, we proposed an analytical model, named as cross interface model (CIM), to accurately reveal the essential mechanism of the two-dimensional thermal transport at the cross interface. The applicability of CIM is validated through the comparison of the analytical results with molecular dynamics simulations for a typical cross interface of two overlapped boron nitride nanoribbons. Besides, it is figured out that the factor ({\eta}) has important influence on the thermal transport besides the thermal resistance inside and between the materials, which is found to be determined by two dimensionless parameters from its expression. Our investigations deepen the understanding of the thermal transport at the cross interface and also facilitate to guide the applications of low-dimensional materials in thermal management.

In last decades, solar stills, as one of the solar desalination technologies, have been well studied in terms of their productivity, efficiency and economics. Recently, to overcome the bottleneck of traditional solar still, improving solar still by optimizing the solar evaporation process based on micro/nanomaterials have been proposed as a promising strategy. In this review, the recent development for achieving high-performance of solar still and solar evaporation are discussed, including materials as well as system configurations. Meanwhile, machine learning was used to analyze the importance of different factors on solar evaporation, where thermal design was founded to be the most significant parameter that contributes in high-efficiency solar evaporation. Moreover, several important points for the further investigations of solar still and solar evaporation were also discussed, including the temperature of the air-water interface, salt rejecting and durability, the effect of solid-liquid interaction on water phase change., Comment: 53pages

Single layer molybdenum disulfide (SLMoS2), a semiconductor possesses intrinsic bandgap and high electron mobility, has attracted great attention due to its unique electronic, optical, mechanical and thermal properties. Although thermal conductivity of SLMoS2 has been widely investigated recently, less studies focus on molybdenum disulfide nanotube (MoS2NT). Here, the comprehensive temperature, size and strain effect on thermal conductivity of MoS2NT are investigated. A chirality-dependent strain effect is identified in thermal conductivity of zigzag nanotube, in which the phonon group velocity can be significantly reduced by strain. Besides, results show that thermal conductivity has a ~T-1 and a ~L\b{eta} relation with temperature from 200 to 400 K and length from 10 to 320 nm, respectively. This work not only provides feasible strategies to modulate the thermal conductivity of MoS2NT, but also offers useful insights into the fundamental mechanisms that govern the thermal conductivity, which can be used for the thermal management of low dimensional materials in optical, electronic and thermoelectrical devices. Introduction.

Understanding thermal transport through nanoscale van der Waals interfaces is vital for addressing thermal management challenges in nanoelectronic devices. In this work, the interfacial thermal conductance (GCA) between copper phthalocyanine (CuPc) nanoribbons is reported to be on the order of 10^5 Wm-2K-1 at 300 K, which is over two orders of magnitude lower than the value predicted by molecular dynamics (MD) simulations for a perfectly smooth interface between two parallelly aligned CuPc nanoribbons. Further MD simulations and contact mechanics analysis reveal that surface roughness can significantly reduce the adhesion energy and effective contact area between CuPc nanoribbons, and thus result in an ultralow GCA. In addition, the adhesion energy at the interface also depends on the stacking configuration of two CuPc nanoribbons, which may also contribute to the observed ultralow GCA., Comment: Main manuscript: 21 pages, 3 figures. Supplementary material: 13 pages, 5 figures and 1 table

The convection plays a very important role in heat transfer when MEMS work under air environment. However, traditional measurements of convection heat transfer coefficient require the knowledge of thermal conductivity, which makes measurements complex. In this work, a modified steady state "hot wire" (MSSHW) method is proposed, which can measure the heat transfer coefficient of microwires' convection without the knowledge of thermal conductivity. To verify MSSHW method, the convection heat transfer coefficient of platinum microwires was measured in the atmosphere, whose value is in good agreement with values by both traditional measurement methods and empirical equations. Then, the convection heat transfer coefficient of microwires with different materials and diameters were measured by MSSHW. It is found that the convection heat transfer coefficient of microwire is not sensitive on materials, while it increases from 86 W/(m$^2$K) to 427 W/(m$^2$K) with the diameter of microwires decreasing from 120 ${\mu}$m to 20 ${\mu}$m. Without knowing thermal conductivity of microwires, the MSSHW method provides a more convenient way to measure the convective effect., Comment: 16 pages, 6 figures

Hotspot is a ubiquitous phenomenon in microdevices/chips. In homogeneous nanoscale graphene disk with a hotspot, a graded thermal conductivity is observed previously even when the system size is fixed. However, the underlying physical mechanism is not clear. In this work, the hotspots in homogeneous 2D disk/3D ball and graphene disk are studied based on phonon Boltzmann transport equation. The mechanisms of phonon scattering are analyzed. It is found that for a system with fixed size, the graded thermal conductivity is predictable as long as there is not sufficient phonon scattering, which is independent on material properties, dimensions or system size. This work may shed light on both theoretical and experimental studies on heat dissipation of microelectronics., Comment: 5 pages, 3 figures, 48 references

Silicon Carbide (SiC) is a typical material for third-generation semiconductor. The thermal boundary resistance (TBR) of 4H-SiC/SiO2 interface, was investigated by both experimental measurements and theoretical calculations. The structure of 4H-SiC/SiO2 was characterized by using transmission electron microscopy and X-ray diffraction. The TBR is measured as 8.11*10-8 m2K/W by 3-omega method. Furthermore, the diffuse mismatch model was employed to predict the TBR of different interfaces which is in good agreement with measurements. Heat transport behavior based on phonon scattering perspective was also discussed to understand the variations of TBR across different interfaces. Besides, the intrinsic thermal conductivity of SiO2 thin films (200~1,500 nm in thickness) on 4H-SiC substrates was measured by 3 omega procedure, as 1.42 W/mK at room temperature. It is believed the presented results could provide useful insights on the thermal management and heat dissipation for SiC devices.

There has been an increasing demand for materials with special thermal properties, whereas experimental discovery is high-cost and time-consuming. The emerging discipline `Materials Informatics' is an effective approach that can accelerate materials development by combining material science and big data technique. Recently materials informatics has been applied to the design of novel materials such as thermal interface materials for heat-dissipation, and thermoelectric materials for power generation. This mini-review summarized the research progress on the applications of materials informatics for the thermal transport properties prediction and discovery of materials with special thermal properties, including optimal thermal conductivity, interfacial thermal conductance and thermoelectricity efficiency. In addition, some perspectives are given for the outlook of materials informatics in the field of thermal transport.

Solar steam generation technology can utilize abundant and renewable solar energy for many applications. In this work, we proposed a solar steam generator using wick material with wrapped graphene nanoparticles, and the energy efficiency can reaches up to 80%. Instead of traditional smearing method, the chemical wrapping method was used to better adhere the graphene nanoparticles on the wick materials. Through the SEM morphological results, the graphene nanoparticles are shown to be evenly wrapped across the fibres of the wick material, which have better dispersity and stability. The evaporation rate, instantaneous energy efficiency and the absorptivity of three wick materials with/without nanoparticles under different conditions were compared and analyzed. Among the three different wick materials, the flannel cloth with dense fine hairs can provide three-dimensional contact area for wrapping graphene nanoparticles and thus contribute to better evaporation. Additionally, the influence of two different reduction methods and different concentrations of graphene oxide solution on the energy efficiency was also discussed. Our work offers a simple and effective way of using nanotechnology in practical application for solar steam generation.

Solar evaporation is important for many applications such as desalination, power generation and industrial drying. Recently, some studies on evaporation reported obtaining high energy efficiency and evaporation rate, which are based on floating evaporation setup (FES) with nanomaterials. Here, we proposed a new cheap and simple setup, named as airing evaporation setup (AES). It shows that the energy efficiency of AES reaches up to 87 % under 1 kW/m2 of solar irradiation, which is 14% higher than that of FES. Meanwhile, the total evaporation rate of AES is about 20% higher than that of FES. The theoretical analysis reveals that the main reason for a better performance of AES is the increasing evaporation area. More interesting, AES could be used for designing portable systems due to its simplicity and flexibility. Furthermore, we show that AES and the corresponding wick material can be used in solar desalination, textile quick-drying and warm-keeping., Comment: 21pages, 5 figures

Thermal transport in nanoribbon based nanostructures is critical to advancing its applications. Wave effects of phonons can give rise to controllability of heat conduction in nanostructures beyond that by particle scattering. In this paper, by introducing pillars to form structural resonance, we systematically studied the thermal conductivity of graphene nanoribbon based phononic metamaterials (GNPM) through non-equilibrium molecular dynamical simulation. Interestingly, it is found that the thermal conductivity of GNPM is counter intuitively enhanced by isotope engineering, which is strong contrast to the common notion that isotope engineering reduces thermal conductivity.Further mode analysis and atomic Green function calculation reveal that the unexpected increasing in thermal conductivity originates from the breaking of the resonant hybridization wave effect between the resonant modes and the propagating modes induced by isotope engineering. Besides, factors including the system width and pillar height can also efficiently tune the thermal conductivity of GNPM. This abnormal mechanism provides a new dimension to manipulate phonon transport in nanoribbon based nanostructures through wave effect.