Related papers: Exploration of Order in Chaos with Replica Exchang…
We present an algorithm to sample stochastic differential equations conditioned on rather general constraints, including integral constraints, endpoint constraints, and stochastic integral constraints. The algorithm is a pathspace…
We study classically the problem of two relativistic particles with an invariant Duffing-like potential which reduces to the usual Duffing form in the nonrelativistic limit. We use a special relativistic generalization (RGEM) of the…
Sequential Monte Carlo methods, also known as particle methods, are a popular set of techniques for approximating high-dimensional probability distributions and their normalizing constants. These methods have found numerous applications in…
The inverse problem for a disordered system involves determining the interparticle interaction parameters consistent with a given set of experimental data. Recently, Rutledge has shown (Phys. Rev. E63, 021111 (2001)) that such problems can…
It is shown that a coupled map model for open flow may exhibit spatial chaos and spatial quasiperiodicity with temporal periodicity. The locations of these patterns, which cover a substantial part of parameter space, are indicated in a…
Rare events are events that are expected to occur infrequently, or more technically, those that have low probabilities (say, order of $10^{-3}$ or less) of occurring according to a probability model. In the context of uncertainty…
Partial differential equation is a powerful tool to characterize various physics systems. In practice, measurement errors are often present and probability models are employed to account for such uncertainties. In this paper, we present a…
We present novel roulette schemes for rare-event sampling that are both structure-preserving and unbiased. The boundaries where Monte Carlo markers are split and deleted are placed automatically and adapted during runtime. Extending…
The bifurcation method is a way to do rare event sampling -- to estimate the probability of events that are too rare to be found by direct simulation. We describe the bifurcation method and use it to estimate the transition rate of a double…
The behaviour of the one--dimensional random--forced Burgers equation is investigated in the path integral formalism, using a discrete space--time lattice. We show that by means of Monte Carlo methods one may evaluate observables, such as…
A method to reduce or enhance chaos in Hamiltonian flows with two degrees of freedom is discussed. This method is based on finding a suitable perturbation of the system such that the stability of a set of periodic orbits changes (local…
This paper proposes strategies for designing a system whose computational model is subject to aleatory and epistemic uncertainty. Aleatory variables, which are caused by randomness in physical parameters, are draws from a possibly unknown…
Macroscopically heterogeneous materials, characterised mostly by comparable heterogeneity lengthscale and structural sizes, can no longer be modelled by deterministic approach instead. It is convenient to introduce stochastic approach with…
We propose a novel targeted exploration strategy designed specifically for uncertain linear time-invariant systems with energy-bounded disturbances, i.e., without any assumptions on the distribution of the disturbances. We use classical…
We apply the diagrammatic Monte Carlo approach to three-dimensional Fermi-polaron systems with mass-imbalance, where an impurity interacts resonantly with a noninteracting Fermi sea whose atoms have a different mass. This method allows to…
We propose a flat-histogram Monte Carlo method to efficiently sample fractal landscapes such as escape time functions of open chaotic systems. This is achieved by using a random-walk step which depends on the height of the landscape via the…
We propose a new implementation of the replica-exchange method (REM) in which replicas follow a pre-planned route in temperature space instead of a random walk. Our method satisfies the detailed balance condition in the proposed route. The…
The parallel computational complexity of the quadratic map is studied. A parallel algorithm is described that generates typical pseudotrajectories of length t in a time that scales as log t and increases slowly in the accuracy demanded of…
Monte Carlo is famous for accepting model extensions and model refinements up to infinite dimension. However, this powerful incremental design is based on a premise which has severely limited its application so far: a state-variable can…
In Hamiltonian systems subjected to periodic perturbations the stable and unstable manifolds of the unstable periodic orbits provide the dynamical "skeleton" that drives the mixing process and bounds the chaotic regions of the phase space.…