The random field Ising model with an asymmetric trimodal probability distribution

The Ising model in the presence of a random field is investigated within the mean field approximation based on Landau expansion. The random field is drawn from the trimodal probability distribution P(h_i)=pδ(h_i−h₀)+qδ(h_i+h₀)+rδ(h_i), where the probabilities p,q,r take on values within the interval 0,1] consistent with the constraint p+q+r=1 (asymmetric distribution), h_i is the random field variable and h₀ the respective strength. This probability distribution is an extension of the bimodal one allowing for the existence in the lattice of non magnetic particles or vacant sites. The current random field Ising system displays second order phase transitions, which, for some values of p,q and h₀, are followed by first order phase transitions, thus confirming the existence of a tricritical point and in some cases two tricritical points. Also, reentrance can be seen for appropriate ranges of the aforementioned variables. Using the variational principle, we determine the equilibrium equation for magnetization, solve it for both transitions and at the tricritical point in order to determine the magnetization profile with respect to h₀.