Researcher profile

Abhishek Dhar

Abhishek Dhar contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
17topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

20 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.

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.

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.

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&#39;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.

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.

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.

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&#39;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&#39;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.

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.

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.

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&#39; 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.

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&#39; 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.

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&#39;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.

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 &#34;active&#34; and &#34;passive&#34; 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.

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.

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&#39;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.

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.

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.

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.

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&#39;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&#39;s law. Here we provide the first numerical verification of Fourier&#39;s law in three dimensions. In the two dimensional pinned case we find that localization of phonon modes leads to a heat insulator.

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.