Logarithmic divergent thermal conductivity in two-dimensional nonlinear lattices

Quasiballistic thermal transport in submicron-scale graphene nanoribbons at room-temperature

Phonon transport in two-dimensional materials has been the subject of intensive studies both theoretically and experimentally. Recently observed unique phenomena such as Poiseuille flow at low temperature in graphene nanoribbons (GNRs) initiated strong interest in similar effects at higher temperatures. Here, we carry out massive molecular dynamics simulations to examine thermal transport in GNRs at room temperature (RT) and demonstrate that non-diffusive behaviors including Poiseuille-like local thermal conductivity and second sound are obtained, indicating quasiballistic thermal transport. For narrow GNRs, a Poiseuille-like thermal conductivity profile develops across the nanoribbon width, and wider GNRs exhibit a mixed nature of phonon transport in that diffusive transport is dominant in the middle region whereas non-uniform behavior is observed near lateral GNR boundaries. In addition, transient heating simulations reveal that the driftless second sound can propagate through GNRs regardless of the GNR width. By decomposing the atomic motion into out-of-plane and in-plane modes, it is further shown that the observed quasiballistic thermal transport is primarily contributed by the out-of-plane motion of C atoms in GNRs.

Heat transport in an angular-momentum-conserving lattice

It is expected that the energy-diffusion propagator in a one-dimensional nonlinear lattice with three conserved quantities: energy, momentum, and stretch, consists of a central heat mode and two sound modes. The heat mode follows a Lévy distribution. Consequently, the heat diffusion is super, i.e., the second moment of the diffusion propagator diverges as t^{β} with β>1; and the heat conduction is anomalous, i.e., the heat conductivity is size dependent and diverges with size N by N^{α}, with α>0. In this paper, we study a one-dimensional lattice with two-dimensional transverse motions, in which the total angular momentum also conserves. More importantly, the diffusion of this conserved quantity is ballistic. Surprisingly, the above pictures and the values of the mentioned power exponents keep unchanged. The universality of the scalings is then further extended. On the other hand, the detailed strengths of heat transports are largely enhanced. Such a counterintuitive finding can be explained by the change of the phonon mean-free path of the lattices.

Energy diffusion in two-dimensional momentum-conserving nonlinear lattices: Lévy walk and renormalized phonon

The energy diffusion process in a few two-dimensional Fermi-Pasta-Ulam-type lattices is numerically simulated via the equilibrium local energy spatiotemporal correlation. Just as the nonlinear fluctuating hydrodynamic theory suggested, the diffusion propagator consists of a bell-shaped central heat mode and a sound mode extending with a constant speed. The profiles of the heat and sound modes satisfy the scaling properties from a random-walk-with-velocity-fluctuation process very well. An effective phonon approach is proposed, which expects the frequencies of renormalized phonons as well as the sound speed with quite good accuracy. Since many existing analytical and numerical studies indicate that heat conduction in such two-dimensional momentum-conserving lattices is divergent and the thermal conductivity κ increases logarithmically with lattice length, it is expected that the mean-square displacement of energy diffusion grows as tlnt. Discrepancies, however, are noticeably observed.

Temperature Dependence of Thermal Conductivity of Giant-Scale Supported Monolayer Graphene

Past work has focused on the thermal properties of microscale/nanoscale suspended/supported graphene. However, for the thermal design of graphene-based devices, the thermal properties of giant-scale (~mm) graphene, which reflects the effect of grains, must also be investigated and are critical. In this work, the thermal conductivity variation with temperature of giant-scale chemical vapor decomposition (CVD) graphene supported by poly(methyl methacrylate) (PMMA) is characterized using the differential transient electrothermal technique (diff-TET). Compared to the commonly used optothermal Raman technique, diff-TET employs joule heating as the heating source, a situation under which the temperature difference between optical phonons and acoustic phonons is eased. The thermal conductivity of single-layer graphene (SLG) supported by PMMA was measured as 743 ± 167 W/(m·K) and 287 ± 63 W/(m·K) at 296 K and 125 K, respectively. As temperature decreased from 296 K to 275 K, the thermal conductivity of graphene was decreased by 36.5%, which can be partly explained by compressive strain buildup in graphene due to the thermal expansion mismatch.

Heat current flows across an interface in two-dimensional lattices

Heat current J that flows through a few typical two-dimensional nonlinear lattices is systematically studied. Each lattice consists of two identical segments that are coupled by an interface with strength k_{int}. It is found that the two-universality-class scenario that is revealed in one-dimensional systems is still valid in the two-dimensional systems. Namely, J may follow k_{int} in two entirely different ways, depending on whether or not the interface potential energy decays to zero. We also study the dependence of J on lattice width N_{Y} and transverse interaction strength k_{Y}. Universal power-law decay or divergence is observed. Finally, we check the equipartition theorem in the systems since it is the basis of all our theoretical analyses. Surprisingly, it holds perfectly even at the interface where there is a finite temperature jump, which makes the system far from equilibrium. However, the equipartition of potential energy, which is observed in one-dimensional systems, is no longer satisfied due to the interaction between different dimensions.

Magnetic field-dependent resistance crossover and anomalous magnetoresistance in topological insulator Bi2Te3

We report a metal-insulator like transition in single-crystalline 3D topological insulator Bi2Te3at a temperature of 230 K in the presence of an external magnetic field applied normal to the surface. This transition becomes more prominent at larger magnetic field strength with the residual resistance value increasing linearly with the magnetic field. At low temperature, the magnetic field dependence of the magnetoresistance shows a transition from logarithmic to linear behavior and the onset magnetic field value for this transition decreases with increasing temperature. The logarithmic magnetoresistance indicates the weak anti-localization of the surface Dirac electrons while the high temperature behavior originates from the bulk carriers due to intrinsic impurities. At even higher temperatures beyond ∼230 K, a completely classical Lorentz model type quadratic behavior of the magnetoresistance is observed. We also show that the experimentally observed anomalies at ∼230 K in the magneto-transport properties do not originate from any stacking fault in Bi2Te3.

Thermal transport of carbon nanomaterials

The diversity of thermal transport properties in carbon nanomaterials enables them to be used in different thermal fields such as heat dissipation, thermal management, and thermoelectric conversion. In the past two decades, much effort has been devoted to study the thermal conductivities of different carbon nanomaterials. In this review, different theoretical methods and experimental techniques for investigating thermal transport in nanosystems are first summarized. Then, the thermal transport properties of various pure carbon nanomaterials including 1D carbon nanotubes, 2D graphene, 3D carbon foam, are reviewed in details and the associated underlying physical mechanisms are presented. Meanwhile, we discuss several important influences on the thermal conductivities of carbon nanomaterials, including size, structural defects, chemisorption and strain. Moreover, we introduce different nanostructuring pathways to manipulate the thermal conductivities of carbon-based nanocomposites and focus on the wave nature of phonons for controlling thermal transport. At last, we briefly review the potential applications of carbon nanomaterials in the fields of thermal devices and thermoelectric conversion.

Anomalous energy diffusion in two-dimensional nonlinear lattices

Heat transport in one-dimensional (1D) momentum-conserved lattices is generally assumed to be anomalous, thus yielding a power-law divergence of thermal conductivity with system length. However, whether heat transport in a two-dimensional (2D) system is anomalous or not is still up for debate because of the difficulties involved in experimental measurements or due to the insufficiently large simulation cell size. Here we simulate energy and momentum diffusion in the 2D nonlinear lattices using the method of fluctuation correlation functions. Our simulations confirm that energy diffusion in the 2D momentum-conserved lattices is anomalous and can be well described by the Lévy-stable distribution. As is expected, we verify that 2D nonlinear lattices with on-site potentials exhibit normal energy diffusion, independent of the dimension. Contrary to the hypothesis of a 1D system, we further clarify that anomalous heat transport in the 2D momentum-conserved system cannot be corroborated by the momentum superdiffusion any longer. Our findings offer some valuable insights into mechanisms of thermal transport in 2D system.

Energy and spin diffusion in the one-dimensional classical Heisenberg spin chain at finite and infinite temperatures

The energy and spin diffusion behaviors in the one-dimensional classical Heisenberg spin chain have been systematically investigated using the equilibrium diffusion method. The spatiotemporal autocorrelation functions for energy and spin are calculated at finite and infinite temperatures. As conserved quantities, the spreading of excess energy and spin can be used to determine their actual diffusion behaviors. At low temperatures, the energy diffusion shows almost ballistic behavior, and spin shows superdiffusion behavior for finite chain size. For energy diffusion, normal diffusion behavior can be obtained when the temperature is higher than 0.75. For spin diffusion, normal diffusion behavior is observed at infinite temperature.

Disorder limits the coherent phonon transport in two-dimensional phononic crystal structures

Recently, increasing efforts are being made to control thermal transport via coherent phonons in periodic phononic structures; however, the direct observation of coherent phonon transport is experimentally very difficult at ambient temperature, and the importance of coherent phonons to the total thermal conductivity has not been critically assessed to date. In this study, using the non-equilibrium molecular dynamics simulations, we studied coherent phonon transport in a C3N phononic crystal (CNPnC) structure at room temperature by changing the porosity. When the holes were randomly distributed to construct the disordered C3N (D-C3N) structure, the localization of the coherent phonons was revealed by the phonon transmission coefficient, phonon wave packet simulation, phonon participation ratio and spatial energy density, which led to a significant reduction in the thermal conductivity. Finally, the effects of the length, temperature and strain on the thermal conductivity of CNPnC and D-C3N have also been discussed. Our study provides a solid understanding of the coherent phonon transport behavior, which will be beneficial for phononic-related control based on coherent phonons.

Influence of the degree of a complex network on heat conduction

Devices made of nanotubes and nanowires networks are of great interest for applications and have caught increasing attention in recent years. In this work, we study heat conduction in a network model with nodes being atoms and links being one-dimensional chains of atoms. It is found that heat conduction in the complex network is fundamentally different from that of regular lattices. It depends very sensitively on the average degrees of complex networks and the degrees of nodes that are attached to the two heat baths. For example, when the two heat source nodes have the same degree k_{0}, the heat flux reaches a maximum at an optimized value of k_{0} and decreases with the increase of the average degree 〈k〉. In other words, the source nodes with optimal degree k_{0} and the sparse network are more favorable to heat flux. Thermal rectification effect is found when the two heat source nodes have different degrees or the network model has multiple heat source nodes. Theoretical analysis is provided to explain the numerical results.

Randomness-Induced Phonon Localization in Graphene Heat Conduction

Through nonequilibrium molecular dynamics simulations, we report the direct numerical evidence of the coherent phonons participating in thermal transport at room temperature in graphene phononic crystal (GPnC) structure and evaluate their contribution to thermal conductivity based on the two-phonon model. With decreasing period length in GPnC, the transition from the incoherent to coherent phonon transport is clearly observed. When a random perturbation to the positions of holes is introduced in a graphene sheet, the phonon wave-packet simulation reveals the presence of notable localization of coherent phonons, leading to the significant reduction of thermal conductivity and suppressed length dependence. Finally, the effects of period length and temperature on the coherent phonon contribution to thermal conductivity are also discussed. Our work establishes a deep understanding of the coherent phonons transport behavior in periodic phononic structures, which provides effective guidance for engineering thermal transport based on a new path via phonon localization.

Renormalized vibrations and normal energy transport in 1d FPU-like discrete nonlinear Schrödinger equations

For one-dimensional (1d) nonlinear atomic lattices, the models with on-site nonlinearities such as the Frenkel-Kontorova (FK) and ϕ⁴ lattices have normal energy transport while the models with inter-site nonlinearities such as the Fermi-Pasta-Ulam-β (FPU-β) lattice exhibit anomalous energy transport. The 1d Discrete Nonlinear Schrödinger (DNLS) equations with on-site nonlinearities has been previously studied and normal energy transport has also been found. Here, we investigate the energy transport of 1d FPU-like DNLS equations with inter-site nonlinearities. Extended from the FPU-β lattice, the renormalized vibration theory is developed for the FPU-like DNLS models and the predicted renormalized vibrations are verified by direct numerical simulations same as the FPU-β lattice. However, the energy diffusion processes are explored and normal energy transport is observed for the 1d FPU-like DNLS models, which is different from their atomic lattice counterpart of FPU-β lattice. The reason might be that, unlike nonlinear atomic lattices where models with on-site nonlinearities have one less conserved quantities than the models with inter-site nonlinearities, the DNLS models with on-site or inter-site nonlinearities have the same number of conserved quantities as the result of gauge transformation.

Multiparticle collision simulations of two-dimensional one-component plasmas: Anomalous transport and dimensional crossovers

By means of hybrid multiparticle collsion-particle-in-cell (MPC-PIC) simulations we study the dynamical scaling of energy and density correlations at equilibrium in moderately coupled two-dimensional (2D) and quasi-one-dimensional (1D) plasmas. We find that the predictions of nonlinear fluctuating hydrodynamics for the structure factors of density and energy fluctuations in 1D systems with three global conservation laws hold true also for 2D systems that are more extended along one of the two spatial dimensions. Moreover, from the analysis of the equilibrium energy correlators and density structure factors of both 1D and 2D neutral plasmas, we find that neglecting the contribution of the fluctuations of the vanishing self-consistent electrostatic fields overestimates the interval of frequencies over which the anomalous transport is observed. Such violations of the expected scaling in the currents correlation are found in different regimes, hindering the observation of the asymptotic scaling predicted by the theory.

Role of the range of the interactions in thermal conduction

We investigate thermal transport along a one-dimensional lattice of classical inertial rotators, with attractive couplings that decrease with distance as r^{-α} (α≥0), subject at its ends to Brownian heat reservoirs at different temperatures with average value T. By means of numerical integration of the equations of motion, we show the effects of the range of the interactions in the temperature profile and energy transport and determine the domain of validity of Fourier's law in this context. We find that Fourier's law, as signaled by a finite κ in the thermodynamic limit, holds only for sufficiently short-range interactions, with α>α_{c}(T). For α<α_{c}(T), a kind of insulator behavior emerges at any T.

Stretch diffusion and heat conduction in one-dimensional nonlinear lattices

For heat conduction in one-dimensional (1D) nonlinear Hamiltonian lattices, it has been known that conserved quantities play an important role in determining the actual heat conduction behavior. In closed or microcanonical Hamiltonian systems, the total energy and stretch are always conserved. Depending on the existence of external on-site potential, the total momentum can be conserved or not. All the momentum-conserving lattices have anomalous heat conduction except the 1D coupled rotator lattice. It was recently claimed that "whenever stretch (momentum) is not conserved in a 1D model, the momentum (stretch) and energy fields exhibit normal diffusion." The stretch in a coupled rotator lattice was also argued to be nonconserved due to the requirement of a finite partition function, which enables the coupled rotator lattice to fulfill this claim. In this work, we will systematically investigate stretch diffusion and heat conduction in terms of energy diffusion for typical 1D nonlinear lattices. Contrary to what was claimed, no clear connection between conserved quantities and heat conduction can be established. The actual situation might be more complicated than what was proposed.

Heat conduction and energy diffusion in momentum-conserving one-dimensional full-lattice ding-a-ling model

The ding-a-ling model is a kind of half lattice and half hard-point-gas (HPG) model. The original ding-a-ling model proposed by Casati et al. does not conserve total momentum and has been found to exhibit normal heat conduction behavior. Recently, a modified ding-a-ling model which conserves total momentum has been studied and normal heat conduction has also been claimed. In this work, we propose a full-lattice ding-a-ling model without hard point collisions where total momentum is also conserved. We investigate the heat conduction and energy diffusion of this full-lattice ding-a-ling model with three different nonlinear inter-particle potential forms. For symmetrical potential lattices, the thermal conductivities diverges with lattice length and their energy diffusions are superdiffusive signaturing anomalous heat conduction. For asymmetrical potential lattices, although the thermal conductivity seems to converge as the length increases, the energy diffusion is definitely deviating from normal diffusion behavior indicating anomalous heat conduction as well. No normal heat conduction behavior can be found for the full-lattice ding-a-ling model.

Heat-current correlation loss induced by finite-size effects in a one-dimensional nonlinear lattice

The Green-Kubo formula provides a mathematical expression for heat conductivity in terms of integrals of the heat-current correlation function, which should be calculated in the thermodynamic limit. In finite systems this function generally decreases, i.e., it decays faster than it does in infinite systems. We compared the values of the correlation function in a one-dimensional purely quartic lattice with various lengths, and found that this loss is much smaller than is conventionally estimated. By studying the heat diffusion process in this lattice, we found that, in contrast to the conventional belief, the collisions between sound modes do not noticeably affect the current correlation function. Therefore, its loss being surprisingly small can be well understood. This finding allows one to calculate the heat conductivity in a very large system with desirable accuracy by performing simulations in a system with much smaller size, and thus greatly broadens the application of the Green-Kubo method.

Temperature-dependent thermal conductivities of one-dimensional nonlinear Klein-Gordon lattices with a soft on-site potential

The temperature-dependent thermal conductivities of one-dimensional nonlinear Klein-Gordon lattices with soft on-site potential (soft-KG) are investigated systematically. Similarly to the previously studied hard-KG lattices, the existence of renormalized phonons is also confirmed in soft-KG lattices. In particular, the temperature dependence of the renormalized phonon frequency predicted by a classical field theory is verified by detailed numerical simulations. However, the thermal conductivities of soft-KG lattices exhibit the opposite trend in temperature dependence in comparison with those of hard-KG lattices. The interesting thing is that the temperature-dependent thermal conductivities of both soft- and hard-KG lattices can be interpreted in the same framework of effective phonon theory. According to the effective phonon theory, the exponents of the power-law dependence of the thermal conductivities as a function of temperature are only determined by the exponents of the soft or hard on-site potentials. These theoretical predictions are consistently verified very well by extensive numerical simulations.

Heat conduction in nanoscale materials: a statistical-mechanics derivation of the local heat flux

We derive a coarse-grained model for heat conduction in nanoscale mechanical systems. Starting with an all-atom description, this approach yields a reduced model, in the form of conservation laws of momentum and energy. The model closure is accomplished by introducing a quasilocal thermodynamic equilibrium, followed by a linear response approximation. Of particular interest is the constitutive relation for the heat flux, which is expressed nonlocally in terms of the spatial and temporal variation of the temperature. Nanowires made of copper and silicon are presented as examples.

Length-dependent thermal conductivity in suspended single-layer graphene

Graphene exhibits extraordinary electronic and mechanical properties, and extremely high thermal conductivity. Being a very stable atomically thick membrane that can be suspended between two leads, graphene provides a perfect test platform for studying thermal conductivity in two-dimensional systems, which is of primary importance for phonon transport in low-dimensional materials. Here we report experimental measurements and non-equilibrium molecular dynamics simulations of thermal conduction in suspended single-layer graphene as a function of both temperature and sample length. Interestingly and in contrast to bulk materials, at 300 K, thermal conductivity keeps increasing and remains logarithmically divergent with sample length even for sample lengths much larger than the average phonon mean free path. This result is a consequence of the two-dimensional nature of phonons in graphene, and provides fundamental understanding of thermal transport in two-dimensional materials.

Thermal transport in graphene oxide--from ballistic extreme to amorphous limit

Graphene oxide is being used in energy, optical, electronic and sensor devices due to its unique properties. However, unlike its counterpart - graphene - the thermal transport properties of graphene oxide remain unknown. In this work, we use large-scale molecular dynamics simulations with reactive potentials to systematically study the role of oxygen adatoms on the thermal transport in graphene oxide. For pristine graphene, highly ballistic thermal transport is observed. As the oxygen coverage increases, the thermal conductivity is significantly reduced. An oxygen coverage of 5% can reduce the graphene thermal conductivity by ~90% and a coverage of 20% lower it to ~8.8 W/mK. This value is even lower than the calculated amorphous limit (~11.6 W/mK for graphene), which is usually regarded as the minimal possible thermal conductivity of a solid. Analyses show that the large reduction in thermal conductivity is due to the significantly enhanced phonon scattering induced by the oxygen defects which introduce dramatic structural deformations. These results provide important insight to the thermal transport physics in graphene oxide and offer valuable information for the design of graphene oxide-based materials and devices.

Validity of Fourier's law in one-dimensional momentum-conserving lattices with asymmetric interparticle interactions

We have numerically studied heat conduction in a few one-dimensional momentum-conserving lattices with asymmetric interparticle interactions by the nonequilibrium heat bath method, the equilibrium Green-Kubo method, and the heat current power spectra analysis. Very strong finite-size effects are clearly observed. Such effects make the heat conduction obey a Fourier-like law in a wide range of lattice lengths. However, in yet longer lattice lengths, the heat conductivity regains its power-law divergence. Therefore, the power-law divergence of the heat conductivity in the thermodynamic limit is verified, as is expected by many existing theories.