Related papers: Cayley modification for strongly stable path-integ…
Molecular dynamics (MD) simulations provide atomistic insights into the structure, dynamics, and function of biomolecules by generating time-resolved, high-dimensional trajectories. Analyzing such data benefits from estimating the minimal…
The theoretical justification of the Hybrid Monte Carlo algorithm depends upon the molecular dynamics trajectories within it being exactly reversible. If computations were carried out with exact arithmetic then it would be easy to ensure…
The paper deals with the problem of integration of equations of motion in nonholonomic systems. By means of well-known theory of the differential equations with an invariant measure the new integrable systems are discovered. Among them…
A Hamiltonian formulation of generic many-particle systems with space-dependent balanced loss and gain coefficients is presented. It is shown that the balancing of loss and gain necessarily occurs in a pair-wise fashion. Further, using a…
A general and rigorous methodology to compute the quantum equilibrium isotope effect is described. Unlike standard approaches, ours does not assume separability of rotational and vibrational motions and does not make the harmonic…
The integration time step is a critical determinant of performance in molecular dynamics simulations, governing the trade-off between speed and fidelity. Although 2 fs remains the standard in atomistic biomolecular simulations, the push for…
Internal coordinate molecular dynamics (ICMD) is an efficient method for studying biopolymers, but it is readily applicable only to molecules with tree topologies, that is with no internal flexible rings. Common examples violating this…
We consider the class of stationary-increment harmonizable stable processes with infinite control measure, which most notably includes real harmonizable fractional stable motions. We give conditions for the integrability of the paths of…
We review and compare numerical methods that simultaneously control temperature while preserving the momentum, a family of particle simulation methods commonly used for the modelling of complex fluids and polymers. The class of methods…
Due to particle conservation, Canonical Molecular Dynamics (MD) simulations fail in the description of surface phase transitions involving coverage or lateral density changes. However, a step on the surface can act effectively as a source…
The magnetohydrodynamic-kinetic (MHD-kinetic) hybrid model [Park et. al., 1992] has been widely applied in studying energetic particles (EPs) problems in fusion plasmas for past decades. The pressure-coupling scheme or the current-coupling…
This paper is devoted to the numerical symplectic approximation of the charged-particle dynamics (CPD) with arbitrary electromagnetic fields. By utilizing continuous-stage methods and exponential integrators, a general class of symplectic…
Predictive modeling of the phonon/thermal transport properties of materials is vital to rational design for a diverse spectrum of engineering applications. Classical Molecular Dynamics (MD) simulations serve as a tool to simulate the time…
Synthetic molecular dynamics (synMD) trajectories from learned generative models have been proposed as a useful addition to the biomolecular simulation toolbox. The computational expense of explicitly integrating the equations of motion in…
Molecular dynamics (MD) simulation is a widely used technique to simulate molecular systems, most commonly at the all-atom resolution where equations of motion are integrated with timesteps on the order of femtoseconds…
Chemical space which encompasses all stable compounds is unfathomably large and its dimension scales linearly with the number of atoms considered. The success of machine learning methods suggests that many physical quantities exhibit…
A proper treatment of electrostatic interactions is crucial for the accurate calculation of forces in computer simulations. Electrostatic interactions are typically modeled using Ewald based methods, which have become one of the…
In this paper we present a Path Integral Monte Carlo (PIMC) simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles. This work represents a quantum…
Attochemistry aims to exploit the properties of coherent electronic wavepackets excited via attosecond pulses, to control the formation of photoproducts. Such molecular processes can in principle be simulated with various nonadiabatic…
We employ a linear stability analysis and direct numerical simulations to study the characteristics of wall-modes in thermal convection in a rectangular box under strong and inclined magnetic fields. The walls of the convection cell are…