Related papers: A New Time-Reversible Integrator for Molecular Dyn…
Molecular-scale computation is crucial for smart materials and nanoscale devices, yet creating single-molecule systems capable of complex computations remains challenging. We present a theoretical framework for a single-molecule computer…
This is the second and final article on the tutorial on manipulator differential kinematics. In Part 1, we described a method of modelling kinematics using the elementary transform sequence (ETS), before formulating forward kinematics and…
We propose a time-adaptive predictor/multi-corrector method to solve hyperbolic partial differential equations, based on the generalized-$\alpha$ scheme that provides user-control on the numerical dissipation and second-order accuracy in…
A transient magneto-quasistatic vector potential formulation involving nonlinear material is spatially discretized using the finite element method of first and second polynomial order. By applying a generalized Schur complement the…
The depth of networks plays a crucial role in the effectiveness of deep learning. However, the memory requirement for backpropagation scales linearly with the number of layers, which leads to memory bottlenecks during training. Moreover,…
We consider a new splitting based on the Sherman-Morrison-Woodbury formula, which is particularly effective with iterative methods for the numerical solution of large linear systems. These systems involve matrices that are perturbations of…
Newton-type solvers have been extensively employed for solving a variety of nonlinear system of algebraic equations. However, for some complex nonlinear system of algebraic equations, efficiently solving these systems remains a challenging…
We consider the problem of motion of several rigid bodies immersed in a perfect compressible fluid. Using the method of convex integration we establish the existence of infinitely many weak solutions with {\it a priori} prescribed motion of…
This paper invites the reader to experiment with an easy-to-use MATLAB implementation of Metropolis integrators for Molecular Dynamics (MD) simulation. These integrators are analysis-based, in the sense that they can rigorously simulate…
A novel approach is introduced for deriving exact solutions to nonlinear systems of ordinary differential equations. This method consists of four parts. In the initial part, the examined nonlinear differential equation system is transformed…
The time irreversibility problem is the dichotomy of the reversible microscopic dynamics and the irreversible macroscopic physics. This problem was considered by Boltzmann, Poincar\'e, Bogolyubov and many other authors and though some…
We derive a new criterion for estimating characteristic dynamical timescales in N-body simulations. The criterion uses the second, third, and fourth derivatives of particle positions: acceleration, jerk, and snap. It can be used for…
Molecular dynamics simulations are a cornerstone in science, allowing to investigate from the system's thermodynamics to analyse intricate molecular interactions. In general, to create extended molecular trajectories can be a…
A new algorithm is introduced to integrate the equations of rotational motion. The algorithm is derived within a leapfrog framework and the quantities involved into the integration are mid-step angular momenta and on-step orientational…
We derived a number of numerical methods to treat biomolecular systems with multiple time scales. Based on the splitting of the operators associated with the slow-varying and fast-varying forces, new multiple time-stepping (MTS) methods are…
Explicit step-truncation tensor methods have recently proven successful in integrating initial value problems for high-dimensional partial differential equations (PDEs). However, the combination of non-linearity and stiffness may introduce…
Quantum mechanics based ab-initio molecular dynamics (MD) simulation schemes offer an accurate and direct means to monitor the time-evolution of materials. Nevertheless, the expensive and repetitive energy and force computations required in…
A very simple and accurate numerical method which is applicable to systems of differentio-integral equations with quite general boundary conditions has been devised. Although the basic idea of this method stems from the Keller Box method,…
In this paper, we extend several time reversible numerical integrators to solve the Lorentz force equations from second order accuracy to higher order accuracy for relativistic charged particle tracking in electromagnetic fields. A fourth…
The choice of numerical integrator in approximating solutions to dynamic partial differential equations depends on the smallest time-scale of the problem at hand. Large-scale deformations in elastic solids contain both shear waves and bulk…