Testing of the GROMOS Force-Field Parameter Set 54A8: Structural Properties of Electrolyte Solutions, Lipid Bilayers, and Proteins

We demonstrate the accuracy of the hypernetted chain closure and of the mean-field approximation for the calculation of the fluid-state properties of systems interacting by means of bounded and positive-definite pair potentials with oscillating Fourier transforms. Subsequently, we prove the validity of a bilinear, random-phase density functional for arbitrary inhomogeneous phases of the same systems. On the basis of this functional, we calculate analytically the freezing parameters of the latter. We demonstrate explicitly that the stable crystals feature a lattice constant that is independent of density and whose value is dictated by the position of the negative minimum of the Fourier transform of the pair potential. This property is equivalent with the existence of clusters, whose population scales proportionally to the density. We establish that regardless of the form of the interaction potential and of the location on the freezing line, all cluster crystals have a universal Lindemann ratio L = 0.189 at freezing. We further make an explicit link between the aforementioned density functional and the harmonic theory of crystals. This allows us to establish an equivalence between the emergence of clusters and the existence of negative Fourier components of the interaction potential. Finally, we make a connection between the class of models at hand and the system of infinite-dimensional hard spheres, when the limits of interaction steepness and space dimension are both taken to infinity in a particularly described fashion. Show more

A coarse-grained simulation method to predict NMR spectra of ions diffusing in porous carbons is proposed. The coarse-grained model uses input from molecular dynamics simulations such as the free-energy profile for ionic adsorption, and density-functional theory calculations are used to predict the NMR chemical shift of the diffusing ions. The approach is used to compute NMR spectra of ions in slit pores with pore widths ranging from 2 to 10 nm. As diffusion inside pores is fast, the NMR spectrum of an ion trapped in a single mesopore will be a sharp peak with a pore size dependent chemical shift. To account for the experimentally observed NMR line shapes, our simulations must model the relatively slow exchange between different pores. We show that the computed NMR line shapes depend on both the pore size distribution and the spatial arrangement of the pores. The technique presented in this work provides a tool to extract information about the spatial distribution of pore sizes from NMR spectra. Such information is diffcult to obtain from other characterisation techniques. Show more