The breakthrough of a quantum chemist by classical dynamics: Martin Karplus and the birth of computer simulations of chemical reactions

Abstract1964–1965 was an early, crucial period in Martin Karplus' research—a time when, rather unexpectedly, he approached the problem of reactive collisions using a quasiclassical approximation with the aid of computer technologies. This marked a substantial departure from the quantum-chemical studies of nuclear magnetic resonance that had, until then, dominated his work. The historical perspective outlined by George Schatz, as well Karplus' own biography, partly frames the contours of this remarkable period in the history of theoretical chemistry. Yet, the available historical literature is not sufficiently complete to allow us to understand Karplus' transition from nuclear magnetic resonance to reaction dynamics. In this article, we discuss the intellectual ground on which Karplus operated around 1964, further commenting on the relevance of his quantum and quasiclassical studies and pondering how Karplus' approach eventually led to his interest in the simulation of complex biomolecules

Combining rare events techniques: phase change in Si nanoparticles

We introduce a combined Restrained MD/Parallel Tempering approach to study the difference in free energy as a function of a set of collective variables between two states in presence of unknown slow degrees of freedom. We applied this method to study the relative stability of the amorphous vs crystalline nanoparticles of size ranging between 0.8 and 1.8 nm as a function of the temperature. We found that, at variance with bulk systems, at low T small nanoparticles are amorphous and undergo an amorphous-to-crystalline phase transition at higher T. On the contrary, large nanoparticles recover the bulk-like behavior: crystalline at low $T$ and amorphous at high T

Bulk viscosity of the Lennard-Jones system at the triple point by dynamical Non Equilibrium Molecular Dynamics

Non-equilibrium Molecular Dynamics (NEMD) calculations of the bulk viscosity of the triple point Lennard-Jones fluid are performed with the aim of investigating the origin of the observed disagreement between Green-Kubo estimates and previous NEMD data. We show that a careful application of the Doll's perturbation field, the dynamical NEMD method, the instantaneous form of the perturbation and the "subtraction technique" provides a NEMD estimate of the bulk viscosity at zero field in full agreement with the value obtained by the Green-Kubo formula. As previously reported for the shear viscosity, we find that the bulk viscosity exhibits a large linear regime with the field intensity which confirms the Lennard-Jones fluid as a genuine Newtonian fluid even at triple point.Comment: 27 pages, 11 figure

Partitioning a macroscopic system into independent subsystems

We discuss the problem of partitioning a macroscopic system into a collection of independent subsystems. The partitioning of a system into replica-like subsystems is nowadays a subject of major interest in several field of theoretical and applied physics, and the thermodynamic approach currently favoured by practitioners is based on a phenomenological definition of an interface energy associated with the partition, due to a lack of easily computable expressions for a microscopic (i.e.~particle-based) interface energy. In this article, we outline a general approach to derive sharp and computable bounds for the interface free energy in terms of microscopic statistical quantities. We discuss potential applications in nanothermodynamics and outline possible future directions.Comment: This is an author-created, un-copyedited version of an article accepted for publication in JSTA

Approximating Time-Dependent Quantum Statistical Properties

Computing quantum dynamics in condensed matter systems is an open challenge due to the exponential scaling of exact algorithms with the number of degrees of freedom. Current methods try to reduce the cost of the calculation using classical dynamics as the key ingredient of approximations of the quantum time evolution. Two main approaches exist, quantum classical and semi-classical, but they suffer from various difficulties, in particular when trying to go beyond the classical approximation. It may then be useful to reconsider the problem focusing on statistical time-dependent averages rather than directly on the dynamics. In this paper, we discuss a recently developed scheme for calculating symmetrized correlation functions. In this scheme, the full (complex time) evolution is broken into segments alternating thermal and real-time propagation, and the latter is reduced to classical dynamics via a linearization approximation. Increasing the number of segments systematically improves the result with respect to full classical dynamics, but at a cost which is still prohibitive. If only one segment is considered, a cumulant expansion can be used to obtain a computationally efficient algorithm, which has proven accurate for condensed phase systems in moderately quantum regimes. This scheme is summarized in the second part of the paper. We conclude by outlining how the cumulant expansion formally provides a way to improve convergence also for more than one segment. Future work will focus on testing the numerical performance of this extension and, more importantly, on investigating the limit for the number of segments that goes to infinity of the approximate expression for the symmetrized correlation function to assess formally its convergence to the exact result

Jarzynski on work and free energy relations: The case of variable volume

Derivations of the Jarzynski equality (JE) appear to be quite general, and applicable to any particle system, whether deterministic or stochastic, under equally general perturbations of an initial equilibrium state at given temperature T. At the same time, the definitions of the quantities appearing in the JE, in particular the work, have been questioned. Answers have been given, but a deeper understanding of the range of phenomena to which the JE applies is necessary, both conceptually and in order to interpret the experiments in which it is used. In fact, domains in which the JE is not applicable have been identified. To clarify the issue, we scrutinize the applicability of the JE to a Hamiltonian particle system in a variable volume. We find that, in this case, the standard interpretation of the terms appearing in the JE is not adequate

Dynamical Non-Equilibrium Molecular Dynamics

In this review, we discuss the Dynamical approach to Non-Equilibrium Molecular Dynamics (D-NEMD), which extends stationary NEMD to time-dependent situations, be they responses or relaxations. Based on the original Onsager regression hypothesis, implemented in the nineteen-seventies by Ciccotti, Jacucci and MacDonald, the approach permits one to separate the problem of dynamical evolution from the problem of sampling the initial condition. D-NEMD provides the theoretical framework to compute time-dependent macroscopic dynamical behaviors by averaging on a large sample of non-equilibrium trajectories starting from an ensemble of initial conditions generated from a suitable (equilibrium or non-equilibrium) distribution at time zero. We also discuss how to generate a large class of initial distributions. The same approach applies also to the calculation of the rate constants of activated processes. The range of problems treatable by this method is illustrated by discussing applications to a few key hydrodynamic processes (the “classical” flow under shear, the formation of convective cells and the relaxation of an interface between two immiscible liquids)

On the force-velocity relationship of a bundle of rigid bio-filaments

In various cellular processes, bio-filaments like F-actin and F-tubulin are able to exploit chemical energy associated with polymerization to perform mechanical work against an obstacle loaded with an external force. The force-velocity relationship quantitatively summarizes the nature of this process. By a stochastic dynamical model, we give, together with the evolution of a staggered bundle of Nfrigid living filaments facing a loaded wall, the corresponding force-velocity relationship. We compute the evolution of the model in the infinite wall diffusion limit and in supercritical conditions (monomer density reduced by critical density ρ^1>1), and we show that this solution remains valid for moderate non-zero values of the ratio between the wall diffusion and the chemical time scales. We consider two classical protocols: the bundle is opposed either to a constant load or to an optical trap setup, characterized by a harmonic restoring force. The constant load case leads, for each F value, to a stationary velocity Vstat(F;Nf,ρ^1) after a relaxation with characteristic time τmicro(F). When the bundle (initially taken as an assembly of filament seeds) is subjected to a harmonic restoring force (optical trap load), the bundle elongates and the load increases up to stalling over a characteristic time τOT. Extracted from this single experiment, the force-velocity VOT(F;Nf,ρ^1) curve is found to coincide with Vstat(F;Nf,ρ^1), except at low loads. We show that this result follows from the adiabatic separation between τmicroand τOT, i.e., τmicro≈ τOT

An observable for vacancy characterization and diffusion in crystals

To locate the position and characterize the dynamics of a vacancy in a crystal, we propose to represent it by the ground state density of a quantum probe quasi-particle for the Hamiltonian associated to the potential energy field generated by the atoms in the sample. In this description, the h^2/2mu coefficient of the kinetic energy term is a tunable parameter controlling the density localization in the regions of relevant minima of the potential energy field. Based on this description, we derive a set of collective variables that we use in rare event simulations to identify some of the vacancy diffusion paths in a 2D crystal. Our simulations reveal, in addition to the simple and expected nearest neighbor hopping path, a collective migration mechanism of the vacancy. This mechanism involves several lattice sites and produces a long range migration of the vacancy. Finally, we also observed a vacancy induced crystal reorientation process