Using thermodynamic fluctuation theory, we study the finite-size rounding of anomalies occurring at first-order phase transitions of the corresponding infinite system. Explicit expressions for thermodynamic functions are derived both for "symmetric transitions" (such as the jump of the spontaneous magnetization in the Ising model from
+Msp to $\ensuremath{-}{M}_{\mathrm{sp}}$ as the field changes from
0+ to ${0}^{\ensuremath{-}}$) as well as for asymmetric cases, but restricting attention to (hyper)cubic system shapes. As an explicit example for the usefulness of these considerations in Monte Carlo simulations where it may be a problem to (i) locate a phase transition and (ii) distinguish first-order from second-order transitions, we present numerical results for the two-dimensional nearest-neighbor Ising ferromagnet in a field, both below the critical temperature
Tc and at
Tc. The numerical results are found to be in very good agreement with the phenomenological theory and it is shown that one may extract the magnitudes of jumps occurring at first-order phase transitions in a well-defined and accurate way.
The first-order phase transition in the one-dimensional
q-state Potts model with long-range interactions decaying with distance as
1/r1+σ has been studied by Monte Carlo numerical simulations for
0<σ≤1 and integer values of
q>2. On the basis of finite-size scaling analysis of interface free energy
ΔFL, specific heat and Binder's fourth order cumulant, we obtain the first-order transition which occurs for
σ below a threshold value
σc(q).