https://doi.org/10.1140/epjb/s10051-021-00165-0
Topical Review - Computational Methods
Mass-Zero constrained dynamics and statistics for the shell model in magnetic field
1
Centre Européen de Calcul Atomique et Moléculaire (CECAM), Ecole Polytechnique Fédérale de Lausanne, 1015, Lausanne, Switzerland
2
Department of Mathematical Sciences, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Torino, Italy
3
Institute for Applied Computing “Mauro Picone” (IAC), CNR Via dei Taurini 19, 00185, Rome, Italy
4
Università di Roma La Sapienza, Ple. A. Moro 5, 00185, Roma, Italy
5
School of Physics, University College of Dublin UCD-Belfield, Dublin 4, Ireland
b
alessandro.coretti@polito.it
Received:
22
April
2021
Accepted:
13
July
2021
Published online:
3
August
2021
In several domains of physics, including first principle simulations and classical models for polarizable systems, the minimization of an energy function with respect to a set of auxiliary variables must be performed to define the dynamics of physical degrees of freedom. In this paper, we discuss a recent algorithm proposed to efficiently and rigorously simulate this type of systems: the Mass-Zero (MaZe) Constrained Dynamics. In MaZe, the minimum condition is imposed as a constraint on the auxiliary variables treated as degrees of freedom of zero inertia driven by the physical system. The method is formulated in the Lagrangian framework, enabling the properties of the approach to emerge naturally from a fully consistent dynamical and statistical viewpoint. We begin by presenting MaZe for typical minimization problems where the imposed constraints are holonomic and summarizing its key formal properties, notably the exact Born–Oppenheimer dynamics followed by the physical variables and the exact sampling of the corresponding physical probability density. We then generalize the approach to the case of conditions on the auxiliary variables that linearly involve their velocities. Such conditions occur, for example, when describing systems in external magnetic field and they require to adapt MaZe to integrate semiholonomic constraints. The new development is presented in the second part of this paper and illustrated via a proof-of-principle calculation of the charge transport properties of a simple classical polarizable model of NaCl.
© The Author(s) 2021
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
