BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

T. J. Oliveira

T. J. Oliveira contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechcond-mat.softcond-mat.mes-hallcond-mat.mtrl-sciphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging
24works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2020arXiv

Height fluctuations in homoepitaxial thin film growth: A numerical study

We report on the investigation of height distributions (HDs) and spatial covariances of two-dimensional surfaces obtained from extensive numerical simulations of the celebrated Clarke-Vvedensky (CV) model for homoepitaxial thin film growth. In this model, the effect of temperature, deposition flux, and strengths of atom-atom interactions are encoded in two parameters: the diffusion to deposition ratio $R=D/F$ and $\varepsilon$, which is related to the probability of an adatom "breaking" a lateral bond. We demonstrate that the HDs present a strong dependence on both $R$ and $\varepsilon$, and even after the deposition of $10^5$ monolayers (MLs) they are still far from the asymptotics in some cases. For instance, the temporal evolution of the HDs' skewness (kurtosis) displays a pronounced minimum (maximum), for small $R$ and $\varepsilon$, and only at long times it passes to increase (decrease) toward its asymptotic value. However, it is hard to determine whether they converge to a single value or different nonuniversal ones. For large $R$ and/or $\varepsilon$, on the other hand, these quantities clearly converge to the values expected for the Villain-Lai-Das Sarma (VLDS) universality class. A similar behavior is observed in the spatial covariances, but with weaker finite-time effects, so that rescaled curves of them collapse quite well with the one for the VLDS class at long times. Simulations of a model with limited mobility of particles, which captures some essential features of the CV model in the limit of irreversible aggregation ($\varepsilon=0$), reveal a similar scenario. Overall, these results point out that the study of fluctuations in homoepitaxial thin films' surfaces can be a very difficult task and shall be performed very carefully, once typical experimental films have $\lesssim 10^4$ MLs, so that their HDs and covariances can be in the realm of transient regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Permeability and kinetic coefficients for mesoscale BCF surface step dynamics: discrete 2D deposition-diffusion equation analysis

A discrete version of deposition-diffusion equations appropriate for description of step flow on a vicinal surface is analyzed for a two-dimensional grid of adsorption sites representing the stepped surface and explicitly incorporating kinks along the step edges. Model energetics and kinetics appropriately account for binding of adatoms at steps and kinks, distinct terrace and edge diffusion rates, and possibly asymmetric barriers for attachment to steps. Analysis of adatom attachment fluxes as well as limiting values of adatom densities at step edges for non-uniform deposition scenarios allows determination of both permeability and kinetic coefficients. Behavior of these quantities is assessed as a function of key system parameters including kink density, step attachment barriers, and the step edge diffusion rate.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Universality and dependence on initial conditions in the class of the nonlinear molecular beam epitaxy equation

We report extensive numerical simulations of growth models belonging to the nonlinear molecular beam epitaxy (nMBE) class, on flat (fixed-size) and expanding substrates (ES). In both $d=1+1$ and $2+1$, we find that growth regime height distributions (HDs), and spatial and temporal covariances are universal, but are dependent on the initial conditions, while the critical exponents are the same for flat and ES systems. Thus, the nMBE class does split into subclasses, as does the Kardar-Parisi-Zhang (KPZ) class. Applying the "KPZ ansatz" to nMBE models, we estimate the cumulants of the $1+1$ HDs. Spatial covariance for the flat subclass is hallmarked by a minimum, which is not present in the ES one. Temporal correlations are shown to decay following well-known conjectures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Substrate effects and diffusion dominated roughening in Cu2O electrodeposition

Cuprous oxide (Cu2O) films from 25 nm to 1500 nm were electrodeposited on n-Si(100) and Ni/n-Si(100) substrates from aqueous solution at room temperature. X-ray diffraction and transmission electron microscopy imaging show that the Cu2O structure and morphology is strongly affected by the substrate choice, with V shape and U shape columnar growth on n-Si(100) and Ni/n-Si(100), respectively. Atomic force microscopy reveals the presence of rounded grains at the surface in both cases. Anomalous and normal roughening are observed in films grown on n-Si and Ni, respectively, but estimates of scaling exponents are not conclusive. On the other hand, the distributions of local heights, roughness, and extremal heights show good agreement with those of the fourth order linear stochastic equation of Mullins and Herring (MH). Thus, surface dynamics in both systems is dominated by diffusion of adsorbed molecules, with no large scale effect of possible inhomogeneities in mass flux from the solution or in reaction and adsorption rates. In growth on n-Si substrates, the noise amplitude of the MH equation increases in time as t^{0.8}, while the coefficient of the curvature-related term is time-independent. Step edge energy barriers restrict the mass flux across grain boundaries, thus a broad size distribution of initial grains leads to coarsening of the larger ones. This explains their V shape in the thickest films and establishes a connection with the anomalous roughening. These effects are reduced in films grown on Ni/n-Si, which initially have much larger grains with narrower size distributions and, consequently, smaller fluctuations in coarse grained growth rates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Temperature effect on (2+1) experimental Kardar-Parisi-Zhang growth

We report on the effect of substrate temperature (T) on both local structure and long-wavelength fluctuations of polycrystalline CdTe thin films deposited on Si(001). A strong T-dependent mound evolution is observed and explained in terms of the energy barrier to inter-grain diffusion at grain boundaries, as corroborated by Monte Carlo simulations. This leads to transitions from uncorrelated growth to a crossover from random-to-correlated growth and transient anomalous scaling as T increases. Due to these finite-time effects, we were not able to determine the universality class of the system through the critical exponents. Nevertheless, we demonstrate that this can be circumvented by analyzing height, roughness and maximal height distributions, which allow us to prove that CdTe grows asymptotically according to the Kardar-Parisi-Zhang (KPZ) equation in a broad range of T. More important, one finds positive (negative) velocity excess in the growth at low (high) T, indicating that it is possible to control the KPZ non-linearity by adjusting the temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Width and extremal height distributions of fluctuating interfaces with window boundary conditions

We present a detailed study of squared local roughness (SLRDs) and local extremal height distributions (LEHDs), calculated in windows of lateral size $l$, for interfaces in several universality classes, in substrate dimensions $d_s = 1$ and $d_s = 2$. We show that their cumulants follow a Family-Vicsek type scaling, and, at early times, when $ξ\ll l$ ($ξ$ is the correlation length), the rescaled SLRDs are given by log-normal distributions, with their $n$th cumulant scaling as $(ξ/l)^{(n-1)d_s}$. This give rise to an interesting temporal scaling for such cumulants $\left\langle w_n \right\rangle_c \sim t^{γ_n}$, with $γ_n = 2 n β+ {(n-1)d_s}/{z} = \left[ 2 n + {(n-1)d_s}/α \right] β$. This scaling is analytically proved for the Edwards-Wilkinson (EW) and Random Deposition interfaces, and numerically confirmed for other classes. In general, it is featured by small corrections and, thus, it yields exponents $γ_n$'s (and, consequently, $α$, $β$ and $z$) in nice agreement with their respective universality class. Thus, it is an useful framework for numerical and experimental investigations, where it is, usually, hard to estimate the dynamic $z$ and mainly the (global) roughness $α$ exponents. The stationary (for $ξ\gg l$) SLRDs and LEHDs of Kardar-Parisi-Zhang (KPZ) class are also investigated and, for some models, strong finite-size corrections are found. However, we demonstrate that good evidences of their universality can be obtained through successive extrapolations of their cumulant ratios for long times and large $l$'s. We also show that SLRDs and LEHDs are the same for flat and curved KPZ interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interface fluctuations for deposition on enlarging flat substrates

We investigate solid-on-solid models that belong to the Kardar-Parisi-Zhang (KPZ) universality class on substrates that expand laterally at a constant rate by duplication of columns. Despite the null global curvature, we show that all investigated models have asymptotic height distributions and spatial covariances in agreement with those expected for the KPZ subclass for curved surfaces. In $1+1$ dimensions, the height distribution and covariance are given by the GUE Tracy-Widom distribution and the Airy$_2$ process, instead of the GOE and Airy$_1$ foreseen for flat interfaces. These results imply that, when the KPZ class splits into the curved and flat subclasses, as conventionally considered, the expanding substrate may play a role equivalent to, or perhaps more important than the global curvature. Moreover, the translational invariance of the interfaces evolving on growing domains allowed us to accurately determine, in $2+1$ dimensions, the analogue of the GUE Tracy-Widom distribution for height distribution and that of the Airy$_2$ process for spatial covariance. Temporal covariance is also calculated and shown to be universal in each dimension and in each of the two subclasses. A logarithmic correction associated to the duplication of column is observed and theoretically elucidated. Finally, crossover between regimes with fixed-size and enlarging substrates is also investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating the initial growth of a deposit from colloidal suspensions

We study the short time properties of a two-dimensional film growth model in which incident particles execute advective-diffusive motion with a vertical step followed by $D$ horizontal steps. The model represents some features of the deposition of anisotropic colloidal particles of the experiment in Phys. Rev. Lett. {\bf 110}, 035501 (2013), in which wandering particles are attracted to particle-rich regions in the deposit. Height profiles changing from rough to columnar structure are observed as $D$ increases from $0$ (ballistic deposition) to $8$, with striking similarity to the experimental ones. The effective growth exponents matches the experimental estimates and the scaling of those exponents on $D$ show a remarkable effect of the range of the particle-deposit interaction. The nearly ellipsoidal shape of colloidal particles is represented for the calculation of roughness exponents in conditions that parallel the experimental ones, giving a range of estimates that also include the experimental values. The effective dynamic exponents calculated from the autocorrelation function are shown to be suitable to decide between a true dynamic scaling or transient behavior, particularly because the latter leads to deviations in an exponent relation. These results are consistent with arguments on short time unstable (columnar) growth of Phys. Rev. Lett {\bf 111}, 209601 (2013), indicating that critical quenched KPZ dynamics does not explain that colloidal particle deposition problem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scaling in reversible submonolayer deposition

The scaling of island and monomer density, capture zone distributions (CZDs), and island size distributions (ISDs) in reversible submonolayer growth was studied using the Clarke-Vvedensky model. An approach based on rate-equation results for irreversible aggregation (IA) models is extended to predict several scaling regimes in square and triangular lattices, in agreement with simulation results. Consistently with previous works, a regime I with fractal islands is observed at low temperatures, corresponding to IA with critical island size i=1, and a crossover to a second regime appears as the temperature is increased to εR^{2/3} ~ 1, where εis the single bond detachment probability and R is the diffusion-to-deposition ratio. In the square (triangular) lattice, a regime with scaling similar to IA with i=3 (i=2) is observed after that crossover. In the triangular lattice, a subsequent crossover to an IA regime with i=3 is observed, which is explained by the recurrence properties of random walks in two dimensional lattices, which is beyond the mean-field approaches. At high temperatures, a crossover to a fully reversible regime is observed, characterized by a large density of small islands, a small density of very large islands, and total island and monomer densities increasing with temperature, in contrast to IA models. CZDs and ISDs with Gaussian right tails appear in all regimes for R ~ 10^7 or larger, including the fully reversible regime, where the CZDs are bimodal. This shows that the Pimpinelli-Einstein (PE) approach for IA explains the main mechanisms for the large islands to compete for free adatom aggregation in the reversible model, and may be the reason for its successful application to a variety of materials and growth conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Surface and bulk properties of ballistic deposition models with bond breaking

We introduce a new class of growth models, with a surface restructuring mechanism in which impinging particles may dislodge suspended particles, previously aggregated on the same column in the deposit. The flux of these particles is controlled through a probability $p$. These systems present a crossover, for small values of $p$, from random to correlated (KPZ) growth of surface roughness, which is studied through scaling arguments and Monte Carlo simulations on one- and two-dimensional substrates. We show that the crossover characteristic time $t_{\times}$ scales with $p$ according to $t_{\times}\sim p^{-y}$ with $y=(n+1)$ and that the interface width at saturation $W_{sat}$ scales as $W_{sat}\sim p^{-δ}$ with $δ= (n+1)/2$, where $n$ is either the maximal number of broken bonds or of dislodged suspended particles. This result shows that the sets of exponents $y=1$ and $δ=1/2$ or $y=2$ and $δ=1$ found in all previous works focusing on systems with this same type of crossover are not universal. Using scaling arguments, we show that the bulk porosity $P$ of the deposits scales as $P\sim p^{y-δ}$ for small values of $p$. This general scaling relation is confirmed by our numerical simulations and explains previous results present in literature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Universal fluctuations in the growth of semiconductor thin films

Scaling of surface fluctuations of polycrystalline CdTe/Si(100) films grown by hot wall epitaxy are studied. The growth exponent of surface roughness and the dynamic exponent of the auto-correlation function in the mound growth regime agree with the values of the Kardar-Parisi-Zhang (KPZ) class. The scaled distributions of heights, local roughness, and extremal heights show remarkable collapse with those of the KPZ class, giving the first experimental observation of KPZ distributions in $2+1$ dimensions. Deviations from KPZ values in the long-time estimates of dynamic and roughness exponents are explained by spurious effects of multi-peaked coalescing mounds and by effects of grain shapes. Thus, this scheme for investigating universality classes of growing films advances over the simple comparison of scaling exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Crossover in the scaling of island size and capture zone distributions

Simulations of irreversible growth of extended (fractal and square) islands with critical island sizes i=1 and 2 are performed in broad ranges of coverage θand diffusion-to-deposition ratios R in order to investigate scaling of island size and capture zone area distributions (ISD, CZD). Large θand small R lead to a crossover from the CZD predicted by the theory of Pimpinelli and Einstein (PE), with Gaussian right tail, to CZD with simple exponential decays. The corresponding ISD also cross over from Gaussian or faster decays to simple exponential ones. For fractal islands, these features are explained by changes in the island growth kinetics, from a competition for capture of diffusing adatoms (PE scaling) to aggregation of adatoms with effectively irrelevant diffusion, which is characteristic of random sequential adsorption (RSA) without surface diffusion. This interpretation is confirmed by studying the crossover with similar CZ areas (of order 100 sites) in a model with freezing of diffusing adatoms that corresponds to i=0. For square islands, deviations from PE predictions appear for coverages near θ=0.2 and are mainly related to island coalescence. Our results show that the range of applicability of the PE theory is narrow, thus observing the predicted Gaussian tail of CZD may be difficult in real systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic modelling of epitaxial film growth with up- and downward step barriers

The formation of three-dimensional structures during the epitaxial growth of films is associated to the reflection of diffusing particles in descending terraces due to the presence of the so-called Ehrlich-Schwoebel (ES) barrier. We generalize this concept in a solid-on-solid growth model, in which a barrier dependent on the particle coordination (number of lateral bonds) exists whenever the particle performs an interlayer diffusion. The rules do not distinguish explicitly if the particle is executing a descending or an ascending interlayer diffusion. We show that the usual model, with a step barrier in descending steps, produces spurious, columnar, and highly unstable morphologies if the growth temperature is varied in a usual range of mound formation experiments. Our model generates well-behaved mounded morphologies for the same ES barriers that produce anomalous morphologies in the standard model. Moreover, mounds are also obtained when the step barrier has an equal value for all particles independently if they are free or bonded. Kinetic roughening is observed at long times, when the surface roughness w and the characteristic length $ξ$ scale as $w ~ t^β$ and $ξ~ t^ζ$ where $β\approx 0.31$ and $ζ\approx 0.22$, independently of the growth temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Roughness exponents and grain shapes

In surfaces with grainy features, the local roughness $w$ shows a crossover at a characteristic length $r_c$, with roughness exponent changing from $α_1\approx 1$ to a smaller $α_2$. The grain shape, the choice of $w$ or height-height correlation function (HHCF) $C$, and the procedure to calculate root mean-square averages are shown to have remarkable effects on $α_1$. With grains of pyramidal shape, $α_1$ can be as low as 0.71, which is much lower than the previous prediction 0.85 for rounded grains. The same crossover is observed in the HHCF, but with initial exponent $χ_1\approx 0.5$ for flat grains, while for some conical grains it may increase to $χ_1\approx 0.7$. The universality class of the growth process determines the exponents $α_2=χ_2$ after the crossover, but has no effect on the initial exponents $α_1$ and $χ_1$, supporting the geometric interpretation of their values. For all grain shapes and different definitions of surface roughness or HHCF, we still observe that the crossover length $r_c$ is an accurate estimate of the grain size. The exponents obtained in several recent experimental works on different materials are explained by those models, with some surface images qualitatively similar to our model films.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Scaling of island size and capture zone distributions in submonolayer growth

Island size and capture zone distributions (ISDs, CZDs) are studied numerically in submonolayer growth with various critical island sizes and shapes. CZDs scaled by the variance show excellent agreement with the Wigner surmise, confirming the Pimpinelli-Einstein approach for large CZs / large island dynamics. The ISD decay as exp(-s^x), with x=4, x approx 2.4 and 2 for point, fractal, and square islands, respectively. A scaling approach explains the values of x from the Gaussian decay of CZD and the efficiency of islands to capture diffusing adatoms.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Solution on the Bethe lattice of a hard core athermal gas with two kinds of particles

Athermal lattice gases of particles with first neighbor exclusion have been studied for a long time as simple models exhibiting a fluid-solid transition. At low concentration the particles occupy randomly both sublattices, but as the concentration is increased one of the sublattices is occupied preferentially. Here we study a mixed lattice gas with excluded volume interactions only in the grand-canonical formalism with two kinds of particles: small ones, which occupy a single lattice site and large ones, which occupy one site and its first neighbors. We solve the model on a Bethe lattice of arbitrary coordination number $q$. In the parameter space defined by the activities of both particles. At low values of the activity of small particles ($z_1$) we find a continuous transition from the fluid to the solid phase as the activity of large particles ($z_2$) is increased. At higher values of $z_1$ the transition becomes discontinuous, both regimes are separated by a tricritical point. The critical line has a negative slope at $z_1=0$ and displays a minimum before reaching the tricritical point, so that a reentrant behavior is observed for constant values of $z_2$ in the region of low density of small particles. The isobaric curves of the total density of particles as a function of $z_1$ (or $z_2$) show a minimum in the fluid phase.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Universal fluctuations in KPZ growth on one-dimensional flat substrates

We present a numerical study of the evolution of height distributions (HDs) obtained in interface growth models belonging to the Kardar-Parisi-Zhang (KPZ) universality class. The growth is done on an initially flat substrate. The HDs obtained for all investigated models are very well fitted by the theoretically predicted Gaussian Orthogonal Ensemble (GOE) distribution. The first cumulant has a shift that vanishes as $t^{-1/3}$, while the cumulants of order $2\le n\le 4$ converge to GOE as $t^{-2/3}$ or faster, behaviors previously observed in other KPZ systems. These results yield a new evidence for the universality of the GOE distribution in KPZ growth on flat substrates. Finally, we further show that the surfaces are described by the Airy$_{1}$ process.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Universal fluctuations in radial growth models belonging to the KPZ universality class

We investigate the radius distributions (RD) of surfaces obtained with large-scale simulations of radial clusters that belong to the KPZ universality class. For all investigated models, the RDs are given by the Tracy-Widom distribution of the Gaussian unitary ensemble, in agreement with the conjecture of the KPZ universality class for curved surfaces. The quantitative agreement was also confirmed by two-point correlation functions asymptotically given by the covariance of the Airy$_2$ process. Our simulation results fill the last lacking gap of the conjecture that had been recently verified analytically and experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Bethe lattice solution of a model of SAW's with up to 3 monomers per site and no restriction

In the multiple monomers per site (MMS) model, polymeric chains are represented by walks on a lattice which may visit each site up to K times. We have solved the unrestricted version of this model, where immediate reversals of the walks are allowed (RA) for K = 3 on a Bethe lattice with arbitrary coordination number in the grand-canonical formalism. We found transitions between a non-polymerized and two polymerized phases, which may be continuous or discontinuous. In the canonical situation, the transitions between the extended and the collapsed polymeric phases are always continuous. The transition line is partly composed by tricritical points and partially by critical endpoints, both lines meeting at a multicritical point. In the subspace of the parameter space where the model is related to SASAW's (self-attracting self-avoiding walks), the collapse transition is tricritical. We discuss the relation of our results with simulations and previous Bethe and Husimi lattice calculations for the MMS model found in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Solution of an associating lattice gas model with density anomaly on a Husimi lattice

We study a model of a lattice gas with orientational degrees of freedom which resemble the formation of hydrogen bonds between the molecules. In this model, which is the simplified version of the Henriques-Barbosa model, no distinction is made between donors and acceptors in the bonding arms. We solve the model in the grand-canonical ensemble on a Husimi lattice built with hexagonal plaquettes with a central site. The ground-state of the model, which was originally defined on the triangular lattice, is exactly reproduced by the solution on this Husimi lattice. In the phase diagram, one gas and two liquid (high density-HDL and low density-LDL) phases are present. All phase transitions (GAS-LDL, GAS-HDL, and LDL-HDL) are discontinuous, and the three phases coexist at a triple point. A line of temperatures of maximum density (TMD) in the isobars is found in the metastable GAS phase, as well as another line of temperatures of minimum density (TmD) appears in the LDL phase, part of it in the stable region and another in the metastable region of this phase. These findings are at variance with simulational results for the same model on the triangular lattice, which suggested a phase diagram with two critical points. However, our results show very good quantitative agreement with the simulations, both for the coexistence loci and the densities of particles and of hydrogen bonds. We discuss the comparison of the simulations with our results.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Grand canonical and canonical solution of self-avoiding walks with up to three monomers per site on the Bethe lattice

We solve a model of polymers represented by self-avoiding walks on a lattice which may visit the same site up to three times in the grand-canonical formalism on the Bethe lattice. This may be a model for the collapse transition of polymers where only interactions between monomers at the same site are considered. The phase diagram of the model is very rich, displaying coexistence and critical surfaces, critical, critical endpoint and tricritical lines, as well as a multicritical point. From the grand-canonical results, we present an argument to obtain the properties of the model in the canonical ensemble, and compare our results with simulations in the literature. We do actually find extended and collapsed phases, but the transition between them, composed by a line of critical endpoints and a line of tricritical points, separated by the multicritical point, is always continuous. This result is at variance with the simulations for the model, which suggest that part of the line should be a discontinuous transition. Finally, we discuss the connection of the present model with the standard model for the collapse of polymers (self-avoiding self-attracting walks), where the transition between the extended and collapsed phases is a tricritical point.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Universal and nonuniversal features in the crossover from linear to nonlinear interface growth

We study a restricted solid-on-solid (RSOS) model involving deposition and evaporation with probabilities p and 1-p, respectively, in one-dimensional substrates. It presents a crossover from Edwards-Wilkinson (EW) to Kardar-Parisi-Zhang (KPZ) scaling for p~0.5. The associated KPZ equation is analytically derived, exhibiting a coefficient lambda of the nonlinear term proportional to q=p-1/2, which is confirmed numerically by calculation of tilt-dependent growth velocities for several values of p. This linear λ-q relation contrasts to the apparently universal parabolic law obtained in competitive models mixing EW and KPZ components. The regions where the interface roughness shows pure EW and KPZ scaling are identified for 0.55<=p<=0.8, which provides numerical estimates of the crossover times t_c. They scale as t_c ~ lambda^(-phi) with phi=4.1+-0.1, which is in excellent agreement with the theoretically predicted universal value phi=4 and improves previous numerical estimates, which suggested phi~3.

0 citations0 reviewsNo code linkedOpen access