Source author record

Arnab Sen

Arnab Sen 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

40works
18topics
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

40 published item(s)

preprint2024arXiv

Weak universality, quantum many-body scars and anomalous infinite-temperature autocorrelations in a one-dimensional spin model with duality

We study a one-dimensional spin-$1/2$ model with three-spin interactions and a transverse magnetic field $h$. The model has a $Z_2 \times Z_2$ symmetry, and a duality between $h$ and $1/h$. The self-dual point at $h=1$ is a quantum critical point with a continuous phase transition. We compute the critical exponents $z$, $β$, $γ$ and $ν$, and the central charge $c$ numerically using exact diagonalization (ED) for systems with periodic boundary conditions. We find that both $z$ and $c$ are equal to $1$, implying that the critical point is governed by a conformal field theory. The values obtained for $β/ν$, $γ/ν$, and $ν$ from ED suggest that the model exhibits Ashkin-Teller criticality with an effective coupling that is intermediate between the four-state Potts model and two decoupled transverse field Ising models. An analysis on larger systems but with open boundaries using density-matrix renormalization group calculations, however, shows that the self-dual point may be in the same universality class as the four-state Potts model. An energy level spacing analysis shows that the model is not integrable. For a system with periodic boundary conditions, there are an exponentially large number of exact mid-spectrum zero-energy eigenstates. A subset of these eigenstates have wave functions which are independent of $h$ and have unusual entanglement structure, suggesting that they are quantum many-body scars. The number of such states scales at least linearly with system size. Finally, we study the infinite-temperature autocorrelation functions close to one end of an open system. We find that some of the autocorrelators relax anomalously in time, with pronounced oscillations and very small decay rates if $h \gg 1$ or $h \ll 1$. If $h$ is close to the critical point, the autocorrelators decay quickly to zero except for an autocorrelator at the end site.

preprint2022arXiv

Quantum order-by-disorder induced phase transition in Rydberg ladders with staggered detuning

$^{87}{\rm Rb}$ atoms are known to have long-lived Rydberg excited states with controllable excitation amplitude (detuning) and strong repulsive van der Waals interaction $V_{{\bf r} {\bf r'}}$ between excited atoms at sites ${\bf r}$ and ${\bf r'}$. Here we study such atoms in a two-leg ladder geometry in the presence of both staggered and uniform detuning with amplitudes $Δ$ and $λ$ respectively. We show that when $V_{\bf r r'} \gg(\ll) Δ, λ$ for $|{\bf r}-{\bf r'}|=1(>1)$, these ladders host a plateau for a wide range of $λ/Δ$ where the ground states are selected by a quantum order-by-disorder mechanism from a macroscopically degenerate manifold of Fock states with fixed Rydberg excitation density $1/4$. Our study further unravels the presence of an emergent Ising transition stabilized via the order-by-disorder mechanism inside the plateau. We identify the competing terms responsible for the transition and estimate a critical detuning $λ_c/Δ=1/3$ which agrees well with exact-diagonalization based numerical studies. We also study the fate of this transition for a realistic interaction potential $V_{{\bf r} {\bf r'}} = V_0 /|{\bf r}-{\bf r'}|^6$, demonstrate that it survives for a wide range of $V_0$, and provide analytic estimate of $λ_c$ as a function of $V_0$. This allows for the possibility of a direct verification of this transition in standard experiments which we discuss.

preprint2022arXiv

Scars from protected zero modes and beyond in $U(1)$ quantum link and quantum dimer models

We demonstrate the presence of anomalous high-energy eigenstates, or many-body scars, in $U(1)$ quantum link and quantum dimer models on square and rectangular lattices. In particular, we consider the paradigmatic Rokhsar-Kivelson Hamiltonian $H=\mathcal{O}_{\mathrm{kin}} + λ\mathcal{O}_{\mathrm{pot}}$ where $\mathcal{O}_{\mathrm{pot}}$ ($\mathcal{O}_{\mathrm{kin}}$) is defined as a sum of terms on elementary plaquettes that are diagonal (off-diagonal) in the computational basis. Both these interacting models possess an exponentially large number of mid-spectrum zero modes in system size at $λ=0$ that are protected by an index theorem preventing any mixing with the nonzero modes at this coupling. We classify different types of scars for $|λ| \lesssim \mathcal{O}(1)$ both at zero and finite winding number sectors complementing and significantly generalizing our previous work [Banerjee and Sen, Phys. Rev. Lett. 126, 220601 (2021)]. The scars at finite $λ$ show a rich variety with those that are composed solely from the zero modes of $\mathcal{O}_{\mathrm{kin}}$, those that contain an admixture of both the zero and the nonzero modes of $\mathcal{O}_{\mathrm{kin}}$, and finally those composed solely from the nonzero modes of $\mathcal{O}_{\mathrm{kin}}$. We give analytic expressions for certain "lego scars" for the quantum dimer model on rectangular lattices where one of the linear dimensions can be made arbitrarily large, with the building blocks (legos) being composed of emergent singlets and other more complicated entangled structures.

preprint2021arXiv

Coincident angle-resolved state-selective photoelectron spectroscopy of acetylene molecules: a candidate system for time-resolved dynamics

The acetylene-vinylidene system serves as a benchmark for investigations of ultrafast dynamical processes where the coupling of the electronic and nuclear degrees of freedom provides a fertile playground to explore the femto- and sub-femto-second physics with coherent extreme-ultraviolet (EUV) photon sources both on the table-top as well as free-electron lasers. We focus on detailed investigations of this molecular system in the photon energy range $19...40$ eV where EUV pulses can probe the dynamics effectively. We employ photoelectron-photoion coincidence (PEPICO) spectroscopy to uncover hitherto unrevealed aspects of this system. In this work, the role of excited states of the $C_{2}H_{2}^{+}$ cation, the primary photoion, is specifically addressed. From photoelectron energy spectra and angular distributions, the nature of the dissociation and isomerization channels is discerned. Exploiting the $4π$-collection geometry of velocity map imaging spectrometer, we not only probe pathways where the efficiency of photoionization is inherently high but also perform PEPICO spectroscopy on relatively weak channels.

preprint2021arXiv

Dynamical relaxation of correlators in periodically driven integrable quantum systems

We show that the correlation functions of a class of periodically driven integrable closed quantum systems approach their steady state value as $n^{-(α+1)/β}$, where $n$ is the number of drive cycles and $α$ and $β$ denote positive integers. We find that generically $β=2$ within a dynamical phase characterized by a fixed $α$; however, its value can change to $β=3$ or $β=4$ either at critical drive frequencies separating two dynamical phases or at special points within a phase. We show that such decays are realized in both driven Su-Schrieffer-Heeger (SSH) and one-dimensional (1D) transverse field Ising models, discuss the role of symmetries of the Floquet spectrum in determining $β$, and chart out the values of $α$ and $β$ realized in these models. We analyze the SSH model for a continuous drive protocol using a Floquet perturbation theory which provides analytical insight into the behavior of the correlation functions in terms of its Floquet Hamiltonian. This is supplemented by an exact numerical study of a similar behavior for the 1D Ising model driven by a square pulse protocol. For both models, we find a crossover timescale $n_c$ which diverges at the transition. We also unravel a long-time oscillatory behavior of the correlators when the critical drive frequency, $ω_c$, is approached from below ($ω< ω_c$). We tie such behavior to the presence of multiple stationary points in the Floquet spectrum of these models and provide an analytic expression for the time period of these oscillations.

preprint2021arXiv

Periodically driven Rydberg chains with staggered detuning

We study the stroboscopic dynamics of a periodically driven finite Rydberg chain with staggered ($Δ$) and time-dependent uniform ($λ(t)$) detuning terms using exact diagonalization (ED). We show that at intermediate drive frequencies ($ω_D$), the presence of a finite $Δ$ results in violation of the eigenstate thermalization hypothesis (ETH) via clustering of Floquet eigenstates. Such clustering is lost at special commensurate drive frequencies for which $\hbar ω_d=n Δ$ ($n \in Z$) leading to restoration of ergodicity. The violation of ETH in these driven finite-sized chains is also evident from the dynamical freezing displayed by the density-density correlation function at specific $ω_D$. Such a correlator exhibits stable oscillations with perfect revivals when driven close to the freezing frequencies for initial all spin-down ($|0\rangle$) or Neel ($|{\mathbb Z}_2\rangle$, with up-spins on even sites) states. The amplitudes of these oscillations vanish at the freezing frequencies and reduces upon increasing $Δ$; their frequencies, however, remains pinned to $Δ/\hbar$ in the large $Δ$ limit. In contrast, for the $|{\bar {\mathbb Z}_2}\rangle$ (time-reversed partner of $|{\mathbb Z}_2\rangle$) initial state, we find complete absence of such oscillations leading to freezing for a range of $ω_D$; this range increases with $Δ$. We also study the properties of quantum many-body scars in the Floquet spectrum of the model as a function of $Δ$ and show the existence of novel mid-spectrum scars at large $Δ$. We supplement our numerical results with those from an analytic Floquet Hamiltonian computed using Floquet perturbation theory (FPT) and also provide a semi-analytic computation of the quantum scar states within a forward scattering approximation (FSA).

preprint2021arXiv

Signatures of multifractality in a periodically driven interacting Aubry-André model

We study the many-body localization (MBL) transition of Floquet eigenstates in a driven, interacting fermionic chain with an incommensurate Aubry-André potential and a time-periodic hopping amplitude as a function of the drive frequency $ω_D$ using exact diagonalization (ED). We find that the nature of the Floquet eigenstates change from ergodic to Floquet-MBL with increasing frequency; moreover, for a significant range of intermediate $ω_D$, the Floquet eigenstates exhibit non-trivial fractal dimensions. We find a possible transition from the ergodic to this multifractal phase followed by a gradual crossover to the MBL phase as the drive frequency is increased. We also study the fermion auto-correlation function, entanglement entropy, normalized participation ratio (NPR), fermion transport and the inverse participation ratio (IPR) as a function of $ω_D$. We show that the auto-correlation, fermion transport and NPR displays qualitatively different characteristics (compared to their behavior in the ergodic and MBL regions) for the range of $ω_D$ which supports multifractal eigenstates. In contrast, the entanglement growth in this frequency range tend to have similar features as in the MBL regime; its rate of growth is controlled by $ω_D$. Our analysis thus indicates that the multifractal nature of Floquet-MBL eigenstates can be detected by studying auto-correlation function and fermionic transport of these driven chains. We support our numerical results with a semi-analytic expression of the Floquet Hamiltonian obtained using Floquet perturbation theory (FPT) and discuss possible experiments which can test our predictions.

preprint2020arXiv

Dynamics of the vacuum state in a periodically driven Rydberg chain

We study the dynamics of the periodically driven Rydberg chain starting from the state with zero Rydberg excitations (vacuum state denoted by $|0\rangle$) using a square pulse protocol in the high drive amplitude limit. We show, using exact diagonalization for finite system sizes ($L\le 26$), that the Floquet Hamiltonian of the system, within a range of drive frequencies which we chart out, hosts a set of quantum scars which have large overlap with the $|0\rangle$ state. These scars are distinct from their counterparts having high overlap with the maximal Rydberg excitation state ($|\mathbb{Z}_2\rangle$); they coexist with the latter class of scars and lead to persistent coherent oscillations of the density-density correlator starting from the $|0\rangle$ state. We also identify special drive frequencies at which the system undergoes perfect dynamic freezing and provide an analytic explanation for this phenomenon. Finally, we demonstrate that for a wide range of drive frequencies, the system reaches a steady state with sub-thermal values of the density-density correlator. The presence of such sub-thermal steady states, which are absent for dynamics starting from the $|\mathbb{Z}_2\rangle$ state, imply a weak violation of the eigenstate thermalization hypothesis in finite sized Rydberg chains distinct from that due to the scar-induced persistent oscillations reported earlier. We conjecture that in the thermodynamic limit such states may exist as pre-thermal steady states that show anomalously slow relaxation. We supplement our numerical results by deriving an analytic expression for the Floquet Hamiltonian using a Floquet perturbation theory in the high amplitude limit which provides an analytic, albeit qualitative, understanding of these phenomena at arbitrary drive frequencies. We discuss experiments which can test our theory.

preprint2020arXiv

Entanglement measures and non-equilibrium dynamics of quantum many-body systems: a path integral approach

We present a path integral formalism for expressing matrix elements of the density matrix of a quantum many-body system between any two coherent states in terms of standard Matsubara action with periodic(anti-periodic) boundary conditions on bosonic(fermionic) fields. We show that this enables us to express several entanglement measures for bosonic/fermionic many-body systems described by a Gaussian action in terms of the Matsubara Green function. We apply this formalism to compute various entanglement measures for the two-dimensional Bose-Hubbard model in the strong-coupling regime, both in the presence and absence of Abelian and non-Abelian synthetic gauge fields, within a strong coupling mean-field theory. In addition, our method provides an alternative formalism for addressing time evolution of quantum-many body systems, with Gaussian actions, driven out of equilibrium without the use of Keldysh technique. We demonstrate this by deriving analytical expressions of the return probability and the counting statistics of several operators for a class of integrable models represented by free Dirac fermions subjected to a periodic drive in terms of the elements of their Floquet Hamiltonians. We provide a detailed comparison of our method with the earlier, related, techniques used for similar computations, discuss the significance of our results, and chart out other systems where our formalism can be used.

preprint2020arXiv

Restoring coherence via aperiodic drives in a many-body quantum system

We study the unitary dynamics of randomly or quasi-periodically driven tilted Bose-Hubbard (tBH) model in one dimension deep inside its Mott phase starting from a $\mathbb{Z}_2$ symmetry-broken state. The randomness is implemented via a telegraph noise protocol in the drive period while the quasi-periodic drive is chosen to correspond to a Thue-Morse sequence. The periodically driven tBH model (with a square pulse protocol characterized by a time period $T$) is known to exhibit transitions from dynamical regimes with long-time coherent oscillations to those with rapid thermalization. Here we show that starting from a regime where the periodic drive leads to rapid thermalization, a random drive, which consists of a random sequence of square pulses with period $T+αdT$, where $α=\pm 1$ is a random number and $dT$ is the amplitude of the noise, restores long-time coherent oscillations for special values of $dT$. A similar phenomenon can be seen for a quasi-periodic drive following a Thue-Morse sequence where such coherent behavior is shown to occur for a larger number of points in the $(T, dT)$ plane due to the additional structure of the drive protocol. We chart out the dynamics of the system in the presence of such aperiodic drives, provide a qualitative analytical understanding of this phenomenon, point out the role of quantum scars behind it, and discuss experiments which can test our theory.

preprint2019arXiv

Collapse and revival of quantum many-body scars via Floquet engineering

The presence of quantum scars, athermal eigenstates of a many-body Hamiltonian with finite energy density, leads to absence of ergodicity and long-time coherent dynamics in closed quantum systems starting from simple initial states. Such non-ergodic coherent dynamics, where the system does not explore its entire phase space, has been experimentally observed in a chain of ultracold Rydberg atoms. We show, via study of a periodically driven Rydberg chain, that the drive frequency acts as a tuning parameter for several reentrant transitions between ergodic and non-ergodic regimes. The former regime shows rapid thermalization of correlation functions and absence of scars in the spectrum of the system's Floquet Hamiltonian. The latter regime, in contrast, has scars in its Floquet spectrum which control the long-time coherent dynamics of correlation functions. Our results open a new possibility of drive frequency-induced tuning between ergodic and non-ergodic dynamics in experimentally realizable disorder-free quantum many-body systems.

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.

preprint2016arXiv

Entanglement generation in periodically driven integrable systems: dynamical phase transitions and steady state

We study a class of periodically driven $d-$dimensional integrable models and show that after $n$ drive cycles with frequency $ω$, pure states with non-area-law entanglement entropy $S_n(l) \sim l^{α(n,ω)}$ are generated, where $l$ is the linear dimension of the subsystem, and $d-1 \le α(n,ω) \le d$. We identify and analyze the crossover phenomenon from an area ($S \sim l^{ d-1}$ for $d\geq1$) to a volume ($S \sim l^{d}$) law and provide a criterion for their occurrence which constitutes a generalization of Hastings' theorem to driven integrable systems in one dimension. We also find that $S_n$ generically decays to $S_{\infty}$ as $(ω/n)^{(d+2)/2}$ for fast and $(ω/n)^{d/2}$ for slow periodic drives; these two dynamical phases are separated by a topological transition in the eigensprectrum of the Floquet Hamiltonian. This dynamical transition manifests itself in the temporal behavior of all local correlation functions and does not require a critical point crossing during the drive. We find that these dynamical phases show a rich re-entrant behavior as a function of $ω$ for $d=1$ models, and also discuss the dynamical transition for $d>1$ models. Finally, we study entanglement properties of the steady state and show that singular features (cusps and kinks in $d=1$) appear in $S_{\infty}$ as a function of $ω$ whenever there is a crossing of the Floquet bands. We discuss experiments which can test our theory.

preprint2016arXiv

Level spectroscopy in a two-dimensional quantum magnet: Linearly dispersing spinons at the deconfined quantum critical point

We study the level structure of excitations at the "deconfined" critical point separating antiferromagnetic and valence-bond-solid phases in two-dimensional quantum spin systems using the $J$-$Q$ model as an example. Energy gaps in different spin ($S$) and momentum (${\bf k}$) sectors are extracted from imaginary-time correlation functions obtained in quantum Monte Carlo simulations. We find strong quantitative evidence for deconfined linearly dispersing spinons with gapless points at ${\bf k}=(0,0)$, $(π,0)$, $(0,π)$, and $(π,π)$, as inferred from two-spinon excitations ($S=0$ and $S=1$ states) around these points. We also observe a duality between singlet and triplet excitations at the critical point and inside the ordered phases, in support of an enhanced symmetry, possibly SO(5).

preprint2016arXiv

Parisi formula, disorder chaos and fluctuation for the ground state energy in the spherical mixed p-spin models

We show that the limiting ground state energy of the spherical mixed p-spin model can be identified as the infimum of certain variational problem. This complements the well-known Parisi formula for the limiting free energy in the spherical model. As an application, we obtain explicit formulas for the limiting ground state energy in the replica symmetry, one level of replica symmetry breaking and full replica symmetry breaking phases at zero temperature. In addition, our approach leads to new results on disorder chaos in spherical mixed even p-spin models. In particular, we prove that when there is no external field, the location of the ground state energy is chaotic under small perturbations of the disorder. We also establish that in the spherical mixed even p-spin model, the ground state energy superconcentrates in the absence of external field, while it obeys a central limit theorem if the external field is present.

preprint2016arXiv

Signatures and conditions for phase band crossings in periodically driven integrable systems

We present generic conditions for phase band crossings for a class of periodically driven integrable systems represented by free fermionic models subjected to arbitrary periodic drive protocols characterized by a frequency $ω_D$. These models provide a representation for the Ising and $XY$ models in $d=1$, the Kitaev model in $d=2$, several kinds of superconductors, and Dirac fermions in graphene and atop topological insulator surfaces. Our results demonstrate that the presence of a critical point/region in the system Hamiltonian (which is traversed at a finite rate during the dynamics) may change the conditions for phase band crossings that occur at the critical modes. We also show that for $d>1$, phase band crossings leave their imprint on the equal-time off-diagonal fermionic correlation functions of these models; the Fourier transforms of such correlation functions, $F_{\vec k_0}( ω_0)$, have maxima and minima at specific frequencies which can be directly related to $ω_D$ and the time at which the phase bands cross at $\vec k = \vec k_0$. We discuss the significance of our results in the contexts of generic Hamiltonians with $N>2$ phase bands and the underlying symmetry of the driven Hamiltonian.

preprint2015arXiv

Classical spin-liquid on the maximally frustrated honeycomb lattice

We show that the honeycomb Heisenberg antiferromagnet with $J_1/2=J_2=J_3$, where $J_{1/2/3}$ are first-, second- and third-neighbour couplings respectively, forms a classical spin liquid with pinch-point singularities in the structure factor at the Brillouin zone corners. Upon dilution with non-magnetic ions, fractionalised degrees of freedom carrying $1/3$ of the free moment emerge. Their effective description in the limit of low-temperature is that of spins randomly located on a triangular lattice, with a frustrated interaction of long-ranged logarithmic form. The XY version of this magnet exhibits nematic thermal order by disorder, which comes with a clear experimental diagnostic.

preprint2015arXiv

Coalescing Brownian flows: A new approach

The coalescing Brownian flow on $\mathbb{R}$ is a process which was introduced by Arratia [Coalescing Brownian motions on the line (1979) Univ. Wisconsin, Madison] and Tóth and Werner [Probab. Theory Related Fields 111 (1998) 375-452], and which formally corresponds to starting coalescing Brownian motions from every space-time point. We provide a new state space and topology for this process and obtain an invariance principle for coalescing random walks. This result holds under a finite variance assumption and is thus optimal. In previous works by Fontes et al. [Ann. Probab. 32 (2004) 2857-2883], Newman et al. [Electron. J. Probab. 10 (2005) 21-60], the topology and state-space required a moment of order $3-\varepsilon$ for this convergence to hold. The proof relies crucially on recent work of Schramm and Smirnov on scaling limits of critical percolation in the plane. Our approach is sufficiently simple that we can handle substantially more complicated coalescing flows with little extra work - in particular similar results are obtained in the case of coalescing Brownian motions on the Sierpinski gasket. This is the first such result where the limiting paths do not enjoy the noncrossing property.

preprint2015arXiv

Extinction window of mean field branching annihilating random walk

We study a model of growing population that competes for resources. At each time step, all existing particles reproduce and the offspring randomly move to neighboring sites. Then at any site with more than one offspring, the particles are annihilated. This is a nonmonotone model, which makes the analysis more difficult. We consider the extinction window of this model in the finite mean-field case, where there are $n$ sites but movement is allowed to any site (the complete graph). We show that although the system survives for exponential time, the extinction window is logarithmic.

preprint2015arXiv

Juntas in the $\ell^{1}$-grid and Lipschitz maps between discrete tori

We show that if $A \subset [k]^n$, then $A$ is $ε$-close to a junta depending upon at most $\exp(O(|\partial A|/(k^{n-1}ε)))$ coordinates, where $\partial A$ denotes the edge-boundary of $A$ in the $\ell^1$-grid. This is sharp up to the value of the absolute constant in the exponent. This result can be seen as a generalisation of the Junta theorem for the discrete cube, from [E. Friedgut, Boolean functions with low average sensitivity depend on few coordinates, Combinatorica 18 (1998), 27-35], or as a characterization of large subsets of the $\ell^1$-grid whose edge-boundary is small. We use it to prove a result on the structure of Lipschitz functions between two discrete tori; this can be seen as a discrete, quantitative analogue of a recent result of Austin [T. Austin, On the failure of concentration for the $\ell^{\infty}$-ball, preprint]. We also prove a refined version of our junta theorem, which is sharp in a wider range of cases.

preprint2015arXiv

Random Coulomb antiferromagnets: from diluted spin liquids to Euclidean random matrices

We study a disordered classical Heisenberg magnet with uniformly antiferromagnetic interactions which are frustrated on account of their long-range Coulomb form, {\em i.e.} $J(r)\sim -A\ln r$ in $d=2$ and $J(r)\sim A/r$ in $d=3$. This arises naturally as the $T\rightarrow 0$ limit of the emergent interactions between vacancy-induced degrees of freedom in a class of diluted Coulomb spin liquids (including the classical Heisenberg antiferromagnets on checkerboard, SCGO and pyrochlore lattices) and presents a novel variant of a disordered long-range spin Hamiltonian. Using detailed analytical and numerical studies we establish that this model exhibits a very broad paramagnetic regime that extends to very large values of $A$ in both $d=2$ and $d=3$. In $d=2$, using the lattice-Green function based finite-size regularization of the Coulomb potential (which corresponds naturally to the underlying low-temperature limit of the emergent interactions between orphan-spins), we only find evidence that freezing into a glassy state occurs in the limit of strong coupling, $A=\infty$, while no such transition seems to exist at all in $d=3$. We also demonstrate the presence and importance of screening for such a magnet. We analyse the spectrum of the Euclidean random matrices describing a Gaussian version of this problem, and identify a corresponding quantum mechanical scattering problem.

preprint2015arXiv

Velocity of excitations in ordered, disordered and critical antiferromagnets

We test three different approaches, based on quantum Monte Carlo simulations, for computing the velocity $c$ of triplet excitations in antiferromagnets. We consider the standard $S=1/2$ one- and two-dimensional Heisenberg models, as well as a bilayer Heisenberg model at its critical point. Computing correlation functions in imaginary time and using their long-time behavior, we extract the lowest excitation energy versus momentum using improved fitting procedures and a generalized moment method. The velocity is then obtained from the dispersion relation. We also exploit winding numbers to define a cubic space-time geometry, where the velocity is obtained as the ratio of the spatial and temporal lengths of the system when all winding number fluctuations are equal. The two methods give consistent results for both ordered and critical systems, but the winding-number estimator is more precise. For the Heisenberg chain, we accurately reproduce the exactly known velocity. For the two-dimensional Heisenberg model, our results are consistent with other recent calculations, but with an improved statistical precision; $c=1.65847(4)$. We also use the hydrodynamic relation $c^2=ρ_s/χ_\perp(q\to 0)$ between $c$, the spin stiffness $ρ_s$, and the transversal susceptibility $χ_\perp$, using the smallest non-zero momentum $q=2π/L$. This method also is well controlled in two dimensions, but the cubic criterion for winding numbers delivers better numerical precision. In one dimension the hydrodynamic relation is affected by logarithmic corrections which make accurate extrapolations difficult. As an application of the winding-number method, for the quantum-critical bilayer model our high-precision determination of the velocity enables us to quantitatively test, at an unprecedented level, field-theoretic predictions for low-temperature scaling forms where $c$ enters.

preprint2014arXiv

A topological spin glass in diluted spin ice

It is a salient experimental fact that a large fraction of candidate spin liquid materials freeze as the temperature is lowered. The question naturally arises whether such freezing is intrinsic to the spin liquid ("disorder-free glassiness") or extrinsic, in the sense that a topological phase simply coexists with standard freezing of impurities. Here, we demonstrate a surprising third alternative, namely that freezing and topological liquidity are inseparably linked. The topological phase reacts to the introduction of disorder by generating degrees of freedom of a new type (along with interactions between them), which in turn undergo a freezing transition while the topological phase supporting them remains intact.

preprint2014arXiv

On Testing Independence and Goodness-of-fit in Linear Models

We consider a linear regression model and propose an omnibus test to simultaneously check the assumption of independence between the error and the predictor variables, and the goodness-of-fit of the parametric model. Our approach is based on testing for independence between the residual obtained from the parametric fit and the predictor using the Hilbert--Schmidt independence criterion (Gretton et al. (2008)). The proposed method requires no user-defined regularization, is simple to compute, based merely on pairwise distances between points in the sample, and is consistent against all alternatives. We develop distribution theory for the proposed test statistic, both under the null and the alternative hypotheses, and devise a bootstrap scheme to approximate its null distribution. We prove the consistency of the bootstrap scheme. A simulation study shows that our method has better power than its main competitors. Two real datasets are analyzed to demonstrate the scope and usefulness of our method.

preprint2012arXiv

Ashkin-Teller criticality and pseudo first-order behavior in a frustrated Ising model on the square lattice

We study the challenging thermal phase transition to stripe order in the frustrated square-lattice Ising model with couplings J1<0 (nearest-neighbor, ferromagnetic) and J2>0 (second-neighbor, antiferromagnetic) for g=J2/|J1|>1/2. Using Monte Carlo simulations and known analytical results, we demonstrate Ashkin-Teller criticality for g>= g*, i.e., the critical exponents vary continuously between those of the 4-state Potts model at g=g* and the Ising model for g -> infinity. Thus, stripe transitions offer a route to realizing a related class of conformal field theories with conformal charge c=1 and varying exponents. The transition is first-order for g<g*= 0.67(1), much lower than previously believed, and exhibits pseudo first-order behavior for g* < g < 1.

preprint2012arXiv

Coalescing systems of non-Brownian particles

A well-known result of Arratia shows that one can make rigorous the notion of starting an independent Brownian motion at every point of an arbitrary closed subset of the real line and then building a set-valued process by requiring particles to coalesce when they collide. Arratia noted that the value of this process will be almost surely a locally finite set at all positive times, and a finite set almost surely if the initial value is compact: the key to both of these facts is the observation that, because of the topology of the real line and the continuity of Brownian sample paths, at the time when two particles collide one or the other of them must have already collided with each particle that was initially between them. We investigate whether such instantaneous coalescence still occurs for coalescing systems of particles where either the state space of the individual particles is not locally homeomorphic to an interval or the sample paths of the individual particles are discontinuous. We give a quite general criterion for a coalescing system of particles on a compact state space to coalesce to a finite set at all positive times almost surely and show that there is almost sure instantaneous coalescence to a locally finite set for systems of Brownian motions on the Sierpinski gasket and stable processes on the real line with stable index greater than one.

preprint2012arXiv

Geometric influences

We present a new definition of influences in product spaces of continuous distributions. Our definition is geometric, and for monotone sets it is identical with the measure of the boundary with respect to uniform enlargement. We prove analogs of the Kahn-Kalai-Linial (KKL) and Talagrand's influence sum bounds for the new definition. We further prove an analog of a result of Friedgut showing that sets with small "influence sum" are essentially determined by a small number of coordinates. In particular, we establish the following tight analog of the KKL bound: for any set in $\mathbb{R}^n$ of Gaussian measure $t$, there exists a coordinate $i$ such that the $i$th geometric influence of the set is at least $ct(1-t)\sqrt{\log n}/n$, where $c$ is a universal constant. This result is then used to obtain an isoperimetric inequality for the Gaussian measure on $\mathbb{R}^n$ and the class of sets invariant under transitive permutation group of the coordinates.

preprint2012arXiv

Geometric Influences II: Correlation Inequalities and Noise Sensitivity

In a recent paper, we presented a new definition of influences in product spaces of continuous distributions, and showed that analogues of the most fundamental results on discrete influences, such as the KKL theorem, hold for the new definition in Gaussian space. In this paper we prove Gaussian analogues of two of the central applications of influences: Talagrand's lower bound on the correlation of increasing subsets of the discrete cube, and the Benjamini-Kalai-Schramm (BKS) noise sensitivity theorem. We then use the Gaussian results to obtain analogues of Talagrand's bound for all discrete probability spaces and to reestablish analogues of the BKS theorem for biased two-point product spaces.

preprint2012arXiv

On reverse hypercontractivity

We study the notion of reverse hypercontractivity. We show that reverse hypercontractive inequalities are implied by standard hypercontractive inequalities as well as by the modified log-Sobolev inequality. Our proof is based on a new comparison lemma for Dirichlet forms and an extension of the Strook-Varapolos inequality. A consequence of our analysis is that {\em all} simple operators $L=Id-\E$ as well as their tensors satisfy uniform reverse hypercontractive inequalities. That is, for all $q<p<1$ and every positive valued function $f$ for $t \geq \log \frac{1-q}{1-p}$ we have $\| e^{-tL}f\|_{q} \geq \| f\|_{p}$. This should be contrasted with the case of hypercontractive inequalities for simple operators where $t$ is known to depend not only on $p$ and $q$ but also on the underlying space. The new reverse hypercontractive inequalities established here imply new mixing and isoperimetric results for short random walks in product spaces, for certain card-shufflings, for Glauber dynamics in high-temperatures spin systems as well as for queueing processes. The inequalities further imply a quantitative Arrow impossibility theorem for general product distributions and inverse polynomial bounds in the number of players for the non-interactive correlation distillation problem with $m$-sided dice.

preprint2012arXiv

Phase transitions in the frustrated Ising model on the square lattice

We consider the thermal phase transition from a paramagnetic to stripe-antiferromagnetic phase in the frustrated two-dimensional square-lattice Ising model with competing interactions J1<0 (nearest neighbor, ferromagnetic) and J2 >0 (second neighbor, antiferromagnetic). The striped phase breaks a Z4 symmetry and is stabilized at low temperatures for g=J2/|J1|>1/2. Despite the simplicity of the model, it has proved difficult to precisely determine the order and the universality class of the phase transitions. This was done convincingly only recently by Jin et al. [PRL 108, 045702 (2012)]. Here, we further elucidate the nature of these transitions and their anomalies by employing a combination of cluster mean-field theory, Monte Carlo simulations, and transfer-matrix calculations. The J1-J2 model has a line of very weak first-order phase transitions in the whole region 1/2<g<g*, where g* = 0.67(1). Thereafter, the transitions from g above g* are continuous and can be fully mapped, using universality arguments, to the critical line of the well known Ashkin-Teller model from its 4-state Potts point to the decoupled Ising limit. We also comment on the pseudo-first-order behavior at the Potts point and its neighborhood in the Ashkin-Teller model on finite lattices, which in turn leads to the appearance of similar effects in the vicinity of the multicritical point g* in the J1-J2 model. The continuous transitions near g* can therefore be mistaken to be first-order transitions, and this realization was the key to understanding the paramagnetic-striped transition for the full range of g>1/2. Most of our results are based on Monte Carlo calculations, while the cluster mean-field and transfer-matrix results provide useful methodological bench-marks for weakly first-order behaviors and Ashkin-Teller criticality.

preprint2012arXiv

Stochastic population growth in spatially heterogeneous environments

Classical ecological theory predicts that environmental stochasticity increases extinction risk by reducing the average per-capita growth rate of populations. To understand the interactive effects of environmental stochasticity, spatial heterogeneity, and dispersal on population growth, we study the following model for population abundances in $n$ patches: the conditional law of $X_{t+dt}$ given $X_t=x$ is such that when $dt$ is small the conditional mean of $X_{t+dt}^i-X_t^i$ is approximately $[x^iμ_i+\sum_j(x^j D_{ji}-x^i D_{ij})]dt$, where $X_t^i$ and $μ_i$ are the abundance and per capita growth rate in the $i$-th patch respectivly, and $D_{ij}$ is the dispersal rate from the $i$-th to the $j$-th patch, and the conditional covariance of $X_{t+dt}^i-X_t^i$ and $X_{t+dt}^j-X_t^j$ is approximately $x^i x^j σ_{ij}dt$. We show for such a spatially extended population that if $S_t=(X_t^1+...+X_t^n)$ is the total population abundance, then $Y_t=X_t/S_t$, the vector of patch proportions, converges in law to a random vector $Y_\infty$ as $t\to\infty$, and the stochastic growth rate $\lim_{t\to\infty}t^{-1}\log S_t$ equals the space-time average per-capita growth rate $\sum_iμ_i\E[Y_\infty^i]$ experienced by the population minus half of the space-time average temporal variation $\E[\sum_{i,j}σ_{ij}Y_\infty^i Y_\infty^j]$ experienced by the population. We derive analytic results for the law of $Y_\infty$, find which choice of the dispersal mechanism $D$ produces an optimal stochastic growth rate for a freely dispersing population, and investigate the effect on the stochastic growth rate of constraints on dispersal rates. Our results provide fundamental insights into "ideal free" movement in the face of uncertainty, the persistence of coupled sink populations, the evolution of dispersal rates, and the single large or several small (SLOSS) debate in conservation biology.

preprint2012arXiv

Vacancy-induced spin textures and their interactions in a classical spin liquid

Motivated by experiments on the archetypal frustrated magnet SrCr$_{9p}$Ga$_{12-9p}$O$_{19}$ (SCGO), we study the classical Heisenberg model on the pyrochlore slab (Kagomé bilayer) lattice with site-dilution $x=1-p$. This allows us to address generic aspects of the physics of non-magnetic vacancies in a classical spin liquid. We explicitly demonstrate that the pure ($x=0$) system remains a spin-liquid down to the lowest temperatures, with an unusual {\em non-monotonic} temperature dependence of the susceptibility, which even turns diamagnetic for the apical spins between the two kagome layers. For $x> 0$ but small, the low temperature magnetic response of the system is most naturally described in terms of the properties of spatially extended spin textures that cloak an "orphan" $S=3/2$ Cr$^{3+}$ spin in direct proximity to a pair of missing sites belonging to the same triangular simplex. In the $T \rightarrow 0$ limit, these orphan-texture complexes each carry a net magnetization that is exactly half the magnetic moment of an individual spin of the undiluted system. Furthermore, we demonstrate that they interact via an entropic {\em temperature dependent} {\em pair-wise exchange interaction} $J_{eff}(T,\vec{r}) \sim T {\mathcal J} (\vec{r} \sqrt{T})$ that has a logarithmic form at short-distances and decays exponentially beyond a thermal correlation length $ξ(T) \sim 1/\sqrt{T}$. The sign of $J_{eff}$ depends on whether the two orphan spins belong to the same Kagome layer or not. We provide a detailed analytical account of these properties using an effective field theory approach specifically tailored for the problem at hand. These results are in quantitative agreement with large-scale Monte Carlo numerics.

preprint2011arXiv

Absolute continuity of the limiting eigenvalue distribution of the random Toeplitz matrix

We show that the limiting eigenvalue distribution of random symmetric Toeplitz matrices is absolutely continuous with density bounded by 8, partially answering a question of Bryc, Dembo and Jiang (2006). The main tool used in the proof is a spectral averaging technique from the theory of random Schrödinger operators. The similar question for Hankel matrices remains open.

preprint2011arXiv

Spectra of large random trees

We analyze the eigenvalues of the adjacency matrices of a wide variety of random trees. Using general, broadly applicable arguments based on the interlacing inequalities for the eigenvalues of a principal submatrix of a Hermitian matrix and a suitable notion of local weak convergence for an ensemble of random trees, we show that the empirical spectral distributions for each of a number of random tree models converge to a deterministic (model dependent) limit as the number of vertices goes to infinity. We conclude for ensembles such as the linear preferential attachment models, random recursive trees, and the uniform random trees that the limiting spectral distribution has a set of atoms that is dense in the real line. We obtain precise asymptotics on the mass assigned to zero by the empirical spectral measures via the connection with the cardinality of a maximal matching. Moreover, we show that the total weight of a weighted matching is asymptotically equivalent to a constant multiple of the number of vertices when the edge weights are independent, identically distributed, non-negative random variables with finite expected value. We greatly extend a celebrated result obtained by Schwenk for the uniform random trees by showing that, under mild conditions, with probability converging to one, the spectrum of a realization is shared by at least one other tree. For the the linear preferential attachment model with parameter $a > -1$, we show that the suitably rescaled $k$ largest eigenvalues converge jointly.

preprint2010arXiv

Example of a first-order Néel to Valence-Bond-Solid transition in two-dimensions

We consider the $S=1/2$ Heisenberg model with nearest-neighbor interaction $J$ and an additional multi-spin interaction $Q_3$ on the square lattice. The $Q_3$ term consists of three bond-singlet projectors and is chosen to favor the formation of a valence-bond solid (VBS) where the valence bonds (singlet pairs) form a staggered pattern. The model exhibits a quantum phase transition from the Néel state to the VBS as a function of $Q_3/J$. We study the model using quantum Monte Carlo (stochastic series expansion) simulations. The Néel-VBS transition in this case is strongly first-order in nature, in contrast to similar previously studied models with continuous transitions into columnar VBS states. The qualitatively different transitions illustrate the important role of an emerging U(1) symmetry in the latter case, which is not possible in the present model due to the staggered VBS pattern (which does not allow local fluctuations necessary to rotate the local VBS order parameter).

preprint2010arXiv

Fractional spin textures in the frustrated magnet SrCr$_{9p}$Ga$_{12-9p}$O$_{19}$

We consider the archetypal frustrated antiferromagnet SrCr$_{9p}$Ga$_{12-9p}$O$_{19}$ in its well-known spin-liquid state, and demonstrate that a Cr$^{3+}$ spin $S=3/2$ ion in direct proximity to a pair of vacancies (in disordered $p<1$ samples) is cloaked by a spatially extended spin texture that encodes the correlations of the parent spin-liquid. In this spin-liquid regime, the combined object has a magnetic response identical to a classical spin of length $S/2 = 3/4$, which dominates over the small intrinsic susceptibility of the pure system. This fractional-spin texture leaves an unmistakable imprint on the measured $^{71}$Ga nuclear magnetic resonance (NMR) lineshapes, which we compute using Monte-Carlo simulations and compare with experimental data.

preprint2008arXiv

Collinear ordering of easy-axis triangular lattice antiferromagnets

Antiferromagnetically coupled moments on the frustrated triangular lattice typically order in a coplanar state at low temperature. Here, we demonstrate that the presence of not-very-large easy axis single ion anisotropy leads to an interesting orientationally ordered collinear state in triangular lattice antiferromagnets with moments $S \ge 3/2$. This ordered state breaks the symmetry of $π/3$ rotations about a lattice site, while leaving intact the translational symmetry of the lattice.

preprint2008arXiv

Semiclassical spin liquid state of easy axis Kagome antiferromagnets

Motivated by recent experiments on Nd-langasite, we consider the effect of strong easy axis single-ion anisotropy $D$ on $S > 3/2$ spins interacting with antiferromagnetic exchange $J$ on the Kagome lattice. When $T \ll DS^2$, the collinear low energy states selected by the anisotropy map on to configurations of the classical Kagome lattice Ising antiferromagnet. However, the low temperature limit is quite different from the cooperative Ising paramagnet that obtains classically for $T \ll JS^2$. We find that sub-leading ${\mathcal O}(J^3S/D^2)$ multi-spin interactions arising from the transverse quantum dynamics result in a crossover from an intermediate temperature classical cooperative Ising paramagnet to a semiclassical spin liquid with distinct short-ranged correlations for $T \ll J^3S/D^2$.