Polymer adsorption on curved surfaces

The conformational behavior of a coarse-grained finite polymer chain near an attractive spherical surface was investigated by means of multicanonical Monte Carlo computer simulations. In a detailed analysis of canonical equilibrium data over a wide range of sphere radius and temperature, we have constructed entire phase diagrams both for nongrafted and end-grafted polymers. For the identification of the conformational phases, we have calculated several energetic and structural observables such as gyration tensor based shape parameters and their fluctuations by canonical statistical analysis. Despite the simplicity of our model, it qualitatively represents in the considered parameter range real systems that are studied in experiments. The work discussed here could have experimental implications from protein-ligand interactions to designing nanosmart materials.

Subphase transitions in first-order aggregation processes

In this paper, we investigate the properties of aggregation transitions in the context of generic coarse-grained homopolymer systems. By means of parallel replica-exchange Monte Carlo methods, we perform extensive simulations of systems consisting of up to 20 individual oligomer chains with five monomers each. Using the tools of the versatile microcanonical inflection-point analysis, we show that the aggregation transition is a first-order process consisting of a sequence of subtransitions between intermediate structural phases. We unravel the properties of these intermediate phases by collecting and analyzing their individual contributions towards the density of states of the system. The central theme of this systematic study revolves around translational entropy and its role in the striking phenomena of missing intermediate phases. We conclude with a brief discussion of the scaling properties of the transition temperature and the latent heat.

Interlocking order parameter fluctuations in structural transitions between adsorbed polymer phases

By means of contact-density chain-growth simulations of a simple coarse-grained lattice model for a polymer grafted at a solid homogeneous substrate, we investigate the complementary behavior of the numbers of surface-monomer and monomer-monomer contacts under various solvent and thermal conditions. This pair of contact numbers represents an appropriate set of order parameters that enables the distinct discrimination of significantly different compact phases of polymer adsorption. Depending on the transition scenario, these order parameters can interlock in perfect cooperation. The analysis helps understand the transitions from compact filmlike adsorbed polymer conformations into layered morphologies and dissolved adsorbed structures, respectively, in more detail.

Confinement effects upon the separation of structural transitions in linear systems with restricted bond fluctuation ranges

By means of advanced parallel replica-exchange Monte Carlo methods we examine the influence of elasticity and confinement on the structural transitions of linear systems with restricted bonded interaction. For this purpose, we adopt a model for coarse-grained flexible polymers of finite length in the dilute regime. Hyperphase diagrams are constructed using energy-dependent canonical quantities to demonstrate the effects of the changes in the range of the confined interaction on the liquid and solid structural phases. With increasing bonded interaction range we observe the disappearance of the liquid phase and the fusion of the gas-liquid (or Θ) and the liquid-solid transitions. One of the most remarkable features, the liquid-gas transition, changes from second to first order if the confined interaction range exceeds a threshold that separates polymeric from nonpolymeric systems. The notoriously difficult sampling of the entropically suppressed conformations in the region of very strong first-order transitions is improved by using multiple Gaussian modified ensembles.

Interplay of Coil-Globule Transition and Surface Adsorption of a Lattice HP Protein Model

An end-grafted hydrophobic-polar (HP) model protein chain with alternating H and P monomers is studied to examine interactions between the critical adsorption transition due to surface attraction and the collapse transition due to pairwise attractive H–H interactions. We find that the critical adsorption phenomenon can always be observed; however, the critical adsorption temperature T_CAP is influenced by the attractive H–H interactions in some cases. When the collapse temperature T_c is lower than T_CAP, the critical adsorption of the HP chain is similar to that of a homopolymer without intrachain attractions and T_CAP remains unchanged, whereas the collapse transition is suppressed by the adsorption. In contrast, for cases where T_c is close to or higher than T_CAP, T_CAP of the HP chain is increased, indicating that a collapsed chain is more easily adsorbed on the surface. The strength of the H–H attraction also influences the statistical size and shape of the polymer, with strong H–H attractions resulting in adsorbed and collapsed chains adopting two-dimensional, circular conformations.

Adsorption and pattern recognition of polymers at complex surfaces with attractive stripelike motifs

We construct the complete structural phase diagram of polymer adsorption at substrates with attractive stripelike patterns in the parameter space spanned by the adsorption affinity of the stripes and temperature. Results were obtained by extensive generalized-ensemble Monte Carlo simulations of a generic model for the hybrid organic-inorganic system. By comparing with adhesion properties at homogeneous substrates, we find substantial differences in the formation of adsorbed polymer structures if translational invariance at the surface is broken by a regular pattern. Beside a more specific understanding of polymer adsorption processes, our results are potentially relevant for the design of macromolecular pattern recognition devices such as sensors.

Generic folding and transition hierarchies for surface adsorption of hydrophobic-polar lattice model proteins

The thermodynamic behavior and structural properties of hydrophobic-polar (HP) lattice proteins interacting with attractive surfaces are studied by means of Wang-Landau sampling. Three benchmark HP sequences (48mer, 67mer, and 103mer) are considered with different types of surfaces, each of which attract either all monomers, only hydrophobic (H) monomers, or only polar (P) monomers, respectively. The diversity of folding behavior in dependence of surface strength is discussed. Analyzing the combined patterns of various structural observables, such as, e.g., the derivatives of the numbers of surface contacts, together with the specific heat, we are able to identify generic categories of folding and transition hierarchies. We also infer a connection between these transition categories and the relative surface strengths, i.e., the ratio of the surface attractive strength to the interchain attraction among H monomers. The validity of our proposed classification scheme is reinforced by the analysis of additional benchmark sequences. We thus believe that the folding hierarchies and identification scheme are generic for HP proteins interacting with attractive surfaces, regardless of chain length, sequence, or surface attraction.

"Pull moves" for rectangular lattice polymer models are not fully reversible

"Pull moves" is a popular move set for lattice polymer model simulations. We show that the proof given for its reversibility earlier is flawed, and some moves are irreversible, which leads to biases in the parameters estimated from the simulations. We show how to make the move set fully reversible.

Selective adsorption of lattice peptides on patterned surfaces

To study the adsorption of individual peptides in implicit solvent, we propose a version of the Wang-Landau Monte Carlo algorithm that uses a single surface, with no need for a confining wall or grafting. Our "wall-free" method is both more efficient than the traditional ones and free of additional assumptions or approximations. We illustrate it by simulating an HP-model lattice peptide on planar surfaces with a variety of patterns of adsorption sites, discovering a temperature-induced switch of surface selection which is due to a balance of energetic and entropic effects.

Phase diagrams of knotted and unknotted ring polymers

The phase diagram for a lattice ring polymer under applied force, with variable solvent quality, for different topological knot states, is determined for the first time. In addition to eliminating pseudophases where the polymer is flattened into a single layer, it is found that nontrivial knots result in additional pseudophases under tensile force conditions.

A new Monte Carlo method for getting the density of states of atomic cluster systems

A novel Monte Carlo flat histogram algorithm is proposed to get the classical density of states in terms of the potential energy, g(E(p)), for systems with continuous variables such as atomic clusters. It aims at avoiding the long iterative process of the Wang-Landau method and controlling carefully the convergence, but keeping the ability to overcome energy barriers. Our algorithm is based on a preliminary mapping in a series of points (called a σ-mapping), obtained by a two-parameter local probing of g(E(p)), and it converges in only two subsequent reweighting iterations on large intervals. The method is illustrated on the model system of a 432 atom cluster bound by a Rydberg type potential. Convergence properties are first examined in detail, particularly in the phase transition zone. We get g(E(p)) varying by a factor 10(3700) over the energy range [0.01 < E(p) < 6000 eV], covered by only eight overlapping intervals. Canonical quantities are derived, such as the internal energy U(T) and the heat capacity C(V)(T). This reveals the solid to liquid phase transition, lying in our conditions at the triple point. This phase transition is further studied by computing a Lindemann-Berry index, the atomic cluster density n(r), and the pressure, demonstrating the progressive surface melting at this triple point. Some limited results are also given for 1224 and 4044 atom clusters.

Monte Carlo simulations of the HP model (the "Ising model" of protein folding)

Using Wang-Landau sampling with suitable Monte Carlo trial moves (pull moves and bond-rebridging moves combined) we have determined the density of states and thermodynamic properties for a short sequence of the HP protein model. For free chains these proteins are known to first undergo a collapse "transition" to a globule state followed by a second "transition" into a native state. When placed in the proximity of an attractive surface, there is a competition between surface adsorption and folding that leads to an intriguing sequence of "transitions". These transitions depend upon the relative interaction strengths and are largely inaccessible to "standard" Monte Carlo methods.