No Arabic abstract
We provide a non-trivial test of supersymmetry in the random-field Ising model at five spatial dimensions, by means of extensive zero-temperature numerical simulations. Indeed, supersymmetry relates correlation functions in a D-dimensional disordered system with some other correlation functions in a D-2 clean system. We first show how to check these relationships in a finite-size scaling calculation, and then perform a high-accuracy test. While the supersymmetric predictions are satisfied even to our high-accuracy at D=5, they fail to describe our results at D=4.
We revisit perturbative RG analysis in the replicated Landau-Ginzburg description of the Random Field Ising Model near the upper critical dimension 6. Working in a field basis with manifest vicinity to a weakly-coupled Parisi-Sourlas supersymmetric fixed point (Cardy, 1985), we look for interactions which may destabilize the SUSY RG flow and lead to the loss of dimensional reduction. This problem is reduced to studying the anomalous dimensions of leaders -- lowest dimension parts of $S_n$-invariant perturbations in the Cardy basis. Leader operators are classified as non-susy-writable, susy-writable or susy-null depending on their symmetry. Susy-writable leaders are additionally classified as belonging to superprimary multiplets transforming in particular $textrm{OSp}(d | 2)$ representations. We enumerate all leaders up to 6d dimension $Delta = 12$, and compute their perturbative anomalous dimensions (up to two loops). We thus identify two perturbations (with susy-null and non-susy-writable leaders) becoming relevant below a critical dimension $d_c approx 4.2$ - $4.7$. This supports the scenario that the SUSY fixed point exists for all $3 < d leq 6$, but becomes unstable for $d < d_c$.
The one-parametric Wang-Landau (WL) method is implemented together with an extrapolation scheme to yield approximations of the two-dimensional (exchange-energy, field-energy) density of states (DOS) of the 3D bimodal random-field Ising model (RFIM). The present approach generalizes our earlier WL implementations, by handling the final stage of the WL process as an entropic sampling scheme, appropriate for the recording of the required two-parametric histograms. We test the accuracy of the proposed extrapolation scheme and then apply it to study the size-shift behavior of the phase diagram of the 3D bimodal RFIM. We present a finite-size converging approach and a well-behaved sequence of estimates for the critical disorder strength. Their asymptotic shift-behavior yields the critical disorder strength and the associated correlation lengths exponent, in agreement with previous estimates from ground-state studies of the model.
In the two-dimensional Ising model weak random surface field is predicted to be a marginally irrelevant perturbation at the critical point. We study this question by extensive Monte Carlo simulations for various strength of disorder. The calculated effective (temperature or size dependent) critical exponents fit with the field-theoretical results and can be interpreted in terms of the predicted logarithmic corrections to the pure systems critical behaviour.
We present a complementary estimation of the critical exponent $alpha$ of the specific heat of the 5D random-field Ising model from zero-temperature numerical simulations. Our result $alpha = 0.12(2)$ is consistent with the estimation coming from the modified hyperscaling relation and provides additional evidence in favor of the recently proposed restoration of dimensional reduction in the random-field Ising model at $D = 5$.
Two numerical strategies based on the Wang-Landau and Lee entropic sampling schemes are implemented to investigate the first-order transition features of the 3D bimodal ($pm h$) random-field Ising model at the strong disorder regime. We consider simple cubic lattices with linear sizes in the range $L=4-32$ and simulate the system for two values of the disorder strength: $h=2$ and $h=2.25$. The nature of the transition is elucidated by applying the Lee-Kosterlitz free-energy barrier method. Our results indicate that, despite the strong first-order-like characteristics, the transition remains continuous, in disagreement with the early mean-field theory prediction of a tricritical point at high values of the random-field.