Catalog footprint

What is connected

110works
28topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

110 published item(s)

preprint2022arXiv

A unified theory of second sound in two dimensional 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.

preprint2021arXiv

Acoustically Manipulating Internal Structure of Disk-in-Sphere Endoskeletal Droplets

Manipulation of micro/nano particles has been well studied and demonstrated by optical, electromagnetic, and acoustic approaches, or their combinations. Manipulation of internal structure of droplet/particle is rarely explored and remains challenging due to its complicated nature. Here we demonstrated the manipulation of internal structure of disk-in-sphere endoskeletal droplets using acoustic wave for the first time. We developed a model to investigate the physical mechanisms behind this novel phenomenon. Theoretical analysis of the acoustic interactions indicated that these assembly dynamics arise from a balance of the primary and secondary radiation forces. Additionally, the disk orientation was found to change with acoustic driving frequency, which allowed on-demand, reversible adjusting disk orientations with respect to the substrate. This novel dynamic behavior leads to unique reversible arrangements of the endoskeletal droplets and their internal architecture, which may provide a new avenue for directed assembly of novel hierarchical colloidal architectures and intracellular organelles or intra-organoid structures.

preprint2021arXiv

Reciprocity of thermal diffusion in time-modulated systems

The reciprocity principle governs the symmetry in transmission of electromagnetic and acoustic waves, as well as the diffusion of heat between two points in space, with important consequences for thermal management and energy harvesting. There has been significant recent interest in materials with time-modulated properties, which have been shown to efficiently break reciprocity for light, sound, and even charge diffusion. Quite surprisingly, here we show that, from a practical point of view, time modulation cannot generally be used to break reciprocity for thermal diffusion. We establish a theoretical framework to accurately describe the behavior of diffusive processes under time modulation, and prove that thermal reciprocity in dynamic materials is generally preserved by the continuity equation, unless some external bias or special material is considered. We then experimentally demonstrate reciprocal heat transfer in a time-modulated device. Our findings correct previous misconceptions regarding reciprocity breaking for thermal diffusion, revealing the generality of symmetry constraints in heat transfer, and clarifying its differences from other transport processes in what concerns the principles of reciprocity and microscopic reversibility.

preprint2020arXiv

Effect of Interfacial Thermal Resistance in Thermal Cloak

When heat transfers through interface between two different materials, it will encounter an interfacial thermal resistance (ITR) that makes the temperature discontinuous. This effect has been totally neglected so far in the research of thermal cloak, in particular when the thermal cloak is built with multilayer structures. In this paper, we investigate the effect of ITR on the performance of the thermal cloak by using both analytical and numerical method. Our results show that the existence of ITR will distort the external field, thus destroy the cloak. Moreover, we found that the effect of ITR can be quantified by a parameter called characteristic length.

preprint2020arXiv

Quantum energy transfer between nonlinearly-coupled bosonic bath and a fermionic chain: an exactly solvable model

The evolution of a quantum system towards thermal equilibrium is usually studied by approximate methods, which have their limits of validity and should be checked against analytically solvable models. In this paper, we propose an analytically solvable model to investigate the energy transfer between a bosonic bath and a fermionic chain which are nonlinearly-coupled to each other. The bosonic bath consists of an infinite collection of non-interacting bosonic modes, while the fermionic chain is represented by a chain of interacting fermions with nearest-neighbor interactions. We compare behaviors of the temperature-dependent energy current $J_{T}$ and temperature-independent energy current $J_{TI}$ for different bath configurations. With respect to the bath spectrum, $J_{T}$ decays exponentially for Lorentz-Drude type bath, which is the same as the conventional approximations. On the other hand, the decay rate is $1/t^{3}$ for Ohmic type and $1/t$ for white noise, which doesn't have conventional counterparts. For the temperature-independent current $J_{TI}$, the decay rate is divergent for the Lorentz-Drude type bath, $1/t^{4}$ for the Ohmic bath, and $1/t$ for the white noise. When further considering the dynamics of the fermionic chain, the current will be modulated based on the envelope from the bath. As an example, for a bosonic bath with Ohmic spectrum, when the fermionic chain is uniformly-coupled, we have $J_{T}\propto1/t^{6}$ and $J_{TI}\propto1/t^{3}$. Remarkably, for perfect state transfer (PST) couplings, there always exists an oscillating quantum energy current $J_{TI}$. Moreover, it is interesting that $J_{T}$ is proportional to $(N-1)^{1/2}$ at certain times for PST couplings under Lorentz-Drude or Ohmic bath.

preprint2019arXiv

Phonon renormalization induced by electric field in ferroelectric P(VDF-TrFE) nanofibers

We report phonon renormalization induced by an external electric field E in ferroelectric poly(vinylidene fluoride-trifluoroethylene) [P(VDF-TrFE)] nanofibers through measuring the E-dependent thermal conductivity. Our experimental results are in excellent agreement with the theoretical ones derived from the lattice dynamics. The renormalization is attributed to the anharmonicity that modifies the phonon spectrum when the atoms are pulled away from their equilibrium positions by the electric field. Our finding provides an efficient way to manipulate the thermal conductivity by tuning external fields in ferroelectric materials.

preprint2018arXiv

Dimensional crossover of heat conduction in amorphous Polyimide nanofibers

The mechanism of thermal conductivity in amorphous polymers, especially polymer fibers, is unclear in comparison with that in inorganic materials. Here, we report the observation of across over of heat conduction behavior from three dimensions (3D) to quasi-one dimension (1D) in Polyimide(PI) nanofibers at a given temperature. A theoretical model based on the random walk theory has been proposed to quantitatively describe the interplay between the inter-chain hopping and the intra-chain hopping in nanofibers. This model explains well the diameter dependence of thermal conductivity and also speculates the upper limit of thermal conductivity of amorphous polymers in the quasi-1D limit.

preprint2016arXiv

A series circuit of thermal rectifiers: an effective way to enhance rectification ratio

The low rectification ratio limits the application of thermal rectifiers. In this paper, it is found that a series circuit of thermal rectifiers in asymmetric graphene/graphene phononic crystal (GPnC) structures can improve rectification significantly. The simulations by using non-equilibrium molecular dynamics (NEMD) are performed. Moreover, the size effect on thermal rectification is also studied and discussed.

preprint2016arXiv

Cargo Towing by Artificial Swimmers

An active swimmer can tow a passive cargo by binding it to form a self-propelling dimer. The orientation of the cargo relative to the axis of the active dimer's head is determined by the hydrodynamic interactions associated with the propulsion mechanism of the latter. We show how the tower-cargo angular configuration greatly influences the dimer's diffusivity and, therefore, the efficiency of the active swimmer as a micro-towing motor.

preprint2016arXiv

Interfacial thermal conductance across metal-insulator/semiconductor interfaces due to surface states

We point out that the effective channel for the interfacial thermal conductance, the inverse of Kapitza resistance, of metal-insulator/semiconductor interfaces is governed by the electron-phonon interaction mediated by the surface states allowed in a thin region near the interface. Our detailed calculations demonstrate that the interfacial thermal conductance across Pb/Pt/Al/Au-diamond interfaces are only slightly different among these metals, and reproduce well the experimental results of the interfacial thermal conductance across metal-diamond interfaces observed by Stoner et al. [Phys. Rev. Lett. 68, 1563 (1992)] and most recently by Hohensee et al. [Nature Commun. 6, 6578 (2015)].

preprint2016arXiv

Layer thickness-dependent phonon properties and thermal conductivity of MoS2

For conventional materials, the thermal conductivity of thin film is usually suppressed when the thickness decreases due to phonon-boundary scattering. However, this is not necessarily true for the van der Waals solids if the thickness is reduced to only a few layers. In this letter, the layer thickness-dependent phonon properties and thermal conductivity in the few-layer MoS2 are studied using the first-principles-based Peierls-Boltzmann transport equation approach. The basal-plane thermal conductivity is found to monotonically reduce from 138 W/mK to 98 W/mK for naturally occurring MoS2, and from 155 W/mK to 115 W/mK for isotopically pure MoS2, when its thickness increases from one layer to three layers. The thermal conductivity of tri-layer MoS2 approaches to that of bulk MoS2. Both the change of phonon dispersion and the thickness-induced anharmonicity are important for explaining such a thermal conductivity reduction. The increased anharmonicity in bi-layer MoS2 results in stronger phonon scattering for ZAi modes, which is linked to the breakdown of the symmetry in single-layer MoS2.

preprint2016arXiv

Manipulate Temperature Dependence of Thermal Conductivity of Graphene Phononic Crystal

By using non-equilibrium molecular dynamics simulations(NEMD), the modulation on temperature dependence of thermal conductivity of graphene phononic crystals (GPnCs) are investigated. It is found that the temperature dependence of thermal conductivity of GPnCs follows $T^{-α}$ behavior. The power exponents ($α$) can be efficiently tuned by changing the characteristic size of GPnCs. The phonon participation ratio spectra and dispersion relation reveal that the long-range phonon modes are more affected in GPnCs with larger size of holes ($L_0$). Our results suggest that constructing GPnCs is an effective method to manipulate the temperature dependence of thermal conductivity of graphene, which would be beneficial for developing GPnCs-based thermal management and signal processing devices.

preprint2016arXiv

Quantum thermal transport through anharmonic systems: A self-consistent approach

We propose a feasible and effective approach to study quantum thermal transport through anharmonic systems. The main idea is to obtain an {\it effective} harmonic Hamiltonian for the anharmonic system by applying the self-consistent phonon theory. Using the effective harmonic Hamiltonian we study thermal transport within the framework of nonequilibrium Green's function method using the celebrated Caroli formula. We corroborate our quantum self-consistent approach using the quantum master equation that can deal with anharmonicity exactly, but is limited to the weak system-bath coupling regime. Finally, in order demonstrate its strength we apply the quantum self-consistent approach to study thermal rectification in a weakly coupled two segment anharmonic system.

preprint2016arXiv

Spin-dependent Seebeck Effect in Aharonov-Bohm Rings with Rashba and Dresselhaus Spin-orbit Interactions

We theoretically investigate the spin-dependent Seebeck effect in an Aharonov-Bohm mesoscopic ring in the presence of both Rashba and Dresselhaus spin-orbit interactions under magnetic flux perpendicular to the ring. We apply the Green's function method to calculate the spin Seebeck coefficient employing the tight-binding Hamiltonian. It is found that the spin Seebeck coefficient is proportional to the slope of the energy-dependent transmission coefficients. We study the strong dependence of spin Seebeck coefficient on the Fermi energy, magnetic flux, strength of spin-orbit coupling, and temperature. Maximum spin Seebeck coefficients can be obtained when the strengths of Rashba and Dresselhaus spin-orbit couplings are slightly different. The spin Seebeck coefficient can be reduced by increasing temperature and disorder.

preprint2016arXiv

Superior thermal conductivity in suspended bilayer hexagonal boron nitride

We reported the basal-plane thermal conductivity in exfoliated bilayer hexagonal boron nitride h-BN that was measured using suspended prepatterned microstructures. The h-BN sample suitable for thermal measurements was fabricated by dry-transfer method, whose sample quality, due to less polymer residues on surfaces, is believed to be superior to that of PMMA-mediated samples. The measured room temperature thermal conductivity is around 484 Wm-1K-1(+141 Wm-1K-1/ -24 Wm-1K-1) which exceeds that in bulk h-BN, providing experimental observation of the thickness-dependent thermal conductivity in suspended few-layer h-BN.

preprint2016arXiv

Variational approach to renormalized phonon in momentum-nonconserving nonlinear lattices

A previously proposed variational approach for momentum-conserving systems [J. Liu et.al., Phys. Rev. E 91, 042910 (2015)] is extended to systematically investigate general momentum-nonconserving nonlinear lattices. Two intrinsic identities characterizing optimal reference systems are revealed, which enables us to derive explicit expressions for optimal variational parameters. The resulting optimal harmonic reference systems provide information for the band gap as well as the dispersion of renormalized phonons in nonlinear lattices. As a demonstration, we consider the one-dimensional ϕ^?4 lattice. By combining the transfer integral operator method, we show that the phonon band gap endows a simple power-law temperature dependence in the weak stochasticity regime where predicted dispersion is reliable by comparing with numerical results. In addition, an exact relation between ensemble averages of the ϕ^?4 lattice in the whole temperature range is found, regardless of the existence of the strong stochasticity threshold.

preprint2015arXiv

1D momentum-conserving systems: the conundrum of anomalous versus normal heat transport

Transport and diffusion of heat in one dimensional (1D) nonlinear systems which {\it conserve momentum} is typically thought to proceed anomalously. Notable exceptions, however, exist of which the rotator model is a prominent case. Therefore, the quest arises to identify the origin of manifest anomalous transport in those low dimensional systems. Here, we develop the theory for both, momentum/heat diffusion and its corresponding momentum/heat transport. We demonstrate that the second temporal derivative of the mean squared deviation of the momentum spread is proportional to the equilibrium correlation of the total momentum flux. This result in turn relates, via the integrated momentum flux correlation, to an effective viscosity, or equivalently, to the underlying momentum diffusivity. We put forward the intriguing hypothesis that a fluid-like momentum dynamics with a {\it finite viscosity} causes {\it normal} heat transport; its corollary being that superdiffusive momentum diffusion with an intrinsic {\it diverging viscosity} in turn yields {\it anomalous} heat transport. This very hypothesis is corroborated over wide extended time scales by use of precise molecular dynamics simulations. The numerical validation of the hypothesis involves three distinct archetype classes of nonlinear pair-interaction potentials: (i) a globally bounded pair interaction (the noted rotator model), (ii) unbounded interactions acting at large distances (the Fermi-Pasta-Ulam $β$ model, or the rotator model amended with harmonic pair interactions) and (iii), a pair interaction potential being unbounded at short distances while displaying an asymptotic free part (Lennard-Jones model).

preprint2015arXiv

Boosting thermoelectric efficiency using time-dependent control

Thermoelectric efficiency is defined as the ratio of power delivered to the load of a device to the rate of heat flow from the source. Till date, it has been studied in presence of thermodynamic constraints set by the Onsager reciprocal relation and the second law of thermodynamics that severely bottleneck the thermoelectric efficiency. In this study, we propose a pathway to bypass these constraints using a time-dependent control and present a theoretical framework to study dynamic thermoelectric transport in the far from equilibrium regime. The presence of a control yields the sought after substantial efficiency enhancement and importantly a significant amount of power supplied by the control is utilised to convert the wasted-heat energy into useful-electric energy. Our findings are robust against nonlinear interactions and suggest that external time-dependent forcing, which can be incorporated with existing devices, provides a beneficial scheme to boost thermoelectric efficiency.

preprint2015arXiv

Direction dependent thermal conductivity of monolayer phosphorene: parameterization of Stillinger-Weber potential and molecular dynamics study

A Stillinger-Weber interatomic potential is parameterized for phosphorene. It well reproduces the crystal structure, cohesive energy and phonon dispersion predicted by first-principles calculations. The thermal conductivity of phosphorene is further explored by equilibrium molecular dynamics simulations adopting the optimal set of potential parameters. At room temperature, the intrinsic thermal conductivities along zigzag and armchair directions are about 152.7 and 33.0 W/mK, respectively, with a large anisotropy ratio of five. The remarkably directional dependence of thermal conductivity in phosphorene, consistent with previous reports, is mainly due to the strong anisotropy of phonon group velocities, and weak anisotropy of phonon lifetimes as revealed by lattice dynamics calculations. Moreover, the effective phonon mean free paths at zigzag and armchair directions are about 141.4 and 43.4nm, respectively.

preprint2015arXiv

Heat conduction and energy diffusion in momentum-conserving 1D 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 {\it 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.

preprint2015arXiv

Invisible sensor: Simultaneous sensing and camouflaging in multiphysical fields

To manipulate various types of physical signals in one single device has long captivated the attention of scientists and engineers. This however is very challenging, if not impossible, even for emerging metamaterials. Up to date, many artificial materials have been proposed, theoretically and (or) experimentally, for manipulating various waves/signals on a one-function-one-device basis. In this work, for the very first time, we employ undecorated natural materials to experimentally demonstrate a simultaneous camouflage for thermal current and electric dc current on the same device. It demonstrates how ingenuity can overcome the limitations of natural material systems without the need for complex decoration to impart inhomogeneous and (or) anisotropic properties, which was previously considered impossible to accomplish except by using metamaterials.

preprint2015arXiv

Renormalized phonons in nonlinear lattices: A variational approach

We propose a variational approach to study renormalized phonons in momentum conserving nonlinear lattices with either symmetric or asymmetric potentials. To investigate the influence of pressure to phonon properties, we derive an inequality which provides both the lower and upper bound of the Gibbs free energy as the associated variational principle. This inequality is a direct extension to the Gibbs-Bogoliubov inequality. Taking the symmetry effect into account, the reference system for the variational approach is chosen to be harmonic with an asymmetric quadratic potential which contains variational parameters. We demonstrate the power of this approach by applying it to one dimensional nonlinear lattices with a symmetric or asymmetric Fermi-Pasta- Ulam type potential. For a system with a symmetric potential and zero pressure, we recover existing results. For other systems which beyond the scope of existing theories, including those having the symmetric potential and pressure, and those having the asymmetric potential with or without pressure, we also obtain accurate sound velocity.

preprint2015arXiv

Stretch diffusion and heat conduction in 1D nonlinear lattices

In the study of 1D nonlinear Hamiltonian lattices, the conserved quantities play an important role in determining the actual behavior of heat conduction. Besides the total energy, total momentum and total stretch could also be conserved quantities. In microcanonical Hamiltonian dynamics, the total energy is always conserved. It was recently argued by Das and Dhar that whenever stretch (momentum) is not conserved in a 1D model, the momentum (stretch) and energy fields exhibit normal diffusion. In this work, we will systematically investigate the stretch diffusions for typical 1D nonlinear lattices. No clear connection between the conserved quantities and heat conduction can be established. The actual situation is more complicated than what Das and Dhar claimed.

preprint2015arXiv

Temperature and Frequency Dependent Mean Free Paths of Renormalized Phonons in Nonlinear Lattices

In the regime of strong nonlinearity, the validity of conventional perturbation based phonon transport theories is questionable. In particular, the renormalized phonons instead of phonons are responsible for heat transport in nonlinear lattices. In this work, we directly study the temperature and frequency dependent Mean Free Path (MFP) of renormalized phonons with the newly developed numerical tuning fork method. The typical 1D nonlinear lattices such as Fermi-Pasta-Ulam $β$ (FPU-$β$) lattice and $ϕ^4$ lattice are investigated in details. It is found that the MFPs are inversely proportional to the frequencies of renormalized phonons rather than the square of phonon frequencies predicted by existing phonon scattering theory.

preprint2015arXiv

Temperature dependence of thermal conductivities of coupled rotator lattice and the momentum diffusion in standard map

In contrary to other 1D momentum-conserving lattices such as the Fermi-Pasta-Ulam $β$ (FPU-$β$) lattice, the 1D coupled rotator lattice is a notable exception which conserves total momentum while exhibits normal heat conduction behavior. The temperature behavior of the thermal conductivities of 1D coupled rotator lattice had been studied in previous works trying to reveal the underlying physical mechanism for normal heat conduction. However, two different temperature behaviors of thermal conductivities have been claimed for the same coupled rotator lattice. These different temperature behaviors also intrigue the debate whether there is a phase transition of thermal conductivities as the function of temperature. In this work, we will revisit the temperature dependent thermal conductivities for the 1D coupled rotator lattice. We find that the temperature dependence follows a power law behavior which is different with the previously found temperature behaviors. Our results also support the claim that there is no phase transition for 1D coupled rotator lattice. We also give some discussion about the similarity of diffusion behaviors between the 1D coupled rotator lattice and the single kicked rotator also called the Chirikov standard map.

preprint2015arXiv

Thermoelectric transport through a quantum nanoelectromechanical system and its backaction

We present a comprehensive study of thermoelectric transport properties of a quantum nanoelectromechanical system (NEMS) described by a single-electron-transistor (SET) coupled to a quantum nanomechanical resonator (NR). The effects of a quantum NR on the electronic current are investigated with special emphasis on how the SET-NR coupling strength plays a role in such a NEMS. We find that the SET-NR coupling is not only able to suppress or enhance the thermoelectric current but can also switch its direction. The effect of the NR on the thermoelectric coefficients of the SET are studied and we find that even a small SET-NR coupling could dramatically suppress the figure of merits ZT . Lastly, we investigate the backaction of electronic current on the NR and possible routes of heating or cooling the NR are discussed. We find that by appropriately tuning the gate voltage the backaction can be eliminated, which could find possible applications to enhance the sensitivity of detection devices.

preprint2015arXiv

Transient unidirectional energy flow and diode-like phenomenon induced by non-Markovian environments

Relying on an exact time evolution scheme, we identify a novel transient energy transfer phe- nomenon in an exactly-solvable quantum microscopic model consisting of a three-level system coupled to two non-Markovian zero-temperature bosonic baths through two separable quantum channels. The dynamics of this model can be solved exactly using the quantum-state-diffusion equation formalism, demonstrating finite intervals of unidirectional energy flow across the system, typically, from the non-Markovian environment towards the more Markovian bath. Furthermore, when introducing a spatial asymmetry into the system, an analogue of the rectification effect is realized. In the long time limit, the dynamics arrives at a stationary state and the effects recede. Understanding temporal characteristics of directional energy flow will aid in designing microscopic energy transfer devices.

preprint2014arXiv

Co-existence of size-dependent and size-independent thermal conductivities in single layer black phosphorus

Thermal conductivity of single layer black phosphorus (BP) is investigated by combining density functional calculations and Peierls-Boltzmann transport equation. Differing from isotropic and divergent thermal conductivities in two-dimensional graphene and MoS$_2$, an compelling co-existence of size-dependent and size-independent thermal conductivities are discovered for single layer BP along zigzag (ZZ) and armchair (AM) direction, respectively. Besides, thermal conductivities in single layer BP are found to be highly anisotropic because of orientation dependent group velocities, e.g., thermal conductivities at 300 K are 83.5 and 24.3 W/m-K along ZZ and AM directions for single layer BP with a size of 10 $μm$, respectively.

preprint2014arXiv

Competing for Attention in Social Media under Information Overload Conditions

Although the many forms of modern social media have become major channels for the dissemination of information, they are becoming overloaded because of the rapidly-expanding number of information feeds. We analyze the expanding user-generated content in Sina Weibo, the largest micro-blog site in China, and find evidence that popular messages often follow a mechanism that differs from that found in the spread of disease, in contrast to common believe. In this mechanism, an individual with more friends needs more repeated exposures to spread further the information. Moreover, our data suggest that in contrast to epidemics, for certain messages the chance of an individual to share the message is proportional to the fraction of its neighbours who shared it with him/her. Thus the greater the number of friends an individual has the greater the number of repeated contacts needed to spread the message, which is a result of competition for attention. We model this process using a fractional susceptible infected recovered (FSIR) model, where the infection probability of a node is proportional to its fraction of infected neighbors. Our findings have dramatic implications for information contagion. For example, using the FSIR model we find that real-world social networks have a finite epidemic threshold. This is in contrast to the zero threshold that conventional wisdom derives from disease epidemic models. This means that when individuals are overloaded with excess information feeds, the information either reaches out the population if it is above the critical epidemic threshold, or it would never be well received, leading to only a handful of information contents that can be widely spread throughout the population.

preprint2014arXiv

Exchange fluctuation theorem for heat transport between multi-terminal harmonic systems

We study full counting statistics for transferred heat and entropy production between multi-terminal systems in absence of a finite junction. The systems are modelled as collections of coupled harmonic oscillators which are kept at different equilibrium temperatures and are connected via arbitrary time dependent couplings. Following consistent quantum framework and two-time measurement concept we obtain analytical expressions for the generalized cumulant generating function. We discuss transient and steady-state fluctuation theorems for the transferred quantities. We also address the effect of coupling strength on the exchange fluctuation theorem.

preprint2014arXiv

Inhomogeneous Thermal Conductivity Enhances Thermoelectric Cooling

We theoretically investigate the enhancement of thermoelectric cooling performance in thermoelectric devices made of materials with inhomogeneous thermal conductivity, beyond the usual practice of enhancing thermoelectric figure of merit ZT. The dissipation of Joule heat in such thermoelectric devices is asymmetric which can give rise to better thermoelectric cooling performance. Although the thermoelectric figure of merit and the coefficient-of-performance are only slightly enhanced, both the maximum cooling power and the maximum cooling temperature difference can be enhanced significantly. This finding can be used to increase the heat absorption at the cold end. The asymmetric dissipation of Joule heat also leads to thermal rectification.

preprint2014arXiv

Invisibility's Flicker: Detecting Thermal Cloaks via Transient Effects

Recent research on the development of a thermal cloak has concentrated on engineering an inhomogeneous thermal conductivity and homogeneous volumetric heat capacity. While the perfect cloak of inhomogeneous $κ$ and $ρc_p$ is known to be exact (no signals scattering or penetrating to the cloak's interior), no such analysis has been considered for this case. Using analytic, computational, and experimental techniques, we demonstrate that these approximate cloaks are detectable. Although they work as perfect cloaks in the steady-state, their transient (time-dependent) response is imperfect and a detectable amount of heat is scattered. This is sufficient to determine the presence of a cloak and any heat source it contains, but the material composition hidden within the cloak is not detectable in practice.

preprint2014arXiv

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, when temperature at 300K, thermal conductivity keeps increasing and remains logarithmic divergence 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 into thermal transport in two-dimensional materials.

preprint2014arXiv

Manipulating Steady Heat Conduction by Sensu-shaped Thermal Metamaterials

The ability to design the control of heat flow has innumerable benefits in the design of electronic systems such as thermoelectric energy harvesters, solid-state lighting, and thermal imagers, where the thermal design plays a key role in performance and device reliability. However, to realize one advanced control function of thermal flux, one needs to design one sophisticated, multilayered and inhomogeneous thermal structure with different composition/shape at different regions of one device. In this work, we employ one identical sensu-unit with facile natural composition to experimentally realize a new class of thermal metamaterials for controlling thermal conduction (e.g., thermal concentrator, focusing/resolving, uniform heating), only resorting to positioning and locating the same unit element of sensu-shape structure. The thermal metamaterial unit and the proper arrangement of multiple identical units are capable of transferring, redistributing and managing thermal energy in a versatile fashion. It is also shown that our sensu-shape unit elements can be used in manipulating dc currents without any change in the layout for the thermal counterpart. The proposed scheme can also be applied to control dc electric currents and dc magnetic fields that governed by Laplace equation. These could markedly enhance the capabilities in thermal sensing, thermal imaging, thermal-energy storage, thermal packaging, thermal therapy, and more domains beyond.

preprint2014arXiv

Nanoscale Graphene Disk: A Natural Functionally Graded Material --The Thermal Conductivity of Nanoscale Graphene Disk by Molecular Dynamics Simulation

In this letter, we investigate numerically (by non-equilibrium molecular dynamics) and analytically the thermal conductivity of nanoscale graphene disks (NGDs), and discussed the possibility to realize FGM with only one material, NGDs. We found that the NGD has a graded thermal conductivity and can be used as FGM in a large temperature range. Moreover, we show the dependent of NGDs' thermal conductivity on radius and temperature. Our study may inspire experimentalists to develop NGD based FGMs and help heat removal of hot spots on chips by graphene.

preprint2014arXiv

Phonon nanocapacitor for storage and lasing of terahertz lattice waves

We introduce a novel ultra-compact nanocapacitor of coherent phonons formed by high-finesse interference mirrors based on atomic-scale semiconductor metamaterials. Our molecular dynamics simulations show that the nanocapacitor stores THz monochromatic lattice waves, which can be used for phonon lasing - the emission of coherent phonons. Either one- or two-color phonon lasing can be realized depending on the geometry of the nanodevice. The two color regimes of the capacitor originates from the distinct transmittance dependance on the phonon wave packet incident angle for the two phonon polarizations at their respective resonances. Phonon nanocapacitor can be charged by cooling the sample equilibrated at room temperature or by the pump-probe technique. The nanocapacitor can be discharged by applying tunable reversible strain, resulting in the emission of coherent THz acoustic beams.

preprint2014arXiv

Predicting market instability: New dynamics between volume and volatility

Econophysics and econometrics agree that there is a correlation between volume and volatility in a time series. Using empirical data and their distributions, we further investigate this correlation and discover new ways that volatility and volume interact, particularly when the levels of both are high. We find that the distribution of the volume-conditional volatility is well fit by a power-law function with an exponential cutoff. We find that the volume-conditional volatility distribution scales with volume, and collapses these distributions to a single curve. We exploit the characteristics of the volume-volatility scatter plot to find a strong correlation between logarithmic volume and a quantity we define as local maximum volatility (LMV), which indicates the largest volatility observed in a given range of trading volumes. This finding supports our empirical analysis showing that volume is an excellent predictor of the maximum value of volatility for both same-day and near-future time periods. We also use a joint conditional probability that includes both volatility and volume to demonstrate that invoking both allows us to better predict the largest next-day volatility than invoking either one alone.

preprint2014arXiv

Robust acoustic cloaking with density-near-zero materials

Using the phenomenon of extraordinary sound transmission (EST) through artificially constructed density-near-zero (DNZ) materials, we propose an alternative approach to realize a very robust acoustic cloaking, which includes one absorber and one projector connected by an energy channel. The elementary unit cell is made of one single-piece homogeneous copper, but it well maintains both the original planar wavefront and the nearly perfect one-dimensional transmission, in the presence of any inserted objects. Furthermore, the overall cloaked space can be increased, designed, and distributed by arbitrarily assembling additional DNZ cells without any upper limit in the total cloaking volume. We demonstrate acoustic cloaking by different combinations of DNZ cells in free space as well as in acoustic waveguides, which enables any objects inside the cells imperceptible along undistorted sound paths.

preprint2014arXiv

Significant Reduction of Graphene Thermal Conductivity by Phononic Crystal Structure

We studied the thermal conductivity of graphene phononic crystal (GPnC), also named as graphene nanomesh, by molecular dynamics simulations. The dependences of thermal conductivity of GPnCs on both length and temperature are investigated. It is found that the thermal conductivity of GPnCs is significantly lower than that of graphene and can be efficiently tuned by changing the porosity and period length. For example, the ratio of thermal conductivity of GPnC to thermal conductivity of graphene can be changed from 0.1 to 0.01 when the porosity is changed from about 21% to 65%. The phonon participation ratio spectra reveal that more phonon modes are localized in GPnCs with larger porosity. Our results suggest that creating GPnCs is a valuable method to efficiently manipulate the thermal conductivity of graphene.

preprint2014arXiv

Spin Seebeck Effect in Asymmetric Four-Terminal Systems with Rashba Spin-Orbit Coupling

We propose a new type of the spin Seebeck effect (SSE) emerging from the Rashba spin-orbit coupling in asymmetric four-terminal electron systems. This system generates spin currents or spin voltages along the longitudinal direction parallel to the temperature gradient in the absence of magnetic fields. The remarkable result arises from the breaking of reflection symmetry along the transverse direction. In the meantime, the SSE along the transverse direction, so-called the spin Nernst effect, with spin currents or spin voltages perpendicular to the temperature gradient can be simultaneously realized in our system. We further find that it is possible to use the temperature differences between four leads to tune the spin Seebeck coefficients.

preprint2014arXiv

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

The temperature-dependent thermal conductivities of one-dimensional nonlinear Klein-Gordon lattices with soft on-site potential (soft-KG) have been investigated systematically. Similar to the previously studied hard-KG lattices, the existence of renormalized phonons has also been confirmed in soft-KG lattices. In particular, the temperature-dependence of renormalized phonon frequency predicted by a classical field theory has been verified by detailed numerical simulations. However, the thermal conductivities of soft-KG lattices exhibit opposite trend in the temperature dependence in comparison with the hard-KG lattices. The interesting thing is that both the temperature-dependent thermal conductivities of 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 the function of temperature are only determined by the exponents of the soft or hard on-site potentials. These theoretical predictions have been consistently verified very well by extensive numerical simulations.

preprint2014arXiv

Thermal Boundary Conductance Across Metal-Nonmetal Interfaces: Effects of Electron-Phonon Coupling both in Metal and at Interface

We theoretically investigate the thermal boundary conductance across metal-nonmetal interfaces in the presence of the electron-phonon coupling not only in metal but also at interface. The thermal energy can be transferred from metal to nonmetal via three channels: (1) the phonon-phonon coupling at interface; (2) the electron-phonon coupling at interface; and (3) the electron-phonon coupling within metal and then subsequently the phonon-phonon coupling at interface. We find that these three channels can be described by an equivalent series-parallel thermal resistor network, based on which we derive out the analytic expression of the thermal boundary conductance. We then exemplify different contributions from each channel to the thermal boundary conductance in three typical interfaces: Pb-diamond, Ti-diamond, and TiN-MgO. Our results reveal that the competition among above channels determines the thermal boundary conductance.

preprint2014arXiv

Thermal Interface Conductance between Aluminum and Silicon by Molecular Dynamics Simulations

The thermal interface conductance between Al and Si was simulated by a non-equilibrium molecular dynamics method. In the simulations, the coupling between electrons and phonons in Al are considered by using a stochastic force. The results show the size dependence of the interface thermal conductance and the effect of electron-phonon coupling on the interface thermal conductance. To understand the mechanism of interface resistance, the vibration power spectra are calculated. We find that the atomic level disorder near the interface is an important aspect of interfacial phonon transport, which leads to a modification of the phonon states near the interface. There, the vibrational spectrum near the interface greatly differs from the bulk. This change in the vibrational spectrum affects the results predicted by AMM and DMM theories and indicates new physics is involved with phonon transport across interfaces. Keywords:

preprint2014arXiv

Thermoelectric Properties of Nanoscale three dimensional Si Phononic Crystal

The thermoelectric properties of n type nanoscale three dimensional (3D) Si phononic crystals (PnCs) with spherical pores are studied. Density functional theory and Boltzmann transport equation under the relaxation time approximation are applied to study the electronic transport coefficients, electrical conductivity, Seebeck coefficient and electronic thermal conductivity. We found that the electronic transport coefficients in 3D Si PnC at room temperature (300 K) change very little compared with that of Si, for example, electrical conductivity and electronic thermal conductivity is decreased by 0.26 to 0.41 and 0.39 to 0.55 depending on carrier concentration, respectively, and the Seebeck coefficient is similar to that of bulk Si. However, the lattice thermal conductivity of 3D Si PnCs with spherical pores is decreased by a factor of 500 calculated by molecular dynamics methods, leading to the ZT of 0.76, which is about 30 times of that of porous Si. This work indicates that 3D Si PnC is a promising candidate for high efficiency thermoelectric materials.

preprint2014arXiv

Triggering waves in nonlinear lattices: Quest for anharmonic phonons and corresponding mean free paths

Guided by a stylized experiment we develop a self-consistent anharmonic phonon concept for nonlinear lattices which allows for explicit "visualization." The idea uses a small external driving force which excites the front particles in a nonlinear lattice slab and subsequently one monitors the excited wave evolution using molecular dynamics simulations. This allows for a simultaneous, direct determination of the existence of the phonon mean free path with its corresponding anharmonic phonon wavenumber as a function of temperature. The concept for the mean free path is very distinct from known prior approaches: the latter evaluate the mean free path only indirectly, via using both, a scale for the phonon relaxation time and yet another one for the phonon velocity. Notably, the concept here is neither limited to small lattice nonlinearities nor to small frequencies. The scheme is tested for three strongly nonlinear lattices of timely current interest which either exhibit normal or anomalous heat transport.

preprint2014arXiv

Wave-packet rectification in nonlinear electronic systems: A tunable Aharonov-Bohm diode

Rectification of electron wave-packets propagating along a quasi-one dimensional chain is commonly achieved via the simultaneous action of nonlinearity and longitudinal asymmetry, both confined to a limited portion of the chain termed wave diode. However, it is conceivable that, in the presence of an external magnetic field, spatial asymmetry perpendicular to the direction of propagation suffices to ensure rectification. This is the case of a nonlinear ring-shaped lattice with different upper and lower halves (diode), which is attached to two elastic chains (leads). The resulting device is mirror symmetric with respect to the ring vertical axis, but mirror asymmetric with respect to the chain direction. Wave propagation along the two diode paths can be modeled for simplicity by a discrete Schrödinger equation with cubic nonlinearities. Numerical simulations demonstrate that, thanks to the Aharonov-Bohm effect, such a diode can be operated by tuning the magnetic flux across the ring.

preprint2013arXiv

Acoustic Wavefront Manipulation: Impedance Inhomogeneity and Extraordinary Reflection

Optical wavefront can be manipulated by interfering elementary beams with phase inhomogeneity. Therefore a surface allowing huge, abrupt and position-variant phase change would enable all possibilities of wavefront engineering. However, one may not have the luxury of efficient abrupt-phase-changing materials in acoustics. This motivates us to establish a counterpart mechanism for acoustics, in order to empower the wide spectrum of novel acoustic applications. Remarkably, the proposed impedance-governed generalized Snell's law (IGSL) of reflection is distinguished from that in optics. Via the manipulation of inhomogeneous acoustic impedance, extraordinary reflection can be tailored for unprecedented wavefront manipulation while ordinary reflection can be surprisingly switched on or off. Our results may power the acoustic-wave manipulation and engineering. We demonstrate novel acoustic applications by planar surfaces designed with IGSL.

preprint2013arXiv

Anomalous Heat Diffusion

Consider anomalous energy spread in solid phases, i.e., $MSD= \int (x -{< x >}_E)^2 ρ_E(x,t)dx \propto t^β$, as induced by a small initial excess energy perturbation distribution $ρ_{E}(x,t=0)$ away from equilibrium. The associated total thermal equilibrium heat flux autocorrelation function $C_{JJ}(t)$ is shown to obey rigorously the intriguing relation, $d^2 MSD/dt^2 = 2C_{JJ}(t)/(k_BT^2c)$, where $c$ is the specific volumetric heat capacity. Its integral assumes a time-local Helfand-moment relation; i.e. $ dMSD/dt|_{t=t_s} = 2/(k_BT^2c)\int_0^{t_s} C_{JJ}(s)ds$, where the chosen cut-off time $t_s$ is determined by the maximal signal velocity for heat transfer. Given the premise that the averaged nonequilibrium heat flux is governed by an anomalous heat conductivity, energy diffusion scaling necessarily determines a corresponding anomalous thermal conductivity scaling behavior.

preprint2013arXiv

Bilayer Isotropic Thermal Cloak

Invisibility has attracted intensive research in various communities, e.g., optics, electromagnetics, acoustics, thermodynamics, etc. However, the most of them have only been experimentally achieved by virtue of simplified approaches due to their inhomogeneous and extreme parameters imposed by transformation-optic method, and usually require challenging realization with metamaterials. In this paper, we demonstrate an advanced bilayer thermal cloak with naturally available materials first time. This scheme, directly from thermal conduction equation, has been validated as an exact cloak rather than a reduced one, and we experimentally confirmed its perfect performance (heat-front maintenance and heat protection) in an actual setup. The proposed scheme may open a new avenue to control the diffusive heat flow in ways inconceivable with phonons.

preprint2013arXiv

Classical Heat Transport in Anharmonic Molecular Junctions: Exact Solutions

We study full counting statistics for classical heat transport through anharmonic/nonlinear molecular junctions formed by interacting oscillators. Analytical result of the steady state heat flux for an overdamped anharmonic junction with arbitrary temperature bias is obtained. It is found that the thermal conductance can be expressed in terms of temperature dependent effective force constant. The role of anharmonicity is identified. We also give the general formula for the second cumulant of heat in steady state, as well as the average geometric heat flux when two temperatures are modulated adiabatically. We present an anharmonic example for which all cumulants for heat can be obtained exactly. For a bounded single oscillator model with mass we found that the cumulants are independent of the nonlinear potential.

preprint2013arXiv

Creation of Tunable Homogeneous Thermal Cloak with Constant Conductivity

Invisible cloak has long captivated the popular conjecture and attracted intensive research in various communities of wave dynamics, e.g., optics, electromagnetics, acoustics, etc. However, their inhomogeneous and extreme parameters imposed by transformation-optic method will usually require challenging realization with metamaterials, resulting in narrow bandwidth, loss, polarization-dependence, etc. On the contrary, we demonstrate that tunable thermodynamic cloak can be achieved with homogeneous and finite conductivity only employing naturally available materials. The controlled localization of thermal distribution inside the coating layer has been presented, and it shows that an incomplete cloak can function perfectly. Practical realization of such homogeneous thermal cloak has been suggested by using two naturally occurring conductive materials, which provides an unprecedentedly plausible way to flexibly realize flexible thermal cloak and manipulate thermal flow.

preprint2013arXiv

Cumulants of heat transfer across nonlinear quantum systems

We consider thermal conduction across a general nonlinear phononic junction. Based on two-time observation protocol and the field theoretical/algebraic method, the cumulants of the heat transferred in both transient and steady-state regimes are studied on an equal footing, and practical formulae for the calculation of the cumulant generating function of heat transfer are obtained. As an application, the developed general formalism is used to study anharmonic effects on fluctuation of steady-state heat transfer across a single-site junction with a quartic nonlinear on-site pinning potential. An explicit nonlinear modification to cumulant generating function exact up to the first order is given, in which Gallavotti-Cohen fluctuation symmetry is verified. Numerically a self-consistent procedure is introduced, which works well for strong nonlinearity.

preprint2013arXiv

Extreme Low Thermal Conductivity in Nanoscale 3D Si Phononic Crystal with Spherical Pores

Thermoelectric material provides a high hope for converting harmful and useless heat into useful energy, electricity. It can also be used as solid state Peltier coolers in integrated circuits, an outstanding challenge for electronic engineers. The desire for a high efficient thermoelectric material has never been so keen. Although many works have been done, we are still far from having a recipe for thermoelectric materials. Nano-structuring provides an effective way to increase figure of merit (ZT) by reducing the thermal conductivity without affecting electronic property.1 Here, we propose a novel nanoscale three-dimensional (3D) Si phononic crystal (PnC) with spherical pores, which can reduce the thermal conductivity of bulk Si by a factor up to 10,000 times at room temperature. The extreme-low thermal conductivity could lead to a larger value of ZT than unity. The thermal conductivity changes little when temperature increases from room temperature to 1100 K. The phonon participation ratio spectra show there are more phonon localizations as the porosity of PnC increases.

preprint2013arXiv

Nonlinearity enhanced interfacial thermal conductance and rectification

We study the nonlinear interfacial thermal transport across atomic junctions by the quantum self-consistent mean field (QSCMF) theory based on nonequilibrium Green's function approach; the QSCMF theory we propose is very precise and matches well with the exact results from quantum master equations. The nonlinearity at the interface is studied by effective temperature dependent interfacial coupling calculated from the QSCMF theory. We find that nonlinearity can provide an extra channel for phonon transport in addition to the phonon scattering which usually blocks heat transfer. For weak linearly coupled interface, the nonlinearity can enhance the interfacial thermal transport; with increasing nonlinearity or temperature, the thermal conductance shows nonmonotonical behavior. The interfacial nonlinearity also induces thermal rectification, which depends on the mismatch of the two leads and also the interfacial linear coupling.

preprint2013arXiv

Reduction of Thermal Conductivity by Nanoscale 3D Phononic Crystal

The thermal conductivity of nanostructures needs to be as small as possible so that it will have a greater efficiency for solid-state electricity generation/refrigeration by thermoelectrics. We studied how the period length and the mass ratio affect the thermal conductivity of isotopic nanoscale 3D phononic crystal of Si. Simulation results by equilibrium molecular dynamics show isotopic nanoscale 3D phononic crystal can make a significance reduction on the thermal conductivity of bulk Si at high temperature,1000 K. Size and mass effects are obvious in manipulating thermal conductivity. The thermal conductivity decreases as the period length and mass ratio increases. The phonon dispersion curves show the decrease of group velocities in 3D phononic crystals. The phonon's localization and band gap is clearly shown in spectra of normalized inverse participation ratio in nanoscale 3D phononic crystal structure.

preprint2013arXiv

Substrate Coupling Suppresses Size Dependence of Thermal Conductivity in Supported Graphene

Thermal conductivity $κ$ of both suspended and supported graphene has been studied by using molecular dynamics simulations. Obvious length dependence is observed in $κ$ of suspended single-layer graphene (SLG), while $κ$ of supported SLG is insensitive to the length. The simulation result of room temperature $κ$ of supported SLG is in good agreement with experimental value. In contrast to the decrease in $κ$ induced by inter-layer interaction in suspended few-layer graphene (FLG), $κ$ of supported FLG is found to increase rapidly with the layer thickness, reaching about 90% of that of bulk graphite at six layers, and eventually saturates at the thickness of 13.4 nm. More interestingly, unlike the remarkable substrate dependent $κ$ in SLG, the effect of substrate on thermal transport is much weaker in FLG. The underlying physics is investigated and presented.

preprint2013arXiv

Thermal transport across metal-insulator interface via electron-phonon interaction

The thermal transport across metal-insulator interface can be characterized by electron-phonon interaction through which an electron lead is coupled to a phonon lead if phonon-phonon coupling at the interface is very weak. We investigate the thermal conductance and rectification flowing between the electron part and the phonon part using nonequilibrium Green's function method. It is found that the thermal conductance has a nonmonotonic behavior as a function of average temperature or the coupling strength between the phonon leads in the metal part and the insulator one. The metal-insulator interface shows evident thermal rectification effect, which can reverse with changing of average temperature or the electron-phonon coupling.

preprint2013arXiv

Thermal transport in nanostructures

This review summarizes recent studies of thermal transport in nanoscaled semiconductors. Different from bulk materials, new physics and novel thermal properties arise in low dimensional nanostructures, such as the abnormal heat conduction, the size dependence of thermal conductivity, phonon boundary/edge scatterings. It is also demonstrated that phonons transport super-diffusively in low dimensional structures, in other words, Fourier's law is not applicable. Based on manipulating phonons, we also discuss envisioned applications of nanostructures in a broad area, ranging from thermoelectrics, heat dissipation to phononic devices.

preprint2013arXiv

Topological Magnon Insulator in Insulating Ferromagnet

In the ferromagnetic insulator with Dzyaloshinskii-Moriya interaction, we theoretically predict and numerically verify a topological magnon insulator, where the charge-free magnon is topologically protected to transport along the edge while it is insulating in the bulk. Within the bulk band gaps, edge states form a connected loop as a $4π$- or $8π$-period Möbius strip in the wave vector space. As a consequence, the chiral energy current traveling along the corresponding edge is topologically protected from defects or disorders. Using the nonequilibrium Green's function method, we demonstrate that the energy current carried by magnons with energy in the bulk gap localizes at edges and prefers to travel along one edge in only one direction at nonequilibrium steady state. Our prediction about topological magnon insulator could be observed at a wide energy range in the thin film of the insulating ferromagnet, such as Lu$_2$V$_2$O$_7$.

preprint2013arXiv

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.

preprint2012arXiv

Anomalous Heat Conduction and Anomalous Diffusion in Low Dimensional Nanoscale Systems

Thermal transport is an important energy transfer process in nature. Phonon is the major energy carrier for heat in semiconductor and dielectric materials. In analogy to Ohm's law for electrical conductivity, Fourier's law is a fundamental rule of heat transfer in solids. It states that the thermal conductivity is independent of sample scale and geometry. Although Fourier's law has received great success in describing macroscopic thermal transport in the past two hundreds years, its validity in low dimensional systems is still an open question. Here we give a brief review of the recent developments in experimental, theoretical and numerical studies of heat transport in low dimensional systems, include lattice models, nanowires, nanotubes and graphenes. We will demonstrate that the phonon transports in low dimensional systems super-diffusively, which leads to a size dependent thermal conductivity. In other words, Fourier's law is breakdown in low dimensional structures.

preprint2012arXiv

Controlling complex networks: How much energy is needed?

The outstanding problem of controlling complex networks is relevant to many areas of science and engineering, and has the potential to generate technological breakthroughs as well. We address the physically important issue of the energy required for achieving control by deriving and validating scaling laws for the lower and upper energy bounds. These bounds represent a reasonable estimate of the energy cost associated with control, and provide a step forward from the current research on controllability toward ultimate control of complex networked dynamical systems.

preprint2012arXiv

Geometric Heat Flux for Classical Thermal Transport in Interacting Open Systems

We study classical heat conduction in a dissipative open system composed of interacting oscillators. By exactly solving a twisted Fokker-Planck equation which describes the full counting statistics of heat flux flowing through the system, we identify the geometric-phase-like effect and examine its impact on the classical heat transport. Particularly, we find that the nonlinear interaction as well as the closely related temperature-dependence of system-parameters are crucial in manifesting the geometric-phase contribution of heat flux. Finally, we propose an electronic experiment based on RC circuits to verify our theoretical predictions.

preprint2012arXiv

How Does Folding Modulate Thermal Conductivity of Graphene

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

preprint2012arXiv

Impacts of Atomistic Coating on Thermal Conductivity of Germanium Nanowires

By using non-equilibrium molecular dynamics simulations, we demonstrated that thermal conductivity of Germanium nanowires can be reduced more than 25% at room temperature by atomistic coating. There is a critical coating thickness beyond which thermal conductivity of the coated nanowire is larger than that of the host nanowire. The diameter dependent critical coating thickness and minimum thermal conductivity are explored. Moreover, we found that interface roughness can induce further reduction of thermal conductivity in coated nanowires. From the vibrational eigen-mode analysis, it is found that coating induces localization for low frequency phonons, while interface roughness localizes the high frequency phonons. Our results provide an available approach to tune thermal conductivity of nanowires by atomic layer coating.

preprint2012arXiv

Phononics: Manipulating heat flow with electronic analogs and beyond

The form of energy termed heat that typically derives from lattice vibrations, i.e. the phonons, is usually considered as waste energy and, moreover, deleterious to information processing. However, with this colloquium, we attempt to rebut this common view: By use of tailored models we demonstrate that phonons can be manipulated like electrons and photons can, thus enabling controlled heat transport. Moreover, we explain that phonons can be put to beneficial use to carry and process information. In a first part we present ways to control heat transport and how to process information for physical systems which are driven by a temperature bias. Particularly, we put forward the toolkit of familiar electronic analogs for exercising phononics; i.e. phononic devices which act as thermal diodes, thermal transistors, thermal logic gates and thermal memories, etc.. These concepts are then put to work to transport, control and rectify heat in physical realistic nanosystems by devising practical designs of hybrid nanostructures that permit the operation of functional phononic devices and, as well, report first experimental realizations. Next, we discuss yet richer possibilities to manipulate heat flow by use of time varying thermal bath temperatures or various other external fields. These give rise to a plenty of intriguing phononic nonequilibrium phenomena as for example the directed shuttling of heat, a geometrical phase induced heat pumping, or the phonon Hall effect, that all may find its way into operation with electronic analogs.

preprint2012arXiv

Quantum transport of double quantum dots coupled to an oscillator in arbitrary strong coupling regime

In this paper, we investigate the quantum transport of a double quantum dot coupled with a nanomechanical resonator at arbitrary strong electron-phonon coupling regimes. We employ the generalized quantum master equation to study full counting statistics of currents. We demonstrate the coherent phonon states method can be applied to decouple the electron-phonon interaction non-perturbatively. With the help of this non-perturbative treatment of electron-phonon couplings, we find that the phonon-assisted resonant tunneling emerges when the excess energy from the left quantum dot to the right one can excite integer number of phonons and multi-phonon excitations can enhance the transport in strong electron-phonon coupling regime. Moreover, we find that as the electron-phonon coupling increases, it first plays a constructive role to assist the transport, and then plays the role of scattering and strongly represses the transport.

preprint2012arXiv

Reverse engineering of complex dynamical networks in the presence of time-delayed interactions based on noisy time series

Reverse engineering of complex dynamical networks is important for a variety of fields where uncovering the full topology of unknown networks and estimating parameters characterizing the network structure and dynamical processes are of interest. We consider complex oscillator networks with time-delayed interactions in a noisy environment, and develop an effective method to infer the full topology of the network and evaluate the amount of time delay based solely on noise- contaminated time series. In particular, we develop an analytic theory establishing that the dynamical correlation matrix, which can be constructed purely from time series, can be manipulated to yield both the network topology and the amount of time delay simultaneously. Extensive numerical support is provided to validate the method. While our method provides a viable solution to the network inverse problem, significant difficulties, limitations, and challenges still remain, and these are discussed thoroughly.

preprint2012arXiv

Thermal conductivities of one-dimensional anharmonic/nonlinear lattices: renormalized phonons and effective phonon theory

Heat transport in low-dimensional systems has attracted enormous attention from both theoretical and experimental aspects due to its significance to the perception of fundamental energy transport theory and its potential applications in the emerging filed of phononics: manipulating heat flow with electronic anologs. We consider the heat conduction of one-dimensional nonlinear lattice models. The energy carriers responsible for the heat transport have been identified as the renormalized phonons. Within the framework of renormalized phonons, a phenomenological theory, effective phonon theory, has been developed to explain the heat transport in general one-dimensional nonlinear lattices. With the help of numerical simulations, it has been verified that this effective phonon theory is able to predict the scaling exponents of temperature-dependent thermal conductivities quantitatively and consistently.

preprint2012arXiv

Thermal Contact Resistance Across Nanoscale Silicon Dioxide and Silicon Interface

Silicon dioxide and silicon (SiO$_{2}$/Si) interface plays a very important role in semiconductor industry. However, at nanoscale, its interfacial thermal properties haven't been well understood so far. In this paper, we systematically study the interfacial thermal resistance (Kapitza resistance) of a heterojunction composed of amorphous silicon dioxide and crystalline silicon by using molecular dynamics simulations. Numerical results have shown that Kapitza resistance at SiO$_{2}$/Si interface depends on the interfacial coupling strength remarkably. In the weak interfacial coupling limit, Kapitza resistance depends on both the detailed interfacial structure and the length of the heterojunction, showing large fluctuation among different samples. In contrast, it is almost insensitive to the detailed interfacial structure or the length of the heterojunction in the strong interfacial coupling limit, giving rise to a nearly constant value around 0.9 $\times10^{-9}$ m$^{2}$KW$^{-1}$ at room temperature. Moreover, the temperature dependent Kapitza resistance in the strong interfacial coupling limit has also been examined. Our study provides useful guidance to the thermal management and heat dissipation across nanoscale SiO$_{2}$/Si interface, in particular for the design of silicon nanowire based nano electronics and photonics devices.

preprint2012arXiv

Thermoelectric transport with electron-phonon coupling and electron-electron interaction in molecular junctions

Within the framework of nonequilibrium Green's functions, we investigate the thermoelectric transport in a single molecular junction with electron-phonon and electron-electron interactions. By transforming into a displaced phonon basis, we are able to deal with these interactions non-perturbatively. Then, by invoking the weak tunneling limit, we are able to calculate the thermoelectricity. Results show that at low temperatures, resonances of the thermoelectric figure of merit ZT occur around the sides of resonances of electronic conductance but drops dramatically to zero at exactly these resonant points. We find ZT can be enhanced by increasing electron-phonon coupling and Coulomb repulsion, and an optimal enhancement is obtained when these two interactions are competing. Our results indicate a great potential for single-molecular-junctions as good thermoelectric devices over a wide range of temperatures.

preprint2011arXiv

A nonequilibrium Green's function study of thermoelectric properties in single-walled carbon nanotubes

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

preprint2011arXiv

A Universal Gauge for Thermal Conductivity of Silicon Nanowires With Different Cross Sectional Geometries

By using molecular dynamics simulations, we study thermal conductivity of silicon nanowires (SiNWs) with different cross sectional geometries. It is found that thermal conductivity decreases monotonically with the increase of surface-to-volume ratio (SVR). More interestingly, a simple universal linear dependence of thermal conductivity on SVR is observed for SiNWs with modest cross sectional area (larger than 20 nm^2), regardless of the cross sectional geometry. As a result, among different shaped SiNWs with the same cross sectional area, the one with triangular cross section has the lowest thermal conductivity. Our study provides not only a universal gauge for thermal conductivity among different cross sectional geometries, but also a designing guidance to tune thermal conductivity by geometry.

preprint2011arXiv

Controllability of Complex Networks with Nonlinear Dynamics

Sun et al. provided an insightful comment arXiv:1108.5739v1 on our manuscript entitled "Controllability of Complex Networks with Nonlinear Dynamics" on arXiv. We agree on their main point that linearization about locally desired states can be violated in general by the breakdown of local control of the linearized complex network with nonlinear state. Therefore, we withdraw our manuscript. However, other than nonlinear dynamics, our claim that a single-node-control can fully control the general bidirectional/undirected linear network with 1D self-dynamics is still valid, which is similar to (but different from) the conclusion of arXiv:1106.2573v3 that all-node-control with a single signal can fully control any direct linear network with nodal-dynamics (1D self-dynamics).

preprint2011arXiv

Counting statistics of heat transport in harmonic junctions -- transient and steady states

We study the statistics of heat transferred in a given time interval $t_M$, through a finite harmonic chain, called the center $(C)$, which is connected with two heat baths, the left $(L)$ and the right $(R)$, that are maintained at two different temperatures. The center atoms are driven by an external time-dependent force. We calculate the cumulant generating function (CGF) for the heat transferred out of the left lead, $Q_L$, based on two-time measurement concept and using nonequilibrium Green's function (NEGF) method. The CGF can be concisely expressed in terms of Green's functions of the center and an argument-shifted self-energy of the lead. The expression of CGF is valid in both transient and steady state regimes. We consider three different initial conditions for the density operator and show numerically, for one-dimensional (1D) linear chains, how transient behavior differs from each other but finally approaches the same steady state, independent of the initial distributions. We also derive the CGF for the joint probability distribution $P(Q_L,Q_R)$, and discuss the correlations between $Q_L$ and $Q_R$. We calculate the total entropy flow to the reservoirs. In the steady state we explicitly show that the CGF obeys steady state fluctuation theorem (SSFT). Classical results are obtained by taking $\hbar \to 0$. The method is also applied to the counting of the electron number and electron energy, for which the associated self-energy is obtained from the usual lead self-energy by multiplying a phase or shifting the contour time, respectively.

preprint2011arXiv

Heat generation and transport due to time-dependent forces

We study heat transport for solids in the presence of arbitrary time-dependent force. Using nonequilibrium Green's function (NEGF) approach we present an exact analytical expression of current for the linear system. We found that the heat current can be expressed in terms of the displacement of the atoms in the center and the self energy of the baths. We carry out the calculation for the oscillatory driven force and study the steady state properties for one-dimensional linear chain and two-dimensional square lattice. We found that the heat current is related to the density of states of the system and is independent of the bath temperature in ballistic transport. The baths absorb energy only when their intrinsic frquency resonates with the applied frequency We also generalize the problem for multiple heat baths with different temperatures. We also discuss the effect due to nonlinear interactions in the center.

preprint2011arXiv

Interfacial thermal transport in atomic junctions

We study ballistic interfacial thermal transport across atomic junctions. Exact expressions for phonon transmission coefficients are derived for thermal transport in one-junction and two-junction chains, and verified by numerical calculation based on a nonequilibrium Green's function method. For a single-junction case, we find that the phonon transmission coefficient typically decreases monotonically with increasing freqency. However, in the range between equal frequency spectrum and equal acoustic impedance, it increases first then decreases, which explains why the Kapitza resistance calculated from the acoustic mismatch model is far larger than the experimental values at low temperatures. The junction thermal conductance reaches a maximum when the interfacial coupling equals the harmonic average of the spring constants of the two semi-infinite chains. For three-dimensional junctions, in the weak coupling limit, we find that the conductance is proportional to the square of the interfacial coupling, while for intermediate coupling strength the conductance is approximately proportional to the interfacial coupling strength. For two-junction chains, the transmission coefficient oscillates with the frequency due to interference effects. The oscillations between the two envelop lines can be understood analytically, thus providing guidelines in designing phonon frequency filters.

preprint2011arXiv

Localizations on Complex Networks

We study the structural characteristics of complex networks using the representative eigenvectors of the adjacent matrix. The probability distribution function of the components of the representative eigenvectors are proposed to describe the localization on networks where the Euclidean distance is invalid. Several quantities are used to describe the localization properties of the representative states, such as the participation ratio, the structural entropy, and the probability distribution function of the nearest neighbor level spacings for spectra of complex networks. Whole-cell networks in the real world and the Watts-Strogatz small-world and Barabasi-Albert scale-free networks are considered. The networks have nontrivial localization properties due to the nontrivial topological structures. It is found that the ascending-order-ranked series of the occurrence probabilities at the nodes behave generally multifractally. This characteristic can be used as a structural measure of complex networks.

preprint2011arXiv

Multi-resonance of energy transport and absence of heat pump in a force-driven lattice

Energy transport control in low dimensional nano-scale systems has attracted much attention in recent years. In this paper, we investigate the energy transport properties of Frenkel-Kontorova lattice subject to a periodic driving force, in particular, the resonance behavior of the energy current by varying the external driving frequency. It is discovered that in certain parameter ranges, multiple resonance peaks, instead of a single resonance, emerge. By comparing the nonlinear lattice model with a harmonic chain, we unravel the underlying physical mechanism for such resonance phenomenon. Other parameter dependencies of the resonance behavior are examined as well. Finally, we demonstrate that heat pumping is actually absent in this force-driven model.

preprint2011arXiv

Phonon Coherent Resonance and Its Effect on Thermal Transport In Core-Shell Nanowires

We study heat current autocorrelation function and thermal conductivity in core-shell nanowires by using molecular dynamics simulations. Interestingly, a remarkable oscillation effect in heat current autocorrelation function is observed in core-shell NWs, while the same effect is absent in pure silicon nanowires, nanotube structures and random doped nanowires. Detailed characterizations of the oscillation signal reveal that this intriguing oscillation is caused by the coherent resonance effect of the transverse and longitudinal phonon modes. This phonon resonance results in the localization of the longitudinal modes, which leads to the reduction of thermal conductivity in core-shell nanowires. Our study reveals that a coherent mechanism can be used to tune thermal conductivity in core-shell nanowires.

preprint2011arXiv

Phonon Hall effect in ionic crystals in the presence of static magnetic field

We study phonon Hall effect (PHE) for ionic crystals in the presence of static magnetic field. Using Green-Kubo formula, we present an exact calculation of thermal conductivity tensor by considering both positive and negative frequency phonons. Numerical results are shown for some lattices such as hexagonal lattices, triangular lattices, and square lattices. We find that the PHE occurs on the nonmagnetic ionic crystal NaCl, although the magnitude is very small which is due to the tiny charge-to-mass ratio of the ions. The off-diagonal thermal conductivity is finite for nonzero magnetic field and changes sign for high value of magnetic field at high temperature. We also found that the off-diagonal thermal conductivity diverges as $\pm{1/T}$ at low temperature.

preprint2011arXiv

Spectral Properties of Directed Random Networks with Modular Structure

We study spectra of directed networks with inhibitory and excitatory couplings. We investigate in particular eigenvector localization properties of various model networks for different value of correlation among their entries. Spectra of random networks, with completely uncorrelated entries show a circular distribution with delocalized eigenvectors, where as networks with correlated entries have localized eigenvectors. In order to understand the origin of localization we track the spectra as a function of connection probability and directionality. As connections are made directed, eigenstates start occurring in complex conjugate pairs and the eigenvalue distribution combined with the localization measure shows a rich pattern. Moreover, for a very well distinguished community structure, the whole spectrum is localized except few eigenstates at boundary of the circular distribution. As the network deviates from the community structure there is a sudden change in the localization property for a very small value of deformation from the perfect community structure. We search for this effect for the whole range of correlation strengths and for different community configurations. Furthermore, we investigate spectral properties of a metabolic network of zebrafish, and compare them with those of the model networks.

preprint2011arXiv

The phonon Hall effect: theory and application

We present a systematic theory of the phonon Hall effect in a ballistic crystal lattice system, and apply it on the kagome lattice which is ubiquitous in various real materials. By proposing a proper second quantization for the non-Hermite Hamiltonian in the polarization-vector space, we obtain a new heat current density operator with two separate contributions: the normal velocity responsible for the longitudinal phonon transport, and the anomalous velocity manifesting itself as the Hall effect of transverse phonon transport. As exemplified in kagome lattices, our theory predicts that the direction of Hall conductivity at low magnetic field can be reversed by tuning temperatures, which we hope can be verified by experiments in the future. Three phonon-Hall-conductivity singularities induced by phonon-band-topology change are discovered as well, which correspond to the degeneracies at three different symmetric center points, Γ, K, X, in the wave-vector space of the kagome lattice.

preprint2011arXiv

Topological effect on thermal conductivity in graphene

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

preprint2010arXiv

Ballistic thermal rectification in nanoscale three-terminal junctions

We study ballistic thermal transport in three-terminal atomic nanojunctions by the nonequilibrium Green's function method. We find that there is ballistic thermal rectification in asymmetric three-terminal structures because of the incoherent phonon scattering from the control terminal. With spin-phonon interaction, we also find the ballistic thermal rectification even in symmetric three-terminal paramagnetic structures.

preprint2010arXiv

Berry-Phase induced Heat Pumping and its Impact on the Fluctuation Theorem

Applying adiabatic, cyclic two parameter modulations we investigate quantum heat transfer across an anharmonic molecular junction contacted with two heat baths. We demonstrate that the pumped heat typically exhibits a Berry phase effect in providing an additional geometric contribution to heat flux. Remarkably, a robust fractional quantized geometric phonon response is identified as well. The presence of this geometric phase contribution in turn causes a breakdown of the fluctuation theorem of the Gallavotti-Cohen type for quantum heat transfer. This can be restored only if (i) the geometric phase contribution vanishes and if (ii) the cyclic protocol preserves the detailed balance symmetry.

preprint2010arXiv

Emergence and control of heat current from strict zero thermal bias

It is an ever-growing challenge to develop nano-devices for controlling energy transport. An open question is whether we can create and control heat current at {\it strict zero thermal bias}, and if yes, how to do it. In this paper, we demonstrate that a nonlinear asymmetric system, when pushed out of equilibrium, can produce heat current in the absence of a thermal bias. The emergence and control of heat current over a broad range of parameters are studied. Our results reveal the following three necessary conditions: non-equilibrium source, symmetry breaking, and nonlinearity. We also demonstrate that when heat baths are correlated, symmetry breaking is sufficient to generate heat current.

preprint2010arXiv

Energy carriers in the Fermi-Pasta-Ulam $β$ lattice: Solitons or Phonons?

We investigate anomalous energy transport processes in the Fermi-Pasta-Ulam $β$ lattice. They are determined by the maximum sound velocity of the relevant weakly damped energy carriers. That velocity can be numerically resolved by measuring the propagating fronts of the correlation functions of energy/momentum fluctuations at different times. The numerical results are compared with the predictions for solitons and effective (renormalized) phonons, respectively. Excellent agreement has been found for the prediction of effective long wavelength phonons, giving strong evidence that the energy carriers should be effective phonons rather than solitons.

preprint2010arXiv

How to Improve The Accuracy of Equilibrium Molecular Dynamics For Computation of Thermal Conductivity?

Equilibrium molecular dynamics (EMD) simulations through Green-Kubo formula (GKF) have been widely used in the study of thermal conductivity of various materials. However, there exist controversial simulation results which have huge discrepancies with experimental ones in literatures. In this paper, we demonstrate that the fluctuation in calculated thermal conductivity is due to the uncertainty in determination of the truncation time, which is related to the ensemble and size dependent phonon relaxation time. We thus propose a new scheme in the direct integration of heat current autocorrelation function (HCACF) and a nonzero correction in the double-exponential-fitting of HCACF to describe correctly the contribution to thermal conductivity from low frequency phonons. By using crystalline Silicon (Si) and Germanium (Ge) as examples, we demonstrate that our method can give rise to the values of thermal conductivity in an excellent agreement with experimental ones.

preprint2010arXiv

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

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

preprint2010arXiv

Molecular Dynamics Simulations of Heat Conduction in Nanostructures: Effect of Heat Bath

We investigate systematically the impacts of heat bath used in molecular dynamics simulations on heat conduction in nanostructures exemplified by Silicon Nanowires (SiNWs) and Silicon/Germanium nano junction. It is found that multiple layers of Nosé-Hoover heat bath are required to reduce the temperature jump at the boundary, while only a single layer of Langevin heat bath is sufficient to generate a linear temperature profile with small boundary temperature jump. Moreover, an intermediate value of heat bath parameter is recommended for both Nosé-Hoover and Langevin heat bath in order to achieve correct temperature profile and thermal conductivity in homogeneous materials. Furthermore, the thermal rectification ratio in Si/Ge thermal diode depends on the choice of Nosé-Hoover heat bath parameter remarkably, which may lead to non-physical results. In contrast, Langevin heat bath is recommended because it can produce consistent results with experiment in large heat bath parameter range.

preprint2010arXiv

Noise bridges dynamical correlation and topology in coupled oscillator networks

We study the relationship between dynamical properties and interaction patterns in complex oscillator networks in the presence of noise. A striking finding is that noise leads to a general, one-to-one correspondence between the dynamical correlation and the connections among oscillators for a variety of node dynamics and network structures. The universal finding enables an accurate prediction of the full network topology based solely on measuring the dynamical correlation. The power of the method for network inference is demonstrated by the high success rate in identifying links for distinct dynamics on both model and real-life networks. The method can have potential applications in various fields due to its generality, high accuracy and efficiency.

preprint2010arXiv

Phonon Transport in Suspended Single Layer Graphene

We report the first temperature dependent phonon transport measurements in suspended Cu-CVD single layer graphene (SLG) from 15K to 380K using microfabricated suspended devices. The thermal conductance per unit cross section $σ$/A increases with temperature and exhibits a peak near T~280K ($\pm$10K) due to the Umklapp process. At low temperatures (T<140K), the temperature dependent thermal conductivity scales as ~T^{1.5}, suggesting that the main contribution to thermal conductance arises from flexural acoustic (ZA) phonons in suspended SLG. The $σ$/A reaches a high value of 1.7$\times10^5 T^{1.5}$ W/m^2K, which is approaching the expected ballistic phonon thermal conductance for two-dimensional graphene sheets. Our results not only clarify the ambiguity in the thermal conductance, but also demonstrate the potential of Cu-CVD graphene for heat related applications.

preprint2010arXiv

Random matrix analysis of localization properties of Gene co-expression network

We analyze gene co-expression network under the random matrix theory framework. The nearest neighbor spacing distribution of the adjacency matrix of this network follows Gaussian orthogonal statistics of random matrix theory (RMT). Spectral rigidity test follows random matrix prediction for a certain range, and deviates after wards. Eigenvector analysis of the network using inverse participation ratio (IPR) suggests that the statistics of bulk of the eigenvalues of network is consistent with those of the real symmetric random matrix, whereas few eigenvalues are localized. Based on these IPR calculations, we can divide eigenvalues in three sets; (A) The non-degenerate part that follows RMT. (B) The non-degenerate part, at both ends and at intermediate eigenvalues, which deviate from RMT and expected to contain information about {\it important nodes} in the network. (C) The degenerate part with $zero$ eigenvalue, which fluctuates around RMT predicted value. We identify nodes corresponding to the dominant modes of the corresponding eigenvectors and analyze their structural properties.

preprint2010arXiv

Remarkable Reduction of Thermal Conductivity in Silicon Nanotubes

We propose to reduce the thermal conductivity of silicon nanowires (SiNWs) by introducing small hole at the centre, i.e. construct silicon nanotube (SiNT) structures. Our numerical results demonstrate that a very small hole (only 1% reduction in cross section area) can induce a 35% reduction in room temperature thermal conductivity. Moreover, with the same cross section area, thermal conductivity of SiNT is only about 33% of that of SiNW at room temperature. The spatial distribution of vibrational energy reveals that localization modes are concentrated on the inner and outer surfaces of SiNTs. The enhanced surface-to-volume ratio in SiNTs reduces the percentage of delocalized modes, which is believed to be responsible for the reduction of thermal conductivity. Our study suggests SiNT is a promising thermoelectric material with low thermal conductivity.

preprint2010arXiv

Thermal contraction in silicon nanowires at low temperatures

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

preprint2010arXiv

Topological Nature of the Phonon Hall Effect

We provide a topological understanding on phonon Hall effect in dielectrics with Raman spinphonon coupling. A general expression for phonon Hall conductivity is obtained in terms of the Berry curvature of band structures. We find a nonmonotonic behavior of phonon Hall conductivity as a function of magnetic field. Moreover, we observe a phase transition in phonon Hall effect, which corresponds to the sudden change of band topology, characterized by the altering of integer Chern numbers. This can be explained by touching and splitting of phonon bands.

preprint2010arXiv

Violation of Fourier's Law and Anomalous Heat Diffusion in Silicon

We study heat conduction and diffusion in silicon nanowires (SiNWs) systematically by using non-equilibrium molecular dynamics. It is found that the thermal conductivity of SiNWs diverges with the length, even when the length is up to 1,100 nm which is much longer than the phonon mean free path. Moreover, an anomalous heat diffusion is observed which is believed to be responsible for the length dependent thermal conductivity. Our results provide strong evidence that Fourier's law of heat conduction is not valid in low dimensional nanostructures.

preprint2009arXiv

Size dependent thermoelectric properties of silicon nanowires

By using first-principles tight-binding electronic structure calculation and Boltzmann transport equation, we investigate the size dependence of thermoelectric properties of silicon nanowires (SiNWs). With cross section area increasing, the electrical conductivity increases slowly, while the Seebeck coefficient reduces remarkably. This leads to a quick reduction of cooling power factor with diameter. Moreover, the figure of merit also decreases with transverse size. Our results demonstrate that in thermoelectric application, NW with small diameter is preferred. We also predict that isotopic doping can increase the value of ZT significantly. With 50% 29Si doping (28Si0.529Si0.5 NW), the ZT can be increased by 31%.

preprint2009arXiv

Thermal conductivity of one-dimensional lattices with self-consistent heat baths: a heuristic derivation

We derive the thermal conductivities of one-dimensional harmonic and anharmonic lattices with self-consistent heat baths (BRV lattice) from the Single-Mode Relaxation Time (SMRT) approximation. For harmonic lattice, we obtain the same result as previous works. However, our approach is heuristic and reveals phonon picture explicitly within the heat transport process. The results for harmonic and anharmonic lattices are compared with numerical calculations from Green-Kubo formula. The consistency between derivation and simulation strongly supports that effective (renormalized) phonons are energy carriers in anharmonic lattices although there exist some other excitations such as solitons and breathers.

preprint2007arXiv

Stability of quantum motion in regular systems: a uniform semiclassical approach

We study the stability of quantum motion of classically regular systems in presence of small perturbations. Onthe base of a uniform semiclassical theory we derive the fidelity decay which displays a quite complexbehaviour, from Gaussian to power law decay $t^{-α}$ with $1 \le α\le 2$. Semiclassical estimates are given for the time scales separating the different decaying regions and numerical results are presented which confirm our theoretical predictions.

preprint2007arXiv

Ultralow Thermal Conductivity of Isotope-Doped Silicon Nanowires

The thermal conductivity of silicon nanowires (SiNWs) is investigated by molecular dynamics (MD) simulation. It is found that the thermal conductivity of SiNWs can be reduced exponentially by isotopic defects at room temperature. The thermal conductivity reaches the minimum, which is about 27% of that of pure 28Si NW, when doped with fifty percent isotope atoms. The thermal conductivity of isotopic-superlattice structured SiNWs depends clearly on the period of superlattice. At a critical period of 1.09 nm, the thermal conductivity is only 25% of the value of pure Si NW. An anomalous enhancement of thermal conductivity is observed when the superlattice period is smaller than this critical length. The ultra-low thermal conductivity of superlattice structured SiNWs is explained with phonon spectrum theory.

preprint2006arXiv

Current behavior of a quantum Hamiltonian ratchet in resonance

We investigate the ratchet current that appears in a kicked Hamiltonian system when the period of the kicks corresponds to the regime of quantum resonance. In the classical analogue, a spatial-temporal symmetry should be broken to obtain a net directed current. It was recently discovered that in quantum resonance the temporal symmetry can be kept, and we prove that breaking the spatial symmetry is a necessary condition to find this effect. Moreover, we show numerically and analytically how the direction of the motion is dramatically influenced by the strength of the kicking potential and the value of the period. By increasing the strength of the interaction this direction changes periodically, providing us with a non-expected source of current reversals in this quantum model. These reversals depend on the kicking period also, though this behavior is theoretically more difficult to analyze. Finally, we generalize the discussion to the case of a non-uniform initial condition.

preprint2006arXiv

Negative differential thermal resistance and thermal transistor

We report on the first model of a thermal transistor to control heat flow. Like its electronic counterpart, our thermal transistor is a three-terminal device with the important feature that the current through the two terminals can be controlled by small changes in the temperature or in the current through the third terminal. This control feature allows us to switch the device between "off" (insulating) and "on" (conducting) states or to amplify a small current. The thermal transistor model is possible because of the negative differential thermal resistance.

preprint2005arXiv

Modification to the pre-factor of the semiclassical propagator

We modify the pre-factor of the semiclassical propagator to improve its efficiency in practical implementations. The new pre-factor represents the smooth portion of an orbit's contribution, and leads to fast convergence in numerical calculations. As an illustration of the accuracy and efficiency of the resultant propagator, we numerically calculate overlaps between quantum and semiclassical wave functions, as well as low-lying spectrum density in a 10-dimensional system contains unstable classical orbits. This sheds light on applying semiclassical propagator to high dimensional systems.

preprint2005arXiv

Uniform semiclassical approach to fidelity decay in the deep Lyapunov regime

We use the uniform semiclassical approximation in order to derive the fidelity decay in the regime of large perturbations. Numerical computations are presented which agree with our theoretical predictions. Moreover, our theory allows to explain previous findings, such as the deviation from the Lyapunov decay rate in cases where the classical finite-time instability is non-uniform in phase space.

preprint2003arXiv

Crossover of quantum Loschmidt echo from golden rule decay to perturbation-independent decay

We study the crossover of the quantum Loschmidt echo (or fidelity) from the golden rule regime to the perturbation-independent exponential decay regime by using the kicked top model. It is shown that the deviation of the perturbation-independent decay of the averaged fidelity from the Lyapunov decay results from quantum fluctuations in individual fidelity, which are caused by the coherence in the initial coherent states. With an averaging procedure suppressing the quantum fluctuations effectively, the perturbation-independent decay is found to be close to the Lyapunov decay. We also show that the Fourier transform of the fidelity is determined directly by the initial state and the eigenstates of the Floquet operators of the two classically chaotic systems concerned. The absolute value part and the phase part of the Fourier transform of the fidelity are found to be divided into several correlated parts, which is a manifestation of the coherence of the initial coherent state. In the whole crossover region, some important properties of the fidelity, such as the exponent of its exponential decay and the short initial time within which the fidelity almost does not change, are found to be closely related to properties of the central part of its Fourier transform.

preprint1999arXiv

Quantum chaos of a kicked particle in a 1D infinite square potential well

We study quantum chaos in a non-KAM system, i.e. a kicked particle in a one-dimensional infinite square potential well. Within the perturbative regime the classical phase space displays stochastic web structures, and the diffusion coefficient D in the regime increases with the perturbative strength K giving a scaling $D \propto K^{2.5}$, and in the large K regime D goes as K^2. Quantum mechanically, we observe that the level spacing statistics of the quasi eigenenergies changes from Poisson to Wigner distribution as the kick strength increases. The quasi eigenstates show power-law localization in the small K region, which become extended one at large K. Possible experimental realization of this model is also discussed.