Related papers: Reaction rate calculation by parallel path swappin…
We present a novel approach to the parallelization of the parabolic fast multipole method for a space-time boundary element method for the heat equation. We exploit the special temporal structure of the involved operators to provide an…
The accuracy of rate constants calculated using transition state theory depends crucially on the correct identification of a recrossing--free dividing surface. We show here that it is possible to define such optimal dividing surface in…
In this paper we demonstrate that tempering Markov chain Monte Carlo samplers for Bayesian models by recursively subsampling observations without replacement can improve the performance of baseline samplers in terms of effective sample size…
The impact of various types of external potentials to the Peyrard-Bishop DNA denaturation is investigated through statistical mechanics approach. The partition function is obtained using transfer integral method, and further the stretching…
While gradient-based discrete samplers are effective in sampling from complex distributions, they are susceptible to getting trapped in local minima, particularly in high-dimensional, multimodal discrete distributions, owing to the…
We present an extension to the time-dependent photo-ionization code C$^2$-Ray to calculate photo-heating in an efficient and accurate way. In C$^2$-Ray, the thermal calculation demands relatively small time-steps for accurate results. We…
The Random Phase Approximation (RPA) is a widely employed post Hartree-Fock or DFT method, capable of capturing van der Waal interactions and other dynamic correlation effects at relatively low costs of $\mathcal O(N^3)$ in time and…
Frequently, the design of physicochemical processes requires screening of large numbers of alternative designs with complex geometries. These geometries may result in conformal meshes which introduce stability issues, significant…
Ab initio instanton rate theory is a computational method for rigorously including tunnelling effects into calculations of chemical reaction rates based on a potential-energy surface computed on the fly from electronic-structure theory.…
Spontaneous structural rearrangements play a central role in the organization and function of complex biomolecular systems. In principle, physics-based computer simulations like Molecular Dynamics (MD) enable us to investigate these…
RNA polymerase (RNAP) elongates RNA by walking along a DNA template and selectively incorporating ribonucleoside triphosphates (rNTPs). Rather than mechanically replicating the base sequence, RNAP conditions binding and chemistry on the…
Bistable biochemical switches are ubiquitous in gene regulatory networks and signal transduction pathways. Their switching dynamics, however, are difficult to study directly in experiments or conventional computer simulations, because…
Simulating rare events, such as the transformation of a reactant into a product in a chemical reaction typically requires enhanced sampling techniques that rely on heuristically chosen collective variables (CVs). We propose using…
Characterizing thermally activated transitions in high-dimensional rugged energy surfaces is a very challenging task for classical computers. Here, we develop a quantum annealing scheme to solve this problem. First, the task of finding the…
We extend replica exchange simulation in two ways, and apply our approaches to biomolecules. The first generalization permits exchange simulation between models of differing resolution -- i.e., between detailed and coarse-grained models.…
Several different methods exist for efficient approximation of paths in multiscale stochastic chemical systems. Another approach is to use bursts of stochastic simulation to estimate the parameters of a stochastic differential equation…
We recently introduced a novel replica-exchange scheme in which an individual replica can sample from states encountered by other replicas at any previous time by way of a global configuration database, enabling the fast propagation of…
Molecular Dynamics (MD) simulations are fundamental computational tools for the study of proteins and their free energy landscapes. However, sampling protein conformational changes through MD simulations is challenging due to the relatively…
We show that the acceptance probability for swaps in the parallel tempering Monte Carlo method for classical canonical systems is given by a universal function that depends on the average statistical fluctuations of the potential and on the…
We propose Parallelised Diffeomorphic Sampling-based Motion Planning (PDMP). PDMP is a novel parallelised framework that uses bijective and differentiable mappings, or diffeomorphisms, to transform sampling distributions of sampling-based…