Source author record

Ramses van Zon

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

ResearcherUnclaimed source record

cond-mat.soft cond-mat.stat-mech cond-mat.dis-nn physics.chem-ph quant-ph

Catalog footprint

What is connected

5works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2012arXiv

Free Energy Landscape of Protein-like Chains with Discontinuous Potentials

In this article the configurational space of two simple protein models consisting of polymers composed of a periodic sequence of four different kinds of monomers is studied as a function of temperature. In the protein models, hydrogen bond interactions, electrostatic repulsion, and covalent bond vibrations are modeled by discontinuous step, shoulder and square-well potentials, respectively. The protein-like chains exhibit a secondary alpha helix structure in their folded states at low temperatures, and allow a natural definition of a configuration by considering which beads are bonded. Free energies and entropies of configurations are computed using the parallel tempering method in combination with hybrid Monte Carlo sampling of the canonical ensemble of the discontinuous potential system. The probability of observing the most common configuration is used to analyze the nature of the free energy landscape, and it is found that the model with the least number of possible bonds exhibits a funnel-like free energy landscape at low enough temperature for chains with fewer than 30 beads. For longer proteins, the landscape consists of several minima, where the configuration with the lowest free energy changes significantly by lowering the temperature and the probability of observing the most common configuration never approaches one due to the degeneracy of the lowest accessible potential energy.

preprint2012arXiv

Mapping quantum-classical Liouville equation: projectors and trajectories

The evolution of a mixed quantum-classical system is expressed in the mapping formalism where discrete quantum states are mapped onto oscillator states, resulting in a phase space description of the quantum degrees of freedom. By defining projection operators onto the mapping states corresponding to the physical quantum states, it is shown that the mapping quantum-classical Liouville operator commutes with the projection operator so that the dynamics is confined to the physical space. It is also shown that a trajectory-based solution of this equation can be constructed that requires the simulation of an ensemble of entangled trajectories. An approximation to this evolution equation which retains only the Poisson bracket contribution to the evolution operator does admit a solution in an ensemble of independent trajectories but it is shown that this operator does not commute with the projection operators and the dynamics may take the system outside the physical space. The dynamical instabilities, utility and domain of validity of this approximate dynamics are discussed. The effects are illustrated by simulations on several quantum systems.

preprint2011arXiv

Free Energy Landscape of a Protein-Like Chain in a Fluid with Discontinuous Potentials

The free energy landscape of a protein-like chain in a fluid was studied by combining discontinuous molecular dynamics and parallel tempering. The model protein is a repeating sequence of four different beads, with interactions mimicking those in real proteins. Neighbor distances and angles are restricted to physical ranges and one out of the four kinds of beads can form hydrogen bonds with each other, except if they are too close in the chain. In contrast to earlier studies of this model, an explicit square-well solvent is included. Beads that can form intra-chain hydrogen bonds, can also form (weaker) hydrogen bonds with solvent molecules, while other beads are insoluble. By categorizing the protein configurations according to their intra-chain bonds, one can distinguish unfolded, helical, and collapsed helical structures. Simulations for chains of 15, 20 and 25 beads show that at low temperatures, the most likely structures are helical or collapsed helical, despite the low entropy of these structures. The temperature at which helical structures become dominant is higher than in the absence of a solvent. The cooperative effect of the solvent is attributed to the presence of hydrophobic beads. A phase transition of the solvent prevented the simulations of the 20-bead and 25-bead chains of reaching low enough temperatures to confirm whether the free energy landscape is funnel-shaped, although the results do not contradict that possibility.

preprint2010arXiv

Constructing smooth potentials of mean force, radial, distribution functions and probability densities from sampled data

In this paper a method of obtaining smooth analytical estimates of probability densities, radial distribution functions and potentials of mean force from sampled data in a statistically controlled fashion is presented. The approach is general and can be applied to any density of a single random variable. The method outlined here avoids the use of histograms, which require the specification of a physical parameter (bin size) and tend to give noisy results. The technique is an extension of the Berg-Harris method [B.A. Berg and R.C. Harris, Comp. Phys. Comm. 179, 443 (2008)], which is typically inaccurate for radial distribution functions and potentials of mean force due to a non-uniform Jacobian factor. In addition, the standard method often requires a large number of Fourier modes to represent radial distribution functions, which tends to lead to oscillatory fits. It is shown that the issues of poor sampling due to a Jacobian factor can be resolved using a biased resampling scheme, while the requirement of a large number of Fourier modes is mitigated through an automated piecewise construction approach. The method is demonstrated by analyzing the radial distribution functions in an energy-discretized water model. In addition, the fitting procedure is illustrated on three more applications for which the original Berg-Harris method is not suitable, namely, a random variable with a discontinuous probability density, a density with long tails, and the distribution of the first arrival times of a diffusing particle to a sphere, which has both long tails and short-time structure. In all cases, the resampled, piecewise analytical fit outperforms the histogram and the original Berg-Harris method.

preprint2004arXiv

Glassy dynamics and domains: exact results for the East model

A general matrix-based scheme for analyzing the long-time dynamics in kinetically constrained models such as the East model is presented. The treatment developed here is motivated by the expectation that slowly-relaxing spin domains of arbitrary size govern the highly cooperative events that lead to spin relaxation at long times. To account for the role of large spin domains in the dynamics, a complete basis expressed in terms of domains of all sizes is introduced. It is first demonstrated that accounting for single domains of all possible sizes leads to a simple analytical result for the two-time single-spin correlation function in the East model that is in excellent quantitative agreement with simulation data for equilibrium spin up density values c greater or equal to 0.6. It is then shown that including also two neighboring domains leads to a closed expression that describes the slow relaxation of the system down to c approximately 0.3. Ingredients of generalizing the method to lower values of c are also provided, as well as to other models. The main advantage of this approach is that it gives explicit analytical results and that it requires neither an arbitrary closure for the memory kernel nor the construction of an irreducible memory kernel. It also allows one to calculate quantities that measure heterogeneity in the same framework, as is illustrated on the neighbor-pair correlation function and the distribution of relaxation times.