Search

The Open Databases Integration for Materials Design (OPTIMADE) application programming interface (API) empowers users with holistic access to a growing federation of databases, enhancing the accessibility and discoverability of materials and chemical data. Since the first release of the OPTIMADE specification (v1.0), the API has undergone significant development, leading to the upcoming v1.2 release, and has underpinned multiple scientific studies. In this work, we highlight the latest features of the API format, accompanying software tools, and provide an update on the implementation of OPTIMADE in contributing materials databases. We end by providing several use cases that demonstrate the utility of the OPTIMADE API in materials research that continue to drive its ongoing development.

Interactive notebooks are a precious tool for creating graphical user interfaces and teaching materials. Python and Jupyter are becoming increasingly popular in this context, with Jupyter widgets at the core of the interactive functionalities. However, while packages and libraries which offer a broad range of general-purpose widgets exist, there is limited development of specialized widgets for computational physics, chemistry and materials science. This deficiency implies significant time investments for the development of effective Jupyter notebooks for research and education in these domains. Here, we present custom Jupyter widgets that we have developed to target the needs of these communities. These widgets constitute high-quality interactive graphical components and can be employed, for example, to visualize and manipulate data, or to explore different visual representations of concepts, clarifying the relationships existing between them. In addition, we discuss with one example how similar functionality can be exposed in the form of JupyterLab extensions, modifying the JupyterLab interface for an enhanced user experience when working with applications within the targeted scientific domains., Comment: 19 pages, 8 figures

A new algorithm for efficient and fully time-reversible integration of first-principles molecular dynamics based on orbital-free density functional theory (OFDFT) is presented. The algorithm adapts to this nontrivial case the recently introduced Mass-Zero (MaZe) constrained dynamics. The formalism ensures that full adiabatic separation is enforced between nuclear and electronic degrees of freedom and, consequently, that the exact Born-Oppenheimer probability for the nuclei is sampled. Numerical integration of the MaZe dynamics combines standard molecular dynamics algorithms, e.g. Verlet or velocity Verlet, with the SHAKE method to impose the minimum conditions on the electronic degrees of freedom as a set of constraints. The developments presented in this work, that include a bespoke adaptation of the standard SHAKE algorithm, ensure that the quasilinear scaling of OFDFT is preserved by the new method for a broad range of kinetic and exchange-correlation functionals, including nonlocal ones. The efficiency and accuracy of the approach is demonstrated via calculations of static and dynamic properties of liquid sodium in the constant energy and constant temperature ensembles.

In this paper we present the Open Software Services for Classrooms and Research (OSSCAR) platform. OSSCAR provides an open collaborative environment to develop, deploy and access educational resources in the form of web applications. To minimize efforts in the creation and use of new educational material, it combines software tools that have emerged as standards with custom domain-specific ones. The technical solutions adopted to create and distribute content are described and motivated on the basis of reliability, sustainability, ease of uptake and use. Examples from courses in the domains of physics, chemistry, and materials science are shown to demonstrate the style and level of interactivity of typical applications. The tools presented are easy to use, and create a uniform and open environment exploitable by a large community of teachers, students, and researchers with the goal of facilitating learning and avoiding, when possible, duplication of efforts in creating teaching material. Contributions to expand the educational content of the OSSCAR project are welcome., Comment: 15 pages, 8 figures

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., Comment: 21 pages, 7 figures

We demonstrate the accuracy and efficiency of a recently introduced approach to account for nuclear quantum effects (NQE) in molecular simulations: the adaptive Quantum Thermal Bath (adQTB). In this method, zero point energy is introduced through a generalized Langevin thermostat designed to precisely enforce the quantum fluctuation-dissipation theorem. We propose a refined adQTB algorithm with improved accuracy and we report adQTB simulations of liquid water. Through extensive comparison with reference path integral calculations, we demonstrate that it provides excellent accuracy for a broad range of structural and thermodynamic observables as well as infrared vibrational spectra. The adQTB has a computational cost comparable to classical molecular dynamics, enabling simulations of up to millions of degrees of freedom.

It has recently been pointed out that Hamiltonian particle systems in constant magnetic fields satisfy generalized time-reversal symmetries that enable to prove useful statistical relationships based on equilibrium phase-space probability distributions without the need to invert, as commonly considered necessary, the magnetic field. Among these relations, that hold without need of Casimir modifications, one finds the standard linear response Green-Kubo relations, and consequently the Onsager reciprocal relations. Going beyond linear response is also possible, for instance in terms of transient and steady state Fluctuation Relations (FRs). Here we highlight how the generalized time-reversal symmetries ensure that the (transient) FRs theory directly applies also for systems in external magnetic fields. Furthermore we show that transient FR can indeed be verified in nonequilibrium molecular dynamics simulations, for systems subjected to magnetic and electric fields, which are thermostatted \`a la Nos\'e-Hoover. The result is nontrivial because, since it is not immediate within which sizes and time scales the effects can actually be observable, it is not obvious what one may obtain by real molecular dynamics simulations., Comment: Submitted to Entropy, 17 pages, 4 figures, 1 table

The validity of the Fluctuation Relations (FR) for systems in a constant magnetic field is investigated. Recently introduced time-reversal symmetries that hold in presence of static electric and magnetic fields and of deterministic thermostats are used to prove the transient FR without invoking, as commonly done, inversion of the magnetic field. Steady-state FR are also derived, under the t-mixing condition. These results extend the predictive power of important statistical mechanics relations. We illustrate this via the non-linear response for the cumulants of the dissipation, showing how the new FR enable to determine analytically null cumulants also for systems in a single magnetic field., Comment: 1 figure, added references

Classical molecular dynamics simulations have recently become a standard tool for the study of electrochemical systems. State-of-the-art approaches represent the electrodes as perfect conductors, modelling their responses to the charge distribution of electrolytes via the so-called fluctuating charge model. These fluctuating charges are additional degrees of freedom that, in a Born-Oppenheimer spirit, adapt instantaneously to changes in the environment to keep each electrode at a constant potential. Here we show that this model can be treated in the framework of constrained molecular dynamics, leading to a symplectic and time-reversible algorithm for the evolution of all the degrees of freedom of the system. The computational cost and the accuracy of the new method are similar to current alternative implementations of the model. The advantage lies in the accuracy and long term stability guaranteed by the formal properties of the algorithm and in the possibility to systematically introduce additional kinematic conditions of arbitrary number and form. We illustrate the performance of the constrained dynamics approach by enforcing the electroneutrality of the electrodes in a simple capacitor consisting of two graphite electrodes separated by a slab of liquid water.

In recent work [Coretti et al., The Journal of Chemical Physics, 2018, 149, 191102], a new algorithm to solve numerically the dynamics of the shell model for polarization was presented. The approach, broadly applicable to systems involving adiabatically separated dynamical variables, employs constrained molecular dynamics to strictly enforce the condition that the force on the fast degrees of freedom, modeled as having zero mass, is null at each time step. The algorithm is symplectic and fully time reversible, and results in stable and efficient propagation. In this paper we complete the discussion of the mechanics of mass zero constrained dynamics by showing how to adapt it to problems where the fast degrees of freedom must satisfy additional conditions. This extension includes, in particular, the important case of first principles molecular dynamics. We then consider the statistical mechanics of the mass zero constrained dynamical system demonstrating that the marginal probability sampled by the dynamics in the physical phase space recovers the form of the Born-Oppenheimer probability density., Comment: 10 pages, no figures, submitted to PCCP (Physical Chemistry Chemical Physics)

We revisit the statistical mechanics of charge fluctuations in capacitors. In constant-potential classical molecular simulations, the atomic charge of electrode atoms are treated as additional degrees of freedom which evolve in time so as to satisfy the constraint of fixed electrostatic potential for each configuration of the electrolyte. The present work clarifies the role of the overall electroneutrality constraint, as well as the link between the averages computed within the Born-Oppenheimer approximation and that of the full constant-potential ensemble. This allows us in particular to derive a complete fluctuation-dissipation relation for the differential capacitance, that includes a contribution from the charge fluctuations (around the charges satisfying the constant-potential and electroneutrality constraints) also present in the absence of an electrolyte. We provide a simple expression for this contribution from the elements of the inverse of the matrix defining the quadratic form of the fluctuating charges in the energy. We then illustrate numerically the validity of our results, and recover the expected result for an empty capacitor with structureless electrodes at large inter-electrode distances. By considering a variety of liquids between graphite electrodes, we confirm that this contribution to the total differential capacitance is small compared to that induced by the thermal fluctuations of the electrolyte.

The quantum-classical Liouville equation provides a description of the dynamics of a quantum subsystem coupled to a classical environment. Representing this equation in the mapping basis leads to a continuous description of discrete quantum states of the subsystem and may provide an alternate route to the construction of simulation schemes. In the mapping basis the quantum-classical Liouville equation consists of a Poisson bracket contribution and a more complex term. By transforming the evolution equation, term-by-term, back to the subsystem basis, the complex term (excess coupling term) is identified as being due to a fraction of the back reaction of the quantum subsystem on its environment. A simple approximation to quantum-classical Liouville dynamics in the mapping basis is obtained by retaining only the Poisson bracket contribution. This approximate mapping form of the quantum-classical Liouville equation can be simulated easily by Newtonian trajectories. We provide an analysis of the effects of neglecting the presence of the excess coupling term on the expectation values of various types of observables. Calculations are carried out on nonadiabatic population and quantum coherence dynamics for curve crossing models. For these observables, the effects of the excess coupling term enter indirectly in the computation and good estimates are obtained with the simplified propagation.