Related papers: Pychastic: Precise Brownian Dynamics using Taylor-…
We present an algorithm for the simulation of the exact real-time dynamics of classical many-body systems with discrete energy levels. In the same spirit of kinetic Monte Carlo methods, a stochastic solution of the master equation is found,…
Brownian dynamics algorithms integrate numerically Langevin equations and allow to probe long time scales in simulations. A common requirement for such algorithms is that interactions in the system should vary little during an integration…
Hamiltonian systems are differential equations which describe systems in classical mechanics, plasma physics, and sampling problems. They exhibit many structural properties, such as a lack of attractors and the presence of conservation…
Brownian motion in one or more dimensions is extensively used as a stochastic process to model natural and engineering signals, as well as financial data. Most works dealing with multidimensional Brownian motion consider the different…
We present a major update to QuSpin, SciPostPhys.2.1.003 -- an open-source Python package for exact diagonalization and quantum dynamics of arbitrary boson, fermion and spin many-body systems, supporting the use of various (user-defined)…
The connection between fundamental interactions acting in molecules in a fluid and macroscopically measured properties, such as the viscosity between colloidal particles coated with polymers, is studied here. The role that hydrodynamic and…
Multiple software packages currently exist for the computation of bulk topological invariants in both idealized tight-binding models and realistic Wannier tight-binding models derived from density functional theory. Currently, only one…
We introduce a framework for approximate dynamic programming that we apply to discrete time chains on $\mathbb{Z}_+^d$ with countable action sets. Our approach is grounded in the approximation of the (controlled) chain's generator by that…
We present the Python Tree Tensor Network package (pyTTN) for the evaluation of dynamical properties of closed and open quantum systems that makes use of Tree Tensor Network (TTN), or equivalently the multi-layer multiconfiguration…
We investigate the feasibility of modelling turbulence via numeric functional integration. By transforming the Burgers' equation into a functional integral we are able to calculate equal-time spatial correlation of system variables using…
In these lecture notes, we explore the mathematical preliminaries and foundational concepts that connect stochastic processes with partial differential equations. We begin by investigating Brownian motion, which serves as a model for random…
The position-based dynamics (PBD) algorithm is a popular and versatile technique for real-time simulation of deformable bodies, but is only applicable to forces that can be expressed as linearly compliant constraints. In this work, we…
A novel approach to account for hard-body interactions in (overdamped) Brownian dynamics simulations is proposed for systems with non-vanishing force fields. The scheme exploits the analytically known transition probability for a Brownian…
We develop a new algorithm for the Brownian dynamics of soft matter systems that evolves time by spatially correlated Monte Carlo moves. The algorithm uses vector wavelets as its basic moves and produces hydrodynamics in the low Reynolds…
It is shown that under a certain condition on a semimartingale and a time-change, any stochastic integral driven by the time-changed semimartingale is a time-changed stochastic integral driven by the original semimartingale. As a direct…
A software package, called DFTBaby, is published, which provides the electronic structure needed for running non-adiabatic molecular dynamics simulations at the level of charge-consistent tight-binding DFT. A long-range correction is…
In the Ermak-McCammon algorithm for Brownian Dynamics, the hydrodynamic interactions (HI) between N spherical particles are described by a 3N x 3N diffusion tensor. This tensor has to be factorized at each timestep with a runtime of O(N^3),…
Brownian dynamics simulations of fairly long, highly detailed polymer chains, at the resolution of a single Kuhn step, remains computationally prohibitive even on the modern processors. This is especially true when the beads on the chain…
The solution of a parabolic stochastic partial differential equation (SPDE) driven by an infinite-dimensional Brownian motion is in general not a semi-martingale anymore and does in general not satisfy an It\^{o} formula like the solution…
Automated data-driven modeling, the process of directly discovering the governing equations of a system from data, is increasingly being used across the scientific community. PySINDy is a Python package that provides tools for applying the…