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

Abhishek Dhar

Abhishek Dhar appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softcond-mat.mes-hallquant-phcond-mat.dis-nncond-mat.supr-conmath-phmath.MPPopulations and EvolutionBiological Physicscond-mat.mtrl-scicond-mat.quant-gaseess.SYhep-thmath.HOmath.OCnlin.CDphysics.class-phphysics.comp-phphysics.soc-phSystems and Control

Catalog footprint

What is connected

56works
21topics
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

56 published item(s)

preprint2026arXiv

Impurity dynamics in a zero-temperature gas

If energy is suddenly released in a localized region of space uniformly filled with identical stationary hard spheres, the outcome is a blast with an asymptotically spherical shock wave separating moving and stationary hard spheres. The radius $R(t)$ of the region filled with the moving spheres grows as $t^{2/(d+2)}$, where $d$ is the spatial dimension. The simplest way to inject energy is to kick a few `impurity' particles. Using hydrodynamics and kinetic theory, we argue that the typical displacement of an impurity scales as $R_{\rm imp} \sim λ(R/λ)^{(4+3d^2)/(8+3d^2)}$, where $λ$ is the mean-free path in the initial state. The number of collisions experienced by each impurity grows as $(R/λ)^{(8+2d^2)/(8+3d^2)}$, while its average speed decreases as $t^{-d(8-2d+3d^2)/[(2+d)(8+3d^2)]}$. In $2D$, the predictions for impurity displacement, collision numbers, and speed are $t^{2/5},~t^{2/5}$ and $t^{-2/5}$, respectively. These predictions are in reasonable agreement with the results of molecular dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Adaptive Output Feedback Model Predictive Control

Model predictive control (MPC) for uncertain systems in the presence of hard constraints on state and input is a non-trivial problem, and the challenge is increased manyfold in the absence of state measurements. In this paper, we propose an adaptive output feedback MPC technique, based on a novel combination of an adaptive observer and robust MPC, for single-input single-output discrete-time linear time-invariant systems. At each time instant, the adaptive observer provides estimates of the states and the system parameters that are then leveraged in the MPC optimization routine while robustly accounting for the estimation errors. The solution to the optimization problem results in a homothetic tube where the state estimate trajectory lies. The true state evolves inside a larger outer tube obtained by augmenting a set, invariant to the state estimation error, around the homothetic tube sections. The proof for recursive feasibility for the proposed `homothetic and invariant' two-tube approach is provided, along with simulation results on an academic system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fundamental limitations in Lindblad descriptions of systems weakly coupled to baths

It is very common in the literature to write down a Markovian quantum master equation in Lindblad form to describe a system with multiple degrees of freedom and weakly connected to multiple thermal baths which can, in general, be at different temperatures and chemical potentials. However, the microscopically derived quantum master equation up to leading order in system-bath coupling is of the so-called Redfield form which is known to not preserve complete positivity in most cases. Additional approximations to the Redfield equation are required to obtain a Lindblad form. We lay down some fundamental requirements for any further approximations to Redfield equation, which, if violated, leads to physical inconsistencies like inaccuracies in the leading order populations and coherences in the energy eigenbasis, violation of thermalization, violation of local conservation laws at the non-equilibrium steady state (NESS). We argue that one or more of these conditions will generically be violated in all the weak system-bath-coupling Lindblad descriptions existing in literature to our knowledge. As an example, we study the recently derived Universal Lindblad Equation (ULE) and use these conditions to show violation of local conservation laws due to inaccurate coherences but accurate populations in energy eigenbasis. Finally, we exemplify our analytical results numerically in an interacting open quantum spin system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gap Statistics for Confined Particles with Power-Law Interactions

We consider the $N$ particle classical Riesz gas confined in a one-dimensional external harmonic potential with power law interaction of the form $1/r^k$ where $r$ is the separation between particles. As special limits it contains several systems such as Dyson's log-gas ($k\to 0^+$), Calogero-Moser model ($k=2$), 1d one component plasma ($k=-1$) and the hard-rod gas ($k\to \infty$). Despite its growing importance, only large-$N$ field theory and average density profile are known for general $k$. In this Letter, we study the fluctuations in the system by looking at the statistics of the gap between successive particles. This quantity is analogous to the well-known level spacing statistics which is ubiquitous in several branches of physics. We show that the variance goes as $N^{-b_k}$ and we find the $k$ dependence of $b_k$ via direct Monte Carlo simulations. We provide supporting arguments based on microscopic Hessian calculation and a quadratic field theory approach. We compute the gap distribution and study its system size scaling. Except in the range $-1<k<0$, we find scaling for all $k>-2$ with both Gaussian and non-Gaussian scaling forms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Run-and-tumble particle in one-dimensional confining potential: Steady state, relaxation and first passage properties

We study the dynamics of a one-dimensional run and tumble particle subjected to confining potentials of the type $V(x) = α\, |x|^p$, with $p>0$. The noise that drives the particle dynamics is telegraphic and alternates between $\pm 1$ values. We show that the stationary probability density $P(x)$ has a rich behavior in the $(p, α)$-plane. For $p>1$, the distribution has a finite support in $[x_-,x_+]$ and there is a critical line $α_c(p)$ that separates an active-like phase for $α> α_c(p)$ where $P(x)$ diverges at $x_\pm$, from a passive-like phase for $α< α_c(p)$ where $P(x)$ vanishes at $x_\pm$. For $p<1$, the stationary density $P(x)$ collapses to a delta function at the origin, $P(x) = δ(x)$. In the marginal case $p=1$, we show that, for $α< α_c$, the stationary density $P(x)$ is a symmetric exponential, while for $α> α_c$, it again is a delta function $P(x) = δ(x)$. For the special cases $p=2$ and $p=1$, we obtain exactly the full time-dependent distribution $P(x,t)$, that allows us to study how the system relaxes to its stationary state. In addition, in these two cases, we also study analytically the full distribution of the first-passage time to the origin. Numerical simulations are in complete agreement with our analytical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

COVID-19: Analytic results for a modified SEIR model and comparison of different intervention strategies

The Susceptible-Exposed-Infected-Recovered (SEIR) epidemiological model is one of the standard models of disease spreading. Here we analyse an extended SEIR model that accounts for asymptomatic carriers, believed to play an important role in COVID-19 transmission. For this model we derive a number of analytic results for important quantities such as the peak number of infections, the time taken to reach the peak and the size of the final affected population. We also propose an accurate way of specifying initial conditions for the numerics (from insufficient data) using the fact that the early time exponential growth is well-described by the dominant eigenvector of the linearized equations. Secondly we explore the effect of different intervention strategies such as social distancing (SD) and testing-quarantining (TQ). The two intervention strategies (SD and TQ) try to reduce the disease reproductive number, $R_0$, to a target value $R^{\rm target}_0 < 1$, but in distinct ways, which we implement in our model equations. We find that for the same $R^{\rm target}_0 < 1$, TQ is more efficient in controlling the pandemic than SD. However, for TQ to be effective, it has to be based on contact tracing and our study quantifies the required ratio of tests-per-day to the number of new cases-per-day. Our analysis shows that the largest eigenvalue of the linearised dynamics provides a simple understanding of the disease progression, both pre- and post- intervention, and explains observed data for many countries. We apply our results to the COVID data for India to obtain heuristic projections for the course of the pandemic, and note that the predictions strongly depend on the assumed fraction of asymptomatic carriers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Equivalence of NEGF and scattering approaches to electron transport in the Kitaev chain

We consider electron transport in a Kitaev chain connected at its two ends to normal metallic leads kept at different temperatures and chemical potentials. Transport in this set-up is usually studied using two frameworks -- the nonequilibrium Green's function (NEGF) approach or the scattering approach. In the NEGF approach the current and other steady state properties of a system are expressed in terms of Green's functions that involve the wire properties and self-energy corrections arising from the leads. In the scattering approach, transport is studied in terms of the scattering amplitudes of plane waves incident on the wire from the reservoirs. Here we show explicitly that these two approaches produce identical results for the conductance of the Kitaev chain. Further we show that the NEGF expression for conductance can be written in such a way that there is a one-to-one correspondence of the various terms in the NEGF expression to the amplitudes for normal transmission, Andreev transmission and Andreev reflection in the scattering approach. Thereby, we obtain closed form expressions for these. We obtain the wavefunctions of zero energy Majorana bound states(MBS) of the wire connected to leads and prove that they are present in the same parameter regime in which they occur for an isolated wire. These bound states give rise to perfect Andreev reflection responsible for zero bias quantized conductance peak. We discuss the dependence of the width of this peak on different parameters of the Hamiltonian and relate it to the MBS wavefunction properties. We find that the peak broadens if the weight of the MBS in the reservoirs increases and vice versa.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Revisiting the Mazur bound and the Suzuki equality

Among the few known rigorous results for time-dependent equilibrium correlations, important for understanding transport properties, are the Mazur bound and the Suzuki equality. The Mazur inequality gives a lower bound, on the long-time average of the time-dependent auto-correlation function of observables, in terms of equilibrium correlation functions involving conserved quantities. On the other hand, Suzuki proposes an exact equality for quantum systems. In this paper, we discuss the relation between the two results and in particular, look for the analogue of the Suzuki result for classical systems. This requires us to examine as to what constitutes a complete set of conserved quantities required to saturate the Mazur bound. We present analytic arguments as well as illustrative numerical results from a number of different systems. Our examples include systems with few degrees of freedom as well as many-particle integrable models, both free and interacting.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A critique of the Covid-19 analysis for India by Singh and Adhikari

In a recent paper (arXiv:2003.12055), Singh and Adhikari present results of an analysis of a mathematical model of epidemiology based on which they argue that a $49$ day lockdown is required in India for containing the pandemic in India. We assert that as a model study with the stated assumptions, the analysis presented in the paper is perfectly valid, however, any serious comparison with real data and attempts at prediction from the model are highly problematic. The main point of the present paper is to convincingly establish this assertion while providing a warning that the results and analysis of such mathematical models should not be taken at face value and need to be used with great caution by policy makers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gap statistics of two interacting run and tumble particles in one dimension

We study the dynamics of the separation (gap) between a pair of interacting run and tumble particles (RTPs) moving in one dimension in the presence of additional thermal noise. On a ring geometry the distribution of the gap approaches a steady state. We analytically compute this distribution and find that this is exponentially localised in space, in contrast to the `jammed' configuration, seen earlier in the absence of thermal noise. We also study the relaxation which is an exponential, characterised by a time scale $τ_r$. We observe that this time scale undergoes a crossover from a size independent value to a size dependent form with increasing size $l$ of the ring. We study the full eigenvalue spectrum of the evolution operator $\mathcal{L}$ and find that the spectrum can be classified into four sectors depending on the symmetries of $\mathcal{L}$. For large $l$, we find explicit expressions for the low lying eigenvalues in each of the sectors. On infinite line the separation does not reach a steady state. In the long times we find that the particles behave as interacting Brownian particles, except for the presence of a peak in the distribution at small separation which is a remnant of activity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How Archimedes showed that $π$ is approximately equal to 22/7

The ratio of the circumference, C, of a circle to its diameter, D, is a constant number denoted by $π$ and is independent of the size of the circle. It is known that $π$ is an irrational number and therefore cannot be expressed as a common fraction. Its value is approximately equal to 3.141592. Since Archimedes was one of the first persons to suggest a rational approximation of 22/7 for $π$, it is sometimes referred to as Archimedes' constant. In this article, we discuss how Archimedes came up with his formula. Archimedes in fact proved that 223/71 < $π$ < 22/7. Here we provide an improved lower bound.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nonlinear Fluctuating Hydrodynamics for the Classical XXZ Spin Chain

Using the framework of nonlinear fluctuating hydrodynamics (NFH), we examine equilibrium spatio-temporal correlations in classical ferromagnetic spin chains with nearest neighbor interactions. In particular, we consider the classical XXZ-Heisenberg spin chain (also known as Lattice Landau Lifshitz or LLL model) evolving deterministically and chaotically via Hamiltonian dynamics, for which energy and $z$-magnetization are the only locally conserved fields. For the easy-plane case, this system has a low-temperature regime in which the difference between neighboring spin's angular orientations in the XY plane is an \textit{almost conserved} field. According to the predictions of NFH, the dynamic correlations in this regime exhibit a heat peak and propagating sound peaks, all with anomalous broadening. We present a detailed molecular dynamics test of these predictions and find a reasonably accurate verification. We find that, in a suitable intermediate temperature regime, the system shows two sound peaks with Kardar-Parisi-Zhang (KPZ) scaling and a heat peak where the expected anomalous broadening is less clear. In high temperature regimes of both easy plane and easy axis case of LLL, our numerics show clear diffusive spin and energy peaks and absence of any sound modes, as one would expect. We also simulate an integrable version of the XXZ-model, for which the ballistic component instead moves with a broad range of speeds rather than being concentrated in narrower peaks around the sound speed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Steady State of an Active Brownian Particle in Two-Dimensional Harmonic Trap

We find an exact series solution for the steady-state probability distribution of a harmonically trapped active Brownian particle in two dimensions, in the presence of translational diffusion. This series solution allows us to efficiently explore the behavior of the system in different parameter regimes. Identifying "active" and "passive" regimes, we predict a surprising re-entrant active-to-passive transition with increasing trap stiffness. Our numerical simulations validate this finding. We discuss various interesting limiting cases wherein closed form expressions for the distributions can be obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermalization of local observables in the $α$-FPUT chain

Most studies on the problem of equilibration of the Fermi-Pasta-Ulam-Tsingou (FPUT) system have focused on equipartition of energy being attained amongst the normal modes of the corresponding harmonic system. In the present work, we instead discuss the equilibration problem in terms of local variables, and consider initial conditions corresponding to spatially localized energy. We estimate the time-scales for equipartition of space localized degrees of freedom and find significant differences with the times scales observed for normal modes. Measuring thermalization in classical systems necessarily requires some averaging, and this could involve one over initial conditions or over time or spatial averaging. Here we consider averaging over initial conditions chosen from a narrow distribution in phase space. We examine in detail the effect of the width of the initial phase space distribution, and of integrability and chaos, on the time scales for thermalization. We show how thermalization properties of the system, quantified by its equilibration time, defined in this work, can be related to chaos, given by the maximal Lyapunov exponent. Somewhat surprisingly we also find that the ensemble averaging can lead to thermalization of the integrable Toda chain, though on much longer time scales.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Transport in spinless superconducting wires

We consider electron transport in a model of a spinless superconductor described by a Kitaev type lattice Hamiltonian where the electron interactions are modelled through a superconducting pairing term. The superconductor is sandwiched between two normal metals kept at different temperatures and chemical potentials and are themselves modelled as non-interacting spinless fermions. For this set-up we compute the exact steady state properties of the system using the quantum Langevin equation approach. Closed form exact expressions for particle current, energy current and other two-point correlations are obtained in the Landauer-type forms and involve two nonequilibrium Green's functions. The current expressions are found out to be sum of three terms having simple physical interpretations. We then discuss a numerical approach where we construct the time-evolution of the two point correlators of the system from the eigenspectrum of the complete quadratic Hamiltonian describing the system and leads. By starting from an initial state corresponding to the leads in thermal equilibrium and the system in an arbitrary state, the long time solution for the correlations, before recurrence times, gives us steady state properties. We use this independent numerical method for verifying the results of the exact solution. We also investigate analytically the presence of high energy bound states and obtain expressions for their contributions to two point correlators. As applications of our general formalism we present results on thermal conductance and on the conductance of a Kitaev chain with next nearest neighbour interactions which allows topological phases with different winding numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Transport, correlations, and chaos in a classical disordered anharmonic chain

We explore transport properties in a disordered nonlinear chain of classical harmonic oscillators and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $N$, disorder strength $Δ$, and temperature $T$. The conductivity $κ_N$, obtained for finite length ($N$) systems, saturates to a value $κ_\infty >0$ in the large $N$ limit, for all values of disorder strength $Δ$ and temperature $T>0$. We show evidence that for any $Δ>0$ the conductivity goes to zero faster than any power of $T$ in the $(T/Δ) \to 0$ limit, and find that the form $κ_\infty \sim e^{-B |\ln(C Δ/T)|^3}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multi-oscillator chaotic islands. The finite-size effect can be $κ_N < κ_{\infty}$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_{\infty}$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the Out-of-Time-Ordered-Commutator.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universal scaling in active single-file dynamics

We study the single-file dynamics of three classes of active particles: run-and-tumble particles, active Brownian particles and active Ornstein-Uhlenbeck particles. At high activity values, the particles, interacting via purely repulsive and short-ranged forces, aggregate into several motile and dynamical clusters of comparable size, and do not display bulk phase-segregation. In this dynamical steady-state, we find that the cluster size distribution of these aggregates is a scaled function of the density and activity parameters across the three models of active particles with the same scaling function. The velocity distribution of these motile clusters is non-Gaussian. We show that the effective dynamics of these clusters can explain the observed emergent scaling of the mean-squared displacement of tagged particles for all the three models with identical scaling exponents and functions. Concomitant with the clustering seen at high activities, we observe that the static density correlation function displays rich structures, including multiple peaks that are reminiscent of particle clustering induced by effective attractive interactions, while the dynamical variant shows non-diffusive scaling. Our study reveals a universal scaling behavior in the single-file dynamics of interacting active particles.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Thermal Conductivity of Glass-Forming Liquids

Thermal conductivity of a model glass-forming system in the liquid and glass states is studied using extensive numerical simulations. We show that near the glass transition temperture, where the structural relaxation time becomes very long, the measured thermal conductivity decreases with increasing age. Secondly the thermal conductivity of the disordered solid obtained at low temperatures depends on the cooling rate with which it was prepared, with lower cooling rates leading to lower thermal conductivity. Our analysis links this decrease of the thermal conductivity with increased exploration of lower-energy inherent structures of the underlying potential energy landscape. Further we show that the lowering of conductivity for lower-energy inherent structures is related to the high frequency harmonic modes associated with the inherent structure being less extended.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Eigenstate Gibbs Ensemble in Integrable Quantum Systems

The Eigenstate Thermalization Hypothesis implies that for a thermodynamically large system in one of its eigenstates, the reduced density matrix describing any finite subsystem is determined solely by a set of {\it relevant} conserved quantities. In a generic system, only the energy plays that role and hence eigenstates appear locally thermal. Integrable systems, on the other hand, possess an extensive number of such conserved quantities and hence the reduced density matrix requires specification of an infinite number of parameters (Generalized Gibbs Ensemble). However, here we show by unbiased statistical sampling of the individual eigenstates with a given finite energy density, that the local description of an overwhelming majority of these states of even such an integrable system is actually Gibbs-like, i.e. requires only the energy density of the eigenstate. Rare eigenstates that cannot be represented by the Gibbs ensemble can also be sampled efficiently by our method and their local properties are then shown to be described by appropriately {\it truncated} Generalized Gibbs Ensembles. We further show that the presence of these rare eigenstates differentiates the model from the generic (non-integrable) case and leads to the system being described by a Generalized Gibbs Ensemble at long time under a unitary dynamics following a sudden quench, even when the initial state is a Gibbs-like eigenstate of the pre-quench Hamiltonian.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equilibrium dynamical correlations in the Toda chain and other integrable models

We investigate the form of equilibrium spatio-temporal correlation functions of conserved quantities, and of energy transport in the Toda lattice and in other integrable models. From numerical simulations we find that the correlations satisfy ballistic scaling with a remarkable collapse of data from different times. We examine special limiting choices of parameter values, for which the Toda lattice tends to either the harmonic chain or the equal mass hard-particle gas. In both these limiting cases, one can obtain the correlations exactly and we find excellent agreement with the direct Toda simulation results. We also discuss a transformation to "normal mode" variables, as commonly done in hydrodynamic theory of non-integrable systems, and find that this is useful, to some extent, even for the integrable system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Non-linear transport in an out-of-equilibrium single-site Bose Hubbard model: scaling, rectification and time dynamics

Recent experiments in hybrid-quantum systems facilitate the potential realization of one of the most fundamental interacting Hamiltonian-Reservoir system, namely, the single-site Bose-Hubbard model coupled to two reservoirs at different temperatures. Using Redfield equations in Born-Markov approximation, we compute non-equilibrium average particle number, energy and currents beyond linear response regime, both time-dynamics and steady state and investigate its dependence on various tunable parameters analytically. We find interesting scaling laws in high temperature regimes that are independent of choice of bath spectral functions. We also demonstrate that the system shows very interesting particle and energy current rectification properties which can be controlled via the relative strength of interaction and temperatures, as well as via the degree of asymmetry in system-bath coupling. Specifically, we find inversion of direction of thermal rectification as a function of the relative strength of the interaction strength and the temperatures. We also show that, in the limit of low-temperature and high interaction strength our results are consistent with the non-equilibrium spin-Boson model. Our results are experimentally relevant not only to hybrid quantum systems, but also in other areas such as molecular junctions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Out of equilibrium open quantum systems: a comparison of approximate Quantum Master Equation approaches with exact results

We present the Born-Markov approximated Redfield quantum master equation (RQME) description for an open system of non-interacting particles (bosons or fermions) on an arbitrary lattice of $N$ sites in any dimension and weakly connected to multiple reservoirs at different temperatures and chemical potentials. The RQME can be reduced to the Lindblad equation, of various forms, by making further approximations. By studying the $N=2$ case, we show that RQME gives results which agree with exact analytical results for steady state properties and with exact numerics for time-dependent properties, over a wide range of parameters. In comparison, the Lindblad equations have a limited domain of validity in non-equilibrium. We conclude that it is indeed justified to use microscopically derived full RQME to go beyond the limitations of Lindblad equations in out-of-equilibrium systems. We also derive closed form analytical results for out-of-equilibrium time dynamics of two-point correlation functions. These results explicitly show the approach to steady state and thermalization. These results are experimentally relevant for cold atoms, cavity QED and far-from-equilibrium quantum dot experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detection of a quantum particle on a lattice under repeated projective measurements

We consider a quantum particle, moving on a lattice with a tight-binding Hamiltonian, which is subjected to measurements to detect it's arrival at a particular chosen set of sites. The projective measurements are made at regular time intervals $τ$, and we consider the evolution of the wave function till the time a detection occurs. We study the probabilities of its first detection at some time and conversely the probability of it not being detected (i.e., surviving) up to that time. We propose a general perturbative approach for understanding the dynamics which maps the evolution operator, consisting of unitary transformations followed by projections, to one described by a non-Hermitian Hamiltonian. For some examples, of a particle moving on one and two-dimensional lattices with one or more detection sites, we use this approach to find exact expressions for the survival probability and find excellent agreement with direct numerical results. A mean field model with hopping between all pairs of sites and detection at one site is solved exactly. For the one- and two-dimensional systems, the survival probability is shown to have a power-law decay with time, where the power depends on the initial position of the particle. Finally, we show an interesting and non-trivial connection between the dynamics of the particle in our model and the evolution of a particle under a non-Hermitian Hamiltonian with a large absorbing potential at some sites.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Driven inelastic Maxwell gases

We consider the inelastic Maxwell model, which consists of a collection of particles that are characterized by only their velocities, and evolving through binary collisions and external driving. At any instant, a particle is equally likely to collide with any of the remaining particles. The system evolves in continuous time with mutual collisions and driving taken to be point processes with rates $τ_c^-{1}$ and $τ_w^{-1}$ respectively. The mutual collisions conserve momentum and are inelastic, with a coefficient of restitution $r$. The velocity change of a particle with velocity $v$, due to driving, is taken to be $Δv=-(1+r_w) v+η$, mimicking the collision with a vibrating wall, where $r_w$ the coefficient of restitution of the particle with the "wall" and $η$ is Gaussian white noise. The Ornstein-Uhlenbeck driving mechanism given by $\frac{dv}{dt}=-Γv+η$ is found to be a special case of the driving modeled as a point process. Using both the continuum and discrete versions we show that while the equations for the one-particle and the two-particle velocity distribution functions do not close, the joint evolution equations of the variance and the two-particle velocity correlation functions close. With the exact formula for the variance we find that, for $r_w\ne-1$, the system goes to a steady state. On the other hand, for $r_w=-1$, the system does not have a steady state. Similarly, the system goes to a steady state for the Ornstein-Uhlenbeck driving with $Γ\not=0$, whereas for the purely diffusive driving ($Γ=0$), the system does not have a steady state.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exact probability distribution for the two-tag displacement in single-file motion

We consider a gas of point particles moving on the one-dimensional line with a hard-core inter-particle interaction that prevents particle crossings --- this is usually referred to as single-file motion. The individual particle dynamics can be arbitrary and they only interact when they meet. Starting from initial conditions such that particles are uniformly distributed, we observe the displacement of a tagged particle at time $t$, with respect to the initial position of another tagged particle, such that their tags differ by $r$. For $r=0$, this is the usual well studied problem of the tagged particle motion. Using a mapping to a non-interacting particle system we compute the exact probability distribution function for the two-tagged particle displacement, for general single particle dynamics. As by-products, we compute the large deviation function, various cumulants and, for the case of Hamiltonian dynamics, the two-particle velocity auto-correlation function.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pumping single-file colloids: Absence of current reversal

We consider the single-file motion of colloidal particles interacting via short-ranged repulsion and placed in a traveling wave potential, that varies periodically in time and space. Under suitable driving conditions, a directed time-averaged flow of colloids is generated. We obtain analytic results for the model using a perturbative approach to solve the Fokker-Planck equations. The predictions show good agreement with numerical simulations. We find peaks in the time-averaged directed current as a function of driving frequency, wavelength and particle density and discuss possible experimental realizations. Surprisingly, unlike a closely related exclusion dynamics on a lattice, the directed current in the present model does not show current reversal with density. A linear response formula relating current response to equilibrium correlations is also proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum time of arrival distribution in a simple lattice model

Imagine an experiment where a quantum particle inside a box is released at some time in some initial state. A detector is placed at a fixed location inside the box and its clicking signifies arrival of the particle at the detector. What is the \emph{time of arrival} (TOA) of the particle at the detector ? Within the paradigm of the measurement postulate of quantum mechanics, one can use the idea of projective measurements to define the TOA. We consider the setup where a detector keeps making instantaneous measurements at regular finite time intervals {\emph{till}} it detects the particle at some time $t$, which is defined as the TOA. This is a stochastic variable and, for a simple lattice model of a free particle in a one-dimensional box, we find interesting features such as power-law tails in its distribution and in the probability of survival (non-detection). We propose a perturbative calculational approach which yields results that compare very well with exact numerics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Role of conserved quantities in normal heat transport in one dimenison

Although one-dimensional systems that exhibit translational symmetry are generally believed to exhibit anomalous heat transport, previous work has shown that the model of coupled rotators on a one-dimensional lattice constitute a possible exception. We investigate the equilibrium spatiotemporal correlations of energy and momentum of the rotator model, and find that both these fields diffuse normally. The normal diffusion is explained within the framework of stochastic hydrodynamics by observing that the angle variables do not constitute a conserved field, which leads to the absence of long-wavelength currents in the system. As an outcome of our analysis, we propose some general criteria for normal transport based on the existence of conservation laws.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tagged particle diffusion in one-dimensional systems with Hamiltonian dynamics - II

We study various temporal correlation functions of a tagged particle in one-dimensional systems of interacting point particles evolving with Hamiltonian dynamics. Initial conditions of the particles are chosen from the canonical thermal distribution. The correlation functions are studied in finite systems, and their forms examined at short and long times. Various one-dimensional systems are studied. Results of numerical simulations for the Fermi-Pasta-Ulam chain are qualitatively similar to results for the harmonic chain, and agree unexpectedly well with a simple description in terms of linearized equations for damped fluctuating sound waves. Simulation results for the alternate mass hard particle gas reveal that - in contradiction to our earlier results [1] with smaller system sizes - the diffusion constant slowly converges to a constant value, in a manner consistent with mode coupling theories. Our simulations also show that the behaviour of the Lennard-Jones gas depends on its density. At low densities, it behaves like a hard-particle gas, and at high densities like an anharmonic chain. In all the systems studied, the tagged particle was found to show normal diffusion asymptotically, with convergence times depending on the system under study. Finite size effects show up at time scales larger than sound traversal times, their nature being system-specific.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Waiting for rare entropic fluctuations

Non-equilibrium fluctuations of various stochastic variables, such as work and entropy production, have been widely discussed recently in the context of large deviations, cumulants and fluctuation relations. Typically, one looks at the distribution of these observables, at large fixed time. To characterize the precise stochastic nature of the process, we here address the distribution in the time domain. In particular, we focus on the first passage time distribution (FPTD) of entropy production, in several realistic models. We find that the fluctuation relation symmetry plays a crucial role in getting the typical asymptotic behavior. Similarities and differences to the simple random walk picture are discussed. For a driven particle in the ring geometry, the mean residence time is connected to the particle current and the steady state distribution, and it leads to a fluctuation relation-like symmetry in terms of the FPTD.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-Energy Tail of the Velocity Distribution of Driven Inelastic Maxwell Gases

A model of homogeneously driven dissipative system, consisting of a collection of $N$ particles that are characterized by only their velocities, is considered. Adopting a discrete time dynamics, at each time step, a pair of velocities is randomly selected. They undergo inelastic collision with probability $p$. With probability $(1-p)$, energy of the system is changed by changing the velocities of both the particles independently according to $v\rightarrow -r_w v +η$, where $η$ is a Gaussian noise drawn independently for each particle as well as at each time steps. For the case $r_w=- 1$, although the energy of the system seems to saturate (indicating a steady state) after time steps of $O(N)$, it grows linearly with time after time steps of $O(N^2)$, indicating the absence of a eventual steady state. For $ -1 <r_w \leq 1$, the system reaches a steady state, where the average energy per particle and the correlation of velocities are obtained exactly. In the thermodynamic limit of large $N$, an exact equation is obtained for the moment generating function. In the limit of nearly elastic collisions and weak energy injection, the velocity distribution is shown to be a Gaussian. Otherwise, for $|r_w| < 1$, the high-energy tail of the velocity distribution is Gaussian, with a different variance, while for $r_w=+1$ the velocity distribution has an exponential tail.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical test of hydrodynamic fluctuation theory in the Fermi-Pasta-Ulam chain

Recent work has developed a nonlinear hydrodynamic fluctuation theory for a chain of coupled anharmonic oscillators governing the conserved fields, namely stretch, momentum, and energy. The linear theory yields two propagating sound modes and one diffusing heat mode. In contrast, the nonlinear theory predicts that, at long times, the sound mode correlations satisfy Kardar-Parisi-Zhang (KPZ) scaling, while the heat mode correlations satisfies Levy-walk scaling. In the present contribution we report on molecular dynamics simulations of Fermi-Pasta-Ulam chains to compute various spatiotemporal correlation functions and compare them with the predictions of the theory. We find very good agreement in many cases, but also some deviations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universal Large Deviations for the Tagged Particle in Single File Motion

We consider a gas of point particles moving in a one-dimensional channel with a hard-core inter-particle interaction that prevents particle crossings --- this is called single-file motion. Starting from equilibrium initial conditions we observe the motion of a tagged particle. It is well known that if the individual particle dynamics is diffusive, then the tagged particle motion is sub-diffusive, while for ballistic particle dynamics, the tagged particle motion is diffusive. Here we compute exactly the large deviation function for the tagged particle displacement and show that this is universal, independent of the individual dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anomalous transport and current fluctuations in a model of diffusing Levy walkers

A Levy walk is a non-Markovian stochastic process in which the elementary steps of the walker consist of motion with constant speed in randomly chosen directions and for a random period of time. The time of flight is chosen from a long-tailed distribution with a finite mean but an infinite variance. Here we consider an open system with boundary injection and removal of particles, at prescribed rates, and study the steady state properties of the system. In particular, we compute density profiles, current and current fluctuations in this system. We also consider the case of a finite density of Levy walkers on the ring geometry. Here we introduce a size dependent cut-off in the time of flight distribution and consider properties of current fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exact solution of a Levy walk model for anomalous heat transport

The Levy walk model is studied in the context of the anomalous heat conduction of one dimensional systems. In this model the heat carriers execute Levy-walks instead of normal diffusion as expected in systems where Fourier's law holds. Here we calculate exactly the average heat current, the large deviation function of its fluctuations and the temperature profile of the Levy-walk model maintained in a steady state by contact with two heat baths (the open geometry). We find that the current is non-locally connected to the temperature gradient. As observed in recent simulations of mechanical models, all the cumulants of the current fluctuations have the same system-size dependence in the open geometry. For the ring geometry, we argue that a size dependent cut-off time is necessary for the Levy walk model to behave as mechanical models. This modification does not affect the results on transport in the open geometry for large enough system sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Heat conduction in the α-β-Fermi-Pasta-Ulam chain

Recent simulation results on heat conduction in a one-dimensional chain with an asymmetric inter-particle interaction potential and no onsite potential found non-anomalous heat transport in accordance to Fourier's law. This is a surprising result since it was long believed that heat conduction in one-dimensional systems is in general anomalous in the sense that the thermal conductivity diverges as the system size goes to infinity. In this paper we report on detailed numerical simulations of this problem to investigate the possibility of a finite temperature phase transition in this system. Our results indicate that the unexpected results for asymmetric potentials is a result of insufficient chain length, and does not represent the asymptotic behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Landauer formula for phonon heat conduction: relation between energy transmittance and transmission coefficient

The heat current across a quantum harmonic system connected to reservoirs at different temperatures is given by the Landauer formula, in terms of an integral over phonon frequencies ω, of the energy transmittance T(ω). There are several different ways to derive this formula, for example using the Keldysh approach or the Langevin equation approach. The energy transmittance T(ω) is usually expressed in terms of nonequilibrium phonon Green's function and it is expected that it is related to the transmission coefficient τ(ω) of plane waves across the system. In this paper, for a one-dimensional set-up of a finite harmonic chain connected to reservoirs which are also semi-infinite harmonic chains, we present a simple and direct demonstration of the relation between T(ω) and τ(ω). Our approach is easily extendable to the case where both system and reservoirs are in higher dimensions and have arbitrary geometries, in which case the meaning of τ and its relation to T are more non-trivial.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-equilibrium stationary state of a harmonic crystal with alternating masses

We analyze the non-equilibrium steady states (NESS) of a one dimensional harmonic chain of $N$ atoms with alternating masses connected to heat reservoirs at unequal temperatures. We find that the temperature profile defined through the local kinetic energy $T(j) \equiv {<p_j^2>}/{m_j}$, oscillates with period two in the bulk of the system. Depending on boundary conditions, either the heavier or the lighter particles in the bulk are hotter. We obtain exact expressions for the bulk temperature profile and steady state current in the limit $N \rightarrow \infty$. These depend on whether $N$ is odd or even. We also study similar temperature oscillations in the NESS of systems with noise in the dynamics. These die out as $N \rightarrow \infty$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonequilibrium density matrix description of steady state quantum transport

With this work we investigate the stationary nonequilibrium density matrix of current carrying nonequilibrium steady states of in-between quantum systems that are connected to reservoirs. We describe the analytical procedure to obtain the explicit result for the reduced density matrix of quantum transport when the system, the connecting reservoirs and, as well, the system-reservoir interactions are described by quadratic Hamiltonians. Our procedure is detailed for both, electronic transport described by the tight-binding Hamiltonian and for phonon transport described by harmonic Hamiltonians. For the special case of weak system-reservoir couplings, a more detailed description of the steady-state density matrix is obtained. Several paradigm transport setups for inter-electrode electron transport and low-dimensional phonon heat flux are elucidated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Scaling behavior in the convection-driven Brazil-nut effect

The Brazil-nut effect is the phenomenon in which a large intruder particle immersed in a vertically shaken bed of smaller particles rises to the top, even when it is much denser. The usual practice, while describing these experiments, has been to use the dimensionless acceleration Γ=a ω^2/g, where a and ωare respectively the amplitude and the angular frequency of vibration and g is the acceleration due to gravity. Considering a vibrated quasi-two-dimensional bed of mustard seeds, we show here that the peak-to-peak velocity of shaking v= aω, rather than Γ, is the relevant parameter in the regime where boundary-driven granular convection is the main driving mechanism. We find that the rise-time τof an intruder is described by the scaling law τ~ (v-v_c)^{-α}, where v_c is identified as the critical vibration velocity for the onset of convective motion of the mustard seeds. This scaling form holds over a wide range of (a,ω), diameter and density of the intruder.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tagged particle diffusion in one-dimensional gas with Hamiltonian dynamics

We consider a one-dimensional gas of hard point particles in a finite box that are in thermal equilibrium and evolving under Hamiltonian dynamics. Tagged particle correlation functions of the middle particle are studied. For the special case where all particles have the same mass, we obtain analytic results for the velocity auto-correlation function in the short time diffusive regime and the long time approach to the saturation value when finite-size effects become relevant. In the case where the masses are unequal, numerical simulations indicate sub-diffusive behaviour with mean square displacement of the tagged particle growing as t/ln(t) with time t. Also various correlation functions, involving the velocity and position of the tagged particle, show damped oscillations at long times that are absent for the equal mass case.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Additivity Principle in High-dimensional Deterministic Systems

The additivity principle (AP), conjectured by Bodineau and Derrida [Phys. Rev. Lett. vol.92, 180601 (2004)], is discussed for the case of heat conduction in three-dimensional disordered harmonic lattices to consider the effects of deterministic dynamics, higher dimensionality, and different transport regimes, i.e., ballistic, diffusive, and anomalous transport. The cumulant generating function (CGF) for heat transfer is accurately calculated, and compared with the one given by the AP. In the diffusive regime, we find a clear agreement with the conjecture even if the system is high-dimensional. Surprisingly even in the anomalous regime the CGF is also well fitted by the AP. Lower dimensional systems are also studied and the importance of three-dimensionality for the validity is stressed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Application of importance sampling to the computation of large deviations in non-equilibrium processes

We present an algorithm for finding the probabilities of rare events in nonequilibrium processes. The algorithm consists of evolving the system with a modified dynamics for which the required event occurs more frequently. By keeping track of the relative weight of phase-space trajectories generated by the modified and the original dynamics one can obtain the required probabilities. The algorithm is tested on two model systems of steady-state particle and heat transport where we find a huge improvement from direct simulation methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Generating Function Formula of Heat Transfer in Harmonic Networks

We consider heat transfer across an arbitrary harmonic network connected to two heat baths at different temperatures. The network has $N$ positional degrees of freedom, of which $N_L$ are connected to a bath at temperature $T_L$ and $N_R$ are connected to a bath at temperature $T_R$. We derive an exact formula for the cumulant generating function for heat transfer between the two baths. The formula is valid even for $N_L \ne N_R$ and satisfies the Gallavotti-Cohen fluctuation symmetry. Since harmonic crystals in three dimensions are known to exhibit different regimes of transport such as ballistic, anomalous and diffusive, our result implies validity of the fluctuation theorem in all regimes. Our exact formula provides a powerful tool to study other properties of nonequilibrium current fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Heat conduction in disordered harmonic lattices with energy conserving noise

We study heat conduction in a harmonic crystal whose bulk dynamics is supplemented by random reversals (flips) of the velocity of each particle at a rate $λ$. The system is maintained in a nonequilibrium stationary state(NESS) by contacts with Langevin reservoirs at different temperatures. We show that the one-body and pair correlations in this system are the same (after an appropriate mapping of parameters) as those obtained for a model with self-consistent reservoirs. This is true both for the case of equal and random(quenched) masses. While the heat conductivity in the NESS of the ordered system is known explicitly, much less is known about the random mass case. Here we investigate the random system, with velocity flips. We improve the bounds on the Green-Kubo conductivity obtained by C.Bernardin. The conductivity of the 1D system is then studied both numerically and analytically. This sheds some light on the effect of noise on the transport properties of systems with localized states caused by quenched disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Large deviations of heat flow in harmonic chains

We consider heat transport across a harmonic chain connected at its two ends to white-noise Langevin reservoirs at different temperatures. In the steady state of this system the heat $Q$ flowing from one reservoir into the system in a finite time $τ$ has a distribution $P(Q,τ)$. We study the large time form of the corresponding moment generating function $<e^{-λQ}>\sim g(λ) e^{τμ(λ)}$. Exact formal expressions, in terms of phonon Green's functions, are obtained for both $μ(λ)$ and also the lowest order correction $g(λ)$. We point out that, in general a knowledge of both $μ(λ)$ and $g(λ)$ is required for finding the large deviation function associated with $P(Q,τ)$. The function $μ(λ)$ is known to be the largest eigenvector of an appropriate Fokker-Planck type operator and our method also gives the corresponding eigenvector exactly.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tunable thermal conductivity in defect engineered nanowires at low temperatures

We measure the thermal conductivity ($κ$) of individual InAs nanowires (NWs), and find that it is 3 orders of magnitude smaller than the bulk value in the temperature range of 10 to 50 K. We argue that the low $κ$ arises from the strong localization of phonons in the random superlattice of twin-defects oriented perpendicular to the axis of the NW. We observe significant electronic contribution arising from the surface accumulation layer which gives rise to tunability of $κ$ with the application of electrostatic gate and magnetic field. Our devices and measurements of $κ$ at different carrier concentrations and magnetic field without introducing structural defects, offer a means to study new aspects of nanoscale thermal transport.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Energy current magnification in coupled oscillator loops

Motivated by studies on current magnification in quantum mesoscopic systems we consider sound and heat transmission in classical models of oscillator chains. A loop of coupled oscillators is connected to two leads through which one can either transmit monochromatic waves or white noise signal from heat baths. We look for the possibility of current magnification in this system due to some asymmetry introduced between the two arms in the loop. We find that current magnification is indeed obtained for particular frequency ranges. However the integrated current shows the effect only in the presence of a pinning potential for the atoms in the leads. We also study the effect of anharmonicity on current magnification.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Heat conduction and phonon localization in disordered harmonic crystals

We investigate the steady state heat current in two and three dimensional isotopically disordered harmonic lattices. Using localization theory as well as kinetic theory we estimate the system size dependence of the current. These estimates are compared with numerical results obtained using an exact formula for the current given in terms of a phonon transmission function, as well as by direct nonequilibrium simulations. We find that heat conduction by high-frequency modes is suppressed by localization while low-frequency modes are strongly affected by boundary conditions. Our {\color{black}heuristic} arguments show that Fourier's law is valid in a three dimensional disordered solid except for special boundary conditions. We also study the pinned case relevant to localization in quantum systems and often used as a model system to study the validity of Fourier's law. Here we provide the first numerical verification of Fourier's law in three dimensions. In the two dimensional pinned case we find that localization of phonon modes leads to a heat insulator.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Heat conduction in a three dimensional anharmonic crystal

We perform nonequilibrium simulations of heat conduction in a three dimensional anharmonic lattice. By studying slabs of length N and width W, we examine the cross-over from one-dimensional to three dimensional behavior of the thermal conductivity. We find that for large N, the cross-over takes place at a small value of the aspect ratio W/N. From our numerical data we conclude that the three dimensional system has a finite non-diverging thermal conductivity and thus provide the first verification of Fourier's law in a system without pinning.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Linear response formula for finite frequency thermal conductance of open systems

An exact linear response expression is obtained for the heat current in a classical Hamiltonian system coupled to heat baths with time-dependent temperatures. The expression is equally valid at zero and finite frequencies. We present numerical results on the frequency dependence of the response function for three different one-dimensional models of coupled oscillators connected to Langevin baths with oscillating temperatures. For momentum conserving systems, a low frequency peak is seen that, is higher than the zero frequency response for large systems. For momentum non-conserving systems, there is no low frequency peak. The momentum non-conserving system is expected to satisfy Fourier's law, however, at the single bond level, we do not see any clear agreement with the predictions of the diffusion equation even at low frequencies. We also derive an exact analytical expression for the response of a chain of harmonic oscillators to a (not necessarily small) temperature difference; the agreement with the linear response simulation results for the same system is excellent.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stochastic pump of interacting particles

We consider the overdamped motion of Brownian particles, interacting via particle exclusion, in an external potential that varies with time and space. We show that periodic potentials that maintain specific position-dependent phase relations generate time-averaged directed current of particles. We obtain analytic results for a lattice version of the model using a recently developed perturbative approach. Many interesting features like particle-hole symmetry, current reversal with changing density, and system-size dependence of current are obtained. We propose possible experiments to test our predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Heat transport and phonon localization in mass-disordered harmonic crystals

We investigate the steady state heat current in two and three dimensional disordered harmonic crystals in a slab geometry, connected at the boundaries to stochastic white noise heat baths at different temperatures.The disorder causes short wavelength phonon modes to be localized so the heat current in this system is carried by the extended phonon modes which can be either diffusive or ballistic. Using ideas both from localization theory and from kinetic theory we estimate the contribution of various modes to the heat current and from this we obtain the asymptotic system size dependence of the current. These estimates are compared with results obtained from a numerical evaluation of an exact formula for the current, given in terms of a frequency transmission function, as well as from direct nonequilibrium simulations. These yield a strong dependence of the heat flux on boundary conditions. Our analytical arguments show that for realistic boundary conditions the conductivity is finite in three dimensions but we are not able to verify this numerically, except in the case where the system is subjected to an external pinning potential. This case is closely related to the problem of localization of electrons in a random potential and here we numerically verify that the pinned three dimensional system satisfies Fourier's law while the two dimensional system is a heat insulator. We also investigate the inverse participation ratio of different normal modes.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Fluctuation theorem in quantum heat conduction

We consider steady state heat conduction across a quantum harmonic chain connected to reservoirs modelled by infinite collection of oscillators. The heat, $Q$, flowing across the oscillator in a time interval $τ$ is a stochastic variable and we study the probability distribution function $P(Q)$. In the large $τ$ limit we use the formalism of full counting statistics (FCS) to compute the generating function of $P(Q)$ exactly. We show that $P(Q)$ satisfies the steady state fluctuation theorem (SSFT) regardless of the specifics of system, and it is nongaussian with clear exponential tails. The effect of finite $τ$ and nonlinearity is considered in the classical limit through Langevin simulations. We also obtain predictions of universal heat current fluctuations at low temperatures in clean wires.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Heat conduction in the disordered Fermi-Pasta-Ulam chain

We address the question of the effect of disorder on heat conduction in an anharmonic chain with interactions given by the Fermi-Pasta-Ulam (FPU) potential. In contrast to the conclusions of an earlier paper [Phys. Rev. Lett. 86, 63 (2001)] which found that disorder could induce a finite thermal conductivity at low temperatures, we find no evidence of a finite temperature transition in conducting properties. Instead, we find that at low temperatures, small system size transport properties are dominated by disorder but the asymptotic system size dependence of current is given by the usual FPU result J ~ 1/N^{2/3}. We also present new interesting results on the binary-mass ordered FPU chain.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Ergodicity properties of energy conserving single spin flip dynamics in the XY model

A single spin flip stochastic energy conserving dynamics for the XY model is considered. We study the ergodicity properties of the dynamics. It is shown that phase space trajectories densely fill the geometrically connected parts of the energy surface. We also show that while the dynamics is discrete and the phase point jumps around, it cannot make transitions between closed disconnected parts of the energy surface. Thus the number of distinct sectors depends on the number of geometrically disconnected parts of the energy surface. Information on the connectivity of the surfaces is obtained by studying the critical points of the energy function. We study in detail the case of two spins and find that the number of sectors can be either one or two, depending on the external fields and the energy. For a periodic lattice in $d$ dimensions, we find regions in phase space where the dynamics is non-ergodic and obtain a lower bound on the number of disconnected sectors. We provide some numerical evidence which suggests that such regions might be of small measure so that the dynamics is effectively ergodic.

0 citations0 reviewsNo code linkedOpen access