Source author record

Kabir Ramola

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

17 published item(s)

preprint2023arXiv

Enhanced Vibrational Stability in Glass Droplets

We show through simulations of amorphous solids prepared in open boundary conditions that they possess significantly fewer low-frequency vibrational modes compared to their periodic boundary counterparts. Specifically, using measurements of the vibrational density of states, we find that the $D(ω) \sim ω^4$ law changes to $D(ω) \sim ω^δ$ with $δ\approx 5$ in two dimensions and $δ\approx 4.5$ in three dimensions. Crucially, this enhanced stability is achieved when utilizing slow annealing protocols to generate solid configurations. We perform an anharmonic analysis of the minima corresponding to the lowest-frequency modes in such open-boundary systems and discuss their correlation with the density of states. A study of various system sizes further reveals that small systems display a higher degree of localization in vibrations. Lastly, we confine open-boundary solids in order to introduce macroscopic stresses in the system which are absent in the unconfined system, and find that the $D(ω) \sim ω^4$ behavior is recovered.

preprint2022arXiv

First contact breaking distributions in strained disordered crystals

We derive exact probability distributions for the strain ($ε$) at which the first stress drop event occurs in uniformly strained disordered crystals, with quenched disorder introduced through polydispersity in particle sizes. We characterize these first stress drop events numerically as well as theoretically, and identify them with the first contact breaking event in the system. Our theoretical results are corroborated with numerical simulations of quasistatic volumetric strain applied to disordered near-crystalline configurations of athermal soft particles. We develop a general technique to determine the $distribution$ of strains at which the first stress drop events occur, through an exact mapping between the cumulative distribution of first contact breaking events and the volume of a convex polytope whose dimension is determined by the number of defects $N_d$ in the system. An exact numerical computation of this polytope volume for systems with small numbers of defects displays a remarkable match with the distribution of strains generated through direct numerical simulations. Finally, we derive the distribution of strains at which the first stress drop occurs, assuming that individual contact breaking events are uncorrelated, which accurately reproduces distributions obtained from direct numerical simulations.

preprint2022arXiv

Tensor Electromagnetism and Emergent Elasticity in Jammed Solids

The theory of mechanical response and stress transmission in disordered, jammed solids poses several open questions of how non-periodic networks -- apparently indistinguishable from a snapshot of a fluid -- sustain shear. We present a stress-only theory of emergent elasticity for a non-thermal amorphous assembly of grains in a jammed solid, where each grain is subjected to mechanical constraints of force and torque balance. These grain-level constraints lead to the Gauss's law of an emergent $U(1)$ tensor electromagnetism, which then accounts for the mechanical response of such solids. This formulation of amorphous elasticity has several immediate consequences. The mechanical response maps exactly to the static, dielectric response of this tensorial electromagnetism with the polarizability of the medium mapping to emergent elastic moduli. External forces act as vector electric charges whereas the tensorial magnetic fields are sourced by momentum density. The dynamics in the electric and magnetic sectors, naturally translate into the dynamics of the rigid jammed network and ballistic particle motion respectively. The theoretical predictions for both stress-stress correlations and responses are borne out by the results of numerical simulations of frictionless granular packings in the static limit of the theory in both 2D and 3D.

preprint2021arXiv

Universal non-Debye low-frequency vibrations in sheared amorphous solids

We study energy minimized configurations of amorphous solids with a simple shear degree of freedom. We show that the low-frequency regime of the vibrational density of states of structural glass formers is crucially sensitive to the stress-ensemble from which the configurations are sampled. In both two and three dimensions, a shear-stabilized ensemble displays a $D(ω_{\min}) \sim ω^{5}_{\min}$ regime, as opposed to the $ω^{4}_{\min}$ regime observed under unstrained conditions. We also study an ensemble of two dimensional, strained amorphous solids near a plastic event. We show that the minimum eigenvalue distribution at a strain $γ$ near the plastic event occurring at $γ_{P}$, displays a collapse when scaled by $\sqrt{γ_P - γ}$, and with the number of particles as $N^{-0.22}$. Notably, at low-frequencies, this scaled distribution displays a robust $D(ω_{\min}) \sim ω^{6}_{\min}$ power-law regime, which survives in the large $N$ limit. Finally, we probe the universal properties of this ensemble through a characterization of the second and third eigenvalues of the Hessian matrix near a plastic event.

preprint2020arXiv

Arrested States in Persistent Active Matter: Gelation without Attraction

We explore phase separation and kinetic arrest in a model active colloidal system consisting of self-propelled, hard-core particles with nonconvex shapes. The passive limit of the model, namely cross-shaped particles on a square lattice, exhibits a first-order transition from a fluid phase to a solid phase with increasing density. Quenches into the two-phase coexistence region exhibit an aging regime. The nonconvex shape of the particles eases jamming in the passive system and leads to strong inhibition of rotations of the active particles. Using numerical simulations and analytical modeling, we quantify the nonequilibrium phase behavior as a function of density and activity. If we view activity as the analog of attraction strength, the phase diagram exhibits strong similarities to that of attractive colloids, exhibiting both aging, glassy states and gel-like arrested states. The two types of dynamically arrested states, glasses and gels, are distinguished by the appearance of density heterogenities in the latter. In the infinitely persistent limit, we show that a coarse-grained model based on the asymmetric exclusion process quantitatively predicts the density profiles of the gel states. The predictions remain qualitatively valid for finite rotation rates. Using these results, we classify the activity-driven phases and identify the boundaries separating them.

preprint2020arXiv

Athermal fluctuations in disordered crystals

We analyze the fluctuations in particle positions and inter-particle forces in disordered jammed crystals in the limit of weak disorder. We demonstrate that such athermal systems are fundamentally different from their thermal counterparts, characterized by constrained fluctuations of forces perpendicular to the lattice directions. We develop a disorder perturbation expansion in polydispersity about the crystalline state, which we use to derive exact results to linear order. We show that constrained fluctuations result as a consequence of local force balance conditions, and are characterized by non-Gaussian distributions which we derive exactly. We analytically predict several properties of such systems, including the scaling of the average coordination with polydispersity and packing fraction, which we verify with numerical simulations using soft disks with one-sided harmonic interactions.

preprint2020arXiv

Timescale divergence at the shear jamming transition

We find that in simulations of quasi-statically sheared frictional disks, the shear jamming transition can be characterized by an abrupt jump in the number of force bearing contacts between particles. This mechanical coordination number increases discontinuously from $Z = 0$ to $Z \gtrsim d +1$ at a critical shear value $γ_c$, as opposed to a smooth increase in the number of geometric contacts. This is accompanied by a diverging timescale $τ^*$ that characterizes the time required by the system to attain force balance when subjected to a perturbation. As the global shear $γ$ approaches the critical value $γ_c$ from below, one observes the divergence of the time taken to relax to a state where all the inter-particle contacts have uniformly zero force. Above $γ_{c}$, the system settles into a state characterized by finite forces between particles, with the timescale also increasing as $γ\to γ_{c}^{+}$. By using two different protocols to generate force balanced configurations, we show that this timescale divergence is a robust feature that accompanies the shear jamming transition.

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.

preprint2015arXiv

Spatial Extent of Branching Brownian Motion

We study the one dimensional branching Brownian motion starting at the origin and investigate the correlation between the rightmost ($X_{\max}\geq 0$) and leftmost ($X_{\min} \leq 0$) visited sites up to time $t$. At each time step the existing particles in the system either diffuse (with diffusion constant $D$), die (with rate $a$) or split into two particles (with rate $b$). We focus on the regime $b \leq a$ where these two extreme values $X_{\max}$ and $X_{\min}$ are strongly correlated. We show that at large time $t$, the joint probability distribution function (PDF) of the two extreme points becomes stationary $P(X,Y,t \to \infty) \to p(X,Y)$. Our exact results for $p(X,Y)$ demonstrate that the correlation between $X_{\max}$ and $X_{\min}$ is nonzero, even in the stationary state. From this joint PDF, we compute exactly the stationary PDF $p(ζ)$ of the (dimensionless) span $ζ= {(X_{\max} - X_{\min})}/{\sqrt{D/b}}$, which is the distance between the rightmost and leftmost visited sites. This span distribution is characterized by a linear behavior ${p}(ζ) \sim \frac{1}{2} \left(1 + Δ\right) ζ$ for small spans, with $Δ= \left(\frac{a}{b} -1\right)$. In the critical case ($Δ= 0$) this distribution has a non-trivial power law tail ${p}(ζ) \sim 8 π\sqrt{3} /ζ^3$ for large spans. On the other hand, in the subcritical case ($Δ> 0$), we show that the span distribution decays exponentially as ${p}(ζ) \sim (A^2/2) ζ\exp \left(- \sqrtΔ~ζ\right)$ for large spans, where $A$ is a non-trivial function of $Δ$ which we compute exactly. We show that these asymptotic behaviors carry the signatures of the correlation between $X_{\max}$ and $X_{\min}$. Finally we verify our results via direct Monte Carlo simulations.

preprint2014arXiv

Branching Brownian Motion Conditioned on Particle Numbers

We study analytically the order and gap statistics of particles at time $t$ for the one dimensional branching Brownian motion, conditioned to have a fixed number of particles at $t$. The dynamics of the process proceeds in continuous time where at each time step, every particle in the system either diffuses (with diffusion constant $D$), dies (with rate $d$) or splits into two independent particles (with rate $b$). We derive exact results for the probability distribution function of $g_k(t) = x_k(t) - x_{k+1}(t)$, the distance between successive particles, conditioned on the event that there are exactly $n$ particles in the system at a given time $t$. We show that at large times these conditional distributions become stationary $P(g_k, t \to \infty|n) = p(g_k|n)$. We show that they are characterised by an exponential tail $p(g_k|n) \sim \exp[-\sqrt{\frac{|b - d|}{2 D}} ~g_k]$ for large gaps in the subcritical ($b < d$) and supercritical ($b > d$) phases, and a power law tail $p(g_k) \sim 8\left(\frac{D}{b}\right){g_k}^{-3}$ at the critical point ($b = d$), independently of $n$ and $k$. Some of these results for the critical case were announced in a recent letter [K. Ramola, S. N. Majumdar and G. Schehr, Phys. Rev. Lett. 112, 210602 (2014)].

preprint2014arXiv

Fluctuations of random matrix products and 1D Dirac equation with random mass

We study the fluctuations of certain random matrix products $Π_N=M_N\cdots M_2M_1$ of $\mathrm{SL}(2,\mathbb{R})$, describing localisation properties of the one-dimensional Dirac equation with random mass. In the continuum limit, i.e. when matrices $M_n$'s are close to the identity matrix, we obtain convenient integral representations for the variance $Γ_2=\lim_{N\to\infty}\mathrm{Var}(\ln||Π_N||)/N$. The case studied exhibits a saturation of the variance at low energy $\varepsilon$ along with a vanishing Lyapunov exponent $Γ_1=\lim_{N\to\infty}\ln||Π_N||/N$, leading to the behaviour $Γ_2/Γ_1\sim\ln(1/|\varepsilon|)\to\infty$ as $\varepsilon\to0$. Our continuum description sheds new light on the Kappus-Wegner (band center) anomaly.

preprint2014arXiv

Universal Order and Gap Statistics of Critical Branching Brownian Motion

We study the order statistics of one dimensional branching Brownian motion in which particles either diffuse (with diffusion constant $D$), die (with rate $d$) or split into two particles (with rate $b$). At the critical point $b=d$ which we focus on, we show that, at large time $t$, the particles are collectively bunched together. We find indeed that there are two length scales in the system: (i) the diffusive length scale $\sim \sqrt{Dt}$ which controls the collective fluctuations of the whole bunch and (ii) the length scale of the gap between the bunched particles $\sim \sqrt{D/b}$. We compute the probability distribution function $P(g_k,t|n)$ of the $k$th gap $g_k = x_k - x_{k+1}$ between the $k$th and $(k+1)$th particles given that the system contains exactly $n>k$ particles at time $t$. We show that at large $t$, it converges to a stationary distribution $P(g_k,t\to \infty|n) = p(g_k|n)$ with an algebraic tail $p(g_k|n) \sim 8(D/b) g_k^{-3}$, for $g_k \gg 1$, independent of $k$ and $n$. We verify our predictions with Monte Carlo simulations.

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}$.

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

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.

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.

preprint2008arXiv

Fragment Formation in Biased Random Walks

We analyse a biased random walk on a 1D lattice with unequal step lengths. Such a walk was recently shown to undergo a phase transition from a state containing a single connected cluster of visited sites to one with several clusters of visited sites (fragments) separated by unvisited sites at a critical probability p_c, [PRL 99, 180602 (2007)]. The behaviour of rho(l), the probability of formation of fragments of length l is analysed. An exact expression for the generating function of rho(l) at the critical point is derived. We prove that the asymptotic behaviour is of the form rho(l) ~ 3/[l(log l)^2].