Internal coordinate density of state from molecular dynamics simulation

Pure and Sc-doped diopside (CaMgSi2O6) vibrational spectra: modelling and experiments

We investigated the structure of pure and Sc-doped synthetic diopside (a monoclinic single-chain silicate nominally CaMgSi2O6); in Sc-doped diopside, Sc3+ substitutes Mg2+ in the structure and, to achieve charge balance, vacancies form at the expense of Ca2+. We compared the structure obtained from ab initio modelling techniques at the density functional theory (DFT) level with the structure solved by employing single crystal X-ray diffraction. Furthermore, we compared IR and Raman spectroscopy experiments with vibrational density of states (VDOS) calculated from the Fourier transform of the velocity autocorrelation function obtained using ab initio (DFT) molecular dynamics simulations. In this framework, we developed a computational tool to assign the vibrational mode associated with a specific frequency. This method consists of projecting velocities along a specific set of internal coordinates such as stretching or bending, in cases involving defects or vacancies, to calculate a partial VDOS (pVDOS) that takes into account only the vibrations associated with selected internal modes, aiding the interpretation of the total VDOS and the experimental spectra in a relevant way. The computed data were validated with the experiments and we observed that doping the diopside structure with Sc produces peak broadening and the occurrence of new peaks in the Raman spectra and that site vacancies are associated with the nearby Sc site. The present work constitutes an interesting starting point to exploit the calculated VDOS/pVDOS to characterize experimental vibrational spectra of complex systems containing local vacancies, substitutions or defects as the Sc-doped diopside.

High-performance transformation of protein structure representation from internal to Cartesian coordinates

We present a highly parallel algorithm to convert internal coordinates of a polymeric molecule into Cartesian coordinates. Traditionally, converting the structures of polymers (e.g., proteins) from internal to Cartesian coordinates has been performed serially, due to an inherent linear dependency along the polymer chain. We show this dependency can be removed using a tree-based concatenation of coordinate transforms between segments, and then parallelized efficiently on graphics processing units (GPUs). The conversion algorithm is applicable to protein engineering and fitting protein structures to experimental data, and we observe an order of magnitude speedup using parallel processing on a GPU compared to serial execution on a CPU.

Molecular mechanism of water permeation in a helium impermeable graphene and graphene oxide membrane

Layers of graphene oxide (GO) are found to be good for the permeation of water but not for helium (Science, 2012, 335(6067), 442-444) suggesting that the GO layers are dynamic in the formation of a permeation route depending on the environment they are in (i.e., water or helium). To probe the microscopic origin of this observation we calculate the potential of mean force (PMF) of GO sheets (with oxidized and reduced parts), with the inter-planar distance as a reaction coordinate in helium and water. Our PMF calculation shows that the equilibrium interlayer distance between the oxidized part of the GO sheets in helium is at 4.8 Å leaving no space for helium permeation. In contrast, the PMF of the oxidized part of the GO in water shows two minima, one at 4.8 Å and another at 6.8 Å, corresponding to no water and a water filled region, thus giving rise to a permeation path. The increased electrostatic interaction between water with the oxidized part of the sheet helps the sheet open up and pushes water inside. Based on the entropy calculations for water trapped between graphene sheets and oxidized graphene sheets at different inter-sheet spacings, we also show the thermodynamics of filling.