# Exploring different regimes in finite-size scaling of the droplet condensation-evaporation transition

###### Abstract

We present a finite-size scaling analysis of the droplet condensation-evaporation transition of a lattice gas (in two and three dimensions) and a Lennard-Jones gas (in three dimensions) at fixed density. Parallel multicanonical simulations allow sampling of the required system sizes with precise equilibrium estimates. In the limit of large systems, we verify the theoretical leading-order scaling prediction for both the transition temperature and the finite-size rounding. In addition, we present an emerging intermediate scaling regime, consistent in all considered cases and with similar recent observations for polymer aggregation. While the intermediate regime locally may show a different effective scaling, we show that it is a gradual crossover to the large-system scaling behavior by including empirical higher-order corrections. This implies that care has to be taken when considering scaling ranges, possibly leading to completely wrong predictions for the thermodynamic limit. In this study, we consider a crossing of the phase boundary orthogonal to the usual fixed temperature studies. We show that this is an equivalent approach and, under certain conditions, may show smaller finite-size corrections.

## I Introduction

The limit of large systems, or the thermodynamic limit, is of general interest when studying finite systems. It is this limit, that is commonly accessible to experiments. Recent developments in experimental and simulation techniques open the door to mesoscopic length scales MuellerPlathe2002 (). Here, general concepts of finite-size scaling may be tested and verified. Considering, for example, the equilibrium properties of several homopolymers, the aggregation transition temperature was shown to exhibit systematic finite-size effects that deviate from the expected behavior for particle gas systems Zierenberg2014JCP (). It was argued that the observed scaling is valid for intermediate system sizes, where the aggregate includes most of the polymers, but in the limit of increasing polymer number the particle picture should be recovered. On the other hand, it is a reasonable assumption that the intermediate regime should be also apparent for particle gas condensation, which will be a main focus of the present paper.

First-order phase transitions in spin systems are usually separating homogeneous phases, which leads in general to finite-size scaling corrections on the order of the inverse system volume. Exceptions occur for example if the low-temperature phase is exponentially degenerated Mueller2014 (). The general situation changes when we consider a phase transition between a homogeneous phase and a mixed phase – as in the case of the droplet condensation-evaporation transition separating a supersaturated gas phase and a mixed phase consisting of a single droplet in equilibrium with surrounding vapor. There exists a large amount of theoretical literature on this topic about the leading-order scaling behavior of large systems Biskup2002 (); Neuhaus2003 (); Binder2003 () with origins already in the 80s Binder1980 (). Usually, a system at fixed temperature is considered, estimating the transition density at which a single macroscopic droplet forms. In the limit of very large systems, this transition density is supposed to have finite-size corrections of the order , while the rounding of the transition should scale with . The former was numerically verified to leading-order at fixed temperature for a lattice gas Neuhaus2003 (); Nussbaumer2006 (); Nussbaumer2008 (); Zierenberg2014JPCS (). Several studies of a three-dimensional Lennard-Jones gas Binder2004 (); MacDowell2006 (); Binder2009 () also verified the finite-size corrections of the transition density. For a fixed density, the scaling of the transition temperature was demonstrated for the two- and three-dimensional lattice gas Malakis2007 (), while the scaling exponent of the finite-size rounding could not be verified. There exist complementary microcanonical studies of the Ising model at fixed density (or magnetization) Pleimling2001 (); Pleimling2009 (), demonstrating that the occurring transition shows signatures of a first-order transition.

The present study aims to fill the gap with respect to the finite-size scaling of the droplet condensation-evaporation transition at fixed density and the gradual crossover from an effective intermediate regime to the asymptotic large-system limit. In order to draw general conclusions, we consider a lattice-gas model in two and three dimensions as well as a Lennard-Jones gas in three dimensions. We will show that the proposed scaling behavior of the transition temperature and rounding is numerically recaptured in the limit of large (but achievable) system sizes for all dimensions and models considered. Where possible, we will directly compare to analytic solutions or low-temperature series expansions. In addition, we discuss an emerging intermediate scaling regime, best visible in the rounding of the transition and consistent with recent observations of finite-size effects in polymer aggregation Zierenberg2014JCP (). This regime with effective local scaling behavior and the gradual crossover to the large-system regime may be described by the finite-size scaling behavior when including empirical higher-order correction terms.

The remaining part of this paper is organized as follows. In Sec. II we briefly recapture the main leading-order results from the literature at fixed temperature and convert them to the scenario of fixed density. After describing our models and methods in Sec. III, we will present our results in Sec. IV and finish with the conclusions in Sec. V.

## Ii Theory

We consider a -dimensional liquid-vapor system of particles in a (periodic) box of volume , where for sufficiently dilute systems the condensation-evaporation transition separates a homogeneous supersaturated gas phase from a mixed phase of a droplet in equilibrium with surrounding bulk gas. At fixed density, this occurs at the condensation-evaporation transition temperature . However, the usual description of this problem is in the language of a grand-canonical scheme, considering a fixed temperature and variable density (or in case of a lattice gas equivalently the magnetization of the Ising model) with the finite-size transition density . This is orthogonal to the scenario of a fixed density with variable temperature. Figure 1 shows a sketch of the infinite-size phase boundary or and the different scenarios crossing it for finite system sizes. It needs to be mentioned that both schemes are working in the canonical ensemble for each point . So in fact, any finite-size transition point belongs to one fixed- and one fixed- scheme, simultaneously. The same holds for any canonical function . Thus, the orthogonal crossing schemes are equivalent and we may translate a functional dependence from one scheme to another by a Taylor series expansion. Then, expanding around some yields

(1) |

which may be solved for . The remaining task is to identify and evaluate suitable functional dependencies.

Using this procedure, we will recapture for a fixed density the finite-size correction to the transition temperature and the scaling of the finite-size rounding (or width of the finite-system transition region) . This is directly related to the linear extension of the droplet , which becomes the relevant length scale for large system sizes.

### ii.1 Finite-size scaling of the droplet condensation-evaporation transition

Biskup et al. Biskup2002 () showed for a -dimensional liquid-vapor system in equilibrium a vanishing probability of intermediate-sized droplets. This was proven rigorously for the lattice-gas interpretation of the 2D Ising model and justifies the restriction of the discussion to two relevant phases: the homogeneous gas phase and the mixed phase of a single macroscopic droplet in equilibrium with a bulk gas surrounding it. At coexistence both phases are equally probable and in general (finite volume) the dominance of either phase is the solution of a saddle point problem (see also Ref. Binder2003 ()). Choosing a fixed temperature , we consider a system at supersaturated density with particle excess , where is the equilibrium density of the infinite-system gas phase at this temperature. In fact, the particle excess is temperature dependent, see Fig. 1. A fraction of this particle excess is considered to be in the single macroscopic droplet, i.e., , while the remaining excess is in the gas phase such that . Up to higher-order corrections from surface and density fluctuations as well as translational entropy of the droplet motion, Biskup et al. Biskup2002 () introduced a rescaled (size-independent) density

(2) |

with the infinite-system liquid and gas densities and , the Wulff-shape surface free-energy per unit volume and the (reduced) isothermal compressibility , see also Ref. Zierenberg2014JPCS (). Here, only is variable and the remaining parameters are constants encoding all system-specific details, like the droplet shape. In the limit of large system sizes, they calculate the fraction of particles in the largest droplet as a “universal” function of . At the transition density the system changes between the vapor phase () and the mixed phase ().

Now consider a fixed density with varying temperature . Then the system-specific constants become functions of the temperature, namely , , and . The temperature-dependent particle excess may be understood by the fact that for decreasing temperature (fixed density) the difference increases and with it the supersaturation, see also Fig. 1. We may rewrite Eq. (2) as

(3) |

where we identify

(4) |

At the condensation-evaporation transition, is constant and the left-hand side of Eq. (3) is depending only on the system size . Then, for a fixed finite system size, a suitable combination of and solves Eq. (3) yielding the finite-size transition point at in Fig. 1. This transition point may be obtained either numerically exact or by a Taylor expansion. Keeping constant, we proceed by expanding in Eq. (3) around the infinite-system transition temperature , where and thus . Solving for the finite-size transition temperature yields

(5) |

In terms of the number of particles this means to first order

(6) |

We need to emphasize that Eq. (3) is only the leading-order result. Numeric tests at fixed temperature already showed apparent higher-order corrections Nussbaumer2006 (); Nussbaumer2008 (); Zierenberg2014JPCS ().

For the Ising model, we can go a little further and evaluate the leading finite-size scaling explicitly. To this end, we relate the magnetization and the density via and identify the spontaneous magnetization as well as, from the symmetry of the Ising model, . In addition, we find the magnetic susceptibility and , due to the shift in the energy scale when exploiting the equivalence of the Ising model and the lattice gas model Zierenberg2014JPCS (). Thus, we obtain

(7) |

with when rewriting the Hamiltonian according to the definition in Sec. III and Ref. Zierenberg2014JPCS ().

In two dimensions, the involved quantities are known analytically or up to arbitrary precision: is described by the Onsager-Yang equation Onsager (), is obtained from sufficiently long series expansions Orrick2001 (); Nickelchi2 (); chi2new () and can be obtained from the volume of the Wulff plot Zia1990 (). For a collection of equations we refer to Ref. Nussbaumer2008 (). We numerically evaluated Eq. (7) in Eq. (3), fixing the density and volume and solving with a bisection algorithm for the corresponding transition temperature. We will compare this result to the numerical finite-size results later in Fig. 6 (top) and refer to it as the full solution of Eq. (3).

In three dimensions, we may make use of low-temperature series expansions of the spontaneous magnetization 3DLowTempM (). This allows us to estimate the infinite-system transition temperature by solving . We will compare these results to the finite-size scaling fits of our numerical data and refer to it as the solution from low-temperature series expansion in Fig. 6 (center).

The analytic results may be used to recapture the scaling of the droplet size at condensation, where the fraction will be in the largest droplet. For the volume of the droplet at condensation it follows Biskup2002 (), where the total particle excess at condensation encodes the shape of the droplet (see Eq. (2)). In general, the volume of an ideal droplet may be expressed by , where is a geometric shape factor that allows to describe both spherical and cubic droplets, where the latter may occur in lattice systems below the roughening transition Zierenberg2014JPCS (). Equating both droplet volumes, inserting from Eq. (2), and solving for the radius yields

(8) |

consistent with results in the literature Binder1980 (); Binder2003 (). For fixed density, this is equivalent to . Moreover, the leading finite-size scaling corrections in Eq. (5) may be expressed in powers of , i.e., .

### ii.2 Rounding of the droplet condensation-evaporation transition

In analogy, one may argue that the rounding of the transition at fixed density should scale with the system size in the same way as the rounding of the transition density (or magnetization) at fixed temperature. The latter was derived by Binder Binder2003 () using a phenomenological theory. He predicted a finite-size rounding proportional to . Starting from a two-state approximation the discussion may be reduced to a homogeneous and an inhomogeneous phase. Then, the rounding of the transition may be related to a characteristic width of order unity, where is the free-energy difference between both phases.

Expanding the free-energy difference around the finite-size transition temperature , recalling that , yields

(9) |

The free-energy difference vanishes at in the limit of large system sizes, considering that both phases contribute with equal probability. In the gas phase, the particles may be considered non-interacting. Thus, the energy difference is dominated by the droplet energy, which depends on the droplet volume . At , Eq. (8) relates the droplet radius to the system volume and the energy difference may be approximated to leading order as . The finite-size corrections from the transition temperature appear merely as corrections to the energy difference, such that Eq. (9) simplifies to in the limit of large system sizes. The condition yields to leading order the rounding width and in terms of the particle number

(10) |

Notice that the radius of the droplet as relevant length scale was used in the argumentation such that for large system sizes the rounding may be identified as .

## Iii Models and Methods

In order to investigate the universal aspects of the condensation-evaporation transition, we employ two different particle gas models, for a sketch see Fig. 2. Bridging the gap to analytic solutions Biskup2002 (), we consider a lattice gas model equivalent to the Ising model at fixed magnetization Pleimling2001 (); Nussbaumer2006 (). We add to our investigation a Lennard-Jones gas model, which — in contrast to the lattice gas — is not symmetric with respect to particle-hole exchange (see also Binder2004 (); Binder2009 ()). The boundary conditions are periodic, which is common for a finite-size scaling analysis. That way, non-intended interactions with the boundaries are avoided. For both models, we consider a finite interaction range smaller than the linear extension of the system, which safely allows the use of periodic boundary conditions.

### iii.1 Lattice gas

In the discrete lattice gas (DLG) model, particles are described by occupied sites on a two-dimensional (2D) square or three-dimensional (3D) cubic lattice of size . The interaction is purely short-range, with attraction between nearest neighbors and hard-core repulsion by the condition that each site can be occupied by only one particle. Introducing a state-variable for each site, , the Hamiltonian may be written as

(11) |

equivalent to the Ising model () at fixed magnetization if . Here, the density can only be adjusted approximately by occupying an integer number of sites .

### iii.2 Lennard-Jones gas

The Lennard-Jones gas (LJG) model is a continuous particle gas model, where in principle all particles interact with each other via the potential

(12) |

with and such that the potential minimum is at . Being consistent with literature and reducing the computational demand, we introduce a cutoff radius above which the potential is zero. The remaining potential is then shifted by in order to be continuous:

(13) |

That way, we may use a domain decomposition, only calculating interactions with particles in neighboring domains. This reduces the computational demand especially in the gas phase. The Hamiltonian is given by

(14) |

Here, we consider only the three-dimensional case of particles in a cubic box of size with periodic boundary conditions and .

For the present study, we considered the Lennard-Jones gas with the same density as for the lattice gas. The Lennard-Jones gas may be reasonably applied to non-polar gases, for example Argon (Ar). For this case, molecular dynamics simulations were matched with experiments already in the 60s Rahman1964 (). Corresponding parameters are and , which serve well for an order-of-magnitude comparison. This would lead to a real temperature . In order to compare to the literature boiling point, the density would have to be adjusted accordingly, e.g., for Argon, Tegeler1999 () with a gas density at atmospheric pressure. The unit length in our system is related to the parametrized length scale as . Assuming a molar weight of for Argon, this yields the conversion for the density . In dimensionless units, the experimental boiling (evaporation) temperature and gas density yield and , respectively.

### iii.3 Parallel multicanonical simulations

The droplet condensation-evaporation transition shows a phase coexistence between a supersaturated gas phase and a mixed phase consisting of a droplet and a remaining gas phase. The resulting free-energy barrier Nussbaumer2010 (); Neuhaus2003 () decreases the probability to pass from one phase to the other in the canonical ensemble which leads to effects like hysteresis and other effects characteristic for first-order like phase transitions. In order to overcome these barriers and to sample the transition point with high accuracy, we apply multicanonical simulations Berg1992 (); JankeMUCA (). This is a natural choice for the fixed-density scheme, while for fixed-temperature a similar approach is not straight forward. The Boltzmann weight , where with , is replaced by an a priori unknown weight function which is iteratively modified in order to yield a flat histogram JankeMUCA (). This allows the sampling of suppressed states between the coexisting phases and produces accurate estimates of the observables when reweighting around the transition temperature. As a result, we obtain a full temperature range from a single simulation applying standard time-series and histogram reweighting techniques.

However, this requires to sample a broad energy range with sufficiently many tunnel events across the transition point from high energy to low energy and vice versa. We make use of a parallel implementation of the multicanonical method Zierenberg2013CPC (), which allows a reduction of computational time. In the cases of lattice gas condensation, it was tested and shown to speed up the simulation by the number of cores used Zierenberg2014JPCS (). We considered lattice sizes up to for the 2D and for the 3D lattice gas as well as Lennard-Jones particle systems with up to particles with a fixed density . The Monte Carlo updates in the lattice systems include local particle shifts to a nearest-neighbor site combined with particle displacements to a random new site. For the Lennard-Jones gas, we restrict ourselves to local random particle displacements.

### iii.4 Metropolis simulations

In addition, we applied standard Metropolis simulations in order to sample the droplet phase at coexistence. We can use the finite-size scaling corrections from the fit to the transition temperatures obtained with the parallel multicanonical simulations. Then, we make use of the first-order nature of the transition causing a suppression of the intermediate states between the vapor and the mixed phase. Preparing a system in the droplet phase, the Metropolis simulation thermalizes and only samples the droplet phase for a finite time. With increasing system size, this finite time gets larger and we may safely sample the size of the critical droplet with a lot less statistics than would be required for the multicanonical simulation. However, we need the multicanonical data of suitable system sizes in order to obtain valid estimates of the transition temperature. We do not use the Metropolis data to draw conclusions about the transition temperature or the finite-size rounding but instead directly measure the critical droplet size. This allows to reach an additional order of magnitude in . An example for the 3D cases is shown in Fig. 3.

### iii.5 Observables

The focus of this study lies on two observables and their thermal derivatives, namely the energy and the fraction of particles in the largest cluster . The latter is determined by identifying the largest cluster as the maximal number of connected particles . In case of the lattice gas, the connected particles are simply nearest neighbors. For the Lennard-Jones gas particles are defined as connected if . Thus, is a measure of the “mass” of the largest droplet and allows one to study the homogeneous nucleation transition efficiently. An example for the 3D Lennard-Jones gas is given in Fig. 4 (left). The consideration of these types of geometric clusters is a safe choice for dilute systems. Note that for rather dense lattice systems also the stochastic Swendsen-Wang cluster definition has been considered Schmitz2013 ().

Anticipating a first-order phase transition, we expect in the thermodynamic limit a discontinuity in the energy as well as in . Hence, the thermal derivatives of the observables will show a pronounced peak at the transition temperature for finite systems, see Fig. 5. The determination of the peak may be done very precisely due to the full temperature range from multicanonical simulations. The error is estimated by using jackknife error analysis Efron1982 (), where we make use of the independent parallel production runs combining all but one time series for each jackknife bin. In the case of the energy, the thermal derivative is related to the specific heat .

Moreover, we will use the mass-fraction in order to directly verify the predictions about the scaling of the critical droplet size that was assumed and used in the theory. We want to focus on the droplet phase and estimate the average fraction of particles in the largest cluster at the transition temperature . To this end, we reweight our multicanonical data to droplet size distributions by time-series reweighting, adding each reweighting factor to the histogram bin corresponding to the observable. An example is shown in Fig. 4 (right). Assuming that both phases are at equal weight, we identify such that . Then

(15) |

is a robust estimator for sufficiently large systems. In fact, we directly deal with the number of particles in a droplet such that the integrals become discrete sums. Errors are again estimated using the jackknife error analysis. Systematic errors may arise from imperfect estimates of the coexistence temperature, the definition of the droplet boundary and, especially for small systems, if the distribution does not show a sufficiently pronounced dip between the two phases.

## Iv Results

The results are presented for a fixed particle density , up to one exemplary case when comparing the Lennard-Jones transition temperature to the Argon boiling point. This choice of a dilute system ensures that the observed transition is really that associated with the formation of a droplet and not of a cylinder or slap MacDowell2006 (). We have to emphasize in the beginning that the lattice model obviously shows discretization effects for small systems because this explicit density may not be realized for any system size and in general remains approximate. We will see, however, that this does not influence our main conclusions.

We consider only a single density, because we want to test the leading-order finite-size corrections and scaling exponents. The density will, however, influence the finite-size scaling limit. In order to illustrate this, consider the low-temperature series expansion of the 3D Ising model. Here, the spontaneous magnetization is given to first order as . For or , this yields for the density and to first order for the inverse temperature . This result may be similarly deduced from microcanonical arguments in continuous systems (see Ref. Zierenberg2014JCP ()).

### iv.1 Finite-size scaling of the droplet condensation-evaporation transition temperature

The transition temperature, where the transition-related observables show a discontinuity, may be estimated by the peak locations of their thermal derivatives, for an example see Fig. 5. In the limit of large systems, this transition temperature is predicted to scale as

(16) |

see Sec. II.1. Figure 6 shows numerical results for the three considered situations, the lattice gas in two and three dimensions and the Lennard-Jones gas in three dimensions. The estimated transition temperatures are obtained from the maxima of the specific heat (red pluses) as well as the thermal derivative of the fraction of particles in the largest cluster (green crosses). Both estimates are remarkably similar as expected for first-order phase transitions, and hence we only show local fits to the estimated transition temperature derived from the specific heat.

Model | Ref | |||||||||
---|---|---|---|---|---|---|---|---|---|---|

Range | Range | Range | ||||||||

2D DLG | [2500: ] | 0.39884(3) | 0.4 | [ 400: ] | 0.3982(1) | 1.8 | [ 324:900] | 0.3903(2) | 8.0 | |

3D DLG | [2160: ] | 0.62341(4) | 0.6 | [1663: ] | 0.6229(4) | 1.9 | [ 68:243] | 0.5840(3) | 1.6 | |

3D LJG | [ 160: ] | 0.7106(4) | 0.8 | [ 10: ] | 0.7011(4) | 1.1 | [ 12: 48] | 0.6597(4) | 1.0 |

In all three cases, we observe for large system sizes a proper finite-size scaling behavior according to the predicted scaling, shown by a good quality of the least-squares fits (a reduced per degree of freedom of about ). In accordance with the literature, we verify that this requires rather large system sizes for the leading order fit (dashed dark blue fit). The intermediate sized systems may be included in an empirical fit with the next higher orders (dotted light blue fit), i.e., . In addition, we want to mention that for intermediate system sizes an effective different scaling behavior may be observed. While not completely apparent, it may be justified by sufficiently good , that the transition temperature can be locally describes by a behavior (dash-dotted orange fit). For details of the individual fits see Table 1 and the following discussion. The scaling of the intermediate regime is consistent with studies of flexible homopolymer aggregation Zierenberg2014JCP (), where a large fraction of the system is involved in the formation of the droplet/aggregate. If almost all constituents are involved in the transition, then the linear extension of the homogeneous, isotropic condensate is just – which justifies the assumed scaling behavior. This leads to two possible conclusions: Either the local scaling function is a polymer property or it is a generic property for a small number of polymers, where the latter seems to be in agreement with the presented results. In order to explore this observation further, we will investigate the size of the critical droplet and the rounding of the transition (which in fact provides suitable fitting ranges for the transition temperature) in the following subsections.

In the case of the 2D lattice gas, we may compare directly to the analytic solution of the infinite system and to the full solution of Eq. (3) for finite systems derived in Sec. II.1. The full solution shows large deviations for small systems as expected. With increasing system size, however, it starts to describe the finite-size scaling approximately. This already gives a hint to the choice of a proper fitting-range with a leading-order scaling behavior. A least-square fit of the leading-order for the largest system sizes yields an adequate and an infinite-size transition temperature . This is consistent with the analytic result inverting Onsager’s solution for the magnetization Onsager (), with ,

(17) |

The resulting is shown in Fig. 6 (top) by the arrow. Including the next order and fitting yields a with , which slightly deviates from the exact solution. This may be taken as a hint that our next-order term is only an effective correction and additional corrections of the same order may be apparent. Interesting to compare is also the amplitude of the leading-order correction . The linear and higher order fits yield and respectively. This may be compared to the power-series expansion of the full solution (see Eq. (5) and Eq. (7)). Making use of the analytic solution for , the series expansion for up to order and the integral solution of (for a list of equations see Ref. Nussbaumer2008 ()), we may numerically differentiate Eq. (7) and calculate in lattice gas units in decent agreement with the leading-order fit. The previously mentioned intermediate scaling regime is not very prominent for the 2D lattice gas where a least-square fit to in the (already small) range still yields a . Moreover, the infinite-size extrapolation is obviously wrong. It is worth noting, that the intermediate scaling regime in polymer aggregation Zierenberg2014JCP () was observed in three dimensions, which suggests that the prominence of this regime may depend on the dimension.

For the 3D lattice gas, we may compare the finite-size scaling results of the infinite-size transition temperature to low-temperature expansions, see Sec. II.1. A fit of the leading-order scaling behavior to sufficiently large systems in the range yields with which is in the vicinity of the (not exact) low-temperature expansion 3DLowTempM (), again shown by the arrow. Increasing the fit range still yields reasonably good fits with higher of the same quality as the fit including the next order , see also Table 1. Notice that in the case of the 3D lattice gas also the small system sizes seem to coincide with the leading order fit. However, considering only an intermediate regime allows to fit the behavior with qualitatively good local agreement. If the larger system sizes were not present, this could be interpreted as the leading order scaling behavior. Comparing the fit to this effective ansatz with the low-temperature series expansion shows, however, strong deviations.

In the case of the 3D Lennard-Jones gas one may best see the arising peculiarities. The leading-order fit for to yields with but shows a clear deviation for small system sizes. Including the next order for yields with , in rough agreement with the leading-order result, and recaptures the deviation of the small system sizes. This is consistent with results for the same Lennard-Jones model Trokhymchuk1999 (); Chapela1977 (). However, considering an intermediate regime with the ansatz yields a qualitatively good fit with and , which deviates strongly from the fit. Again, if the largest system sizes were not present this could be interpreted as the finite-size scaling corrections, especially if no reference temperature is available. Locally it seems that this is an intermediate scaling regime which is, however, already covered by the theoretically predicted scaling behavior including empirical higher-order corrections. This shows the necessity to control that the available data is really in the expected leading-order scaling regime.

Comparing to the boiling temperature of Argon, we have to exemplary consider a density . Then, the leading-order fit yields (not shown here) which differs from the experimental result (see Sec. III.2) by . While the order-of-magnitude is comparable, the difference is expected for the truncated Lennard-Jones potential Trokhymchuk1999 ().

The largest systems considered for the Lennard-Jones gas included 512 particles in a box of length . This is a lot smaller than the system sizes considered in Ref. Binder2009 () () at fixed temperature (for their parameterization) with typical particle numbers . Still, in this fixed- approach they did not see the predicted scaling behavior of the transition density shift but needed to extrapolate an effective exponent (smaller than ) in order to recover the theoretical prediction . Similarly, a direct fit of the leading-order power-low exponent to the present data yields already for smaller system sizes . This implies that an orthogonal phase boundary crossing may lead in certain situations to reduced finite-size corrections and serves as a useful, complementary approach.

### iv.2 Finite-size rounding of the droplet condensation-evaporation transition

In the limit of large systems, the rounding of the transition is predicted to scale as for fixed density, see Sec. II.2. Figure 7 shows the rescaled specific heat from Fig. 5 for the 3D Lennard-Jones gas. At a first-order transition, the maximum of the specific heat should scale inverse proportional to the rounding of the transition Imry1980 (). The rounding of the transition is here estimated as the half-width of the specific-heat peak, defined as the width where . Errors are obtained by jackknife error analysis.

Figure 8 shows the finite-size scaling of the rounding for all three considered systems. In all cases, we could clearly identify two different scaling regimes: For intermediate system sizes one may observe (dashed line) and for large system sizes (solid line). The latter is the expected scaling of the rounding as predicted by theory Binder2003 (); Binder2009 () and reformulated in Sec. II.2. Thus, having only data available in the intermediate regime would suggest a wrong finite-size scaling behavior consistent with the intermediate regime for the transition temperature. However, the intermediate regime is nicely fitted by considering the predicted scaling behavior for large system sizes and empirically including the next higher-order corrections (dash-dotted light blue fit), i.e., .

In the case of the 2D lattice gas, smallest systems show no systematic behavior which explains the large deviations in Fig. 6. For the intermediate regime one can see a direct particle or volume dependence ( fixed). The onset of the large-system regime is consistent with the scaling range in the transition temperature that approximately coincides with the full solution (Fig. 6). A fit including higher order corrections yields a including already system sizes and thus including the intermediate regime.

For the 3D lattice gas and the Lennard-Jones gas, the intermediate regime is apparent for quite small systems and also shows the -proportionality up to the crossover to the large-system regime. Again, the crossover is consistent with a good choice of a leading-order fitting range for the finite-size scaling of the transition temperature. The fit to the rounding of the transition, including higher-order corrections, allows to include system sizes with for the lattice gas and with for the Lennard-Jones gas case. Again, this includes (parts of) the intermediate regime by considering empirical higher-order corrections.

Previous studies of the lattice gas in two and three dimensions at fixed- showed significant deviations from the predicted exponents for the transition rounding Malakis2007 (), using average densities of states from Wang-Landau simulations. For the 3D case they found an effective scaling of the rounding with . The present results on the other hand clearly confirm the “large”-system scaling behavior. Direct fits of the power-law behavior to the largest system sizes of the 3D lattice gas model yield effective exponents and for and respectively, close to the predicted scaling .

The finite-size rounding shows to be a good observable to identify the previously noticed intermediate regime. The width of the transition is associated to the fluctuations of the system Imry1980 (), which should depend on the inverse volume of the relevant system size. The relevant system size remains the droplet, which scales for large systems as . On the other hand, for small systems the droplet includes a large fraction of the system, see the following discussion. Thus the local scaling of the intermediate regime is not surprising. Knowing the scaling behavior for large systems, however, allows to anticipate this gradual crossover with empirical higher-order corrections.

### iv.3 Scaling regimes of the critical droplet size

We identified the leading-order finite-size scaling corrections as powers of the critical droplet radius. This requires that the mass of the critical droplet scales in the limit of large systems as . We want to test this relation by considering the fraction of particles in the largest cluster — or in other words the mass of the droplet divided by the total number of particles — that should scale as .

Figure 9 shows that this assumption seems to be valid for all three considered models. In each case, we show the expectation value in the droplet phase at the transition temperature. The horizontal dashed line at corresponds to all particles being in the droplet. For one, the results from our multicanonical simulation (red pluses) are measured as described in Sec. III.5 with jackknife errors. In addition, there are statistical averages from canonical Metropolis simulations (green crosses) prepared in the droplet phase at as extrapolated in Sec. IV.1. To be precise, we considered the following pairs : for the 3D lattice gas as well as for the 3D Lennard-Jones gas. In this case, the errors are statistical errors including the integrated autocorrelation times. The straight line is the expected scaling behavior shifted as a guide to the eye in order to be compared with the numerical data.

The 2D lattice gas shows the expected scaling behavior quite early, already for roughly 100 particles. In the case of the 3D lattice gas the situation already changes and the expected scaling behavior only starts for quite large system sizes of the order of 2000 particles. This is consistent with the 3D Lennard-Jones gas, where, however, the available system sizes are much smaller. Still the expected scaling may be anticipated. Including the Metropolis data, it shows that the leading-order extrapolation of the available transition temperatures is not precise enough and seems to overestimate the transition temperature for larger systems, consistent with the observations in Sec. IV.1.

All cases show that for small system sizes, the majority of particles is included in the single macroscopic droplet. A corresponding constant behavior (the horizontal dashed line) would explain an intermediate scaling regime where the relevant length scale depends on the system volume. From Fig. 9, this may be at most guessed locally. However, regarding the measurement of the average critical droplet we have to mention that the small systems including only a few particles also show narrow distributions. These are discretized by nature and only defined in the range . Thus, the observable may not provide the required measurement sensitivity. Still, having a large fraction of particles in the droplet, especially compared to the vanishing fraction in the limit of large systems, suffices for the argument that locally an effective intermediate scaling becomes visible. Moreover, it becomes clear that there are corrections to the size of the critical droplet for small and intermediate systems. This should also have an effect on the scaling of the transition temperature and rounding.

## V Conclusions

We presented results on the canonical finite-size scaling of the droplet condensation-evaporation transition. To this end, we performed a comparative study in which we considered two- and three-dimensional lattice gas as well as three-dimensional off-lattice Lennard-Jones gas models in order to draw general conclusions. In all cases, we could verify the theoretical predictions about the finite-size corrections to the transition temperature and the finite-size rounding in the limit of large systems Biskup2002 (); Binder2003 (). In addition, we identified an intermediate regime with an effective different scaling behavior in all three cases, most apparent in the finite-size rounding. This regime is consistent with results from flexible polymers Zierenberg2014JCP (), where the aggregate included most of the polymers. Our measurement of the average droplet size in equilibrium with vapor at coexistence showed that indeed for small system sizes a large fraction of the system is included in the single macroscopic droplet. The scaling of the droplet size, however, is nicely described by the theoretical predictions Binder1980 (); Binder2003 () in the limit of large systems.

Without knowledge of the larger system sizes, the intermediate regime could give way to a wrong finite-size scaling behavior with qualitatively good fit results for both the transition temperature and rounding. However, we were able to show that there is a gradual crossover to the large-system scaling regime which can be approximately described by including empirical higher-order corrections to the theoretical scaling predictions, namely polynomial orders of the critical droplet size. If only the leading-order corrections are considered, the rounding is a good source getting a lower fit limit for a qualitatively good fit. This was compared to analytic and low-temperature results for the lattice cases. It is thus also suitable in order to identify the finite-size scaling regime for phase transitions with mixed phases, like the condensation-evaporation transition but also more complex systems like polymer aggregation. Considering only the local, intermediate regime and applying standard finite-size scaling approaches (with possibly wrong corrections) would then lead to wrong estimates of the infinite-system limit. A similar note of caution has been recently demonstrated in the finite-size scaling of self-avoiding walks on percolation clusters Fricke2014 ().

An intuitive approach to the leading-order finite-size scaling corrections is the competition of volume () and surface () contributions. For a first-order phase transition, this gives rise to a finite-size correction of the order , where is the relevant length scale of the system Privman1990 (); Borgs1995 (); Borgs2002 (); Zierenberg2014JCP (). We argue that the linear extension of the droplet at coexistence plays this dominant role. The condensation-evaporation transition clearly connects a homogeneous and a mixed phase in the canonical ensemble. However, within the canonical approach, we may consider a (virtual) subsystem with the volume of the transition droplet. By translational invariance, this subsystem may be always constructed around the largest droplet. Above the transition temperature, this subsystem includes a homogeneous gas phase while at and below the transition it is filled by the largest droplet and hence shows a homogeneous liquid phase. Thus, this may be interpreted as a grand-canonical transition between homogeneous phases in the virtual system spanned by the volume of the critical droplet. By construction, this virtual volume would have open boundary conditions yielding a finite-size shift of order and a finite-size rounding of order . This picture is consistent with rigorous results for non-periodic first-order phase transitions Borgs2002 (). However, this argument relies on the finite-size scaling of , which was shown to be already non-trivial.

The present study shows that considering an orthogonal crossing of the phase boundary still yields the same finite-size corrections and serves as a complementary tool. On the example of the transition rounding, the fixed-density approach was shown to be closer to the expected large-system scaling behavior already for smaller system sizes. In general, both directions have their advantages and drawbacks both numerically and systematically. This may be exploited for one’s benefit by choosing the suitable direction for the problem at hand, as was also recently demonstrated for the Blume-Capel model Zierenberg2015 ().

###### Acknowledgements.

We thank A. Malakis, M. Mueller and P. Schierz for useful discussions. The project was funded by the European Union and the Free State of Saxony. Computing time provided by the John von Neumann Institute for Computing (NIC) under grant No. HLZ21 on the supercomputer JUROPA at Jülich Supercomputing Centre (JSC) is gratefully acknowledged. Part of this work has been financially supported by the Deutsche Forschungsgemeinschaft (DFG) through the Leipzig Graduate School of Excellence GSC185 “BuildMoNa” and by the Deutsch-Französische Hochschule (DFH-UFA) through the German-French Graduate School under Grant No. CDFA-02-07.## References

- (1) F. Müller-Plathe, ChemPhysChem. 3, 754 (2002).
- (2) J. Zierenberg, M. Mueller, P. Schierz, M. Marenz, and W. Janke, J. Chem. Phys. 141, 114908 (2014).
- (3) M. Mueller, W. Janke, and D. A. Johnston, Phys. Rev. Lett. 112, 200601 (2014).
- (4) M. Biskup, L. Chayes, and R. Kotecký, Europhys. Lett. 60, 21 (2002); Commun. Math. Phys. 242, 137 (2003).
- (5) T. Neuhaus and J. S. Hager, J. Stat. Phys. 113, 47 (2003).
- (6) K. Binder, Physica A 319, 99 (2003).
- (7) K. Binder and M. H. Kalos, J. Stat. Phys. 22, 363 (1980).
- (8) A. Nußbaumer, E. Bittner, T. Neuhaus, and W. Janke, Europhys. Lett. 75, 716 (2006).
- (9) A. Nußbaumer, E. Bittner, and W. Janke, Phys. Rev. E 77, 041109 (2008).
- (10) J. Zierenberg, M. Wiedenmann, and W. Janke, J. Phys.: Conf. Ser. 510, 012017 (2014).
- (11) L. G. MacDowell, P. Virnau, M. Müller, and K. Binder, J. Chem. Phys. 120, 5293 (2004).
- (12) L. G. MacDowell, V. K. Shen, and J. R. Errington, J. Chem. Phys. 125, 034705 (2006).
- (13) M. Schrader, P. Virnau, and K. Binder, Phys. Rev. E 79, 061104 (2009).
- (14) S. S. Martinos, A. Malakis, and I. Hadjiagapiou, Physica A 384, 368 (2007).
- (15) M. Pleimling and A. Hüller, J. Stat. Phys. 104, 971 (2001).
- (16) M. Kastner and M. Pleimling, Phys. Rev. Lett. 102, 240604 (2009).
- (17) L. Onsager, Nuovo Cimento (Suppl.) 6, 261 (1949); C. N. Yang, Phys. Rev. 85, 808 (1952).
- (18) W. P. Orrick, B. G. Nickel, A. J. Guttmann, and J. H. H. Perk, Phys. Rev. Lett. 86, 4120 (2001); J. Stat. Phys. 102, 795 (2001).
- (19) B. Nickel, J. Phys. A: Math. Gen. 32, 3889 (1999); J. Phys. A: Math. Gen. 33, 1693 (2000).
- (20) S. Boukraa, A. J. Guttmann, S. Hassani, I. Jensen, J.-M. Maillard, B. Nickel, and N. Zenine, J. Phys. A: Math. Theor. 41, 455202 (2008); Y. Chan, A. J. Guttmann, B. G. Nickel, and J. H. H. Perk, J. Stat. Phys. 145, 549 (2011).
- (21) K. Leung and R. K. P. Zia, J. Phys. A: Math. Gen. 23, 4593 (1990).
- (22) M. F. Sykes, J. W. Essam, and D. S. Gaunt, J. Math. Phys. 6, 283 (1965); G. Bhanot, M. Creutz, and J. Lacki, Phys. Rev. Lett. 69, 1841 (1992).
- (23) A. Rahman, Phys. Rev. 136, A405 (1964).
- (24) Ch. Tegeler, R. Span, and W. Wagner, J. Phys. Chem. Ref. Data 28, 779 (1999).
- (25) A. Nußbaumer, E. Bittner, and W. Janke, Prog. Theor. Phys. Suppl. 184, 400 (2010).
- (26) B. A. Berg and T. Neuhaus, Phys. Lett. B 267, 249 (1991); Phys. Rev. Lett. 68, 9 (1992).
- (27) W. Janke, Int. J. Mod. Phys. C 03, 1137 (1992); Physica A 254, 164 (1998).
- (28) J. Zierenberg, M. Marenz, and W. Janke, Comput. Phys. Comm. 184, 1155 (2013).
- (29) F. Schmitz, P. Virnau, and K. Binder, Phys. Rev. E 87, 053302 (2013).
- (30) B. Efron, The Jackknife, the Bootstrap and other Resampling Plans (Society for Industrial and Applied Mathematics, Philadelphia, 1982).
- (31) A. Trokhymchuk and J. Alejandre, J. Chem. Phys. 111, 8510 (1999).
- (32) G. A. Chapela, G. Saville, S. M. Thompson, and J. S. Rowlinson, J. Chem. Soc., Faraday Trans. 2 73, 1133 (1977).
- (33) Y. Imry, Phys. Rev. B 21, 2042 (1980).
- (34) N. Fricke and W. Janke, Phys. Rev. Lett. 113, 255701 (2014).
- (35) V. Privman and J. Rudnick, J. Stat. Phys. 60, 551 (1990).
- (36) C. Borgs and R. Kotecký, J. Stat. Phys. 79, 43 (1995).
- (37) C. Borgs, R. Kotecký, and I. Medved, J. Stat. Phys. 109, 67 (2002).
- (38) J. Zierenberg, N. G. Fytas, and W. Janke, Phys. Rev. E 91, 032126 (2015).