Dimensional crossover in Kardar-Parisi-Zhang growth.

Two-dimensional (2D) Kardar-Parisi-Zhang (KPZ) growth is usually investigated on substrates of lateral sizes L_{x}=L_{y}, so that L_{x} and the correlation length (ξ) are the only relevant lengths determining the scaling behavior. However, in cylindrical geometry, as well as in flat rectangular substrates L_{x}≠L_{y} and, thus, the surfaces can become correlated in a single direction, when ξ∼L_{x}≪L_{y}. From extensive simulations of several KPZ models, we demonstrate that this yields a dimensional crossover in their dynamics, with the roughness scaling as W∼t^{ß_{2D}} for t≪t_{c} and W∼t^{ß_{1D}} for t≫t_{c}, where t_{c}∼L_{x}^{1/z_{2D}}. The height distributions (HDs) also cross over from the 2D flat (cylindrical) HD to the asymptotic Tracy-Widom Gaussian orthogonal ensemble (Gaussian unitary ensemble) distribution. Moreover, 2D to one-dimensional (1D) crossovers are found also in the asymptotic growth velocity and in the steady-state regime of flat systems, where a family of universal HDs exists, interpolating between the 2D and 1D ones as L_{y}/L_{x} increases. Importantly, the crossover scalings are fully determined and indicate a possible way to solve 2D KPZ models.

Universality class of the special adsorption point of two-dimensional lattice polymers.

In recent work [Rodrigues et al., Phys. Rev. E 100, 022121 (2019)10.1103/PhysRevE.100.022121], evidence was found that the surface adsorption transition of interacting self-avoiding trails (ISATs) placed on the square lattice displays a nonuniversal behavior at the special adsorption point (SAP) where the collapsing polymers adsorb. In fact, different surface exponents ϕ^{(s)} and 1/δ^{(s)} were found at the SAP depending on whether the surface orientation is horizontal (HS) or diagonal (DS). Here, we revisit these systems and study other ones, through extensive Monte Carlo simulations, considering much longer trails than previous works. Importantly, we demonstrate that the different exponents observed in the reference above are due to the presence of a surface-attached-globule (SAG) phase in the DS system, which changes the multicritical nature of the SAP and is absent in the HS case. By considering a modified horizontal surface (mHS), on which the trails are forbidden from having two consecutive steps, resembling the DS situation, a stable SAG phase is found in the phase diagram, and both DS and mHS systems present similar 1/δ^{(s)} exponents at the SAP, namely, 1/δ^{(s)}≈0.44, whereas 1/δ^{(s)}≈0.34 in the HS case. Intriguingly, while ϕ^{(s)}≈1/δ^{(s)} is found for the DS and HS scenarios, as expected, in the mHS case ϕ^{(s)} is about 10% smaller than 1/δ^{(s)}. These results strongly indicate that at least two universality classes exist for the SAPs of adsorbing ISATs on the square lattice.

One-point height fluctuations and two-point correlators of (2+1) cylindrical KPZ systems.

While the one-point height distributions (HDs) and two-point covariances of (2+1) Kardar-Parisi-Zhang (KPZ) systems have been investigated in several recent works for flat and spherical geometries, for the cylindrical one the HD was analyzed for few models and nothing is known about the spatial and temporal covariances. Here, we report results for these quantities, obtained from extensive numerical simulations of discrete KPZ models, for three different setups yielding cylindrical growth. Beyond demonstrating the universality of the HD and covariances, our results reveal other interesting features of this geometry. For example, the spatial covariances measured along the longitudinal and azimuthal directions are different, with the former being quite similar to the curve for flat (2+1) KPZ systems, while the latter resembles the Airy_{2} covariance of circular (1+1) KPZ interfaces. We also argue (and present numerical evidence) that, in general, the rescaled temporal covariance A(t/t_{0}) decays asymptotically as A(x)∼x^{-λ[over ¯]} with an exponent λ[over ¯]=ß+d^{*}/z, where d^{*} is the number of interface sides kept fixed during the growth (being d^{*}=1 for the systems analyzed here). Overall, these results complete the picture of the main statistics for the (2+1) KPZ class.

Semianalytical solutions of Ising-like and Potts-like magnetic polymers on the Bethe lattice.

We study magnetic polymers, defined as self-avoiding walks where each monomer i carries a "spin" s_{i} and interacts with its first neighbor monomers, let us say j, via a coupling constant J(s_{i},s_{j}). Ising-like [s_{i}=±1, with J(s_{i},s_{j})=ɛs_{i}s_{j}] and Potts-like [s_{i}=1,...,q, with J(s_{i},s_{j})=ɛ(s_{i})δ(s_{i},s_{j})] models are investigated. Some particular cases of these systems have recently been studied in the continuum and on regular lattices and are related to interesting applications. Here, we solve these models on Bethe lattices of branching number σ, focusing on the ferromagnetic case in zero external magnetic field. In most cases, the phase diagrams present a nonpolymerized (NP) and two polymerized phases: a paramagnetic (PP) and a ferromagnetic (FP) one. However, quite different thermodynamic properties are found depending on q in the Potts-like polymers and on whether one uses the Ising or Potts coupling in the two-state systems. For instance, for the standard Potts model [where ɛ(1)=⋯=ɛ(q)] with q=2, beyond the θ-point (where the critical and discontinuous NP-PP transition lines meet), a second tricritical point connecting a critical and a discontinuous transition line between the PP-FP phases is found in the system. A triple point where the NP-PP, NP-FP, and PP-FP first-order transition lines meet is also present in the phase diagram. For q≥3 the PP-FP transition is always discontinuous, and the scenario with the triple and θ points appears for q≤6. Interestingly, for q≥7, as well as for the Ising-like model the θ point becomes metastable and the critical NP-PP transition line ends at a critical end point (CEP), where it meets the NP-FP and PP-FP coexistence lines. Importantly, these results indicate that when q≤6 the spin ordering transition is preceded by the polymer collapse transition, whereas for q≥7 and in the Ising case these transitions happen together at the CEPs. Some interesting nonstandard Potts models are also studied, such as the lattice version of the model for epigenetic marks in the chromatin introduced by Michieletto et al. [Phys. Rev. X 6, 041047 (2016)2160-330810.1103/PhysRevX.6.041047]. In addition, the solution of the dilute Ising and dilute Potts models on the Bethe lattice are also presented here, since they are important to understand the PP-FP transitions.

Height distributions in interface growth: The role of the averaging process.

Height distributions (HDs) are key quantities to uncover universality and geometry-dependence in evolving interfaces. To quantitatively characterize HDs, one uses adimensional ratios of their first central moments (m_{n}) or cumulants (κ_{n}), especially the skewness S and kurtosis K, whose accurate estimate demands an averaging over all L^{d} points of the height profile at a given time, in translation-invariant interfaces, and over N independent samples. One way of doing this is by calculating m_{n}(t) [or κ_{n}(t)] for each sample and then carrying out an average of them for the N interfaces, with S and K being calculated only at the end. Another approach consists in directly calculating the ratios for each interface and, then, averaging the N values. It turns out, however, that S and K for the growth regime HDs display strong finite-size and -time effects when estimated from these "interface statistics," as already observed in some previous works and clearly shown here, through extensive simulations of several discrete growth models belonging to the EW and KPZ classes on one- and two-dimensional substrates of sizes L=const. and L∼t. Importantly, I demonstrate that with "1-point statistics," i.e., by calculating m_{n}(t) [or κ_{n}(t)] once for all NL^{d} heights together, these corrections become very weak, so that S and K attain values very close to the asymptotic ones already at short times and for small L's. However, I find that this "1-point" (1-pt) approach fails in uncovering the universality of the HDs in the steady-state regime (SSR) of systems whose average height, h[over ¯], is a fluctuating variable. In fact, as demonstrated here, in this regime the 1-pt height evolves as h(t)=h[over ¯](t)+s_{λ}A^{1/2}L^{α}ζ+⋯-where P(ζ) is the underlying SSR HD-and the fluctuations in h[over ¯] yield S_{1-pt}∼t^{-1/2} and K_{1-pt}∼t^{-1}. Nonetheless, by analyzing P(h-h[over ¯]), the cumulants of P(ζ) can be accurately determined. I also show that different, but universal, asymptotic values for S and K (related, so, to different HDs) can be found from the "interface statistics" in the SSR. This reveals the importance of employing the various complementary approaches to reliably determine the universality class of a given system through its different HDs.

Kardar-Parisi-Zhang growth on square domains that enlarge nonlinearly in time.

Fundamental properties of an interface evolving on a domain of size L, such as its height distribution (HD) and two-point covariances, are known to assume universal but different forms depending on whether L is fixed (flat geometry) or expands linearly in time (radial growth). The interesting situation where L varies nonlinearly, however, is far less explored and it has never been tackled for two-dimensional (2D) interfaces. Here, we study discrete Kardar-Parisi-Zhang (KPZ) growth models deposited on square lattice substrates, whose (average) lateral size enlarges as L=L_{0}+ωt^{γ}. Our numerical simulations reveal that the competition between the substrate expansion and the increase of the correlation length parallel to the substrate, ξ≃ct^{1/z}, gives rise to a number of interesting results. For instance, when γ<1/z the interface becomes fully correlated, but its squared roughness, W_{2}, keeps increasing as W_{2}∼t^{2αγ}, as previously observed for one-dimensional (1D) systems. A careful analysis of this scaling, accounting for an intrinsic width on it, allows us to estimate the roughness exponent of the 2D KPZ class as α=0.387(1), which is very accurate and robust, once it was obtained averaging the exponents for different models and growth conditions (i.e., for various γ^{'}s and ω^{'}s). In this correlated regime, the HDs and covariances are consistent with those expected for the steady-state regime of the 2D KPZ class for flat geometry. For γ≈1/z, we find a family of distributions and covariances continuously interpolating between those for the steady-state and the growth regime of radial KPZ interfaces, as the ratio ω/c augments. When γ>1/z the system stays forever in the growth regime and the HDs always converge to the same asymptotic distribution, which is the one for the radial case. The spatial covariances, on the other hand, are (γ,ω)-dependent, showing a trend towards the covariance of a random deposition in enlarging substrates as the expansion rate increases. These results considerably generalize our understanding of the height fluctuations in 2D KPZ systems, revealing a scenario very similar to the one previously found in the 1D case.

Entropy of fully packed rigid rods on generalized Husimi trees: A route to the square-lattice limit.

Although hard rigid rods (k-mers) defined on the square lattice have been widely studied in the literature, their entropy per site, s(k), in the full-packing limit is only known exactly for dimers (k=2) and numerically for trimers (k=3). Here, we investigate this entropy for rods with k≤7, by defining and solving them on Husimi lattices built with diagonal and regular square-lattice clusters of effective lateral size L, where L defines the level of approximation to the square lattice. Due to an L-parity effect, by increasing L we obtain two systematic sequences of values for the entropies s_{L}(k) for each type of cluster, whose extrapolations to L→∞ provide estimates of these entropies for the square lattice. For dimers, our estimates for s(2) differ from the exact result by only 0.03%, while that for s(3) differs from best available estimates by 3%. In this paper, we also obtain a new estimate for s(4). For larger k, we find that the extrapolated results from the Husimi tree calculations do not lie between the lower and upper bounds established in the literature for s(k). In fact, we observe that, to obtain reliable estimates for these entropies, we should deal with levels L that increase with k. However, it is very challenging computationally to advance to solve the problem for large values of L and for large rods. In addition, the exact calculations on the generalized Husimi trees provide strong evidence for the fully packed phase to be disordered for k≥4, in contrast to the results for the Bethe lattice wherein it is nematic.

Kardar-Parisi-Zhang universality class in (d+1)-dimensions.

The determination of the exact exponents of the KPZ class in any substrate dimension d is one of the most important open issues in Statistical Physics. Based on the behavior of the dimensional variation of some exact exponent differences for other growth equations, I find here that the KPZ growth exponents (related to the temporal scaling of the fluctuations) are given by ß_{d}=7/8d+13. These exponents present an excellent agreement with the most accurate estimates for them in the literature. Moreover, they are confirmed here through extensive Monte Carlo simulations of discrete growth models and real-space renormalization group (RG) calculations for directed polymers in random media (DPRM), up to d=15. The left-tail exponents of the probability density functions for the DPRM energy provide another striking verification of the analytical result above.

Husimi-lattice solutions and the coherent-anomaly-method analysis for hard-square lattice gases.

Although lattice gases composed of particles preventing up to their kth nearest neighbors from being occupied (the kNN models) have been widely investigated in the literature, the location and the universality class of the fluid-columnar transition in the 2NN model on the square lattice are still a topic of debate. Here, we present grand-canonical solutions of this model on Husimi lattices built with diagonal square lattices, with 2L(L+1) sites, for L⩽7. The systematic sequence of mean-field solutions confirms the existence of a continuous transition in this system, and extrapolations of the critical chemical potential µ_{2,c}(L) and particle density ρ_{2,c}(L) to L→∞ yield estimates of these quantities in close agreement with previous results for the 2NN model on the square lattice. To confirm the reliability of this approach, we employ it also for the 1NN model, where very accurate estimates for the critical parameters µ_{1,c} and ρ_{1,c}-for the fluid-solid transition in this model on the square lattice-are found from extrapolations of data for L⩽6. The nonclassical critical exponents for these transitions are investigated through the coherent anomaly method (CAM), which in the 1NN case yields ß and ν differing by at most 6% from the expected Ising exponents. For the 2NN model, the CAM analysis is somewhat inconclusive, because the exponents sensibly depend on the value of µ_{2,c} used to calculate them. Notwithstanding, our results suggest that ß and ν are considerably larger than the Ashkin-Teller exponents reported in numerical studies of the 2NN system.

Fluid-fluid demixing and density anomaly in a ternary mixture of hard spheres.

We report the grand-canonical solution of a ternary mixture of discrete hard spheres defined on a Husimi lattice built with cubes, which provides a mean-field approximation for this system on the cubic lattice. The mixture is composed of pointlike particles (0NN) and particles which exclude up to their first (1NN) and second neighbors (2NN), with activities z_{0}, z_{1}, and z_{2}, respectively. Our solution reveals a very rich thermodynamic behavior, with two solid phases associated with the ordering of 1NN (S1) or 2NN particles (S2), and two fluid phases, one being regular (RF) and the other characterized by a dominance of 0NN particles (F0 phase). However, in most part of the phase diagram these fluid (F) phases are indistinguishable. Discontinuous transitions are observed between all the four phases, yielding several coexistence surfaces in the system, among which a fluid-fluid and a solid-solid demixing surface. The former one is limited by a line of critical points and a line of triple points (where the phases RF-F0-S2 coexist), both meeting at a special point, after which the fluid-fluid coexistence becomes metastable. Another line of triple points is found, connecting the F-S1, F-S2, and S1-S2 coexistence surfaces. A critical F-S1 surface is also observed meeting the F-S1 coexistence one at a line of tricritical points. Furthermore, a thermodynamic anomaly characterized by minima in isobaric curves of the total density of particles is found, yielding three surfaces of minimal density in the activity space, depending on which activity is kept fixed during its calculation.

Order-disorder transition in a two-dimensional associating lattice gas.

We study an associating lattice gas (ALG) using Monte Carlo simulation on the triangular lattice and semianalytical solutions on Husimi lattices. In this model, the molecules have an orientational degree of freedom and the interactions depend on the relative orientations of nearest-neighbor molecules, mimicking the formation of hydrogen bonds. We focus on the transition between the high-density liquid (HDL) phase and the isotropic phase in the limit of full occupancy, corresponding to chemical potential µâ†’∞, which has not yet been studied systematically. Simulations yield a continuous phase transition at τ_{c}=k_{B}T_{c}/γ=0.4763(1) (where -γ is the bond energy) between the low-temperature HDL phase, with a nonvanishing mean orientation of the molecules, and the high-temperature isotropic phase. Results for critical exponents and the Binder cumulant indicate that the transition belongs to the three-state Potts model universality class, even though the ALG Hamiltonian does not have the full permutation symmetry of the Potts model. In contrast with simulation, the Husimi lattice analyses furnish a discontinuous phase transition, characterized by a discontinuity of the nematic order parameter. The transition temperatures (τ_{c}=0.51403 and 0.51207 for trees built with triangles and hexagons, respectively) are slightly higher than that found via simulation. Since the Husimi lattice studies show that the ALG phase diagram features a discontinuous isotropic-HDL line for finite µ, three possible scenarios arise for the triangular lattice. The first is that in the limit µâ†’∞ the first-order line ends in a critical point; the second is a change in the nature of the transition at some finite chemical potential; the third is that the entire line is one of continuous phase transitions. Results from other ALG models and the fact that mean-field approximations show a discontinuous phase transition for the three-state Potts model (known to possess a continuous transition) lends some weight to the third alternative.

Thermodynamic behavior of binary mixtures of hard spheres: Semianalytical solutions on a Husimi lattice built with cubes.

We study binary mixtures of hard particles, which exclude up to their kth nearest neighbors (kNN) on the simple cubic lattice and have activities z_{k}. In the first model analyzed, point particles (0NN) are mixed with 1NN ones. The grand-canonical solution of this model on a Husimi lattice built with cubes unveils a phase diagram with a fluid and a solid phase separated by a continuous and a discontinuous transition line which meet at a tricritical point. A density anomaly, characterized by minima in isobaric curves of the total density of particles against z_{0} (or z_{1}), is also observed in this system. Overall, this scenario is identical to the one previously found for this model when defined on the square lattice. The second model investigated consists of the mixture of 1NN particles with 2NN ones. In this case, a very rich phase behavior is found in its Husimi lattice solution, with two solid phases-one associated with the ordering of 1NN particles (S1) and the other with the ordering of 2NN ones (S2)-beyond the fluid (F) phase. While the transitions between F-S2 and S1-S2 phases are always discontinuous, the F-S1 transition is continuous (discontinuous) for small (large) z_{2}. The critical and coexistence F-S1 lines meet at a tricritical point. Moreover, the coexistence F-S1,F-S2, and S1-S2 lines meet at a triple point. Density anomalies are absent in this case.

Three stable phases and thermodynamic anomaly in a binary mixture of hard particles.

While the realistic modeling of the thermodynamic behavior of fluids usually demands elaborated atomistic models, much has been learned from simplified ones. Here, we investigate a model where pointlike particles (with activity z0) are mixed with molecules that exclude their first and second neighbors (i.e., cubes of lateral size λ=3a, with activity z2), both placed on the sites of a simple cubic lattice with parameter a. Only hard-core interactions exist among the particles so that the model is athermal. Despite its simplicity, the grand-canonical solution of this model on a Husimi lattice built with cubes revels a fluid-fluid demixing, yielding a phase diagram with two fluid phases (one of them dominated by small particles-F0) and a solidlike phase coexisting at a triple-point. Moreover, the fluid-fluid coexistence line ends at a critical point. An anomaly in the total density (ρT) of particles is also found, which is hallmarked by minima in the isobaric curves of ρT vs z0 (or z2). Interestingly, the line of minimum density crosses the phase diagram starting inside the region where both fluid phases are stable, passing through the F0 one and ending deep inside its metastable region, in a point where the spinodals of both fluid phases cross each other.

Adsorption of two-dimensional polymers with two- and three-body self-interactions.

Using extensive Monte Carlo simulations, we investigate the surface adsorption of self-avoiding trails on the triangular lattice with two- and three-body on-site monomer-monomer interactions. In the parameter space of two-body, three-body, and surface interaction strengths, the phase diagram displays four phases: swollen (coil), globule, crystal, and adsorbed. For small values of the surface interaction, we confirm the presence of swollen, globule, and crystal bulk phases. For sufficiently large values of the surface interaction, the system is in an adsorbed state, and the adsorption transition can be continuous or discontinuous, depending on the bulk phase. As such, the phase diagram contains a rich phase structure with transition surfaces that meet in multicritical lines joining in a single special multicritical point. The adsorbed phase displays two distinct regions with different characteristics, dominated by either single- or double-layer adsorbed ground states. Interestingly, we find that there is no finite-temperature phase transition between these two regions though rather a smooth crossover.

Initial pseudo-steady state & asymptotic KPZ universality in semiconductor on polymer deposition.

The Kardar-Parisi-Zhang (KPZ) class is a paradigmatic example of universality in nonequilibrium phenomena, but clear experimental evidences of asymptotic 2D-KPZ statistics are still very rare, and far less understanding stems from its short-time behavior. We tackle such issues by analyzing surface fluctuations of CdTe films deposited on polymeric substrates, based on a huge spatio-temporal surface sampling acquired through atomic force microscopy. A pseudo-steady state (where average surface roughness and spatial correlations stay constant in time) is observed at initial times, persisting up to deposition of ~104 monolayers. This state results from a fine balance between roughening and smoothening, as supported by a phenomenological growth model. KPZ statistics arises at long times, thoroughly verified by universal exponents, spatial covariance and several distributions. Recent theoretical generalizations of the Family-Vicsek scaling and the emergence of log-normal distributions during interface growth are experimentally confirmed. These results confirm that high vacuum vapor deposition of CdTe constitutes a genuine 2D-KPZ system, and expand our knowledge about possible substrate-induced short-time behaviors.

Grand-canonical solution of semiflexible self-avoiding trails on the Bethe lattice.

We consider a model of semiflexible interacting self-avoiding trails (sISATs) on a lattice, where the walks are constrained to visit each lattice edge at most once. Such models have been studied as an alternative to the self-attracting self-avoiding walks (SASAWs) to investigate the collapse transition of polymers, with the attractive interactions being on site as opposed to nearest-neighbor interactions in SASAWs. The grand-canonical version of the sISAT model is solved on a four-coordinated Bethe lattice, and four phases appear: non-polymerized (NP), regular polymerized (P), dense polymerized (DP), and anisotropic nematic (AN), the last one present in the phase diagram only for sufficiently stiff chains. The last two phases are dense, in the sense that all lattice sites are visited once in the AN phase and twice in the DP phase. In general, critical NP-P and DP-P transition surfaces meet with a NP-DP coexistence surface at a line of bicritical points. The region in which the AN phase is stable is limited by a discontinuous critical transition to the P phase, and we study this somewhat unusual transition in some detail. In the limit of rods, where the chains are totally rigid, the P phase is absent and the three coexistence lines (NP-AN, AN-DP, and NP-DP) meet at a triple point, which is the endpoint of the bicritical line.

Collapse transition in polymer models with multiple monomers per site and multiple bonds per edge.

We present results from extensive Monte Carlo simulations of polymer models where each lattice site can be visited by up to K monomers and no restriction is imposed on the number of bonds on each lattice edge. These multiple monomer per site (MMS) models are investigated on the square and cubic lattices, for K=2 and 3, by associating Boltzmann weights ω_{0}=1, ω_{1}=e^{ß_{1}}, and ω_{2}=e^{ß_{2}} to sites visited by 1, 2, and 3 monomers, respectively. Two versions of the MMS models are considered for which immediate reversals of the walks are allowed (RA) or forbidden (RF). In contrast to previous simulations of these models, we find the same thermodynamic behavior for both RA and RF versions. In three dimensions, the phase diagrams, in space ß_{2}×ß_{1}, are featured by coil and globule phases separated by a line of Θ points, as thoroughly demonstrated by the metric ν_{t}, crossover ϕ_{t}, and entropic γ_{t} exponents. The existence of the Θ lines is also confirmed by the second virial coefficient. This shows that no discontinuous collapse transition exists in these models, in contrast to previous claims based on a weak bimodality observed in some distributions, which indeed exists in a narrow region very close to the Θ line when ß_{1}<0. Interestingly, in two dimensions, only a crossover is found between the coil and globule phases.

Nature of the collapse transition in interacting self-avoiding trails.

We study the interacting self-avoiding trail (ISAT) model on a Bethe lattice of general coordination q and on a Husimi lattice built with squares and coordination q=4. The exact grand-canonical solutions of the model are obtained, considering that up to K monomers can be placed on a site and associating a weight ω_{i} with an i-fold visited site. Very rich phase diagrams are found with nonpolymerized, regular polymerized, and dense polymerized phases separated by lines (or surfaces) of continuous and discontinuous transitions. For a Bethe lattice with q=4 and K=2, the collapse transition is identified with a bicritical point and the collapsed phase is associated with the dense polymerized (solidlike) phase instead of the regular polymerized (liquidlike) phase. A similar result is found for the Husimi lattice, which may explain the difference between the collapse transition for ISATs and for interacting self-avoiding walks on the square lattice. For q=6 and K=3 (studied on the Bethe lattice only), a more complex phase diagram is found, with two critical planes and two coexistence surfaces, separated by two tricritical and two critical end-point lines meeting at a multicritical point. The mapping of the phase diagrams in the canonical ensemble is discussed and compared with simulational results for regular lattices.

Point island models for nucleation and growth of supported nanoclusters during surface deposition.

Point island models (PIMs) are presented for the formation of supported nanoclusters (or islands) during deposition on flat crystalline substrates at lower submonolayer coverages. These models treat islands as occupying a single adsorption site, although carrying a label to track their size (i.e., they suppress island structure). However, they are particularly effective in describing the island size and spatial distributions. In fact, these PIMs provide fundamental insight into the key features for homogeneous nucleation and growth processes on surfaces. PIMs are also versatile being readily adapted to treat both diffusion-limited and attachment-limited growth and also a variety of other nucleation processes with modified mechanisms. Their behavior is readily and precisely assessed by kinetic Monte Carlo simulation.

Transfer-matrix study of a hard-square lattice gas with two kinds of particles and density anomaly.

Using transfer matrix and finite-size scaling methods, we study the thermodynamic behavior of a lattice gas with two kinds of particles on the square lattice. Only excluded volume interactions are considered, so that the model is athermal. Large particles exclude the site they occupy and its four first neighbors, while small particles exclude only their site. Two thermodynamic phases are found: a disordered phase where large particles occupy both sublattices with the same probability and an ordered phase where one of the two sublattices is preferentially occupied by them. The transition between these phases is continuous at small concentrations of the small particles and discontinuous at larger concentrations, both transitions are separated by a tricritical point. Estimates of the central charge suggest that the critical line is in the Ising universality class, while the tricritical point has tricritical Ising (Blume-Emery-Griffiths) exponents. The isobaric curves of the total density as functions of the fugacity of small or large particles display a minimum in the disordered phase.