Monte Carlo simulations of 4d simplicial quantum gravity

Dynamical triangulations of four-dimensional Euclidean quantum gravity give rise to an interesting, numerically accessible model of quantum gravity. We give a simple introduction to the model and discuss two particularly important issues. One is that contrary to recent claims there is strong analytical and numerical evidence for the existence of an exponential bound that makes the partition function well-defined. The other is that there may be an ambiguity in the choice of the measure of the discrete model which could even lead to the existence of different universality classes.Comment: 16 pages, LaTeX, epsf, 4 uuencoded figures; contribution to the JMP special issue on "Quantum Geometry and Diffeomorphism-Invariant Quantum Field Theory

Interdisciplinary Monte Carlo Simulations

Biological, linguistic, sociological and economical applications of statistical physics are reviewed here. They have been made on a variety of computers over a dozen years, not only at the NIC computers. A longer description can be found in our new book, an emphasis on teaching in Eur.J.Phys. 26, S 79 and AIP Conf. Proc. 779, 49, 56, 69 and 75.Comment: 11 pages including many Figs.; for 3rd NIC Symposium, Julich, 3/0

Effect of atomic scale plasticity on hydrogen diffusion in iron: Quantum mechanically informed and on-the-fly kinetic Monte Carlo simulations

We present an off-lattice, on-the-fly kinetic Monte Carlo (KMC) model for simulating stress-assisted diffusion and trapping of hydrogen by crystalline defects in iron. Given an embedded atom (EAM) potential as input, energy barriers for diffusion are ascertained on the fly from the local environments of H atoms. To reduce computational cost, on-the-fly calculations are supplemented with precomputed strain-dependent energy barriers in defect-free parts of the crystal. These precomputed barriers, obtained with high-accuracy density functional theory calculations, are used to ascertain the veracity of the EAM barriers and correct them when necessary. Examples of bulk diffusion in crystals containing a screw dipole and vacancies are presented. Effective diffusivities obtained from KMC simulations are found to be in good agreement with theory. Our model provides an avenue for simulating the interaction of hydrogen with cracks, dislocations, grain boundaries, and other lattice defects, over extended time scales, albeit at atomistic length scales

Reaction-controlled diffusion: Monte Carlo simulations

We study the coupled two-species non-equilibrium reaction-controlled diffusion model introduced by Trimper et al. [Phys. Rev. E 62, 6071 (2000)] by means of detailed Monte Carlo simulations in one and two dimensions. Particles of type A may independently hop to an adjacent lattice site provided it is occupied by at least one B particle. The B particle species undergoes diffusion-limited reactions. In an active state with nonzero, essentially homogeneous B particle saturation density, the A species displays normal diffusion. In an inactive, absorbing phase with exponentially decaying B density, the A particles become localized. In situations with algebraic decay rho_B(t) ~ t^{-alpha_B}, as occuring either at a non-equilibrium continuous phase transition separating active and absorbing states, or in a power-law inactive phase, the A particles propagate subdiffusively with mean-square displacement ~ t^{1-alpha_A}. We find that within the accuracy of our simulation data, \alpha_A = \alpha_B as predicted by a simple mean-field approach. This remains true even in the presence of strong spatio-temporal fluctuations of the B density. However, in contrast with the mean-field results, our data yield a distinctly non-Gaussian A particle displacement distribution n_A(x,t) that obeys dynamic scaling and looks remarkably similar for the different processes investigated here. Fluctuations of effective diffusion rates cause a marked enhancement of n_A(x,t) at low displacements |x|, indicating a considerable fraction of practically localized A particles, as well as at large traversed distances.Comment: Revtex, 19 pages, 27 eps figures include

Monte carlo simulations of parapatric speciation

Parapatric speciation is studied using an individual--based model with sexual reproduction. We combine the theory of mutation accumulation for biological ageing with an environmental selection pressure that varies according to the individuals geographical positions and phenotypic traits. Fluctuations and genetic diversity of large populations are crucial ingredients to model the features of evolutionary branching and are intrinsic properties of the model. Its implementation on a spatial lattice gives interesting insights into the population dynamics of speciation on a geographical landscape and the disruptive selection that leads to the divergence of phenotypes. Our results suggest that assortative mating is not an obligatory ingredient to obtain speciation in large populations at low gene flow.Comment: submitted to Phys.Rev.

MONTE CARLO SIMULATIONS OF MUON PRODUCTION

Muon production requirements for a muon collider are presented. Production of muons from pion decay is studied. Lithium lenses and solenoids are considered for focussing pions from a target, and for matching the pions into a decay channel. Pion decay channels of alternating quadrupoles and long solenoids are compared. Monte Carlo simulations are presented for production of $\pi \rightarrow \mu$ by protons over a wide energy range, and criteria for choosing the best proton energy are discussed.Comment: Latex uses mu95.sty, 19 pages, 5 postscript figures. A postscript file can be seen at URL http://www.cap.bnl.gov/~cap/mumu/important.html Search for Publication