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

• We compute annihilation rates of metastable magnetic skyrmions using a form of Langer's theory in the intermediate-to-high damping (IHD) regime. We look at three possible paths to annihilation: isotropic collapse of an isolated skyrmion, isotropic collapse induced by another skyrmion and annihilation at a boundary. We find that the skyrmion's internal modes play a dominant role in the thermally activated transitions compared to the spin-wave excitations and that the relative contribution of internal modes depends on the nature of the transition process. Additionally, the eigenmodes at saddle point configurations are characterized by broken symmetries. Our calculations for a small skyrmion stabilized at zero-field show that the annihilation is largely dominated by the mechanism at the boundary, even though in this case the activation energy is higher than that of isotropic collapses. The potential source of stability of metastable skyrmions is therefore found not to lie in high activation energies, nor in the dynamics at the transition state, but comes from entropic narrowing in the saddle point region which leads to low attempt frequencies. This narrowing effect is found to be primarily associated with the skyrmion's internal modes. Isotropic collapse induced by another skyrmion exhibits the same internal energy barrier as a single skyrmion, but with a different entropic barrier. The probability of induced isotropic collapse is expected to increase with the number of skyrmions present on a racetrack.
• We develop a highly efficient method to numerically simulate thermal fluctuations and correlations in non-relativistic continuous bosonic one-dimensional systems. The method is suitable for arbitrary local interactions as long as the system remains dynamically stable. We start by proving the equivalence of describing the systems through the transfer matrix formalism and a Fokker-Planck equation for a distribution evolving in space. The Fokker-Planck equation is known to be equivalent to a stochastic differential (Itō) equation. The latter is very suitable for computer simulations, allowing the calculation of any desired correlation function. As an illustration, we apply our method to the case of two tunnel-coupled quasi-condensates of bosonic atoms. The results are compared to the predictions of the sine-Gordon model for which we develop analytic expression directly from the transfer matrix formalism.
• The synchronous dynamics of an array of excitable oscillators, coupled via a generic graph, is studied. Non homogeneous perturbations can grow and destroy synchrony, via a self-consistent instability which is solely instigated by the intrinsic network dynamics. By acting on the characteristic time-scale of the network modulation, one can make the examined system to behave as its (partially) averaged analog. This result if formally obtained by proving an extended version of the averaging theorem, which allows for partial averages to be carried out. Oscillation death are reported to follow the onset of the network driven instability.
• We present an algorithm to correct molecular dynamics non-equilibrium simulations based on the knowledge of time-dependent experimental data. This is inspired by a similar scheme used for equilibrium simulations, and is based on the principle of maximum caliber that guarantees that as little subjective information as possible is added to a model, besides the available experimental data. We also show that the same algorithm, in the case of simplified models, can be used to speed up the computational time needed by a simulation by one to two orders of magnitude.
• Speed limit for classical stochastic Markov processes with discrete states is studied. We find that a trade-off inequality exists between the speed of the state transformation and the entropy production. The dynamical activity determines the time scale and plays a crucial role in the inequality. For systems with stationary current, a similar trade-off inequality with the Hatano-Sasa entropy production gives a much better bound on the speed of the state transformation. Our inequalities contain only physically well-defined quantities, and thus the physical picture of these inequalities is clear.
• We provide a thermodynamically consistent description of energy, charge and spin transfers in a thermoelectric quantum-dot spin valve in the collinear configuration based on nonequilibrium Green's function and full counting statistics. We use the fluctuation theorem symmetry and the concept of entropy production to characterize the efficiency with which thermal gradients can be used to pump charges or spins against their chemical potentials, arbitrary far from equilibrium. Close to equilibrium, we recover the Onsager reciprocal relations and the connection to linear response notions of performance such as the figure of merit. We also identify regimes where work extraction is more efficient far then close from equilibrium.
• We study the minimally nonlinear irreversible heat engines in which the time-reversal symmetry for the systems may b e broken. The expressions for the power and the efficiency are derived, in which the effects of the nonlinear terms due to dissipations are included. We show that, as within the linear responses, the minimally nonlinear irreversible heat engines enable attainment of Carnot efficiency at positive power. We also find that the Curzon-Ahlborn limit imposed on the efficiency at maximum power can be overcomed if the time-reversal symmetry is broken.
• Tracking tracer particles in heterogeneous environments plays an important role in unraveling the material properties. These heterogeneous structures are often static and depend on the sample realizations. Sample-to-sample fluctuations of such disorder realizations sometimes become considerably large. When we investigate the sample-to-sample fluctuations, fundamental averaging procedures are a thermal average for a single disorder realization and the disorder average for different disorder realizations. Here, we report on the self-averaging and non-self-averaging transitions in quenched trap models with finite system sizes, where we consider the periodic and the reflecting boundary conditions. Sample-to-sample fluctuations of diffusivity greatly exceeds trajectory-to-trajectory fluctuations of diffusivity in the corresponding annealed model. For a single disorder realization, the time-averaged mean square displacement and position-dependent observables converge to constants with the aid of the existence of the equilibrium distribution. This is a manifestation of ergodicity. As a result, the time-averaged quantities do not depend on the initial condition nor on the thermal histories but depend crucially on the disorder realization.
• Experimental realizations of a variety of atomic binary Bose-Fermi mixtures have brought opportunities for studying composite quantum systems with different spin-statistics. The binary atomic mixtures can exhibit a structural transition from a mixture into phase separation as the boson-fermion interaction increases. By using a path-integral formalism to evaluate the grand partition function and thermodynamic grand potential, we obtain the effective potential of binary Bose-Fermi mixtures. Thermodynamic quantities in a broad range of temperatures and interactions are also derived. The structural transition can be identified as a loop of the effective potential curve, and the volume fraction of phase separation can be determined by the lever rule. For $^6$Li-$^7$Li and $^6$Li-$^{41}$K mixtures, we present the phase diagrams of the mixtures in a box potential at zero and finite temperatures. Due to the flexible densities of atomic gases, the construction of phase separation is more complicated when compared to conventional liquid or solid mixtures where the individual densities are fixed. For harmonically trapped mixtures, we use the local density approximation to map out the finite-temperature density profiles and present typical trap structures, including the mixture, partially separated phases, and fully separated phases.
• Using an operator approach based on a commutator scheme that has been previously applied to Richardson's reduced BCS model and the inhomogeneous Dicke model, we obtain general exact solvability requirements for an anisotropic central spin model with $XXZ$-type hyperfine coupling between the central spin and the spin bath, without any prior knowledge of integrability of the model. We outline the basic steps of the usage of the operator approach, and pedagogically summarize them into two \emphLemmas and two \emphConstraints. Through a step-by-step construction of the eigen-problem, we show that the condition $g'^2_j-g_j^2=c$ naturally arises for the model to be exactly solvable, where $c$ is a constant independent of the bath-spin index $j$, and $\{g_j\}$ and $\{g'_j\}$ are the longitudinal and transverse hyperfine interactions, respectively. The obtained conditions and the resulting Bethe ansatz equations are consistent with that in previous literature.
• We investigate the swim pressure exerted by non-chiral and chiral active particles on convex or concave circular boundaries. Active particles are modeled as non-interacting and non-aligning self-propelled Brownian particles. The convex and concave circular boundaries are used as models representing a fixed inclusion immersed in an active bath and a cavity (or container) enclosing the active particles, respectively. We first present a detailed analysis of the role of convex versus concave boundary curvature and of the chirality of active particles on their spatial distribution, chirality-induced currents, and the swim pressure they exert on the bounding surfaces. The results will then be used to predict the mechanical equilibria of suspended fluid enclosures (generically referred to as 'droplets') in a bulk with active particles being present either inside the bulk fluid or within the suspended droplets. We show that, while droplets containing active particles and suspended in a normal bulk behave in accordance with standard capillary paradigms, those containing a normal fluid exhibit anomalous behaviors when suspended in an active bulk. In the latter case, the excess swim pressure results in non-monotonic dependence of the inside droplet pressure on the droplet radius. As a result, we find a regime of anomalous capillarity for a single droplet, where the inside droplet pressure increases upon increasing the droplet size. In the case of two interconnected droplets, we show that mechanical equilibrium can occur also when they have different sizes. We further identify a regime of anomalous ripening, where two unequal-sized droplets can reach a final state of equal sizes upon interconnection, in stark contrast with the standard Ostwald ripening phenomenon, implying shrinkage of the smaller droplet in favor of the larger one.
• We compute the joint statistics of the momenta $p_i$ of $N$ non-interacting fermions in a trap, near the Fermi edge, with a particular focus on the largest one $p_{\max}$. For a $1d$ harmonic trap, momenta and positions play a symmetric role and hence, the joint statistics of momenta is identical to that of the positions. In particular, $p_{\max}$, as $x_{\max}$, is distributed according to the Tracy-Widom distribution. Here we show that novel "momentum edge statistics" emerge when the curvature of the potential vanishes, i.e. for "flat traps" near their minimum, with $V(x) \sim x^{2n}$ and $n>1$. These are based on generalisations of the Airy kernel that we obtain explicitly. The fluctuations of $p_{\max}$ are governed by new universal distributions determined from the $n$-th member of the second Painlevé hierarchy of non-linear differential equations, with connections to multicritical random matrix models. Finite temperature extensions and possible experimental signatures in cold atoms are discussed.

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)