Related papers: Pychastic: Precise Brownian Dynamics using Taylor-…
We develop efficient numerical methods for performing many-body Brownian dynamics simulations of a recently-observed fingering instability in an active suspension of colloidal rollers sedimented above a wall [M. Driscoll, B. Delmotte, M.…
Including the effect of thermal fluctuations in traditional computational fluid dynamics requires developing numerical techniques for solving the stochastic partial differential equations of fluctuating hydrodynamics. These Langevin…
We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust,…
We present a Python package for ground-state preparation based on the probabilistic imaginary-time evolution algorithm, with particular focus on its state-vector-based implementation. A standard shot-based simulation is also supported, and…
We present a differentiable soft-body physics simulator that can be composed with neural networks as a differentiable layer. In contrast to other differentiable physics approaches that use explicit forward models to define state…
Dynamics simulation with frictional contacts is important for a wide range of applications, from cloth simulation to object manipulation. Recent methods using smoothed lagged friction forces have enabled robust and differentiable simulation…
The solutions of parabolic and hyperbolic stochastic partial differential equations (SPDEs) driven by an infinite dimensional Brownian motion, which is a martingale, are in general not semi-martingales any more and therefore do not satisfy…
In this article we show in details the derivation of an integration scheme for the dissipative particle dynamic model (DPD) using the stochastic Trotter formula [De Fabritiis et al., Physica A, 361, 429 (2006)]. We explain some subtleties…
Recent developments in path integral methodology have significantly reduced the computational expense of including quantum mechanical effects in the nuclear motion in ab initio molecular dynamics simulations. However, the implementation of…
Numerical stability is a crucial requirement of reliable scientific computing. However, despite the pervasiveness of Python in data science, analyzing large Python programs remains challenging due to the lack of scalable numerical analysis…
We propose a numerical method for approximate calculations of the time evolution of point particle systems given only the system's Hamiltonian function and initial conditions. The method both generates and solves the equations of motion…
A Python package for post-processing of plane two-dimensional data from computational fluid dynamics simulations is presented. The package, called turbulucid, provides means for scripted, reproducible analysis of large simulation campaigns…
Probabilistic ordinary differential equation (ODE) solvers have been introduced over the past decade as uncertainty-aware numerical integrators. They typically proceed by assuming a functional prior to the ODE solution, which is then…
We study a class of semi-implicit Taylor-type numerical methods that are easy to implement and designed to solve multidimensional stochastic differential equations driven by a general rough noise, e.g. a fractional Brownian motion. In the…
We present a family of Python modules for the numerical integration of ordinary, delay, or stochastic differential equations. The key features are that the user enters the derivative symbolically and it is just-in-time-compiled, allowing…
PHYSBO (optimization tools for PHYSics based on Bayesian Optimization) is a Python library for fast and scalable Bayesian optimization. It has been developed mainly for application in the basic sciences such as physics and materials…
This paper studies the existence and uniqueness of solution of It\^o type stochastic differential equation $dx(t)=b(t, x(t), \om)dt+\si(t,x(t), \om) d B(t)$, where $B(t)$ is a fractional Brownian motion of Hurst parameter $H>1/2$ and…
We derive an exact Markovian kinetic equation for an oscillator linearly coupled to a heat bath, describing quantum Brownian motion. Our work is based on the subdynamics formulation developed by Prigogine and collaborators. The space of…
We present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point…
Simulation of contact and friction dynamics is an important basis for control- and learning-based algorithms. However, the numerical difficulties of contact interactions pose a challenge for robust and efficient simulators. A…