Related papers: Hessian-free force-gradient integrators and their …
In this paper, we derive a variational integrator for certain highly oscillatory problems in mechanics. To do this, we take a new approach to the splitting of fast and slow potential forces: rather than splitting these forces at the level…
We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian…
Master-field simulations offer an approach to lattice QCD in which calculations are performed on a small number of large-volume gauge-field configurations. The latter is advantageous for simulations in which the global topological charge is…
Progress in computing the spectrum of excited baryons and mesons in lattice QCD is described. Large sets of spatially-extended hadron operators are used. The need for multi-hadron operators in addition to single-hadron operators is…
We propose an efficient approach for time integration of Klein-Gordon equations with highly oscillatory in time input terms. The new methods are highly accurate in the entire range, from slowly varying up to highly oscillatory regimes. Our…
This article presents a Hamiltonian lattice formulation of static Casimir systems at a level of generality appropriate for an introductory investigation. Background structure - represented by a lattice potential V(x) - is introduced along…
A practical and simple stable method for calculating Fourier integrals is proposed, effective both at low and at high frequencies. An approach based on the fruitful idea of Levin, to use of the collocation method to approximate the slowly…
We introduce a novel technique for constructing higher-order variational integrators for Hamiltonian systems of ODEs. In particular, we are concerned with generating globally smooth approximations to solutions of a Hamiltonian system. Our…
A novel finite element formulation for gradient-regularized damage models is presented which allows for the robust, efficient, and mesh-independent simulation of damage phenomena in engineering and biological materials. The paper presents a…
The rise of exascale supercomputers has fueled competition among GPU vendors, driving lattice QCD developers to write code that supports multiple APIs. Moreover, new developments in algorithms and physics research require frequent updates…
Stable partitioned techniques for simulating unsteady fluid-structure interaction (FSI) are known to be computationally expensive when high added-mass is involved. Multiple coupling strategies have been developed to accelerate these…
One of the major frontiers of lattice field theory is the inclusion of light fermions in simulations, particularly in pursuit of accurate, first principles predictions from lattice QCD. With dedicated Teraflops-scale computers currently…
Lattice field theory is a non-perturbative tool for studying properties of strongly interacting field theories, which is particularly amenable to numerical calculations and has quantifiable systematic errors. In these lectures we apply…
Second-order methods for neural network optimization have several advantages over methods based on first-order gradient descent, including better scaling to large mini-batch sizes and fewer updates needed for convergence. But they are…
Lattice models are valuable tools to gain insight into the statistical physics of heteropolymers. We rigorously map the partition function of these models into a vacuum expectation value of a $\mathbb{Z}_2$ lattice gauge theory (LGT), with…
The Hybrid Integral Gain Systems (HIGS) has recently gained a lot of attention in the control of precision motion systems. HIGS is a nonlinear low pass filter with a 52 degree phase advantage over its linear counterpart. This property…
Differentiable simulators promise faster computation time for reinforcement learning by replacing zeroth-order gradient estimates of a stochastic objective with an estimate based on first-order gradients. However, it is yet unclear what…
Four-Fermi quantum field theories in (2+1) dimensions lie among the simplest models in high-energy physics, the understanding of which requires a non-perturbative lattice formulation addressing their strongly-coupled fixed points. These…
Many applications in computational physics that use numerical integrators based on splitting and composition can benefit from the development of optimized algorithms and from choosing the best ordering of terms. The cost in programming and…
We propose new semi-implicit numerical methods for the integration of the stochastic Landau-Lifshitz equation with built-in angular momentum conservation. The performance of the proposed integrators is tested on the 1D Heisenberg chain. For…