Abstract:The proper choice of collective variables (CVs) is central to biased-sampling free energy reconstruction methods in molecular dynamics simulations. The PLUMED 2 library, for instance, provides several sophisticated CV choices, implemented in a C++ framework; however, developing new CVs is still time consuming due to the need to provide code for the analytical derivatives of all functions with respect to atomic coordinates. We present two solutions to this problem, namely (a) symbolic differentiation and code generation, and (b) automatic code differentiation, in both cases leveraging open-source libraries (SymPy and Stan Math, respectively). The two approaches are demonstrated and discussed in detail implementing a realistic example CV, the local radius of curvature of a polymer. Users may use the code as a template to streamline the implementation of their own CVs using high-level constructs and automatic gradient computation. Program Program Title: Practical approaches to the differentiation of collective variables in free energy codes: computer-algebra code generation and automatic differentiation Program Files doi: http://dx.doi.org/10.17632/r4r67bvkdn.1 Licensing provisions: GNU Lesser General Public License Version 3 (LGPL-3) Programming languages: C++, Python Nature of problem: The C++ implementation of collective variables (CVs, functions of atomic coordinates to be used in biased sampling applications) in biasing libraries for atomistic simulations, such as PLUMED [1], requires computation of both the variable and its gradient with respect to the atomic coordinates; coding and testing the analytical derivatives complicate the implementation of new CVs. Solution method: The paper shows two approaches to automate the computation of CV gradients, namely, symbolic differentiation with code generation and automatic code differentiation, demonstrating their implementation entirely with open-source software (respectively, SymPy and the Stan Math Library). Additional comments: The paper’s accompanying code serves as an example and template for the methods described in the paper; it is distributed as the two modules curvature_codegen and curvature_autodiff integrated in PLUMED 2’s source tree; the latest version is available at https://github.com/tonigi/plumed2-automatic-gradients . [1] Tribello GA, Bonomi M, Branduardi D, Camilloni C, Bussi G. PLUMED 2: New feathers for an old bird. Computer Physics Communications. 2014 Feb;185(2):604-13.

Code and Data Repository for Convergence Rates of Zeroth Order Gradient Descent for Łojasiewicz Functions

Convergence Rates of Zeroth Order Gradient Descent for Łojasiewicz Functions

Convergence Rates of Stochastic Zeroth-order Gradient Descent for \L ojasiewicz Functions

The Exponential Convergence of the CR Yamabe Flow

Convergence Analysis of Gradient Algorithms on Riemannian Manifolds Without Curvature Constraints and Application to Riemannian Mass

Laplacian Smoothing Gradient Descent

Lévy Gradient Descent: Augmented Random Search for Geophysical Inverse Problems

Convergence Analysis of a Finite Difference Scheme for the Gradient Flow associated with the ROF Model

Linear Convergence of Inexact Descent Method and Inexact Proximal Gradient Algorithms for Lower-Order Regularization Problems

Linear Convergence of Subgradient Algorithm for Convex Feasibility on Riemannian Manifolds

Code and Data Repository for A FAST Method for Nested Estimation

Inexact Riemannian Gradient Descent Method for Nonconvex Optimization

Convergence error analysis of reflected gradient Langevin dynamics for non-convex constrained optimization

On the Convergence of AdaGrad(Norm) on $\R^{d}$: Beyond Convexity, Non-Asymptotic Rate and Acceleration

How to differentiate collective variables in free energy codes: Computer-algebra code generation and automatic differentiation

Convergence Error Analysis of Reflected Gradient Langevin Dynamics for Globally Optimizing Non-Convex Constrained Problems

HPC-driven computational reproducibility in numerical relativity codes: A use case study with IllinoisGRMHD

Statistical Inference for Polyak-Ruppert Averaged Zeroth-order Stochastic Gradient Algorithm

Preconditioned Descent Algorithms for P-Laplacian.

The Average Rate of Convergence of the Exact Line Search Gradient Descent Method

Convergence Rate Bounds for the Mirror Descent Method: IQCs, Popov Criterion and Bregman Divergence