Publications 2014
Publications de l'année 2014 :

A new dipolar potential for numerical simulations of polar fluids on the 4D hypersphere
JeanMichel Caillol ^{1} Martin Trulsson ^{2}
The Journal of Chemical Physics, 2014, 141, pp.124111
We present a new method for Monte Carlo or Molecular Dynamics numerical simulations of three dimensional polar fluids. The simulation cell is defined to be the surface of the northern hemisphere of a fourdimensional (hyper)sphere. The point dipoles are constrained to remain tangent to the sphere and their interactions are derived from the basic laws of electrostatics in this geometry. The dipoledipole potential has two singularities which correspond to the following boundary conditions : when a dipole leaves the northern hemisphere at some point of the equator, it reappears at the antipodal point bearing the same dipole moment. We derive all the formal expressions needed to obtain the thermodynamic and structural properties of a polar liquid at thermal equilibrium in actual numerical simulation. We notably establish the expression of the static dielectric constant of the fluid as well as the behavior of the pair correlation at large distances. We report and discuss the results of extensive numerical Monte Carlo simulations for two reference states of a fluid of dipolar hard spheres and compare these results with previous methods with a special emphasis on finite size effects.
 1. LPT  Laboratoire de Physique Théorique d'Orsay [Orsay]
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

A singularpotential random matrix model arising in meanfield glassy systems
Gernot Akemann ^{1} Dario Villamaina ^{2} Pierpaolo Vivo ^{3}
Physical Review E, American Physical Society, 2014, 89, pp.062146
We consider an invariant random matrix model where the standard Gaussian potential is distorted by an additional single pole of order $m$. We compute the average or macroscopic spectral density in the limit of large matrix size, solving the loop equation with the additional constraint of vanishing trace on average. The density is generally supported on two disconnected intervals lying on the two sides of the pole. In the limit of having no pole, we recover the standard semicircle. Obtained in the planar limit, our results apply to matrices with orthogonal, unitary or symplectic symmetry alike. The orthogonal case with $m=2$ is motivated by an application to spin glass physics. In the SherringtonKirkpatrick meanfield model, in the paramagnetic phase and for sufficiently large systems the spin glass susceptibility is a random variable, depending on the realization of disorder. It is essentially given by a linear statistics on the eigenvalues of the coupling matrix. As such its large deviation function can be computed using standard Coulomb fluid techniques. The resulting free energy of the associated fluid precisely corresponds to the partition function of our random matrix model. Numerical simulations provide an excellent confirmation of our analytical results.
 1. Fakultät für Physik
 2. LPTENS  Laboratoire de Physique Théorique de l'ENS
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

A universal set of qubit quantum channels
Daniel Braun ^{1} Olivier Giraud ^{2} Ion Nechita ^{1} Clement Pellegrini ^{3} Marko Znidaric ^{4}
Journal of Physics A: Mathematical and Theoretical, Institute of Physics: Hybrid Open Access, 2014, 47 (13), pp.135302. <10.1088/17518113/47/13/135302>
We investigate the set of quantum channels acting on a single qubit. We provide a compact generalization of the FujiwaraAlgoet conditions for complete positivity to nonunital qubit channels, which we then use to characterize the possible geometric forms of the pure output of the channel. We provide universal sets of quantum channels for all unital qubit channels as well as for all extremal (not necessarily unital) qubit channels, in the sense that all qubit channels in these sets can be obtained by concatenation of channels in the corresponding universal set. We also show that our universal sets are essentially minimal.
 1. Information et Chaos Quantiques (LPT)
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. IMT  Institut de Mathématiques de Toulouse
 4. Department of Physics, Faculty of Mathematics and Physics

Absence of a gap in the Gaffnian state
Thierry Jolicoeur ^{1} Takahiro Mizusaki ^{2} Philippe Lecheminant ^{3}
Physical Review B (Condensed Matter), American Physical Society, 2014, 90, pp.075116
We study the Gaffnian trial wavefunction proposed to describe fractional quantum Hall correlations at Bose filling factor $\nu=2/3$ and Fermi filling $\nu=2/5$. A family of Hamiltonians interpolating between a hardcore interaction for which the physics is known and a projector whose ground state is the Gaffnian is studied in detail. We give evidence for the absence of a gap by using largescale exact diagonalizations in the spherical geometry. This is in agreement with recent arguments based on the fact that this wavefunction is constructed from a nonunitary conformal field theory. By using the cylinder geometry, we discuss in detail the nature of the underlying minimal model and we show the appearance of heterotic conformal towers in the edge energy spectrum where left and right movers are generated by distinct primary operators.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Institute of Natural Sciences
 3. LPTM  Laboratoire de Physique Théorique et Modélisation

Archive ouverte HAL – Repulsive Casimir Effect with Chern Insulators
Pablo RodriguezLopez ^{1} Adolfo Grushin ^{2}
Physical Review Letters, American Physical Society, 2014, 112 (5), pp.056804. 〈10.1103/PhysRevLett.112.056804〉
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. ICMM  Instituto de Ciencia de Materiales de Madrid

Area distribution of twodimensional random walks and non Hermitian Hofstadter quantum mechanics
Sergey Matveenko ^{1, 2} Stephane Ouvry ^{2}
Journal of Physics A: Mathematical and Theoretical, Institute of Physics: Hybrid Open Access, 2014, 47, pp.185001
When random walks on a square lattice are biased horizontally to move solely to the right, the probability distribution of their algebraic area can be exactly obtained. We explicitly map this biased classical random system on a non hermitian Hofstadterlike quantum model where a charged particle on a square lattice coupled to a perpendicular magnetic field hopps only to the right. In the commensurate case when the magnetic flux per unit cell is rational, an exact solution of the quantum model is obtained. Periodicity on the lattice allows to relate traces of the Nth power of the Hamiltonian to probability distribution generating functions of biased walks of length N.
 1. Landau Institute for Theoretical Physics
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Bound states and entanglement in the excited states of quantum spin chains
Jan Mölter ^{1} Thomas Barthel ^{2, 1} Ulrich Schollwöck ^{1} Vincenzo Alba ^{1}
Journal of Statistical Mechanics, 2014, pp.P10029
We investigate entanglement properties of the excited states of the spin1/2 Heisenberg (XXX) chain with isotropic antiferromagnetic interactions, by exploiting the Bethe ansatz solution of the model. We consider eigenstates obtained from both real and complex solutions ("strings") of the Bethe equations. Physically, the former are states of interacting magnons, whereas the latter contain bound states of groups of particles. We first focus on the situation with few particles in the chain. Using exact results and semiclassical arguments, we derive an upper bound S_MAX for the entanglement entropy. This exhibits an intermediate behavior between logarithmic and extensive, and it is saturated for highlyentangled states. As a function of the eigenstate energy, the entanglement entropy is organized in bands. Their number depends on the number of blocks of contiguous BetheTakahashi quantum numbers. In presence of bound states a significant reduction in the entanglement entropy occurs, reflecting that a group of bound particles behaves effectively as a single particle. Interestingly, the associated entanglement spectrum shows edgerelated levels. At finite particle density, the semiclassical bound S_MAX becomes inaccurate. For highlyentangled states S_A\propto L_c, with L_c the chord length, signaling the crossover to extensive entanglement. Finally, we consider eigenstates containing a single pair of bound particles. No significant entanglement reduction occurs, in contrast with the fewparticle case.
 1. Department of Physics and Arnold Sommerfeld Center for Theoretical Physics
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Calculation of multifractal dimensions in spin chains
Y. Y. Atas ^{1} E. Bogomolny ^{1}
Philosophical transactions of the Royal Society of London. Series A: Mathematical and physical sciences, Royal Society, The, 2014, 372, pp.20120520
It was demonstrated in [Phys. Rev. E 86, 021104, (2012)], that the groundstate wave functions for a large variety of onedimensional spin1/2 models are multifractals in the natural spinz basis. We present here the details of analytical derivations and numerical confirmations of these results.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

CasimirPolder interaction for gently curved surfaces
Giuseppe Bimonte ^{1} Thorsten Emig ^{2} Mehran Kardar ^{3}
Physical Review D, American Physical Society, 2014, 90, pp.081702
We use a derivative expansion for gently curved surfaces to compute the leading and the nexttoleading curvature corrections to the CasimirPolder interaction between a polarizable small particle and a nonplanar surface. While our methods apply to any homogeneous and isotropic surface, explicit results are presented here for perfect conductors. We show that the derivative expansion of the CasimirPolder potential follows from a resummation of its perturbative series, for small inplane momenta. We consider the retarded, nonretarded and classical high temperature limits.
 1. INFN, Sezione di Napoli  Istituto Nazionale di Fisica Nucleare, Sezione di Napoli
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. Department of Physics

Cauchy’s formulas for random walks in bounded domains
Mazzolo Alain ^{1} De Mulatier Clélia ^{2, 3} Zoia Andrea ^{4}
Journal of Mathematical Physics, American Institute of Physics (AIP), 2014, 55, pp.083308
Cauchy's formula was originally established for random straight paths crossing a body $B \subset \mathbb{R}^{n}$ and basically relates the average chord length through $B$ to the ratio between the volume and the surface of the body itself. The original statement was later extended in the context of transport theory so as to cover the stochastic paths of Pearson random walks with exponentially distributed flight lengths traversing a bounded domain. Some heuristic arguments suggest that Cauchy's formula may also hold true for Pearson random walks with arbitrarily distributed flight lengths. For such a broad class of stochastic processes, we rigorously derive a generalized Cauchy's formula for the average length travelled by the walkers in the body, and show that this quantity depends indeed only on the ratio between the volume and the surface, provided that some constraints are imposed on the entrance step of the walker in $B$. Similar results are obtained also for the average number of collisions performed by the walker in $B$, and an extension to absorbing media is discussed.
 1.
 2.
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 4.

Classical Heisenberg spins with longrange interactions: relaxation to equilibrium for finite systems
Julien Barré ^{1} Shamik Gupta ^{2}
Journal of Statistical Mechanics: Theory and Experiment, IOP Science, 2014, 2014 (P02017), pp.20. <10.1088/17425468/2014/02/P02017>
Systems with longrange interactions often relax towards statistical equilibrium over timescales that diverge with N , the number of particles. A recent work [S. Gupta and D. Mukamel, J. Stat. Mech.: Theory Exp. P03015 (2011)] analyzed a model system comprising N globally coupled classical Heisenberg spins and evolving under classical spin dynamics. It was numerically shown to relax to equilibrium over a time that scales superlinearly with N . Here, we present a detailed study of the LenardBalescu operator that accounts at leading order for the finiteN effects driving this relaxation. We demonstrate that corrections at this order are identically zero, so that relaxation occurs over a time longer than of order N , in agreement with the reported numerical results.
 1. JAD  Laboratoire Jean Alexandre Dieudonné
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Clustering of branching Brownian motions in confined geometries
Zoia AndreaDumonteil EricMazzolo AlainDe Mulatier CléliaRosso Alberto
Physical Review E, American Physical Society, 2014, 90, pp.042118
We study the evolution of a collection of individuals subject to Brownian diffusion, reproduction and disappearance. In particular, we focus on the case where the individuals are initially prepared at equilibrium within a confined geometry. Such systems are widespread in physics and biology and apply for instance to the study of neutron populations in nuclear reactors and the dynamics of bacterial colonies, only to name a few. The fluctuations affecting the number of individuals in space and time may lead to a strong patchiness, with particles clustered together. We show that the analysis of this peculiar behaviour can be rather easily carried out by resorting to a backward formalism based on the Green's function, which allows the key physical observables, namely, the particle concentration and the pair correlation function, to be explicitly derived.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Condensation transition in joint large deviations of linear statistics
Juraj SzavitsNossan ^{1} Martin R. Evans ^{1} Satya N. Majumdar ^{2}
Journal of Physics A: Mathematical and Theoretical, Institute of Physics: Hybrid Open Access, 2014, 47, pp.455004
Real space condensation is known to occur in stochastic models of mass transport in the regime in which the globally conserved mass density is greater than a critical value. It has been shown within models with factorised stationary states that the condensation can be understood in terms of sums of independent and identically distributed random variables: these exhibit condensation when they are conditioned to a large deviation of their sum. It is well understood that the condensation, whereby one of the random variables contributes a finite fraction to the sum, occurs only if the underlying probability distribution (modulo exponential) is heavytailed, i.e. decaying slower than exponential. Here we study a similar phenomenon in which condensation is exhibited for nonheavytailed distributions, provided random variables are additionally conditioned on a large deviation of certain linear statistics. We provide a detailed theoretical analysis explaining the phenomenon, which is supported by Monte Carlo simulations (for the case where the additional constraint is the sample variance) and demonstrated in several physical systems. Our results suggest that the condensation is a generic phenomenon that pertains to both typical and rare events.
 1. SUPA, School of Physics, University of Edinburgh
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Constraint driven condensation in large fluctuations of linear statistics
Juraj SzavitsNossan ^{1} Martin R. Evans ^{1} Satya N. Majumdar ^{2}
Physical Review Letters, American Physical Society, 2014, 112, pp.020602
Condensation is the phenomenon whereby one of a sum of random variables contributes a finite fraction to the sum. It is manifested as an aggregation phenomenon in diverse physical systems such as coalescence in granular media, jamming in traffic and gelation in networks. We show here that the same condensation scenario, which normally happens only if the underlying probability distribution has tails heavier than exponential, can occur for lighttailed distributions in the presence of additional constraints. We demonstrate this phenomenon on the sample variance, whose probability distribution conditioned on the particular value of the sample mean undergoes a phase transition. The transition is manifested by a change in behavior of the large deviation rate function.
 1. SUPA, School of Physics, University of Edinburgh
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Counterions between or at asymmetrically charged walls: 2D freefermion point
L. Samaj ^{1, 2} E. Trizac ^{1}
Journal of Statistical Physics, Springer Verlag (Germany), 2014, 156, pp.932
This work contributes to the problem of determining effective interaction between asymmetrically (likely or oppositely) charged objects whose total charge is neutralized by mobile pointlike counterions of the same charge, the whole system being in thermal equilibrium. The problem is formulated in two spatial dimensions with logarithmic Coulomb interactions. The charged objects correspond to two parallel lines at distance d, with fixed line charge densities. Two versions of the model are considered: the standard "unconstrained" one with particles moving freely between the lines and the "constrained" one with particles confined to the lines. We solve exactly both systems at the freefermion coupling and compare the results for the pressure (i.e. the force between the lines per unit length of one of the lines) with the meanfield PoissonBoltzmann solution. For the unconstrained model, the larged asymptotic behaviour of the freefermion pressure differs from that predicted by the meanfield theory. For the constrained model, the asymptotic pressure coincides with the attractive van der WaalsCasimir fluctuational force. For both models, there are fundamental differences between the cases of likelycharged and oppositelycharged lines, the latter case corresponding at large distances d to a capacitor.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Institute of Physics

Critical dynamics and domain motion from permittivity of the electronic ferroelectric (TMTTF)2AsF6
S. Brazovskii ^{1, 2} Pierre Monceau ^{3} F. Ya. Nad ^{4}
International schoolconference on electronic crystals: ECRYS2014, Aug 2014, Cargèse, France. 460, pp.79, Physica B: Condensed Matter
The quasi onedimensional organic conductor (TMTTF)2AsF6 shows the charge ordering transition at Tc101K to a state of the ferroelectric Mott insulator which is still well conducting. We present and interpret the experimental data on the gigantic dielectric response in the vicinity of TCO, concentrating on the frequency dependence of the inverse $1/\epsilon$ of the complex permittivity $\epsilon=\epsilon^\prime+i\epsilon^{\prime\prime}$. Surprisingly for a ferroelectric, we could closely approach the 2nd order phase transition and to deeply reach the critical dynamics of the polarization. We could analyse the critical slowingdown when approaching Tc from both sides and to extract the anomalous power law for the frequency dependence of the order parameter viscosity. Moreover, below Tc we could extract a sharp absorption feature coming from a motion of domain walls which shows up at a frequency well below the relaxation rate.
 1. International Institute of Physics
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. MagSup
 4. Kotelnikov Institute of Radio Engineering and Electronics

Delocalization of weakly interacting bosons in a 1D quasiperiodic potential
V. P. Michal ^{1} B. L. Altshuler ^{2} G. V. Shlyapnikov ^{1, 3, 4}
Physical Review Letters, American Physical Society, 2014, 113, pp.045304
We consider weakly interacting bosons in a 1D quasiperiodic potential (AubryAzbelHarper model) in the regime where all singleparticle states are localized. We show that the interparticle interaction may lead to the manybody delocalization and we obtain the finitetemperature phase diagram. Counterintuitively, in a wide range of parameters the delocalization requires stronger cou pling as the temperature increases. This means that the system of bosons can undergo a transition from a fluid to insulator (glass) state under heating.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Physics Dept.
 3. VAN DER WAALSZEEMAN INSTITUTE  University of Amsterdam Van der WaalsZeeman Institute
 4. Russian Quantum Center

Dielectric Square Resonator Investigated with Microwave Experiments
S. Bittner ^{1} E. Bogomolny ^{2} B. Dietz ^{3} M. MiskiOglu ^{3} A. Richter ^{3}
Physical Review E, American Physical Society, 2014, 90, pp.052909
We present a detailed experimental study of the symmetry properties and the momentum space representation of the field distributions of a dielectric square resonator as well as the comparison with a semiclassical model. The experiments have been performed with a flat ceramic microwave resonator. Both the resonance spectra and the field distributions were measured. The momentum space representations of the latter evidenced that the resonant states are each related to a specific classical torus, leading to the regular structure of the spectrum. Furthermore, they allow for a precise determination of the refractive index. Measurements with different arrangements of the emitting and the receiving antennas were performed and their influence on the symmetry properties of the field distributions was investigated in detail, showing that resonances with specific symmetries can be selected purposefully. In addition, the length spectrum deduced from the measured resonance spectra and the trace formula for the dielectric square resonator are discussed in the framework of the semiclassical model.
 1. LPQM  Laboratoire de Photonique Quantique et Moléculaire
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. Institut für Kernphysik

Diffusion in periodic, correlated random forcing landscapes
David S Dean ^{1} Shamik Gupta ^{2, *} Gleb Oshanin ^{3} Alberto Rosso ^{2} Grégory Schehr ^{2}
Journal of Physics A: Mathematical and Theoretical, IOP Publishing, 2014, 47 (37), pp.372001. <10.1088/17518113/47/37/372001>
We study the dynamics of a Brownian particle in a strongly correlated quenched random potential defined as a periodicallyextended (with period $L$) finite trajectory of a fractional Brownian motion with arbitrary Hurst exponent $H \in (0,1)$. While the periodicity ensures that the ultimate longtime behavior is diffusive, the generalised Sinai potential considered here leads to a strong logarithmic confinement of particle trajectories at intermediate times. These two competing trends lead to dynamical frustration and result in a rich statistical behavior of the diffusion coefficient $D_L$: Although one has the typical value $D^{\rm typ}_L \sim \exp(\beta L^H)$, we show via an exact analytical approach that the positive moments ($k>0$) scale like $\langle D^k_L \rangle \sim \exp{[c' (k \beta L^{H})^{1/(1+H)}]}$, and the negative ones as $\langle D^{k}_L \rangle \sim \exp(a' (k \beta L^{H})^2)$, $c'$ and $a'$ being numerical constants and $\beta$ the inverse temperature. These results demonstrate that $D_L$ is strongly nonselfaveraging. We further show that the probability distribution of $D_L$ has a lognormal left tail and a highly singular, onesided logstable right tail reminiscent of a Lifshitz singularity.
 1. LOMA  Laboratoire Ondes et Matière d'Aquitaine
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. LPTMC  Laboratoire de Physique Théorique de la Matière Condensée

Diffusion with resetting in arbitrary spatial dimension
Martin R. Evans ^{1} Satya N. Majumdar ^{2}
Journal of Physics A: Mathematical and Theoretical, Institute of Physics: Hybrid Open Access, 2014, 47, pp.285001
We consider diffusion in arbitrary spatial dimension d with the addition of a resetting process wherein the diffusive particle stochastically resets to a fixed position at a constant rate $r$. We compute the nonequilibrium stationary state which exhibits nonGaussian behaviour. We then consider the presence of an absorbing target centred at the origin and compute the survival probability and mean time to absorption of the diffusive particle by the target. The mean absorption time is finite and has a minimum value at an optimal resetting rate $r^*$ which depends on dimension. Finally we consider the problem of a finite density of diffusive particles, each resetting to its own initial position. While the typical survival probability of the target at the origin decays exponentially with time regardless of spatial dimension, the average survival probability decays asymptotically as $\exp A (\log t)^d$ where $A$ is a constant. We explain these findings using an interpretation as a renewal process and arguments invoking extreme value statistics.
 1. SUPA, School of Physics, University of Edinburgh
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Domainwall melting in ultracold boson systems with holes and spinflip defects
Jad C. Halimeh ^{1} Anton Wöllert ^{1} Ian P. McCulloch ^{2} Ulrich Schollwöck ^{1} Thomas Barthel ^{1, 3}
Physical Review A, American Physical Society, 2014, 89, pp.063603
Quantum magnetism is a fundamental phenomenon of nature. As of late, it has garnered a lot of interest because experiments with ultracold atomic gases in optical lattices can be used as a simulator for phenomena of magnetic systems. A paradigmatic example is the timeevolution of a domainwall state of a spin1/2 Heisenberg chain, the socalled domainwall melting. The model can be implemented by having two species of bosonic atoms with unity filling and strong onsite repulsion U in an optical lattice. In this paper, we study the domainwall melting in such a setup on the basis of the timedependent density matrix renormalization group (tDMRG). We are particularly interested in the effects of defects that originate from an imperfect preparation of the initial state. Typical defects are holes (empty sites) and flipped spins. We show that the dominating effects of holes on observables like the spatially resolved magnetization can be taken account of by a linear combination of spatially shifted observables from the clean case. For sufficiently large $U$, further effects due to holes become negligible. In contrast, the effects of spin flips are more severe as their dynamics occur on the same timescale as that of the domainwall melting itself. It is hence advisable to avoid preparation schemes that are based on spinflips.
 1. Department of Physics and Arnold Sommerfeld Center for Theoretical Physics
 2. Centre for Engineered Quantum Systems
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Dynamical patterns of phase transformations from selftrapping of quantum excitons
Tianyou Yi ^{1} Natasha Kirova ^{2} Serguei Brazovskii ^{1}
International schoolconference on electronic crystals: ECRYS2014, Aug 2014, Cargèse, France. 460, pp.73, Physica B: Condensed Matter
Phase transitions induced by short optical pulses is a new mainstream in studies of cooperative electronic states. Its special realization in systems with neutralionic transformations stands out in a way that the optical pumping goes to excitons rather than to electronic bands. We present a semiphenomenological modeling of spaciotemporal effects applicable to any system where the optical excitons are coupled to a symmetry breaking order parameter. In our scenario, after a short initial pulse of photons, a quasicondensate of excitons appears as a macroscopic quantum state which then evolves interacting with other degrees of freedom prone to instability. This coupling leads to selftrapping of excitons; that locally enhances their density which can surpass a critical value to trigger the phase transformation, even if the mean density is below the required threshold. The system is stratified in domains which evolve through dynamical phase transitions and may persist even after the initiating excitons have recombined. We recover dynamic interplays of fields such as the excitons wave function, electronic charge transfer and polarization, lattice deformations.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. LPS  Laboratoire de Physique des Solides

Dynamics and termination cost of spatially coupled meanfield models
Francesco Caltagirone ^{1} Silvio Franz ^{2} Richard Morris ^{1} Lenka Zdeborová ^{1}
Physical Review E, American Physical Society, 2014, 89, pp.012102
This work is motivated by recent progress in information theory and signal processing where the socalled 'spatially coupled' design of systems leads to considerably better performance. We address relevant open questions about spatially coupled systems through the study of a simple Ising model. In particular, we consider a chain of CurieWeiss models that are coupled by interactions up to a certain range. Indeed, it is well known that the pure (uncoupled) CurieWeiss model undergoes a first order phase transition driven by the magnetic field, and furthermore, in the spinodal region such systems are unable to reach equilibrium in subexponential time if initialized in the metastable state. By contrast, the spatially coupled system is, instead, able to reach the equilibrium even when initialized to the metastable state. The equilibrium phase propagates along the chain in the form of a travelling wave. Here we study the speed of the wavefront and the socalled 'termination cost' \textit{i.e.}, the conditions necessary for the propagation to occur. We reach several interesting conclusions about optimization of the speed and the cost.
 1. IPHT  Institut de Physique Théorique (ex SPhT)
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Elastic MultiBody Interactions on a Lattice
D. S. Petrov ^{1}
Physical Review A, American Physical Society, 2014, 90, pp.021601
We show that by coupling two hyperfine states of an atom in an optical lattice one can independently control two, three, and fourbody onsite interactions in a nonperturbative manner. In particular, under typical conditions of current experiments one can have a purely three or fourbody interacting gas of $^{39}$K atoms characterized by onsite interaction shifts of several 100Hz.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Ensemble inequivalence in a meanfield XY model with ferromagnetic and nematic couplings
Arkady Pikovsky ^{1, 2} Shamik Gupta ^{3} Tarcisio N. Teles ^{4} Fernanda P. C. Benetti ^{4} Renato Pakter ^{4} Yan Levin ^{4} Stefano Ruffo ^{5}
Physical Review E, American Physical Society, 2014, pp.062141
We explore ensemble inequivalence in longrange interacting systems by studying an XY model of classical spins with ferromagnetic and nematic coupling. We demonstrate the inequivalence by mapping the microcanonical phase diagram onto the canonical one, and also by doing the inverse mapping. We show that the equilibrium phase diagrams within the two ensembles strongly disagree within the regions of firstorder transitions, exhibiting interesting features like temperature jumps. In particular, we discuss the coexistence and forbidden regions of different macroscopic states in both the phase diagrams.
 1. UNIPOTSDAM  Department of Physics and Astronomie, Univ. Potsdam
 2. Department of Control Theory
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 4. IFUFRGS  Instituto de Física
 5. Dipartimento di Fisica e Astronomia and CSDC

Entanglement Generation Using Discrete Solitons in Coulomb Crystals
H. Landa ^{1, 2} A. Retzker ^{3} T. Schaetz ^{4} B. Reznik ^{1}
Physical Review Letters, American Physical Society, 2014, 113, pp.053001
Laser cooled and trapped ions can crystallize and feature discrete solitons, that are nonlinear, topologicallyprotected configurations of the Coulomb crystal. Such solitons, as their continuum counterparts, can move within the crystal, while their discreteness leads to the existence of a gapseparated, spatiallylocalized motional mode of oscillation above the spectrum. Suggesting that these unique properties of discrete solitons can be used for generating entanglement between different sites of the crystal, we study a detailed proposal in the context of stateoftheart experimental techniques. We analyze the interaction of periodicallydriven planar ion crystals with optical forces, revealing the effects of micromotion in radiofrequency traps inherent to such structures, as opposed to linear ion chains. The proposed method requires Doppler cooling of the crystal and sideband cooling of the soliton's localized modes alone. Since the gap separation of the latter is nearly independent of the crystal size, this approach could be particularly useful for producing entanglement and studying systemenvironment interactions in large, two and possibly threedimensional systems.
 1. Raymond and Beverly Sackler Faculty of Exact Sciences
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. HUJ  Hebrew University of Jerusalem
 4. University of Freiburg [Freiburg]

Exact nonequilibrium solutions of the Boltzmann equation under a timedependent external force
David GuéryOdelin ^{1} J. G. Muga ^{2, 3} M. J. RuizMontero ^{4} E. Trizac ^{5}
Physical Review Letters, American Physical Society, 2014, 112, pp.180602
We construct a novel class of exact solutions to the Boltzmann equation, in both its classical and quantum formulation, for arbitrary collision laws. When the system is subjected to a specific external forcing, the precise form of which is worked out, non equilibrium dampingless solutions are admissible. They do not contradict the $H$theorem, but are constructed from its requirements. Interestingly, these solutions hold for timedependent confinement. We exploit them, in a reverseengineering perspective, to work out a protocol that shortcuts any adiabatic transformation between two equilibrium states in an arbitrarily short timespan, for an interacting system. Particle simulations of the direct Monte Carlo type fully corroborate the analytical predictions.
 1. Atomes Froids (LCAR)
 2. Departamento de QumicaFsica, UPVEHU
 3. Department of Physics, Shanghai University
 4. Fisica Teorica, Universidad de Sevilla
 5. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Excitonic Mechanism of Local Phase Transformations by Optical Pumping
Serguei Brazovskii ^{1} Natasha Kirova ^{2}
Journal of Superconductivity and Novel Magnetism, Springer Verlag (Germany), 2014, 27, pp.1009  1013. <10.1007/s1094801324329>
Transformations of cooperative electronic states by impacts of optical pumping or/and electrostatic doping is a new mainstream in physics of correlated systems. Here we present a semiphenomenological modelling of spaciotemporal effects in a system where the light absorption goes through a channel creating the excitons – intramolecular ones or bound electronhole pairs and finally the condensate of optical excitons feeds and stimulates phase transformations. Interacting with a nearcritical order parameter and deformations, the excitons are subject to selftrapping. That locally enhances their density which can surpass a critical value to trigger the phase transformation, even if the mean density is below the required threshold. The model can be used e.g. as a simplified version of optically induced neutralionic transitions in organic chain compounds.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. LPS  Laboratoire de Physique des Solides

Experimental and theoretical investigation of angular dependence of the Casimir force between sinusoidally corrugated surfaces
A. A. Banishev ^{1} J. Wagner ^{1} T. Emig ^{2} R. Zandi ^{1} U. Mohideen ^{1}
Physical Review B (Condensed Matter), American Physical Society, 2014, 89, pp.235436
In the current work we present the complete results for the measurement of normal Casimir force between shallow and smooth sinusoidally corrugated gold coated sphere and a plate at various angles between the corrugations using an atomic force microscope. All measured data were compared with the theoretical approach using the proximity force approximation and theory based on derivative expansion. In both cases real material properties of the surfaces and nonzero temperature were taken into account. Special attention is paid to the description of electrostatic interactions between corrugated surfaces at different angels between corrugations and samples reparation and characterization. The measured forces are found to be in good agreement with the theory including correlation effects of geometry and material properties and deviate significantly from the predictions of the proximity force approximation approach. This provides the quantitative confirmation for the observation of diffractiontype effects that are disregarded within the PFA approach. The obtained results open new opportunities for control of the Casimir effect in micromechanical systems.
 1. Department of Physics and Astronomy
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

First order transition for the optimal search time of Lévy flights with resetting
Lukasz Kusmierz ^{1, 2} Satya N. Majumdar ^{3} Sanjib Sabhapandit ^{4} Gregory Schehr ^{3}
Physical Review Letters, American Physical Society, 2014, 113, pp.220602
We study analytically an intermittent search process in one dimension. There is an immobile target at the origin and a searcher undergoes a discrete time jump process starting at $x_0\geq0$, where successive jumps are drawn independently from an arbitrary jump distribution $f(\eta)$. In addition, with a probability $0\leq r \leq1$ the position of the searcher is reset to its initial position $x_0$. The efficiency of the search strategy is characterized by the mean time to find the target, i.e., the mean first passage time (MFPT) to the origin. For arbitrary jump distribution $f(\eta)$, initial position $x_0$ and resetting probability $r$, we compute analytically the MFPT. For the heavytailed L\'evy stable jump distribution characterized by the L\'evy index $0<\mu < 2$, we show that, for any given $x_0$, the MFPT has a global minimum in the $(\mu,r)$ plane at $(\mu^*(x_0),r^*(x_0))$. We find a remarkable firstorder phase transition as $x_0$ crosses a critical value $x_0^*$ at which the optimal parameters change discontinuously. Our analytical results are in good agreement with numerical simulations.
 1. Institute of Physics of Krakow
 2. Department of Automatics and Biomedical Engineering
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 4. Raman Research Institute

Fluctuating interfaces subject to stochastic resetting
Shamik Gupta ^{1} Satya N. Majumdar ^{1} Gregory Schehr ^{1}
Physical Review Letters, American Physical Society, 2014, 112, pp.220601
We study onedimensional fluctuating interfaces of length $L$ where the interface stochastically resets to a fixed initial profile at a constant rate $r$. For finite $r$ in the limit $L \to \infty$, the system settles into a nonequilibrium stationary state with nonGaussian interface fluctuations, which we characterize analytically for the KardarParisiZhang and EdwardsWilkinson universality class. Our results are corroborated by numerical simulations. We also discuss the generality of our results for a fluctuating interface in a generic universality class.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Fluctuations in partitioning systems with few degrees of freedom
L. Cerino ^{1} G. Gradenigo ^{2, 3} A. Sarracino ^{2, 4, 5} D. Villamaina ^{2} A. Vulpiani ^{6, 7}
Physical Review E, American Physical Society, 2014, 89, pp.042105
We study the behavior of a moving wall in contact with a particle gas and subjected to an external force. We compare the fluctuations of the system observed in the microcanonical and canonical ensembles, at varying the number of particles. Static and dynamic correlations signal significant differences between the two ensembles. Furthermore, velocityvelocity correlations of the moving wall present a complex twotime relaxation which cannot be reproduced by a standard Langevinlike description. Quite remarkably, increasing the number of gas particles in an elongated geometry, we find a typical timescale, related to the interaction between the partitioning wall and the particles, which grows macroscopically.
 1. Dipartimento di Fisica [Rome]
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. CEA  CEA Saclay
 4. Dipartimento di Fisica
 5. Kavli Institute for Theoretical Physics China
 6. INFN  Istituto Nazionale di Fisica Nucleare [Sezione di Roma 1]
 7. Kavli Institute for Theoretical Physics China,

Fluctuations of random matrix products and 1D Dirac equation with random mass
Kabir Ramola ^{1} Christophe Texier ^{1}
Journal of Statistical Physics, Springer Verlag (Germany), 2014, 157, pp.497514
We study the fluctuations of certain random matrix products $\Pi_N=M_N\cdots M_2M_1$ of $\mathrm{SL}(2,\mathbb{R})$, describing localisation properties of the onedimensional 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 $\Gamma_2=\lim_{N\to\infty}\mathrm{Var}(\ln\Pi_N)/N$. The case studied exhibits a saturation of the variance at low energy $\varepsilon$ along with a vanishing Lyapunov exponent $\Gamma_1=\lim_{N\to\infty}\ln\Pi_N/N$, leading to the behaviour $\Gamma_2/\Gamma_1\sim\ln(1/\varepsilon)\to\infty$ as $\varepsilon\to0$. Our continuum description sheds new light on the KappusWegner (band center) anomaly.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Fractal free energy landscapes in structural glasses
Patrick Charbonneau ^{1, 2} Jorge Kurchan ^{3} Giorgio Parisi ^{4, 5} Pierfrancesco Urbani ^{6, 7} Francesco Zamponi ^{2}
Nature Communications, Nature Publishing Group, 2014, 5, pp.3725. 〈10.1038/ncomms4725〉
Glasses are amorphous solids whose constituent particles are caged by their neighbors and thus cannot flow. This sluggishness is often ascribed to the free energy landscape containing multiple minima (basins) separated by high barriers. Here we show, using theory and numerical simulation, that the landscape is much rougher than is classically assumed. Deep in the glass, it undergoes a "roughness transition" to fractal basins. This brings about isostaticity at jamming and marginality of glassy states near jamming. Critical exponents for the basin width, the weak force distribution, and the spatial spread of quasicontacts at jamming can be analytically determined. Their value is found to be compatible with numerical observations. This advance therefore incorporates the jamming transition of granular materials into the framework of glass theory. Because temperature and pressure control which features of the landscape are experienced, glass mechanics and transport are expected to reflect the features of the topology we discuss here. Hitherto mysterious properties of lowtemperature glasses could be explained by this approach.
 1. Duke university [Durham]
 2. LPTENS  Laboratoire de Physique Théorique de l'ENS
 3. LPS  Laboratoire de Physique Statistique de l'ENS
 4. SMC/INFM  Department of Physics, Sapienza University of Rome
 5. INFN  Istituto Nazionale di Fisica Nucleare [Sezione di Roma 1]
 6. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 7. IPHT  Institut de Physique Théorique  UMR CNRS 3681

From generalized directed animals to the asymmetric simple exclusion process
Nils Haug ^{1} Sergei Nechaev ^{2, 3, 4} Mikhail Tamm ^{4, 5}
Journal of Statistical Mechanics, 2014, pp.P10013
Using the generalized normally ordered form of words in a locallyfree group of $n$ generators, we show that in the limit $n\to\infty$, the partition function of weighted directed lattice animals on a semiinfinite strip coincides with the partition function of stationary configurations of the asymmetric simple exclusion process (ASEP) with arbitrary entry/escape rates through open boundaries. We relate the features of the ASEP in the different regimes of the phase diagram to the geometric features of the associated generalized directed animals by showing the results of numerical simulations. In particular, we show how the presence of shocks at the first order transition line translates into the directed animal picture. Using the evolution equation for generalized, weighted Lukasiewicz paths, we also provide a straightforward calculation of the known ASEP generating function.
 1. School of Mathematical Sciences
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. P. N. Lebedev Physical Institute
 4. Physics Department
 5. Department of Applied Mathematics

From microstructural features to effective toughness in disordered brittle solids
Vincent Démery ^{1} Laurent Ponson ^{1} Alberto Rosso ^{2}
Europhysics Letters, EDP Science, 2014, 105, pp.34003
The relevant parameters at the microstructure scale that govern the macroscopic toughness of disordered brittle materials are investigated theoretically. We focus on planar crack propagation and describe the front evolution as the propagation of a longrange elastic line within a plane with random distribution of toughness. Our study reveals two regimes: in the collective pinning regime, the macroscopic toughness can be expressed as a function of a few parameters only, namely the average and the standard deviation of the local toughness distribution and the correlation lengths of the heterogeneous toughness field; in the individual pinning regime, the passage from micro to macroscale is more subtle and the full distribution of local toughness is required to be predictive. Beyond the failure of brittle solids, our findings illustrate the complex filtering process of microscale quantities towards the larger scales into play in a broad range of systems governed by the propagation of an elastic interface in a disordered medium.
 1. IJLRA  Institut Jean Le Rond d'Alembert
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Generalized Onsager theory for strongly anisometric patchy colloids
H. H. Wensink ^{1} E. Trizac ^{2}
Journal of Chemical Physics, American Institute of Physics (AIP), 2014, 140, pp.024901
The implications of soft 'patchy' interactions on the orientational disorderorder transition of strongly elongated colloidal rods and flat disks is studied within a simple Onsagervan der Waals density functional theory. The theory provides a generic framework for studying the liquid crystal phase behaviour of highly anisometric cylindrical colloids, which carry a distinct geometrical pattern of repulsive or attractive soft interactions localised on the particle surface. In this paper, we apply our theory to the case of charged rods and disks for which the local electrostatic interactions can be described by a screenedCoulomb potential. We consider infinitely thin rodlike cylinders with a uniform line charge and infinitely thin discotic cylinders with several distinctly different surface charge patterns. Irrespective of the backbone shape, the isotropicnematic phase diagrams of charged colloids feature a generic destabilization of nematic order, a dramatic narrowing of the biphasic density region and a reentrant phenomenon upon reducing the electrostatic screening. At higher particle density the electrostatic repulsion leads to a complete suppression of nematic order in favour of spatially inhomogeneous liquid crystals.
 1. LPS  Laboratoire de Physique des Solides
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Geometrical origins of contractility in disordered actomyosin networks
Lenz Martin ^{1}
Physical Review X, 2014, 4, pp.041002
Movement within eukaryotic cells largely originates from localized forces exerted by myosin motors on scaffolds of actin filaments. Although individual motors locally exert both contractile and extensile forces, large actomyosin structures at the cellular scale are overwhelmingly contractile, suggesting that the scaffold serves to favor contraction over extension. While this mechanism is well understood in highly organized striated muscle, its origin in disordered networks such as the cell cortex is unknown. Here we develop a mathematical model of the actin scaffold's local two or threedimensional mechanics and identify four competing contraction mechanisms. We predict that one mechanism dominates, whereby local deformations of the actin break the balance between contraction and extension. In this mechanism, contractile forces result mostly from motors plucking the filaments transversely rather than buckling them longitudinally. These findings sheds light on recent $\textit{in vitro}$ experiments, and provides a new geometrical understanding of contractility in the myriad of disordered actomyosin systems found $\textit{in vivo}$.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Heat transfer between anisotropic nanopartricles: Enhancement and switching
Roberta Incardone ^{1, 2} Thorsten Emig ^{3} Matthias Krüger ^{1, 2}
Europhysics Letters, EDP Science, 2014, 106, pp.41001
We theoretically study heat transfer between two anisotropic nanoparticles in vacuum, and derive closed expressions in terms of the anisotropic dipole polarizabilities. We show that transfer between two small spheroids can be many times as large as the one for two spheres of same volumes. Such increase with anisotropy is also found for the heat emission of an isolated small spheroid. Furthermore, we observe a strong dependence of transfer on the relative orientation, yielding the interpretation as a heat transfer switch. The switch quality, given as the ratio of transfer in the ''on'' and ''off'' positions, is observed to be as large as $10^3$ in the near field and even larger in the far field.
 1. MaxPlanckInstitut fur Intelligente Systeme
 2. Institut fur Theoretische Physik IV
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Index Distribution of Cauchy Random Matrices
Ricardo Marino ^{1} Satya N. Majumdar ^{1} Grégory Schehr ^{1} Pierpaolo Vivo ^{1}
Journal of Physics A: Mathematical and Theoretical, Institute of Physics: Hybrid Open Access, 2014, 47, pp.055001
Using a Coulomb gas technique, we compute analytically the probability $\mathcal{P}_\beta^{(C)}(N_+,N)$ that a large $N\times N$ Cauchy random matrix has $N_+$ positive eigenvalues, where $N_+$ is called the index of the ensemble. We show that this probability scales for large $N$ as $\mathcal{P}_\beta^{(C)}(N_+,N)\approx \exp\left[\beta N^2 \psi_C(N_+/N)\right]$, where $\beta$ is the Dyson index of the ensemble. The rate function $\psi_C(\kappa)$ is computed in terms of single integrals that are easily evaluated numerically and amenable to an asymptotic analysis. We find that the rate function, around its minimum at $\kappa=1/2$, has a quadratic behavior modulated by a logarithmic singularity. As a consequence, the variance of the index scales for large $N$ as $\mathrm{Var}(N_+)\sim \sigma_C\ln N$, where $\sigma_C=2/(\beta\pi^2)$ is twice as large as the corresponding prefactor in the Gaussian and Wishart cases. The analytical results are checked by numerical simulations and against an exact finite $N$ formula which, for $\beta=2$, can be derived using orthogonal polynomials.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Inferring the origin of an epidemic with a dynamic messagepassing algorithm
Andrey Y. Lokhov ^{1} Marc Mézard ^{1, 2} Hiroki Ohta ^{1} Lenka Zdeborová ^{3}
Physical Review E, American Physical Society, 2014, 90, pp.012801
We study the problem of estimating the origin of an epidemic outbreak  given a contact network and a snapshot of epidemic spread at a certain time, determine the infection source. Finding the source is important in different contexts of computer or social networks. We assume that the epidemic spread follows the most commonly used susceptibleinfectedrecovered model. We introduce an inference algorithm based on dynamic messagepassing equations, and we show that it leads to significant improvement of performance compared to existing approaches. Importantly, this algorithm remains efficient in the case where one knows the state of only a fraction of nodes.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. ENS  Ecole Normale Supérieure
 3. IPHT  Institut de Physique Théorique (ex SPhT)

Inhomogeneous and nonstationary Hall states of the CDW with quantized normal carriers
Serguei Brazovskii ^{1}
International schoolconference on electronic crystals: ECRYS2014, Aug 2014, Cargèse, France. 460, pp.236, Physica B: Condensed Matter
We suggest a theory for a deformable and sliding charge density wave (CDW) in the Hall bar geometry for the quantum limit when the carriers in remnant small pockets are concentrated at lowest Landau levels (LL) forming a fractionally ($\nu<1$) filled quantum Hall state. The gigantic polarizability of the CDW allows for a strong redistribution of electronic densities up to a complete charge segregation when all carriers occupy, with the maximum filling, a fraction $\nu$ of the chain length  thus forming the integer quantum Hall state, while leaving the fraction $(1\nu)$ of the chain length unoccupied. The electric field in charged regions easily exceeds the pinning threshold of the CDW, then the depinning propagates into the nominally pinned central region via sharp domain walls. Resulting picture is that of compensated collective and normal pulsing countercurrents driven by the Hall voltage. This scenario is illustrated by numerical modeling for nonstationary distributions of the current and the electric field. This picture can interpret experiments in mesajunctions showing depinning by the Hall voltage and the generation of voltagecontrolled high frequency oscillations (Yu.I. Latyshev, P. Monceau, A.A. Sinchenko, et al, presented at ECRYS2011, unpublished).
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Instanton Approach to Large $N$ HarishChandraItzyksonZuber Integrals
Joel Bun ^{1, 2, 3} JeanPhilippe Bouchaud ^{1} Satya N. Majumdar ^{2} Marc Potters ^{1}
Physical Review Letters, American Physical Society, 2014, 113, pp.070201
We reconsider the large $N$ asymptotics of HarishChandraItzyksonZuber integrals. We provide, using Dyson's Brownian motion and the method of instantons, an alternative, transparent derivation of the Matytsin formalism for the unitary case. Our method is easily generalized to the orthogonal and symplectic ensembles. We obtain an explicit solution of Matytsin's equations in the case of Wigner matrices, as well as a general expansion method in the dilute limit, when the spectrum of eigenvalues spreads over very wide regions.
 1. Science et Finance
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. Pole Univ Leonard de Vinci

Interplay of anisotropy in shape and interactions in charged platelet suspensions
Sara JabbariFarouji ^{1, 2} JeanJacques Weis ^{3} Patrick Davidson ^{4} Pierre Levitz ^{5} Emmanuel Trizac ^{2}
Journal of Chemical Physics, American Institute of Physics (AIP), 2014, 141, pp.224510
Motivated by the intriguing phase behavior of charged colloidal platelets, we investigate the structure and dynamics of charged repulsive disks by means of MonteCarlo simulations. The electrostatic interactions are taken into account through an effective twobody potential, obtained within the nonlinear PoissonBoltzmann formalism, which has the form of anisotropic screened Coulomb potential. Recently, we showed that the original intrinsic anisotropy of the electrostatic potential in competition with excluded volume effects leads to a rich phase behavior that not only includes various liquidcrsytalline phases but also predicts the existence of novel structures composed of alternating nematicantinematic sheets. Here, we examine the structural and dynamical signatures of each of the observed structures for both translational and rotational degrees of freedom. Finally, we discuss the influence of effective charge value and our results in relation to experimental findings on charged platelet suspensions.
 1. LIPhy, Laboratoire Interdisciplinaire de Physique
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. Laboratoire de Physique Théorique
 4. LPS  Laboratoire de Physique des Solides
 5. PECSA  Physicochimie des Electrolytes, Colloïdes et Sciences Analytiques

Islands of stability in motif distributions of random networks
M. V. Tamm ^{1, 2} A. B. Shkarin ^{3} V. A. Avetisov ^{2, 4} O. V. Valba ^{2, 5, 6} S. K. Nechaev ^{2, 6, 7}
Physical Review Letters, American Physical Society, 2014, 113, pp.095701
We consider random nondirected networks subject to dynamics conserving vertex degrees and study analytically and numerically equilibrium threevertex motif distributions in the presence of an external field, $h$, coupled to one of the motifs. For small $h$ the numerics is well described by the "chemical kinetics" for the concentrations of motifs based on the law of mass action. For larger $h$ a transition into some trapped motif state occurs in Erd\H{o}sRényi networks. We explain the existence of the transition by employing the notion of the entropy of the motif distribution and describe it in terms of a phenomenological Landautype theory with a nonzero cubic term. A localization transition should always occur if the entropy function is nonconvex. We conjecture that this phenomenon is the origin of the motifs' pattern formation in real evolutionary networks.
 1. Physics Department
 2. Department of Applied Mathematics
 3. Department of Physics
 4. The Semenov Institute of Chemical Physics
 5. MIPT  Moscow Institute of Physics and Technology
 6. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 7. P. N. Lebedev Physical Institute

Kovacslike memory effect in driven granular gases
A. Prados ^{1, 2} E. Trizac ^{1}
Physical Review Letters, American Physical Society, 2014, 112, pp.198001
While memory effects have been reported for dense enough disordered systems such as glasses, we show here by a combination of analytical and simulation techniques that they are also intrinsic to the dynamics of dilute granular gases. By means of a certain driving protocol, we prepare the gas in a state where the granular temperature $T$ coincides with its long time limit. However, $T$ does not subsequently remain constant, but exhibits a nonmonotonic evolution before reaching its nonequilibrium steady value. The corresponding socalled Kovacs hump displays a normal behavior for weak dissipation (as observed in molecular systems), but is reversed under strong dissipation, where it thus becomes anomalous.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Fisica Teorica

Kuramoto model of synchronization: Equilibrium and nonequilibrium aspects
Gupta Shamik ^{1} Campa Alessandro ^{2} Ruffo Stefano ^{3}
Journal of Statistical Mechanics: Theory and Experiment, IOP Science, 2014, pp.R08001
Recently, there has been considerable interest in the study of spontaneous synchronization, particularly within the framework of the Kuramoto model. The model comprises oscillators with distributed natural frequencies interacting through a meanfield coupling, and serves as a paradigm to study synchronization. In this review, we put forward a general framework in which we discuss in a unified way known results with more recent developments obtained for a generalized Kuramoto model that includes inertial effects and noise. We describe the model from a different perspective, highlighting the longrange nature of the interaction between the oscillators, and emphasizing the equilibrium and outofequilibrium aspects of its dynamics from a statistical physics point of view. In the first main section, we introduce the model and discuss for the noiseless and noisy dynamics and unimodal frequency distributions the synchronization transition that occurs in the stationary state. In the second section, we introduce the generalized model, and discuss its synchronization phase diagram for unimodal frequency distributions. In the third section, we describe deviations from the meanfield setting of the Kuramoto model by considering the generalized dynamics on a onedimensional periodic lattice on the sites of which the oscillators reside and interact with a coupling that decays as an inverse powerlaw of their separation. For two specific cases, namely, in the absence of noise and inertia, and in the case when the natural frequencies are the same for all the oscillators, we discuss how the longtime transition to synchrony is governed by the dynamics of the meanfield mode (zero Fourier mode) of the spatial distribution of the oscillator phases.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Department of Technology and Health
 3. Dipartimento di Fisica e Astronomia and CSDC

Lamplighter model of a random copolymer adsorption on a line
L. I. Nazarov ^{1} S. K. Nechaev ^{2, 3} M. V. Tamm ^{4, 5}
Condensed Matter Physics, Institute for Condensed Matter Physics, 2014, 17, pp.33002
We present a model of an ABdiblock random copolymer sequential selfpackaging with local quenched interactions on a onedimensional infinite sticky substrate. It is assumed that the AA and BB contacts are favorable, while AB are not. The position of a newly added monomer is selected in view of the local contact energy minimization. The model demonstrates a selforganization behavior with the nontrivial dependence of the total energy, $E$ (the number of unfavorable contacts), on the number of chain monomers, $N$: $E\sim N^{3/4}$ for quenched random equally probable distribution of A and Bmonomers along the chain. The model is treated by mapping it onto the "lamplighter" random walk and the diffusioncontrolled chemical reaction of $X+X\to 0$ type with the subdiffusive motion of reagents.
 1. Institute of Mining, Siberian Branch
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. P. N. Lebedev Physical Institute
 4. Physics Department
 5. Department of Applied Mathematics

Local logarithmic correlators as limits of Coulomb gas integrals
Raoul Santachiara ^{1} Jacopo Viti ^{2}
Nuclear Physics B, Elsevier, 2014, 882, pp.229262
We will describe how logarithmic singularities arise as limits of Coulomb Gas integrals. Our approach will combine analytic properties of the timelike Liouville structure constants, together with the recursive formula of the Virasoro conformal blocks. Although the Coulomb Gas formalism forces a diagonal coupling between the chiral and antichiral sectors of the Conformal Field Theory (CFT), we present new results for the multiscreening integrals which are potentially interesting for applications to critical statistical systems described by Logarithmic CFTs. In particular our findings extend and complement previous results, derived with Coulomb Gas methods, at $c=0$ and $c=2$.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. LPTENS  Laboratoire de Physique Théorique de l'ENS

Longrange correlations in a locally driven exclusion process
Tridib Sadhu ^{1} Satya N. Majumdar ^{2} David Mukamel ^{3}
Physical Review E, American Physical Society, 2014, 90, pp.012109
We show that the presence of a driven bond in an otherwise diffusive lattice gas with simple exclusion interaction results in longrange densitydensity correlation in its stationary state. In dimensions $d>1$ we show that in the thermodynamic limit this correlation decays as $C(r,s)\sim (r^2+s^2)^{d}$ at large distances $r$ and $s$ away from the drive with $rs>>1$. This is derived using an electrostatic analogy whereby $C(r,s)$ is expressed as the potential due to a configuration of electrostatic charges distributed in $2d$dimension. At bulk density $\rho=1/2$ we show that the potential is that of a localized quadrupolar charge. At other densities the same is correct in leading order in the strength of the drive and is argued numerically to be valid at higher orders.
 1. IPHT  Institut de Physique Théorique (ex SPhT)
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. Weizmann Institute

Maximal distance travelled by N vicious walkers till their survival
Anupam Kundu ^{1} Satya N. Majumdar ^{1} Gregory Schehr ^{1}
Journal of Statistical Physics, Springer Verlag (Germany), 2014, 157, pp.124157
We consider $N$ Brownian particles moving on a line starting from initial positions ${\bf{u}}\equiv \{u_1,u_2,\dots u_N\}$ such that $0
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Measurements of fiducial and differential cross sections for Higgs boson production in the diphoton decay channel at $\sqrt{s}=8$ TeV with ATLAS
S. Franz ^{1} S. Gupta ^{1} and all
Journal of High Energy Physics (JHEP), 2014, 09, pp.112. <10.1007/JHEP09(2014)112>
Measurements of fiducial and differential cross sections are presented for Higgs boson production in protonproton collisions at a centreofmass energy of $\sqrt{s}=8$ TeV. The analysis is performed in the $H \rightarrow \gamma\gamma$ decay channel using 20.3 fb$^{1}$ of data recorded by the ATLAS experiment at the CERN Large Hadron Collider. The signal is extracted using a fit to the diphoton invariant mass spectrum assuming that the width of the resonance is much smaller than the experimental resolution. The signal yields are corrected for the effects of detector inefficiency and resolution. The $pp\rightarrow H \rightarrow \gamma\gamma$ fiducial cross section is measured to be $43.2 \pm 9.4 (stat) {}^{+3.2}_{2.9} (syst) \pm 1.2 (lumi)$ fb for a Higgs boson of mass 125.4 GeV decaying to two isolated photons that have transverse momentum greater than 35% and 25% of the diphoton invariant mass and each with absolute pseudorapidity less than 2.37. Four additional fiducial cross sections and two crosssection limits are presented in phase space regions that test the theoretical modelling of different Higgs boson production mechanisms, or are sensitive to physics beyond the Standard Model. Differential cross sections are also presented, as a function of variables related to the diphoton kinematics and the jet activity produced in the Higgs boson events. The observed spectra are statistically limited but broadly in line with the theoretical expectations.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Memory effect in uniformly heated granular gases
E. Trizac ^{1} A. Prados ^{1, 2}
Physical Review E, American Physical Society, 2014, 90, pp.012204
We evidence a Kovacslike memory effect in a uniformly driven granular gas. A system of inelastic hard particles, in the low density limit, can reach a nonequilibrium steady state when properly forced. By following a certain protocol for the drive time dependence, we prepare the gas in a state where the granular temperature coincides with its long time value. The temperature subsequently does not remain constant, but exhibits a nonmonotonic evolution with either a maximum or a minimum, depending on the dissipation, and on the protocol. We present a theoretical analysis of this memory effect, at BoltzmannFokkerPlanck equation level, and show that when dissipation exceeds a threshold, the response can be coined anomalous. We find an excellent agreement between the analytical predictions and direct Monte Carlo simulations.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Fisica Teorica

Memory effects in avalanche dynamics: a key to the statistical properties of earthquakes
E. A. Jagla ^{1} François P. Landes ^{2} Alberto Rosso ^{2}
Physical Review Letters, American Physical Society, 2014, 112, pp.174301
Many complex systems respond to continuous input of energy by accumulation of stress over time and sudden energy releases in the form of avalanches. Avalanches are paradigmatic nonequilibrium phenomena displaying power law size distribution and involving all the length scales in the system. Conventional avalanche models disregard memory effects and thus miss basic features observed in real systems. Notable examples are aftershocks and the anomalous exponent of the GutenbergRichter law which characterize earthquake statistics. We propose a model which accounts for memory effects through the introduction of viscoelastic relaxation at an intermediate time scale. We demonstrate that in the resulting dynamics, coherent oscillations of the stress field emerge spontaneously without fine tuning of any parameter. Remarkably, in two dimensions, which is relevant in seismicity, these oscillations generate instability patterns that produce realistic earthquake dynamics with the correct GutenbergRichter exponent.
 1. Centro Atómico Bariloche and Instituto Balseiro
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Modeling of dislocations in a CDW junction: interference of the CDW and the normal carriers
Alvaro Rojo Bravo ^{1} Tianyou Yi ^{2} Natasha Kirova ^{3, 4} Serguei Brazovskii ^{1, 3}
International schoolconference on electronic crystals: ECRYS2014, Aug 2014, Cargèse, France. 460, pp.16, Physica B: Condensed Matter
We derive and study equations for dissipative transient processes in a constraint incommensurate charge density wave (CDW) with remnant pockets or a thermal population of normal carriers. The attention was paid to give the correct conservation of condensed and normal electrons, which was problematic at presence of moving dislocation cores if working within an intuitive GinzburgLandau like model. We performed a numeric modelling for stationary and transient states in a rectangular geometry when the voltage V or the normal current are applied across the conducting chains. We observe creation of an array of electronic vortices, the dislocations, at or close to the junction surface; their number increases stepwise with increasing V. The dislocation core strongly concentrates the normal carriers but the CDW phase distortions almost neutralize the total charge. At other regimes, the lines of the zero CDW amplitude flash across the sample working as phase slips. The studies were inspired by, and can be applied to experiments on mesajunctions in NbSe3 and TaS3 (Yu.I. Latyshev et al in proceedings of ECRYS 2008 and 2011).
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. South University of Science and Technology of China
 3. International Institute of Physics
 4. LPS  Laboratoire de Physique des Solides

Move ordering and communities in complex networks describing the game of go
Vivek Kandiah ^{1} Bertrand Georgeot ^{1} Olivier Giraud ^{2}
The European Physical Journal B, Springer, 2014, 87, pp.246. <10.1140/epjb/e2014504971>
We analyze the game of go from the point of view of complex networks. We construct three different directed networks of increasing complexity, defining nodes as local patterns on plaquettes of increasing sizes, and links as actual successions of these patterns in databases of real games. We discuss the peculiarities of these networks compared to other types of networks. We explore the ranking vectors and community structure of the networks and show that this approach enables to extract groups of moves with common strategic properties. We also investigate different networks built from games with players of different levels or from different phases of the game. We discuss how the study of the community structure of these networks may help to improve the computer simulations of the game. More generally, we believe such studies may help to improve the understanding of human decision process.
 1. LPT  Laboratoire de Physique Théorique  IRSAMC
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Multiple ordering transitions in a chiral liquid
Pierre Ronceray ^{1, 2} Peter Harrowell ^{3}
Physical Review Letters, American Physical Society, 2014, 112, pp.017801
We present here a numerical study of a lattice model of a chiral liquid. The low symmetry of the favoured local structure depresses the freezing point to reveal an exotic liquidliquid transition characterised by the appearance of an extended chirality,prior to freezing. What mechanisms impede crystallisation in liquids with low molecular symmetry ? The ordered liquid can be readily supercooled to zero temperature, as the combination of critical slowing down and competing crystal polymorphs results in a dramatically slow crystallisation process.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. ENS Paris  École normale supérieure  Paris
 3. Faculty of Sciences

Nearextreme eigenvalues and the first gap of Hermitian random matrices
Anthony Perret ^{1} Gregory Schehr ^{1}
Journal of Statistical Physics, Springer Verlag (Germany), 2014, 156, pp.843876
We study the phenomenon of "crowding" near the largest eigenvalue $\lambda_{\max}$ of random $N \times N$ matrices belonging to the Gaussian Unitary Ensemble (GUE) of random matrix theory. We focus on two distinct quantities: (i) the density of states (DOS) near $\lambda_{\max}$, $\rho_{\rm DOS}(r,N)$, which is the average density of eigenvalues located at a distance $r$ from $\lambda_{\max}$ and (ii) the probability density function of the gap between the first two largest eigenvalues, $p_{\rm GAP}(r,N)$. In the edge scaling limit where $r = {\cal O}(N^{1/6})$, which is described by a double scaling limit of a system of unconventional orthogonal polynomials, we show that $\rho_{\rm DOS}(r,N)$ and $p_{\rm GAP}(r,N)$ are characterized by scaling functions which can be expressed in terms of the solution of a Lax pair associated to the Painlevé XXXIV equation. This provides an alternative and simpler expression for the gap distribution, which was recently studied by Witte, Bornemann and Forrester in Nonlinearity 26, 1799 (2013). Our expressions allow to obtain precise asymptotic behaviors of these scaling functions both for small and large arguments.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Nonequilibrium statistical mechanics of the heat bath for two Brownian particles
Caterina De Bacco ^{1} Fulvio Baldovin ^{2} Enzo Orlandini ^{2} Ken Sekimoto ^{3, 4}
Physical Review Letters, American Physical Society, 2014, 112, pp.180605
We propose a new look at the heat bath for two Brownian particles, in which the heat bath as a 'system' is both perturbed and sensed by the Brownian particles. Nonlocal thermal fluctuation give rise to bathmediated static forces between the particles. Based on the general sumrule of the linear response theory, we derive an explicit relation linking these forces to the friction kernel describing the particles' dynamics. The relation is analytically confirmed in the case of two solvable models and could be experimentally challenged. Our results point out that the inclusion of the environment as a part of the whole system is important for micron or nanoscale physics.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Dipartimento di Fisica e Astronomia G. Galilei and Sezione I NFN
 3. MSC  Matière et Systèmes Complexes
 4. LPCT  Laboratoire de PhysicoChimie Théorique

Nonlocal response in a lattice gas under a shear drive
Tridib Sadhu ^{1} Satya N. Majumdar ^{2} David Mukamel ^{3}
Journal of Physics A: Mathematical and Theoretical, Institute of Physics: Hybrid Open Access, 2014, 47, pp.505005
In equilibrium, the effect of a spatially localised perturbation is typically confined around the perturbed region. Quite contrary to this, in a nonequilibrium stationary state often the entire system is affected. This appears to be a generic feature of nonequilibrium. We study such nonlocal response in the stationary state of a lattice gas with a shear drive at the boundary which keeps the system out of equilibrium. We show that a perturbation in the form of a localised blockage at the boundary, induces algebraically decaying density and current profile. In two examples, noninteracting particles and particles with simple exclusion, we analytically derive the powerlaw tail of the profiles.
 1. IPHT  Institut de Physique Théorique (ex SPhT)
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. Weizmann Institute

Nonequilibrium firstorder transition in coupled oscillator systems with inertia and noise
Shamik Gupta ^{1} Alessandro Campa ^{2} Stefano Ruffo ^{3}
Physical Review E, American Physical Society, 2014, 89, pp.022123
We study the dynamics of a system of coupled oscillators of distributed natural frequencies, by including the features of both thermal noise, parametrized by a temperature, and inertial terms, parametrized by a moment of inertia. For a general unimodal frequency distribution, we report here the complete phase diagram of the model in the space of dimensionless moment of inertia, temperature, and width of the frequency distribution. We demonstrate that the system undergoes a nonequilibrium firstorder phase transition from a synchronized phase at low parameter values to an incoherent phase at high values. We provide strong numerical evidence for the existence of both the synchronized and the incoherent phase, treating the latter analytically to obtain the corresponding linear stability threshold that bounds the firstorder transition point from below. In the limit of zero noise and inertia, when the dynamics reduces to the one of the Kuramoto model, we recover the associated known continuous transition. At finite noise and inertia but in the absence of natural frequencies, the dynamics becomes that of a wellstudied model of longrange interactions, the Hamiltonian meanfield model. Close to the firstorder phase transition, we show that the escape time out of metastable states scales exponentially with the number of oscillators, which we explain to be stemming from the longrange nature of the interaction between the oscillators.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Complex Systems and Theoretical Physic Unit
 3. Dipartimento di Fisica e Astronomia and CSDC

Nonlinear waves of polarization in twocomponent BoseEinstein condensates
A. M. Kamchatnov ^{1} Y. V. Kartashov ^{1} P. É. Larré ^{2} N. Pavloff ^{2}
Physical Review A, American Physical Society, 2014, 89, pp.033618
Waves with different symmetries exist in twocomponent BoseEinstein condensates (BECs) whose dynamics is described by a system of coupled GrossPitaevskii (GP) equations. A first type of waves corresponds to excitations for which the motion of both components is locally in phase. In the second type of waves the two components have a counterphase local motion. In the case of different values of inter and intracomponent interaction constants, the long wavelength behavior of these two modes corresponds to two types of sound with different velocities. In the limit of weak nonlinearity and small dispersion the first mode is described by the wellknown Kortewegde Vries (KdV) equation. We show that in the same limit the second mode can be described by the Gardner (modified KdV) equation, if the intracomponent interaction constants have close enough values. This leads to a rich phenomenology of nonlinear excitations (solitons, kinks, algebraic solitons, breathers) which does not exist in the KdV description.
 1. Institute of Spectroscopy
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Observation of a disordered bosonic insulator from weak to strong interactions
Chiara D Errico ^{1, 2} Eleonora Lucioni ^{1, 2} Luca Tanzi ^{1} Lorenzo Gori ^{1} Guillaume Roux ^{3} Ian P. McCulloch ^{4} Thierry Giamarchi ^{5} Massimo Inguscio ^{1, 6} Giovanni Modugno ^{1, 2}
Physical Review Letters, American Physical Society, 2014, 113, pp.095301
We employ ultracold atoms with controllable disorder and interaction to study the paradigmatic problem of disordered bosons in the full disorderinteraction plane. Combining measurements of coherence, transport and excitation spectra, we get evidence of an insulating regime extending from weak to strong interaction and surrounding a superfluidlike regime, in general agreement with the theory. For strong interaction, we reveal the presence of a stronglycorrelated Bose glass coexisting with a Mott insulator.
 1. European Laboratory for Nonlinear Spectroscopy (LENS)
 2. Istituto Nazionale di Ottica
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 4. Centre for Engineered Quantum Systems
 5. DPMC  Département de Physique de la Matière Condensée
 6. INRIM

Observation of a Strong AtomDimer Attraction in a MassImbalanced FermiFermi Mixture
Michael Jag ^{1, 2} Matteo Zaccanti ^{1, 3} Marko Cetina ^{1} Rianne S. Lous ^{1, 4} Florian Schreck ^{1} Rudolf Grimm ^{1, 4} Dmitry S. Petrov ^{5} Jesper Levinsen ^{6, 7}
Physical Review Letters, American Physical Society, 2014, 112, pp.075302
We investigate a mixture of ultracold fermionic $^{40}$K atoms and weakly bound $^{6}$Li$^{40}$K dimers on the repulsive side of a heteronuclear atomic Feshbach resonance. By radiofrequency spectroscopy we demonstrate that the normally repulsive atomdimer interaction is turned into a strong attraction. The phenomenon can be understood as a threebody effect in which two heavy $^{40}$K fermions exchange the light $^{6}$Li atom, leading to attraction in odd partialwave channels (mainly pwave). Our observations show that mass imbalance in a fermionic system can profoundly change the character of interactions as compared to the wellestablished massbalanced case.
 1. Institut für Quantenoptik und Quanteninformation
 2. Institut fur Experimentalphysik
 3. CNR  Istituto Nazionale Ottica
 4. Institut fur Experimentalphysik
 5. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 6. TCM Group, Cavendish Laboratory
 7. Aarhus Institute of Advanced Studies

On the density of shear transformation zones in amorphous solids
Jie Lin ^{1} Alaa Saade ^{1} Edan Lerner ^{1} Alberto Rosso ^{2} Matthieu Wyart ^{1}
Europhysics Letters, EDP Science, 2014, 105, pp.26003
We study the stability of amorphous solids, focusing on the distribution P(x) of the local stress increase x that would lead to an instability. We argue that this distribution is singular P(x)x^{\theta}, where the exponent {\theta} is nonzero if the elastic interaction between rearranging regions is nonmonotonic, and increases with the interaction range. For a class of finite dimensional models we show that stability implies a lower bound on {\theta}, which is found to lie near saturation. For quadrupolar interactions these models yield {\theta} ~ 0.6 for d=2 and \theta ~ 0.4 in d=3 where d is the spatial dimension, accurately capturing previously unresolved observations in atomistic models, both in quasistatic flow and after a fast quench.
 1. New York University, Center for Soft Matter Research
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

On the Gap and Time Interval between the First Two Maxima of Long Random Walks
Satya N. Majumdar ^{1} Philippe Mounaix ^{2} Gregory Schehr ^{1}
Journal of Statistical Mechanics, 2014, pp.P09013
In the context of order statistics of discrete time random walks (RW), we investigate the statistics of the gap, $G_n$, and the number of time steps, $L_n$, between the two highest positions of a Markovian onedimensional random walker, starting from $x_0 = 0$, after $n$ time steps (taking the $x$axis vertical). The jumps $\eta_i = x_i  x_{i1}$ are independent and identically distributed random variables drawn from a symmetric probability distribution function (PDF), $f(\eta)$, the Fourier transform of which has the small $k$ behavior $1  \hat f(k) \propto k^\mu$, with $0 < \mu \leq 2$. For $\mu=2$, the variance of the jump distribution is finite and the RW (properly scaled) converges to a Brownian motion. For $0<\mu<2$, the RW is a L\'evy flight of index $\mu$. We show that the joint PDF of $G_n$ and $L_n$ converges to a well defined stationary bivariate distribution $p(g,l)$ as the RW duration $n$ goes to infinity. We present a thorough analytical study of the limiting joint distribution $p(g,l)$, as well as of its associated marginals $p_{\rm gap}(g)$ and $p_{\rm time}(l)$, revealing a rich variety of behaviors depending on the tail of $f(\eta)$ (from slow decreasing algebraic tail to fast decreasing superexponential tail). We also address the problem for a random bridge where the RW starts and ends at the origin after $n$ time steps. We show that in the large $n$ limit, the PDF of $G_n$ and $L_n$ converges to the {\it same} stationary distribution $p(g,l)$ as in the case of the freeend RW. Finally, we present a numerical check of our analytical predictions. Some of these results were announced in a recent letter [S. N. Majumdar, Ph. Mounaix, G. Schehr, Phys. Rev. Lett. {\bf 111}, 070601 (2013)].
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. CPHT  Centre de Physique Théorique [Palaiseau]

Onedimensional disordered quantum mechanics and Sinai diffusion with random absorbers
Aurélien Grabsch ^{1} Christophe Texier ^{2} Yves Tourigny ^{3}
Journal of Statistical Physics, Springer Verlag (Germany), 2014, 155, pp.237276
We study the onedimensional Schrödinger equation with a disordered potential of the form $V (x) = \phi(x)^2+\phi'(x) + \kappa(x) $ where $\phi(x)$ is a Gaussian white noise with mean $\mu g$ and variance $g$, and $\kappa(x)$ is a random superposition of delta functions distributed uniformly on the real line with mean density $\rho$ and mean strength $v$. Our study is motivated by the close connection between this problem and classical diffusion in a random environment (the Sinai problem) in the presence of random absorbers~: $\phi(x)$ models the force field acting on the diffusing particle and $\kappa(x)$ models the absorption properties of the medium in which the diffusion takes place. The focus is on the calculation of the complex Lyapunov exponent $ \Omega(E) = \gamma(E)  \mathrm{i} \pi N(E) $, where $N$ is the integrated density of states per unit length and $\gamma$ the reciprocal of the localisation length. By using the continuous version of the DysonSchmidt method, we find an exact formula, in terms of a Hankel function, in the particular case where the strength of the delta functions is exponentiallydistributed with mean $v=2g$. Building on this result, we then solve the general case  in the lowenergy limit  in terms of an infinite sum of Hankel functions. Our main result, valid without restrictions on the parameters of the model, is that the integrated density of states exhibits the power law behaviour $$ N(E) \underset{E\to0+}{\sim} E^\nu \hspace{0.5cm} \mbox{where } \nu=\sqrt{\mu^2+2\rho/g}\:. $$ This confirms and extends several results obtained previously by approximate methods.
 1. ENS Cachan  École normale supérieure  Cachan
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. School of Mathematics [Bristol]

Outage Capacity for the Optical MIMO Channel
Apostolos Karadimitrakis ^{1} Aris L. Moustakas ^{1, 2} Pierpaolo Vivo ^{3}
France. 60, pp.4370, 2014, Information theory IEEE Transactions on
MIMO processing techniques in fiber optical communications have been proposed as a promising approach to meet increasing demand for information throughput. In this context, the multiple channels correspond to the multiple modes and/or multiple cores in the fiber. In this paper we characterize the distribution of the mutual information with Gaussian input in a simple channel model for this system. Assuming significant cross talk between cores, negligible backscattering and nearlossless propagation in the fiber, we model the transmission channel as a random complex unitary matrix. The loss in the transmission may be parameterized by a number of unutilized channels in the fiber. We analyze the system in a dual fashion. First, we evaluate a closedform expression for the outage probability, which is handy for small matrices. We also apply the asymptotic approach, in particular the Coulomb gas method from statistical mechanics, to obtain closedform results for the ergodic mutual information, its variance as well as the outage probability for Gaussian input in the limit of large number of cores/modes. By comparing our analytic results to simulations, we see that, despite the fact that this method is nominally valid for large number of modes, our method is quite accurate even for small to modest number of channels.
 1. UNIV. ATHENS  University of Athens, Greece
 2. E3S  Supélec Sciences des Systèmes  EA4454
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Phase transitions and edge scaling of number variance in Gaussian random matrices
Ricardo Marino ^{1} Satya N. Majumdar ^{1} Grégory Schehr ^{1} Pierpaolo Vivo ^{1}
Physical Review Letters, American Physical Society, 2014, 112, pp.254101
We consider $N\times N$ Gaussian random matrices, whose average density of eigenvalues has the Wigner semicircle form over $[\sqrt{2},\sqrt{2}]$. For such matrices, using a Coulomb gas technique, we compute the large $N$ behavior of the probability $\mathcal{P}_{\scriptscriptstyle N,L}(N_L)$ that $N_L$ eigenvalues lie within the box $[L,L]$. This probability scales as $\mathcal{P}_{\scriptscriptstyle N,L}(N_L=\kappa_L N)\approx\exp\left({\beta} N^2 \psi_L(\kappa_L)\right)$, where $\beta$ is the Dyson index of the ensemble and $\psi_L(\kappa_L)$ is a $\beta$independent rate function that we compute exactly. We identify three regimes as $L$ is varied: (i) $\, N^{1}\ll L<\sqrt{2}$ (bulk), (ii) $\ L\sim\sqrt{2}$ on a scale of $\mathcal{O}(N^{{2}/{3}})$ (edge) and (iii) $\ L > \sqrt{2}$ (tail). We find a dramatic nonmonotonic behavior of the number variance $V_N(L)$ as a function of $L$: after a logarithmic growth $\propto \ln (N L)$ in the bulk (when $L \sim {\cal O}(1/N)$), $V_N(L)$ decreases abruptly as $L$ approaches the edge of the semicircle before it decays as a stretched exponential for $L > \sqrt{2}$. This "dropoff" of $V_N(L)$ at the edge is described by a scaling function $\tilde V_{\beta}$ which smoothly interpolates between the bulk (i) and the tail (iii). For $\beta = 2$ we compute $\tilde V_2$ explicitly in terms of the Airy kernel. These analytical results, verified by numerical simulations, directly provide for $\beta=2$ the full statistics of particlenumber fluctuations at zero temperature of 1d spinless fermions in a harmonic trap.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Phase transitions in the condition number distribution of Gaussian random matrices
Isaac Pérez Castillo ^{1} Eytan Katzav ^{2} Pierpaolo Vivo ^{3, 4}
Physical Review E, American Physical Society, 2014, 90, pp.050103
We study the statistics of the condition number $\kappa=\lambda_{\mathrm{max}}/\lambda_{\mathrm{min}}$ (the ratio between largest and smallest squared singular values) of $N\times M$ Gaussian random matrices. Using a Coulomb fluid technique, we derive analytically and for large $N$ the cumulative $\mathcal{P}[\kappa
 1. Departamento de Sistemas Complejos, Instituto de Física, UNAM
 2. Racah Institute of Physics  Racah Institute of Physics
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 4. Department of Mathematics, King's College London

Polydispersity and optimal relaxation in the hard sphere fluid
Matthieu Barbier ^{1} Emmanuel Trizac ^{1}
Journal of Statistical Physics, Springer Verlag (Germany), 2014, 154, pp.13651390
We consider the mass heterogeneity in a gas of polydisperse hard particles as a key to optimizing a dynamical property: the kinetic relaxation rate. Using the framework of the Boltzmann equation, we study the long time approach of a perturbed velocity distribution toward the equilibrium Maxwellian solution. We work out the cases of discrete as well as continuous distributions of masses, as found in dilute fluids of mesoscopic particles such as granular matter and colloids. On the basis of analytical and numerical evidence, we formulate a dynamical equipartition principle that leads to the result that no such continuous dispersion in fact minimizes the relaxation time, as the global optimum is characterized by a finite number of species. This optimal mixture is found to depend on the dimension d of space, ranging from five species for d=1 to a single one for d>=4. The role of the collisional kernel is also discussed, and extensions to dissipative systems are shown to be possible.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Quantum Flutter: Signatures and Robustness
Michael Knap ^{1, 2, 3} Charles J. M. Mathy ^{3, 4} Martin Ganahl ^{2} Mikhail B. Zvonarev ^{5} Eugene Demler ^{4}
Physical Review Letters, American Physical Society, 2014, 112, pp.015302
We investigate the motion of an impurity particle injected with finite velocity into an interacting onedimensional quantum gas. Using largescale numerical simulations based on matrix product states, we observe and quantitatively analyze longlived oscillations of the impurity momentum around a nonzero saturation value, called quantum flutter. We show that the quantum flutter frequency is equal to the energy difference between two branches of collective excitations of the model. We propose an explanation of the finite saturation momentum of the impurity based on the properties of the edge of the excitation spectrum. Our results indicate that quantum flutter exists away from integrability, and provide parameter regions in which it could be observed in experiments with ultracold atoms using currently available technology.
 1. Department of Physics
 2. Institute of Theoretical and Computational Physics
 3. CFA  HarvardSmithsonian Center for Astrophysics
 4. Department of Physics
 5. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Quasiparticle spectroscopy as a probe of the topological phase in graphene with heavy adatoms
Paul Soulé ^{1, 2} Marcel Franz ^{2, 3}
4 pages, 2 figures; Version 2 : references added, minor changes. 2014
Electrons in graphene with heavy adatoms (such as In or Tl) have been predicted to form a 2D topological insulator phase with a substantial spectral gap potentially suitable for future practical applications. In order to facilitate the ongoing experimental efforts to identify this phase we perform a theoretical study of its spectral properties in a model graphene system with randomly distributed adatoms. Our extensive modeling shows that random heavy adatoms produce a full spectral gap (as opposed to a mobility gap) accompanied by distinctive quasiparticle interference patterns observable by means of Fouriertransform scanning tunneling spectroscopy.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Department of Physics and Astronomy
 3. Quantum Matter Institute

Relevance of electron spin dissipative processes on Dynamic Nuclear Polarization via Thermal Mixing
Sonia Colombo Serra ^{1} Marta Filibian ^{2} Pietro Carretta ^{3} Alberto Rosso ^{4} Fabio Tedoldi ^{1}
Physical Chemestry Chemical Physics, 2014, 16, pp.753764
The available theoretical approaches aiming at describing Dynamic Nuclear spin Polarization (DNP) in solutions containing molecules of biomedical interest and paramagnetic centers are not able to model the behaviour observed upon varying the concentration of trityl radicals or the polarization enhancement caused by moderate addition of gadolinium complexes. In this manuscript, we first show experimentally that the nuclear steady state polarization reached in solutions of pyruvic acid with 15 mM trityl radicals is substantially independent from the average internuclear distance. This evidences a leading role of electron (over nuclear) spin relaxation processes in determining the ultimate performances of DNP. Accordingly, we have devised a variant of the Thermal Mixing model for inhomogenously broadened electron resonance lines which includes a relaxation term describing the exchange of magnetic anisotropy energy of the electron spin system with the lattice. Thanks to this additional term, the dependence of the nuclear polarization on the electron concentration can be properly accounted for. Moreover, the model predicts a strong increase of the final polarization on shortening the electron spinlattice relaxation time, providing a possible explanation for the effect of gadolinium doping.
 1. Centro Ricerche Bracco
 2. CNISM  Dipartimento di Fisica and Unità
 3. CNISM  Dipartimento di Fisica and Unità
 4. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Reweighted belief propagation and quiet planting for random KSAT
Florent Krzakala ^{1} Marc Mézard ^{2, 3} Lenka Zdeborová ^{4}
Journal on Satisfiability, Boolean Modeling and Computation, IOS Press, 2014, 8, pp.149
We study the random Ksatisfiability problem using a partition function where each solution is reweighted according to the number of variables that satisfy every clause. We apply belief propagation and the related cavity method to the reweighted partition function. This allows us to obtain several new results on the properties of random Ksatisfiability problem. In particular the reweighting allows to introduce a planted ensemble that generates instances that are, in some region of parameters, equivalent to random instances. We are hence able to generate at the same time a typical random SAT instance and one of its solutions. We study the relation between clustering and belief propagation fixed points and we give a direct evidence for the existence of purely entropic (rather than energetic) barriers between clusters in some region of parameters in the random Ksatisfiability problem. We exhibit, in some large planted instances, solutions with a nontrivial whitening core; such solutions were known to exist but were so far never found on very large instances. Finally, we discuss algorithmic hardness of such planted instances and we determine a region of parameters in which planting leads to satisfiable benchmarks that, up to our knowledge, are the hardest known.
 1. LPCT  Laboratoire de PhysicoChimie Théorique
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 3. ENS  Ecole Normale Supérieure
 4. IPHT  Institut de Physique Théorique (ex SPhT)

Scaling description of the yielding transition in soft amorphous solids at zero temperature
Jie Lin ^{1, 2} Edan Lerner ^{3} Alberto Rosso ^{4} Matthieu Wyart ^{3}
Proceedings of the National Academy of Sciences U S A, National Academy of Sciences, 2014, 111, pp.14382
Yield stress materials flow if a sufficiently large shear stress is ap plied. Although such materials are ubiquitous and relevant for indus try, there is no accepted microscopic description of how they yield, even in the simplest situations where temperature is negligible and where flow inhomogeneities such as shear bands or fractures are ab sent. Here we propose a scaling description of the yielding transition in amorphous solids made of soft particles at zero temperature. Our description makes a connection between the HerschelBulkley expo nent characterizing the singularity of the flow curve near the yield stress {\Sigma}c, the extension and duration of the avalanches of plasticity observed at threshold, and the density P(x) of soft spots, or shear transformation zones, as a function of the stress increment x be yond which they yield. We argue that the critical exponents of the yielding transition can be expressed in terms of three independent exponents {\theta}, df and z, characterizing respectively the density of soft spots, the fractal dimension of the avalanches, and their duration. Our description shares some similarity with the depinning transition that occurs when an elastic manifold is driven through a random potential, but also presents some striking differences. We test our arguments in an elastoplastic model, an automaton model similar to those used in depinning, but with a different interaction kernel, and find satisfying agreement with our predictions both in two and three dimensions.
 1. Division of Parasitology
 2. Parasitology, Center of Infectious Diseases
 3. New York University, Center for Soft Matter Research
 4. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Scaling properties of fieldinduced superdiffusion in Continous Time Random Walks
R. Burioni ^{1, 2} G. Gradenigo ^{3, 4} A. Sarracino ^{3} A. Vezzani ^{5, 6} A. Vulpiani ^{7}
Small system nonequilibrium fluctuations, dynamics and stochastics, and anomalous behavior, China. 62, pp.514, 2014, Communications in Theoretical Physics
We consider a broad class of Continuous Time Random Walks with large fluctuations effects in space and time distributions: a random walk with trapping, describing subdiffusion in disordered and glassy materials, and a L\'evy walk process, often used to model superdiffusive effects in inhomogeneous materials. We derive the scaling form of the probability distributions and the asymptotic properties of all its moments in the presence of a field by two powerful techniques, based on matching conditions and on the estimate of the contribution of rare events to powerlaw tails in a field.
 1. Fisica
 2. INFN
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 4. CEA  CEA Saclay
 5. Fisica
 6. Fisica
 7. INFN  Istituto Nazionale di Fisica Nucleare [Sezione di Roma 1]

Shortest nodedisjoint paths on random graphs
Caterina De Bacco ^{1} Silvio Franz ^{1} David Saad ^{2} Chi Ho Yeung ^{2, 3, 4}
Journal of Statistical Mechanics, 2014, pp.P07009
A localized method to distribute paths on random graphs is devised, aimed at finding the shortest paths between given source/destination pairs while avoiding path overlaps at nodes. We propose a method based on messagepassing techniques to process global information and distribute paths optimally. Statistical properties such as scaling with system size and number of paths, average pathlength and the transition to the frustrated regime are analysed. The performance of the suggested algorithm is evaluated through a comparison against a greedy algorithm.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Aston University
 3. Department of Physics
 4. Department of Science and Environmental Studies

SO(5) symmetry in the quantum Hall effect in graphene
Fengcheng Wu ^{1} Inti Sodemann ^{1} Yasufumi Araki ^{1} Allan H. MacDonald ^{1} Thierry Jolicoeur ^{2}
Physical Review B (Condensed Matter), American Physical Society, 2014, pp.235432
Electrons in graphene have four flavors associated with lowenergy spin and valley degrees of freedom. The fractional quantum Hall effect in graphene is dominated by longrange Coulomb interactions which are invariant under rotations in spinvalley space. This SU(4) symmetry is spontaneously broken at most filling factors, and also weakly broken by atomic scale valleydependent and valleyexchange interactions with coupling constants $g_{z}$ and $g_{\perp}$. In this paper we demonstrate that when $g_{z}=g_{\perp}$ an exact SO(5) symmetry survives which unifies the N\'eel spin order parameter of the antiferromagnetic state and the $XY$ valley order parameter of the Kekul\'e distortion state into a single fivecomponent order parameter. The proximity of the highly insulating quantum Hall state observed in graphene at $\nu=0$ to an ideal SO(5) symmetric quantum Hall state remains an open experimental question. We illustrate the physics associated with this SO(5) symmetry by studying the multiplet structure and collective dynamics of filling factor $\nu=0$ quantum Hall states based on exactdiagonalization and lowenergy effective theory approaches. This allows to illustrate how manifestations of the SO(5) symmetry would survive even when it is weakly broken.
 1. University of Texas at Austin
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Spectral density of the quantum Ising model in two fields: Gaussian and multiGaussian approximations
Y. Y. Atas ^{1} E. Bogomolny ^{1}
Journal of Physics A: Mathematical and Theoretical, Institute of Physics: Hybrid Open Access, 2014, 47, pp.335201
Spectral density of quantum Ising model in two fields for large but finite number of spins $N$, is discussed in detail. When all coupling constants are of the same order, spectral densities in the bulk are well approximated by a Gaussian function which is typical behaviour for manybody models with shortrange interactions. The main part of the paper is devoted to the investigation of a different characteristic case when spectral densities have peaks related with strong degeneracies of unperturbed states in certain limits of coupling constants. In the strict limit $N\to\infty$, peaks overlap and disappear but for values of $N$ accessible in numerical calculations they often strongly influence spectral densities and other quantities as well. A simple method is developed which permits to find general approximation formulae for multipeak structure of spectral density in good agreement with numerics.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Spinon excitation spectra of the $J_1$$J_2$ chain from analytical calculations in the dimer basis and exact diagonalization
Lavarélo Arthur ^{1} Roux Guillaume ^{1}
The European Physical Journal B, Springer, 2014, 87, pp.229
The excitation spectrum of the frustrated spin$1/2$ Heisenberg chain is reexamined using variational and exact diagonalization calculations. We show that the overlap matrix of the shortrange resonating valence bond states basis can be inverted which yields tractable equations for single and two spinons excitations. Older results are recovered and new ones, such as the bondstate dispersion relation and its size with momentum at the MajumdarGhosh point are found. In particular, this approach yields a gap opening at $J_2=0.25J_1$ and an onset of incommensurability in the dispersion relation at $J_2=9/17J_1$ [as in S. Brehmer \emph{et al.}, J. Phys.: Condens. Matter \textbf{10}, 1103 (1998)]. These analytical results provide a good support for the understanding of exact diagonalization spectra, assuming an independent spinons picture.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Statistics of the first passage time of Brownian motion conditioned by maximum value or area
Michael J. Kearney ^{1} Satya N. Majumdar ^{2}
Journal of Physics A: Mathematical and Theoretical, Institute of Physics: Hybrid Open Access, 2014, 47, pp.465001
We derive the moments of the first passage time for Brownian motion conditioned by either the maximum value or the area swept out by the motion. These quantities are the natural counterparts to the moments of the maximum value and area of Brownian excursions of fixed duration, which we also derive for completeness within the same mathematical framework. Various applications are indicated.
 1. Faculty of Engineering and Physical Sciences
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Strong pinning of propagation fronts in adverse flow
Thomas Gueudré ^{1} Awadhesh Kumar Dubey ^{2} Laurent Talon ^{2} Alberto Rosso ^{3}
Physical Review E, American Physical Society, 2014, 89, pp.041004
Reaction fronts evolving in a porous medium exhibit a rich dynamical behaviour. In presence of an adverse flow, experiments show that the front slows down and eventually gets pinned, displaying a particular sawtooth shape. Extensive numerical simulations of the hydrodynamic equations confirm the experimental observations. Here we propose a stylized model, predicting two possible outcomes of the experiments for large adverse flow: either the front develops a sawtooth shape, or it acquires a complicated structure with islands and overhangs. A simple criterion allows to distinguish between the two scenarios and its validity is reproduced by direct hydrodynamical simulations. Our model gives a better understanding of the transition and is relevant in a variety of domains, when the pinning regime is strong and only relies on a small number of sites.
 1. LPTENS  Laboratoire de Physique Théorique de l'ENS
 2. FAST  Fluides, automatique, systèmes thermiques
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Synchronization transitions in globally coupled rotors in presence of noise and inertia: Exact results
Maxim Komarov ^{1, 2} Shamik Gupta ^{3} Arkady Pikovsky ^{1, 2}
Europhysics Letters, EDP Science, 2014, 106, pp.40003
We study a generic model of globally coupled rotors that includes the effects of noise, phase shift in the coupling, and distributions of moments of inertia and natural frequencies of oscillation. As particular cases, the setup includes previously studied SakaguchiKuramoto, Hamiltonian and Brownian meanfield, and TanakaLichtenbergOishi and AcebrónBonillaSpigler models. We derive an exact solution of the selfconsistent equations for the order parameter in the stationary state, valid for arbitrary parameters in the dynamics, and demonstrate nontrivial phase transitions to synchrony that include reentrant synchronous regimes.
 1. UNIPOTSDAM  Department of Physics and Astronomie, Univ. Potsdam
 2. Department of Control Theory
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Taking one charge off a twodimensional Wigner crystal
Moritz Antlanger ^{1, 2} Martial Mazars ^{1} Ladislav Šamaj ^{3, 4} Gerhard Kahl ^{2} Emmanuel Trizac ^{4}
LPTORSAY 1410. 17 pages, 11 figures. 2014
A planar array of identical charges at vanishing temperature forms a Wigner crystal with hexagonal symmetry. We take off one (reference) charge in a perpendicular direction, hold it fixed, and search for the ground state of the whole system. The planar projection of the reference charge should then evolve from a sixfold coordination (center of a hexagon) for small distances to a threefold arrangement (center of a triangle), at large distances $d$ from the plane. The aim of this paper is to describe the corresponding nontrivial lattice transformation. For that purpose, two numerical methods (direct energy minimization and Monte Carlo simulations), together with an analytical treatment, are presented. Our results indicate that the $d=0$ and $d\to\infty$ limiting cases extend for finite values of $d$ from the respective starting points into two sequences of stable states, with intersecting energies at some value $d_t$; beyond this value the branches continue as metastable states.
 1. LPT  Laboratoire de Physique Théorique d'Orsay [Orsay]
 2. Institut für Theoretische Physik and Center for Computational Materials Science
 3. Institute of Physics
 4. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

The effect of disorder geometry on the critical force in disordered elastic systems
V. Démery ^{1} V. Lecomte ^{2} A. Rosso ^{3}
Journal of Statistical Mechanics : Theory and Experiment, 2014, P03009, 20 p. <10.1088/17425468/2014/03/P03009>
 1. IJLRA  Institut Jean Le Rond d'Alembert
 2. LPMA  Laboratoire de Probabilités et Modèles Aléatoires
 3. Laboratoire Physique Théorique et Modèles Statistiques (UMR CNRS 8626), Université de ParisSud, Orsay Cedex,

The role of the glassy dynamics and thermal mixing in the dynamic nuclear polarization and relaxation mechanisms of pyruvic acid
Marta Filibian ^{1} Sonia Colombo Serra ^{2} Marco Moscardini ^{1} Alberto Rosso ^{3} Fabio Tedoldi ^{2} Pietro Carretta ^{1}
Physical Chemistry Chemical Physics, Royal Society of Chemistry, 2014, 16, pp.27025
The temperature dependence of $^1$H and $^{13}$C nuclear spinlattice relaxation rate $1/T_1$ has been studied between 4.2 K and 1.6 K in pure pyruvic acid and in pyruvic acid containing trityl radicals at a concentration of 15 mM. The temperature dependence of $1/T_1$ is found to follow a quadratic power law for both nuclei in the two samples. Remarkably that is the same temperature dependence found also for the electron spinlattice relaxation rate in the sample containing radicals. Dynamical nuclear polarization experiments on this latter sample show that below 4 K $^{13}$C build up rate and $1/T_1$ scale with $1/T_{1e}$ and their values are quantitatively consistent with the presence of a thermal mixing regime between the nuclear and electron spin reservoirs. These results are explained by considering the effect of the pyruvic acid glassy dynamics on the relaxation rates and by assuming that below 4 K dynamical nuclear polarization is driven by a very good thermal contact between the nuclear and electron spin reservoirs.
 1. CNISM  Dipartimento di Fisica and Unità
 2. Centro Ricerche Bracco
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

The threebody parameter for Efimov states in lithium6
Huang Bo ^{1} O'Hara Kenneth M. ^{2, 3} Grimm Rudolf ^{1, 3} Hutson Jeremy M. ^{4} Petrov Dmitry S. ^{5}
Physical Review A, American Physical Society, 2014, 90, pp.043636
We present a stateoftheart reanalysis of experimental results on Efimov resonances in the threefermion system of $^6$Li. We discuss different definitions of the 3body parameter (3BP) for Efimov states, and adopt a definition that excludes effects due to deviations from universal scaling for lowlying states. We develop a finitetemperature model for the case of three distinguishable fermions and apply it to the excitedstate Efimov resonance to obtain the most accurate determination to date of the 3BP in an atomic threebody system. Our analysis of groundstate Efimov resonances in the same system yields values for the threebody parameter that are consistent with the excitedstate result. Recent work has suggested that the reduced 3BP for atomic systems is a nearuniversal quantity, almost independent of the particular atom involved. However, the value of the 3BP obtained for $^6$Li is significantly ($\sim 20$%) different from that previously obtained from the excitedstate resonance in Cs. The difference between these values poses a challenge for theory.
 1. Institut für Experimentalphysik
 2. Pennsylvania State University
 3. ICOQI  Institut für Quantenoptik und Quanteninformation
 4. JQC  Joint Quantum Centre
 5. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

The two scenarios for quantum multifractality breakdown
Rémy Dubertrand ^{1} Ignacio GarcíaMata ^{2} Bertrand Georgeot ^{1} Olivier Giraud ^{3} Gabriel Lemarié ^{1} J. Martin ^{4}
Physical Review Letters, American Physical Society, 2014, 112, pp.234101. <10.1103/PhysRevLett.112.234101>
We expose two scenarios for the breakdown of quantum multifractality under the effect of perturbations. In the first scenario, multifractality survives below a certain scale of the quantum fluctuations. In the other one, the fluctuations of the wave functions are changed at every scale and each multifractal dimension smoothly goes to the ergodic value. We use as generic examples a onedimensional dynamical system and the threedimensional Anderson model at the metalinsulator transition. Our results suggest that quantum multifractality breakdown is universal and follows one of these two scenarios depending on the perturbation. We also discuss the experimental implications.
 1. Information et Chaos Quantiques (LPT)
 2. Instituto de Investigaciones Físicas de Mar del Plata
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 4. Institut de Physique Nucléaire, Atomique et de Spectroscopie

Thinking outside the box: fluctuations and finite size effects
Dario Villamaina ^{1} Emmanuel Trizac ^{1}
European Journal of Physics, Institute of Physics: Hybrid Open Access, 2014, 35, pp.3
The isothermal compressibility of an interacting or non interacting system may be extracted from the fluctuations of the number of particles in a well chosen control volume. Finite size effects are prevalent and should then be accounted for to obtain a meaningful, thermodynamic compressibility. In the traditional computational setup where a given simulation box is replicated with periodic boundary conditions, we study particle number fluctuations outside the box (i.e. when the control volume exceeds the box itself), which bear relevant thermodynamic information. We also investigate the related problem of extracting the compressibility from the structure factor in the small wavevector limit ($k\to 0$). The calculation should be restricted to the discrete set of wavevectors $k$ that are compatible with the periodicity of the system, and we assess the consequences of considering other $k$ values, a widespread error among beginners.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Thomas precession, persistent spin currents and quantum forces
Stephane Ouvry ^{1} Leonid Pastur ^{2} Andrey Yanovsky ^{2}
Europhysics Letters, EDP Science, 2014, 106, pp.18001
We consider Tinvariant spin currents induced by spinorbit interactions which originate from the confined motion of spin carriers in nanostructures. The resulting Thomas spin precession is a fundamental and purely kinematic relativistic effect occurring when the acceleration of carriers is not parallel to their velocity. In the case, where the carriers (e.g. electrons) have magnetic moment the forces due to the electric field of the spin current can, in certain conditions, exceed the van der WaalsCasimir forces by several orders of magnitude. We also discuss a possible experimental setup tailored to use these forces for checking the existence of a nonzero anomalous magnetic moment of the photon.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. Low Temperature Physics Institute, Kharkiv

ThreeBody Interacting Bosons in Free Space
D. S. Petrov ^{1}
Physical Review Letters, American Physical Society, 2014, 112, pp.103201
We propose a method of controlling two and threebody interactions in an ultracold Bose gas in any dimension. The method requires us to have two coupled internal singleparticle states split in energy such that the upper state is occupied virtually but amply during collisions. By varying system parameters one can switch off the twobody interaction while maintaining a strong threebody one. The mechanism can be implemented for dipolar bosons in the bilayer configuration with tunnelling or in an atomic system by using radiofrequency fields to couple two hyperfine states. One can then aim to observe a purely threebodyinteracting gas, dilute selftrapped droplets, the paired superfluid phase, Pfaffian state, and other exotic phenomena.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Top eigenvalue of a random matrix: large deviations and third order phase transition
Satya N. Majumdar ^{1} Gregory Schehr ^{1}
Journal of Statistical Mechanics: Theory and Experiment, IOP Science, 2014, pp.P01012
We study the fluctuations of the largest eigenvalue $\lambda_{\max}$ of $N \times N$ random matrices in the limit of large $N$. The main focus is on Gaussian $\beta$ensembles, including in particular the Gaussian orthogonal ($\beta=1$), unitary ($\beta=2$) and symplectic ($\beta = 4$) ensembles. The probability density function (PDF) of $\lambda_{\max}$ consists, for large $N$, of a central part described by TracyWidom distributions flanked, on both sides, by two large deviations tails. While the central part characterizes the typical fluctuations of $\lambda_{\max}$  of order ${\cal O}(N^{2/3})$ , the large deviations tails are instead associated to extremely rare fluctuations  of order ${\cal O}(1)$. Here we review some recent developments in the theory of these extremely rare events using a Coulomb gas approach. We discuss in particular the thirdorder phase transition which separates the left tail from the right tail, a transition akin to the socalled GrossWittenWadia phase transition found in 2d lattice quantum chromodynamics. We also discuss the occurrence of similar thirdorder transitions in various physical problems, including nonintersecting Brownian motions, conductance fluctuations in mesoscopic physics and entanglement in a bipartite system.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Topological transition in disordered planar matching: combinatorial arcs expansion
Andrey Y. Lokhov ^{1} Olga V. Valba ^{2, 3} Sergei K. Nechaev ^{1, 4, 3} Mikhail V. Tamm ^{5, 3}
Journal of Statistical Mechanics: Theory and Experiment, IOP Science, 2014, pp.P12004
In this paper, we investigate analytically the properties of the disordered Bernoulli model of planar matching. This model is characterized by a topological phase transition, yielding complete planar matching solutions only above a critical density threshold. We develop a combinatorial procedure of arcs expansion that explicitly takes into account the contribution of short arcs, and allows to obtain an accurate analytical estimation of the critical value by reducing the global constrained problem to a set of local ones. As an application to a toy representation of the RNA secondary structures, we suggest generalized models that incorporate a onetoone correspondence between the contact matrix and the RNAtype sequence, thus giving sense to the notion of effective noninteger alphabets.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. N. Semenov Institute of Chemical Physics
 3. Department of Applied Mathematics
 4. P. N. Lebedev Physical Institute
 5. Physics Department

Ultrafast switching to a stable hidden topologically protected quantum state in an electronic crystal
L. Stojchevska ^{1, 2} I. Vaskivskyi ^{1} T. Mertelj ^{1} P. Kusar ^{1} D. Svetin ^{1} S. Brazovskii ^{3, 4} D. Mihailovic ^{1, 2, 5}
Science, American Association for the Advancement of Science, 2014, 344, pp.177180
Hidden states of matter with novel and unusual properties may be created if a system out of equilibrium can be induced to follow a trajectory to a state which is inaccessible or does not even exist under normal equilibrium conditions. Here we report on the discovery of a hidden (H) topologically protected electronic state in a layered dichalcogenide 1TTaS2 crystal reached as a result of a quench caused by a single 35 fs laser pulse. The properties of the H state are markedly different from any other state of the system: it exhibits a large drop of electrical resistance, strongly modified single particle and collective mode spectra and a marked change of optical reflectivity. Particularly important and unusual, the H state is stable for an arbitrarily long time until a laser pulse, electrical current or thermal erase procedure is applied, causing it to revert to the thermodynamic ground state. Major observed events can be reproduced by a kinetic model describing the conversion of photo excited electrons and holes into an electronically ordered crystal, thus converting a Mott insulator to a conducting H state. Its longtime stability follows from the topological protection of the number of periods in the electronic crystal.
 1. Department of Complex Matter
 2. Jozef Stefan International Postgraduate School
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 4. International Institute of Physics
 5. CENN Nanocenter

Universal covariance formula for linear statistics on random matrices
Fabio Deelan Cunden ^{1, 2} Pierpaolo Vivo ^{3}
Physical Review Letters, American Physical Society, 2014, 113, pp.070202
We derive an analytical formula for the covariance $\mathrm{Cov}(A,B)$ of two smooth linear statistics $A=\sum_i a(\lambda_i)$ and $B=\sum_i b(\lambda_i)$ to leading order for $N\to\infty$, where $\{\lambda_i\}$ are the $N$ real eigenvalues of a general onecut randommatrix model with Dyson index $\beta$. The formula, carrying the universal $1/\beta$ prefactor, depends on the randommatrix ensemble only through the edge points $[\lambda_,\lambda_+]$ of the limiting spectral density. For $A=B$, we recover in some special cases the classical variance formulas by Beenakker and DysonMehta, clarifying the respective ranges of applicability. Some choices of $a(x)$ and $b(x)$ lead to a striking \emph{decorrelation} of the corresponding linear statistics. We provide two applications  the joint statistics of conductance and shot noise in ideal chaotic cavities, and some new fluctuation relations for traces of powers of random matrices.
 1. Dipartimento di Matematica
 2. INFN, Sezione di Bari  Istituto Nazionale di Fisica Nucleare, Sezione di Bari
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Universal Order and Gap Statistics of Critical Branching Brownian Motion
Kabir Ramola ^{1} Satya N. Majumdar ^{1} Gregory Schehr ^{1}
Physical Review Letters, American Physical Society, 2014, 112, pp.210602
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,tn)$ 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 \inftyn) = p(g_kn)$ with an algebraic tail $p(g_kn) \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.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Universal properties of branching random walks in confined geometries
Clélia De Mulatier ^{1, 2} Alain Mazzolo ^{2} Andrea Zoia ^{2}
Europhysics Letters, EDP Science, 2014, 107, pp.30001
Characterizing the occupation statistics of a radiation flow through confined geometries is key to such technological issues as nuclear reactor design and medical diagnosis. This amounts to assessing the distribution of the travelled length $\ell$ and the number of collisions $n$ performed by the underlying stochastic transport process, for which remarkably simple Cauchylike formulas were established in the case of branching Pearson random walks with exponentially distributed jumps. In this Letter, we show that such formulas strikingly carry over to the much broader class of branching processes with arbitrary jumps, provided that scattering is isotropic and the average jump size is finite.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. DM2S  Département de Modélisation des Systèmes et Structures

Universal statistics of longest lasting records of random walks and Lévy flights
Claude Godreche ^{1} Satya N. Majumdar ^{2} Gregory Schehr ^{2}
t14/073. 23 pages, 4 figures, Typos corrected. 2014
We study the record statistics of random walks after $n$ steps, $x_0, x_1,\ldots, x_n$, with arbitrary symmetric and continuous distribution $p(\eta)$ of the jumps $\eta_i = x_i  x_{i1}$. We consider the age of the records, i.e. the time up to which a record survives. Depending on how the age of the current last record is defined, we propose three distinct sequences of ages (indexed by $\alpha$ = I, II, III) associated to a given sequence of records. We then focus on the longest lasting record, which is the longest element among this sequence of ages. To characterize the statistics of these longest lasting records, we compute: (i) the probability that the record of the longest age is broken at step $n$, denoted by $Q^{\alpha}(n)$, which we call the probability of record breaking and: (ii) the duration of the longest lasting record, $\ell_{\max}^{\alpha}(n)$. We show that both $Q^{\alpha}(n)$ and the full statistics of $\ell_{\max}^{\alpha}(n)$ are universal, i.e. independent of the jump distribution $p(\eta)$. We compute exactly the large $n$ asymptotic behaviors of $Q^{\alpha}(n)$ as well as $\langle \ell_{\max}^{\alpha}(n)\rangle$ (when it exists) and show that each case gives rise to a different universal constant associated to random walks (including Lévy flights). While two of them appeared before in the excursion theory of Brownian motion, for which we provide here a simpler derivation, the third case gives rise to a nontrivial new constant $C^{\rm III} = 0.241749 \ldots$ associated to the records of random walks. Other observables characterizing the ages of the records, exhibiting an interesting universal behavior, are also discussed.
 1. Institut de Physique Théorique CEA Saclay
 2. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Velocity filtration and temperature inversion in a system with longrange interactions
Lapo Casetti ^{1, 2} Shamik Gupta ^{3}
6 pages, 9 figures, epj macros. 2014
Temperature inversion due to velocity filtration, a mechanism originally proposed to explain the heating of the solar corona, is demonstrated to occur also in a simple paradigmatic model with longrange interactions, the Hamiltonian meanfield model. Using molecular dynamics simulations, we show that when the system settles into an inhomogeneous quasistationary state in which the velocity distribution has suprathermal tails, the temperature and density profiles are anticorrelated: denser parts of the system are colder than dilute ones. We argue that this may be a generic property of longrange interacting systems.
 1. INAFOsservatorio Astrofisico di Arcetri
 2. Dipartimento di Fisica e Astronomia and CSDC
 3. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Waiting time distribution for trains of quantized electron pulses
Mathias Albert ^{1}
PHYS. REV. B, 2014, 90, pp.10
We consider a sequence of quantized Lorentzian pulses of noninteracting electrons impinging on a quantum point contact and study the waiting time distribution (WTD), for any transmission and any number of pulses. As the degree of overlap between the electronic wave functions is tuned, the WTD reveals how the correlations between particles are modified. In the weak overlap regime, the WTD is made of several equidistant peaks, separated by the same period as the incoming pulses, contained in an almost exponentially decaying envelope. In the other limit, the WTD of a single quantum channel subjected to a constant voltage is recovered. In both cases, the WTD stresses the difference between the fluctuations induced by the scatterer and the ones encoded in the incoming quantum state. A clear crossover between these two situations is studied with numerical and analytical calculations based on scattering theory.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques

Winding statistics of a Brownian particle on a ring
Anupam Kundu ^{1} Alain Comtet ^{1, 2} Satya N. Majumdar ^{1}
Journal of Physics A: Mathematical and Theoretical, Institute of Physics: Hybrid Open Access, 2014, 47 (38), pp.385001
We consider a Brownian particle moving on a ring. We study the probability distributions of the total number of turns and the net number of counterclockwise turns the particle makes till time t. Using a method based on the renewal properties of Brownian walker, we find exact analytical expressions of these distributions. This method serves as an alternative to the standard path integral techniques which are not always easily adaptable for certain observables. For large t, we show that these distributions have Gaussian scaling forms. We also compute large deviation functions associated to these distributions characterizing atypically large fluctuations. We provide numerical simulations in support of our analytical results.
 1. LPTMS  Laboratoire de Physique Théorique et Modèles Statistiques
 2. IHP  Institut Henri Poincaré