content_image/368516.jpg content_image/368514.jpg
of α and β . The bulk densities obtained using the two methods differ by ≤ 1 . 6%. Thus the MC method captures the behavior found using the Langevin equation to good precision. 

<img src='content_image/368516.jpg'>

We perform MC simulations of systems of L = 100, 200 and 500 wells. Far from the phase transition, density profiles depend only weakly on system size, but near the transition there are significant finite-size effects, as illustrated in Fig. 5. In this case, as L increases, the profile tends to a near-constant value except for a sharp increase near the exit. 

To minimize boundary effects we study the bulk density ρ , defined as the mean density over the 10% of sites nearest the center: 

<img src='content_image/368514.jpg'>

