Researcher profile

J. L. Lebowitz

J. L. Lebowitz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

28 published item(s)

preprint2022arXiv

Stationary States of the One-Dimensional Discrete-Time Facilitated Symmetric Exclusion Process

We describe the extremal translation invariant stationary (ETIS) states of the facilitated exclusion process on $\mathbb{Z}$. In this model all particles on sites with one occupied and one empty neighbor jump at each integer time to the empty neighbor site, and if two particles attempt to jump into the same empty site we choose one randomly to succeed. The ETIS states are qualitatively different for densities $ρ<1/2$, $ρ=1/2$, and $1/2<ρ<1$, but in each density region we find states which may be grouped into families, each of which is in natural correspondence with the set of all ergodic measures on $\{0,1\}^{\mathbb{Z}}$. For $ρ<1/2$ there is one such family, containing all the ergodic states in which the probability of two adjacent occupied sites is zero. For $ρ=1/2$ there are two families, in which configurations translate to the left and right, respectively, with constant speed 2. For the high density case there is a continuum of families. We show that all ETIS states at densities $ρ\le1/2$ belong to these families, and conjecture that also at high density there are no other ETIS states. We also study the possible ETIS states which might occur if the conjecture fails.

preprint2020arXiv

The Discrete-Time Facilitated Totally Asymmetric Simple Exclusion Process

We describe the translation invariant stationary states of the one dimensional discrete-time facilitated totally asymmetric simple exclusion process (F-TASEP). In this system a particle at site $j$ in $Z$ jumps, at integer times, to site $j+1$, provided site $j-1$ is occupied and site $j+1$ is empty. This defines a deterministic noninvertible dynamical evolution from any specified initial configuration on $\{0,1\}^{Z}$. When started with a Bernoulli product measure at density $ρ$ the system approaches a stationary state, with phase transitions at $ρ=1/2$ and $ρ=2/3$. We discuss various properties of these states in the different density regimes $0<ρ<1/2$, $1/2<ρ<2/3$, and $2/3<ρ<1$; for example, we show that the pair correlation $g(j)=\langleη(i)η(i+j)\rangle$ satisfies, for all $n\in Z$, $\sum_{j=kn+1}^{k(n+1)}g(j)=kρ^2$, with $k=2$ when $0 \le ρ\le 1/2$ and $k=3$ when $2/3 \le ρ\le 1$, and conjecture (on the basis of simulations) that the same identity holds with $k=6$ when $1/2 \le ρ\le 2/3$. The $ρ<1/2$ stationary state referred to above is also the stationary state for the deterministic discrete-time TASEP at density $ρ$ (with Bernoulli initial state) or, after exchange of particles and holes, at density $1-ρ$.

preprint2019arXiv

Exact Solution of the F-TASEP

We obtain the exact solution of the facilitated totally asymmetric simple exclusion process (F-TASEP) in 1D. The model is closely related to the conserved lattice gas (CLG) model and to some cellular automaton traffic models. In the F-TASEP a particle at site $j$ in $\mathbb{Z}$ jumps, at integer times, to site $j+1$, provided site $j-1$ is occupied and site $j+1$ is empty. When started with a Bernoulli product measure at density $ρ$ the system approaches a stationary state. This non-equilibrium steady state (NESS) has phase transitions at $ρ=1/2$ and $ρ=2/3$. The different density regimes $0<ρ<1/2$, $1/2<ρ<2/3$, and $2/3<ρ<1$ exhibit many surprising properties; for example, the pair correlation $g(j)=\langleη(i)η(i+j)\rangle$ satisfies, for all $n\in\mathbb{Z}$, $\sum_{j=kn+1}^{k(n+1)}g(j)=kρ^2$, with $k=2$ when $0\leρ\le1/2$, $k=6$ when $1/2\leρ\le2/3$, and $k=3$ when $2/3\leρ\le1$. The quantity $\lim_{L\to\infty}V_L/L$, where $V_L$ is the variance in the number of particles in an interval of length L, jumps discontinuosly from $ρ(1-ρ)$ to 0 when $ρ\to1/2$ and when $ρ\to2/3$.

preprint2016arXiv

Translation invariant extensions of finite volume measures

We investigate the following questions: Given a measure $μ_Λ$ on configurations on a subset $Λ$ of a lattice $\mathbb{L}$, where a configuration is an element of $Ω^Λ$ for some fixed set $Ω$, does there exist a measure $μ$ on configurations on all of $\mathbb{L}$, invariant under some specified symmetry group of $\mathbb{L}$, such that $μ_Λ$ is its marginal on configurations on $Λ$? When the answer is yes, what are the properties, e.g., the entropies, of such measures? Our primary focus is the case in which $\mathbb{L}=\mathbb{Z}^d$ and the symmetries are the translations. For the case in which $Λ$ is an interval in $\mathbb{Z}$ we give a simple necessary and sufficient condition, local translation invariance (LTI), for extendibility. For LTI measures we construct extensions having maximal entropy, which we show are Gibbs measures; this construction extends to the case in which $\mathbb{L}$ is the Bethe lattice. On $\mathbb{Z}$ we also consider extensions supported on periodic configurations, which are analyzed using de~Bruijn graphs and which include the extensions with minimal entropy. When $Λ\subset\mathbb{Z}$ is not an interval, or when $Λ\subset\mathbb{Z}^d$ with $d>1$, the LTI condition is necessary but not sufficient for extendibility. For $\mathbb{Z}^d$ with $d>1$, extendibility is in some sense undecidable.

preprint2015arXiv

Central limit theorems, Lee-Yang zeros, and graph-counting polynomials

We consider the asymptotic normalcy of families of random variables $X$ which count the number of occupied sites in some large set. We write $Prob(X=m)=p_mz_0^m/P(z_0)$, where $P(z)$ is the generating function $P(z)=\sum_{j=0}^{N}p_jz^j$ and $z_0>0$. We give sufficient criteria, involving the location of the zeros of $P(z)$, for these families to satisfy a central limit theorem (CLT) and even a local CLT (LCLT); the theorems hold in the sense of estimates valid for large $N$ (we assume that $Var(X)$ is large when $N$ is). For example, if all the zeros lie in the closed left half plane then $X$ is asymptotically normal, and when the zeros satisfy some additional conditions then $X$ satisfies an LCLT. We apply these results to cases in which $X$ counts the number of edges in the (random) set of &#34;occupied&#34; edges in a graph, with constraints on the number of occupied edges attached to a given vertex. Our results also apply to systems of interacting particles, with $X$ counting the number of particles in a box $Λ$ whose size approaches infinity; $P(z)$ is then the grand canonical partition function and its zeros are the Lee-Yang zeros.

preprint2014arXiv

KPZ universality class and the anchored Toom interface

We revisit the anchored Toom interface and use KPZ scaling theory to argue that the interface fluctuations are governed by the Airy_1 process with the role of space and time interchanged. There is no free parameter. The predictions are numerically well confirmed for space-time statistics in the stationary state. In particular the spatial fluctuations of the interface are given by the GOE edge distribution of Tracy and Widom.

preprint2013arXiv

Lyapunov functionals for boundary-driven nonlinear drift-diffusions

We exhibit a large class of Lyapunov functionals for nonlinear drift-diffusion equations with non-homogeneous Dirichlet boundary conditions. These are generalizations of large deviation functionals for underlying stochastic many-particle systems, the zero range process and the Ginzburg-Landau dynamics, which we describe briefly. As an application, we prove linear inequalities between such an entropy-like functional and its entropy production functional for the boundary-driven porous medium equation in a bounded domain with positive Dirichlet conditions: this implies exponential rates of relaxation related to the first Dirichlet eigenvalue of the domain. We also derive Lyapunov functions for systems of nonlinear diffusion equations, and for nonlinear Markov processes with non-reversible stationary measures.

preprint2013arXiv

Propagation of Chaos for a Thermostated Kinetic Model

We consider a system of N point particles moving on a d-dimensional torus. Each particle is subject to a uniform field E and random speed conserving collisions. This model is a variant of the Drude-Lorentz model of electrical conduction. In order to avoid heating by the external field, the particles also interact with a Gaussian thermostat which keeps the total kinetic energy of the system constant. The thermostat induces a mean-field type of interaction between the particles. Here we prove that, starting from a product measure, in the large N limit, the one particle velocity distribution satisfies a self consistent Vlasov-Boltzmann equation.. This is a consequence of &#34;propagation of chaos&#34;, which we also prove for this model.

preprint2013arXiv

Time evolution of electron flow in a model diode: Non-perturbative analysis

Using a combination of Eulerian and Lagrangian variables we obtain some exact results and good approximation schemes for the time evolution of the electron flow from a no-current state to a final stationary current state in a planar one-dimensional diode. The electrons can be injected externally or generated by the cathode via field emission governed by a current-field law. The case of equipotential electrodes and fixed injection is studied along with a positive anode potential. When the current is fixed externally the approach to the stationary state goes without oscillations if the initial electron velocity is high enough and the anode can absorb the injected flow. Otherwise the accumulated space charge creates a potential barrier which reflects the flow and leads to its oscillations, but our method of analysis is invalid in such conditions. In the field emission case the flow goes to its stationary state through a train of decaying oscillations whose period is of the order of the electron transit time, in agreement with earlier studies based on perturbation techniques. Our approximate method does not permit very high cathode emissivity, though the method works when the stationary current density is only about 10% smaller than the Child-Langmuir limit.

preprint2012arXiv

Approach to equilibrium for the stochastic NLS

We study the approach to equilibrium, described by a Gibbs measure, for a system on a $d$-dimensional torus evolving according to a stochastic nonlinear Schrödinger equation (SNLS) with a high frequency truncation. We prove exponential approach to the truncated Gibbs measure both for the focusing and defocusing cases when the dynamics is constrained via suitable boundary conditions to regions of the Fourier space where the Hamiltonian is convex. Our method is based on establishing a spectral gap for the non self-adjoint Fokker-Planck operator governing the time evolution of the measure, which is {\it uniform} in the frequency truncation $N$. The limit $N\to\infty$ is discussed.

preprint2012arXiv

Non-equilibrium stationary state of a harmonic crystal with alternating masses

We analyze the non-equilibrium steady states (NESS) of a one dimensional harmonic chain of $N$ atoms with alternating masses connected to heat reservoirs at unequal temperatures. We find that the temperature profile defined through the local kinetic energy $T(j) \equiv {<p_j^2>}/{m_j}$, oscillates with period two in the bulk of the system. Depending on boundary conditions, either the heavier or the lighter particles in the bulk are hotter. We obtain exact expressions for the bulk temperature profile and steady state current in the limit $N \rightarrow \infty$. These depend on whether $N$ is odd or even. We also study similar temperature oscillations in the NESS of systems with noise in the dynamics. These die out as $N \rightarrow \infty$.

preprint2012arXiv

The blockage problem

We investigate the totally asymmetric exclusion process on Z, with the jump rate at site i given by r_i=1 for i nonzero, r_0=r. It is easy to see that the maximal stationary current j(r) is nondecreasing in r and that j(r)=1/4 for r>=1; it is a long outstanding problem to determine whether or not the critical value r_c of r such that j(r)=1/4 for r>r_c is strictly less than 1. Here we present a heuristic argument, based on the analysis of the first sixteen terms in a formal power series expansion of j(r) obtained from finite volume systems, that r_c=1 and that for r less than 1 and near 1, j(r) behaves as 1/4-γ\exp[-{a/(1-r)}] with a approximately equal to 2. We also give some new exact results about this system; in particular we prove that j(r)=J_max(r), with J_max(r) the hydrodynamic maximal current defined by Seppalainen, and thus establish continuity of j(r). Finally we describe a related exactly solvable model, a semi-infinite system in which the site i=0 is always occupied. For that system, the critical r is 1/2 and the analogue j_s(r) of j(r) satisfies j_s(r)=r(1-r) for r<=1/2; j_s(r) is the limit of finite volume currents inside the curve |r(1-r)|=1/4 in the complex r plane and we suggest that analogous behavior may hold for the original system.

preprint2011arXiv

Analysis of Droplets in Lattice Systems with Long-range Kac Potentials

We investigate the geometry of typical equilibrium configurations for a lattice gas in a finite macroscopic domain with attractive, long range Kac potentials. We focus on the case when the system is below the critical temperature and has a fixed number of occupied sites.We connect the properties of typical configurations to the analysis of the constrained minimizers of a mesoscopic non-local free energy functional, which we prove to be the large deviation functional for a density profile in the canonical Gibbs measure with prescribed global density. In the case in which the global density of occupied sites lies between the two equilibrium densities that one would have without a constraint on the particle number, a &#34;droplet&#34; of the high (low) density phase may or may not form in a background of the low (high) density phase. We determine the critical density for droplet formation, and the nature of the droplet, as a function of the temperature and the size of the system, by combining the present large deviation principle with the analysis of the mesoscopic functional given in CCELM

preprint2010arXiv

A diffusive system driven by a battery or by a smoothly varying field

We consider the steady state of a one dimensional diffusive system, such as the symmetric simple exclusion process (SSEP) on a ring, driven by a battery at the origin or by a smoothly varying field along the ring. The battery appears as the limiting case of a smoothly varying field, when the field becomes a delta function at the origin. We find that in the scaling limit, the long range pair correlation functions of the system driven by a battery turn out to be very different from the ones known in the steady state of the SSEP maintained out of equilibrium by contact with two reservoirs, even when the steady state density profiles are identical in both models.

preprint2010arXiv

Ionization of Coulomb systems in $\RR^3$ by time periodic forcings of arbitrary size

We analyze the long time behavior of solutions of the Schrödinger equation $iψ_t=(-Δ-b/r+V(t,x))ψ$, $x\in\RR^3$, $r=|x|$, describing a Coulomb system subjected to a spatially compactly supported time periodic potential $V(t,x)=V(t+2π/ω,x)$ with zero time average. We show that, for any $V(t,x)$ of the form $2Ω(r)\sin (ωt-θ)$, with $Ω(r)$ nonzero on its support, Floquet bound states do not exist. This implies that the system ionizes, {\em i.e.} $P(t,K)=\int_K|ψ(t,x)|^2dx\to 0$ as $t\to\infty$ for any compact set $K\subset\RR^3$. Furthermore, if the initial state is compactly supported and has only finitely many spherical harmonic modes, then $P(t,K)$ decays like $t^{-5/3}$ as $t \to \infty $. To prove these statements, we develop a rigorous WKB theory for infinite systems of ordinary differential equations.

preprint2009arXiv

Droplet minimizers for the Gates-Lebowitz-Penrose free energy functional

We study the structure of the constrained minimizers of the Gates-Lebowitz-Penrose free-energy functional ${\mathcal F}_{\rm GLP}(m)$, non-local functional of a density field $m(x)$, $x\in {\mathcal T}_L$, a $d$-dimensional torus of side length $L$. At low temperatures, ${\mathcal F}_{\rm GLP}$ is not convex, and has two distinct global minimizers, corresponding to two equilibrium states. Here we constrain the average density $L^{-d}\int_{{\cal T}_L}m(x)\dd x$ to be a fixed value $n$ between the densities in the two equilibrium states, but close to the low density equilibrium value. In this case, a &#34;droplet&#34; of the high density phase may or may not form in a background of the low density phase, depending on the values $n$ and $L$. We determine the critical density for droplet formation, and the nature of the droplet, as a function of $n$ and $L$. The relation between the free energy and the large deviations functional for a particle model with long-range Kac potentials, proven in some cases, and expected to be true in general, then provides information on the structure of typical microscopic configurations of the Gibbs measure when the range of the Kac potential is large enough.

preprint2009arXiv

Phase diagram of the ABC model on an interval

The three species asymmetric ABC model was initially defined on a ring by Evans, Kafri, Koduvely, and Mukamel, and the weakly asymmetric version was later studied by Clincy, Derrida, and Evans. Here the latter model is studied on a one-dimensional lattice of N sites with closed (zero flux) boundaries. In this geometry the local particle conserving dynamics satisfies detailed balance with respect to a canonical Gibbs measure with long range asymmetric pair interactions. This generalizes results for the ring case, where detailed balance holds, and in fact the steady state measure is known only for the case of equal densities of the different species: in the latter case the stationary states of the system on a ring and on an interval are the same. We prove that in the N to infinity limit the scaled density profiles are given by (pieces of) the periodic trajectory of a particle moving in a quartic confining potential. We further prove uniqueness of the profiles, i.e., the existence of a single phase, in all regions of the parameter space (of average densities and temperature) except at low temperature with all densities equal; in this case a continuum of phases, differing by translation, coexist. The results for the equal density case apply also to the system on the ring, and there extend results of Clincy et al.

preprint2007arXiv

Entropy of Open Lattice Systems

We investigate the behavior of the Gibbs-Shannon entropy of the stationary nonequilibrium measure describing a one-dimensional lattice gas, of L sites, with symmetric exclusion dynamics and in contact with particle reservoirs at different densities. In the hydrodynamic scaling limit, L to infinity, the leading order (O(L)) behavior of this entropy has been shown by Bahadoran to be that of a product measure corresponding to strict local equilibrium; we compute the first correction, which is O(1). The computation uses a formal expansion of the entropy in terms of truncated correlation functions; for this system the k-th such correlation is shown to be O(L^{-k+1}). This entropy correction depends only on the scaled truncated pair correlation, which describes the covariance of the density field. It coincides, in the large L limit, with the corresponding correction obtained from a Gaussian measure with the same covariance.

preprint2006arXiv

Time asymptotics of the Schroedinger wave function in time-periodic potentials

We study the transition to the continuum of an initially bound quantum particle in $\RR^d$, $d=1,2,3$, subjected, for $t\ge 0$, to a time periodic forcing of arbitrary magnitude. The analysis is carried out for compactly supported potentials, satisfying certain auxiliary conditions. It provides complete analytic information on the time Laplace transform of the wave function. From this, comprehensive time asymptotic properties (Borel summable transseries) follow. We obtain in particular a criterion for whether the wave function gets fully delocalized (complete ionization). This criterion shows that complete ionization is generic and provides a convenient test for particular cases. When satisfied it implies absence of discrete spectrum and resonances of the associated Floquet operator. As an illustration we show that the parametric harmonic perturbation of a potential chosen to be any nonzero multiple of the characteristic function of a measurable compact set has this property.

preprint2002arXiv

Dynamics of a Massive Piston in an Ideal Gas: Oscillatory Motion and Approach to Equilibrium

We study numerically and theoretically (on a heuristic level) the time evolution of a gas confined to a cube of size $L^3$ divided into two parts by a piston with mass $M_L \sim L^2$ which can only move in the $x$-direction. Starting with a uniform ``double-peaked&#39;&#39; (non Maxwellian) distribution of the gas and a stationary piston, we find that (a) after an initial quiescent period the system becomes unstable and the piston performs a damped oscillatory motion, and (b) there is a thermalization of the system leading to a Maxwellian distribution of the gas velocities. The time of the onset of the instability appears to grow like $L \log L$ while the relaxation time to the Maxwellian grows like $L^{7/2}$.

preprint2002arXiv

Large Deviation of the Density Profile in the Steady State of the Open Symmetric Simple Exclusion Process

We consider an open one dimensional lattice gas on sites $i=1,...,N$, with particles jumping independently with rate 1 to neighboring interior empty sites, the {\it simple symmetric exclusion process}. The particle fluxes at the left and right boundaries, corresponding to exchanges with reservoirs at different chemical potentials, create a stationary nonequilibrium state (SNS) with a steady flux of particles through the system. The mean density profile in this state, which is linear, describes the typical behavior of a macroscopic system, i.e., this profile occurs with probability 1 when $N \to \infty$. The probability of microscopic configurations corresponding to some other profile $ρ(x)$, $x = i/N$, has the asymptotic form $\exp[-N {\cal F}(\{ρ\})]$; $\cal F$ is the {\it large deviation functional}. In contrast to equilibrium systems, for which ${\cal F}_{eq}(\{ρ\})$ is just the integral of the appropriately normalized local free energy density, the $\cal F$ we find here for the nonequilibrium system is a nonlocal function of $ρ$. This gives rise to the long range correlations in the SNS predicted by fluctuating hydrodynamics and suggests similar non-local behavior of $\cal F$ in general SNS, where the long range correlations have been observed experimentally.

preprint2002arXiv

Scaling Dynamics of a Massive Piston in a Cube Filled With Ideal Gas: Exact Results

We continue the study of the time evolution of a system consisting of a piston in a cubical container of large size $L$ filled with an ideal gas. The piston has mass $M\sim L^2$ and undergoes elastic collisions with $N\sim L^3$ gas particles of mass $m$. In a previous paper, Lebowitz, Piasecki and Sinai considered a scaling regime, with time and space scaled by $L$, in which they argued heuristically that the motion of the piston and the one particle distribution of the gas satisfy autonomous coupled differential equations. Here we state exact results and sketch proofs for this behavior.

preprint2000arXiv

Diffusion Effects on the Breakdown of a Linear Amplifier Model Driven by the Square of a Gaussian Field

We investigate solutions to the equation $\partial_t{\cal E} - {\cal D}Δ{\cal E} = λS^2{\cal E}$, where $S(x,t)$ is a Gaussian stochastic field with covariance $C(x-x&#39;,t,t&#39;)$, and $x\in {\mathbb R}^d$. It is shown that the coupling $λ_{cN}(t)$ at which the $N$-th moment $<{\cal E}^N(x,t)>$ diverges at time $t$, is always less or equal for ${\cal D}>0$ than for ${\cal D}=0$. Equality holds under some reasonable assumptions on $C$ and, in this case, $λ_{cN}(t)=Nλ_c(t)$ where $λ_c(t)$ is the value of $λ$ at which $<\exp\lbrack λ\int_0^tS^2(0,s)ds\rbrack>$ diverges. The ${\cal D}=0$ case is solved for a class of $S$. The dependence of $λ_{cN}(t)$ on $d$ is analyzed. Similar behavior is conjectured when diffusion is replaced by diffraction, ${\cal D}\to i{\cal D}$, the case of interest for backscattering instabilities in laser-plasma interaction.

preprint1998arXiv

Liquid-Vapor Phase Transitions for Systems with Finite Range Interactions

We consider particles in $\R^d, d \geq 2$ interacting via attractive pair and repulsive four-body potentials of the Kac type. Perturbing about mean field theory, valid when the interaction range becomes infinite, we prove rigorously the existence of a liquid-gas phase transition, when the interaction range is finite but long compared to the interparticle spacing for a range of temperature.

preprint1998arXiv

Modelling of Phase Separation in Alloys with Coherent Elastic Misfit

Elastic interactions arising from a difference of lattice spacing between two coherent phases can have a strong influence on the phase separation (coarsening) of alloys. If the elastic moduli are different in the two phases, the elastic interactions may accelerate, slow down or even stop the phase separation process. If the material is elastically anisotropic, the precipitates can be shaped like plates or needles instead of spheres and can form regular precipitate superlattices. Tensions or compressions applied externally to the specimen may have a strong effect on the shapes and arrangement of the precipitates. In this paper, we review the main theoretical approaches that have been used to model these effects and we relate them to experimental observations. The theoretical approaches considered are (i) `macroscopic&#39; models treating the two phases as elastic media separated by a sharp interface (ii) `mesoscopic&#39; models in which the concentration varies continuously across the interface (iii) `microscopic&#39; models which use the positions of individual atoms.

preprint1998arXiv

Self-Diffusion in Simple Models: Systems with Long-Range Jumps

We review some exact results for the motion of a tagged particle in simple models. Then, we study the density dependence of the self diffusion coefficient, $D_N(ρ)$, in lattice systems with simple symmetric exclusion in which the particles can jump, with equal rates, to a set of $N$ neighboring sites. We obtain positive upper and lower bounds on $F_N(ρ)=N((1-\r)-[D_N(ρ)/D_N(0)])/(ρ(1-ρ))$ for $ρ\in [0,1]$. Computer simulations for the square, triangular and one dimensional lattice suggest that $F_N$ becomes effectively independent of $N$ for $N\ge 20$.

preprint1997arXiv

Metastability in the two-dimensional Ising model with free boundary conditions

We investigate metastability in the two dimensional Ising model in a square with free boundary conditions at low temperatures. Starting with all spins down in a small positive magnetic field, we show that the exit from this metastable phase occurs via the nucleation of a critical droplet in one of the four corners of the system. We compute the lifetime of the metastable phase analytically in the limit $T\to 0$, $h\to 0$ and via Monte Carlo simulations at fixed values of $T$ and $h$ and find good agreement. This system models the effects of boundary domains in magnetic storage systems exiting from a metastable phase when a small external field is applied.

preprint1997arXiv

Phase Transitions in the Multicomponent Widom-Rowlinson Model and in Hard Cubes on the BCC--Lattice

We use Monte Carlo techniques and analytical methods to study the phase diagram of the M--component Widom-Rowlinson model on the bcc-lattice: there are M species all with the same fugacity z and a nearest neighbor hard core exclusion between unlike particles. Simulations show that for M greater or equal 3 there is a ``crystal phase&#39;&#39; for z lying between z_c(M) and z_d(M) while for z > z_d(M) there are M demixed phases each consisting mostly of one species. For M=2 there is a direct second order transition from the gas phase to the demixed phase while for M greater or equal 3 the transition at z_d(M) appears to be first order putting it in the Potts model universality class. For M large, Pirogov-Sinai theory gives z_d(M) ~ M-2+2/(3M^2) + ... . In the crystal phase the particles preferentially occupy one of the sublattices, independent of species, i.e. spatial symmetry but not particle symmetry is broken. For M to infinity this transition approaches that of the one component hard cube gas with fugacity y = zM. We find by direct simulations of such a system a transition at y_c ~ 0.71 which is consistent with the simulation z_c(M) for large M. This transition appears to be always of the Ising type.