A Tutorial on Advanced Dynamic Monte Carlo Methods for Systems with Discrete State Spaces

Advanced algorithms are necessary to obtain faster-than-real-time dynamic simulations in a number of different physical problems that are characterized by widely disparate time scales. Recent advanced dynamic Monte Carlo algorithms that preserve the dynamics of the model are described. These include the $n$-fold way algorithm, the Monte Carlo with Absorbing Markov Chains (MCAMC) algorithm, and the Projective Dynamics (PD) algorithm. To demonstrate the use of these algorithms, they are applied to some simplified models of dynamic physical systems. The models studied include a model for ion motion through a pore such as a biological ion channel and the metastable decay of the ferromagnetic Ising model. Non-trivial parallelization issues for these dynamic algorithms, which are in the class of parallel discrete event simulations, are discussed. Efforts are made to keep the article at an elementary level by concentrating on a simple model in each case that illustrates the use of the advanced dynamic Monte Carlo algorithm.Comment: 53 pages, 17 figure

Constraints on low energy QCD parameters from η→3π\eta \to 3\pi and ππ\pi\pi scattering

The $\eta \to 3\pi$ decays are a valuable source of information on low energy QCD. Yet they were not used for an extraction of the three flavor chiral symmetry breaking order parameters until now. We use a Bayesian approach in the framework of resummed chiral perturbation theory to obtain constraints on the quark condensate and pseudoscalar decay constant in the chiral limit. We compare our results with recent CHPT and lattice QCD fits and find some tension, as the $\eta \to 3\pi$ data seem to prefer a larger ratio of the chiral order parameters. The results also disfavor a very large value of the pseudoscalar decay constant in the chiral limit, which was found by some recent works. In addition, we present results of a combined analysis including $\eta \to 3\pi$ decays and $\pi\pi$ scattering and though the picture does not changed appreciably, we find some tension between the data we use. We also try to extract information on the mass difference of the light quarks, but the uncertainties prove to be large.Comment: 23 pages, 8 figure

Weak solutions for some compressible multicomponent fluid models

The principle purpose of this work is to investigate a "viscous" version of a "simple" but still realistic bi-fluid model described in [Bresch, Desjardin, Ghidaglia, Grenier, Hillairet] whose "non-viscous" version is derived from physical considerations in \cite[Ishii, Hibiki]{ISHI} as a particular sample of a multifluid model with algebraic closure. The goal is to show existence of weak solutions for large initial data on an arbitrarily large time interval. We achieve this goal by transforming the model to an academic system which resembles to the compressible Navier-Stokes equations, with however two continuity equations and a momentum equation endowed with pressure of complicated structure dependent on two variable densities. The new "academic system" is then solved by an adaptation of the Lions--Feireisl approach for solving compressible Navier--Stokes equation, completed with several observations related to the DiPerna--Lions transport theory inspired by [Maltese, Michalek, Mucha, Novotny, Pokorny, Zatorska] and [Vasseur, Wen, Yu]. We also explain how these techniques can be generalized to a model of mixtures with more then two species. This is the first result on the existence of weak solutions for any realistic multifluid system