BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jian-Sheng Wang

Jian-Sheng Wang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mes-hallcond-mat.stat-mechcond-mat.mtrl-sciquant-phphysics.comp-phcond-mathep-latphysics.gen-ph

Catalog footprint

What is connected

84works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

84 published item(s)

preprint2022arXiv

Anharmonic quantum thermal transport across a van der Waals interface

We investigate the anharmonic phonon scattering across a weakly interacting interface by developing a quantum mechanics-based theory. We find that the contribution from anharmonic three-phonon scatterings to interfacial thermal conductance can be cast into Landauer formula with transmission function being temperature-dependent. Surprisingly, in the weak coupling limit, the transmission due to anharmonic phonon scattering is unbounded with increasing temperature, which is physically impossible for two-phonon processes. We further reveal that the anharmonic contribution in a real heterogeneous interface (e.g., between graphene and monolayer molybdenum disulfide) can dominate over the harmonic process even at room temperature, highlighting the important role of anharmonicity in weakly interacting heterogeneous systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Microscopic theory of photon-induced energy, momentum, and angular momentum transport in the nonequilibrium regime

We set up a general microscopic theory for the transfer of energy, momentum, and angular momentum mediated by photons. Using the nonequilibrium Green's function method, we propose a unified Meir-Wingreen formalism for the energy emitted, force experienced, and torque experienced by the objects due to the fluctuating electromagnetic field. Our theory does not require the local thermal equilibrium that is the central assumption of the conventional theory of fluctuational electrodynamics (FE). The obtained formulas are valid for arbitrary objects as well as the environment without the requirement of reciprocity. To show the capability of our microscopic theory, we apply the general formulas to transport problems of graphene edges in both equilibrium and nonequilibrium situations. We show the local equilibrium energy radiation of graphene obeys the well-known $T^4$ law with a converged theoretical emissivity of 2.058$\%$. In the ballistic nonequilibrium situation driven by chemical potential biases, we observe nonzero results for force and torque from the graphene edges, which go beyond the predictive ability of the FE theory. Our method is general and efficient for large systems, which paves the way for studying more complex transport phenomena in the nonequilibrium regime.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Angular momentum radiation from current-carrying molecular junctions

We consider the radiation of angular momentum (AM) from current-carrying molecular junctions. Using the nonequilibrium Green's function method, we derive a convenient formula for the AM radiation and apply it to a prototypical benzene molecule junction. We discuss the selection rules for inelastic transitions between the molecular angular momentum eigenstates due to a 6-fold rotational symmetry. Our study provides important insights into the generation of light with AM from DC-biased molecular junctions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Caroli formula in near-field heat transfer between parallel graphene sheets

In this work we conduct a close-up investigation into the nature of near-field heat transfer (NFHT) of two graphene sheets in parallel-plate geometry. We develop a fully microscopic and quantum approach using nonequilibrium Green's function method. A Caroli formula for heat flux is proposed and numerically verified. We show our near-field-to-black-body heat flux ratios generally exhibit $1/d^α$ dependence, with an effective exponent $α\approx 2.2$, at long distances exceeding 100 nm and up to one micron; in the opposite $d\rightarrow 0$ limit, the values converge to a range within an order of magnitude. We justify this feature by noting it is owing to the breakdown of local conductivity theory, which predicts a $1/d$ dependence. Furthermore, from the numerical result, we find in addition to thermal wavelength, $λ_{th}$, a shorter distance scale $\sim$ 10 - 100 nm, comparable to the graphene thermal length ($\hbar v_{F}/k_{B} T$) or Fermi wavelength ($k_{F}^{-1}$), marks the transition point between the short- and long-distance transfer behaviors; within that point, relatively large variation of heat flux in response to doping level becomes a typical characteristic. The emergence of such large variation is tied to relative NFHT contributions from the intra- and inter-band transitions. Beyond that point, scaling of thermal flux $\propto 1/d^α$ can be generally observed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

First-principles method to study near-field radiative heat transfer

We present a general and convenient first principle method to study near-field radiative heat transfer. We show that the Landauer-like expression of heat flux can be expressed in terms of a frequency and wave-vector dependent macroscopic dielectric function which can be obtained from the linear response density functional theory. A random phase approximation is used to calculate the response function. We computed the heat transfer in three systems -- graphene, molybdenum disulfide (MoS$_2$), and hexagonal boron nitride (h-BN). Our results show that the near-field heat flux exceeds the blackbody limit up to four orders of magnitude. With the increase of the distances between two parallel sheets, a $1/d^2$ dependence of heat flux is shown, consistent with Coulomb's law. The heat transfer capacity is sensitive to the dielectric properties of materials. Influences from chemical potential and temperature are also discussed. Our method can be applied to a wide range of materials including systems with inhomogeneities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ab initio dipolar electron-phonon interactions in two-dimensional materials

We develop an ab initio formalism for dipolar electron-phonon interactions (EPI) in two-dimensional (2D) materials. Unlike purely longitudinal Fröhlich model, we show that the out-of-plane dipoles also contribute to the long-wavelength non-analytical behavior of EPI. And the 2D dipolar EPI plays an important role not only in the typical polar material MoS$_2$, but also in graphane and fluorinated graphene. By incorporating this formalism into Wannier-Fourier interpolation, we enable accurate EPI calculations for 2D materials and subsequent intrinsic carrier mobility prediction. The results show that Fröhlich model is inadequate for 2D materials and correct long-wavelength interaction must be included for the reliable prediction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Geometric effect on near-field heat transfer analysis using efficient graphene and nanotube models

Following the recent research enthusiasm on the effect of geometry on near-field heat transfer (NFHT) enhancement, we present an analysis based on simplified yet highly efficient graphene and nanotube models. Two geometries are considered: that of two parallel infinite "graphene" surfaces and that of a one-dimensional infinite "nanotube" line in parallel with an infinite surface. Due to its symmetry, the former is in principal simpler to analyze and even so, earlier works suggested that the application of a full model in this problem still demands heavy computations. Among other findings, our simplified computation - having successfully replicated the results of relevant earlier works - suggests a sharper NFHT enhancement dependence on distance for the line-surface system, namely $J\sim d^{-5.1}$ as compared to $J\sim d^{-2.2}$ for the parallel surface. Such comparisons together with applications of our efficient approach would be the important first steps in the attempt to find a general rule describing geometric dependence of NFHT.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nonadiabatic particle and energy pump at strong system-reservoir coupling

We study the dynamics of electron and energy currents in a nonadiabatic pump. The pump is a quantum dot nanojunction with time-varying gate potential and tunnel couplings to the leads. The leads are unbiased and maintained at the same temperature and chemical potential. We find that synchronized variations of the gate and tunnel couplings can pump electrons and energy from the left to the right lead. Inspired by quantum heat engines, we devise a four-stroke operating protocol that can optimally pump energy and hence, we investigate energy transfer and the coefficient of performance of the device. We compare our device to a two-stroke pump and find that the latter's lower performance is due to the bi-directional flow of energy currents resulting in low net energy currents. The performance of our four-stroke pump can be improved, up to a point, by increasing the net energy carried by the pumped electrons through energy charging via the gate potential. This is achieved by increasing the durations of energy charging and discharging strokes in the pump's protocol. However, despite the large energy output for long charging and discharging strokes, the energy required to maintain the strokes become large too resulting in a stagnant pump performance. Our pump operates effectively only in the strong lead coupling regime and becomes a dud in the weak coupling regime due to the net output energy flowing in the reverse direction. We use nonequilibirum Green's functions techniques to calculate the currents and capture the effects of strong lead-channel coupling exactly while simultaneously incorporating three time-varying parameters. Results from our work could aid in the design of high-performance quantum pumps.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermoelectric Penta-Silicene with a High Room-Temperature Figure of Merit

Silicon is one of the most frequently used chemical elements of the periodic table in nanotechnology. Two-dimensional (2D) silicene, a silicon analog of graphene, has been readily obtained to make field-effect transistors since 2015. Recently, as new members of the silicene family, penta-silicene and its nanoribbon have been experimentally grown on Ag(110) surface with exotic electronic properties. However, the thermoelectric performance of penta-silicene has not been so far studied that would hinder its potential applications of electric generation from waste heat and solid-state Peltier coolers. Based on the Boltzmann transport theory and ab initio calculations, we find that penta-silicene shows a remarkable room temperature figure of merit ZT of 3.4 and 3.0 at the reachable hole and electron concentrations, respectively. We attribute this high ZT to the superior pudding-mold electronic band structure and ultralow lattice thermal conductivity. The discovery provides new insight into the transport property of pentagonal nanostructures and highlights the potential applications of thermoelectric materials at room temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Three-terminal interface as a thermoelectric generator beyond Seebeck effect

We investigate thermoelectric transport through interfaces with inelastic scatterings by developing a quantum theory, which has been extensively validated by existing theories. We find that under temperature bias, while a two-terminal conductor-insulator interface behaves only as a thermal resistor, a three-terminal conductor-insulator-conductor interface can function as an electricity generator caused by phonon-mediated electron scatterings with heat-charge current separation. Unlike conventional thermoelectrics which is a property of a bulk caused by the Seebeck effect, this thermoelectric behavior is a property of an interface driven by electron-phonon scatterings.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A microscopic theory for ultra-near-field radiation

Using the nonequilibrium Green's function (NEGF) formalism, we propose a microscopic theory for near-field radiative heat transfer between metal plates. Tight-binding models for the electrons are coupled to the electromagnetic field continuum. Our approach differs from the established ones based on fluctuational electrodynamics, in that it describes truly nonequilibrium steady states, and is nonlocal in system's dielectric properties. For a two quantum-dot model a new length scale emerges at which the heat current shows a peak. This length scale is related to the physics of parallel plate capacitors. The three-dimensional model results are consistent with the theory of Polder and van Hove except at very short distances.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electron and phonon drag in thermoelectric transport through coherent molecular conductors

We study thermoelectric transport through a coherent molecular conductor connected to two electron and two phonon baths using the nonequilibrium Green's function method. We focus on the mutual drag between electron and phonon transport as a result of `momentum' transfer, which happens only when there are at least two phonon degrees of freedom. After deriving expressions for the linear drag coefficients, obeying the Onsager relation, we further investigate their effect on nonequilibrium transport. We show that the drag effect is closely related to two other phenomena: (1) adiabatic charge pumping through a coherent conductor; (2) the current-induced nonconservative and effective magnetic forces on phonons.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phonon transport in a one-dimensional harmonic chain with long-range interaction and mass disorder

Atomic mass and interatomic interaction are the two key quantities that significantly affect the heat conduction carried by phonons. Here, we study the effects of long-range (LR) interatomic interaction and mass-disorder on the phonon transport in a one-dimensional harmonic chain withup to 10^5 atoms. We find that while LR interaction reduces the transmission of low frequency phonons, it enhances the transmission of high frequency phonons by suppressing the localization effects caused by mass disorder. Therefore, long-range interaction is able to boost heat conductance in the high temperature regime or in the large size regime, where the high frequency modes are important.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Critique of the Gibbs volume entropy and its implication

Dunkel and Hilbert, "Consistent thermostatistics forbids negative absolute temperatures," Nature Physics, {\bf 10}, 67 (2014), and Hilbert, Hänggi, and Dunkel, "Thermodynamic laws in isolated systems," Phys. Rev. E {\bf 90}, 062116 (2014) have presented an unusual view of thermodynamics that sets aside several properties that have traditionally have been assumed to be true. Among other features, their results do not satisfy the postulates of thermodynamics originally proposed by Tisza and Callen. In their theory, differences in the temperatures of two objects cannot determine the direction of heat flow when they are brought into thermal contact. They deny that negative temperatures are possible. We disagree with most of their assertions, and we present arguments in favor of a more traditional interpretation of thermodynamics. We show explicitly that it is possible to deduce the thermodynamic entropy for a paramagnet along the lines of classical thermodynamics and that it agrees with Boltzmann entropy in equilibrium. A Carnot engine with efficiency larger than one is a matter of definition and is possible for inverted energy distributions, regardless of whether negative temperatures are used in the analysis. We elaborate on Penrose's argument that an adiabatic and reversible process connecting systems with Hamiltonian $H$ to $-H$ is possible, thus negative temperatures logically must exist. We give a demonstration that the Boltzmann temperature determines the direction of heat flow while Gibbs temperature does not, for sufficiently large systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Current-induced forces and hot-spots in biased nano-junctions

We investigate theoretically the interplay of current-induced forces (CIF), Joule heating, and heat transport inside a current-carrying nano-conductor. We find that the CIF, due to the electron-phonon coherence, can control the spatial heat dissipation in the conductor. This yields a significant asymmetric concentration of excess heating (hot-spot) even for a symmetric conductor. When coupled to the electrode phonons, CIF drive different phonon heat flux into the two electrodes. First-principles calculations on realistic biased nano-junctions illustrate the importance of the effect.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large adiabatic temperature and magnetic entropy changes in EuTiO3

We have investigated the magnetocaloric effect in single and polycrystalline samples of quantum paraelectric EuTiO3 by magnetization and heat capacity measurements. Single crystalline EuTiO3 shows antiferromagnetic ordering due to Eu2+ magnetic moments below TN = 5.6 K. This compound shows a giant magnetocaloric effect around its Neel temperature. The isothermal magnetic entropy change is 49 Jkg-1K-1, the adiabatic temperature change is 21 K and the refrigeration capacity is 500 JKg-1 for a field change of 7 T at TN. The single crystal and polycrystalline samples show similar values of the magnetic entropy change and adiabatic temperature changes. The large magnetocaloric effect is due to suppression of the spin entropy associated with localized 4f moment of Eu2+ ions. The giant magnetocaloric effect together with negligible hysteresis, suggest that EuTiO3 could be a potential material for magnetic refrigeration below 20 K.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Negative temperatures and the definition of entropy

The concept of negative temperature has recently received renewed interest in the context of debates about the correct definition of the thermodynamic entropy in statistical mechanics. Several researchers have identified the thermodynamic entropy exclusively with the "volume entropy'' suggested by Gibbs, and have further concluded that by this definition, negative temperatures violate the principles of thermodynamics. We disagree with these conclusions. We demonstrate that volume entropy is inconsistent with the postulates of thermodynamics for systems with non-monotonic energy densities, while a definition of entropy based on the probability distributions of macroscopic variables does satisfy the postulates of thermodynamics. Our results confirm that negative temperature is a valid extension of thermodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Gibbs "volume" entropy is incorrect

In recent papers, several authors have claimed that a definition of the thermodynamic entropy in terms of the logarithm of a volume in phase space, originally suggested by Gibbs, is the only valid definition. Arguing from the Gibbs entropy, these authors claim that thermodynamics cannot be extended to include negative temperatures. However, the Gibbs entropy fails to satisfy the postulates of thermodynamics, leading to serious errors. In particular, predictions of the Gibbs entropy for systems with non-monotonic energy densities are incorrect. We show that the correct expression for the equilibrium entropy contains an integral over a surface in phase space, and negative temperature is a valid thermodynamic concept.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improved Dyson series expansion for steady-state quantum transport beyond the weak coupling limit - divergences and resolution

We present a general theory to calculate the steady-state heat and electronic currents for nonlinear systems using a perturbative expansion in the system-bath coupling. We explicitly demonstrate that using the truncated Dyson-series leads to divergences in the steady-state limit, thus making it impossible to be used for actual applications. In order to resolve the divergences we propose a unique choice of initial-condition for the reduced density matrix, which removes the divergences at each order. Our approach not only allows us to use the truncated Dyson-series, with a reasonable choice of initial condition, but also gives the expected result that the steady-state solutions should be independent of initial preparations. Using our improved Dyson series we evaluate the heat and electronic currents upto fourth-order in system-bath coupling, a considerable improvement over the standard quantum master equation techniques. We then numerically corroborate our theory for archetypal settings of linear systems using the exact nonequilibrium Green's function approach. Lastly, to demonstrate the advantage of our approach we deal with the nonlinear spin-boson model to evaluate heat current upto fourth-order and find signatures of cotunnelling process.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Giant magnetocaloric effect in magnetoelectric Eu1-xBaxTiO3

We report the magnetic entropy change (Delta Sm) in magnetoelectric Eu1-xBaxTiO3 for x = 0.1- 0.9. We find - delta Sm = 11 (40) J/kg.K in x = 0.1 for a field change of 1 (5) Tesla respectively, which is the largest value among all Eu-based oxides. Delta Sm arises from the field-induced suppression of the spin entropy of Eu2+:4f7 localized moments. While -delta Sm decreases with increasing x, -DeltaSm = 6.58 J/kg.K observed in the high spin diluted composition x = 0.9 is larger than that in many manganites. Our results indicate that these magnetoelectrics are potential candidates for cryogenic magnetic refrigeration.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonequilibrium Green's function method for quantum thermal transport

This review deals with the nonequilibrium Green's function (NEGF) method applied to the problems of energy transport due to atomic vibrations (phonons), primarily for small junction systems. We present a pedagogical introduction to the subject, deriving some of the well-known results such as the Laudauer-like formula for heat current in ballistic systems. The main aim of the review is to build the machinery of the method so that it can be applied to other situations, which are not directly treated here. In addition to the above, we consider a number of applications of NEGF, not in routine model system calculations, but in a few new aspects showing the power and usefulness of the formalism. In particular, we discuss the problems of multiple leads, coupled left-right-lead system, and system without a center. We also apply the method to the problem of full counting statistics. In the case of nonlinear systems, we make general comments on the thermal expansion effect, phonon relaxation time, and a certain class of mean-field approximations. Lastly, we examine the relationship between NEGF, reduced density matrix, and master equation approaches to thermal transport.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reduced density matrix for nonequilibrium steady states: A modified Redfield solution approach

We describe a method to obtain the reduced density matrix (RDM) correct up to second order in system-bath coupling in \emph{nonequilibrium} steady state situations. The RDM is obtained via a scheme based on analytic continuation, using the time-local Redfield-like quantum master equation, which was earlier used by the same authors [J. Chem. Phys. \textbf{136}, 194110 (2012)] to obtain the correct thermal equilibrium description. This nonequilibrium modified Redfield solution is then corroborated with the exact RDM obtained via the nonequilibrium Green's function technique for the quantum harmonic oscillator. Lastly, the scheme is compared to different quantum master equations (QMEs), namely the time-local Redfield-like and the Lindblad-like QMEs, in order to illustrate the differences between each of these approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spin rectification in thermally driven XXZ spin chain via the spin-Seebeck effect

We study the phenomenon of spin-current rectification in a one-dimensional XXZ spin chain in the presence of a thermal drive. In our model a pure spin current is generated by a temperature difference between two harmonic heat baths which in turn creates a spin voltage via the spin-Seebeck effect. Along with a local spin-current operator definition and the nonequilibrium modified Redfield solution we study the spin-rectification ratio as a function of system size and external magnetic field. Intriguing effects are observed at low temperatures such as oscillations with system size and high range of tunability with external magnetic field making magnetic insulators, which are well described by the XXZ model, ideal candidates to build spin devices based on rectification.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Work distribution under continuous quantum histories

We present a general scheme to obtain work distribution in closed systems under continuous quantum histories of corresponding "power" operator. The scheme is tested by analytically calculating the quantum work distribution for a prototype model of a center-shifted one-dimensional quantum harmonic oscillator. Furthermore, its relationship to the path integral formalism with the Wick rotation is immediately noticed. We find that, based on the continuous measurement, generally Jarzynski equality in quantum regime is not valid, though it is recovered in high-temperature limit. Also we briefly explore the extension to open systems, in which the work done and heat transfer can be calculated in the same set of quantum histories. Our results for this model can be directly verified by single-molecule force spectroscopy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A generalized Caroli formula for transmission coefficient with lead-lead coupling

We present a generalized transmission coefficient formula for the lead-junction-lead system, in which interaction between the leads has been taken into account. Based on it the Caroli formula could be easily recovered and a transmission coefficient formula for interface problem in the ballistic system can be obtained. The condition of validity for the formula is carefully explored. We mainly focus on heat transport. However, the corresponding electrical transport could be similarly dealt with. Also, an illustrative example is given to clarify the precise meaning of the quantities used in the formula, such as the concept of the reduced interacting matrix in different situations. In addition, an explicit transmission coefficient formula for a general one-dimensional interface setup is obtained based on the derived interface formula.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cumulant generating function formula of heat transfer in ballistic system with lead-lead coupling

Based on two-time observation protocol, we consider heat transfer in a given time interval $t_{M}$ in lead-junction-lead system taking coupling between the leads into account. In view of the two-time observation, consistency conditions are carefully verified in our specific family of quantum histories. Furthermore, its implication is briefly explored. Then using nonequilibrium Green's function (NEGF) method, we obtain an exact formula for the cumulant generating function for heat transfer between the two leads valid in both transient and steady-state regimes. Also, a compact formula for the cumulant generating function in the long-time limit is derived, for which the Gallavotti-Cohen (GC) fluctuation symmetry is explicitly verified. In addition, a brief discussion of Di Ventra's trick regarding whether the effect of the repartitioning procedure of the total Hamiltonian on nonequilibium steady-state current fluctuation exists is given. All kinds of properties of nonequilibrium current fluctuations such as the fluctuation theorem in different time regimes could be readily given according to these exact formulas.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dimensional crossover of thermal conductance in graphene nanoribbons: A first-principles approach

First-principles density-functional calculations are performed to investigate the thermal transport properties in graphene nanoribbons (GNRs). The dimensional crossover of thermal conductance from one to two dimensions (2D) is clearly demonstrated with increasing ribbon width. The thermal conductance of GNRs in a few nanometer width already exhibits an approximate low-temperature dependence of $T^{1.5}$, like that of 2D graphene sheet which is attributed to the quadratic nature of dispersion relation for the out-of-plane acoustic phonon modes. Using a zone-folding method, we heuristically derive the dimensional crossover of thermal conductance with the increase of ribbon width. Combining our calculations with the experimental phonon mean-free path, some typical values of thermal conductivity at room temperature are estimated for GNRs and for 2D graphene sheet, respectively. Our findings clarify the issue of low-temperature dependence of thermal transport in GNRs and suggest a calibration range of thermal conductivity for experimental measurements in graphene-based materials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Generalized Gibbs state with modified Redfield solution: Exact agreement up to second order

A novel scheme for the steady state solution of the standard Redfield quantum master equation is developed which yields agreement with the exact result for the corresponding reduced density matrix up to second order in the system-bath coupling strength. We achieve this objective by use of an analytic continuation of the off-diagonal matrix elements of the Redfield solution towards its diagonal limit. Notably, our scheme does not require the provision of yet higher order relaxation tensors. Testing this modified method for a heat bath consisting of a collection of harmonic oscillators we assess that the system relaxes towards its correct coupling-dependent, generalized quantum Gibbs state in second order. We numerically compare our formulation for a damped quantum harmonic system with the nonequilibrium Green's function formalism: we find good agreement at low temperatures for coupling strengths that are even larger than expected from the very regime of validity of the second-order Redfield quantum master equation. Yet another advantage of our method is that it markedly reduces the numerical complexity of the problem; thus allowing to study efficiently large-sized \emph{system} Hilbert spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Role of the on-site pinning potential in establishing quasi-steady-state conditions of heat transport in finite quantum systems

We study the transport of energy in a finite linear harmonic chain by solving the Heisenberg equation of motion, as well as by using nonequilibrium Green's functions to verify our results. The initial state of the system consists of two separate and finite linear chains that are in their respective equilibriums at different temperatures. The chains are then abruptly attached to form a composite chain. The time evolution of the current from just after switch-on to the transient regime and then to later times is determined numerically. We expect the current to approach a steady-state value at later times. Surprisingly, this is possible only if a nonzero quadratic on-site pinning potential is applied to each particle in the chain. If there is no on-site potential a recurrent phenomenon appears when the time scale is longer than the traveling time of sound to make a round trip from the midpoint to a chain edge and then back. Analytic expressions for the transient and steady-state currents are derived to further elucidate the role of the on-site potential.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Steady state thermal transport in anharmonic systems: Application to molecular junctions

We develop a general theory for thermal transport in anharmonic systems under the weak system-bath coupling approximation similar to the quantum master equation formalism. A current operator is derived, which is valid not only in the steady state, but in the transient regime as well. Here we focus on the effects of anharmonicity on the steady-state thermal conductance of a mono and diatomic molecular junctions. We also study molecules being confined in a double-well potential. We find that when the molecules have a non-linear on-site potential the low-temperature thermal conductance is dramatically affected by the strength of non-linearity, whereas for the diatomic molecule connected by an anharmonic spring the strength of anharmonicity plays almost no role in the low-temperature regime. In case of the molecules confined in a double-well potential we find that the height of the barrier greatly affects the thermal conductance; once the molecules can feel the effect of the barrier we observe negative differential thermal conductance at both high and low temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tuning thermal transport in nanotubes with topological defects

Using the atomistic nonequilibrium Green's function, we find that thermal conductance of carbon nanotubes with presence of topological lattice imperfects is remarkably reduced, due to the strong Rayleigh scattering of high-frequency phonons. Phonon transmission across multiple defects behaves as a cascade scattering based with the random phase approximation. We elucidate that phonon scattering by structural defects is related to the spatial fluctuations of local vibrational density of states (LVDOS). An effective method of tuning thermal transport in low-dimensional systems through the modulation of LVDOS has been proposed. Our findings provide insights into experimentally controlling thermal transport in nanoscale devices

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bright and dark modes induced by graphene bubbles

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

First-principles study of heat transport properties of graphene nanoribbons

We combine density-functional theory and the nonequilibrium Green's function method to study the thermal conductance of graphene nanoribbons with armchair and zigzag edges. Zigzag ribbons have higher thermal conductance than armchair ribbons of comparable widths due to an anisotropy in phonon dispersion for graphene nanoribbons: low-frequency bands in zigzag ribbons are more dispersive than those in armchair ribbons. The edges with and without hydrogen-passivation modify the atomic structure and ultimately influence the phonon thermal transport differently for the two ribbon types.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Geometrical effects on spin injection: 3D spin drift diffusion model

We discuss a three-dimensional (3D) spin drift diffusion (SDD) model to inject spin from a ferromagnet (FM) to a normal metal (N) or semiconductor (SC). Using this model we investigate the problem of spin injection into isotropic materials like GaAs and study the effect of FM contact area and SC thickness on spin injection. We find that in order to achieve detectable spin injection a small contact area or thick SC samples are essential for direct contact spin injection devices. We investigate the use of thin metal films (Cu) proposed by S.B. Kumar et al. and show that they are an excellent substitute for tunnelling barriers (TB) in the regime of small contact area. Since most tunnelling barriers are prone to pinhole defects, we study the effect of pinholes in AlO tunnelling barriers and show that the reduction in the spin-injection ratio ($γ$) is solely due to the effective area of the pinholes and there is no correlation between the number of pinholes and the spin injection ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Graphene-based tortional resonator from molecular dynamics simulation

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Minimum thermal conductance in graphene and boron nitride superlattice

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Molecular dynamics simulation for heat transport in thin diamond nanowires

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermal expansion in multiple layers of graphene

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Transient behavior of full counting statistics in thermal transport

The generating function of energy counting statistics is derived for phononic junction systems. It is expressed in terms of the contour-ordered self-energy of the lead with shifted arguments, $Σ^A(τ, τ') = Σ_L\bigl(τ+\hbar x(τ), τ' + \hbar x(τ')\bigr) - Σ_L(τ, τ')$, where $Σ_L(τ,τ')$ is the usual contour-ordered self-energy of the left lead. The cumulants of the energy transferred in a given time $t_M$ from the lead to the center is obtained by taking derivatives. A transient result of the first four cumulants of a graphene junction is presented. It is found that measurements cause the energy to flow into the lead.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mechanism of phonon localized edge modes

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Transient behavior of heat transport in a thermal switch

We study the time-dependent transport of heat in a nanoscale thermal switch. The switch consists of left and right leads that are initially uncoupled. During switch-on the coupling between the leads is abruptly turned on. We use the nonequilibrium Green's function formalism and numerically solve the constructed Dyson equation to determine the nonperturbative heat current. At the transient regime we find that the current initially flows simultaneously into both of the leads and then afterwards oscillates between flowing into and out of the leads. At later times the oscillations decay away and the current settles into flowing from the hotter to the colder lead. We find the transient behavior to be influenced by the extra energy added during switch-on. Such a transient behavior also exists even when there is no temperature difference between the leads. The current at the long-time limit approaches the steady-state value independently calculated from the Landauer formula.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tunable heat pump by modulating the coupling to the leads

We follow the nonequilibrium Green's function formalism to study time-dependent thermal transport in a linear chain system consisting of two semi-infinite leads connected together by a coupling that is harmonically modulated in time. The modulation is driven by an external agent that can absorb and emit energy. We determine the energy current flowing out of the leads exactly by solving numerically the Dyson equation for the contour-ordered Green's function. The amplitude of the modulated coupling is of the same order as the interparticle coupling within each lead. When the leads have the same temperature, our numerical results show that modulating the coupling between the leads may direct energy to either flow into the leads simultaneously or flow out of the leads simultaneously, depending on the values of the driving frequency and temperature. A special combination of values of the driving frequency and temperature exists wherein no net energy flows into or out of the leads, even for long times. When one of the leads is warmer than the other, net energy flows out of the warmer lead. For the cooler lead, however, the direction of the energy current flow depends on the values of the driving frequency and temperature. In addition, we find transient effects to become more pronounced for higher values of the driving frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Quantum transport in honeycomb lattice ribbons with armchair and zigzag edges coupled to semi-infinite linear chain leads

We study quantum transport in honeycomb lattice ribbons with either armchair or zigzag edges. The ribbons are coupled to semi-infinite linear chains serving as the input and output leads and we use a tight-binding Hamiltonian with nearest-neighbor hops. For narrow ribbons we find transmission gaps for both types of edges. The center of the gap is at the middle of the band in ribbons with armchair edges. This symmetry is due to a property satisfied by the matrices in the resulting linear problem. In ribbons with zigzag edges the gap center is displaced to the right of the middle of the band. We also find transmission oscillations and resonances within the transmitting region of the band for both types of edges. Extending the length of a ribbon does not affect the width of the transmission gap, as long as the ribbon's length is longer than a critical value when the gap can form. Increasing the width of the ribbon, however, changes the width of the gap. In armchair edges the gap is not well-defined because of the appearance of transmission resonances while in zigzag edges the gap width systematically shrinks as the width of the ribbon is increased. We also find only evanescent waves within the gap and both evanescent and propagating waves in the transmitting regions.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Nonequilibrium Green's function method for thermal transport in junctions

We present a detailed treatment of the nonequilibrium Green's function method for thermal transport due to atomic vibrations in nanostructures. Some of the key equations, such as self-energy and conductance with nonlinear effect, are derived. A self-consistent mean-field theory is proposed. Computational procedures are discussed. The method is applied to a number of systems including one-dimensional chains, a benzene ring junction, and carbon nanotubes. Mean-field calculations of the Fermi-Pasta-Ulam model are compared with classical molecular dynamics simulations. We find that nonlinearity suppresses thermal transport even at moderately high temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A Worm Algorithm for Two-Dimensional Spin Glasses

A worm algorithm is proposed for the two-dimensional spin glasses. The method is based on a low-temperature expansion of the partition function. The low-temperature configurations of the spin glass on square lattice can be viewed as strings connecting pairs of frustrated plaquettes. The worm algorithm directly manipulates these strings. It is shown that the worm algorithm is as efficient as any other types of cluster or replica-exchange algorithms. The worm algorithm is even more efficient if free boundary conditions are used. We obtain accurate low-temperature specific heat data consistent with a form c = T^{-2} exp(-2J/(k_BT)), where T is temperature and J is coupling constant, for the +/-J two-dimensional spin glass.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Transition Matrix Monte Carlo Method

We present a formalism of the transition matrix Monte Carlo method. A stochastic matrix in the space of energy can be estimated from Monte Carlo simulation. This matrix is used to compute the density of states, as well as to construct multi-canonical and equal-hit algorithms. We discuss the performance of the methods. The results are compared with single histogram method, multi-canonical method, and other methods. In many aspects, the present method is an improvement over the previous methods. PACS numbers: 02.70.Tt, 05.10.Ln, 05.50.+q. Keywords: Monte Carlo method, flat histogram, multi-canonical ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Anisotropic finite-size scaling analysis of a three-dimensional driven-diffusive system

We study the standard three-dimensional driven diffusive system on a simple cubic lattice where particle jumps along a given lattice direction are biased by an infinitely strong field, while those along other directions follow the usual Kawasaki dynamics. Our goal is to determine which of the several existing theories for critical behavior is valid. We analyze finite-size scaling properties using a range of system shapes and sizes far exceeding previous studies. Four different analytic predictions are tested against the numerical data. Binder and Wang's prediction does not fit the data well. Among the two slightly different versions of Leung, the one including the effects of a dangerous irrelevant variable appears to be better. Recently proposed isotropic finite-size scaling is inconsistent with our data from cubic systems, where systematic deviations are found, especially in scaling at the critical temperature.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Transition Matrix Monte Carlo

Although histogram methods have been extremely effective for analyzing data from Monte Carlo simulations, they do have certain limitations, including the range over which they are valid and the difficulties of combining data from independent simulations. In this paper, we describe an complementary approach to extracting information from Monte Carlo simulations that uses the matrix of transition probabilities. Combining the Transition Matrix with an N-fold way simulation technique produces an extremely flexible and efficient approach to rather general Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Random Sequential Adsorption, Series Expansion and Monte Carlo Simulation

Random sequential adsorption is an irreversible surface deposition of extended objects. In systems with continuous degrees of freedom coverage follows a power law, theta(t) = theta_J - c t^{-alpha}, where the exponent alpha depends on the geometric shape (symmetry) of the objects. Lattice models give typically exponential saturation to jamming coverage. We discuss how such function theta(t) can be computed by series expansions and analyzed with Pade approximations. We consider the applications of efficient Monte Carlo computer simulation method (event-driven method) to random sequential adsorptions with high precision and at very long-time scale.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Locally Frozen Defects in Random Sequential Adsorption with Diffusional Relaxation

Random sequential adsorption with diffusional relaxation, of two by two square objects on the two-dimensional square lattice is studied by Monte Carlo computer simulation. Asymptotically for large lattice sizes, diffusional relaxation allows the deposition process to reach full coverage. The coverage approaches the full occupation value, 1, as a power-law with convergence exponent near 1/2. For a periodic lattice of finite (even) size $L$, the final state is a frozen random rectangular grid of domain walls connecting single-site defects. The domain sizes saturate at L**0.8. Prior to saturation, i.e., asymptotically for infinite lattice, the domain growth is power-law with growth exponent near, or possibly somewhat smaller than, 1/2.

0 citations0 reviewsNo code linkedOpen access