87 research outputs found
CHAMPION: Chalmers Hierarchical Atomic, Molecular, Polymeric & Ionic Analysis Toolkit
We present CHAMPION: a software developed to automatically detect
time-dependent bonds between atoms based on their dynamics, classify the local
graph topology around them, and analyze the physicochemical properties of these
topologies by statistical physics. In stark contrast to methodologies where
bonds are detected based on static conditions such as cut-off distances,
CHAMPION considers pairs of atoms to be bound only if they move together and
act as a bound pair over time. Furthermore, the time-dependent global bond
graph is possible to split into dynamically shifting connected components or
subgraphs around a certain chemical motif and thereby allow the physicochemical
properties of each such topology to be analyzed by statistical physics.
Applicable to condensed matter and liquids in general, and electrolytes in
particular, this allows both quantitative and qualitative descriptions of local
structure, as well as dynamical processes such as speciation and diffusion. We
present here a detailed overview of CHAMPION, including its underlying
methodology, implementation and capabilities.Comment: 11 pages, 8 figure
Non-equilibrium statistical mechanics: From a paradigmatic model to biological transport
Unlike equilibrium statistical mechanics, with its well-established
foundations, a similar widely-accepted framework for non-equilibrium
statistical mechanics (NESM) remains elusive. Here, we review some of the many
recent activities on NESM, focusing on some of the fundamental issues and
general aspects. Using the language of stochastic Markov processes, we
emphasize general properties of the evolution of configurational probabilities,
as described by master equations. Of particular interest are systems in which
the dynamics violate detailed balance, since such systems serve to model a wide
variety of phenomena in nature. We next review two distinct approaches for
investigating such problems. One approach focuses on models sufficiently simple
to allow us to find exact, analytic, non-trivial results. We provide detailed
mathematical analyses of a one-dimensional continuous-time lattice gas, the
totally asymmetric exclusion process (TASEP). It is regarded as a paradigmatic
model for NESM, much like the role the Ising model played for equilibrium
statistical mechanics. It is also the starting point for the second approach,
which attempts to include more realistic ingredients in order to be more
applicable to systems in nature. Restricting ourselves to the area of
biophysics and cellular biology, we review a number of models that are relevant
for transport phenomena. Successes and limitations of these simple models are
also highlighted.Comment: 72 pages, 18 figures, Accepted to: Reports on Progress in Physic
Integrating sequence and array data to create an improved 1000 Genomes Project haplotype reference panel
A major use of the 1000 Genomes Project (1000GP) data is genotype imputation in genome-wide association studies (GWAS). Here we develop a method to estimate haplotypes from low-coverage sequencing data that can take advantage of single-nucleotide polymorphism (SNP) microarray genotypes on the same samples. First the SNP array data are phased to build a backbone (or 'scaffold') of haplotypes across each chromosome. We then phase the sequence data 'onto' this haplotype scaffold. This approach can take advantage of relatedness between sequenced and non-sequenced samples to improve accuracy. We use this method to create a new 1000GP haplotype reference set for use by the human genetic community. Using a set of validation genotypes at SNP and bi-allelic indels we show that these haplotypes have lower genotype discordance and improved imputation performance into downstream GWAS samples, especially at low-frequency variants. © 2014 Macmillan Publishers Limited. All rights reserved
Dynamics and Energetics of Permeation Through Aquaporins. What Do We Learn from Molecular Dynamics Simulations?
A generalized reactive force field for nonlinear hydrogen bonds : hydrogen dynamics and transfer in malonaldehyde
Using molecular dynamics (MD) simulations, the spectroscopy and dynamics of malonaldehyde is investigated. To this end, the recently proposed molecular mechanics with proton transfer (MMPT) potential is generalized to nonlinear hydrogen bonds. The calculated properties for malonaldehyde in both gas and condensed phases, including equilibrium geometries, infrared spectra, tunneling splittings, and hydrogen transfer rates, compare well with previous experimental and computational works. In particular, by using a harmonic bath averaged (HBA) Hamiltonian, which is based on a reaction path Hamiltonian, it is possible to estimate the tunneling splitting in an efficient manner. It is found that a zero point corrected barrier of 6.7 kcal/mol and effective masses of 1.234 (i.e., 23.4% larger than the mass of a physical H-atom) and 1.117 (for the physical D-atom) are consistent with the measured splittings of 21.6 and 2.9 cm(-1), respectively. The HBA Hamiltonian also yields a pair of hydrogen transfer fundamentals at 1573 and 1267 cm(-1), similar to results obtained with a reaction surface Hamiltonian on a MP2/6-31G(**) potential energy surface. This amounts to a substantial redshift of more than 1000 cm(-1) which can be rationalized by comparison with weakly (HCO(+): rare gas) and strongly (H(2)O-H(+)-OH(2)) proton-bound systems. Hydrogen transfer rates in vacuum and water were determined from the validated MMPT potential and it is found that the solvent enhances the rate by a factor of 5 at 300 K. The rates of 2.4/ns and 10/ns are commensurate with previous density functional tight binding ab initio MD studies
- âŠ