Related papers: A new Monte Carlo method for dynamical evolution o…
In this work a replica exchange Monte Carlo scheme which considers an extended isobaric-isothermal ensemble with respect to pressure is applied to study hard spheres (HS). The idea behind the proposal is expanding volume instead of…
These lectures will cover methods for studying the evolution of galaxies since their formation. Because the properties of a galaxy depend on its history, an understanding of galaxy evolution requires that we understand the dynamical…
We present a set of new numerical methods that are relevant to calculating radiation pressure terms in hydrodynamics calculations, with a particular focus on massive star formation. The radiation force is determined from a Monte Carlo…
We study an induced dynamics in the space of energy of single-spin-flip Monte Carlo algorithm. The method gives an efficient reweighting technique. This dynamics is shown to have relaxation times proportional to the specific heat. Thus, it…
This paper, and its companion, investigate the evolution of dense stellar systems due to the influence of two-body gravitational encounters, physical collisions and stellar evolution. Our goal is the simulation of the densest centers of…
We present a radiation hydrodynamics simulation of the formation of a massive star using a Monte Carlo treatment for the radiation field. We find that strong, high speed bipolar cavities are driven by the radiation from the protostar, and…
We study nonequilibrium dynamical properties of inhomogeneous systems, in particular at a free surface or at a defect plane. Thereby we consider nonconserved (model-A) dynamics of a system which is prepared in the high-temperature phase and…
Many problems in materials science and biology involve particles interacting with strong, short-ranged bonds, that can break and form on experimental timescales. Treating such bonds as constraints can significantly speed up sampling their…
We have coupled a fast, parametrized star cluster evolution code to a Markov Chain Monte Carlo code to determine the distribution of probable initial conditions of observed star clusters, which may serve as a starting point for future…
The numerically exact path integral Monte Carlo approach for the real-time evolution of dissipative quantum systems (PIMC), particularly suited for systems with discrete configuration space (tight-binding systems), is extended to treat…
As a planet transits the face of a star, it accelerates along the line-of-sight. The changing delay in the propagation of photons produces an apparent deceleration of the planet across the sky throughout the transit. This persistent…
We review what has been learned recently using N-body simulations about the evolution of globular clusters. While simulations of star clusters have become more realistic, and now include the evolution of single and binary stars, the…
We study the long-term dynamics of a planetary system composed of a star and a planet. Both bodies are considered as extended, non-spherical, rotating objects. There are no assumptions made on the relative angles between the orbital angular…
A real-time path integral Monte Carlo approach is developed to study the dynamics in a many-body quantum system until reaching a nonequilibrium stationary state. The approach is based on augmenting an exact reduced equation for the…
The density relaxation phenomenon is modeled using both Monte Carlo and dissipative MD simulations to investigate the effects of regular taps applied to a vessel having a planar floor filled with monodisperse spheres. Results suggest the…
We present a series of simple, largely analytical models to compute the effects of disruption on the mass function of star clusters. Our calculations include evaporation by two-body relaxation and gravitational shocks and mass loss by…
A continuous-time formulation of the Diffusion Monte Carlo method for lattice models is presented. In its simplest version, without the explicit use of trial wavefunctions for importance sampling, the method is an excellent tool for…
Zero- and two-dimensional crystal defects form in open statistical ensembles, such as the grand canonical, that are usually inaccessible with conventional simulation techniques. This longstanding challenge is overcome with a new Hamiltonian…
The merger of two neutron stars has been proposed as a source of gamma-ray bursts, r-process elements, and detectable gravitational waves. Extracting information from observations of these phenomena requires fully relativistic simulations.…
A new algorithm for implementing the adaptive Monte Carlo method is given. It is used to solve the relativistic Boltzmann equations that describe the time evolution of a nonequilibrium electron-positron pair plasma containing high-energy…