Statistical mechanics of coupled supercooled liquids in finite dimensions
Benjamin Guiselin, Ludovic Berthier, Gilles Tarjus
Abstract
We study the statistical mechanics of supercooled liquids when the system evolves at a temperature T <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>T</mml:mi> </mml:math> with a field \epsilon <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>ϵ</mml:mi> </mml:math> linearly coupled to its overlap with a reference configuration of the same liquid sampled at a temperature T_0 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:msub> <mml:mi>T</mml:mi> <mml:mn>0</mml:mn> </mml:msub> </mml:math> . We use mean-field theory to fully characterize the influence of the reference temperature T_0 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:msub> <mml:mi>T</mml:mi> <mml:mn>0</mml:mn> </mml:msub> </mml:math> , and we mainly study the case of a fixed, low- T_0 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:msub> <mml:mi>T</mml:mi> <mml:mn>0</mml:mn> </mml:msub> </mml:math> value in computer simulations. We numerically investigate the extended phase diagram in the (\epsilon,T) <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mo stretchy="false" form="prefix">(</mml:mo> <mml:mi>ϵ</mml:mi> <mml:mo>,</mml:mo> <mml:mi>T</mml:mi> <mml:mo stretchy="false" form="postfix">)</mml:mo> </mml:mrow> </mml:math> plane of model glass-forming liquids in spatial dimensions d=2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>d</mml:mi> <mml:mo>=</mml:mo> <mml:mn>2</mml:mn> </mml:mrow> </mml:math> and d=3 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>d</mml:mi> <mml:mo>=</mml:mo> <mml:mn>3</mml:mn> </mml:mrow> </mml:math> , relying on umbrella sampling and reweighting techniques. For both 2d <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mn>2</mml:mn> <mml:mi>d</mml:mi> </mml:mrow> </mml:math> and 3d <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mn>3</mml:mn> <mml:mi>d</mml:mi> </mml:mrow> </mml:math> cases, a similar phenomenology with nontrivial thermodynamic fluctuations of the overlap is observed at low temperatures, but a detailed finite-size analysis reveals qualitatively distinct behaviors. We establish the existence of a first-order transition line for nonzero \epsilon <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mi>ϵ</mml:mi> </mml:math> ending in a critical point in the universality class of the random-field Ising model (RFIM) in d=3 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>d</mml:mi> <mml:mo>=</mml:mo> <mml:mn>3</mml:mn> </mml:mrow> </mml:math> . In d=2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mi>d</mml:mi> <mml:mo>=</mml:mo> <mml:mn>2</mml:mn> </mml:mrow> </mml:math> instead, no phase transition is found in large enough systems at least down to temperatures below the extrapolated calorimetric glass transition temperature T_g <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:msub> <mml:mi>T</mml:mi> <mml:mi>g</mml:mi> </mml:msub> </mml:math> . Our results confirm that glass-forming liquid samples of limited size display the thermodynamic fluctuations expected for finite systems undergoing a random first-order transition. They also support the relevance of the physics of the RFIM for supercooled liquids, which may then explain the qualitative difference between 2d <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"> <mml:mrow> <mml:mn>2</mml:mn> <mml:mi>d</mml:mi> </mml:mrow> </mml:math> and 3d <mml:math xmlns:mml="http://www.