Enhanced Sampling in the Well-Tempered Ensemble

We introduce the well-tempered ensemble (WTE) which is the biased ensemble sampled by well-tempered metadynamics when the energy is used as collective variable. WTE can be designed so as to have approximately the same average energy as the canonical ensemble but much larger fluctuations. These two properties lead to an extremely fast exploration of phase space. An even greater efficiency is obtained when WTE is combined with parallel tempering. Unbiased Boltzmann averages are computed on the fly by a recently developed reweighting method [M. Bonomi et al. J. Comput. Chem. 30, 1615 (2009)]. We apply WTE and its parallel tempering variant to the 2d Ising model and to a Go-model of HIV protease, demonstrating in these two representative cases that convergence is accelerated by orders of magnitude.Comment: 7 pages, 5 figure

Entropy based fingerprint for local crystalline order

We introduce a new fingerprint that allows distinguishing between liquid-like and solid-like atomic environments. This fingerprint is based on an approximate expression for the entropy projected on individual atoms. When combined with a local enthalpy, this fingerprint acquires an even finer resolution and it is capable of discriminating between different crystal structures.Comment: 6 pages, 4 figure

Escaping free-energy minima

We introduce a novel and powerful method for exploring the properties of the multidimensional free energy surfaces of complex many-body systems by means of a coarse-grained non-Markovian dynamics in the space defined by a few collective coordinates.A characteristic feature of this dynamics is the presence of a history-dependent potential term that, in time, fills the minima in the free energy surface, allowing the efficient exploration and accurate determination of the free energy surface as a function of the collective coordinates. We demonstrate the usefulness of this approach in the case of the dissociation of a NaCl molecule in water and in the study of the conformational changes of a dialanine in solution.Comment: 3 figure

Evolution of the structure of amorphous ice - from low-density amorphous (LDA) through high-density amorphous (HDA) to very high-density amorphous (VHDA) ice

We report results of molecular dynamics simulations of amorphous ice for pressures up to 22.5 kbar. The high-density amorphous ice (HDA) as prepared by pressure-induced amorphization of Ih ice at T=80 K is annealed to T=170 K at various pressures to allow for relaxation. Upon increase of pressure, relaxed amorphous ice undergoes a pronounced change of structure, ranging from the low-density amorphous ice (LDA) at p=0, through a continuum of HDA states to the limiting very high-density amorphous ice (VHDA) regime above 10 kbar. The main part of the overall structural change takes place within the HDA megabasin, which includes a variety of structures with quite different local and medium-range order as well as network topology and spans a broad range of densities. The VHDA represents the limit to densification by adapting the hydrogen-bonded network topology, without creating interpenetrating networks. The connection between structure and metastability of various forms upon decompression and heating is studied and discussed. We also discuss the analogy with amorphous and crystalline silica. Finally, some conclusions concerning the relation between amorphous ice and supercooled water are drawn.Comment: 11 pages, 12 postscript figures. To be published in The Journal of Chemical Physic

Coexistence of tetrahedral and octahedral-like sites in amorphous phase change materials

Chalcogenide alloys are materials of interest for optical recording and non-volatile memories. We perform ab-initio molecular dynamics simulations aiming at shading light onto the structure of amorphous Ge2Sb2Te5 (GST), the prototypical material in this class. First principles simulations show that amorphous GST obtained by quenching from the liquid phase displays two types of short range order. One third of Ge atoms are in a tetrahedral environment while the remaining Ge, Sb and Te atoms display a defective octahedral environment, reminiscent of cubic crystalline GST.Comment: 6 pages, 10 figures including supplementary materia

Enhancing entropy and enthalpy fluctuations to drive crystallization in atomistic simulations

Crystallization is a process of great practical relevance in which rare but crucial fluctuations lead to the formation of a solid phase starting from the liquid. Like in all first order first transitions there is an interplay between enthalpy and entropy. Based on this idea, to drive crystallization in molecular simulations, we introduce two collective variables, one enthalpic and the other entropic. Defined in this way, these collective variables do not prejudge the structure the system is going to crystallize into. We show the usefulness of this approach by studying the case of sodium and aluminum that crystallize in the bcc and fcc crystalline structure, respectively. Using these two generic collective variables, we perform variationally enhanced sampling and well tempered metadynamics simulations, and find that the systems transform spontaneously and reversibly between the liquid and the solid phases.Comment: 4 pages, 2 figure

Naphthalene crystal shape prediction from molecular dynamics simulations

We used molecular dynamics simulations to predict the steady state crystal shape of naphthalene grown from ethanol solution. The simulations were performed at constant supersaturation by utilizing a recently proposed algorithm [Perego et al., J. Chem. Phys., 142, 2015, 144113]. To bring the crystal growth within the timescale of a molecular dynamics simulation we applied Well-Tempered Metadynamics with a spatially constrained collective variable, which focuses the sampling on the growing layer. We estimated that the resulting steady state crystal shape corresponds to a rhombic prism, which is in line with experiments. Further, we observed that at the investigated supersaturations, the $\{00\bar{1}\}$ face grows in a two step two dimensional nucleation mechanism while the considerably faster growing faces $\{1\bar{1}0\}$ and $\{20\bar{1}\}$ grow new layers with a one step two dimensional nucleation mechanism

The Ciona intestinalis immune-related galectin genes (CiLgals-a and CiLgals-b) are expressed by the gastric epithelium

The transcription of two Ciona intestinalis galectin genes (CiLgals-a and CiLgalseb) is uparegulated by LPS in the pharynxis (hemocytes, vessel epithelium, endostilar zones) which is retained the main organ of the immunity. In this ascidian, for the first time we show, by immunohistochemistry and in situ hybridization methods, that these two immune-related genes are expressed in the gastric epithelium of na\uefve ascidians, whereas the galectins appear to be only contained in the intestine columnar epithelium. In addition, according to previous results on the pharynx, the genes are also expressed and galectins produced by hemocytes scattered in the connective tissue surrounding the gut. The genes expression and galectin localization in several tissues, including the previous findings on the transcription upregulation, the constitutive expression of these genes by endostylar zones and by the gastric epithelium suggest a potential multifunctional role of these galectins. In this respect, it is of interest to define where the CiLgals are normally found as related to the tissue functions. Such an approach should be a starting point for further investigations