The non-equilibrium dynamics of condensation phenomena in nano-pores is studied via Monte Carlo simulation of a lattice gas model. Hysteretic behavior of the particle density as a function of the density of a reservoir is obtained for various pore geometries in two and three dimensions. The shape of the hysteresis loops depend on the characteristics of the pore geometry. The evaporation of particles from a pore can be fitted to a stretched exponential decay of the particle density $rho_f(t) sim exp [ -(t/tau)^beta]$. Phase separation dynamics inside the pore is effectively described by a random walk of the non-wetting phases. Domain evolution is significantly slowed down in presence of random wall-particle potential and gives rise to a temperature dependent growth exponent. On the other hand roughness of the pore wall only delays the onset of a pure domain growth.
The amount and type of self-entanglement of DNA filaments is significantly affected by spatial confinement, which is ubiquitous in biological systems. Motivated by recent advancements in single DNA molecule experiments based on nanofluidic devices, and by the introduction of algorithms capable of detecting knots in open chains, we investigate numerically the entanglement of linear, open DNA chains confined inside nano-slits. The results regard the abundance, type and length of occurring knots and are compared with recent findings for DNA inside nano-channels. In both cases, the width of the confining region, D, spans the 30nm- 1mu m range and the confined DNA chains are 1 to 4mu m long. It is found that the knotting probability is maximum for slit widths in the 70-100nm range. However, over the considered DNA contour lengths, the maximum incidence of knots remains below 20%, while for channel confinement it tops 50%. Further differences of the entanglement are seen for the average contour length of the knotted region which drops significantly below D ~100nm for channel-confinement, while it stays approximately constant for slit-like confinement. These properties ought to reverberate in different kinetic properties of linear DNA depending on confinement and could be detectable experimentally or exploitable in nano-technological applications.
We study non-equilibrium phases for interacting two-dimensional self-propelled particles with isotropic pair-wise interactions using a persistent kinetic Monte Carlo (MC) approach. We establish the quantitative phase diagram, including the motility-induced phase separation (MIPS) that is a commonly observed collective phenomena in active matter. In addition, we demonstrate for several different potential forms the presence of two-step melting, with an intermediate hexatic phase, in regions far from equilibrium. Increased activity can melt a two-dimensional solid and the melting lines remain disjoint from MIPS. We establish this phase diagram for a range of the inter-particle potential stiffnesses, and identify the MIPS phase even in the hard-disk limit. We establish that the full description of the phase behavior requires three independent control parameters.
A new Monte Carlo approach is proposed to investigate the fluid-solid phase transition of the polydisperse system. By using the extended ensemble, a reversible path was constructed to link the monodisperse and corresponding polydisperse system. Once the fluid-solid coexistence point of the monodisperse system is known, the fluid-solid coexistence point of the polydisperse system can be obtained from the simulation. The validity of the method is checked by the simulation of the fluid-solid phase transition of a size-polydisperse hard sphere colloid. The results are in agreement with the previous studies.
Extensive Monte Carlo results are presented for a lattice model of a bottle-brush polymer under good solvent or Theta solvent conditions. Varying the side chain length, backbone length, and the grafting density for a rigid straight backbone, both radial density profiles of monomers and side chain ends are obtained, as well as structure factors describing the scattering from a single side chain and from the total bottle-brush polymer. To describe the structure in the interior of a very long bottle-brush, a periodic boundary condition in the direction along the backbone is used, and to describe effects due to the finiteness of the backbone length, a second set of simulations with free ends of the backbone is performed. In the latter case, the inhomogeneity of the structure in the direction along the backbone is carefully investigated. We use these results to test various phenomenological models that have been proposed to interpret experimental scattering data for bottle-brush macromolecules. These models aim to extract information on the radial density profile of a bottle-brush from the total scattering via suitable convolution approximations. Possibilities to improve such models, guided by our simulation results, are discussed.
A Grand-canonical Monte-Carlo simulation method extended to simulate a mixture of salts is presented. Due to charge neutrality requirement of electrolyte solutions, ions must be added to or removed from the system in groups. This leads to some complications compared to regular Grand Canonical simulation. Here, a recipe for simulation of electrolyte solution of salt mixture is presented. It is then implemented to simulate solution of 1:1, 2:1 and 2:2 salts or their mixtures at different concentrations using the primitive ion model. The osmotic pressures of the electrolyte solutions are calculated and shown to depend linearly on the salt concentrations within the concentration range simulated. We also show that at the same concentration of divalent anions, the presence of divalent cations make it easier to insert monovalent cations into the system. This can explain some quantitative differences observed in experiments of the MgCl$_2$ salt mixture and MgSO$_4$ salt mixture.