Source author record

Giovanni Ciccotti

Giovanni Ciccotti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2020arXiv

Adiabatic motion and statistical mechanics via mass zero constrained dynamics

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.

preprint2016arXiv

On the Properties of a Bundle of Flexible Actin Filaments in an Optical Trap

We establish the Statistical Mechanics framework for a bundle of Nf living and uncrosslinked actin filaments in a supercritical solution of free monomers pressing against a mobile wall. The filaments are anchored normally to a fixed planar surface at one of their ends and, because of their limited flexibility, they grow almost parallel to each other. Their growing ends hit a moving obstacle, depicted as a second planar wall, parallel to the previous one and subjected to a harmonic compressive force. The force constant is denoted as trap strength while the distance between the two walls as trap length to make contact with the experimental optical trap apparatus. For an ideal solution of reactive filaments and free monomers at fixed free monomers chemical potential, we obtain the general expression for the grand potential from which we derive averages and distributions of relevant physical quantities, namely the obstacle position, the bundle polymerization force and the number of filaments in direct contact with the wall. The grafted living filaments are modeled as discrete Wormlike chains, with Factin persistence length, subject to discrete contour length variations to model single monomer (de)polymerization steps. Rigid filaments, either isolated or in bundles, all provide average values of the stalling force in agreement with Hill's predictions, independent of the average trap length. Flexible filaments instead, for values of the trap strength suitable to prevent their lateral escape, provide an average bundle force and an average trap length slightly larger than the corresponding rigid cases (few percents). Still the stalling force remains nearly independent on the average trap length, but results from the product of two strongly L dependent contributions: the fraction of touching filaments and the single filament buckling force.

preprint2015arXiv

A semi-flexible model prediction for the polymerization force exerted by a living F-actin filament on a fixed wall

We consider a single living semi-flexible filament with persistence length l_p in chemical equilibrium with a solution of free monomers at fixed monomer chemical potential mu_1 and fixed temperature T. While one end of the filament is chemically active with single monomer (de)polymerization steps, the other end is grafted normally to a rigid wall to mimick a rigid network from which the filament under consideration emerges. A second rigid wall, parallel to the grafting wall, is fixed at distance L<<l_p from the filament seed. In supercritical conditions the filament tends to grow and impinges onto the second surface which, in suitable conditions (non-escaping filament regime) stops the filament growth. We first establish the grand-potential and derive some general properties, in particular the filament size distribution and the force exerted by the living filament on the obstacle wall. We apply this formalism to the semi-flexible, living, discrete Wormlike chain (d-WLC) model with step size d and persistence length l_p, hitting a hard wall. By original Monte-Carlo calculations we justify the use of the weak bending universal expressions of Gholami et al. (Phys.Rev.E. 74,(2006), 041803) over the whole non escaping filament regime. Employing this universal form for living filaments, we find that the average force exerted by a living filament on a wall at distance L is in practice L independent and very close to the value predicted by Hill, his expression being strictly valid in the rigid filament limit. The average filament force results from the product of the cumulative size fraction x, where the filament is in contact with the wall, times the buckling force on a filament of size L_c ~ L. We discuss several consequences of the L independence of the stalling force for our specific filament model.

preprint2015arXiv

Maximum probability domains for the analysis of the microscopic structure of liquids

We introduce the concept of maximum probability domains (MPDs), developed in the context of the analysis of electronic densities, in the study of the microscopic spatial structures of liquids. The idea of locating a particle in a three dimensional region, by determining the domain where the probability of finding that, and only that, particle is maximum, gives an interesting characterization of the local structure of the liquid. The optimization procedure, required for the search of the domain of maximum probability, is carried out by the implementation of the level set method. Results for a couple of case studies are presented, to illustrate the structure of liquid water at ambient conditions and upon increasing pressure from the point of view of MPDs and to compare the information encoded in the solvation shells of sodium in water with, once again, that extracted from the MPDs.

preprint2014arXiv

Exploring the conformational dynamics of alanine dipeptide in solution subjected to an external electric field: A nonequilibrium molecular dynamics simulation

In this paper, we investigate the conformational dynamics of alanine dipeptide under an external electric field by nonequilibrium molecular dynamics simulation. We consider the case of a constant and of an oscillatory field. In this context we propose a procedure to implement the temperature control, which removes the irrelevant thermal effects of the field. For the constant field different time-scales are identified in the conformational, dipole moment, and orientational dynamics. Moreover, we prove that the solvent structure only marginally changes when the external field is switched on. In the case of oscillatory field, the conformational changes are shown to be as strong as in the previous case, and non-trivial nonequilibrium circular paths in the conformation space are revealed by calculating the integrated net probability fluxes.

preprint2014arXiv

Relaxation of a steep density gradient in a simple fluid: comparison between atomistic and continuum modeling

We compare dynamical nonequilibrium molecular dynamics and continuum simulations of the dynamics of relaxation of a fluid system characterized by a non uniform density profile. Results match quite well as long as the lengthscale of density nonuniformities are greater than the molecular scale (10 times the molecular size). In presence of molecular scale features some of the continuum fields (e.g. density and momentum) are in good agreement with atomistic counterparts, but are smoother. On the contrary, other fields, such at the temperature field, present very large difference with respect to reference (atomistic) ones. This is due to the limited accuracy of some of the empirical relations used in continuum models, the equation of state of the fluid in the example considered.

preprint2014arXiv

Spontaneous chiral symmetry breaking in model bacterial suspensions

Chiral symmetry breaking is ubiquitous in biological systems, from DNA to bacterial suspensions. A key unresolved problem is how chiral structures may spontaneously emerge from achiral interactions. We study a simple model of bacterial suspensions in three dimensions that effectively incorporates active motion and hydrodynamic interactions. We perform large-scale molecular dynamics simulations (up to $10^6$ particles) and describe stable (or long-lived metastable) collective states that exhibit chiral organization although the interactions are achiral. We elucidate under which conditions these chiral states will emerge and grow to large scales. We also study a related equilibrium model that clarifies the role of orientational fluctuations.

preprint2012arXiv

Adaptive resolution molecular dynamics simulation through coupling to an internal particle reservoir

For simulation studies of (macro) molecular liquids it would be of significant interest to be able to adjust or increase the level of resolution within one region of space, while allowing for the free exchange of molecules between open regions of different resolution or representation. We generalize the adaptive resolution idea and suggest an interpretation in terms of an effective generalized grand canonical approach. The method is applied to liquid water at ambient conditions.

preprint2012arXiv

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.

preprint2011arXiv

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.

preprint2011arXiv

Hydrodynamics from dynamical non-equilibrium MD

We review a dynamical approach to non-equilibrium MD (D-NEMD). We show how, using a proper simulation setup, is possible to treat interesting cases in which the initial condition is a stationary non-equilibrium state produced by a suitable dynamical system. We then extend the class of non-equilibrium phenomena that can be studied by atomistic simulations to the case of complex initial conditions consisting in assigning a macroscopic value of a scalar or vector observable or a field. We illustrate the functioning of this method by applying it to the relaxation of an interface between two immiscible liquids. We have shown that our method generate unbiased results while this might not be the case for the often used short time average approach.