Analysis of the static properties of cluster formations in symmetric linear multiblock copolymers

Thermal properties of knotted block copolymer rings with charged monomers subjected to short-range interactions

The thermal properties of coarse-grained knotted copolymer rings fluctuating in a highly screening solution are investigated on a simple cubic lattice using the Wang-Landau Monte Carlo algorithm. The rings contain two kinds of monomers A and B with opposite charges that are subjected to short-range interactions. In view of possible applications in medicine and the construction of intelligent materials, it is shown that the behavior of copolymer rings can be tuned by changing both their monomer configuration and topology. We find several phase transitions depending on the monomer distribution. They include the expansion and collapse of the knotted polymer as well as rearrangements leading to metastable states. The temperatures at which these phase transitions are occurring and other features can be tuned by changing the topology of the system. The processes underlying the observed transitions are identified. In knots formed by diblock copolymers, two different classes of behaviors are detected depending on whether there is an excess of monomers of one kind or not. Moreover, we find that the most stable compact states are formed by copolymers in which units of two A monomers are alternated by units of two B monomers. Remarkably, these compact states are in a lamellar phase. The transition from the lamellar to the expanded state produces in the specific heat capacity a narrow and high peak that is centered at temperatures that are much higher than those of the peaks observed in all other monomer distributions.

Diagrams of States of Single Flexible-Semiflexible Multi-Block Copolymer Chains: A Flat-Histogram Monte Carlo Study

The combination of flexibility and semiflexibility in a single molecule is a powerful design principle both in nature and in materials science. We present results on the conformational behavior of a single multiblock-copolymer chain, consisting of equal amounts of Flexible (F) and Semiflexible (S) blocks with different affinity to an implicit solvent. We consider a manifold of macrostates defined by two terms in the total energy: intermonomer interaction energy and stiffness energy. To obtain diagrams of states (pseudo-phase diagrams), we performed flat-histogram Monte Carlo simulations using the Stochastic Approximation Monte Carlo algorithm (SAMC). We have accumulated two-Dimensional Density of States (2D DoS) functions (defined on the 2D manifold of macrostates) for a SF-multiblock-copolymer chain of length N=64 with block lengths b = 4, 8, 16, and 32 in two different selective solvents. In an analysis of the canonical ensemble, we calculated the heat capacity and determined its maxima and the most probable morphologies in different regions of the state diagrams. These are rich in various, non-trivial morphologies, which are formed without any specific interactions, and depend on the block length and the type of solvent selectivity (preferring S or F blocks, respectively). We compared the diagrams with those for the non-selective solvent and reveal essential changes in some cases. Additionally, we implemented microcanonical analysis in the “conformational” microcanonical (NVU, where U is the potential energy) and the true microcanonical (NVE, where E is the total energy) ensembles with the aim to reveal and classify pseudo-phase transitions, occurring under the change of temperature.

Efficient estimation of contact probabilities from inter-bead distance distributions in simulated polymer chains

The estimation of contact probabilities (CP) from conformations of simulated bead-chain polymer models is a key step in methods that aim to elucidate the spatial organization of chromatin from analysis of experimentally determined contacts between different genomic loci. Although CPs can be estimated simply by counting contacts between beads in a sample of simulated chain conformations, reliable estimation of small CPs through this approach requires a large number of conformations, which can be computationally expensive to obtain. Here we describe an alternative computational method for estimating relatively small CPs without requiring large samples of chain conformations. In particular, we estimate the CPs from functional approximations to the cumulative distribution function (cdf) of the inter-bead distance for each pair of beads. These cdf approximations are obtained by fitting the extended generalized lambda distribution (EGLD) to inter-bead distances determined from a sample of chain conformations, which are in turn generated by Monte Carlo simulations. We find that CPs estimated from fitted EGLD cdfs are significantly more accurate than CPs estimated using contact counts from samples of limited size, and are more precise with all sample sizes, permitting as much as a tenfold reduction in conformation sample size for chains of 200 beads and samples smaller than 10(5) conformations. This method of CP estimation thus has potential to accelerate computational efforts to elucidate the spatial organization of chromatin.

Molecular dynamics simulations of single-component bottle-brush polymers with flexible backbones under poor solvent conditions

Conformations of a single-component bottle-brush polymer with a fully flexible backbone under poor solvent conditions are studied by molecular dynamics simulations, using a coarse-grained bead-spring model with side chains of up to N = 40 effective monomers. By variation of the solvent quality and the grafting density σ with which side chains are grafted onto the flexible backbone, we study for backbone lengths of up to Nb = 100 the crossover from the brush/coil regime to the dense collapsed state. At lower temperatures, where collapsed chains with a constant monomer density are observed, the choice of the above parameters does not play any role and it is the total number of monomers that defines the dimensions of the chains. Furthermore, bottle-brush polymers with longer side chains possess higher spherical symmetry compared to polymers with lower side-chain lengths in contrast to what one may intuitively expect, as the stretching of the side chains is less important than the increase of their length. At higher temperatures, always below the Theta (Θ) temperature, coil-like configurations, similar to a single polymer chain, or brush-like configurations, similar to a homogeneous cylindrical bottle-brush polymer with a rigid backbone, are observed, depending on the choice of the particular parameters N and σ. In the crossover regime between the collapsed state (globule) and the coil/brush regime the acylindricity increases, whereas for temperatures outside of this range, bottle-brush polymers maintain a highly cylindrical symmetry in all configurational states.

A coarse-grained model for DNA-functionalized spherical colloids, revisited: effective pair potential from parallel replica simulations

We discuss a coarse-grained model recently proposed by Starr and Sciortino [J. Phys.: Condens. Matter 18, L347 (2006)] for spherical particles functionalized with short single DNA strands. The model incorporates two key aspects of DNA hybridization, i.e., the specificity of binding between DNA bases and the strong directionality of hydrogen bonds. Here, we calculate the effective potential between two DNA-functionalized particles of equal size using a parallel replica protocol. We find that the transition from bonded to unbonded configurations takes place at considerably lower temperatures compared to those that were originally predicted using standard simulations in the canonical ensemble. We put particular focus on DNA-decorations of tetrahedral and octahedral symmetry, as they are promising candidates for the self-assembly into a single-component diamond structure. Increasing colloid size hinders hybridization of the DNA strands, in agreement with experimental findings.