Multidimensional stochastic approximation Monte Carlo

Crystallisation in Melts of Short, Semi-Flexible Hard-Sphere Polymer Chains: The Role of the Non-Bonded Interaction Range

A melt of short semi-flexible polymers with hard-sphere-type non-bonded interaction undergoes a first-order crystallisation transition at lower density than a melt of hard-sphere monomers or a flexible hard-sphere chain. In contrast to the flexible hard-sphere chains, the semi-flexible ones have an intrinsic stiffness energy scale, which determines the natural temperature scale of the system. In this paper, we investigate the effect of weak additional non-bonded interaction on the phase transition temperature. We study the system using the stochastic approximation Monte Carlo (SAMC) method to estimate the micro-canonical entropy of the system. Since the density of states in the purely hard-sphere non-bonded interaction case already covers 5600 orders of magnitude, we consider the effect of weak interactions as a perturbation. In this case, the system undergoes the same ordering transition with a temperature shift non-uniformly depending on the additional interaction. Short-range attractions impede ordering of the melt of semi-flexible polymers and decrease the transition temperature, whereas relatively long-range attractions assist ordering and shift the transition temperature to higher values, whereas weak repulsive interactions demonstrate an opposite effect on the transition temperature.

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.

Control of accuracy in the Wang-Landau algorithm

The Wang-Landau (WL) algorithm has been widely used for simulations in many areas of physics. Our analysis of the WL algorithm explains its properties and shows that the difference of the largest eigenvalue of the transition matrix in the energy space from unity can be used to control the accuracy of estimating the density of states. Analytic expressions for the matrix elements are given in the case of the one-dimensional Ising model. The proposed method is further confirmed by numerical results for the one-dimensional and two-dimensional Ising models and also the two-dimensional Potts model.

Macroscopically constrained Wang-Landau method for systems with multiple order parameters and its application to drawing complex phase diagrams

A generalized approach to Wang-Landau simulations, macroscopically constrained Wang-Landau, is proposed to simulate the density of states of a system with multiple macroscopic order parameters. The method breaks a multidimensional random-walk process in phase space into many separate, one-dimensional random-walk processes in well-defined subspaces. Each of these random walks is constrained to a different set of values of the macroscopic order parameters. When the multivariable density of states is obtained for one set of values of fieldlike model parameters, the density of states for any other values of these parameters can be obtained by a simple transformation of the total system energy. All thermodynamic quantities of the system can then be rapidly calculated at any point in the phase diagram. We demonstrate how to use the multivariable density of states to draw the phase diagram, as well as order-parameter probability distributions at specific phase points, for a model spin-crossover material: an antiferromagnetic Ising model with ferromagnetic long-range interactions. The fieldlike parameters in this model are an effective magnetic field and the strength of the long-range interaction.