Jul 2026· Journal of Chemical Physics· Vol 165 2· 0 citations· 14 references
Medicine
Abstract
Rare-event molecular dynamics simulations are limited by the separation between atomic and activated timescales. Hyperdynamics accelerates these processes via bias potentials while preserving correct transition kinetics. The minimum-mode following (MMF) formulation provides a ridge-based bias construction but typically employs an identity approximation for the Jacobian together with local force evaluation, which introduces inconsistency in high-dimensional systems. Here, we develop a Jacobian-corrected MMF (J-MMF) method that improves both accuracy and consistency. The method introduces a path-informed semi-analytical Jacobian based on the sequential evolution of minimum-mode vectors and an orthogonal projection that removes spurious force components along the reaction coordinate, without additional force evaluations. Benchmark simulations of adatom diffusion on the Cu(100) surface (1-201 mobile atoms) demonstrate improved accuracy and robust, largely size-independent acceleration compared with standard MMF and bond-boost methods. J-MMF provides a scalable and consistent framework for rare-event simulations in complex atomistic systems.
Direct inversion in the iterative subspace (DIIS), introduced in 1980, remains the practical default for accelerating Hartree-Fock and Kohn-Sham self-consistent-field (SCF) calculations. Many later methods improve robustness or iteration count, but the near-zero overhead of DIIS makes a broad reduction in total wall ti...
The framework connects force tolerances to reference-energy transfer, providing a physical basis for potential assessment and adaptive reference allocation, and establishes a directional residual-work coefficient combining directional curvature mismatch with the spatial distribution of the residual response.
In modern density functional theory, effective Hamiltonians are constructed to reproduce densities and forces of electrons at equilibrium. However, their nonlocal potentials leave systematic errors in spectral properties and real-time dynamics. Access to local optimized or inverted effective potentials would remove thi...
Matija Medvidović, Angel Rubio, Juan Carrasquilla· 1 citation
We present a matrix-free residual-minimization method for preparing guiding-center distribution functions on finite grids. The method minimizes a fixed conservative discrete residual in the physical Jacobian metric using Polak-Ribiere+ iterations with exact quadratic line searches. JAX/XLA evaluates the operator and it...
The coupled-cluster (CC) equations are most frequently solved via fixed-point (FP) iterations. However, when formulated in a non-canonical gauge, as in local correlation CC, the FP iteration may converge slowly or even diverge. Practical fixes, such as level-shifting and a direct inversion of iterative subspace (DIIS),...
Chanaka D. M. Mudiyanselage, Kangbo Li, F. Faulstich· 0 citations
The kinetic Monte Carlo (KMC) method is becoming increasingly popular in the research of heterogeneous catalysts, yielding fundamental understanding of, e.g., their activity, selectivity, and rate of poisoning. KMC simulations can provide a link between ab initio electronic structure calculations and experiments by sim...
Alec Elías Sigurðarson, M. Stamatakis· Journal of Chemical Physics· 0 citations
We use cookies to run the site and, with your consent, for analytics and to show ads.
See our Cookie Policy.