Source author record

Joel L. Lebowitz

Joel L. Lebowitz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

49 published item(s)

preprint2026arXiv

Convergent Power Series for Anharmonic Chain with Periodic Forcing

We study the propagation of energy in one-dimensional anharmonic chains subject to a periodic, localized forcing. For the purely harmonic case, forcing frequencies outside the linear spectrum produce exponentially localized responses, preventing equi-distribution of energy per degree of freedom. We extend this result to anharmonic perturbations with bounded second derivatives and boundary dissipation, proving that for small perturbations and non-resonant forcing, the dynamics converges to a periodic stationary state with energy exponentially localized uniformly in the system size. The perturbed periodic state is described by a convergent power type expansion in the strength of the anharmonicity. This excludes chaoticity induced by anharmonicity, independently of the size of the system. Our perturbative scheme can also be applied in higher dimensions.

preprint2022arXiv

Fluctuation and Entropy in Spectrally Constrained random fields

We investigate the statistical properties of translation invariant random fields (including point processes) on Euclidean spaces (or lattices) under constraints on their spectrum or structure function. An important class of models that motivate our study are hyperuniform and stealthy hyperuniform systems, which are characterised by the vanishing of the structure function at the origin (resp., vanishing in a neighbourhood of the origin). We show that many key features of two classical statistical mechanical measures of randomness - namely, fluctuations and entropy, are governed only by some particular local aspects of their structure function. We obtain exponents for the fluctuations of the local mass in domains of growing size, and show that spatial geometric considerations play an important role - both the shape of the domain and the mode of spectral decay. In doing so, we unveil intriguing oscillatory behaviour of spatial correlations of local masses in adjacent box domains. We describe very general conditions under which we show that the field of local masses exhibit Gaussian asymptotics, with an explicitly described limit. We further demonstrate that stealthy hyperuniform systems with joint densities exhibit degeneracy in their asymptotic entropy per site. In fact, our analysis shows that entropic degeneracy sets in under much milder conditions than stealthiness, as soon as the structure function fails to be logarithmically integrable.

preprint2020arXiv

Exact solution of the Schrodinger equation for photoemission from a metal

We solve rigorously the time dependent Schrödinger equation describing electron emission from a metal surface by a laser field perpendicular to the surface. We consider the system to be one-dimensional, with the half-line $x<0$ corresponding to the bulk of the metal and $x>0$ to the vacuum. The laser field is modeled as a classical electric field oscillating with frequency $ω$, acting only at $x>0$. We consider an initial condition which is a stationary state of the system without a field, and, at time $t=0$, the field is switched on. We prove the existence of a solution $ψ(x,t)$ of the Schrödinger equation for $t>0$, and compute the surface current. The current exhibits a complex oscillatory behavior, which is not captured by the "simple" three step scenario. As $t\to\infty$, $ψ(x,t)$ converges with a rate $t^{-\frac32}$ to a time periodic function with period $\frac{2π}ω$ which coincides with that found by Faisal, Kamiński and Saczuk (Phys Rev A 72, 023412, 2015). However, for realistic values of the parameters, we have found that it can take quite a long time (over 50 laser periods) for the system to converge to its asymptote. Of particular physical importance is the current averaged over a laser period $\frac{2π}ω$, which exhibits a dramatic increase when $\hbarω$ becomes larger than the work function of the metal, which is consistent with the original photoelectric effect.

preprint2019arXiv

Diffusive limit for a Boltzmann-like equation with non-conserved momentum

We consider a kinetic model whose evolution is described by a Boltzmann-like equation for the one-particle phase space distribution $f(x,v,t)$. There are hard-sphere collisions between the particles as well as collisions with randomly fixed scatterers. As a result, this evolution does not conserve momentum but only mass and energy. We prove that the diffusively rescaled $f^\varepsilon(x,v,t)=f(\varepsilon^{-1}x,v,\varepsilon^{-2}t)$, as $\varepsilon\to 0$ tends to a Maxwellian $M_{ρ, 0, T}=\fracρ{(2πT)^{3/2}}\exp[{-\frac{|v|^2}{2T}}]$, where $ρ$ and $T$ are solutions of coupled diffusion equations and estimate the error in $L^2_{x,v}$.

preprint2019arXiv

Gibbs and Boltzmann Entropy in Classical and Quantum Mechanics

The Gibbs entropy of a macroscopic classical system is a function of a probability distribution over phase space, i.e., of an ensemble. In contrast, the Boltzmann entropy is a function on phase space, and is thus defined for an individual system. Our aim is to discuss and compare these two notions of entropy, along with the associated ensemblist and individualist views of thermal equilibrium. Using the Gibbsian ensembles for the computation of the Gibbs entropy, the two notions yield the same (leading order) values for the entropy of a macroscopic system in thermal equilibrium. The two approaches do not, however, necessarily agree for non-equilibrium systems. For those, we argue that the Boltzmann entropy is the one that corresponds to thermodynamic entropy, in particular in connection with the second law of thermodynamics. Moreover, we describe the quantum analog of the Boltzmann entropy, and we argue that the individualist (Boltzmannian) concept of equilibrium is supported by the recent works on thermalization of closed quantum systems.

preprint2019arXiv

Heat Conduction in a hard disc system with non-conserved momentum

We describe results of computer simulations of steady state heat transport in a fluid of hard discs undergoing both elastic interparticle collisions and velocity randomizing collisions which do not conserve momentum. The system consists of N discs of radius r in a unit square, periodic in the y-direction and having thermal walls at x = 0 with temperature T0 taking values from 1 to 20 and at x = 1 with T1 = 1. We consider different values of the ratio between randomizing and interparticle collision rates and extrapolate results from different N, to N->infinity, r->0 such that rho=1/2. We find that in the (extrapolated) limit N->infinity, the systems local density and temperature profiles are those of local thermodynamic equilibrium (LTE) and obey Fourier's law. The variance of global quantities, such as the total energy, deviates from its local equilibrium value in a form consistent with macroscopic fluctuation theory.

preprint2016arXiv

A note on Lee-Yang zeros in the negative half-plane

We obtain lower bounds on the inverse compressibility of systems whose Lee-Yang zeros of the grand-canonical partition function lie in the left half of the complex fugacity plane. This includes in particular systems whose zeros lie on the negative real axis such as the monomer-dimer system on a lattice. We also study the virial expansion of the pressure in powers of the density for such systems. We find no direct connection between the positivity of the virial coefficients and the negativity of the L-Y zeros, and provide examples of either one or both properties holding. An explicit calculation of the partition function of the monomer-dimer system on 2 rows shows that there are at most a finite number of negative virial coefficients in this case.

preprint2016arXiv

Fowler-Nordheim emission modified by laser pulses in the adiabatic regime

We investigate enhanced field emission due to a continuous or pulsed oscillating field added to a constant electric field $E$ at the emitter surface. When the frequency of oscillation, field strength, and property of the emitter material satisfy the Keldysh condition $γ<1/2$ one can use the adiabatic approximation for treating the oscillating field, i.e. consider the tunneling through the instantaneous Fowler-Nordheim barrier created by both fields. Due to the great sensitivity of the emission to the field strength the average tunneling current can be much larger than the current produced by only the constant field. We carry out the computations for arbitrary strong constant electric fields, beyond the commonly used Fowler-Nordheim approximation which exhibit in particular an important property of the wave function inside the potential barrier where it is found to be monotonically decreasing without oscillations.

preprint2015arXiv

Thermal equilibrium of a macroscopic quantum system in a pure state

We consider the notion of thermal equilibrium for an individual closed macroscopic quantum system in a pure state, i.e., described by a wave function. The macroscopic properties in thermal equilibrium of such a system, determined by its wave function, must be the same as those obtained from thermodynamics, e.g., spatial uniformity of temperature and chemical potential. When this is true we say that the system is in macroscopic thermal equilibrium (MATE). Such a system may however not be in microscopic thermal equilibrium (MITE). The latter requires that the reduced density matrices of small subsystems be close to those obtained from the microcanonical, equivalently the canonical, ensemble for the whole system. The distinction between MITE and MATE is particularly relevant for systems with many-body localization (MBL) for which the energy eigenfunctions fail to be in MITE while necessarily most of them, but not all, are in MATE. We note however that for generic macroscopic systems, including those with MBL, most wave functions in an energy shell are in both MATE and MITE. For a classical macroscopic system, MATE holds for most phase points on the energy surface, but MITE fails to hold for any phase point.

preprint2015arXiv

Universal Probability Distribution for the Wave Function of a Quantum System Entangled with Its Environment

A quantum system (with Hilbert space $\mathscr{H}_1$) entangled with its environment (with Hilbert space $\mathscr{H}_2$) is usually not attributed a wave function but only a reduced density matrix $ρ_1$. Nevertheless, there is a precise way of attributing to it a random wave function $ψ_1$, called its conditional wave function, whose probability distribution $μ_1$ depends on the entangled wave function $ψ\in\mathscr{H}_1\otimes\mathscr{H}_2$ in the Hilbert space of system and environment together. It also depends on a choice of orthonormal basis of $\mathscr{H}_2$ but in relevant cases, as we show, not very much. We prove several universality (or typicality) results about $μ_1$, e.g., that if the environment is sufficiently large then for every orthonormal basis of $\mathscr{H}_2$, most entangled states $ψ$ with given reduced density matrix $ρ_1$ are such that $μ_1$ is close to one of the so-called GAP (Gaussian adjusted projected) measures, $GAP(ρ_1)$. We also show that, for most entangled states $ψ$ from a microcanonical subspace (spanned by the eigenvectors of the Hamiltonian with energies in a narrow interval $[E,E+δE]$) and most orthonormal bases of $\mathscr{H}_2$, $μ_1$ is close to $GAP(\mathrm{tr}_2 ρ_{mc})$ with $ρ_{mc}$ the normalized projection to the microcanonical subspace. In particular, if the coupling between the system and the environment is weak, then $μ_1$ is close to $GAP(ρ_β)$ with $ρ_β$ the canonical density matrix on $\mathscr{H}_1$ at inverse temperature $β=β(E)$. This provides the mathematical justification of our claim in [J. Statist. Phys. 125:1193 (2006), http://arxiv.org/abs/quant-ph/0309021] that $GAP$ measures describe the thermal equilibrium distribution of the wave function.

preprint2014arXiv

Exponential approach to, and properties of, a non-equilibrium steady state in a dilute gas

We investigate a kinetic model of a system in contact with several thermal reservoirs at different temperatures $T_α$. Our system is a spatially uniform dilute gas whose internal dynamics is described by the nonlinear Boltzmann equation with Maxwellian collisions. Similarly, the interaction with reservoir $α$ is represented by a Markovian process that has the Maxwellian $M_{T_α}$ as its stationary state. We prove existence and uniqueness of a non-equilibrium steady state (NESS) and show exponential convergence to this NESS in a metric on probability measures introduced into the study of Maxwellian collisions by Gabetta, Toscani and Wenberg (GTW). This shows that the GTW distance between the current velocity distribution to the steady-state velocity distribution is a Lyapunov functional for the system. We also derive expressions for the entropy production in the system plus the reservoirs which is always positive.

preprint2014arXiv

Green-Kubo formula for weakly coupled system with dynamical noise

We study the Green-Kubo (GK) formula $κ(\varepsilon, ξ)$ for the heat conductivity of an infinite chain of $d$-dimensional finite systems (cells) coupled by a smooth nearest neighbour potential $\varepsilon V$. The uncoupled systems evolve according to Hamiltonian dynamics perturbed stochastically by an energy conserving noise of strength $ξ$. Noting that $κ(\varepsilon, ξ)$ exists and is finite whenever $ξ> 0$, we are interested in what happens when the strength of the noise $ξ\to 0$. For this, we start in this work by formally expanding $κ(\varepsilon, ξ)$ in a power series in $\varepsilon$, $κ(\varepsilon, ξ) = \varepsilon^2 \sum_{n\ge 2} \varepsilon^{n-2} κ_n (ξ)$ and investigating the (formal) equations satisfied by $κ_n (ξ$. We show in particular that $κ_2 (ξ)$ is well defined when no pinning potential is present, and coincides formally with the heat conductivity obtained in the weak coupling (van Hove) limit, where time is rescaled as $\varepsilon^{-2}t$, for the cases where the latter has been established \cite{LO, DL}. For one-dimensional systems, we investigate $κ_2 (ξ)$ as $ξ\to 0$ in three cases: the disordered harmonic chain, the rotor chain and a chain of strongly anharmonic oscillators. Moreover, we formally identify $κ_2 (ξ)$ with the conductivity obtained by having the chain between two reservoirs at temperature $T$ and $T+δT$, in the limit $δT\to 0$, $N \to \infty$, $\varepsilon \to 0$.

preprint2014arXiv

Local Central Limit Theorem for Determinantal Point Processes

We prove a local central limit theorem (LCLT) for the number of points $N(J)$ in a region $J$ in $\mathbb R^d$ specified by a determinantal point process with an Hermitian kernel. The only assumption is that the variance of $N(J)$ tends to infinity as $|J| \to \infty$. This extends a previous result giving a weaker central limit theorem (CLT) for these systems. Our result relies on the fact that the Lee-Yang zeros of the generating function for $\{E(k;J)\}$ --- the probabilities of there being exactly $k$ points in $J$ --- all lie on the negative real $z$-axis. In particular, the result applies to the scaled bulk eigenvalue distribution for the Gaussian Unitary Ensemble (GUE) and that of the Ginibre ensemble. For the GUE we can also treat the properly scaled edge eigenvalue distribution. Using identities between gap probabilities, the LCLT can be extended to bulk eigenvalues of the Gaussian Symplectic Ensemble (GSE). A LCLT is also established for the probability density function of the $k$-th largest eigenvalue at the soft edge, and of the spacing between $k$-th neigbors in the bulk.

preprint2013arXiv

Harmonic Systems With Bulk Noises

We consider a harmonic chain in contact with thermal reservoirs at different temperatures and subject to bulk noises of different types: velocity flips or self-consistent reservoirs. While both systems have the same covariances in the nonequilibrium stationary state (NESS) the measures are very different. We study hydrodynamical scaling, large deviations, fluctuations, and long range correlations in both systems. Some of our results extend to higher dimensions.

preprint2013arXiv

On an Extension Problem for Density Matrices

We investigate the problem of the existence of a density matrix rho on the product of three Hilbert spaces with given marginals on the pair (1,2) and the pair (2,3). While we do not solve this problem completely we offer partial results in the form of some necessary and some sufficient conditions on the two marginals. The quantum case differs markedly from the classical (commutative) case, where the obvious necessary compatibility condition suffices, namely, trace_1 (rho_{12}) = \trace_3 (rho_{23}).

preprint2012arXiv

Location of the Lee-Yang zeros and absence of phase transitions in some Ising spin systems

We consider a class of Ising spin systems on a set Λof sites. The sites are grouped into units with the property that each site belongs to either one or two units, and the total internal energy of the system is the sum of the energies of the individual units, which in turn depend only on the number of up spins in the unit. We show that under suitable conditions on these interactions none of the |Λ| Lee-Yang zeros in the complex z = exp{2βh} plane, where βis the inverse temperature and h the uniform magnetic field, touch the positive real axis, at least for large values of β. In some cases one obtains, in an appropriately taken βto infinity limit, a gas of hard objects on a set Λ'; the fugacity for the limiting system is a rescaling of z and the Lee-Yang zeros of the new partition function also avoid the positive real axis. For certain forms of the energies of the individual units the Lee-Yang zeros of both the finite- and zero-temperature systems lie on the negative real axis for all β. One zero-temperature limit of this type, for example, is a monomer-dimer system; our results thus generalize, to finite β, a well-known result of Heilmann and Lieb that the Lee-Yang zeros of monomer-dimer systems are real and negative.

preprint2011arXiv

Boltzmann's Entropy and Large Deviation Lyapunov Functionals for Closed and Open Macroscopic Systems

I give a brief overview of the resolution of the apparent problem of reconciling time symmetric microscopic dynamic with time asymmetric equations describing the evolution of macroscopic variables. I then show how the large deviation function of the stationary state of the microscopic system can be used as a Lyapunov function for the macroscopic evolution equations.

preprint2011arXiv

Checkerboards, stripes and corner energies in spin models with competing interactions

We study the zero temperature phase diagram of Ising spin systems in two dimensions in the presence of competing interactions, long range antiferromagnetic and nearest neighbor ferromagnetic of strength J. We first introduce the notion of a "corner energy" which shows, when the antiferromagnetic interaction decays faster than the fourth power of the distance, that a striped state is favored with respect to a checkerboard state when J is close to J_c, the transition to the ferromagnetic state, i.e., when the length scales of the uniformly magnetized domains become large. Next, we perform detailed analytic computations on the energies of the striped and checkerboard states in the cases of antiferromagnetic interactions with exponential decay and with power law decay r^{-p}, p>2, that depend on the Manhattan distance instead of the Euclidean distance. We prove that the striped phase is always favored compared to the checkerboard phase when the scale of the ground state structure is very large. This happens for J\lesssim J_c if p>3, and for J sufficiently large if 2<p<=3. Many of our considerations involving rigorous bounds carry over to dimensions greater than two and to more general short-range ferromagnetic interactions.

preprint2011arXiv

Necessary and sufficient conditions for realizability of point processes

We give necessary and sufficient conditions for a pair of (generalized) functions $ρ_1(\mathbf{r}_1)$ and $ρ_2(\mathbf{r}_1,\mathbf{r}_2)$, $\mathbf{r}_i\in X$, to be the density and pair correlations of some point process in a topological space $X$, for example, $\mathbb {R}^d$, $\mathbb {Z}^d$ or a subset of these. This is an infinite-dimensional version of the classical "truncated moment" problem. Standard techniques apply in the case in which there can be only a bounded number of points in any compact subset of $X$. Without this restriction we obtain, for compact $X$, strengthened conditions which are necessary and sufficient for the existence of a process satisfying a further requirement---the existence of a finite third order moment. We generalize the latter conditions in two distinct ways when $X$ is not compact.

preprint2011arXiv

Nonequilibrium Stationary States of Harmonic Chains with Bulk Noises

We consider a chain composed of $N$ coupled harmonic oscillators in contact with heat baths at temperature $T_\ell$ and $T_r$ at sites 1 and $N$ respectively. The oscillators are also subjected to non-momentum conserving bulk stochastic noises. These make the heat conductivity satisfy Fourier's law. Here we describe some new results about the hydrodynamical equations for typical macroscopic energy and displacement profiles, as well as their fluctuations and large deviations, in two simple models of this type.

preprint2011arXiv

Phase diagram of a generalized ABC model on the interval

We study the equilibrium phase diagram of a generalized ABC model on an interval of the one-dimensional lattice: each site $i=1,...,N$ is occupied by a particle of type $\a=A,B,C,$ with the average density of each particle species $N_\a/N=r_\a$ fixed. These particles interact via a mean field non-reflection-symmetric pair interaction. The interaction need not be invariant under cyclic permutation of the particle species as in the standard ABC model studied earlier. We prove in some cases and conjecture in others that the scaled infinite system $N\rw\infty$, $i/N\rw x\in[0,1]$ has a unique density profile $\p_\a(x)$ except for some special values of the $r_\a$ for which the system undergoes a second order phase transition from a uniform to a nonuniform periodic profile at a critical temperature $T_c=3\sqrt{r_A r_B r_C}/2π$.

preprint2011arXiv

Proof of Rounding by Quenched Disorder of First Order Transitions in Low-Dimensional Quantum Systems

We prove that for quantum lattice systems in d<=2 dimensions the addition of quenched disorder rounds any first order phase transition in the corresponding conjugate order parameter, both at positive temperatures and at T=0. For systems with continuous symmetry the statement extends up to d<=4 dimensions. This establishes for quantum systems the existence of the Imry-Ma phenomenon which for classical systems was proven by Aizenman and Wehr. The extension of the proof to quantum systems is achieved by carrying out the analysis at the level of thermodynamic quantities rather than equilibrium states.

preprint2011arXiv

Spatial Structure of Stationary Nonequilibrium States in the Thermostatted Periodic Lorentz Gas

We investigate analytically and numerically the spatial structure of the non-equilibrium stationary states (NESS) of a point particle moving in a two dimensional periodic Lorentz gas (Sinai Billiard). The particle is subject to a constant external electric field E as well as a Gaussian thermostat which keeps the speed |v| constant. We show that despite the singular nature of the SRB measure its projections on the space coordinates are absolutely continuous. We further show that these projections satisfy linear response laws for small E. Some of them are computed numerically. We compare these results with those obtained from simple models in which the collisions with the obstacles are replaced by random collisions.Similarities and differences are noted.

preprint2010arXiv

Heat conduction and phonon localization in disordered harmonic crystals

We investigate the steady state heat current in two and three dimensional isotopically disordered harmonic lattices. Using localization theory as well as kinetic theory we estimate the system size dependence of the current. These estimates are compared with numerical results obtained using an exact formula for the current given in terms of a phonon transmission function, as well as by direct nonequilibrium simulations. We find that heat conduction by high-frequency modes is suppressed by localization while low-frequency modes are strongly affected by boundary conditions. Our {\color{black}heuristic} arguments show that Fourier's law is valid in a three dimensional disordered solid except for special boundary conditions. We also study the pinned case relevant to localization in quantum systems and often used as a model system to study the validity of Fourier's law. Here we provide the first numerical verification of Fourier's law in three dimensions. In the two dimensional pinned case we find that localization of phonon modes leads to a heat insulator.

preprint2010arXiv

Long-Time Behavior of Macroscopic Quantum Systems: Commentary Accompanying the English Translation of John von Neumann's 1929 Article on the Quantum Ergodic Theorem

The renewed interest in the foundations of quantum statistical mechanics in recent years has led us to study John von Neumann's 1929 article on the quantum ergodic theorem. We have found this almost forgotten article, which until now has been available only in German, to be a treasure chest, and to be much misunderstood. In it, von Neumann studied the long-time behavior of macroscopic quantum systems. While one of the two theorems announced in his title, the one he calls the "quantum H-theorem", is actually a much weaker statement than Boltzmann's classical H-theorem, the other theorem, which he calls the "quantum ergodic theorem", is a beautiful and very non-trivial result. It expresses a fact we call "normal typicality" and can be summarized as follows: For a "typical" finite family of commuting macroscopic observables, every initial wave function $ψ_0$ from a micro-canonical energy shell so evolves that for most times $t$ in the long run, the joint probability distribution of these observables obtained from $ψ_t$ is close to their micro-canonical distribution.

preprint2010arXiv

Nonequilibrium stationary state of a truncated stochastic NLSE: I. Formulation and mean field approximation

We investigate the stationary state of a model system evolving according to a modified focusing truncated nonlinear Schrödinger equation (NLSE) used to describe the envelope of Langmuir waves in a plasma. We restrict the system to have a finite number of normal modes each of which is in contact with a Langevin heat bath at temperature $T$. Arbitrarily large realizations of the field are prevented by restricting each mode to a maximum amplitude. We consider a simple modeling of wave-breaking in which each mode is set equal to zero when it reaches its maximum amplitude. Without wave-breaking the stationary state is given by a Gibbs measure. With wave-breaking the system attains a nonequilibrium stationary state which is the unique invariant measure of the time evolution. A mean field analysis shows that the system exhibits a transition from a regime of low field values at small $|λ|$, to a regime of higher field values at large $|λ|$, where $λ<0$ specifies the strength of the nonlinearity in the focusing case. Field values at large $|λ|$ are significantly smaller with wave-breaking than without wave-breaking.

preprint2010arXiv

Normal Typicality and von Neumann's Quantum Ergodic Theorem

We discuss the content and significance of John von Neumann's quantum ergodic theorem (QET) of 1929, a strong result arising from the mere mathematical structure of quantum mechanics. The QET is a precise formulation of what we call normal typicality, i.e., the statement that, for typical large systems, every initial wave function $ψ_0$ from an energy shell is "normal": it evolves in such a way that $|ψ_t> <ψ_t|$ is, for most $t$, macroscopically equivalent to the micro-canonical density matrix. The QET has been mostly forgotten after it was criticized as a dynamically vacuous statement in several papers in the 1950s. However, we point out that this criticism does not apply to the actual QET, a correct statement of which does not appear in these papers, but to a different (indeed weaker) statement. Furthermore, we formulate a stronger statement of normal typicality, based on the observation that the bound on the deviations from the average specified by von Neumann is unnecessarily coarse and a much tighter (and more relevant) bound actually follows from his proof.

preprint2010arXiv

On Some Classes of Open Two-Species Exclusion Processes

We investigate some properties of the nonequilibrium stationary state (NESS) of a one dimensional open system consisting of first and second class (type 1 and type 2) particles. The dynamics are totally asymmetric but the rates for the different permitted exchanges (10 -> 01, 12 -> 21, and 20 -> 02) need not be equal. The entrance and exit rates of the different species can also be different. We show that for certain classes of rates one can compute the currents and phase diagram, or at least obtain some monotonicity properties. For other classes one can obtain a matrix representation of the NESS; this generalizes previous work in which second class particles can neither enter nor leave the system. We analyze a simple example of this type and establish the existence of a randomly located shock at which the typical density profiles of all three species are discontinuous.

preprint2010arXiv

Phase transitions with four-spin interactions

Using an extended Lee-Yang theorem and GKS correlation inequalities, we prove, for a class of ferromagnetic multi-spin interactions, that they will have a phase transition(and spontaneous magnetization) if, and only if, the external field $h=0$ (and the temperature is low enough). We also show the absence of phase transitions for some nonferromagnetic interactions. The FKG inequalities are shown to hold for a larger class of multi-spin interactions.

preprint2010arXiv

Restricted equilibrium ensembles: Exact equation of state of a model glass

We investigate the thermodynamic properties of a toy model of glasses: a hard-core lattice gas with nearest neighbor interaction in one dimension. The time-evolution is Markovian, with nearest-neighbor and next-nearest neighbor hoppings, and the transition rates are assumed to satisfy detailed balance condition, but the system is non-ergodic below a glass temperature. Below this temperature, the system is in restricted thermal equilibrium, where both the number of sectors, and the number of accessible states within a sector grow exponentially with the size of the system. Using partition functions that sum only over dynamically accessible states within a sector, and then taking a quenched average over the sectors, we determine the exact equation of state of this system.

preprint2010arXiv

The grand canonical ABC model: a reflection asymmetric mean field Potts model

We investigate the phase diagram of a three-component system of particles on a one-dimensional filled lattice, or equivalently of a one-dimensional three-state Potts model, with reflection asymmetric mean field interactions. The three types of particles are designated as $A$, $B$, and $C$. The system is described by a grand canonical ensemble with temperature $T$ and chemical potentials $Tλ_A$, $Tλ_B$, and $Tλ_C$. We find that for $λ_A=λ_B=λ_C$ the system undergoes a phase transition from a uniform density to a continuum of phases at a critical temperature $\hat T_c=(2π/\sqrt3)^{-1}$. For other values of the chemical potentials the system has a unique equilibrium state. As is the case for the canonical ensemble for this $ABC$ model, the grand canonical ensemble is the stationary measure satisfying detailed balance for a natural dynamics. We note that $\hat T_c=3T_c$, where $T_c$ is the critical temperature for a similar transition in the canonical ensemble at fixed equal densities $r_A=r_B=r_C=1/3$.

preprint2009arXiv

Correlation Inequalities for Quantum Spin Systems with Quenched Centered Disorder

It is shown that random quantum spin systems with centered disorder satisfy correlation inequalities previously proved (arXiv:cond-mat/0612371) in the classical case. Consequences include monotone approach of pressure and ground state energy to the thermodynamic limit. Signs and bounds on the surface pressures for different boundary conditions are also derived for finite range potentials.

preprint2009arXiv

Heat transport and phonon localization in mass-disordered harmonic crystals

We investigate the steady state heat current in two and three dimensional disordered harmonic crystals in a slab geometry, connected at the boundaries to stochastic white noise heat baths at different temperatures.The disorder causes short wavelength phonon modes to be localized so the heat current in this system is carried by the extended phonon modes which can be either diffusive or ballistic. Using ideas both from localization theory and from kinetic theory we estimate the contribution of various modes to the heat current and from this we obtain the asymptotic system size dependence of the current. These estimates are compared with results obtained from a numerical evaluation of an exact formula for the current, given in terms of a frequency transmission function, as well as from direct nonequilibrium simulations. These yield a strong dependence of the heat flux on boundary conditions. Our analytical arguments show that for realistic boundary conditions the conductivity is finite in three dimensions but we are not able to verify this numerically, except in the case where the system is subjected to an external pinning potential. This case is closely related to the problem of localization of electrons in a random potential and here we numerically verify that the pinned three dimensional system satisfies Fourier's law while the two dimensional system is a heat insulator. We also investigate the inverse participation ratio of different normal modes.

preprint2009arXiv

Modulated phases of a 1D sharp interface model in a magnetic field

We investigate the ground states of 1D continuum models having short-range ferromagnetic type interactions and a wide class of competing longer-range antiferromagnetic type interactions. The model is defined in terms of an energy functional, which can be thought of as the Hamiltonian of a coarse-grained microscopic system or as a mesoscopic free energy functional describing various materials. We prove that the ground state is simple periodic whatever the prescribed total magnetization might be. Previous studies of this model of frustrated systems assumed this simple periodicity but, as in many examples in condensed matter physics, it is neither obvious nor always true that ground states do not have a more complicated, or even chaotic structure.

preprint2009arXiv

Nonequilibrium Stationary Solutions of Thermostated Boltzmann Equation in a Field

We consider a system of particles subjected to a uniform external force E and undergoing random collisions with "virtual" fixed obstacles, as in the Drude model of conductivity. The system is maintained in a nonequilibrium stationary state by a Gaussian thermostat. In a suitable limit the system is described by a self consistent Boltzmann equation for the one particle distribution function f. We find that after a long time f(v,t) approaches a stationary velocity distribution f(v) which vanishes for large speeds, i.e. f(v)=0 for |v|>vmax(E), with vmax(E)~1/|E| as |E| -> 0. In that limit f(v)~exp(-c|v|^3) for fixed v, where c depends on mean free path of the particle. f(v) is computed explicitly in one dimension.

preprint2009arXiv

On Spin Systems with Quenched Randomness: Classical and Quantum

The rounding of first order phase transitions by quenched randomness is stated in a form which is applicable to both classical and quantum systems: The free energy, as well as the ground state energy, of a spin system on a $d$-dimensional lattice is continuously differentiable with respect to any parameter in the Hamiltonian to which some randomness has been added when $d \leq 2$. This implies absence of jumps in the associated order parameter, e.g., the magnetization in case of a random magnetic field. A similar result applies in cases of continuous symmetry breaking for $d \leq 4$. Some questions concerning the behavior of related order parameters in such random systems are discussed.

preprint2009arXiv

On the Approach to Thermal Equilibrium of Macroscopic Quantum Systems

We consider an isolated, macroscopic quantum system. Let H be a micro-canonical "energy shell," i.e., a subspace of the system's Hilbert space spanned by the (finitely) many energy eigenstates with energies between E and E + delta E. The thermal equilibrium macro-state at energy E corresponds to a subspace H_{eq} of H such that dim H_{eq}/dim H is close to 1. We say that a system with state vector psi in H is in thermal equilibrium if psi is "close" to H_{eq}. We show that for "typical" Hamiltonians with given eigenvalues, all initial state vectors psi_0 evolve in such a way that psi_t is in thermal equilibrium for most times t. This result is closely related to von Neumann's quantum ergodic theorem of 1929.

preprint2009arXiv

Rounding of First Order Transitions in Low-Dimensional Quantum Systems with Quenched Disorder

We prove that the addition of an arbitrarily small random perturbation of a suitable type to a quantum spin system rounds a first order phase transition in the conjugate order parameter in d <= 2 dimensions, or in systems with continuous symmetry in d <= 4. This establishes rigorously for quantum systems the existence of the Imry-Ma phenomenon, which for classical systems was proven by Aizenman and Wehr.

preprint2009arXiv

Space Charge Effects in Field Emission: One Dimensional Theory

The current associated with field emission is greatly dependent on the electric field at the emitting electrode. This field is a combination of the electric field in vacuum and the space charge created by the current. The latter becomes more important as the currentdensity increases. Here, a study is performed using a modified classical 1D Child-Langmuir description that allows for exact solutions in order to characterize the contributions due to space charge. Methods to connect the 1D approach to an array of periodic 3D structures are considered.

preprint2008arXiv

On the Two Species Asymmetric Exclusion Process with Semi-Permeable Boundaries

We investigate the structure of the nonequilibrium stationary state (NESS) of a system of first and second class particles, as well as vacancies (holes), on L sites of a one-dimensional lattice in contact with first class particle reservoirs at the boundary sites; these particles can enter at site 1, when it is vacant, with rate alpha, and exit from site L with rate beta. Second class particles can neither enter nor leave the system, so the boundaries are semi-permeable. The internal dynamics are described by the usual totally asymmetric exclusion process (TASEP) with second class particles. An exact solution of the NESS was found by Arita. Here we describe two consequences of the fact that the flux of second class particles is zero. First, there exist (pinned and unpinned) fat shocks which determine the general structure of the phase diagram and of the local measures; the latter describe the microscopic structure of the system at different macroscopic points (in the limit L going to infinity in terms of superpositions of extremal measures of the infinite system. Second, the distribution of second class particles is given by an equilibrium ensemble in fixed volume, or equivalently but more simply by a pressure ensemble, in which the pair potential between neighboring particles grows logarithmically with distance. We also point out an unexpected feature in the microscopic structure of the NESS for finite L: if there are n second class particles in the system then the distribution of first class particles (respectively holes) on the first (respectively last) n sites is exchangeable.

preprint2007arXiv

Periodic minimizers in 1D local mean field theory

Using reflection positivity techniques we prove the existence of minimizers for a class of mesoscopic free-energies representing 1D systems with competing interactions. All minimizers are either periodic, with zero average, or of constant sign. If the local term in the free energy satisfies a convexity condition, then all minimizers are either periodic or constant. Examples of both phenomena are given. This extends our previous work where such results were proved for the ground states of lattice systems with ferromagnetic nearest neighbor interactions and dipolar type antiferromagnetic long range interactions.

preprint2007arXiv

Striped phases in two dimensional dipole systems

We prove that a system of discrete 2D in-plane dipoles with four possible orientations, interacting via a 3D dipole-dipole interaction plus a nearest neighbor ferromagnetic term, has periodic striped ground states. As the strength of the ferromagnetic term is increased, the size of the stripes in the ground state increases, becoming infinite, i.e., giving a ferromagentic ground state, when the ferromagentic interaction exceeds a certain critical value. We also give a rigorous proof of the reorientation transition in the ground state of a 2D system of discrete dipoles with six possible orientations, interacting via a 3D dipole-dipole interaction plus a nearest neighbor antiferromagnetic term. As the strength of the antiferromagnetic term is increased the ground state flips from being striped and in-plane to being staggered and out-of-plane. An example of a rotator model with a sinusoidal ground state is also discussed.

preprint2006arXiv

Ising models with long-range dipolar and short-range ferromagnetic interactions

We study the ground state of a $d$--dimensional Ising model with both long range (dipole--like) and nearest neighbor ferromagnetic (FM) interactions. The long range interaction is equal to $r^{-p}$, $p>d$, while the FM interaction has strength $J$. If $p>d+1$ and $J$ is large enough the ground state is FM, while if $d<p\le d+1$ the FM state is not the ground state for any choice of $J$. In $d=1$ we show that for any $p>1$ the ground state has a series of transitions from an antiferromagnetic state of period 2 to $2h$--periodic states of blocks of sizes $h$ with alternating sign, the size $h$ growing when the FM interaction strength $J$ is increased (a generalization of this result to the case $0<p\le 1$ is also discussed). In $d\ge 2$ we prove, for $d<p\le d+1$, that the dominant asymptotic behavior of the ground state energy agrees for large $J$ with that obtained from a periodic striped state conjectured to be the true ground state. The geometry of contours in the ground state is discussed.

preprint2002arXiv

Hydrodynamics of binary fluid phase segregation

Starting with the Vlasov-Boltzmann equation for a binary fluid mixture, we derive an equation for the velocity field $\bm{u}$ when the system is segregated into two phases (at low temperatures) with a sharp interface between them. $\bm{u}$ satisfies the incompressible Navier-Stokes equations together with a jump boundary condition for the pressure across the interface which, in turn, moves with a velocity given by the normal component of $\bm{u} $. Numerical simulations of the Vlasov-Boltzmann equations for shear flows parallel and perpendicular to the interface in a phase segregated mixture support this analysis. We expect similar behavior in real fluid mixtures.

preprint1997arXiv

Solutions to the Boltzmann equation in the Boussinesq regime

We consider a gas in a horizontal slab, in which the top and bottom walls are kept at different temperatures. The system is described by the Boltzmann equation (BE) with Maxwellian boundary conditions specifying the wall temperatures. We study the behavior of the system when the Knudsen number $ε$ is small and the temperature difference between the walls as well as the velocity field is of order $ε$, while the gravitational force is of order $ε^2$. We prove that there exists a solution to the BE for which is near a global Maxwellian, and whose moments are close, up to order $ε^2$ to the density, velocity and temperature obtained from the smooth solution of the Oberbeck-Boussinesq equations, up to the time this one stays regular.

preprint1997arXiv

Stationary Nonequilibrium States in Boundary Driven Hamiltonian Systems: Shear Flow

We investigate stationary nonequilibrium states of systems of particles moving according to Hamiltonian dynamics with specified potentials. The systems are driven away from equilibrium by Maxwell demon ``reflection rules'' at the walls. These deterministic rules conserve energy but not phase space volume, and the resulting global dynamics may or may not be time reversible (or even invertible). Using rules designed to simulate moving walls we can obtain a stationary shear flow. Assuming that for macroscopic systems this flow satisfies the Navier-Stokes equations, we compare the hydrodynamic entropy production with the average rate of phase space volume compression. We find that they are equal when the velocity distribution of particles incident on the walls is a local Maxwellian. An argument for a general equality of this kind, based on the assumption of local thermodynamic equilibrium, is given. Molecular dynamic simulations of hard disks in a channel produce a steady shear flow with the predicted behavior.

preprint1993arXiv

Exact Results for the Asymmetric Simple Exclusion Process with a Blockage

We present new results for the current as a function of transmission rate in the one dimensional totally asymmetric simple exclusion process (TASEP) with a blockage that lowers the jump rate at one site from one to r < 1. Exact finite volume results serve to bound the allowed values for the current in the infinite system. This proves the existence of a gap in allowed density corresponding to a nonequilibrium ``phase transition'' in the infinite system. A series expansion in r, derived from the finite systems, is proven to be asymptotic for all sufficiently large systems. Pade approximants based on this series, which make specific assumptions about the nature of the singularity at r = 1, match numerical data for the ``infinite'' system to a part in 10^4.