Discreteness-Induced Slow Relaxation in Reversible Catalytic Reaction Networks

Slowing down of the relaxation of the fluctuations around equilibrium is investigated both by stochastic simulations and by analysis of Master equation of reversible reaction networks consisting of resources and the corresponding products that work as catalysts. As the number of molecules $N$ is decreased, the relaxation time to equilibrium is prolonged due to the deficiency of catalysts, as demonstrated by the amplification compared to that by the continuum limit. This amplification ratio of the relaxation time is represented by a scaling function as $h = N \exp(-\beta V)$, and it becomes prominent as $N$ becomes less than a critical value $h \sim 1$, where $\beta$ is the inverse temperature and $V$ is the energy gap between a product and a resource

Bethe--Salpeter equation in QCD

We extend to regular QCD the derivation of a confining $q \bar{q}$ Bethe--Salpeter equation previously given for the simplest model of scalar QCD in which quarks are treated as spinless particles. We start from the same assumptions on the Wilson loop integral already adopted in the derivation of a semirelativistic heavy quark potential. We show that, by standard approximations, an effective meson squared mass operator can be obtained from our BS kernel and that, from this, by ${1\over m^2}$ expansion the corresponding Wilson loop potential can be reobtained, spin--dependent and velocity--dependent terms included. We also show that, on the contrary, neglecting spin--dependent terms, relativistic flux tube model is reproduced.Comment: 23 pages, revte

Two phase transitions in the fully frustrated XYXY model

The fully frustrated $XY$ model on a square lattice is studied by means of Monte Carlo simulations. A Kosterlitz-Thouless transition is found at $T_{\rm KT} \approx 0.446$, followed by an ordinary Ising transition at a slightly higher temperature, $T_c \approx 0.452$. The non-Ising exponents reported by others, are explained as a failure of finite size scaling due to the screening length associated with the nearby Kosterlitz-Thouless transition.Comment: REVTEX file, 8 pages, 5 figures in uuencoded postscrip

Constraint-preserving boundary conditions in the 3+1 first-order approach

A set of energy-momentum constraint-preserving boundary conditions is proposed for the first-order Z4 case. The stability of a simple numerical implementation is tested in the linear regime (robust stability test), both with the standard corner and vertex treatment and with a modified finite-differences stencil for boundary points which avoids corners and vertices even in cartesian-like grids. Moreover, the proposed boundary conditions are tested in a strong field scenario, the Gowdy waves metric, showing the expected rate of convergence. The accumulated amount of energy-momentum constraint violations is similar or even smaller than the one generated by either periodic or reflection conditions, which are exact in the Gowdy waves case. As a side theoretical result, a new symmetrizer is explicitly given, which extends the parametric domain of symmetric hyperbolicity for the Z4 formalism. The application of these results to first-order BSSN-like formalisms is also considered.Comment: Revised version, with conclusive numerical evidence. 23 pages, 12 figure

Biogeography and Long-Run Economic Development

The transition from a hunter-gather economy to agricultural production, which made possible the endogenous technological progress that ultimately led to the industrial revolution, is one of the most important events in the thousands of years of humankind’s economic development. In this paper we present theory and evidence showing that exogenous geography and initial condition biogeography exerted decisive influence on the location and timing of transitions to sedentary agriculture, to complex social organization and, eventually, to modern industrial production. Evidence from a large cross-section of countries indicates that the effects of geographic and biogeographic endowments on contemporary levels of economic development are remarkably strong.Geography biogeography and growth; Economic development; Agricultural revolution; Institutions and growth; Plants animals and growth

A hybrid simulation model for a stable auroral arc

International audienceWe present a new type of hybrid simulation model, intended to simulate a single stable auroral arc in the latitude/altitude plane. The ionospheric ions are treated as particles, the electrons are assumed to follow a Boltzmann response and the magnetospheric ions are assumed to be so hot that they form a background population unaffected by the electric fields that arise. The system is driven by assumed parallel electron energisation causing a primary negative charge cloud and an associated potential structure to build up. The results show how a closed potential structure and density depletion of an auroral arc build up and how they decay after the driver is turned off. The model also produces upgoing energetic ion beams and predicts strong static perpendicular electric fields to be found in a relatively narrow altitude range (~ 5000?11 000 km)

The current-voltage relationship revisited: exact and approximate formulas with almost general validity for hot magnetospheric electrons for bi-Maxwellian and kappa distributions

International audienceWe derive the current-voltage relationship in the auroral region taking into account magnetospheric electrons for the bi-Maxwellian and kappa source plasma distribution functions. The current-voltage formulas have in principle been well known for a long time, but the kappa energy flux formulas have not appeared in the literature before. We give a unified treatment of the bi-Maxwellian and kappa distributions, correcting some errors in previous work. We give both exact results and two kinds of approximate formulas for the current density and the energy flux. The first approximation is almost generally valid and is practical to compute. The first approximation formulas are therefore suitable for use in simulations. In the second approximation we assume in addition that the thermal energy is small compared to the potential drop. This yields even simpler linear formulas which are suitable for many types of event studies and which have a more transparent physical interpretation than the first approximation formulas. We also show how it is possible to derive the first approximation formulas even for those distributions for which the exact results can not be computed analytically. The kappa field-aligned conductance value turns out always to be smaller than the corresponding Maxwellian conductance. We also verify that the obtained kappa current density and energy flux formulas go to Maxwellian results when ???

Vortex glass transitions in disordered three-dimensional XY models: Simulations for several different sets of parameters

The anisotropic frustrated 3D XY model with strong disorder in the coupling constants is studied as a model of a disordered superconductor in an applied magnetic field. Simulations with the exchange Monte Carlo method are performed for frustrations f=1/5 and f=1/4, corresponding to two different values of the magnetic field along the z direction. The anisotropy is also varied. The determination of the helicity modulus from twist histograms is discussed in some detail and the helicity modulus is used in finite size scaling analyses of the vortex glass transition. The general picture is that the behavior in [Phys. Rev. Lett. 91, 077002 (2003)] is confirmed. For strong (e.g. isotropic) coupling in the z direction the helicity modulus fails to scale and it is argued that this is due to a too small effective randomness of such systems for the accessible system sizes