Astrophysical Fluid Dynamics via Direct Statistical Simulation

In this paper we introduce the concept of Direct Statistical Simulation (DSS) for astrophysical flows. This technique may be appropriate for problems in astrophysical fluids where the instantaneous dynamics of the flows are of secondary importance to their statistical properties. We give examples of such problems including mixing and transport in planets, stars and disks. The method is described for a general set of evolution equations, before we consider the specific case of a spectral method optimised for problems on a spherical surface. The method is illustrated for the simplest non-trivial example of hydrodynamics and MHD on a rotating spherical surface. We then discuss possible extensions of the method both in terms of computational methods and the range of astrophysical problems that are of interest.Comment: 26 pages, 11 figures, added clarifying remarks and references, and corrected typos. This version is accepted for publication in The Astrophysical Journa

New Keynesian versus old Keynesian government spending multipliers

Renewed interest in fiscal policy has increased the use of quantitative models to evaluate policy. Because of modeling uncertainty, it is essential that policy evaluations be robust to alternative assumptions. We find that models currently being used in practice to evaluate fiscal policy stimulus proposals are not robust. Government spending multipliers in an alternative empirically-estimated and widely-cited new Keynesian model are much smaller than in these old Keynesian models; the estimated stimulus is extremely small with GDP and employment effects only one-sixth as large

Symplectic algorithm for constant-pressure molecular dynamics using a Nose-Poincare thermostat

We present a new algorithm for isothermal-isobaric molecular-dynamics simulation. The method uses an extended Hamiltonian with an Andersen piston combined with the Nos'e-Poincar'e thermostat, recently developed by Bond, Leimkuhler and Laird [J. Comp. Phys., 151, (1999)]. This Nos'e-Poincar'e-Andersen (NPA) formulation has advantages over the Nos'e-Hoover-Andersen approach in that the NPA is Hamiltonian and can take advantage of symplectic integration schemes, which lead to enhanced stability for long-time simulations. The equations of motion are integrated using a Generalized Leapfrog Algorithm and the method is easy to implement, symplectic, explicit and time reversible. To demonstrate the stability of the method we show results for test simulations using a model for aluminum.Comment: 7 page

Dissimilar bouncy walkers

We consider the dynamics of a one-dimensional system consisting of dissimilar hardcore interacting (bouncy) random walkers. The walkers' (diffusing particles') friction constants xi_n, where n labels different bouncy walkers, are drawn from a distribution rho(xi_n). We provide an approximate analytic solution to this recent single-file problem by combining harmonization and effective medium techniques. Two classes of systems are identified: when rho(xi_n) is heavy-tailed, rho(xi_n)=A xi_n^(-1-\alpha) (0<alpha<1) for large xi_n, we identify a new universality class in which density relaxations, characterized by the dynamic structure factor S(Q,t), follows a Mittag-Leffler relaxation, and the the mean square displacement of a tracer particle (MSD) grows as t^delta with time t, where delta=alpha/(1+\alpha). If instead rho is light-tailedsuch that the mean friction constant exist, S(Q,t) decays exponentially and the MSD scales as t^(1/2). We also derive tracer particle force response relations. All results are corroborated by simulations and explained in a simplified model.Comment: 11 pages, to appear in Journal of Chemical Physic

Optimizing evacuation flow in a two-channel exclusion process

We use a basic setup of two coupled exclusion processes to model a stylised situation in evacuation dynamics, in which evacuees have to choose between two escape routes. The coupling between the two processes occurs through one common point at which particles are injected, the process can be controlled by directing incoming individuals into either of the two escape routes. Based on a mean-field approach we determine the phase behaviour of the model, and analytically compute optimal control strategies, maximising the total current through the system. Results are confirmed by numerical simulations. We also show that dynamic intervention, exploiting fluctuations about the mean-field stationary state, can lead to a further increase in total current.Comment: 16 pages, 6 figure

General Monogamy Inequality for Bipartite Qubit Entanglement

We consider multipartite states of qubits and prove that their bipartite quantum entanglement, as quantified by the concurrence, satisfies a monogamy inequality conjectured by Coffman, Kundu, and Wootters. We relate this monogamy inequality to the concept of frustration of correlations in quantum spin systems.Comment: Fixed spelling mistake. Added references. Fixed error in transformation law. Shorter and more explicit proof of capacity formula. Reference added. Rewritten introduction and conclusion

Achievable Qubit Rates for Quantum Information Wires

Suppose Alice and Bob have access to two separated regions, respectively, of a system of electrons moving in the presence of a regular one-dimensional lattice of binding atoms. We consider the problem of communicating as much quantum information, as measured by the qubit rate, through this quantum information wire as possible. We describe a protocol whereby Alice and Bob can achieve a qubit rate for these systems which is proportional to N^(-1/3) qubits per unit time, where N is the number of lattice sites. Our protocol also functions equally in the presence of interactions modelled via the t-J and Hubbard models

A heuristic quantum theory of the integer quantum Hall effect

Contrary to common belief, the current emitted by a contact embedded in a two-dimensional electron gas (2DEG) is quantized in the presence of electric and magnetic fields. This observation suggests a simple, clearly defined model for the quantum current through a Hall device that does not invoke disorder or interactions as the cause of the integer quantum Hall effect (QHE), but is based on a proper quantization of the classical electron drift motion. The theory yields a quantitative description of the breakdown of the QHE at high current densities that is in agreement with experimental data. Furthermore, several of its key points are in line with recent findings of experiments that address the dependency of the QHE on the 2DEG bias voltage, results that are not easily explained within the framework of conventional QHE models.Comment: 20 pages, 6 figure