Wetting transition in the two-dimensional Blume-Capel model: a Monte Carlo study

Interfacial segregation of interacting vacancies and their role on the wetting critical properties of the Blume-Emery-Griffiths model

We study the wetting critical behavior of the three-state (s=±1,0) Blume-Emery-Griffiths model using numerical simulations. This model provides a suitable scenario for the study of the role of vacancies on the wetting behavior of a thin magnetic film. To this aim we study a system confined between parallel walls with competitive short-range surface magnetic fields (h_{L}=-|h_{1}|). We locate relevant critical curves for different values of the biquadratic interaction and use a thermodynamic integration method to calculate the surface tension as well as the interfacial excess energy and determine the wetting transition. Furthermore, we also calculate the local position of the interface along the film and its fluctuations (capillary waves), which are a measure of the interface width. To characterize the role played by vacancies on the interfacial behavior we evaluate the excess density of vacancies, i.e., the density difference between a system with and without interface. We also show that the temperature dependence of both the local position of the interface and its width can be rationalized in term of a finite-size scaling description, and we propose and successfully test the same scaling behavior for the average position of the center of mass of the vacancies and its fluctuations. This shows that the excess of vacancies can be associated to the presence of the interface that causes the observed segregation. This segregation phenomena is also evidenced by explicitly evaluating the interfacial free energy.

Higher-Order Cumulants Drive Neuronal Activity Patterns, Inducing UP-DOWN States in Neural Populations

A major challenge in neuroscience is to understand the role of the higher-order correlations structure of neuronal populations. The dichotomized Gaussian model (DG) generates spike trains by means of thresholding a multivariate Gaussian random variable. The DG inputs are Gaussian distributed, and thus have no interactions beyond the second order in their inputs; however, they can induce higher-order correlations in the outputs. We propose a combination of analytical and numerical techniques to estimate higher-order, above the second, cumulants of the firing probability distributions. Our findings show that a large amount of pairwise interactions in the inputs can induce the system into two possible regimes, one with low activity (“DOWN state”) and another one with high activity (“UP state”), and the appearance of these states is due to a combination between the third- and fourth-order cumulant. This could be part of a mechanism that would help the neural code to upgrade specific information about the stimuli, motivating us to examine the behavior of the critical fluctuations through the Binder cumulant close to the critical point. We show, using the Binder cumulant, that higher-order correlations in the outputs generate a critical neural system that portrays a second-order phase transition.

Monte Carlo study of the interfacial adsorption of the Blume-Capel model

We investigate the scaling of the interfacial adsorption of the two-dimensional Blume-Capel model using Monte Carlo simulations. In particular, we study the finite-size scaling behavior of the interfacial adsorption of the pure model at both its first- and second-order transition regimes, as well as at the vicinity of the tricritical point. Our analysis benefits from the currently existing quite accurate estimates of the relevant (tri)critical-point locations. In all studied cases, the numerical results verify to a level of high accuracy the expected scenarios derived from analytic free-energy scaling arguments. We also investigate the size dependence of the interfacial adsorption under the presence of quenched bond randomness at the originally first-order transition regime (disorder-induced continuous transition) and the relevant self-averaging properties of the system. For this ex-first-order regime, where strong transient effects are shown to be present, our findings support the scenario of a non-divergent scaling, similar to that found in the original second-order transition regime of the pure model.

Interfacial adsorption in two-dimensional pure and random-bond Potts models

We use Monte Carlo simulations to study the finite-size scaling behavior of the interfacial adsorption of the two-dimensional square-lattice q-states Potts model. We consider the pure and random-bond versions of the Potts model for q=3,4,5,8, and 10, thus probing the interfacial properties at the originally continuous, weak, and strong first-order phase transitions. For the pure systems our results support the early scaling predictions for the size dependence of the interfacial adsorption at both first- and second-order phase transitions. For the disordered systems, the interfacial adsorption at the (disordered induced) continuous transitions is discussed, applying standard scaling arguments and invoking findings for bulk critical properties. The self-averaging properties of the interfacial adsorption are also analyzed by studying the infinite limit-size extrapolation of properly defined signal-to-noise ratios.

Numerical evidence of hyperscaling violation in wetting transitions of the random-bond Ising model in d=2 dimensions

We performed extensive simulations of the random-bond Ising model confined between walls where competitive surface fields act. By properly taking the thermodynamic limit we unambiguously determined wetting transition points of the system, as extrapolation of localization-delocalization transitions of the interface between domains of different orientation driven by the respective fields. The finite-size scaling theory for wetting with short-range fields [E. V. Albano and K. Binder, Phys. Rev. Lett. 109, 036101 (2012)PRLTAO0031-900710.1103/PhysRevLett.109.036101] establishes that the average magnetization of the sample, with critical exponent β, is the proper order parameter for the study of wetting. While the hyperscaling relationship given by γ+2β=ν_{∥}+ν_{⊥} requires β=1/2 (γ=4, ν_{∥}=3, and ν_{⊥}=2), the thermodynamic scaling establishes that Δ_{s}=γ+β, which in contrast requires β=0 (Δ_{s}=4), where γ, ν_{∥}, ν_{⊥}, and Δ_{s} are the critical exponents of the susceptibility, the correlation lengths parallel and perpendicular to the interface, and the gap exponent, respectively. So, we formulate a finite-size scaling theory for wetting without hyperscaling and perform numerical simulations that provide strong evidence of hyperscaling violation (i.e., β=0) and a direct measurement of the susceptibility critical exponent γ/ν_{⊥}=2.0±0.2, in agreement with theoretical results for the strong fluctuation regime of wetting transitions with quenched noise.

Tricritical wetting in the two-dimensional Ising magnet due to the presence of localized non-magnetic impurities

Fixed vacancies (non-magnetic impurities) are placed along the centre of Ising strips in order to study the wetting behaviour in this confined system, by means of numerical simulations analysed with the aid of finite size scaling and thermodynamic integration methods. By considering strips of size L × M (L << M) where short-range competitive surface fields (H(s)) act along the M-direction, we observe localization-delocalization transitions of the interface between magnetic domains of different orientation (driven by the corresponding surface fields), which are the precursors of the wetting transitions that occur in the thermodynamic limit. By placing vacancies or equivalently non-magnetic impurities along the centre of the sample, we found that for low vacancy densities the wetting transitions are of second order, while by increasing the concentration of vacancies the transitions become of first order. Second- and first-order lines meet in tricritical wetting points (H(tric)(SW), T(tric)(W)), where H(tric)(SW) and T(Tric)(W) are the magnitude of the surface field and the temperature, respectively. In the phase diagram, tricritical points shift from the high temperature and weak surface field regime at large vacancy densities to the T --> 0, H(tric)(SW) --> 1 limit for low vacancy densities. By comparing the locations of the tricritical points with those corresponding to the case of mobile impurities, we conclude that in order to observe similar effects, in the latter the required density of impurities is much smaller (e.g. by a factor 3-5). Furthermore, a proper density of non magnetic impurities placed along the centre of a strip can effectively pin rather flat magnetic interfaces for suitable values of the competing surface fields and temperature.

Far-from-equilibrium growth of magnetic thin films with Blume-Capel impurities

We investigate the irreversible growth of (2+1)-dimensional magnetic thin films. The spin variable can adopt three states (s(I)=±1,0), and the system is in contact with a thermal bath of temperature T. The deposition process depends on the change of the configuration energy, which, by analogy to the Blume-Capel Hamiltonian in equilibrium systems, depends on Ising-like couplings between neighboring spins (J) and has a crystal field (D) term that controls the density of nonmagnetic impurities (s(I)=0). Once deposited, particles are not allowed to flip, diffuse, or detach. By means of extensive Monte Carlo simulations, we obtain the phase diagram in the crystal field vs temperature parameter space. We show clear evidence of the existence of a tricritical point located at D(t)/J=1.145(10) and k(B)T(t)/J=0.425(10), which separates a first-order transition curve at lower temperatures from a critical second-order transition curve at higher temperatures, in analogy with the previously studied equilibrium Blume-Capel model. Furthermore, we show that, along the second-order transition curve, the critical behavior of the irreversible growth model can be described by means of the critical exponents of the two-dimensional Ising model under equilibrium conditions. Therefore, our findings provide a link between well-known theoretical equilibrium models and nonequilibrium growth processes that are of great interest for many experimental applications, as well as a paradigmatic topic of study in current statistical physics.

First-order and tricritical wetting transitions in the two-dimensional Ising model caused by interfacial pinning at a defect line

We present a study of the critical behavior of the Blume-Capel model with three spin states (S=±1,0) confined between parallel walls separated by a distance L where competitive surface magnetic fields act. By properly choosing the crystal field (D), which regulates the density of nonmagnetic species (S=0), such that those impurities are excluded from the bulk (where D=-∞) except in the middle of the sample [where D(M)(L/2)≠-∞], we are able to control the presence of a defect line in the middle of the sample and study its influence on the interface between domains of different spin orientations. So essentially we study an Ising model with a defect line but, unlike previous work where defect lines in Ising models were defined via weakened bonds, in the present case the defect line is due to mobile vacancies and hence involves additional entropy. In this way, by drawing phase diagrams, i.e., plots of the wetting critical temperature (T(w)) versus the magnitude of the crystal field at the middle of the sample (D(M)), we observe curves of (first-) second-order wetting transitions for (small) high values of D(M). Theses lines meet in tricritical wetting points, i.e., (T(w)(tc),D(M)(tc)), which also depend on the magnitude of the surface magnetic fields. It is found that second-order wetting transitions satisfy the scaling theory for short-range interactions, while first-order ones do not exhibit hysteresis, provided that small samples are used, since fluctuations wash out hysteretic effects. Since hysteresis is observed in large samples, we performed extensive thermodynamic integrations in order to accurately locate the first-order transition points, and a rather good agreement is found by comparing such results with those obtained just by observing the jump of the order parameter in small samples.

Influence of nonuniform surface magnetic fields in wetting transitions in a confined two-dimensional Ising ferromagnet

Wetting transitions are studied in the two-dimensional Ising ferromagnet confined between walls where competitive surface fields act. In our finite samples of size L×M, the walls are separated by a distance L, M being the length of the sample. The surface fields are taken to be short-range and nonuniform, i.e., of the form H(1),δH(1),H(1),δH(1),..., where the parameter -1≤δ≤1 allows us to control the nonuniformity of the fields. By performing Monte Carlo simulations we found that those competitive surface fields lead to the occurrence of an interface between magnetic domains of different orientation that runs parallel to the walls. In finite samples, such an interface undergoes a localization-delocalization transition, which is the precursor of a true wetting transition that takes place in the thermodynamic limit. By exactly working out the ground state (T=0), we found that besides the standard nonwet and wet phases, a surface antiferromagnetic-like state emerges for δ<-1/3 and large fields (H(1)>3), H(1)(tr)/J=3, δ(tr)=-1/3,T=0, being a triple point where three phases coexist. By means of Monte Carlo simulations it is shown that these features of the phase diagram remain at higher temperatures; e.g., we examined in detail the case T=0.7×T(cb). Furthermore, we also recorded phase diagrams for fixed values of δ, i.e., plots of the critical field at the wetting transition (H(1w)) versus T showing, on the one hand, that the exact results of Abraham [Abraham, Phys. Rev. Lett. 44, 1165 (1980)] for δ=1 are recovered, and on the other hand, that extrapolations to T→0 are consistent with our exact results. Based on our numerical results we conjectured that the exact result for the phase diagram worked out by Abraham can be extended for the case of nonuniform fields. In fact, by considering a nonuniform surface field of some period λ, with λ<<M, e.g., [H(1)(x,λ)>0], one can obtain the effective field H(eff) at a λ coarse-grained level given by H(eff)=1/λ∑(x=1)(λ)H(1)(x,λ). Then we conjectured that the exact solution for the phase diagram is now given by H(eff)/J=F(T), where F(T) is a function of the temperature T that straightforwardly follows from Abraham's solution. The conjecture was exhaustively tested by means of computer simulations. Furthermore, it is found that for δ≠1 the nonwet phase becomes enlarged, at the expense of the wet one, i.e., a phenomenon that we call "surface nonuniformity-induced nonwetting," similar to the already known case of "roughness-induced nonwetting."

Non-mean-field behavior of critical wetting transition for short-range forces

Critical wetting transition for short-range forces in three dimensions (d=3) is reinvestigated by means of Monte Carlo simulation. Using an anisotropic finite size scaling approach, as well as approaches that do not rely on finite size scaling, we show that the critical wetting transition shows clear deviation from mean-field behavior. We estimate that the effective critical exponent ν_{∥}^{eff}=1.76 ± 0.08 for J/kT=0.35 and ν_{∥}^{eff}=1.85 ± 0.07 for J/kT=0.25. These values are in accord with predictions of Parry et al. [Phys. Rev. Lett. 100, 136105 (2008)]. We also point out that the anisotropic finite size scaling approach in d=3 requires additional modification in order to reach full consistency of simulational results.

Finite-size scaling approach for critical wetting: rationalization in terms of a bulk transition with an order parameter exponent equal to zero

Clarification of critical wetting with short-range forces by simulations has been hampered by the lack of accurate methods to locate where the transition occurs. We solve this problem by developing an anisotropic finite-size scaling approach and show that then the wetting transition is a "bulk" critical phenomenon with order parameter exponent equal to zero. For the Ising model in two dimensions, known exact results are straightforwardly reproduced. In three dimensions, it is shown that previous estimates for the location of the transition need revision, but the conclusions about a slow crossover away from mean-field behavior remain unaltered.