Related papers: Metric-tensor flexible-cell algorithm for isotherm…
In this paper, we present a novel general framework grounded in the factor graph theory to solve kinematic and dynamic problems for multi-body systems. Although the motion of multi-body systems is considered to be a well-studied problem and…
We construct a general metric-tensor framework for treating inhomogenous adiabatic deformations applied to crystalline insulators, by deriving an effective time-dependent Schr\"odinger equation in the undistorted frame. The response can be…
A recent reformulation [1] of the problem of Coulomb gases in the presence of a dynamical dielectric medium showed that finite temperature simulations of such systems can be accomplished on the basis of completely local Hamiltonians on a…
Geometric integrators of the Schr\"{o}dinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement, but, unfortunately, is restricted to…
We present a computational framework for two-scale asymptotic homogenization to determine the intrinsic magnetic permeability of composites. To this end, considering linear magnetostatics, both vector and scalar potential formulations are…
Most numerical integration algorithms are not designed specifically for Hamiltonian systems and do not respect their characteristic properties, which include the preservation of phase space volume with time. This can lead to spurious…
A common approach for generating an anisotropic mesh is the M-uniform mesh approach where an adaptive mesh is generated as a uniform one in the metric specified by a given tensor M. A key component is the determination of an appropriate…
The presented paper concentrates on the boundary element method (BEM) for the heat equation in three spatial dimensions. In particular, we deal with tensor product space-time meshes allowing for quadrature schemes analytic in time and…
We present a new approach to simulating Hamiltonian dynamics based on implementing linear combinations of unitary operations rather than products of unitary operations. The resulting algorithm has superior performance to existing simulation…
The popular method of Nose and Hoover to create canonically distributed positions and momenta in classical molecular dynamics simulations is generalized to a genuine quantum system of infinite dimensionality. We show that for the quantum…
We develop a density matrix formalism to describe coupled electron-nuclear dynamics. To this end we introduce an effective Hamiltonian formalism that describes electronic transitions and small (quantum) nuclear fluctuations along a…
We study asymptotically constrained systems for numerical integration of the Einstein equations, which are intended to be robust against perturbative errors for the free evolution of the initial data. First, we examine the previously…
We present a simple choice of integration variables that can be used to exploit the near-integrable character of problems in celestial mechanics. The approach is based on the well-known principle of variation of parameters: instead of…
Conventional molecular dynamics simulations macromolecules require long computational times because the most interesting motions are very slow compared with the fast oscillations of bond lengths and bond angles that limit the integration…
A new anisotropic mesh adaptation strategy for finite element solution of elliptic differential equations is presented. It generates anisotropic adaptive meshes as quasi-uniform ones in some metric space, with the metric tensor being…
Flows on symplectic, Poisson, contact, and metriplectic manifolds are reviewed in order to describe our main result, which is to associate a natural metriplectic dynamical system on the general one-jet bundle $J^1N=T^*N\times \mathbb{R}$,…
Polynomial approximations to the inverse of the fermion matrix are used to filter the dynamics of the upper energy scales in HMC simulations. The use of a multiple time-scale integration scheme allows the filtered pseudofermions to be…
A new method is proposed to numerically integrate a dynamical system on a manifold such that the trajectory stably remains on the manifold and preserves first integrals of the system. The idea is that given an initial point in the manifold…
This paper presents a Lagrangian approach to simulating multibody dynamics in a tensegrity framework with an ability to tackle holonomic constraint violations in an energy-preserving scheme. Governing equations are described using…
It has previously been shown that varying the numerical timestep during a symplectic orbital integration leads to a random walk in energy and angular momentum, destroying the phase space-conserving property of symplectic integrators. Here…