format: text

text:
# Statistics of lowest droplets in two-dimensional Gaussian Ising spin glasses

## I Introduction.

Despite three decades of work in the field of spin glasses major issues related to their low-temperature behavior still remain unresolved . Although important achievements have been obtained in the understanding of mean-field theory  the appropriate treatment beyond mean-field to include short-range interactions is yet to be found. The absence of a successful analytical approach to deal with this problem corroborates the present state of our knowledge, often misguided by a non-accurate, if not confusing, interpretation of the numerical data. This situation has generated a hot debate about the correct physical interpretation of the available numerical data. Leaving aside the long-standing controversy whether replica symmetry breaking is or not a good description of the spin-glass phase , there are still unresolved issues which are not as striking but evidence our ignorance about some fundamental questions.

One among these problems is the correct value of the thermal exponent in two-dimensional (2d) Gaussian Ising spin glasses (GISG). This question has received attention from time to time during the last two decades but not enough to settle it definitively and explain the origin of some of the reported discrepancies. The study of the low-$`T`$ properties of the 2d GISG starts with the work by McMillan who proposed  that thermal properties in spin glasses are determined by the scaling behavior of the typical largest excitations (commonly referred as droplets) present in the system. This idea has been further elaborated and extended to deal with equilibrium and dynamical properties of spin glasses in a scenario nowadays referred to as droplet model . The low-$`T`$ behavior in spin glasses is determined by a spectrum of large scale gapless droplets with typical length $`L`$ and energy cost $`EL^\theta `$, $`\theta `$ being the thermal exponent. As these droplets correspond to flipping some domains of spins (assumed to be compact clusters), the energy cost of these excitations arises from the set of unsatisfied bonds on their surface. The striking low-$`T`$ behavior in spin glasses arises from multiple energy cancellations occurring at the surface of the droplet. These cancellations can be seen as the result of a competition between energy and entropy effects: as the droplet becomes progressively larger there are more available conformations for the surface to minimize the energy cost of the unsatisfied bonds. In the absence of cancellations one would expect $`\theta =(d1)/2`$. However, as these cancellations are very important, the inequality $`\theta <(d1)/2`$ holds and $`\theta `$ is by far less than the maximum value $`(d1)/2`$. The value of the thermal exponent $`\theta `$ characterizes the low-$`T`$ critical behavior as it is related to the correlation length exponent $`\nu `$ where $`\xi T^\nu `$ by the identity $`\nu =1/\theta `$. McMillan also used domain-wall renormalization group ideas to introduce a practical way to determine the leading energy cost of these low-lying large-scale excitations . The method consists in measuring the energy defect of a domain-wall spanning the whole system obtained by computing the change of the ground state energy when switching from periodic to anti periodic boundary conditions in one direction. Several works have used McMillan’s method to determine the value of $`\theta `$ in two and three dimensions . Hereafter, in order to keep the discussion as clear as possible, we will denote by $`\theta _{DW}`$ the estimate of the exponent $`\theta `$ obtained by domain-wall calculations. The initial value for $`\theta _{DW}`$ reported by McMillan is $`\theta _{DW}=0.281(5)`$ for pretty modest lattice sizes $`L=38`$. Recent numerical results with much more powerful algorithms have reached sizes $`L500`$ and confirmed the initial result with much larger accuracy  $`\theta _{DW}=0.287(4)`$. These studies would definitively close the problem if it were not by the existence of other alternative estimates of the exponent $`\theta `$, largely consistent among them, which yield a quite different value $`\theta 0.47(2)`$. We will denote this estimate by $`\theta _{TF}`$ as several of these methods use transfer matrix . However, a word of caution is necessary here as the Monte Carlo method and other approaches that are not based on transfer matrix methods report values compatible with that estimate. For instance, Kawashima and Aoki used another method to estimate the stiffness exponent . The idea is to generate a droplet inside a box of size $`L\times L`$ that includes a fixed central spin, with the following procedure. First, the ground state is found with a standard algorithms (we will denote it by the reference configuration). Afterwards, the spins at the boundaries of the box are fixed and the central spin is forced to flip respect to the reference configuration. The droplet of minimum energy that includes the central reversed spin and does not touch the boundaries is computed. The spanning length of the droplets generated in this way allows to define the fractal dimension of both the surface (or perimeter for the two dimensional case) and the volume. It is found that these minimum energy droplets have a fractal volume dimension smaller than 2 and the thermal exponent is $`\theta =0.42(5)`$ in agreement with results obtained from MC methods  and heuristic optimization algorithms . A similar study of minimum energy clusters in the three-dimensional Edwards-Anderson model also reports evidence that $`\theta _{DW}`$ is an upper bound to the actual value of the thermal exponent .

The accuracy of previous estimates is poorer than the values obtained through the domain-wall method as they deal, in one way or another, with all possible excitations and not only with the calculation of ground state energies. More recently, another method has been used to estimate the value of $`\theta `$. It consists in perturbing the original Hamiltonian $`_0`$ with a term $`ϵ𝒫`$, where $`𝒫`$ stands for the perturbation and $`ϵ`$ for its intensity. For example, $`P`$ can be the overlap between the actual configuration and the ground state of the original Hamiltonian $`_0`$. As $`ϵ`$ varies the new ground state of the total Hamiltonian $`=_0+ϵ𝒫`$ remains unchanged until a certain value $`ϵ=ϵ_c`$ is reached where a excited energy level of $`_0`$ becomes the new ground state of $``$. The overlap between the old and the new ground states as well as the value of the shifting energy provoked by the perturbation links its energy cost $`E`$ with its size providing another way to estimate $`\theta `$. We will denote by $`\theta _P`$ the estimate obtained in this way. This method has been recently used in the 2d GISG by Hartmann and Young  reporting the following value $`\theta _P0.31`$. Although slightly more negative than $`\theta _{DW}`$, both $`\theta _{DW},\theta _P`$ appear to be statistically compatible. Yet more accurate estimates are needed to confirm whether $`\theta _P=\theta _{DW}`$.

This last method and the domain-wall method have in common the same feature, i.e. they perturb the original Hamiltonian in one way or another to probe the characteristic energy of excitations that are supposed to be the typical ones that determine the low-$`T`$ thermodynamic properties. In fact, the estimate $`\theta _{DW}`$ can be considered as a particular example of $`\theta _P`$, where the perturbation consists in reversing all the bonds in one of the surfaces of the box. This raises the important question whether the different estimates of $`\theta _P`$, obtained by considering different class of perturbations, are different. The question is rather subtle as there are numerical indications that indeed this could be the case. For instance , measurements of $`\theta _P`$ where the perturbation is a uniform magnetic field yield a value $`\theta _P=0.48(1)`$ compatible with the other competing set of values $`\theta _{TF}`$.

How is that the value of the exponent $`\theta _P`$ could depend on the type of perturbation? This is a very difficult question to answer as our present knowledge is inadequate. We can offer only speculative answers. Strong discrepancies among different types of perturbations could arise if a multi-fractal scenario governs the statistics of excitations in spin glasses. By definition, in all perturbation methods the probed large scale droplets are those which minimize the energy cost but constrained to maximize the value of the perturbation for the selected droplets. Therefore, among all possible large-scale low-lying droplets the perturbation method selectively probes those that maximally overlap with the perturbation. A dependence of the value $`\theta _P`$ on a given class of perturbations could arise if the perturbation selectively probes one or another topological property of the droplet. This rather awkward multi-fractal scenario is not new in the field of disordered systems. Multifractality is known to be present in the localization problem in the strongly disordered regime. A multi-fractal scenario would imply the existence of different critical exponents at $`T=0`$ depending on the type of perturbation applied. On the other hand, the fact that the value estimated for $`\theta _P`$ when the perturbation is a uniform magnetic field appears to be consistent with the value $`\theta _{TF}`$, suggests that maybe some types of perturbation can probe the relevant excitations while others may not. These good observables, which probe the typical excitations, could be called neutral observables in the same spirit as this term has been coined to describe observable dependences of the fluctuation-dissipation ratio (i.e. the effective temperature) in glassy systems. Concomitant, this “perturbation class dependence” issue is presently also debated in the different (but related to a certain degree) field of glassy dynamics.

If the hypothetical multi-fractal scenario holds, then we must face the question about what is the correct procedure to determine the thermal exponent $`\theta `$. As $`\theta `$ determines the free energy cost of droplets, the natural answer is that $`\theta `$ is given by the lowest value among all possible estimates,

$$\theta =\mathrm{min}_𝒫\{\theta _𝒫\}.$$
(1)
With the present available data this relation suggests that the estimate $`\theta _{TF}`$ is the correct value of the thermal exponent and that $`\theta _{DW}`$ as well as many other estimates $`\theta _P`$ are only upper bounds to the true value.

The question we want to address in this paper is the following. Is it possible to devise a method that is alternative to current perturbation methods, in which excitations are not selectively probed by the perturbation, but selected only according to the correct balance between energy and entropy? The main purpose of this paper is to show that the analysis of the statistics of the first or lowest excitations gives a positive answer to this question. As we will see, the method we propose in this paper yields a consistent estimate of $`\theta `$ compatible with the value $`\theta _{TF}`$, therefore supports the result that $`\theta _{DW}`$ and many other $`\theta _P`$ are only upper bounds to the actual value of $`\theta `$. A preliminary account of these results has already appeared in .

The paper is divided as follows. Sec. II describes the basis of the lowest droplet approach and introduces the lowest droplet exponents. Sec. III shows the results obtained in the 2d GISG. Sec. IV analyzes a method to extract the value of the thermal exponent $`\theta `$. Sec. V presents a more powerful method to extract the value of the lowest droplet exponents based on an aspect-ratio analysis. Sec. VI discusses the origin of the finite-volume corrections to the value of the lowest droplet exponent $`\theta _l`$ as a problem of corrections in the statistics of extreme values. Sec. VII analyzes some topological properties of the lowest droplets. Finally Sec. VIII presents the conclusions. There are also two technical appendixes: Appendix A presents the heuristic argument that $`\theta _l=d`$ for Gaussian spin glasses, and Appendix B explains the transfer matrix method we used to obtain the lowest droplets.

## II Basis of the lowest droplet approach.

The purpose of this work is to show an alternative approach to determine the low-$`T`$ behavior of spin glasses by studying the size and energy spectrum of the lowest excitations by introducing two exponents ($`\lambda _l`$ and $`\theta _l`$) needed to fully characterize the zero-temperature fixed point. All along the paper we will denote this exponents as lowest energy droplet exponents, or lowest droplet exponents in short, and that we will abbreviate as LDEs. The exponent $`\lambda _l`$ is the most important one and describes the probability to find a large-scale lowest excitation spanning the whole system, while the exponent $`\theta _l`$ describes the system-size dependence of the average energy cost of these lowest excitations.

The underlying theoretical background of the approach is the following. To investigate the leading low-temperature behavior in spin glasses let us consider expectation values for moments of the order parameter by keeping only the ground state and the first or lowest excitation. This approach was introduced in  and can be shown to capture the low-temperature behavior at the leading order. The method that investigates the low-$`T`$ properties based on a restricted analysis of the spectrum to the absolute lowest excitations has been also used for the study of the localized phase in the disordered Anderson model . The present paper can be seen as the applicability of these ideas to the spin-glass case. At the end of the paper (see Sec. VIII) we will give reasons supporting the validity of our approach.

To generate the spectrum of lowest excitations we consider the following procedure. Let us consider a set of $`𝒩_s`$ samples and for each of them we determine both the configurations of the ground state and the lowest excitation. For a spin model the lowest excitation has $`v`$ spins overturned with respect to the ground state (so the overlap between the ground and that excited state is $`q=12v/V`$, $`V`$ being the volume of the system) and with energy cost or gap $`E`$. It can be easily proved that the lowest excitation must be a connected cluster which we will generically call the lowest droplet. If $`v_s`$ and $`E(s)`$ denote the volume and excitation energy of the lowest droplet for sample $`s`$, in the limit where $`𝒩_s`$ is sent to infinity, we can define the following joint probability distribution

$$P(v,E)=\frac{1}{𝒩_s}\underset{s=1}{\overset{𝒩_s}{}}\delta (vv_s)\delta (EE(s)).$$
(2)
Using the Bayes theorem, this joint probability distribution can be written as $`P(v,E)=g_v\widehat{P}_v(E)`$ where

$$\underset{v=1}{\overset{\frac{V}{2}}{}}g_v=1;_0^{\mathrm{}}dE\widehat{P}_v(E)=1v.$$
(3)
$`g_v`$ is the probability to find a sample such that its lowest droplet has volume $`v`$ and $`\widehat{P}_v(E)`$ is the conditioned probability for that droplet to have a gap equal to $`E`$. In what follows, we separately discuss the scaling behavior of both distributions $`g_v,\widehat{P}_v(E)`$.

Before continuing, and for sake of clarity, let us make an important digression about nomenclature. There are two volumes involved in the problem: the volume $`v`$ of the lowest excitation and the volume $`V`$ of the lattice. If not stated otherwise we will refer to the volume $`v`$ as the size of the excitation while volume will generally refer to the lattice volume $`V`$. Thus, when we speak about finite-size excitations we usually refer to excitations with $`v`$ finite, and finite-volume corrections (which we will sometimes abbreviate as FVC) will refer to the corrections affecting the distribution (2) due to the finite volume $`V`$ of the lattice.

### A The lowest droplet exponent $`\lambda _l`$.

The simplest scenario for the size distribution of the lowest droplets is that all sizes occur with uniform probability. The normalization condition (3) imposes $`g_v1/V`$. This situation is encountered in the 1d GISG  with both free and periodic boundary conditions. However, in the most general situation, this does not hold and low energy droplets are found with a probability that depends on their size $`v`$. The simplest and most general way to incorporate such a dependence is to assume an ansatz solution for $`g_v`$ that factorizes into a power law $`A/V^{\lambda _l+1}`$ with $`\lambda _l>0`$ and a coefficient $`AG(q)`$ which depends only on the overlap $`q`$ between the ground state and the lowest droplet,

$$g_v=\frac{G(q)}{V^{\lambda _l+1}}.$$
(4)
The behavior of $`G(q)`$ can be guessed in both limits $`q1`$ (the case $`q1`$ is equivalent in models with time-reversal symmetry which are those we are considering here) and $`q0`$,

$`G(q0)\mathrm{constant}`$ (5)
$`G(q1){\displaystyle \frac{1}{(1q)^{\lambda _l+1}}}.`$ (6)
The first relation describes the scaling behavior for the number of droplets whose size scales with the total volume of the system. As these can only depend on the volume $`V`$, $`G(0)`$ must converge to a constant. The second relation is consequence of the fact that the number of droplets with finite size $`v`$ cannot depend on $`V`$ in the large $`V`$ limit as these are not affected by the boundaries. On the other hand, the distribution of finite size droplets $`g_v`$ is self-similar as can be seen by inserting (6) in (4) and using the relation $`q=12v/V`$. This yields $`g_v1/v^{\lambda _l+1}`$, the same relation as for the large scale limit (5) where $`g_V1/V^{\lambda _l+1}`$. A simple expression that interpolates both limits is given by,

$$G(q)=\left(A+\frac{B}{(1q)^{\lambda _l+1}}\right).$$
(7)
Note however that, despite its simplicity, expression (7) is only an interpolation and the most we can say about $`G(q)`$ concerns its asymptotic behaviors (5-6).

The ansatz (4), applied only to large-scale excitations, was proposed in . Note that although $`g_v`$ is defined for discrete volumes, in the limit $`V1`$, the values of $`q`$ for consecutive droplet sizes $`vv+1`$ become equally spaced by $`\mathrm{\Delta }q=2/V`$. Therefore, in the limit, $`V1`$, the function $`g(q)=\frac{V}{2}g_v`$ becomes a continuous function if expressed in terms of the variable $`q`$ instead of the integer variable $`v`$,

$$g(q)=\frac{1}{2V^{\lambda _l}}G(q).$$
(8)
A word of caution is in order. Although (4) diverges for $`q=1`$, leading apparently to a violation of the normalization condition (3) for $`g_v`$, it must be emphasized that no excitation has $`q=1`$ so there is a maximum cutoff value $`q^{}=12/V`$ corresponding to one-spin excitations. For instance, if we insert (8) into the normalization condition for $`g(q)`$ we get in the large $`V`$ limit,

$$_0^{q^{}=12/V}g(q)𝑑q=1\frac{AB/\lambda _l}{2V^{\lambda _l}}+\frac{B}{2^{\lambda _l+1}\lambda _l}=1,$$
(9)
implying $`\lambda _l0`$ as expected since otherwise the normalization would not be possible in the large-$`V`$ limit. The divergent term ($`q1`$) in (8) shows that for $`\lambda _l>0`$ one-spin excitations are the most large in number among the whole spectrum of sizes. In fact, from (4), $`g(1)O(1)g(V/2)1/V^{\lambda +1}`$, so the majority of excitations have a finite size. However the average excitation size

$$\overline{v}=\underset{v=1}{\overset{V}{}}vg_v_V\mathrm{}V^{1\lambda _l},$$
(10)
diverges in the $`V\mathrm{}`$ limit and differs from the typical excitation volume $`v_{\mathrm{typ}}O(1)`$. Relation (10) provides a way to measure the exponent $`\lambda _l`$ alternative to the use of the scaling behavior (4).

### B The lowest droplet exponent $`\theta _l`$

The analysis of the gap distribution $`\widehat{P}_v(E)`$ goes along the same lines as we did for the distribution $`g_v`$, but with one important difference. As the gap $`E`$ describes the lowest among all possible excitation energies, it has to scale in the same way for all droplet sizes independently on their size (and, in particular, whether these are finite-size or large scale droplets). This statement refers to a scenario which hereafter we will call “the random energy-size droplet scenario” (RESD scenario) to specifically indicate that the distribution of the lowest energies of droplets is independent of their size. Mathematically it can be expressed as,

$$\widehat{P}_v(E)=\widehat{P}(E),v.$$
(11)
In addition, we follow the standard droplet model and assume that the spectrum is gapless and defined by an exponent $`\theta _l`$ which describes the characteristic energy of the lowest droplets whatever their size or overlap $`q`$ with the ground state. If the scaling function $`\widehat{P}_v(E)`$ is independent of $`v`$ it follows immediately that the non-conditioned or size-averaged gap probability distribution

$$P(E)=\underset{v1}{}g_v\widehat{P}_v(E)=\underset{v1}{}g_v\widehat{P}(E)=\widehat{P}(E)$$
(12)
where we used (11) and the normalization condition (3) for $`g_v`$. From now on, if not stated otherwise, we will always refer to the size-averaged probability distribution $`P(E)`$ with the clear understanding that it coincides with any of the conditioned distributions $`\widehat{P}_v(E)`$. As the spectrum of lowest excitations is gapless, the normalized distribution $`P(E)`$ has the following scaling behavior,

$$P(E)=\frac{1}{L^{\theta _l}}𝒫\left(\frac{E}{L^{\theta _l}}\right).$$
(13)
We stress that the exponent $`\theta _l`$ is completely different from the standard thermal exponent (see next section) as they describe totally different excitations. The thermal exponent $`\theta `$ describes the energy-length relation for droplets typically excited at finite temperatures while the lowest energy exponent $`\theta _l`$ describes the droplets that are separated by the smallest gap respectively to the ground state, so that, in general, $`\theta _l\theta `$.

We will argue below in Sec. II C that $`\theta _l=d`$ for a generic class of spin-glass systems with coupling distributions with finite weight at zero gap. In addition, this relation will provide an alternative interpretation of the lower critical dimension in terms of the exponent $`\lambda _l`$ introduced in Sec. II A describing the properties of the spectrum of sizes of the lowest droplets.

### C The standard thermal exponent $`\theta `$

Now we want to show how the exponents $`\lambda _l`$ and $`\theta _l`$ combine to give the usual scaling exponent $`\theta `$ describing the energy cost of typical thermal excitations in droplet theory. There are several ways to show this result. For simplicity, here we exemplify this relation by analyzing the low-$`T`$ behavior of the second moment of the spin-glass order parameter at the order linear in $`T`$ by keeping only the first excitation. If $`q_{\{\sigma ,\tau \}}=\frac{1}{V}_i\sigma _i\tau _i`$ denotes the overlap between two replicas (i.e. configurations of different systems with the same realization of quenched disorder), then the expectation value $`\overline{q^2}`$ can be written as follows ,

$$\overline{<q^2>}=1\frac{2}{V^2}\underset{v}{}_0^{\mathrm{}}𝑑EP(v,E)v(Vv)\mathrm{sech}^2\left(\frac{E}{2T}\right),$$
(14)
where $`P(v,E)`$ is given by (2). A low-temperature expansion of (14)  up to linear order in $`T`$ yields,

$$\overline{<q^2>}=1\frac{4T}{V^2}\underset{v=1}{\overset{V}{}}g_v\widehat{P}_v(0)v(Vv)$$
(15)
which shows that the leading behavior is determined by both $`g_v`$ and the density of states at zero gap $`\widehat{P}_v(0)`$. In the standard droplet model, it is generally assumed that typical low energy droplets have an average size $`\overline{v}=_vvg_vV`$ of the order of the system size (such as those generated by DW perturbation) and finite weight at zero gap $`\widehat{P}_V(0)1/L^\theta `$ where $`\theta `$ is the thermal exponent. In principle, a single exponent $`\theta `$ describes the scaling behavior of typical large-scale droplets with volume $`vV`$ and determines the zero-temperature critical behavior. As these large-scale droplets are typical they occur with finite (therefore independent of $`V`$) probability $`g_VO(1)`$ while small scale droplets are simply irrelevant $`g_{vO(1)}0`$. This yields,

$$\overline{<q^2>}=1c\frac{T}{L^\theta },$$
(16)
where $`c`$ is a non-universal stiffness constant related to the particular model. One of the most relevant results from the ansatz (4) is that both small and large scale excitations contribute to low-temperature properties. In general, let us consider any expression (such as (15)) involving a sum over all possible volume excitations. Restricting the sum to the large-scale droplets ($`v/V`$ finite) the net contribution to such sum is proportional to $`Vg_V\widehat{P}_V(0)L^{\theta _ld\lambda _l}𝒫(0)`$ (where $`𝒫`$ is the scaling function appearing in (13)). Coming back to (15) and using (4) and (13), we note that both small and large-scale excitations yield a contribution to (15) of the same order and given by,

$$\overline{<q^2>}=1c_l\frac{T}{L^{\theta _l+d\lambda _l}},$$
(17)
where $`c_l`$ is another constant (different from the constant $`c`$ appearing in (16)). Identifying both relations (16) and (17) we obtain the general relation,

$$\theta =\theta _l+d\lambda _l.$$
(18)
This relation shows how the value of $`\theta `$ can be computed from $`\lambda _l`$ and $`\theta _l`$. Through the study of a specific example, we will see later that the exponents $`\theta _l`$ and $`\lambda _l`$ have strong finite-volume corrections arising from the corrections present in the statistics of the extreme values. However, we will present alternative routes to overcome this dependence and provide an accurate estimate of $`\theta `$.

Now we come back to the aforementioned argument at the end of Sec. II B claiming that in the large-volume limit $`\theta _l`$ must converge to the value $`d`$ in the case of coupling distributions with finite gap at zero coupling. The details of the argument are shown in appendix A. The argument has two parts. First, it is proved that one-spin excitations provide an upper bound for the LDE $`\theta _l`$. Then it is argued that this upper bound holds also for any finite-size excitations (such as two-spin clusters, three-spin clusters, and so on). We will see below how this result is supported by the numerical analysis of the data. Let us also note that this result, in a RESD scenario (see Sec. II B) can be linked to the linear dependence of the specific at low temperatures, a result widely accepted, but that has been revisited recently in  to show that it has strong FVC due to the systematic FVC present in the value of $`\theta _l`$. Inserting $`\theta _l=d`$, (18) becomes,

$$\theta =d(\lambda _l1).$$
(19)
This relation provides a way to distinguish the lower critical dimension $`d_{lcd}`$ in terms of the average size distribution of the lowest droplets. According to (10) the relation $`\lambda _l(d_{lcd})=1`$ distinguishes a regime where the average size of the lowest droplet grows with the volume of the system to a regime where the average size of the lowest droplet is finite,

$`d<d_{lcd}:\underset{V\mathrm{}}{lim}\overline{v}(V)=\mathrm{};\lambda _l<1,\theta <0`$ (20)
$`d>d_{lcd}:\underset{V\mathrm{}}{lim}\overline{v}(V)=O(1);\lambda _l>1,\theta >0.`$ (21)
The marginal case $`\lambda _l=1,\theta =0`$ is specially interesting as the average size $`\overline{v}`$ could be finite or diverge with the size but slower than a power law. This scenario corresponds to the mean-field behavior as replica-symmetry is broken in both the standard RSB  or in the TNT  (standing for trivial-non trivial) scenarios. Therefore, the study of the size spectrum of the lowest excitations in spin glasses can be very useful to find out the correct value of the thermal exponent in models without a finite-$`T`$ transition (such as the 2d GISG) as well as establishing the correct low-$`T`$ scenario in models with a finite-$`T`$ transition. In the next section we apply all these ideas to evaluate the thermal exponent for the $`2d`$ GISG.

## III Statistics of the lowest energy droplets in the 2d GISG

Several numerical works have recently searched for low-lying excitations in spin glasses using heuristic algorithms . But, to our knowledge, no study has ever presented exact results about the statistics of lowest excitations. We have exactly computed ground states and lowest excitations in two-dimensional Gaussian spin glasses defined by

$$=\underset{i<j}{}J_{ij}\sigma _i\sigma _j,$$
(22)
where the $`\sigma _i`$ are the spins ($`\pm 1)`$ and the $`J_{ij}`$ are quenched random variables extracted from a Gaussian distribution of zero mean and unit variance. These have been computed by using a transfer matrix method working in the spin basis. Representing each spins state by a weight and a graduation in the energy we can build explicitly the ground state by keeping the largest energy and, by subsequent iteration, the first excitation and so on (see Appendix B for the details on how we compute these quantities). The continuous values for the couplings assures that there is no accidental degeneracy in the system (apart from the trivial time-reversal symmetry $`\sigma \sigma `$). Calculations have been done in systems with free boundary conditions in both directions (FF), periodic boundary conditions in both directions (PP) and free boundary conditions in one direction but periodic in the other (FP). In all cases we find the same qualitative and quantitative results indicating that we are seeing the correct critical behavior.

We have found ground states and lowest droplets for systems ranging from $`L=4`$ up to $`L=11`$ for PP and up to $`L=16`$ for FP and FF. The number of samples is very large, typically $`10^6`$ for all sizes. The large number of samples assures us that many samples have large-scale droplets as first excitations. This provides us with good statistics to properly analyze the sector of large-scale excitations. The large number of samples requires a big amount of computational time so that calculations were done in a PC cluster during several months. For each sample we have evaluated the volume of the excitation $`v`$ (and hence the overlap $`q=12v/V`$ between the ground state and the first excitation) and the gap $`E`$. From these quantities we can construct the $`g_v`$ and the $`\widehat{P}_v(E)`$.

In Fig. 1 we show $`g(q)=\frac{V}{2}g_v`$ as function of $`q`$ for different sizes in the PP and FF cases. We can clearly see that there are excitations of all possible sizes but, as discussed in the paragraph following (9), the typical ones which dominate by far are single spin excitations. To have a rough idea of the number of rare samples giving large scale excitations let us say that nearly half of the total number of samples have one-spin lowest excitations, whereas less than $`10\%`$ of the samples have lowest excitations with overlap $`q`$ in the range $`00.5`$. This disparity increases systematically with size. For the lattice sizes explored the typical number of large-scale droplets is in the range $`10^410^5`$ which is, indeed, quite good to have a good sampling of the sector corresponding to large scale excitations. A detailed analysis of the shape of $`g_v`$ reveals that it has a flat tail for large-scale excitations and a power-law divergence for finite-size excitations. The $`g_v`$ can be excellently fitted by the interpolating formula (7,8)),

$$g(q)=\frac{2}{V^{\lambda _l}}\left(A+\frac{B}{(1q)^{\lambda _l+1}}\right).$$
(23)
As shown in the insets of Fig. 1 a good collapse of the scaling function is obtained with the effective exponent $`\lambda _l^{\mathrm{eff}}0.7`$ for both PP and FF cases. We also plot the line resulting from the fit of (23) with numerical data with the following values for $`A`$ and $`B`$: PPBC: $`A=1.55(3)`$ and $`B=0.777(3)`$; FFBC: $`A=2.02(3)`$ and $`B=0.85(1)`$. Note that the fit is excellent and is hardly distinguishable from the points. The value of $`\lambda _l`$ is compatible with the one obtained by fitting the average size with the expression (10) with the addition of a constant term to account for the small-$`V`$ behavior, $`\overline{v}=C_1+C_2V^{1\lambda }`$. The same exponent $`\lambda _l`$ can be estimated by measuring the ratio $`g(V/2)/g(1)D_1+D_2V^{1\lambda }`$. In both cases we get an effective exponent $`\lambda _l^{\mathrm{eff}}=0.70(5)`$ as best fitting value.

However these different estimates of $`\lambda _l`$ are strongly affected by finite-volume corrections (FVC). To evidence them we have estimated an effective $`L`$ dependent $`\lambda _l^{\mathrm{eff}}(L)`$ exponent by relating the average excitation size at consecutive sizes and using relation (10),

$$\lambda _l^{\mathrm{eff}}(L)=1\frac{1}{d}\frac{\mathrm{log}\left(\frac{\overline{v}(L+1)}{\overline{v}(L)}\right)}{\mathrm{log}\left(\frac{L+1}{L}\right)}.$$
(24)
In Fig. 2 we show $`\lambda _l^{\mathrm{eff}}(L)`$ in the range $`L=411`$ for the PP case. As we can appreciate there is a systematic increase of the effective exponent as we go to large volume sizes without any tendency to saturate. This proves that FVC in our measurements are still big and the estimate $`\lambda _l^{\mathrm{eff}}`$ used to collapse the data in Fig. 1 is still far from the asymptotic exact value.

After having discussed the $`g_v`$ we jump now to discuss the scaling behavior of the energy gap distribution $`\widehat{P}_v(E)`$ and its average $`P(E)`$. In Fig. 3 we show $`P(E)`$ (main figure and inset a) and $`\widehat{P}_v(E)`$ (inset b) for the PP case. Similar results are obtained for the FF and FP cases. Quite remarkably, as was already anticipated in (11), the RESD scenario holds as the distribution $`\widehat{P}_v(E)`$ does not depend on the size $`v`$ of the excitation (see inset b in Fig. 3), hence both large and finite-size excitations are described by the same gap distribution.

In the main figure we can see how the width of distribution $`P(E)`$ progressively shrinks to 0 as $`L`$ increases. Moreover, the $`P(E)`$ has an exponential shape. This is shown in the inset a) of Fig. 3 where we plot $`P(E)`$ in log-normal scale. Nonetheless, a detailed examination of the tails of $`P(E)`$ reveals some deviations from linearity. In Sec. VI we discuss the origin of these deviations. We anticipate, though, that they are consequence of the strong FVC in the range of sizes investigated. In that inset we also verify the scaling ansatz (13) by showing the best data collapse for $`P(E)`$ obtained with an effective exponent $`\theta _l^{\mathrm{eff}}1.7(1)`$. This is very far from the expected value $`\theta _l=2`$ discussed in the preceding Sec. II C and in the Appendix A. A calculation of the moments of $`P(E)`$ (13) for different values of $`L`$ shows that there are also strong sub-dominant corrections to the leading scaling (13) that result in corrections as large as the ones affecting the exponent $`\lambda _l`$.

Again, to manifest the magnitude of FVC in $`\theta _l`$ we have evaluated $`\overline{E}(L)`$, the first moment of $`\widehat{P}(E)`$, obtained by averaging the lowest gap over all possible droplet sizes for different lattice sizes in the range $`L=411`$. We have estimated an effective $`L`$-dependent exponent by means of the following expression,

$$\theta _l^{\mathrm{eff}}(L)=\frac{\mathrm{log}\left(\frac{\overline{E}(L+1)}{\overline{E}(L)}\right)}{\mathrm{log}\left(\frac{L+1}{L}\right)}.$$
(25)
The results are shown in Fig. 4 for the PP case. Again, as for $`\lambda _l^{\mathrm{eff}}`$ (see Fig. 2), we observe that the estimated value for $`\theta _l^{\mathrm{eff}}`$ systematically changes with size showing that, for the sizes we have explored, we are still far from the asymptotic regime.

We can summarize the results of this section saying that both lowest droplet exponents (LDEs) $`\lambda _l`$ and $`\theta _l`$ display strong systematic finite-volume corrections (FVC). In principle, without further elaboration, it is difficult to give an accurate estimate for the thermal exponent $`\theta `$ using (18). An alternative estimate for the exponent $`\theta `$ could be defined from the analysis of the fraction of large-scale excitations with $`q1/2`$, $`f(q1/2)`$ , which is given by

$$f(q1/2)Vg(0)1/V^{\lambda _l1}L^{d(\lambda _l1)}1/L^\theta ,$$
(26)
where we have used $`\theta _l=d`$ (19). Although (26) yields estimates for $`\theta `$, again these are affected by strong finite-volume corrections. In the range of sizes studied in this paper, and using (26) we get $`\theta 0.6`$ quite far from the asymptotic value reported later in Secs. IV and V. How can we go further and estimate $`\theta `$ in a safer way? In the next two sections we shall answer this question.

## IV A good estimate of the lowest droplet exponents

An interesting aspect of the effective $`L`$-dependent exponents shown in Figs. 2 and 4 is that, while their FVC are large, their corrections are of opposite sign. While $`\lambda _l^{\mathrm{eff}}(L)`$ increases with $`L`$, $`\theta _l^{\mathrm{eff}}(L)`$ decreases. As they have to be added to get $`\theta `$ according to the relation (18) their finite-volume corrections cancel out to a certain degree. If we combine the two estimates for the best data collapse given in the previous section ($`\lambda _l^{\mathrm{eff}}=0.70(5)`$, $`\theta _l^{\mathrm{eff}}1.7(1)`$) we obtain $`\theta 0.3(2)`$ which is very close to the DW value in average. However, this estimate is too pessimistic. A better route would be to use the two LDEs estimated from (24),(25) and adding them according to (18)

$$\theta ^{\mathrm{eff}}(L)=\theta _l^{\mathrm{eff}}(L)+d\lambda _l^{\mathrm{eff}}(L).$$
(27)
In Fig. 5 (left panel) we show the value of $`\theta `$ obtained in this way. Note that the value of the thermal exponent $`\theta `$ has negligible FVC but relatively large statistical fluctuations with $`L`$.

A better, albeit related, way to estimate $`\theta `$ is the following. Instead of independently finding out $`\lambda _l`$ and $`\theta _l`$ we look for an estimator which depends on the appropriate combination of the two exponents $`\theta =\theta _l+d\lambda _l`$. The simplest quantity which satisfies this requirements is given by the combination,

$$A(L)=L^d\frac{\overline{E}(L)}{\overline{v}(L)}.$$
(28)
Since $`\overline{E}(L)L^{\theta _l}`$ and $`\overline{v}(L)L^{d(1\lambda _l)}`$, using (18) we obtain $`A(L)L^\theta `$. To estimate the value of $`\theta `$ we follow two different routes: 1) We use (25) by replacing $`\theta _l^{\mathrm{eff}}(L)\theta ^{\mathrm{eff}}(L)`$ and $`\overline{E}(L)A(L)`$. By definition, this procedure gives exactly the estimate (27) shown in the left panel in Fig. 5. 2) A more stable estimate can be obtained from a fit of $`A(L)`$ versus $`L`$, with data in the range $`[L,\mathrm{},L_{max}=11]`$ (for the PPBC case). This is shown in the left panel of Fig. 5 together with the previous estimate (27) and also in the right panel of Fig. 5 but there compared with the effective exponent $`\theta _{DW}`$ obtained from domain-wall calculations. Our best value for $`\theta `$ is

$$\theta =0.46(1).$$
(29)
This value is very close to the finite-temperature (Monte Carlo or transfer matrix) estimates $`\theta _{TF}=0.48(1)`$  but certainly smaller than the domain-wall value $`\theta _{DW}=0.285`$ . Our estimate for $`\theta `$ is compatible with the other possible value $`\theta _{TF}`$ obtained by other methods as discussed in Sec. I but is certainly inconsistent with the value obtained with other methods with results closer to the DW estimate.

All these estimates strongly support the inequality $`\theta =\theta _{TF}<\theta _{DW}`$. However, one cannot exclude a situation where the present tendency of the data gets modified and $`\theta \theta _{DW}`$ in the large-$`L`$ limit . We have already explained in Sec. II C that $`\theta _l`$ must converge to $`2`$ in the large volume limit implying the relation (19). Introducing our estimate (29) in (19) we get,

$$\lambda _l=0.770(5).$$
(30)
A convincing proof of the correctness of the values (29,30) requires proving that the estimate (24) converges to the value (30) when $`L\mathrm{}`$. In the next section we present an aspect-ratio analysis to evidence that the estimates (29),(30) are correct in the large $`L`$ limit.

## V Aspect-ratio analysis of the lowest droplet exponents

In this section, we present some additional data obtained via an aspect-ratio analysis (ARA). This analysis has been proved to be very useful to extract the value of the domain-wall exponent $`\theta _{DW}`$ by generating domain walls in rectangular lattices $`M\times L`$ with different aspect ratios $`M/L`$ . It has been found that, in the limit of large aspect ratio, the value of $`\theta _{DW}`$ for Gaussian spin glasses is largely independent of the boundary conditions. We have seen in Sec. III that our measurements on squared lattices of size $`L\times L`$ mix small excitations with large ones so one does not have a clear-cut separation in the statistical distribution between the two different regimes $`v𝒪(1)`$ and $`v/V𝒪(1)`$. Our main motivation here is to show that, by investigating large aspect-ratios, we can separate these two different scaling regimes. We made our measurements on systems of size $`L\times M`$, with $`M=LR>>L`$ where $`R`$ ranges from $`1`$ up to $`10`$. We have investigated different types of boundary conditions: periodic boundary conditions in both directions (PPBC) and periodic boundary conditions in the $`L`$ direction with free boundary conditions in the $`M`$ direction (FPBC).

In Fig. (6), we display the data for $`g(q)`$ versus $`1q`$ (8) for the FPBC case for $`L=8`$ and $`R=1,5`$ and $`10`$. One can clearly see that the behavior of the distribution $`g(q)`$ drastically changes as one increases $`R`$. Indeed, as we have already seen in Sec. III and in (5), (6), (8), for $`R=1`$ it is very difficult to separate the region of small excitations (a scaling region with $`g(q)\frac{1}{(1q)^{\lambda _l+1}}`$) from the one of large excitations (a constant $`q`$-independent contribution $`g(q)\frac{1}{V^{\lambda _l}}`$). The main advantage of separating these two regions is that one can fit directly each of them. This yields two separate measurements of the LDE $`\lambda _l`$ in addition to the estimate (10) obtained from the $`L`$ dependence of the average size of the excitations.

In Fig. (7), we show $`g(q)`$ versus $`1q`$ for $`R=10`$ for various linear sizes $`L`$ and for the PPBC case. These distributions have been obtained by running a large number of samples ranging from 10 million of samples for $`L=4`$ down to 5 million for the largest size $`L=9`$. We have also inserted in the figure two vertical lines which indicate the limits for the range of values we have chosen for the fits of the scaling behavior of the finite-size excitation sector ($`1q0.07`$) and for the constant contribution corresponding to large scale excitations ($`1q0.25`$). We have chosen these values for the following reasons. First, as one can clearly see in Figs. (6,7), the scaling region for small excitations survives up to excitation sizes $`vL\times L`$. This size provides a threshold value for the overlap $`q_{\mathrm{th}}`$ below which the simple scaling $`g(q)\frac{1}{(1q)^{\lambda _l+1}}`$ does not hold anymore,

$$1q_{\mathrm{th}}=1(1\frac{2v}{V})\frac{2L^2}{RL^2}\frac{2}{R}.$$
(31)
An second, there is a crossover region around $`qq_{\mathrm{th}}`$. A careful look at Fig. (7) shows that the scaling region for small excitations ends around $`1q0.07`$. At this value, one observes a change of the slope of the curves just before entering the regime of large excitations where $`g(q)`$ becomes $`q`$ independent. For $`1q0.25`$, the curves are rather constant and the result of a fit does not depend much on the choice $`1q=0.25`$. This second threshold value is indicated as the rightmost vertical bar in Fig. (7).

In Fig. (8), we show the estimated values of effective lowest droplet exponent $`\lambda _l^{\mathrm{eff}}`$ obtained in three different ways. The first estimate has been obtained by averaging the volume of all excitations for different lattice sizes as explained in Sec. III and then taking a logarithmic derivative, see (24). The second estimate has been obtained by considering the large excitation sector ($`1q0.25`$) and its $`L,R`$ dependence:

$$g(q)(RL^2)^{\lambda _l}.$$
(32)
Averaging the excitation volume within this sector ($`1q0.25`$) and using again the corresponding logarithmic derivatives as in (24) yields the second estimate. The third estimate for $`\lambda _l`$ is obtained from a direct fit of $`g(q)`$ for small values of $`1q`$ :

$$g(q)(1q)^{1\lambda _l}.$$
(33)
This third method is in fact the most direct one since it can be done for each size $`L`$ (while the other two estimates require a fit using data from two different lattice sizes $`L`$ and $`L^{}`$). The first conclusion that we learn from Fig. (8) is that the ARA produces a great improvement on the estimated values of the exponent $`\lambda _l`$. The most stable measurement is the third estimate obtained by fitting the small-size spectrum of the excitations. In that case, $`\lambda _l^{\mathrm{eff}}`$ is nearly constant with a value that converges to

$$\lambda _l=0.77(1)$$
(34)
in excellent agreement with the result (30) of the previous section.

Moreover, one also observes in Fig. (8) that the two other estimated values for $`\lambda _l^{\mathrm{eff}}`$, obtained with the first and second methods, are strongly correlated. This shows that finite-volume corrections, which are expected to affect the value of the exponent obtained from the analysis of large-size excitations, does affect also the value of the exponent obtained by averaging over the whole spectrum. In addition, we also observe that the ARA for large $`R`$ strongly decreases the magnitude of finite-volume corrections. While on a square geometry, the effective exponent $`\lambda _l^{\mathrm{eff}}`$ obtained from the average size of excitations took values in the range $`0.520.62`$ (see Fig.(2)), with the ARA, we obtain for the same exponent values in the range $`0.640.72`$, which are much closer to the expected asymptotic value $`0.77(1)`$.

The same conclusion holds for the lowest droplet exponent $`\theta _l`$. In Fig. (9), we show the effective exponent $`\theta _l^{\mathrm{eff}}`$ obtained by evaluating the logarithmic derivative as in (25). Note that finite-volume corrections are much smaller than with the squared lattices and as a result, the value of the effective exponent converges much faster to the expected value $`2`$. Using a fit of the form $`\theta _l^{\mathrm{eff}}(L)=\theta _l^{\mathrm{eff}}(\mathrm{})+\frac{cst}{L^\alpha }`$, one gets $`\theta _l^{\mathrm{eff}}(\mathrm{})=1.96(6)`$ for the PPBC, the best fit being also represented in Fig. (9). In this figure, we also show the same exponent obtained for the FPBC, where the best fit yields the asymptotic value $`\theta _l^{\mathrm{eff}}(\mathrm{})=2.12(11)`$. In both cases the fitting value we obtain for the exponent is $`\alpha 1`$. Note that the asymptotic values for $`\theta _l^{\mathrm{eff}}`$ are well compatible with our prediction of Sec. II C, $`\theta _l=2`$ (see also the heuristic argument in Appendix A).

## VI Finite-volume corrections (FVC) and their relation to the statistics of extreme values

What is the origin of these strong finite-volume corrections? Intuitively it is not difficult to find an explanation for the strong systematic finite-volume corrections in the lowest droplet exponent $`\theta _l`$. As the word lowest indicates, these exponents describe the statistical distribution of droplet excitations which are at the tail of the energy gap distribution that includes all possible high energy levels. As the volume of the system increases there is more available space to find excitations with lower energy gap. This implies that there is more probability to find a lowest droplet with an energy smaller than a given threshold value $`E^{}`$. As this probability systematically increases with $`L`$ for all samples, $`\theta _l^{\mathrm{eff}}(L)`$ must be a decreasing function of $`L`$. This is in agreement with what we have found. The behavior of $`\lambda _l^{\mathrm{eff}}(L)`$ is more difficult to establish. By the same token, although there is more available space for droplets we expect that different droplet sizes increase or decrease their relative probability in a non-trivial way making difficult to guess how the exponent $`\lambda _l^{\mathrm{eff}}(L)`$ systematically changes with $`L`$.

To understand the origin of finite-volume corrections in the value of $`\theta _l`$ we have focused our attention on the behavior of the upper bound exponent $`\theta _l^1`$ describing the statistics of the lowest one-spin excitations as described in the Appendix A. The gap distribution corresponding to these excitations can be obtained from the local-field distribution evaluated at the ground state. We have numerically computed this distribution for different sizes, the results are shown in Fig. 10. As discussed in the appendix A, the local-field distribution has a finite weight at zero field and is a self-averaging quantity. As the local field distribution is self-averaging, the probability distribution for the lowest one-spin excitations corresponds to the extreme value statistics of the local-field distribution $`p_1(h)`$ where $`h`$ stands for the local field which we assume to be positive as the gap is given by its absolute value (the subindex 1 is used to stress that this distribution describes energy gaps for one-spin excitations only). If $`P_1(h)`$ stands for the probability distribution of the smallest local-fields, then $`P_1(h)`$ can be easily related to $`p_1(h)`$ by standard probability arguments (see for instance, ). Although the argument is very general, here we apply it to one-spin excitations. For a given sample, the lowest value $`h`$ is selected as the minimum value among all the possible $`V`$ local fields $`h_i`$ at each lattice site. The probability $`P_1(h)`$ is given by the following expression,

$$P_1(h)=Vp_1(h)(1_h^{\mathrm{}}p_1(h^{})𝑑h^{})^{V1}=\frac{}{h}\left(_h^{\mathrm{}}p(h^{})𝑑h^{}\right)^V$$
(35)
which accounts for all possible ways the value $`h`$ coincides with the minimum value obtained among all different $`V`$ local fields distributed according to the $`p_1(h)`$. The last identity shows that $`P_1(h)`$ is normalized. This probability can be explicitly worked out in the large $`V`$ limit,

$$P_1(h)=\frac{}{h}\mathrm{exp}\left[Vg_1(h)\right]=Vg_1(h)\mathrm{exp}\left[Vg_1(h)\right].$$
(36)
Up to second order in $`h`$ the function $`g_1(h)`$ is given by,

$$g_1(h)=p_1(0)h+\frac{p_1^{}(0)+(p_1(0))^2}{2}h^2.$$
(37)
From (36) we immediately learn that the gap distribution is an exponential with a sub-leading Gaussian correction whose magnitude decreases as $`1/V`$. Actually, plotting $`P_1(h)/V`$ as function of the scaling variable $`x=hV`$ one gets,

$$\frac{P_1(h)}{V}=g_1^{}(x/V)\mathrm{exp}\left[xp_1(0)\frac{p_1^{}(0)+(p_1(0))^2}{2V}x^2\right].$$
(38)
In the large $`V`$ limit $`g_1^{}(x/V)p_1(0)`$ and the coefficient in front of the Gaussian correction goes asymptotically to zero, therefore the distribution $`P_1(h)`$ converges to an exponential as expected,

$$P_1(h)=Vp_1(0)\mathrm{exp}\left(Vp_1(0)h\right)$$
(39)
in agreement with the scaling relation (13). We can now understand the deviations from the pure exponential behavior discussed in Sec. III in the context of the inset a) shown in Fig. 3. They are simply consequence of the finite-volume corrections of the extreme values of the gap distribution for all energy levels (and not only one-spin excitations as we are discussing here). So one could imagine to compute for a given sample the first $`V`$ energy gaps corresponding to the first $`V`$ excitations. We also assume that the resulting distribution $`p_{\mathrm{all}}(E)`$ is self-averaging in the large volume limit (as it is the local-field distribution $`p_1(h)`$). The lowest energy distribution constructed by taking the minimum value of the gap $`E`$ for each sample yields the extremes distribution $`P_{\mathrm{all}}(E)`$ defined in (13) (from now on we will drop the subindex ’all’ in $`P_{\mathrm{all}}(E)`$ as it coincides with the $`P(E)`$ defined in (13). Also the subindex ’all’ for the $`p_{\mathrm{all}}(E)`$ will be dropped). In fact, the parameters $`p(0),p^{}(0)`$ characterizing this distribution can be obtained from the $`P(E)`$’s shown in Fig. 3. To evaluate them, the best way is to analyze the cumulative distribution $`𝒫(E)=_E^{\mathrm{}}𝑑E^{}P(E^{})`$ which from (35) we can assume to be $`𝒫(E)=\mathrm{exp}[Vg_{\mathrm{all}}(E)]`$. Thus we can fit $`𝒫(E)`$ with an exponential with Gaussian corrections $`A\mathrm{exp}[BxCx^2/2]`$ whose fitting parameters are related to $`p(0)`$ and $`p^{}(0)`$. The best fits yield the following values $`p(0)0.2`$ and $`p^{}(0)=0.3`$.

Coming back to our original goal we discuss now the finite-volume corrections for the estimate $`\theta _l^{\mathrm{eff}}`$, as shown in Fig. 4. From the distribution (36) describing the whole spectrum of excitations we can express the effective exponent (25) for $`L>>1`$ as,

$$\theta _l^{\mathrm{eff}}(L)=\frac{\mathrm{log}\left(\overline{E}(L)\right)}{\mathrm{log}(L)}.$$
(40)
The computation of $`\overline{E}(L)`$ is quite straightforward as it is given by the simple relation,

$$\overline{E}(L)=_0^{\mathrm{}}EP(E)dE=_0^{\mathrm{}}\mathrm{exp}(Vg(E))$$
(41)
where we have used (36) plus an integration by parts. The integral, up to second order in $`1/V`$ yields,

$$\overline{E}=\frac{1}{Vp(0)}\frac{p^{}(0)+(p(0))^2}{V^2(p(0))^3}+O\left(\frac{1}{V^3}\right).$$
(42)
Inserting this result in (40) we finally get,

$$\theta _l^{\mathrm{eff}}(L)=d+\frac{d}{V}\left(1+\frac{p^{}(0)}{(p(0))^2}\right)+O\left(\frac{1}{V^3}\right).$$
(43)
This shows that $`\theta _l^{\mathrm{eff}}(L)`$ approaches $`d`$ from below (as $`p^{}(0)`$ is positive). On the other hand the magnitude of the finite-volume corrections can be pretty large if $`\frac{p^{}(0)}{(p(0))^2}>>1`$. For instance, if one takes the results obtained from the analysis of one-spin excitations one gets $`p_1(0)0.069,p_1^{}(0).125`$ yielding $`\frac{p_1^{}(0)}{(p_1(0))^2}27.`$ which is indeed large. Inserting these values in (43) we obtain an estimation for $`\theta _l^{\mathrm{eff}}(L=12)=1.65`$ in good agreement with numerical results (see Fig. 11).

If we insert the previous estimated values for the whole spectrum of excitations extracted from the $`P(E)`$’s in Fig. 3, we obtain $`\frac{p^{}(0)}{(p(0))^2}7.5`$. From (43) it follows that $`\theta _l^{\mathrm{eff}}(L)2(17.5/V)`$, which for $`L=11`$ yields $`\theta _l^{\mathrm{eff}}=1.87`$. All in all, the magnitude of the effective exponent $`\theta _l`$ is well compatible with the reported value $`\theta _l^{\mathrm{eff}}`$ used in the inset a) in Fig. 3 for the PP case. Note that the FVC corrections to $`\theta _l^{\mathrm{eff}}`$ obtained from the local-field distributions in the FF case are much larger than FVC corrections in the PP case in agreement with ARA results (see Fig. 9). From this analysis it becomes clear that to significantly reduce the magnitude of the finite-volume corrections in the value of $`\theta _l`$ (let us say $`\theta _l1.95`$), we would need larger volumes beyond $`20\times 20`$.

## VII Compactness of the lowest energy droplets.

One intriguing question about the droplet excitations concerns their topological properties. Kawashima and Aoki  have argued that droplet excitations are not compact. Instead, their volume has a fractal structure as the number of lattice points included in the droplet scales with its spanning length (which is a measure of the length scale of the droplet) with an exponent smaller (around 1.80(2)) than the dimension of the system (2).

To answer this question we have computed the surface, i.e. the perimeter $`P`$, of all lowest droplets. The relation between the average perimeter as function of the size $`v`$ of the excitation depends on both the fractal dimension of the surface or perimeter $`d_s`$ and the volume $`d_v`$ of the lowest droplets. $`d_s`$ and $`d_v`$ can be defined in terms of the spanning length $`l`$ of the droplet which can be defined in different ways. For example, one could use the gyration radius, the average distance between the sites contained in the cluster, or the maximum distance among the sites of the cluster. As the typical length scale of our lowest droplets is small, $`l10`$, we have not attempted to estimate it as this can strongly depend on the precise definition of the spanning length. Here, we restrict ourselves to investigate the perimeter-volume dependence. In terms of the spanning length $`l`$ the surface fractal and volume fractal dimensions $`d_s,d_v`$ of the droplets are defined as,

$`lP^{\frac{1}{d_s}}`$ (44)
$`lv^{\frac{1}{d_v}}`$ (45)
which combined give,

$$Pv^{\frac{d_s}{d_v}}.$$
(46)
In Fig. 12 we show $`P(v)`$ as a function of $`v`$ for different lattice sizes. As can be seen, FVC are important for large volumes. However, there is an enveloping curve that is independent of $`L`$ for small volumes and spans a progressively increase range of volumes as $`L`$ increases. This enveloping curve is excellently fitted (continuous curve) by the scaling relation (46) and yields an estimate,

$$\frac{d_s}{d_v}0.632(2)$$
(47)
consistent with the results reported by Kawashima and Aoki $`\frac{d_s}{d_v}=0.61(1)`$ obtained with a completely different method.

## VIII Conclusions

We have shown that a proper description of low-temperature properties in two-dimensional Gaussian spin glasses can be done in terms of the lowest droplet exponents (LDEs) $`\lambda _l`$ and $`\theta _l`$ describing the spectrum of lowest excitations. $`\lambda _l`$ describes the spectrum of sizes of the lowest energy droplets, while $`\theta _l`$ describes the typical energy cost of these lowest droplets whatever their size. Assuming that $`\theta _l=d`$ one concludes that the LDE $`\lambda _l`$ fully characterizes the spin-glass phase. Although independent numerical estimates of $`\theta _l`$ and $`\lambda _l`$ show strong finite-volume corrections, the thermal exponent $`\theta =\theta _l+d\lambda _l`$ can be well estimated giving the results (29), (30)

$$\theta =0.46(1)\lambda _l=0.770(5).$$
(48)
unambiguously showing that $`\theta <\theta _{DW}=0.287(4)`$ . Our estimates (48) have been confirmed via an aspect-ratio analysis which provides estimates much less influenced by finite-volume corrections. Moreover, the result $`\theta _l=2`$ (that is believed to be correct for spin glasses with coupling distributions with finite weight at zero coupling, see the Appendix A) has been numerically confirmed by the aspect-ratio analysis. To sum up, McMillan’s excitations are not the typical low-lying excitations and our approach offers a new and independent way to estimate the thermal exponent $`\theta `$ without the need to generate typical low-lying excitations by looking at the new ground state of the system after perturbing it.

We think that discrepancies on the value of the thermal exponent $`\theta `$ reported by comparing non-perturbative methods (such as finite-temperature transfer-matrix calculations and the present lowest droplet analysis) with perturbative methods such as domain-wall calculations (or perturbations induced by introducing a coupling term in the energy function that induces a large-scale excitation) are serious enough to be taken as a clear indication that our knowledge of the low-temperature properties of the 2d GISG is still inadequate. In this direction we want also to recall the issue of multifractality and the possibility that different exponents could describe the zero-temperature critical point. Is this really possible? Well, to our knowledge no exact result precludes this possibility and, although purely speculative at the present stage, one should seriously think about it. Altogether, the present analysis suggests that the excitations in 2d GISG are very different from the compact droplets proposed in the context of the droplet model. If this were true, the implications of the 2d studies in larger dimensions could be important. There are many routes that can be followed to understand better what is going on and the origin of this discrepancy. Certainly, with the outstanding accuracy of present algorithms to compute ground states in 2d, it would be very interesting to revisit again the analysis of the statistics of the large scale excitations generated by imposing a uniform magnetic field. “Old” results by Rieger et al.  give an estimate for $`\theta `$ that is compatible with our estimate rather than to the domain-wall estimate. This would be an independent check of our values, but using a perturbation method with an appropriate neutral observable such as the global magnetization as has been explained in Sec. I before (1).

The proposed method may appear venturesome as, to our present knowledge, there is no numerical study in the field along this line of research. However, as explained in Sec. I, recent studies on the disordered Anderson model  have revealed that the analysis of the lowest excitation provides a good description of the localized phase. More studies are certainly required to understand better the reliability of the present method to investigate the critical properties of spin glasses. One disadvantage of our approach is that a huge number of samples is needed to reasonably sample large-scale excitations. However, as we saw in Sec. V, the behavior of the $`g(q)`$ for small $`1q`$ can be extracted with a modest number of samples. The advantage, as has been already stressed in Sec. I, is that we do not introduce any external perturbation to generate the excitations.

Finally, we want to comment on the extension of this approach to other models. Of course, the immediate extension one could think of is the 2d $`\pm J`$ model. However, the analysis of this model appears quite troublesome. This model does not have a continuous gap distribution but a discrete one that introduces further complications. As the ground state is not unique one has to redefine the full analysis to properly define the spectrum of lowest excitations. The discreteness of variables could have some unexpected effects in the present approach as seems to happen also with domain-wall calculations . It is more natural to extend the research to other models such as 2d ISG with other continuous coupling distributions without gap (e.g. characterized by $`P(J)|J|^\alpha `$ for $`|J|0`$), Migdal-Kadanoff spin glasses (where both the ground state and the first excitation could be feasibly found with an appropriate algorithm), Gaussian spin glasses beyond $`d=2`$ (where unfortunately, algorithms are much less effective than in 2d as the finding of the ground state becomes a NP complete problem) and finally mean-field spin-glass models where the zero temperature exponents are known and maybe the spectrum of lowest excitations could be analytically tackled. Preliminary results in this case  confirm that the present analysis describes pretty well the data for rather small sizes. We are pretty confident that, in the near future, new results and evidences will finally resolve this interesting problem.

## A Heuristic proof of the identity $`\theta _l=d`$

In this appendix we show that $`\theta _l=d`$. In what follows we do not attempt to present a rigorous proof but we content ourselves to present an heuristic argument. The argument has two parts: first we show that $`d`$ is an upper bound, next we show that the upper bound is the exact value. For the upper bound the argument is well known and goes as follows. Consider the ground state and all possible one-spin excitations. Because one-spin excitations are not necessarily the absolute lowest ones, the statistics of the lowest one-spin excitations must yield an upper bound $`\theta _l^1`$ for the value of $`\theta _l`$, $`\theta _l\theta _l^1`$. The statistics of the lowest one-spin excitations is determined by the behavior of the ground-state local field distribution $`p(h)`$ in the limit $`h0`$. If $`p(h)`$ is self-averaging and $`p(0)`$ is finite (in the large-volume limit) then the statistics of the lowest excitations must be governed by the exponent $`\theta _l^1=d`$. Although we do not know a precise mathematical proof of the statement that $`p(0)`$ is finite, it looks quite intuitive . In any short-range system with a frustrated ground state and a coupling distribution with finite density at zero coupling, we may expect a finite probability to find a cage containing a spin coupled to its neighbors by a set of weak bonds which produce a vanishing net local-field acting on that spin. This argument should generally hold for $`d2`$. Moreover, as its name indicates, the local-field distribution is a local observable. An argument à la Brout proves that it should be self-averaging as all possible local field values are realized across the whole lattice <sup>*</sup><sup>*</sup>*Our numerical results in the 2d GISG confirm this conclusion, see Sec. VI. The next part of the argument consists in proving that an identical upper bound is valid by considering excitations with size strictly larger than 1 but finite. The upper bound derived for the one-spin excitations must necessarily hold for finite-size excitations beyond one-spin excitations (for instance two-spins, three-spins and so on) as the gap corresponding to the finite-size excitations can always be written as a linear combination of a finite number of local fields with coefficients which depend on the ground state configuration. It is easy to verify that the aforementioned properties of the local-field distribution $`p(h)`$ imply that the new gap distribution has a finite weight at zero gap and is self-averaging. This argument however cannot be extended to large-scale (with $`vV`$) excitations in a straightforward way because the distribution for the corresponding gap distribution corresponds to an infinite sum of terms in the $`V\mathrm{}`$ limit. However once we argue that $`\theta _l`$ is an upper bound valid for all finite-size excitations it can be concluded that this upper bound must coincide with the exponent $`\theta _l`$ describing the probability of the absolute lowest excitations. From (4) the fraction of large scale excitations $`vV`$ is given by $`Vg_V=\frac{1}{V^\lambda }`$. In general $`\lambda >0`$ so this fraction vanishesThis fractions is finite only in $`d=1`$ where $`\lambda =0`$. But this case is trivial as the surface of large scale droplets in $`d=1`$ only contains a finite number of broken bonds. in the infinite-volume limit and finite-size excitations determine the result $`\theta _l=d`$ as they dominate the spectrum of lowest excitations. Moreover, if large-scale excitations yield a different value for $`\theta _l`$ this would imply that boundary conditions could affect the value of the thermal exponent. That would be quite unusual as this would mean that the exponents of the $`T=0`$ fixed point would depend on the boundary conditions.

## B Transfer Matrix algorithms

In this appendix, we will briefly explain how we determine the ground state and the first excited state. We will work on a square lattice of size $`L\times L`$. The energy associate to a configuration of spins $`S(i,j)`$ with a fixed configuration of disorder $`J^x(i,j)`$ and $`J^y(i,j)`$ is

$`E={\displaystyle \underset{i=1,L1}{}}{\displaystyle \underset{j=1,L}{}}J^x(i,j)S(i,j)S(i+1,j)+{\displaystyle \underset{i=1,L}{}}{\displaystyle \underset{j=1,L1}{}}J^y(i,j)S(i,j)S(i,j+1)`$ (B1)
$`+B_1{\displaystyle \underset{j=1,L}{}}J^x(L,j)S(L,j)S(1,j)+B_2{\displaystyle \underset{i=1,L}{}}J^y(i,L)S(i,L)S(i,1),`$ (B2)
where $`B_1`$ and $`B_2`$ correspond to the choice of boundary conditions. Here we will consider three cases : Periodic-Periodic boundary conditions (PPBC) with $`B_1=B_2=1`$; Free-Periodic boundary conditions (FPBC) with $`B_1=1;B_2=0`$ (or equivalently $`B_1=0;B_2=1`$) and Free-Free boundary conditions (FFBC) $`B_1=B_2=0`$. We will only consider the case with a Gaussian distribution of the bond disorder $`J^x,J^y`$. To determine the ground state and the first excited states, we proceed as follow: we start by associating a weight for each configurations of spins in the first row of the lattice $`S(1,1),S(1,2),\mathrm{},S(1,L)`$:

$$W(S(1,1),S(1,2),\mathrm{},S(1,L))=B_2J^y(1,L)S(1,L)S(1,1)+\underset{i=1,L1}{}J^y(1,i)S(1,i)S(1,i+1).$$
(B3)
Next, we start iterating the transfer matrix using a sparse-matrix factorization . The first iteration gives

$$W^1(S(1,2),\mathrm{},S(1,L),S(2,1))=max_{S(1,1)}\left[J^x(1,1)S(1,1)S(2,1)+W(S(1,1),\mathrm{},S(1,L))\right].$$
(B4)
Since we are also interested in the first excited state, we define the second largest weight:

$$W^2(S(1,2),\mathrm{},S(1,L),S(2,1))=min_{S(1,1)}\left[J^x(1,1)S(1,1)S(2,1)+W(S(1,1),\mathrm{},S(1,L))\right].$$
(B5)
In the following, we will use the simplified notation

$$W(i,j)W(S(i,j),\mathrm{},S(i,L),S(i+1,1),\mathrm{},S(i+1,j1)).$$
(B6)
Thus eqs.(B4,B5) become

$$W^1(1,2)=max_{S(1,1)}\left[J^x(1,1)S(1,1)S(2,1)+W(1,1)\right]$$
(B7)
and

$$W^2(1,2)=min_{S(1,1)}\left[J^x(1,1)S(1,1)S(2,1)+W(1,1)\right].$$
(B8)
At the next iteration, the numbers of possible weight will again be multiplied by two but we will keep only the two largest ones defined as

$`W^1(1,3)=max_{S(1,2)}\left[J^x(1,2)S(1,2)S(2,2)+J^y(2,1)S(2,1)S(2,2)+W^1(1,2)\right]`$ (B9)
$`W^2(1,3)`$ $`=`$ $`max(min_{S(1,2)}[J^x(1,2)S(1,2)S(2,2)+J^y(2,1)S(2,1)S(2,2)+W^1(1,2)],`$ (B11)
$`max_{S(1,2)}[J^x(1,2)S(1,2)S(2,2)+J^y(2,1)S(2,1)S(2,2)+W^2(1,2)]).`$
The general iteration relations are

$`W^1(i,j)`$ $`=`$ $`max_{S(i,j1)}[J^x(i,j1)S(i,j1)S(i+1,j1)`$ (B13)
$`+J^y(i+1,j2)S(i+1,j2)S(i+1,j1)+W^1(i,j1)]`$
$`W^2(i,j)=max`$ $`(`$ $`min_{S(i,j1)}[J^x(i,j1)S(i,j1)S(i+1,j1)`$ (B17)
$`+J^y(i+1,j2)S(i+1,j2)S(i+1,j1)+W^1(i,j1)],`$
$`max_{S(i,j1)}[J^x(i,j1)S(i,j1)S(i+1,j1)`$
$`+J^y(i+1,j2)S(i+1,j2)S(i+1,j1)+W^2(i,j1)]).`$
In addition, each time that we end the construction of a new row, we must add the boundary term

$$W^1(i,1)W^1(i,1)+B_2J^y(i,L)S(i,L)S(i,1);W^2(i,1)W^2(i,1)+B_2J^y(i,L)S(i,L)S(i,1).$$
(B18)
We still have to take in account the boundary condition corresponding to $`B_1`$. The two type of boundary conditions (Free or Periodic) have to be considered separately :

* Free boundary condition, $`B_1=0`$ : In that case, we iterate up to the construction of the weights associated to the configurations of the spins $`S(L,1),\mathrm{},S(L,L)`$. The energy of the ground state ($`E_0`$) is then simply the maximum among all the weights $`W^1(L,1)`$:

$$E_0=max_{\{S(L,1),\mathrm{},S(L,L)\}}[W^1(L,1)].$$
(B19)
We call $`\{S^0(L,1),\mathrm{},S^0(L,L)\}`$ the configuration of spins on the last row for the ground state. The energy of the first excited state is the second largest weight among $`W^1(L,1)`$ and $`W^2(L,1)`$:

$$E_1=max[max_{\{S(L,1),\mathrm{},S(L,L)\}\{S^0(L,1),\mathrm{},S^0(L,L)\}}[W^1(L,1)],max_{\{S(L,1),\mathrm{},S(L,L)\}}[W^2(L,1)]].$$
(B20)
* Periodic boundary case, $`B_1=1`$ : We first choose one configuration of spins on the first row $`S^i(1,1),S^i(1,2),\mathrm{},S^i(1,L)`$. The weight of this configuration is defined as in (B3). The weight of all the other configurations of spins on this first row are fixed to an arbitrary large negative number. Next we iterate the transfer matrix as described above, up to the construction of the weights $`W^1(L,1)`$ and $`W^2(L,1)`$. Finally, we iterate one additional row, with bonds $`J^x(L,i)`$ and $`J^y(L+1,i)=0`$. Next, we store the two weights $`W^1`$ and $`W^2`$ associated to the initial spins configuration $`S^i(1,1),S^i(1,2),\mathrm{},S^i(1,L)`$. We denote these two weights by

$`𝒲^1`$ $`(S^i(1,1),S^i(1,2),\mathrm{},S^i(1,L))`$ (B22)
$`W^1(S(L+1,1)=S^i(1,1),S(L+1,2)=S^i(1,2),\mathrm{},S(L+1,L)=S^i(1,L))`$
$`𝒲^2`$ $`(S^i(1,1),S^i(1,2),\mathrm{},S^i(1,L))`$ (B24)
$`W^2(S(L+1,1)=S^i(1,1),S(L+1,2)=S^i(1,2),\mathrm{},S(L+1,L)=S^i(1,L)).`$
The energy of the ground state is the maximum on all the $`𝒲^1`$:

$$E_0=max_{\{S^i(1,1),\mathrm{},S^i(1,L)\}}[𝒲^1(S^i(1,1),S^i(1,2),\mathrm{},S^i(1,L))],$$
(B25)
and we denote by $`S^0(1,1),\mathrm{},S^0(1,L)`$ the configuration of spins on the first row for the ground state. The energy of the first excited state is the second largest weight among $`𝒲^1(L,1)`$ and $`𝒲^2(L,1)`$:

$`E_1=max[max_{\{S^i(1,1),\mathrm{},S^i(1,L)\}\{S^0(1,1),\mathrm{},S^0(1,L)\}}[𝒲^1(S^i(1,1),S^i(1,2),\mathrm{},S^i(1,L))],`$ (B26)
$`max_{\{S^i(1,1),\mathrm{},S^i(1,L)\}}[𝒲^2(S^i(1,1),S^i(1,2),\mathrm{},S^i(1,L))]].`$ (B27)
The construction of the ground state and of the first excited state is much more costly in computing time for the periodic case, since we have to repeat $`2^L`$ times the iterations, one time for each configuration $`S^i(1,1),S^i(1,2),\mathrm{},S^i(1,L)`$.

So far, we have only described how to compute the value of the energies associated to the ground state and the first excited state. Since we also want to determine the spins configurations for these two states, we have to store, at each iteration of the transfer matrix, the value of the spin on which one sums, as well as the value of the previous spins. Thus, for each $`2^L`$ weights $`W^1(i,j1)`$, we have to store the configuration

$$𝒞_{S(i,j1),\mathrm{},S(i+1,j2)}(S(1,1),\mathrm{},S(1,L),S(2,1),\mathrm{},S(i,j2)).$$
(B28)
At the next iteration, we will build the weight $`W^1(i,j)`$ with the corresponding configuration

$$𝒞_{S(i,j),\mathrm{},S(i+1,j1)}=(𝒞_{S(i,j1),\mathrm{},S(i+1,j2)},S(i,j1))$$
(B29)
with $`S(i,j1)`$ the value of the spin which corresponds to the maximum in (B13). From this construction, one have access to the spins configurations of the ground state and the first excited state.

Finally, one should also add that this construction can be easily extended to the second excited state, etc. After (B11), one can easily define a third weight which would be associated to the second excited state, and so on.

Acknowledgments. We are grateful to J.-P. Bouchaud, D. Huse, O. Martin, M. A. Moore and A. P. Young for useful comments. We have been supported by the Spanish Ministerio de Ciencia y Tecnología, projects BFM2001-3525 (F.R.) and grant AP98-36523875 (M.S.). M. P. and F. R. acknowledge support from the French-Spanish collaboration (Picasso program and Acciones Integradas HF1998-0097). Funding from the European Science Foundation through the SPHINX program is also acknowledged.


text:
# Neutrino-nucleon scattering rate in the relativistic random phase approximation

## 1 Introduction

Neutrino transport is an important ingredient of the numerical simulation of supernova collapse and of the cooling of the resulting protoneutron star. In the conditions of high density and temperature prevailing in this situation, the neutrino opacities will be appreciably modified by medium effects. This has been confirmed by the calculations of several groups (see e.g.  and references therein).

We perform here a calculation of the effect of correlations on this parameter in the random phase approximation. This contribution is an extension of the results presented in  to the more commonly used parametrizations of the nuclear interactions NL3 and TM1 . The reader is referred to  for detailed analytical expressions as well as for a more complete list of bibliographical references.

In the dense and hot nuclear matter present in the protoneutron star, the mean free path of the neutrino is determined mostly by its interactions with the baryons. Here we will consider only the nucleons. The scattering process occurs via the neutral current, while the emission and absorption processes involve the charged current. In both cases, the differential cross section is given by

$`{\displaystyle \frac{d\sigma }{dE_\nu d\mathrm{\Omega }}}={\displaystyle \frac{G_F^2}{64\pi ^3}}{\displaystyle \frac{E_\nu ^{}}{E_\nu }}m(S^{\mu \nu }L_{\mu \nu })`$ (1)
It is proportional to the imaginary part of the lepton current $`L^{\mu \nu }`$ contracted with the structure function of the nucleons $`S^{\mu \nu }`$. When the nucleons are free, the structure function is proportional, with a thermodynamical factor, to the mean field value of the polarization

$`\mathrm{\Pi }_{WW}^{\mu \nu \text{MF}}(q)={\displaystyle d^4p\mathrm{Tr}\left[\mathrm{\Gamma }_W^\mu G(p)\mathrm{\Gamma }_W^\nu G(p+q)\right]}`$ (2)
where $`G`$ is the propagator of the nucleon and $`\mathrm{\Gamma }_W`$ is the vertex for the weak interaction.

The RPA correlations are introduced by replacing the mean field polarization by the solution of the Dyson equation

$`\mathrm{\Pi }_{WW}^{\mu \nu \text{RPA}}=\mathrm{\Pi }_{WW}^{\mu \nu \text{MF}}+\mathrm{\Pi }_{WS}^{\mu \alpha \text{MF}}D_{SS\alpha \beta }^{(0)}\mathrm{\Pi }_{SW}^{\beta \nu \text{RPA}}`$ (3)
with $`D_{SS}^{(0)}`$ representing the nuclear interaction in free space111The indices $`W`$ or $`S`$ stand for a vertex with a weak or strong coupling respectively. As can be easily seen on the diagramatical representation, this corresponds to an infinite series of loop diagrams, which may also be resummed by defining a screened nuclear interaction $`D_{SS}^{\mathrm{RPA}}`$

$`\mathrm{\Pi }_{WW}^{\mu \nu \text{RPA}}=\mathrm{\Pi }_{WW}^{\mu \nu \text{MF}}+\mathrm{\Pi }_{WS}^{\mu \alpha \text{MF}}D_{SS\alpha \beta }^{\text{RPA}}\mathrm{\Pi }_{SW}^{\beta \nu \text{MF}}`$ (4)
## 2 Polarization insertions

We will describe the nuclear interaction by a relativistic field model with meson exchange. Non-linear meson-meson interactions are introduced in order to obtain a better description of the properties of nuclear matter at saturation. We consider here two widely used parametrizations of the nuclear interaction, NL3 and TM1 , which are fitted to the properties of nuclei and bulk nuclear matter. The interaction piece of the corresponding Lagrangian reads

$`_{\mathrm{int}}=\overline{\psi }(g_\sigma \sigma +g_\omega \gamma ^\mu \omega _\mu +g_\rho \gamma ^\mu \stackrel{}{\rho }_\mu .\stackrel{}{\tau })\psi {\displaystyle \frac{1}{3}}bm_N\sigma ^3{\displaystyle \frac{1}{4}}c\sigma ^4+{\displaystyle \frac{1}{4}}d(\omega _\alpha \omega ^\alpha )`$ (5)
The values of the parameters are collected in the following table.

These interactions do not include the $`\delta `$ meson considered in  nor the tensor coupling of the $`\rho `$ meson. The Lagrangian also does not include a pion term; however it was shown in  that the neutrino-nucleon scattering rate does not depend on the pion and only very weakly on the details of the residual contact interaction222 This statement does not hold for the emission and absorption processes which involve charged $`\rho `$ and $`\pi `$ exchange, corrected at short distance by contact terms.. These parametrizations have also been shown to yield good agreement with the available data on the longitudinal response for quasielastic electron scattering on nuclei in the RPA approximation. (cf. e.g. ).

We present in this contribution results for the neutrino-nucleon scattering rate. In this case, we need to replace the piece $`D^{\mathrm{RPA}}`$ describing the screened nuclear interaction by the matrix of propagators for the neutral mesons $`\sigma `$, $`\omega `$ and $`\rho `$ dressed at RPA level in asymmetric nuclear matter. The detailed expression for these propagators with meson mixing in all channels was derived in .

One should be aware that, in models with non linear meson couplings, the meson propagators are not only dressed by particle-hole loops, but also include a contribution of the non linearities. When performing the derivation of the propagators, one finds that this amounts to replacing the masses of the mesons appearing in the expressions of the propagators given in Eqs. (32-47) of  as follows333 The index H indicates that the values of the meson fields are taken in the Hartree approximation, and $`\eta ^2=\eta ^\mu \eta _\mu `$ with $`\eta ^\mu `$ defined below Eq. (6):

$`m_\sigma ^2`$ $``$ $`M_\sigma ^2=m_\sigma ^2+2bm\sigma _H+3c\sigma _H^2`$
$`m_\omega ^2`$ $``$ $`M_{\omega T}^2=m_\omega ^2+d\omega _H^2\text{(trans.) or}M_{\omega L}^2={\displaystyle \frac{(m_\omega ^2+d\omega _H^2)(m_\omega ^2+3d\omega _H^2)}{m_\omega ^2+d\omega _H^2\left(2\eta ^2+3\right)}}(\mathrm{long}.)`$
## 3 Numerical results

We calculated the differential and total cross sections in asymmetric matter with a proton fraction $`Y_p`$. The proton fraction is determined by the condition that $`\beta `$ equilibrium is realized: $`\widehat{\mu }=\mu _n\mu _p=\mu _e\mu _\nu `$. In cold neutron stars, the neutrinos can leave the star unhindered. On the other hand, in the conditions we are considering here they are still trapped inside the star, and the chemical potential of the neutrino therefore takes a finite value $`\mu _\nu =(6\pi ^2\rho Y_\nu )^{1/3}`$. The lepton fraction $`Y_L=Y_e+Y_\nu `$ is determined by neutrino transport (diffusion equation or Boltzman), to which the cross section serves as input. A typical value is $`Y_L0.4`$.

It is convenient to decompose the polarizations and propagators onto orthogonal projectors formed with the vectors and tensor available in the problem, i.e. the metric $`g^{\mu \nu }=\mathrm{diag}`$ (1,-1,-1,-1), the hydrodynamic velocity $`u^\mu `$ and the transferred momentum $`q^\mu `$.

$`\mathrm{\Pi }_{WW}^{\mu \nu \text{RPA}}`$ $`=`$ $`\mathrm{\Pi }_TT^{\mu \nu }+\mathrm{\Pi }_L\mathrm{\Lambda }^{\mu \nu }+\mathrm{\Pi }_QQ^{\mu \nu }+i\mathrm{\Pi }_EE^{\mu \nu }`$ (6)
$`\mathrm{\Lambda }^{\mu \nu }`$ $`=`$ $`{\displaystyle \frac{\eta ^\mu \eta ^\nu }{\eta ^2}};\eta ^\mu =u^\mu {\displaystyle \frac{q.u}{q^2}}q^\mu (\mathrm{longitudinal})`$
$`T^{\mu \nu }`$ $`=`$ $`g^{\mu \nu }{\displaystyle \frac{\eta ^\mu \eta ^\nu }{\eta ^2}}{\displaystyle \frac{q^\mu q^\nu }{q^2}}(\mathrm{transverse})`$
$`E^{\mu \nu }`$ $`=`$ $`ϵ^{\mu \nu \rho \lambda }\eta _\rho q_\lambda (\mathrm{axial});Q^{\mu \nu }={\displaystyle \frac{q^\mu q^\nu }{q^2}}\text{(does not contribute)}`$
The contraction of the lepton current with the polarization can be expressed by means of three structure functions $`R_1`$, $`R_2`$ and $`R_5`$ related to the previous polarizations by

$`R_1`$ $`=`$ $`{\displaystyle \frac{2}{1e^z}}m[{\displaystyle \frac{q^2}{𝐪^\mathrm{𝟐}}}\mathrm{\Pi }_L+{\displaystyle \frac{w^2}{𝐪^\mathrm{𝟐}}}\mathrm{\Pi }_T],R_2={\displaystyle \frac{2}{1e^z}}m\left[\mathrm{\Pi }_T\right]`$
$`R_5`$ $`=`$ $`{\displaystyle \frac{2}{1e^z}}m\left[\mathrm{\Pi }_E\right]\mathrm{with}z={\displaystyle \frac{\omega \mathrm{\Delta }\mu }{k_BT}}\mathrm{and}q^\mu =(\omega ,𝐪)`$ (7)
The structure function $`R_1`$ involves the $`\sigma `$ meson and the longitudinal part of the $`\omega `$ meson, as well as the longitudinal part of the $`\rho `$. On the other hand the structure functions $`R_2`$ and $`R_5`$ involve the transverse part of the $`\omega `$ and $`\rho `$ mesons. The transverse contribution is only moderately modified by RPA correlations whereas the longitudinal contribution is reduced approximately by a factor of two. The contribution of the axial-vector response function is significant only at high density or energy. The relative magnitude of the contributions of the longitudinal, transverse and axial vector polarizations to the differential cross section were compared. The transverse contribution is dominant and will provide for about 60% of the total result. The corrections will therefore arise mostly from the subdominant longitudinal and axial-vector polarizations $`\mathrm{\Pi }_L`$ and $`\mathrm{\Pi }_E`$. As an exemple we show here the transverse structure function for the TM1 and NL3 parametrizations.

An estimate444This is only an estimate since other processes contribute. of the mean free path is given by the inverse of the total cross section obtained by integrating the differential one over the allowed energy and momentum transfer, with a Pauli blocking factor for the outgoing lepton. For the parameter sets NL3 and TM1, the total neutrino-neutron scattering cross section is found to be reduced by RPA correlations at high density by a factor $`20\%`$ to $`25\%`$. At low density and moderate temperature, on the other hand, RPA correlations would yield an enhancement; however the validity of the model becomes questionable in this range. The reduction factor obtained here is somewhat stronger than the one quoted in . Actually we observe that parameter sets adjusted to give a lower effective mass at saturation yield stronger reduction factors555$`m^{}`$=0.63 for TM1 and 0.6 for NL3 while it was adjusted to 0.8 in .