# Statistical Mechanics (cond-mat.stat-mech)

• We study a model of bosons confined to three coupled wells. The model describes interactions between bosons, tunneling of bosons between adjacent wells, and the effect of an external field. We consider both integrable and non-integrable regimes within the model. These regimes are first analyzed by means of the energy level spacing distribution. We then conduct a study of the quantum dynamics of the system to probe the conditions under which switching behavior can occur. Through variation of an external field on one well, we demonstrate how the system can be controlled between various switched-on and switched-off configurations.
• Chiral active particles (or self-propelled circle swimmers) feature a rich collective behavior, comprising rotating macro-clusters and micro-flock patterns which consist of phase-synchronized rotating clusters with a characteristic self-limited size. These patterns emerge from the competition of alignment interactions and rotations suggesting that they might occur generically in many chiral active matter systems. However, although excluded volume interactions occur naturally among typical circle swimmers, it is not yet clear if macro-clusters and micro-flock patterns survive their presence. The present work shows that both types of pattern do survive but feature strongly enhance fluctuations regarding the size and shape of the individual clusters. Despite these fluctuations, we find that the average micro-flock size still follows the same characteristic scaling law as in the absence of excluded volume interactions, i.e. micro-flock sizes scale linearly with the single-swimmer radius.
• We use information geometry, in which the local distance between models measures their distinguishability from data, to quantify the flow of information under the renormalization group. We show that information about relevant parameters is preserved, with distances along relevant directions maintained under flow. By contrast, irrelevant parameters become less distinguishable under the flow, with distances along irrelevant directions contracting according to renormalization group exponents. We develop a covariant formalism to understand the contraction of the model manifold. We then apply our tools to understand the emergence of the diffusion equation and more general statistical systems described by a free energy. Our results give an information-theoretic justification of universality in terms of the flow of the model manifold under coarse graining.
• We find a series of topological phase transitions of increasing order, beyond the more standard second-order phase transition in a one-dimensional topological superconductor. The jumps in the order of the transitions depend on the range of the pairing interaction, which is parametrized by an algebraic decay with exponent $\alpha$. Remarkably, in the limit $\alpha = 1$ the order of the topological transition becomes infinite. We compute the critical exponents for the series of higher-order transitions in exact form and find that they fulfil the hyperscaling relation. We also study the critical behaviour at the boundary of the system and discuss potential experimental platforms of magnetic atoms in superconductors.
• In this thesis, we considered quantum systems coupled to several baths. We supposed that the system state is governed by the quantum master equation (QME). We investigated the quantum pump and the excess entropy production. In the first half of the thesis, we investigated the quantum pump using the full counting statistics with quantum master equation (FCS-QME) approach. In the latter part of the thesis, we investigated the excess entropy production. The average entropy production is composed of the time integral of the instantaneous steady entropy production rate and the excess entropy production. We define average entropy production rate using the average energy and particle currents, which are calculated by using the full counting statistics with QME. The excess entropy production is given by a line integral in the control parameter space and its integrand is called the Berry-Sinitsyn-Nemenman (BSN) vector. In the weakly nonequilibrium regime, we show that BSN vector is described by $\ln \rho_0^{(-1)}$ and $\rho_0$ where $\rho_0$ is the instantaneous steady state of the QME and $\rho_0^{(-1)}$ is that of the QME which is given by reversing the sign of the Lamb shift term. In general, the potential dose not exist. The origins of the non-existence of the potential are a quantum effect (the Lamb shift) and the breaking of the time-reversal symmetry. The non-existence of the potential means that the excess entropy essentially depends on the path of the modulation. If the system Hamiltonian is non-degenerate or the Lamb shift term is negligible, the excess entropy production approximately reduces to the difference between the von Neumann entropies of the system. We pointed out that the expression of the entropy production obtained in the classical Markov jump process is different from our result and showed that these are approximately equivalent only in the weakly nonequilibrium regime.
• We perform large-scale molecular dynamics simulations to study heated granular fluids in three dimensions. Granular particles dissipate their kinetic energy due to solid frictional interaction with other particles. The velocity of each particle is perturbed by a uniformly-distributed random noise, which mimics the heating. At the early stage of evolution, the kinetic energy of the system decays with time and reaches a steady state at a later stage. The velocity distribution in the steady state shows a non-Gaussian distribution. This has been characterized by using the Sonine polynomial expansion for a wide range of densities. Particles show diffusive motion for densities below the \it jamming density $\phi_{\rm J}$
• We study the homogeneous nearest-neighbor Ising ferromagnet on the right half plane with a Dobrushin type boundary condition --- say plus on the top part of the boundary and minus on the bottom. For sufficiently low temperature $T$, we completely characterize the pure (i.e., extremal) Gibbs states, as follows. There is exactly one for each angle $\theta\in\lbrack-\pi/2,+\pi/2]$; here $\theta$ specifies the asymptotic angle of the interface separating regions where the spin configuration looks like that of the plus (respectively, minus) full-plane state. Some of these conclusions are extended all the way to $T=T_{c}$ by developing new Ising exact solution result -- in particular, there is at least one pure state for each $\theta$.
• We consider a collection of self-driven apolar particles on a substrate that organize into an active nematic phase at sufficiently high density or low noise. Using the dynamical renormalization group, we systematically study the 2d fluctuating ordered phase in a coarse-grained hydrodynamic description involving both the nematic director and the conserved density field. In the presence of noise, we show that the system always displays only quasi-long ranged orientational order beyond a crossover scale. A careful analysis of the nonlinearities permitted by symmetry reveals that activity is dangerously irrelevant over the linearized description, allowing giant number fluctuations to persist though now with strong finite-size effects and a non-universal scaling exponent. Nonlinear effects from the active currents lead to power law correlations in the density field thereby preventing macroscopic phase separation in the thermodynamic limit.
• We first present an empirical study of the Belief Propagation (BP) algorithm, when run on the random field Ising model defined on random regular graphs in the zero temperature limit. We introduce the notion of maximal solutions for the BP equations and we use them to fix a fraction of spins in their ground state configuration. At the phase transition point the fraction of unconstrained spins percolates and their number diverges with the system size. This in turn makes the associated optimization problem highly non trivial in the critical region. Using the bounds on the BP messages provided by the maximal solutions we design a new and very easy to implement BP scheme which is able to output a large number of stable fixed points. On one side this new algorithm is able to provide the minimum energy configuration with high probability in a competitive time. On the other side we found that the number of fixed points of the BP algorithm grows with the system size in the critical region. This unexpected feature poses new relevant questions on the physics of this class of models.
• We study the statistics, in stationary conditions, of the work $W_\tau$ done by the active force in different systems of self-propelled particles in a time $\tau$. We show the existence of a critical value $W_\tau ^\dag$ such that fluctuations with $W_\tau >W_\tau ^\dag$ correspond to configurations where interaction between particles plays a minor role whereas those with $W_\tau < W_\tau ^\dag$ represent states with single particles dragged by clusters. This two-fold behavior is fully mirrored by the probability distribution $P(W_\tau)$ of the work, which does not obey the large-deviation principle for $W_\tau <W_\tau ^\dag$. This pattern of behavior can be interpreted as due to a phase transition occurring at the level of fluctuating quantities and an order parameter is correspondingly identified.
• We provide a reply to a comment by I. Goychuk arXiv:1708.04155, version v2 from 27 Aug 2017 (not under active consideration with Phys. Rev. Lett.), on our Letter E. Aghion, D. A. Kessler and E. Barkai, Phys. Rev. Lett., 118, 260601 (2017).
• By studying the quench dynamics in one-dimensional superlattice systems with inversion symmetry, we find robust crossings in single-particle entanglement spectra for quantum quenches between different symmetry-protected topological phases. The physics behind this phenomenon is the emergence of a dynamical Chern number accompanied by unitarily created momentum-time Skyrmions. We also discuss a possible experimental situation based on Bloch-state tomography in ultracold atomic systems. Our work identifies the unique role of topology in quantum dynamics far from equilibrium.
• Matter exhibits phases and their transitions. These transitions are classified as first-order phase transitions (FOPTs) and continuous ones. While the latter has a well-established theory of the renormalization group, the former is only qualitatively accounted for by classical theories of nucleation, since their predictions often disagree with experiments by orders of magnitude. A theory to integrate FOPTs into the framework of the renormalization-group theory has been proposed but seems to contradict with extant wisdom and lacks numerical evidence. Here we show that universal hysteresis scaling as predicted by the renormalization-group theory emerges unambiguously when the theory is combined intimately with the theory of nucleation and growth in the FOPTs of the paradigmatic two-dimensional Ising model driven by a linearly varying externally applied field below its critical point. This not only provides a new method to rectify the nucleation theories, but also unifies the theories for both classes of transitions and FOPTs can be studied using universality and scaling similar to their continuous counterpart.
• The aim of the paper is to propose an algorithm based on the Multispin Coding technique for the Kob-Andersen glassy dynamics. We first give motivations to speed up the numerical simulation in the context of spin glass models [M. Mezard, G. Parisi, M. Virasoro, Spin Glass Theory and Beyond (World Scientific, Singapore, 1987)], after defining the Markovian dynamics as in [W. Kob, H.C. Andersen, Phys. Rev. E 48, 4364 (1993)] as well as the related interesting observables, we extend it to the more general framework of random regular graphs, listing at the same time some known analytical results [C. Toninelli, G. Biroli, D.S. Fisher, J. Stat. Phys. 120, 167 (2005)]. The purpose of this work is a dual one, firstly, we describe how bitwise operators can be used to build up the algorithm by carefully exploiting the way data are stored on a computer. Since it was first introduced [M. Creutz, L. Jacobs, C. Rebbi, Phys. Rev. D 20, 1915 (1979), C. Rebbi, R.H. Swendsen, Phys. Rev. D 21, 4094 (1980)], this technique has been widely used to perform Monte Carlo simulations for Ising and Potts spin systems, however, it can be successfully adapted to more complex systems in which microscopic parameters may assume boolean values. Secondly, we introduce a random graph in which a characteristic para- meter allows to tune the possible transition point. A consistent part is devoted to listing the numerical results obtained by running numerical simulations.
• In two dimensions (2D), dephasing by a bath cuts off Anderson localization that would otherwise occur at any energy density for fermions with disorder. For an isolated system with short-range interactions, the system can be its own bath, exhibiting diffusive (non-Markovian) thermal density fluctuations. We recast the dephasing of weak localization due to a diffusive bath as a self-interacting polymer loop. We investigate the critical behavior of the loop in $d=4-\epsilon$ dimensions, and find a nontrivial fixed point corresponding to temperature $T^* \sim \epsilon >0$ where the dephasing time diverges. Assuming this fixed point survives to $\epsilon=2$, we associate it with a toy version of the ergodic to many-body-localized (MBL) phase transition in 2D.

Thomas Klimpel Apr 20 2017 09:16 UTC

This paper [appeared][1] in February 2016 in the peer reviewed interdisciplinary journal Chaos by the American Institute of Physics (AIP).

It has been reviewed publicly by amateurs both [favorably][2] and [unfavorably][3]. The favorable review took the last sentence of the abstract ("These invalid

...(continued)
Māris Ozols Sep 15 2016 21:30 UTC

Here is a link for those who also haven't heard of SciPost before: https://scipost.org/

Zoltán Zimborás Sep 15 2016 18:12 UTC

This is the very first paper of SciPost, waiting for the first paper of "Quantum" (http://quantum-journal.org). There are radical (and good!) changes going on in scientific publishing.

Stephen Jordan Apr 15 2016 15:02 UTC

This is a beautiful set of lecture notes.

Mile Gu Nov 20 2015 05:04 UTC

Good question! There shouldn't be any contradiction with the correspondence principle. The reason here is that the quantum models are built to simulate the output behaviour of macroscopic, classical systems, and are not necessarily macroscopic themselves. When we compare quantum and classical comple

...(continued)
hong Nov 20 2015 00:40 UTC

Interesting results. But, just wondering, does it contradict to the correspondence principle?

Marco Tomamichel Jul 17 2015 05:12 UTC

I am no expert at all on strongly correlated systems or topological order, but since you refer to information theory in your abstract, let me still ask you: What is the justification for using $I_2(A:B) = H_2(A) + H_2(B) - H_2(AB)$ for the Rényi mutual information? This quantity has no information-t

...(continued)
Jonathan Oppenehim May 06 2015 14:29 UTC

This article has generated a fair bit of discussion. But I found a few of the statements puzzling (Edgar Lozano also). Take for example, Theorem 1 (ii) (reversibility) which appears to contradict a number of previous results. Should we understand your work function as "work in the paradigm where we

...(continued)
Chris Perry Apr 08 2015 09:02 UTC

Hi Felix,

Thanks a lot for your comments. Yes, $X$ can be any state and we've made this clearer now in the paper - we don't need to assume that it is diagonal.

We've also noted more prominently that achievability of our result in the full thermodynamics case is only when the target state is b

...(continued)
Felix Leditzky Apr 03 2015 11:06 UTC

A question to Theorem 1:

In the process description in eq. (2), you mention that $X$ is any arbitrary state. However, in the proof of Theorem 1, in the converse part you assume that $[X,\sigma]=0$, or equivalently that $X$ is diagonal in the eigenbasis of $\sigma$. Similarly, in the achievability

...(continued)