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

Deepak Dhar

Deepak 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-mechBiological Physicscond-mat.str-elphysics.soc-phq-fin.TRCell Behaviorcond-matcond-mat.dis-nncond-mat.softmath-phmath.MPmath.PRphysics.pop-ph

Catalog footprint

What is connected

41works
13topics
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

41 published item(s)

preprint2022arXiv

Finite-size scaling functions of the phase transition in the ferromagnetic Ising model on random regular graphs

We discuss the finite-size scaling of the ferromagnetic Ising model on random regular graphs. These graphs are locally tree-like, and in the limit of large graphs, the Bethe approximation gives the exact free energy per site. In the thermodynamic limit, the Ising model on these graphs show a phase transition. This transition is rounded off for finite graphs. We verify the scaling theory prediction that this rounding off is described in terms of the scaling variable $[T/T_c -1] S^{1/2}$ (where $T$ and $T_c$ are the temperature and the critical temperature respectively, and $S$ is the number of sites in the graph), and $\textit{not}$ in terms of a power of the diameter of the graph, which varies as $\log S$. We determine the theoretical scaling functions for the specific heat capacity and the magnetic susceptibility of the absolute value of the magnetization in closed form and compare them to Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The phase transition from nematic to high-density disordered phase in a system of hard rods on a lattice

A system of hard rigid rods of length $k$ on hypercubic lattices is known to undergo two phases transitions when chemical potential is increased: from a low density isotropic phase to an intermediate density nematic phase, and on further increase to a high-density phase with no orientational order. In this paper, we argue that, for large $k$, the second phase transition is a first order transition with a discontinuity in density in all dimensions greater than $1$. We show the chemical potential at the transition is $\approx k \ln [k /\ln k]$ for large $k$, and that the density of uncovered sites drops from a value $ \approx (\ln k)/k^2$ to a value of order $\exp(-ak)$, where $a$ is some constant, across the transition. We conjecture that these results are asymptotically exact, in all dimensions $d\geq 2$. We also present evidence of coexistence of nematic and disordered phases from Monte Carlo simulations for rods of length $9$ on the square lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Improved Upper Bounds on the Asymptotic Growth Velocity of Eden Clusters

We consider the asymptotic shape of clusters in the Eden model on a d-dimensional hypercubical lattice. We discuss two improvements for the well-known upper bound to the growth velocity in different directions by that of the independent branching process (IBP). In the IBP, each cell gives rise to a daughter cell at a neighboring site at a constant rate. In the first improvement, we do not allow such births along the bond connecting the cell to its mother cell. In the second, we iteratively evolve the system by a growth as IBP for a duration $Δ$ t, followed by culling process in which if any cell produced a descendant within this interval, who occupies the same site as the cell itself, then the descendant is removed. We study the improvement on the upper bound on the velocity for different dimensions d. The bounds are asymptotically exact in the large-d limit. But in $d =2$, the improvement over the IBP approximation is only a few percent.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

TASEP Speed Process: An Effective Medium Approach

We discuss the approximate phenomenological description of the motion of a single second-class particle in a two-species totally asymmetric simple exclusion process (TASEP) on a 1D lattice. Initially, the second class particle is located at the origin and to its left, all sites are occupied with first class particles while to its right, all sites are vacant. Ferrari and Kipnis proved that in any particular realization, the average velocity of the second class particle tends to a constant, but this mean value has a wide variation in different histories. We discuss this phenomena, here called the TASEP Speed Process, in an approximate effective medium description, in which the second class particle moves in a random background of the space-time dependent average density of the first class particles. We do this in three different approximations of increasing accuracy, treating the motion of the second-class particle first as a simple biassed random walk in a continuum Langevin equation, then as a biased Markovian random walk with space and time dependent jump rates, and finally as a Non-Markovian biassed walk with a non-exponential distribution of waiting times between jumps. We find that, when the displacement at time $T$ is $x_0$, the conditional expectation of displacement, at time $zT$ ($z>1$) is $zx_0$, and the variance of the displacement only varies as $z(z-1)T$. We extend this approach to describe the trajectories of a tagged particle in the case of a \emph{finite} lattice, where there are $L$ classes of particles on an $L$-site line, initially placed in the order of increasing class number. Lastly, we discuss a variant of the problem in which the exchanges between adjacent particles happened at rates proportional to the difference in their labels.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Directed Abelian sandpile with multiple downward neighbors

We study the directed Abelian sandpile model on a square lattice, with $K$ downward neighbors per site, $K > 2$. The $K=3$ case is solved exactly, which extends the earlier known solution for the $K=2$ case. For $K>2$, the avalanche clusters can have holes and side-branches and are thus qualitatively different from the $K=2$ case where avalance clusters are compact. However, we find that the critical exponents for $K>2$ are identical with those for the $K=2$ case, and the large scale structure of the avalanches for $K>2$ tends to the $K=2$ case.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stability of columnar order in assemblies of hard rectangles or squares

A system of $2\times d$ hard rectangles on square lattice is known to show four different phases for $d \geq 14$. As the covered area fraction $ρ$ is increased from $0$ to $1$, the system goes from low-density disordered phase, to orientationally-ordered nematic phase, to a columnar phase with orientational order and also broken translational invariance, to a high density phase in which orientational order is lost. For large d, the threshold density for the first transition $ρ_1^*$ tends to $0$, and the critical density for the third transition $ρ_3^*$ tends to $1$. Interestingly, simulations have shown that the critical density for the second transition $ρ_2^*$ tends to a non-trivial finite value $\approx 0.73$, as $d \rightarrow \infty$, and $ρ_2^* \approx 0.93$ for $d=2$. We provide a theoretical explanation of this interesting result. We develop an approximation scheme to calculate the surface tension between two differently ordered columnar phases. The density at which the surface tension vanishes gives an estimate $ρ_2^* = 0.746$, for $d\to \infty$, and $ρ_2^*=0.923$ for $d=2$. For all values of $d$, these estimates are in good agreement with Monte Carlo data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Many-agent limit of the Extreme Introvert-Extrovert model

We consider a toy model of interacting extrovert and introvert agents introduced earlier by Liu et al [Europhys. Lett. {\bf 100} (2012) 66007]. The number of extroverts, and introverts is $N$ each. At each time step, we select an agent at random, and allow her to modify her state. If an extrovert is selected, she adds a link at random to an unconnected introvert. If an introvert is selected, she removes one of her links. The set of $N^2$ links evolves in time, and may be considered as a set of Ising spins on an $N \times N$ square-grid with single-spin-flip dynamics. This dynamics satisfies detailed balance condition, and the probability of different spin configurations in the steady state can be determined exactly. The effective hamiltonian has long-range multi-spin couplings that depend on the row and column sums of spins. If the relative bias of choosing an extrovert over an introvert is varied, this system undergoes a phase transition from a state with very few links to one in which most links are occupied. We show that the behavior of the system can be determined exactly in the limit of large $N$. The behavior of large fluctuations in the total numer of links near the phase transition is determined. We also discuss two variations, called egalitarian and elitist agents, when the agents preferentially add or delete links to their least/ most-connected neighbor. These shows interesting cooperative behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Oslo model, hyperuniformity, and the quenched Edwards-Wilkinson model

We present simulations of the 1-dimensional Oslo rice pile model in which the critical height at each site is randomly reset after each toppling. We use the fact that the stationary state of this sandpile model is hyperuniform to reach system of sizes $> 10^7$. Most previous simulations were seriously flawed by important finite size corrections. We find that all critical exponents have values consistent with simple rationals: $ν=4/3$ for the correlation length exponent, $D =9/4$ for the fractal dimension of avalanche clusters, and $z=10/7 $ for the dynamical exponent. In addition we relate the hyperuniformity exponent to the correlation length exponent $ν$. Finally we discuss the relationship with the quenched Edwards-Wilkinson (qEW) model, where we find in particular that the local roughness exponent is $α_{\rm loc} = 1$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Columnar order and Ashkin-Teller criticality in mixtures of hard-squares and dimers

We show that critical exponents of the transition to columnar order in a {\em mixture} of $2 \times 1$ dimers and $2 \times 2$ hard-squares on the square lattice {\em depends on the composition of the mixture} in exactly the manner predicted by the theory of Ashkin-Teller criticality, including in the hard-square limit. This result settles the question regarding the nature of the transition in the hard-square lattice gas. It also provides the first example of a polydisperse system whose critical properties depend on composition. Our ideas also lead to some interesting predictions for a class of frustrated quantum magnets that exhibit columnar ordering of the bond-energies at low temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fragmentation of a sheet by propagating, branching and merging cracks

We consider a model of fragmentation of sheet by cracks that move with a velocity in preferred direction, but undergo random transverse displacements as they move. There is a non-zero probability of crack-splitting, and the split cracks move independently. If two cracks meet, they merge, and move as a single crack. In the steady state, there is non-zero density of cracks, and the sheet left behind by the moving cracks is broken into a large number of fragments of different sizes. The evolution operator for this model reduces to the Hamiltonian of quantum XY spin chain, which is exactly integrable. This allows us to determine the steady state, and also the distribution of sizes of fragments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Networks with preferred degree: A mini-review and some new results

Since their inception about a decade ago, dynamic networks which adapt to the state of the nodes have attracted much attention. One simple case of such an adaptive dynamics is a model of social networks in which individuals are typically comfortable with a certain number of contacts, i.e., preferred degrees. This paper is partly a review of earlier work of single homogeneous systems and ones with two interacting networks, and partly a presentation of some new results. In general, the dynamics does not obey detailed balance and the stationary distributions are not known analytically. A particular limit of the latter is a system of extreme introverts and extroverts - the XIE model. Remarkably, in this case, the detailed balance condition is satisfied, the exact distribution and an effective Hamiltonian can be found explicitly. Further, the model exhibits a phase transition in which the total number of links in the system - a macroscopically interesting quantity, displays an extreme Thouless effect. We show that in the limit of large populations and away from the transition, the model reduces to one with non-interacting agents of the majority subgroup. We determine the nature of fluctuations near the transition. We also introduce variants of the model where the agents show preferential attachment or detachment. There are significant changes to the degree distributions in the steady state, some of which can be understood by theoretical arguments and some remain to be explored. Many intriguing questions are posed, providing some food for thought and avenues for future research.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strategy switches and co-action equilibria in a minority game

We propose an analytically tractable variation of the minority game in which rational agents use probabilistic strategies. In our model, $N$ agents choose between two alternatives repeatedly, and those who are in the minority get a pay-off 1, others zero. The agents optimize the expectation value of their discounted future pay-off, the discount parameter being $λ$. We propose an alternative to the standard Nash equilibrium, called co-action equilibrium, which gives higher expected pay-off for all agents. The optimal choice of probabilities of different actions are determined exactly in terms of simple self -consistent equations. The optimal strategy is characterized by $N$ real parameters, which are non-analytic functions of $λ$, even for a finite number of agents. The solution for $N \leq 7$ is worked out explicitly indicating the structure of the solution for larger $N$. For large enough future time horizon, the optimal strategy switches from random choice to a win-stay lose-shift strategy, with the shift probability depending on the current state and $λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A class of exactly solved assisted hopping models of active-absorbing state transitions on a line

We construct a class of assisted hopping models in one dimension in which a particle can move only if it does not lie in an otherwise empty interval of length greater than $n+1$. We determine the exact steady state by a mapping to a gas of defects with only on-site interaction. We show that this system undergoes a phase transition as a function of the density $ρ$ of particles, from a low-density phase with all particles immobile for $ρ\le ρ_c = \frac{1}{n+1}$, to an active state for $ρ> ρ_c$. The mean fraction of movable particles in the active steady state varies as $(ρ- ρ_c)^β$, for $ρ$ near $ρ_c$. We show that for the model with range $n$, the exponent $β=n$, and thus can be made arbitrarily large.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A sandpile model for proportionate growth

An interesting feature of growth in animals is that different parts of the body grow at approximately the same rate. This property is called proportionate growth. In this paper, we review our recent work on patterns formed by adding $N$ grains at a single site in the abelian sandpile model. These simple models show very intricate patterns, show proportionate growth, and sometimes having a striking resemblance to natural forms. We give several examples of such patterns. We discuss the special cases where the asymptotic pattern can be determined exactly. The effect of noise in the background or in the rules on the patterns is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Memoryless nonlinear response: A simple mechanism for the 1/f noise

Discovering the mechanism underlying the ubiquity of $"1/f^α"$ noise has been a long--standing problem. The wide range of systems in which the fluctuations show the implied long--time correlations suggests the existence of some simple and general mechanism that is independent of the details of any specific system. We argue here that a {\it memoryless nonlinear response} suffices to explain the observed non--trivial values of $α$: a random input noisy signal $S(t)$ with a power spectrum varying as $1/f^{α'}$, when fed to an element with such a response function $R$ gives an output $R(S(t))$ that can have a power spectrum $1/f^α$ with $α< α'$. As an illustrative example, we show that an input Brownian noise ($α'=2$) acting on a device with a sigmoidal response function $R(S)= \sgn(S)|S|^x$, with $x<1$, produces an output with $α= 3/2 +x$, for $0 \leq x \leq 1/2$. Our discussion is easily extended to more general types of input noise as well as more general response functions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Proportionate growth in patterns formed in the rotor-router model

We study the growing patterns in the rotor-router model formed by adding $N$ walkers at the center of a $L \times L$ two-dimensional square lattice, starting with a periodic background of arrows, and relaxing to a stable configuration. The pattern is made of large number of triangular regions, where in each region all arrows point in the same direction. The square circumscribing the region, where all the arrows have been rotated atleast one full circle, may be considered as made up of smaller squares of different sizes, all of which grow linearly with $N$, for $ 1 \ll N < 2 L$. We use the Brooks-Smith-Stone-Tutte theorem relating tilings of squares by smaller squares to resistor networks, to determine the exact relative sizes of the different elements of the asymptotic pattern. We also determine the scaling limit of the function describing the variation of number of visits to a site with its position in the pattern. We also present evidence that deviations of the sizes of different small squares from the linear growth for finite $N$ are bounded and quasiperiodic functions of $N$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The nematic-disordered phase transition in systems of long rigid rods on two dimensional lattices

We study the phase transition from a nematic phase to a high-density disordered phase in systems of long rigid rods of length $k$ on the square and triangular lattices. We use an efficient Monte Carlo scheme that partly overcomes the problem of very large relaxation times of nearly jammed configurations. The existence of a continuous transition is observed on both lattices for $k=7$. We study correlations in the high-density disordered phase, and we find evidence of a crossover length scale $ξ^* \gtrsim 1400$, on the square lattice. For distances smaller than $ξ^*$, correlations appear to decay algebraically. Our best estimates of the critical exponents differ from those of the Ising model, but we cannot rule out a crossover to Ising universality class at length scales $\gg ξ^*$. On the triangular lattice, the critical exponents are consistent with those of the two dimensional three-state Potts universality class.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High-Activity Perturbation Expansion for the Hard Square Lattice Gas

We study a system of particles with nearest and next-nearest-neighbour exclusion on the square lattice (hard squares). This system undergoes a transition from a fluid phase at low density to a columnar ordered phase at high density. We develop a systematic high-activity perturbation expansion for the free energy per site about a state with perfect columnar order. We show that the different terms of the series can be regrouped to get a Mayer-like series for a polydisperse system of interacting vertical rods in which the $n$-th term is of order $z^{-(n+1)/2}$, where $z$ is the fugacity associated with each particle. We sum this series to get the exact expansion to order $1/z^{3/2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modelling proportionate growth

An important question in biology is how the relative size of different organs is kept nearly constant during growth of an animal. This property, called proportionate growth, has received increased attention in recent years. We discuss our recent work on a simple model where this feature comes out quite naturally from local rules, without fine tuning any parameter. The patterns produced are composed of large distinguishable structures with sharp boundaries, all of which grow at the same rate, keeping their overall shapes unchanged.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-perturbative corrections to mean-field behavior: spherical model on spider-web graph

We consider the spherical model on a spider-web graph. This graph is effectively infinite-dimensional, similar to the Bethe lattice, but has loops. We show that these lead to non-trivial corrections to the simple mean-field behavior. We first determine all normal modes of the coupled springs problem on this graph, using its large symmetry group. In the thermodynamic limit, the spectrum is a set of $δ$-functions, and all the modes are localized. The fractional number of modes with frequency less than $ω$ varies as $\exp (-C/ω)$ for $ω$ tending to zero, where $C$ is a constant. For an unbiased random walk on the vertices of this graph, this implies that the probability of return to the origin at time $t$ varies as $\exp(- C' t^{1/3})$, for large $t$, where $C'$ is a constant. For the spherical model, we show that while the critical exponents take the values expected from the mean-field theory, the free-energy per site at temperature $T$, near and above the critical temperature $T_c$, also has an essential singularity of the type $\exp[ -K {(T - T_c)}^{-1/2}]$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Power spectrum of mass and activity fluctuations in a sandpile

We consider a directed abelian sandpile on a strip of size $2\times n$, driven by adding a grain randomly at the left boundary after every $T$ time-steps. We establish the exact equivalence of the problem of mass fluctuations in the steady state and the number of zeroes in the ternary-base representation of the position of a random walker on a ring of size $3^n$. We find that while the fluctuations of mass have a power spectrum that varies as $1/f$ for frequencies in the range $ 3^{-2n} \ll f \ll 1/T$, the activity fluctuations in the same frequency range have a power spectrum that is linear in $f$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bethe approximation for a system of hard rigid rods: the random locally tree-like layered lattice

We study the Bethe approximation for a system of long rigid rods of fixed length k, with only excluded volume interaction. For large enough k, this system undergoes an isotropic-nematic phase transition as a function of density of the rods. The Bethe lattice, which is conventionally used to derive the self-consistent equations in the Bethe approximation, is not suitable for studying the hard-rods system, as it does not allow a dense packing of rods. We define a new lattice, called the random locally tree-like layered lattice, which allows a dense packing of rods, and for which the approximation is exact. We find that for a 4-coordinated lattice, k-mers with k>=4 undergo a continuous phase transition. For even coordination number q>=6, the transition exists only for k >= k_{min}(q), and is first order.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cooperation amongst competing agents in minority games

We study a variation of the minority game. There are N agents. Each has to choose between one of two alternatives everyday, and there is reward to each member of the smaller group. The agents cannot communicate with each other, but try to guess the choice others will make, based only the past history of number of people choosing the two alternatives. We describe a simple probabilistic strategy using which the agents acting independently, can still maximize the average number of people benefitting every day. The strategy leads to a very efficient utilization of resources, and the average deviation from the maximum possible can be made of order $(N^ε)$, for any $ε>0$. We also show that a single agent does not expect to gain by not following the strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Glassy dynamics and hysteresis in a linear system of orientable hard rods

We study the dynamics of a one-dimensional fluid of orientable hard rectangles with a non-coarse-grained microscopic mechanism of facilitation. The length occupied by a rectangle depends on its orientation, which is coupled to an external field. The equilibrium properties of our model are essentially those of the Tonks gas, but at high densities, the orientational degrees of freedom become effectively frozen due to jamming. This is a simple analytically tractable model of glassy phase. Under a cyclic variation of the pressure, hysteresis is observed. Following a pressure quench, the orientational persistence exhibits a two-stage decay characteristic of glassy systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pattern formation in fast-growing sandpiles

We study the patterns formed by adding $N$ sand-grains at a single site on an initial periodic background in the Abelian sandpile models, and relaxing the configuration. When the heights at all sites in the initial background are low enough, one gets patterns showing proportionate growth, with the diameter of the pattern formed growing as $N^{1/d}$ for large $N$, in $d$-dimensions. On the other hand, if sites with maximum stable height in the starting configuration form an infinite cluster, we get avalanches that do not stop. In this paper, we describe our unexpected finding of an interesting class of backgrounds in two dimensions, that show an intermediate behavior: For any $N$, the avalanches are finite, but the diameter of the pattern increases as $N^α$, for large $N$, with $1/2 < α\leq 1$. Different values of $α$ can be realized on different backgrounds, and the patterns still show proportionate growth. The non-compact nature of growth simplifies their analysis significantly. We characterize the asymptotic pattern exactly for one illustrative example with $α=1$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resonating valence bond wavefunctions and classical interacting dimer models

We relate properties of nearest-neighbour resonating valence bond (nnRVB) wavefunctions for $SU(g)$ spin systems on two dimensional bipartite lattices to those of fully-packed classical dimer models with potential energy $V$ on the same lattice. We define a cluster expansion of $V$ in terms of $n$-body potentials $V_n$, which are recursively determined from the nnRVB wavefunction on {\em finite subgraphs} of the original lattice. The magnitude of the $n$-body interaction $V_n$ ($n>1$) is of order ${\mathcal O}(g^{-(n-1)})$ for small $g^{-1}$, while $V_1$ reduces to a constant due to the fully-packed nature of the model. At leading non-trivial order on the square lattice, the interacting dimer model only has two-body interactions $V_2(g)$ that favour two parallel dimers on elementary plaquettes. Setting $g=2$ and using the results of earlier work on this interacting dimer model, we find that the long-distance behaviour of the bond-energy correlation function is dominated by an oscillatory term that decays as $ 1/|\vec{r}|^α$ with $α\approx 1.22$ for SU(2) spins. This result is in remarkable quantitative agreement with earlier direct numerical studies of the corresponding wavefunction, which give $α\approx 1.20$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Classical Heisenberg spins on a hexagonal lattice with Kitaev couplings

We analyse the low temperature properties of a system of classical Heisenberg spins on a hexagonal lattice with Kitaev couplings. For a lattice of 2N sites with periodic boundary conditions, we show that the ground states form an (N+1) dimensional manifold. We show that the ensemble of ground states is equivalent to that of a solid-on-solid model with continuously variable heights and nearest neighbour interactions, at a finite temperature. For temperature T tending to zero, all ground states have equal weight, and there is no order-by-disorder in this model. We argue that the bond-energy bond-energy correlations at distance R decay as 1/R^2 at zero temperature. This is verified by Monte Carlo simulations. We also discuss the relation to the quantum spin-S Kitaev model for large S, and obtain lower and upper bounds on the ground state energy of the quantum model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effect of Noise on Patterns Formed by Growing Sandpiles

We consider patterns generated by adding large number of sand grains at a single site in an abelian sandpile model with a periodic initial configuration, and relaxing. The patterns show proportionate growth. We study the robustness of these patterns against different types of noise, \textit{viz.}, randomness in the point of addition, disorder in the initial periodic configuration, and disorder in the connectivity of the underlying lattice. We find that the patterns show a varying degree of robustness to addition of a small amount of noise in each case. However, introducing stochasticity in the toppling rules seems to destroy the asymptotic patterns completely, even for a weak noise. We also discuss a variational formulation of the pattern selection problem in growing abelian sandpiles.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Restricted equilibrium ensembles: Exact equation of state of a model glass

We investigate the thermodynamic properties of a toy model of glasses: a hard-core lattice gas with nearest neighbor interaction in one dimension. The time-evolution is Markovian, with nearest-neighbor and next-nearest neighbor hoppings, and the transition rates are assumed to satisfy detailed balance condition, but the system is non-ergodic below a glass temperature. Below this temperature, the system is in restricted thermal equilibrium, where both the number of sectors, and the number of accessible states within a sector grow exponentially with the size of the system. Using partition functions that sum only over dynamically accessible states within a sector, and then taking a quenched average over the sectors, we determine the exact equation of state of this system.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spin-1 Kitaev model in one dimension

We study a one-dimensional version of the Kitaev model on a ring of size N, in which there is a spin S > 1/2 on each site and the Hamiltonian is J \sum_i S^x_i S^y_{i+1}. The cases where S is integer and half-odd-integer are qualitatively different. We show that there is a Z_2 valued conserved quantity W_n for each bond (n,n+1) of the system. For integer S, the Hilbert space can be decomposed into 2^N sectors, of unequal sizes. The number of states in most of the sectors grows as d^N, where d depends on the sector. The largest sector contains the ground state, and for this sector, for S=1, d =(\sqrt{5}+1)/2. We carry out exact diagonalization for small systems. The extrapolation of our results to large N indicates that the energy gap remains finite in this limit. In the ground state sector, the system can be mapped to a spin-1/2 model. We develop variational wave functions to study the lowest energy states in the ground state and other sectors. The first excited state of the system is the lowest energy state of a different sector and we estimate its excitation energy. We consider a more general Hamiltonian, adding a term λ\sum_n W_n, and show that this has gapless excitations in the range λ^c_1 \leq λ\leq λ^c_2. We use the variational wave functions to study how the ground state energy and the defect density vary near the two critical points λ^c_1 and λ^c_2.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Asymptotic shape of the region visited by an Eulerian Walker

We study an Eulerian walker on a square lattice, starting from an initially randomly oriented background using Monte Carlo simulations. We present evidence that, that, for large number of steps $N$, the asymptotic shape of the set of sites visited by the walker is a perfect circle. The radius of the circle increases as $N^{1/3}$, for large $N$, and the width of the boundary region grows as $N^{α/ 3}$, with $α= 0.40 \pm .05$. If we introduce stochasticity in the evolution rules, the mean square displacement of the walker, $<R_{N}^{2}> \sim N^{2ν}$, shows a crossover from the Eulerian ($ν= 1/3$) to a simple random walk ($ν=1/2$) behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Pattern Formation in Growing Sandpiles with Multiple Sources or Sinks

Adding sand grains at a single site in Abelian sandpile models produces beautiful but complex patterns. We study the effect of sink sites on such patterns. Sinks change the scaling of the diameter of the pattern with the number $N$ of sand grains added. For example, in two dimensions, in presence of a sink site, the diameter of the pattern grows as $\sqrt{(N/\log N)}$ for large $N$, whereas it grows as $\sqrt{N}$ if there are no sink sites. In presence of a line of sink sites, this rate reduces to $N^{1/3}$. We determine the growth rates for these sink geometries along with the case when there are two lines of sink sites forming a wedge, and its generalization to higher dimensions. We characterize one such asymptotic patterns on the two-dimensional F-lattice with a single source adjacent to a line of sink sites, in terms of position of different spatial features in the pattern. For this lattice, we also provide an exact characterization of the pattern with two sources, when the line joining them is along one of the axes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Scaling relation for determining the critical threshold for continuum percolation of overlapping discs of two sizes

We study continuum percolation of overlapping circular discs of two sizes. We propose a phenomenological scaling equation for the increase in the effective size of the larger discs due to the presence of the smaller discs. The critical percolation threshold as a function of the ratio of sizes of discs, for different values of the relative areal densities of two discs, can be described in terms of a scaling function of only one variable. The recent accurate Monte Carlo estimates of critical threshold by Quintanilla and Ziff [Phys. Rev. E, 76 051115 (2007)] are in very good agreement with the proposed scaling relation.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Pattern formation in growing sandpiles

Adding grains at a single site on a flat substrate in the Abelian sandpile models produce beautiful complex patterns. We study in detail the pattern produced by adding grains on a two-dimensional square lattice with directed edges (each site has two arrows directed inward and two outward), starting with a periodic background with half the sites occupied. The size of the pattern formed scales with the number of grains added $N$ as $\sqrt{N}$. We give exact characterization of the asymptotic pattern, in terms of the position and shape of different features of the pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Steady state of Stochastic Sandpile Models

We study the steady state of the abelian sandpile models with stochastic toppling rules. The particle addition operators commute with each other, but in general these operators need not be diagonalizable. We use their abelian algebra to determine their eigenvalues, and the Jordan block structure. These are then used to determine the probability of different configurations in the steady state. We illustrate this procedure by explicitly determining the numerically exact steady state for a one dimensional example, for systems of size $\le12$, and also study the density profile in the steady state.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Critical Behavior of Sandpile Models with Sticky Grains

We revisit the question whether the critical behavior of sandpile models with sticky grains is in the directed percolation universality class. Our earlier theoretical arguments in favor, supported by evidence from numerical simulations [ Phys. Rev. Lett., {\bf 89} (2002) 104303], have been disputed by Bonachela et al. [Phys. Rev. E {\bf 74} (2004) 050102] for sandpiles with no preferred direction. We discuss possible reasons for the discrepancy. Our new results of longer simulations of the one-dimensional undirected model fully support our earlier conclusions.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Emergence of quasi-units in the one dimensional Zhang model

We study the Zhang model of sandpile on a one dimensional chain of length $L$, where a random amount of energy is added at a randomly chosen site at each time step. We show that in spite of this randomness in the input energy, the probability distribution function of energy at a site in the steady state is sharply peaked, and the width of the peak decreases as $ {L}^{-1/2}$ for large $L$. We discuss how the energy added at one time is distributed among different sites by topplings with time. We relate this distribution to the time-dependent probability distribution of the position of a marked grain in the one dimensional Abelian model with discrete heights. We argue that in the large $L$ limit, the variance of energy at site $x$ has a scaling form $L^{-1}g(x/L)$, where $g(ξ)$ varies as $\log(1/ξ)$ for small $ξ$, which agrees very well with the results from numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

On the orientational ordering of long rods on a lattice

We argue that a system of straight rigid rods of length k on square lattice with only hard-core interactions shows two phase transitions as a function of density, rho, for k >= 7. The system undergoes a phase transition from the low-density disordered phase to a nematic phase as rho is increased from 0, at rho = rho_c1, and then again undergoes a reentrant phase transition from the nematic phase to a disordered phase at rho = rho_c2 < 1.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Drift and trapping in biased diffusion on disordered lattices

We reexamine the theory of transition from drift to no-drift in biased diffusion on percolation networks. We argue that for the bias field B equal to the critical value B_c, the average velocity at large times t decreases to zero as 1/log(t). For B < B_c, the time required to reach the steady-state velocity diverges as exp(const/|B_c-B|). We propose an extrapolation form that describes the behavior of average velocity as a function of time at intermediate time scales. This form is found to have a very good agreement with the results of extensive Monte Carlo simulations on a 3-dimensional site-percolation network and moderate bias.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Breakdown of Simple Scaling in Abelian Sandpile Models in One Dimension

We study the abelian sandpile model on decorated one dimensional chains. We determine the structure and the asymptotic form of distribution of avalanche-sizes in these models, and show that these differ qualitatively from the behavior on a simple linear chain. We find that the probability distribution of the total number of topplings $s$ on a finite system of size $L$ is not described by a simple finite size scaling form, but by a linear combination of two simple scaling forms $Prob_L(s) = 1/L f_1(s/L) + 1/L^2 f_2(s/L^2)$, for large $L$, where $f_1$ and $f_2$ are some scaling functions of one argument.

0 citations0 reviewsNo code linkedOpen access