Related papers: Hessian-free force-gradient integrators and their …
Differentiable programming is revolutionizing computational science by enabling automatic differentiation (AD) of numerical simulations. While first-order gradients are well-established, second-order derivatives (Hessians) for implicit…
A consequent approach is proposed to construct symplectic force-gradient algorithms of arbitrarily high orders in the time step for precise integration of motion in classical and quantum mechanics simulations. Within this approach the basic…
Long-range interactions in finite density QCD necessitate a non-perturbative approach in order to reliably map out the key features and spectrum of the QCD phase diagram. However, the complex nature of the fermion determinant in this sector…
Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian…
In most of mesh-free methods, the calculation of interactions between sample points or particles is the most time consuming. When we use mesh-free methods with high spatial orders, the order of the time integration should also be high. If…
Heavy-light mesons, heavy quarkonium and doubly heavy baryons are briefly discussed. Effective field theories (EFTs) of QCD based on the heavy quark mass expansion 1/m_Q provide a unified framework to describe all three systems. They…
Simulating noninteracting fermion systems is a common task in computational many-body physics. In absence of translational symmetries, modeling free fermions on $N$ modes usually requires poly$(N)$ computational resources. While often…
High dimensional integrals are abundant in many fields of research including quantum physics. The aim of this paper is to develop efficient recursive strategies to tackle a class of high dimensional integrals having a special product…
Quantum simulations of lattice gauge theories are anticipated to directly probe the real time dynamics of QCD, but scale unfavorably with the required truncation of the gauge fields. Improved Hamiltonians are derived to correct for the…
In this work we develop a stochastic algorithm to integrate the Cahn-Hilliard equations. The algorithm is based on Gillespie's stochastic simulation algorithm, also known as kinetic Monte Carlo. The deterministic integration of the phase…
This work introduces the nested-set Hessian approximation, a second-order approximation method that can be used in any derivative-free optimization routine that requires such information. It is built on the foundation of the generalized…
Since they were introduced in the 1990s, Lie group integrators have become a method of choice in many application areas. These include multibody dynamics, shape analysis, data science, image registration and biophysical simulations. Two…
The supercomputing platforms available for high performance computing based research evolve at a great rate. However, this rapid development of novel technologies requires constant adaptations and optimizations of the existing codes for…
The effective electroweak Hamiltonian in the gradient-flow formalism is constructed for the current-current operators through next-to-next-to-leading order QCD. The results are presented for two common choices of the operator basis. This…
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…
Force-gradient decomposition methods are used to improve the energy preservation of symplectic schemes applied to Hamiltonian systems. If the potential is composed of different parts with strongly varying dynamics, this multirate potential…
Accelerating the convergence of second-order optimization, particularly Newton-type methods, remains a pivotal challenge in algorithmic research. In this paper, we extend previous work on the \textbf{Quadratic Gradient (QG)} and rigorously…
This paper proposes a hybrid learning and optimization framework for mobile manipulators for complex and physically interactive tasks. The framework exploits an admittance-type physical interface to obtain intuitive and simplified human…
Recent results have shown the stability of frustration-free Hamiltonians to weak local perturbations, assuming several conditions. In this paper, we prove the stability of free fermion Hamiltonians which are gapped and local. These free…
We rigorously establish the asymptotic equivalence between the height function of interacting dimers on the square lattice and the massless Gaussian free field. Our theorem explains the microscopic origin of the sine-Gordon field theory…