Search

The transport of water and protons in the cathode catalyst layer (CCL) of proton exchange membrane (PEM) fuel cells is critical for cell performance, but the underlying mechanism is still unclear. Herein, the ionomer structure and the distribution/transport characteristics of water and protons in CCLs are investigated via all-atom molecular dynamics simulations. The results show that at low water contents, isolated water clusters form in ionomer pores, while proton transport is mainly via the charged sites of the ionomer side chains and the Grotthuss mechanism. Moreover, with increasing water content, water clusters are interconnected to form continuous water channels, which provide effective paths for proton transfer via the vehicular and Grotthuss mechanisms. Increasing the ionomer mass content can enhance the dense arrangement of the ionomer, which in turn increases the density of charge sites and improves the proton transport efficiency. When the ionomer mass content is high, the clustering effect reduces the space for water diffusion, increases the proton transport path, and finally decreases the proton transport efficiency. By providing physics insights into the proton transport mechanism, this study is helpful for the structural design and performance improvement of CCLs of PEM fuel cells., Comment: 14 pages, 15 figures

We investigate the coalescence of surfactant-laden water droplets by using several different surfactant types and a wide range of concentrations by means of a coarse-grained model obtained by the statistical associating fluid theory. Our results demonstrate in detail a universal mass transport mechanism of surfactant across many concentrations and several surfactant types during the process. Coalescence initiation is seen to occur via a single pinch due to aggregation of surface surfactant, and its remnants tend to become engulfed in part inside the forming bridge. Across the board we confirm the existence of an initial thermal regime with constant bridge width followed by a later inertial regime with bridge width scaling roughly as the square root of time, but see no evidence of an intermediate viscous regime. Coalescence becomes slower as surfactant concentration grows, and we see evidence of the appearance of a further slowdown of a different nature for several times the critical concentration. We anticipate that our results provide further insights in the mechanisms of coalescence of surfactant-laden droplets., Comment: 13 pages, 11 figures

Durotaxis motion is a spectacular phenomenon manifesting itself by the autonomous motion of a nano-object between parts of a substrate with different stiffness. This motion usually takes place along a stiffness gradient from softer to stiffer parts of the substrate. Here, we propose a new design of a polymer brush substrate that demonstrates antidurotaxis droplet motion, that is, droplet motion from stiffer to softer parts of the substrate. By carrying out extensive molecular dynamics simulation of a coarse-grained model, we find that antidurotaxis is solely controlled by the gradient in the grafting density of the brush and is favorable for fluids with a strong attraction to the substrate (low surface energy). The driving force of the antidurotaxial motion is the minimization of the droplet-substrate interfacial energy, which is attributed to the penetration of the droplet into the brush. Thus, we anticipate that the proposed substrate design offers a new understanding and possibilities in the area of autonomous motion of droplets for applications in microfluidics, energy conservation, and biology., Comment: 12 pages, 9 figures

MARTINI is a popular coarse-grained force-field that is mainly used in molecular dynamics (MD) simulations. It is based on the ``Lego'' approach where intermolecular interactions between coarse-grained beads representing chemical units of different polarity are obtained through water--octanol partition coefficients. This enables the simulation of a wide range of molecules by only using a finite number of parametrized coarse-grained beads, similar to the Lego game, where a finite number of bricks are used to create larger structures. Moreover, the MARTINI force-field is based on the Lennard-Jones potential with the shortest possible cutoff including attractions, thus rendering it very efficient for MD simulations. However, MD simulation is in general a computationally expensive method. Here, we demonstrate that using the MARTINI ``Lego'' approach is suitable for many-body dissipative particle (MDPD) dynamics, a method that can simulate multi-component and multi-phase soft matter systems in a much faster time (about 4--7 times) than MD. In this study, a DPPC lipid bilayer is chosen to provide evidence for the validity of this approach and various properties are compared to highlight the potential of the method. Thus, we anticipate that our study opens new possibilities for faster simulations of a wide range of soft matter systems by using the MDPD method., Comment: 9 pages, 5 figures

Droplet coalescence is ubiquitous in nature and the same time key to various technologies, such as inkjet printing. Here, we report on the coalescence of polymer droplets with different chain lengths coalescing on substrates of different wettability. By means of molecular dynamics simulations of a coarse-grained model, it is found that the rate of bridge growth is higher in the case of droplets with smaller contact angles (more wettable substrates) and decreases with the increase of the chain length of the polymers. Different behavior has also been identified in the dynamics of the approach of the two droplets during coalescence with the substrate wettability playing a more important role compared to the chain length of the polymers. While the dynamics of the droplet are greatly affected by the latter parameters, the density profile and flow patterns remain the same for the different cases. Thus, we anticipate that our work provides further insights into the coalescence of liquid polymer droplets on solid substrates with implications for relevant technologies., Comment: 10 pages, 5 figures

Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between state and flux coupling and discuss the various assumptions required for both. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature for both pressure tensor and coupling constraints. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines., Comment: 21 pages, 8 figures perspective paper submitted to Physical Chemistry Chemical Physics (PCCP)

The breakup of liquid threads into smaller droplets is a fundamental problem in fluid dynamics. In this study, we estimate the characteristic wavelength of the breakup process by means of many-body dissipative particle dynamics. This wavelength shows a power-law dependence on the Ohnesorge number in line with results from stability analysis. We also discover that the number of satellite droplets exhibits a power-law decay with exponent $0.72 \pm 0.04$ in the product of the Ohnesorge and thermal capillary numbers, while the overall size of main droplets is larger than that based on the characteristic wavelength thanks to the asynchronous breakup of the thread. Finally, we show that the formation of satellite droplets is the result of the advection of pinching points towards the main droplets in a remaining thinning neck, when the velocity gradient of the fluid exhibits two symmetric maxima., Comment: 18 pages, 9 figures

Nano-patterned substrates offer possibilities for controlling the motion of fluids without external energy supply in novel technologies in microfluidics, coatings, etc. Here, we report on the rugotaxial motion of droplets on wrinkled substrates with gradient in the wavelength of the wrinkles by exploring a broad range of parameters, such as amplitude of the wrinkles, substrate wettability, droplet size and wavelength gradient. Adopting a theoretical and molecular dynamics approach, we determine the Cassie-Baxter and Wenzel states of the droplets, investigate the efficiency of rugotaxis as a function of different parameters, and discuss additional effects, such as pinning. We find that shallow wrinkles characterised by small wavelength gradients, and moderate adhesion of the droplet to the substrate favour the rugotaxis motion with growing droplet size, when pinning is avoided. We also find that the driving force in rugotaxis is the gain in interfacial energy between the droplet and the substrate as the droplet enters regions of denser wrinkles (smaller wavelengths of the wrinkles)., Comment: 7 pages, 5 figures

Droplet coalescence is an important process in nature and various technologies (e.g. inkjet printing). Here, we unveil the surfactant mass-transport mechanism and report on several major differences in the coalescence of surfactant-laden droplets as compared to pure water droplets by means of molecular dynamics simulation of a coarse-grained model. Large scale changes to bridge growth dynamics are identified, such as the lack of multiple thermally excited precursors, attenuated collective excitations after contact, slowing down in the inertial regime due to aggregate-induced rigidity and reduced water flow, and a slowing down in the coalescence rate (deceleration) when surfactant concentration increases, while at the same time we also confirm the existence of an initial thermal, and a power-law, inertial, regime of the bridge growth dynamics in both the pure and the surfactant-laden droplets. Thus, we unveil the key mechanisms in one of the fundamental topological processes of liquid droplets containing surfactant, which is crucial in relevant technologies., Comment: 23 pages, 7 figures

Using extensive molecular dynamics simulation of a coarse-grained model, we demonstrate the possibility of sustained unidirectional motion (durotaxis) of droplets without external energy supply when placed on a polymer brush substrate with stiffness gradient in a certain direction. The governing key parameters for the specific substrate design studied, which determine the durotaxis efficiency, are found to be the grafting density of the brush and the droplet adhesion to the brush surface, whereas the strength of the stiffness gradient, the viscosity of the droplet, or the length of the polymer chains of the brush have only a minor effect on the process. It is shown that this durotaxial motion is driven by the steady increase of the interfacial energy between droplet and brush as the droplet moves from softer to stiffer parts of the substrate whereby the mean driving force gradually declines with decreasing roughness of the brush surface. We anticipate that our findings indicate further possibilities in the area of nanoscale motion without external energy supply.

As one of the most fascinating confined water/ice phenomena, two-dimensional square ice has been extensively studied and experimentally confirmed in recent years. Apart from the unidirectional homogeneous square icing patterns considered in previous studies, the multidirectional partitioned square icing patterns are discovered in this study and characterized by molecular dynamics (MD) simulations. Square icing parameters are proposed to quantitatively distinguish the partitioned patterns from the homogeneous patterns and the liquid water. The number of graphene monolayers n is varied in this study, and the results show that it is more energetically favorable to form partitioned square icing patterns when the water molecules are confined between graphite sheets (n >= 2) compared to graphene (n = 1). This phenomenon is insensitive to n as long as n >= 2, because of the short-range nature of the interaction between water molecules and the carbon substrate. Moreover, it is energetically unfavorable to form partitioned square icing patterns for a single layer of water molecules even for n >= 2, verifying that the interaction between layers of water molecules is another dominant factor in the formation of partitioned structures. The conversion from partitioned structure to homogenous square patterns is investigated by changing the pressure and the temperature. Based on the comprehensive MD simulations, this study unveils the formation mechanism of the partitioned square icing patterns., Comment: 10 pages, 9 figures

Atmospheric aerosols can consist of inorganic and organic substances, including surfactants at a significant concentration. Importantly, the latter can reduce the surface tension at the liquid-vapor surfaces, where they preferentially adsorb due to their amphiphilic structure. As a result, processes such as droplet coalescence, development of precipitation and ultimately cloud lifetime, may depend on the presence of surfactants in the aerosols. Here, we present a numerical model for cloud droplet formation, which is based on the Lagrangian particle-based microphysics-scheme super-droplet method and takes into account the presence of surfactant in the droplets. Our results show that surfactant facilitates cloud formation by increasing the number and size of activated droplets, which concentrate at the bottom of the cloud, while the largest droplets are concentrated at the top of the cloud. This indicates a circulation of droplets that involves activation and growth processes from the bottom of the cloud towards the top. Moreover, our conclusions are independent of the particular approach used for modeling the diffusion of Eulerian variables due to the subgrid-scale turbulence. We anticipate that our results will enrich our understanding of the role of surfactants in the behavior of atmospheric aerosols and, importantly, will pave the way for further developments in the numerical modeling of systems with surfactants at macroscopic scales., Comment: 25 pages, 9 figures

We investigate the scaling of the interfacial adsorption of the two-dimensional Blume-Capel model using Monte Carlo simulations. In particular, we study the finite-size scaling behavior of the interfacial adsorption of the pure model at both its first- and second-order transition regimes, as well as at the vicinity of the tricritical point. Our analysis benefits from the currently existing quite accurate estimates of the relevant (tri)critical-point locations. In all studied cases, the numerical results verify to a level of high accuracy the expected scenarios derived from analytic free-energy scaling arguments. We also investigate the size dependence of the interfacial adsorption under the presence of quenched bond randomness at the originally first-order transition regime (disorder-induced continuous transition) and the relevant self-averaging properties of the system. For this ex-first-order regime, where strong transient effects are shown to be present, our findings support the scenario of a non-divergent scaling, similar to that found in the original second-order transition regime of the pure model., Comment: 6 pages, 5 figures, version published in Phys. Rev. E. arXiv admin note: text overlap with arXiv:1610.08224

Using high-precision Monte-Carlo simulations based on a parallel version of the Wang-Landau algorithm and finite-size scaling techniques we study the effect of quenched disorder in the crystal-field coupling of the Blume-Capel model on the square lattice. We mainly focus on the part of the phase diagram where the pure model undergoes a continuous transition, known to fall into the universality class of the pure Ising ferromagnet. A dedicated scaling analysis reveals concrete evidence in favor of the strong universality hypothesis with the presence of additional logarithmic corrections in the scaling of the specific heat. Our results are in agreement with an early real-space renormalization-group study of the model as well as a very recent numerical work where quenched randomness was introduced in the energy exchange coupling. Finally, by properly fine tuning the control parameters of the randomness distribution we also qualitatively investigate the part of the phase diagram where the pure model undergoes a first-order phase transition. For this region, preliminary evidence indicate a smoothening of the transition to second-order with the presence of strong scaling corrections., Comment: 11 pages, 13 figures, minor correction to references appearing in Fig. 1, to be published in Phys. Rev. E

Surface nanobubbles (NBs) are stable gaseous phases in liquids that form at the interface with solid substrates. They have been particularly intriguing for their high stability that contradicts theoretical expectations and their potential relevance for many technological applications. Here, we present the current state of the art in this research area by discussing and contrasting main results obtained from theory, simulation and experiment, and presenting their limitations. We also provide future perspectives anticipating that this review will stimulate further studies in the research area of surface NBs., Comment: 47 pages, 6 figures

The surfactant-driven spreading of droplets is an essential process in many applications ranging from coating flow technology to enhanced oil recovery. Despite the significant advancement in describing spreading processes in surfactant-laden droplets, including the exciting phenomena of superspreading, many features of the underlying mechanisms require further understanding. Here, we have carried out molecular dynamics simulations of a coarse-grained model with a force-field obtained from the statistical associating fluid theory to study droplets laden with common and superspreading surfactants. We have confirmed the important elements of the superspreading mechanism, i.e. the adsorption of surfactant at the contact line and the fast replenishment of surfactant from the bulk. Through a detailed comparison of a range of droplets with different surfactants, our analysis has indicated that the ability of surfactant to adsorb at the interfaces is the key feature of the superspreading mechanism. To this end, surfactants that tend to form aggregates and have a strong hydrophobic attraction in the aggregated cores prevent the fast replenishment of the interfaces, resulting in reduced spreading ability. We also show that various surfactant properties, such as diffusion and architecture, play a secondary role in the spreading process. Moreover, we discuss various drop properties such as the height, contact angle, and surfactant surface concentration, highlighting differences between superspreading and common surfactants. We anticipate that our study will provide further insight for applications requiring the control of wetting., Comment: 26 pages, 11 figures

Polymer-grafted nanoparticles (PGNPs) can provide property profiles than cannot be obtained individually by polymers or nanoparticles (NPs). Here, we have studied the mixing--demixing transition of symmetric copolymer melts of polymer-grafted spherical nanoparticles by means of coarse-grained molecular dynamics simulation and a theoretical mean-field model. We find that a larger size of NPs leads to higher stability for given number of grafted chains and chain length reaching a point where demixing is not possible. Most importantly, the increase in the number of grafted chains, $N_g$, can initially favour the phase separation of PGNPs, but further increase can lead to more difficult demixing. The reason is the increasing impact of an effective core that forms as the grafting density of the tethered polymer chains around the NPs increases. The range and exact values of $N_g$ where this change in behaviour takes place depends on the NP size and the chain length of the grafted polymer chains. Our study elucidates the phase behaviour of PGNPs and in particular the influence of the grafting density on the phase behaviour of the systems anticipating that it will open new doors in the understanding of these systems with implications in materials science and medicine., Comment: 6 pages, 4 figures, final version to be published in Soft Matter

We study using Monte Carlo simulations the finite-size scaling behavior of the interfacial adsorption of the two-dimensional square-lattice $q$-states Potts model. We consider the pure and random-bond versions of the Potts model for $q = 3,4,5,8$ and $q = 10$, thus probing the interfacial properties at the originally continuous, weak, and strong first-order phase transitions. For the pure systems our results support the early scaling predictions for the size dependence of the interfacial adsorption at both first- and second-order phase transitions. For the disordered systems, the interfacial adsorption at the (disordered induced) continuous transitions is discussed, applying standard scaling arguments and invoking findings for bulk critical properties. The self-averaging properties of the interfacial adsorption are also analyzed by studying the infinite limit-size extrapolation of properly defined signal-to-noise ratios., Comment: 7 pages, 5 figures, final version to be published in Phys. Rev. E. arXiv admin note: text overlap with arXiv:1504.07425

Surface nanobubbles are stable gaseous phases in liquids that form on solid substrates. While their existence has been confirmed, there are many open questions related to their formation and dissolution processes along with their structures and properties, which are difficult to investigate experimentally. To address these issues, we carried out molecular dynamics simulations based on atomistic force fields for systems comprised of water, air (N2 and O2), and a Highly Oriented Pyrolytic Graphite (HOPG) substrate. Our results provide insights into the formation/dissolution mechanisms of nanobubbles and estimates for their density, contact angle, and surface tension. We found that the formation of nanobubbles is driven by an initial nucleation process of air molecules and the subsequent coalescence of the formed air clusters. The clusters form favorably on the substrate, which provides an enhanced stability to the clusters. In contrast, nanobubbles formed in the bulk either move randomly to the substrate and spread or move to the water--air surface and pop immediately. Moreover, nanobubbles consist of a condensed gaseous phase with a surface tension smaller than that of an equivalent system under atmospheric conditions, and contact angles larger than those in the equivalent nanodroplet case. We anticipate that this study will provide useful insights into the physics of nanobubbles and will stimulate further research in the field by using all-atom simulations.

Colloidal particles grafted with single-stranded DNA (ssDNA) chains can self-assemble into a number of different crystalline structures, where hybridization of the ssDNA chains creates links between colloids stabilizing their structure. Depending on the geometry and the size of the particles, the grafting density of the ssDNA chains, and the length and choice of DNA sequences, a number of different crystalline structures can be fabricated. However, understanding how these factors contribute synergistically to the self-assembly process of DNA-functionalized nano- or micro-sized particles remains an intensive field of research. Moreover, the fabrication of long-range structures due to kinetic bottlenecks in the self-assembly are additional challenges. Here, we discuss the most recent advances from theory and experiment with particular focus put on recent simulation studies., Comment: 24 pages, 2 figures, Review Article

The phase behaviour of two-component bottle-brush polymers with fully flexible backbones under poor solvent conditions is studied via molecular-dynamics simulations, using a coarse-grained bead-spring model and side chains of up to $N=40$ effective monomers. We consider a symmetric model where side chains of type A and B are grafted alternately onto a flexible backbone. The aim of this study to explore the phase behaviour of two-component bottle-brushes depending on parameters, such as as the grafting density $\sigma$, the backbone length $N_b$, the side-chain length $N$, and the temperature $T$. Based on a cluster analysis, we identify for our range of parameters the regimes of fully phase separated systems, i.e., A-type side chains form one cluster and B-type chains another, while the interface that separates these two clusters contains the backbone monomers. We find that pearl-necklace or Janus-like structures, which normally occur for bottle-brush polymers with rigid backbones under poor solvent conditions, are fully attributed to the backbone rigidity, and, therefore, such structures are unlikely in the case of bottle brushes with fully flexible backbones. Also, a comparative discussion with earlier work on the phase behaviour of single-component bottle-brush polymers with flexible backbones is performed., Comment: 16 pages, 5 figures

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 $\sigma$ with which side chains are grafted onto the flexible backbone, we study for backbone lengths of up to $N_b=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 chains 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 ($\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 $\sigma$. 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., Comment: 18 pages, 11 eps figures, to be published in J. Phys.: Condens. Matter

We use a standard bead-spring model and molecular dynamics simulations to study the static properties of symmetric linear multiblock copolymer chains and their blocks under poor solvent conditions in a dilute solution from the regime close to theta conditions, where the chains adopt a coil-like formation, to the poorer solvent regime where the chains collapse obtaining a globular formation and phase separation between the blocks occurs. We choose interaction parameters as is done for a standard model, i.e., the Lennard-Jones fluid and we consider symmetric chains, i.e., the multiblock copolymer consists of an even number $n$ of alternating chemically different A and B blocks of the same length $N_{A}=N_{B}=N$. We show how usual static properties of the individual blocks and the whole multiblock chain can reflect the phase behavior of such macromolecules. Also, how parameters, such as the number of blocks $n$ can affect properties of the individual blocks, when chains are in a poor solvent for a certain range of $n$. A detailed discussion of the static properties of these symmetric multiblock copolymers is also given. Our results in combination with recent simulation results on the behavior of multiblock copolymer chains provide a complete picture for the behavior of these macromolecules under poor solvent conditions, at least for this most symmetrical case. Due to the standard choice of our parameters, our system can be used as a benchmark for related models, which aim at capturing the basic aspects of the behavior of various biological systems., Comment: 13 pages, 11 figures

By Molecular Dynamics simulation of a coarse-grained bead-spring type model for a cylindrical molecular brush with a backbone chain of $N_b$ effective monomers to which with grafting density $\sigma$ side chains with $N$ effective monomers are tethered, several characteristic length scales are studied for variable solvent quality. Side chain lengths are in the range $5 \le N \le 40$, backbone chain lengths are in the range $50 \le N_b \le 200$, and we perform a comparison to results for the bond fluctuation model on the simple cubic lattice (for which much longer chains are accessible, $N_b \le 1027$, and which corresponds to an athermal, very good, solvent). We obtain linear dimensions of side chains and the backbone chain and discuss their $N$-dependence in terms of power laws and the associated effective exponents. We show that even at the Theta point the side chains are considerably stretched, their linear dimension depending on the solvent quality only weakly. Effective persistence lengths are extracted both from the orientational correlations and from the backbone end-to-end distance; it is shown that different measures of the persistence length (which would all agree for Gaussian chains) are not mutually consistent with each other, and depend distinctly both on $N_b$ and the solvent quality. A brief discussion of pertinent experiments is given., Comment: 30 pages, 13 figures, 1 table

We report a simulation study for bottle-brush polymers grafted on a rigid backbone. Using a standard coarse-grained bead-spring model extensive molecular dynamics simulations for such macromolecules under good solvent conditions are performed. We consider a broad range of parameters and present numerical results for the monomer density profile, density of the untethered ends of the grafted flexible backbones and the correlation function describing the range that neighboring grafted bottle-brushes are affected by the presence of the others due to the excluded volume interactions. The end beads of the flexible backbones of the grafted bottle-brushes do not access the region close to the rigid backbone due to the presence of the side chains of the grafted bottle-brush polymers, which stretch further the chains in the radial directions. Although a number of different correlation lengths exist as a result of the complex structure of these macromolecules, their properties can be tuned with high accuracy in good solvents. Moreover, qualitative differences with "typical" bottle-brushes are discussed. Our results provide a first approach to characterizing such complex macromolecules with a standard bead spring model., Comment: To appear in Journal of Physics Condensed Matter (2011)

When two types of side chains (A,B) are densely grafted to a (stiff) backbone and the resulting bottle-brush polymer is in a solution under poor solvent conditions, an incompatibility between A and B leads to microphase separation in the resulting cylindrical brush. The possible types of ordering are reminiscent of the ordering of block copolymers in cylindrical confinement. Starting from this analogy, Leibler's theory of microphase separation in block copolymer melts is generalized to derive a description of the system in the weak segregation limit. Also molecular dynamics simulation results of a corresponding coarse-grained bead-spring model are presented. Using side chain lengths up to N = 50 effective monomers, the ratio of the Lennard-Jones energy parameter between unlike monomers $(\epsilon_{AB})$ and monomers of the same kind $(\epsilon _{AA} = \epsilon_{BB})$ is varied. Various correlation functions are analyzed to study the conditions when (local) Janus cylinder-type ordering and when (local) microphase separation in the direction along the cylinder axis occurs. Both the analytical theory and the simulations give evidence for short range order due to a tendency towards microphase separation in the axial direction, with a wavelength proportional to the side chain gyration radius, irrespective of temperature and grafting density, for a wide range of these parameters., Comment: 26 pages, 19 figures

We study the critical behavior of the $d=3$ Ising model with bond randomness through extensive Monte Carlo simulations and finite-size scaling techniques. Our results indicate that the critical behavior of the random-bond model is governed by the same universality class as the site- and bond-diluted models, clearly distinct from that of the pure model, thus providing a complete set of universality in disordered systems., Comment: 5 pages, 4 figures, to be published in Phys. Rev. E as a brief report

Molecular dynamics simulations are used to study the phase behavior of linear multiblock copolymers with two types of monomers, A and B, where the length of the polymer blocks $N_{A}$ and $N_{B}$ ($N_{A}=N_{B}=N$), the number of the blocks $n_{A}$ and $n_{B}$ ($n_{A}=n_{B}=n$), and the solvent quality varies. The fraction $f$ of A-type monomers is kept constant and equal to 0.5. Whereas at high enough temperatures these macromolecules form coil structures, where each block A or B forms rather individual clusters, at low enough temperatures A and B monomers from different blocks can join together forming clusters of A or B monomers. The dependence of the formation of these clusters on the varied parameters is discussed in detail, providing a full understanding of the phase behavior of linear multiblock copolymers, at least for this symmetrical case., Comment: 8 pages, 12 figures, to be published in Soft Matter

As an example for the interplay of structure, dynamics, and phase behavior of macromolecular systems, this article focuses on the problem of bottle-brush polymers with either rigid or flexible backbones. On a polymer with chain length $N_b$, side-chains with chain length $N$ are endgrafted with grafting density $\sigma$. Due to the multitude of characteristic length scales and the size of these polymers (typically these cylindrical macromolecules contain of the order of 10000 effective monomeric units) understanding of the structure is a challenge for experiment. But due to excessively large relaxation times (particularly under poor solvent conditions) such macromolecules also are a challenge for simulation studies. Simulation strategies to deal with this challenge, both using Monte Carlo and Molecular Dynamics Methods, will be briefly discussed, and typical results will be used to illustrate the insight that can be gained., Comment: 8 pages, 7 figures, to be published in the proceedings of "NIC Symposium 2010", eds. G. Muenster et al. (NIC, Juelich, 2010)