Search

In real crystals and at finite temperatures point defects are inevitable. Under shear their dynamics severely influence the mechanical properties of these crystals, giving rise to non-linear effects, such as ductility. In an effort to elucidate the complex behavior of crystals under plastic deformation it is crucial to explore and to understand the interplay between the timescale related to the equilibrium point-defect diffusion and the shear-induced timescale. Based on extensive non-equilibrium molecular dynamics simulations we present a detailed investigation on the yielding behavior of cluster crystals, an archetypical model for a defect-rich crystal: in such a system clusters of overlapping particles occupy the lattice sites of a regular (FCC) structure. In equilibrium particles diffuse via site-to-site hopping while maintaining the crystalline structure intact. We investigate these cluster crystals at a fixed density and at different temperatures where the system remains in the FCC structure: temperature allows us to vary the diffusion timescale appropriately. We then expose the crystal to shear, thereby choosing shear rates which cover timescales that are both higher and lower than the equilibrium diffusion timescales. We investigate the macroscopic and microscopic response of our cluster crystal to shear and find that the yielding scenario of such a system does not rely on the diffusion of the particles – it is rather related to the plastic deformation of the underlying crystalline structure. The local bond order parameters and the measurement of local angles between neighboring clusters confirm the cooperative movement of the clusters close to the yield point. Performing complementary, related simulations for an FCC crystal formed by harshly repulsive particles reveals similarities in the yielding behavior between both systems. Still we find that the diffusion of particles does influence characteristic features in the cluster crystal, such as a less prominent increase of order parameters close to the yield point. Our simulations provide for the first time an insight into the role of the diffusion of defects in the yielding behavior of a defect-rich crystal under shear. These observations will thus be helpful in the development of theories for the plastic deformation of defect-rich crystals., The yielding of cluster crystals under shear is primarily governed by the deformation of the underlying (FCC) crystalline structure.

We propose a computationally lean, two-stage approach that reliably predicts self-assembly behavior of complex charged molecules on a metallic surfaces under electrochemical conditions. Stage one uses ab initio simulations to provide reference data for the energies (evaluated for archetypical configurations) to fit the parameters of a conceptually much simpler and computationally less expensive model of the molecules: classical, spherical particles, representing the respective atomic entities, a soft but perfectly conductive wall potential represents the metallic surface. Stage two feeds the energies that emerge from this classical model into highly efficient and reliable optimization techniques to identify via energy minimization the ordered ground state configurations of the molecules. We demonstrate the power of our approach by successfully reproducing, on a semi-quantitative level, the intricate supramolecular ordering observed experimentally for PQP$^+$ and ClO$_4^-$ molecules at an Au(111)-electrolyte interface, including the formation of open-porous, self-hosts--guest, and stratified bilayer phases as a function of the electric field at the solid--liquid interface. We also discuss the role of the perchlorate ions in the self-assembly process, whose positions could not be identified in the related experimental investigations.

Various microorganisms and some mammalian cells are able to swim in viscous fluids by performing nonreciprocal body deformations, such as rotating attached flagella or by distorting their entire body. In order to perform chemotaxis, i.e. to move towards and to stay at high concentrations of nutrients, they adapt their swimming gaits in a nontrivial manner. We propose a model how microswimmers are able to autonomously adapt their shape in order to swim in one dimension towards high field concentrations using an internal decision making machinery modeled by an artificial neural network. We present two methods to measure chemical gradients, spatial and temporal sensing, as known for swimming mammalian cells and bacteria, respectively. Using the NEAT genetic algorithm surprisingly simple neural networks evolve which control the shape deformations of the microswimmer and allow them to navigate in static and complex time-dependent chemical environments. By introducing noisy signal transmission in the neural network the well-known biased run-and-tumble motion emerges. Our work demonstrates that the evolution of a simple internal decision-making machinery, which we can fully interpret and is coupled to the environment, allows navigation in diverse chemical landscapes. These findings are of relevance for intracellular biochemical sensing mechanisms of single cells, or for the simple nervous system of small multicellular organisms such as C. elegans., Comment: 9 pages, 4 figures

Direct correlation functions (DCFs), linked to the second functional derivative of the free energy with respect to the one-particle density, play a fundamental role in a statistical mechanics description of matter. This holds, in particular, for the ordered phases: DCFs contain information about the local structure including defects and encode the thermodynamic properties of crystalline solids; they open a route to the elastic constants beyond low temperature expansions. Via a demanding numerical approach, we have explicitly calculated for the first time the DCF of a solid: based on the fundamental measure concept, we provide results for the DCF of a hard sphere crystal. We demonstrate that this function differs at coexistence significantly from its liquid counterpart—both in shape as well as in its order of magnitude—because it is dominated by vacancies. We provide evidence that the traditional use of liquid DCFs in functional Taylor expansions of the free energy is conceptually wrong and show that the emergent elastic constants are in good agreement with simulation-based results. published

We employ genetic algorithms (GA), which allow for an unbiased search for the global minimum of energy landscapes, to identify the ordered equilibrium configurations formed by binary dipolar systems confined on a plane. A large variety of arrangements is identified, the complexity of which grows with increasing asymmetry between the two components and with growing concentration of the small particles. The effects of the density are briefly discussed and a comparison with results obtained via conventional lattice-sum minimization is presented. Our results can be confirmed by experiments involving Langmuir monolayers of polystyrene dipolar spheres or superparamagnetic colloids confined on the air–water interface and polarized by an external, perpendicular magnetic field.

We give evidence that particles interacting via the simple, radially symmetric square-shoulder potential can self-organise in highly complex, low-symmetry lattices, thereby forming clusters, columns, or lamellae; only at high pressure are compact, high-symmetry structures observed. Our search for these ordered equilibrium structures is based on ideas of genetic algorithms, a strategy that is characterised by a high success rate. A simple mean-field type consideration complements these findings and locates in a semi-quantitative way the cross-over between the competing structures.

The emergent fluctuating hydrodynamics of a viscoelastic fluid modeled by the multiparticle collision dynamics (MPC) approach is studied. The fluid is composed of flexible, Gaussian phantom polymers, which interact by local momentum-conserving stochastic MPC collisions. For comparison, the analytical solution of the linearized Navier-Stokes equation is calculated, where viscoelasticity is taken into account by a time-dependent shear relaxation modulus. The fluid properties are characterized by the transverse velocity autocorrelation function in Fourier space as well as in real space. Various polymer lengths are considered---from dumbbells to (near-)continuous polymers. Viscoelasticity affects the fluid properties and leads to strong correlations, which overall decay exponentially in Fourier space. In real space, the center-of-mass velocity autocorrelation function of individual polymers exhibits a long-time tail independent of polymer length, which decays as $t^{-3/2}$, similar to a Newtonian fluid, in the asymptotic limit $t \to \infty$. Moreover, for long polymers an additional power-law decay appears at time scales shorter than the longest polymer relaxation time with the same time dependence, but negative correlations, and the polymer length dependence $L^{-1/2}$. Good agreement is found between the analytical and simulation results.

We have applied the formalism of classical density functional theory to study the shape and the orientation of the density profiles formed by aspherical, ultrasoft particles. For simplicity we have...

With suitably designed Monte Carlo simulations, we have investigated the properties of mobile, impenetrable, yet deformable particles that are immersed into a porous matrix, the latter one realized by a frozen configuration of spherical particles. By virtue of a model put forward by Batista and Miller [Phys. Rev. Lett. 105, 088305 (2010)], the fluid particles can change in their surroundings, formed by other fluid particles or the matrix particles, their shape within the class of ellipsoids of revolution; such a change in shape is related to a change in energy, which is fed into suitably defined selection rules in the deformation "moves" of the Monte Carlo simulations. This concept represents a simple yet powerful model of realistic, deformable molecules with complex internal structures (such as dendrimers or polymers). For the evaluation of the properties of the system, we have used the well-known quenched-annealed protocol (with its characteristic double average prescription) and have analyzed the simulation data in terms of static properties (the radial distribution function and aspect ratio distribution of the ellipsoids) and dynamic features (notably the mean squared displacement). Our data provide evidence that the degree of deformability of the fluid particles has a distinct impact on the aforementioned properties of the system.

Inverse patchy colloids are nano- to micro-scale particles with a surface divided into differently charged regions. This class of colloids combines directional, selective bonding with a relatively simple particle design: owing to the competitive interplay between the orientation-dependent attraction and repulsion — induced by the interactions between like/oppositely charged areas — experimentally accessible surface patterns are complex enough to favor the stabilization of specific structures of interest. Most important, the behavior of heterogeneously charged units can be ideally controlled by means of external parameters, such as the pH and the salt concentration. We present a concise review about this class of systems, spanning the range from the synthesis of model inverse patchy particles to their self-assembly, covering their coarse-grained modeling and the related numerical/analytical treatments.

Star polymers with magnetically functionalized end groups are presented as a novel polymeric system whose morphology, self-aggregation, and orientation can easily be tuned by exposing these macromolecules simultaneously to an external magnetic field and to shear forces. Our investigations are based on a specialized simulation technique which faithfully takes into account the hydrodynamic interactions of the surrounding, Newtonian solvent. We find that the combination of magnetic field (including both strength and direction) and shear rate controls the mean number of magnetic clusters, which in turn is largely responsible for the static and dynamic behavior. While some properties are similar to comparable non-magnetic star polymers, others exhibit novel phenomena; examples of the latter include the breakup and reorganization of the clusters beyond a critical shear rate, and a strong dependence of the efficiency with which shear rate is translated into whole-body rotations on the direction of the magnetic field.

We present a new two dimensional model for elliptic (i.e., shape anisotropic) patchy colloids, where the impenetrable core of the particles is decorated on its co-)vertices by Kern-Frenkel type of patches. Using (i) well-documented criteria for the overlap of the undecorated ellipses and (ii) proposing new criteria that give evidence if the patchy regions of two particles interact we perform extensive Monte Carlo simulations in the NPT ensemble. Considering elliptic particles of aspect ratios \{kappa} = 2, 4, and 6 with two patches located at the co-vertices and choosing representative values for the temperature we study on a semi-quantitative level the emerging disordered phases in terms of snapshots and radial distribution functions and present results for the equation-of-state.

We present a comprehensive theoretical and experimental study of the solution phase properties of DNA-based family of nanoparticles - dendrimer-like DNA molecules (DL-DNA). These charged DNA dendrimers are novel macromolecular aggregates, which hold high promise in targeted self-assembly of soft matter systems in the bulk and at interfaces. To describe the behavior of this family of dendrimers (with generations ranging from G1 to G7), we use a theoretical model in which base-pairs of a single DL-DNA molecule are modeled by charged monomers, whose interactions are chosen to mimic the equilibrium properties of DNA correctly. Experimental results on the sizes and conformations of DL-DNA are based on static and 1dynamic light scattering; at the same time, Molecular Dynamics simulations are employed to model the equilibrium properties of DL-DNA, which compare favorably with the findings from experiments while at the same time providing a host of additional information and insight into the molecular structure of the nanostructures. We also examine the salt-responsiveness of these macromolecules, finding that despite the strong screening of electrostatic interactions, brought about by the added salt, the macromolecules shrink only slightly, their size robustness stemming from the high bending rigidity of the DNA-segments. The study of these charged dendrimer systems is an important field of research in the area of soft matter due to their potential role for various interdisciplinary applications, ranging from molecular cages and carriers for drug delivery in a living organism to the development of dendrimer- and dendron-based ultra-thin films in the area of nanotechnology. These findings are essential to determine if DL-DNA is a viable candidate for the experimental realization of cluster crystals in the bulk, a novel form of solids with multiple site occupancy., 38 pages, 17 figures

We provide rigorous evidence that the ordered ground state configurations of a system of parallel oriented, ellipsoidal particles, interacting via a Gaussian interaction (termed in literature as Gaussian core nematics) {\it must} be infinitely degenerate: we have demonstrated that these configurations originate from the related ground state configuration of the corresponding symmetric Gaussian core system via a suitable stretching operation of this lattice in combination with an arbitrary rotation. These findings explain related observations in former investigations, which then remained unexplained. Our conclusions have far reaching consequences for the search of ground state configurations of other nematic particles.

It is useful to describe the deformation characteristics of long biological macromolecules, such as deoxyribonucleic acid (DNA), by means of terms such as “bending”, “stretching”, or “twisting”. These terms are borrowed from classical beam theory, a traditional and widely known subfield of continuum mechanics, whereas the standard numerical modeling procedure for macromolecules, which is molecular dynamics, does not allow for explicit introduction of the aforementioned deformation modes. This somehow puts some limit to the mechanical understanding of biological macromolecules. As a remedy, we here propose an upscaling (or homogenization) approach, spanning a new conceptual bridge from molecular dynamics to beam theory. Firstly, we apply the principle of virtual power (PVP) to classical continuum beams subjected to stretching and twisting, as well as to atomic compounds represented as discrete systems of mass points in the framework of molecular dynamics. Equating virtual power densities associated with continuum and discrete representations provides homogenization rules from the atomic compounds to the continuum beam line elements. Secondly, the forces acting on the aforementioned mass points are derived from energy potentials associated with bond stretching, valence and torsion angle variations, as well as electrostatic and van der Waals interactions. Application of this strategy to a specific DNA sequence consisting of 20 base pairs reveals deformation-dependent conformational changes, as well as paradox phenomena such as “stretching due to overwinding”, in line with known experimental observations.

Inverse patchy colloids are patchy particles with differently charged surface regions. In this paper we focus on inverse patchy colloids with two different polar patches and an oppositely charged equatorial belt, and we describe a model and a reliable and efficient numerical algorithm that can be applied to investigate the properties of these particles in molecular dynamics simulations.

We present a comprehensive discussion of the so-called asymmetric Wigner bilayer system, where mobile point charges, all of the same sign, are immersed into the space left between two parallel, homogeneously charged plates (with possibly different charge densities). At vanishing temperatures, the particles are expelled from the slab interior; they necessarily stick to one of the two plates and form there ordered sublattices. Using complementary tools (analytic and numerical), we study systematically the self-assembly of the point charges into ordered ground state configurations as the inter-layer separation and the asymmetry in the charge densities are varied. The overwhelming plethora of emerging Wigner bilayer ground states can be understood in terms of the competition of two strategies of the system: net charge neutrality on each of the plates on the one hand and particles' self-organization into commensurate sublattices on the other hand. The emerging structures range from simple, highly commensurate (and thus very stable) lattices (such as staggered structures, built up by simple motives) to structures with a complicated internal structure. The combined application of our two approaches (whose results agree within remarkable accuracy) allows us to study on a quantitative level phenomena such as over- and underpopulation of the plates by the mobile particles, the nature of phase transitions between the emerging phases (which pertain to two different universality classes), and the physical laws that govern the long-range behaviour of the forces acting between the plates. Extensive, complementary Monte Carlo simulations in the canonical ensemble, which have been carried out at small, but finite temperatures along selected, well-defined pathways in parameter space confirm the analytical and numerical predictions within high accuracy. The simple setup of the Wigner bilayer system offers an attractive possibility to study and to control complex scenarios and strategies of colloidal self-assembly, via the variation of two system parameters.

We put forward a model that allows the calculation of the effective potential of two interacting DNA-stars, i.e., three-armed, Y-shaped, charged macromolecules, built up by three intertwined single-stranded DNAs. These particles are assumed to float on a flat interface separating two media with different dielectric properties. As the only input, our model requires the charge density along the branches and the interaction between two infinitesimally short segments, along two interacting rods. With this effective interaction at hand, a detailed investigations of the self-assembly scenarios of these molecules either via computer simulations or via theoretical frameworks comes within reach.

We consider a novel class of patchy particles inspired by polymer-based complex units where the limited valence in bonding is accompanied by soft interactions and incessant fluctuations of the patch positions, possibly leading to reversible modifications of the patch number and size. We introduce a simple model that takes into account the aforementioned features and we focus on the role played by the patch flexibility on the self-organization of our patchy units in the bulk, with particular attention to the connectivity properties and the morphology of the aggregated networks.

A planar array of identical charges at vanishing temperature forms a Wigner crystal with hexagonal symmetry. We take off one (reference) charge in a perpendicular direction, hold it fixed, and search for the ground state of the whole system. The planar projection of the reference charge should then evolve from a six-fold coordination (center of a hexagon) for small distances to a three-fold arrangement (center of a triangle), at large distances $d$ from the plane. The aim of this paper is to describe the corresponding non-trivial lattice transformation. For that purpose, two numerical methods (direct energy minimization and Monte Carlo simulations), together with an analytical treatment, are presented. Our results indicate that the $d=0$ and $d\to\infty$ limiting cases extend for finite values of $d$ from the respective starting points into two sequences of stable states, with intersecting energies at some value $d_t$; beyond this value the branches continue as metastable states., Comment: 17 pages, 11 figures

The performance of the previously validated LUCIO(®)-Direct-enzyme linked immunosorbent assay (direct ELISA) screening tests according to forensic guidelines is compared to that of cloned enzyme donor immunoassays (CEDIA) test for drugs of abuse in urine as defined in the new re-licensing German medical and psychological assessment (MPA) guidelines. The MPA screening cut-offs correspond to 10 ng/ml 11-nor-delta-9-tetrahydrocannabinol-9-carboxylic acid (THC-COOH), 50 ng/ml amphetamine and designer amphetamines, 25 ng/ml morphine, codeine and dihydrocodeine, 30 ng/ml benzoylecgonine, 50 ng/ml methadone metabolite, 2-ethylidene-1,5-dimethyl-3,3-diphenylpyrrolidine (EDDP) and metabolites of diazepam, oxazepam, bromazepam, alprazolam, flunitrazepam and lorazepam at 50 ng/ml. Average relative sensitivities and relative specificities were 99.7 % and 98.4 % for direct ELISA and 66 % and 91.4 % for CEDIA, respectively.

Systems with short-range attractive and long-range repulsive interactions are able to form mesophases at sufficiently low temperatures. In two dimensions, such mesophases emerge as clusters, stripes or bubbles. Using extensive Monte Carlo simulations we investigate the static and the dynamic properties of such a cluster-forming system over a broad temperature range and for different densities. Via the static properties we analyse how ordering into close packed configurations sets in both at the level of the particles as well as at the level of the clusters. The dynamic properties provide information on how, at low temperature, the motion of individual particles is influenced by the dynamic slowing down of the clusters. Finally, we discuss the different diffusion mechanisms at play at low and intermediate densities.

In Germany drink driving offenders lose their license and must prove abstinence for one year in order to regain it. In this paper we assess the newly introduced ethyl glucuronide (EtG) tests in urine and hair in this alcohol abstinence monitoring. 20% (80 out of 386) of the 3cm long hair samples were tested positive for EtG in hair, compared to only 2% (92 out of 4248 samples) in urine in the same time period. Additionally 50% of the samples positive for EtG in hair had EtG values greater than 30pg/mg hair, indicating chronic alcohol consumption in the last three months. This study shows that four EtG tests in 3cm hair lengths reveal a significantly higher percentage of drink driving offenders who fail to be sober in the rehabilitation period, than do six random EtG tests in urine. Presumably, the hair test is more adequate to monitor long term alcohol abstinence than the urine test as defined by the new driving license re-granting medical and psychological assessment (MPA) in Germany.

The purpose of the study was to compare the detection rate of illicit drugs in urine and hair specimens. The samples were taken from subjects trying to regain their revoked driver's license after a drug- or alcohol-related traffic offence. In 2010, we screened 14 000 urine and 3900 hair samples for amphetamines, methamphetamines, cannabinoids, cocaine, opiates, methadone, and benzodiazepines as well as for ethylglucuronide. We used the low threshold values of the new German guidelines for Medical Psychological Assessment (MPA). Positive screening tests were confirmed with gas chromatography-mass spectrometry (GC-MS), gas chromatography-tandem mass spectrometry (GC-MS/MS) or liquid chromatography-tandem mass spectrometry (LC-MS/MS). The results show that positivity rates for methamphetamines, MDMA, cocaine, and monoacetylmorphine were 1.7-, 5.7-, 3.8- and 9.3-fold higher in hair than in urine. In contrast, the detection rate for benzodiazepines was higher in urine than in hair (oxazepam, 0.21% versus 0%, nordiazepam 0.10% versus 0.03%). The positivity rate in hair for ethylglucuronide was 6-fold (12.7%) that for urine testing (2.1%). The study reveals that in the control of abstinence in the context of driving license re-granting there are in part large differences of positivity rates for some drugs or metabolites between hair and urine samples. These differences should be kept in mind by physicians and psychologists in traffic medicine who are ordering the drug testing. Copyright © 2012 John Wiley & Sons, Ltd. Language: en

In this paper we present the first assessment of the new German driving licence re-granting medical and psychological assessment (MPA) guidelines by comparing over 3500 urine samples tested under the old MPA cut-offs to over 5000 samples tested under the new MPA cut-offs. Since the enzyme multiplied immunoassay technique (EMIT) technology used previously was not sensitive enough to screen for drugs at such low concentrations, as suggested by the new MPA guidelines, enzyme-linked immunosorbent assay (ELISA) screening kits were used to screen for the drugs of abuse at the new MPA cut-offs. The above comparison revealed significantly increased detection rates of drug use or exposure during the rehabilitation period as follows: 1.61, 2.33, 3.33, and 7 times higher for 11-nor-delta-9-tetrahydrocannabinol-9-carboxylic acid (THC-COOH), morphine, benzoylecgonine and amphetamine respectively. The present MPA guidelines seem to be more effective to detect non-abstinence from drugs of abuse and hence to detecting drivers who do not yet fulfil the MPA requirements to regain their revoked driving licence.

International audience; Model systems in which fluid particles move in a disordered matrix of immobile obstacles have been found to be a reasonable representation of a colloidal fluid confined in a disordered porous medium. For systems consisting of hard-sphere particles, the obstacle matrix partitions the space available to the fluid particles into voids of finite volume ("traps") and a percolating void that extends over the entire volume. This geometric distinction plays a key role for the dynamic properties of the confined fluid: while its particles are not able to escape from traps, in the percolating void they can propagate infinitely far. We present a geometric method, based on a Delaunay decomposition, to identify the two different kinds of voids in an arbitrary matrix configuration of finite size under periodic boundary conditions. We subsequently apply a rastering technique, which enables us to statistically characterize the structure of the voids. We investigate the specific case of a quenched-annealed mixture of identical hard spheres, for which, among others, we accurately determine the matrix packing fraction at which the percolation transition of the voids takes place.

We introduce a continuum description of the thermodynamics of colloids with a core-corona architecture. In the case of thick coronas, their overlap can be treated approximately by replacing the exact one-particle density distribution by a suitably shaped step profile, which provides a convenient way of modeling the spherical, columnar, lamellar, and inverted cluster morphologies predicted by numerical simulations and the more involved theories. We use the model to study monodisperse particles with the hard-core/square-shoulder pair interaction as the simplest representatives of the core-corona class. We derive approximate analytical expressions for the enthalpies of the cluster morphologies which offer a clear insight into the mechanisms at work, and we calculate the lattice spacing and the cluster size for all morphologies of the phase sequence as well as the phase-transition pressures. By comparing the results with the exact crystalline minimum-enthalpy configurations, we show that the accuracy of the theory increases with shoulder width. We discuss possible extensions of the theory that could account for the finite-temperature effects.

This study attempts to assess the utility of the urine matrix as an alternative to cosmetically treated hair for the detection of drugs and alcohol for driving licence re-granting in 1026 cosmetically treated hair samples and 33 262 urine routine samples. No significant difference was observed between the percentage positive samples in cosmetically treated hair to those in urine at both the 95% and 99% significance level for amphetamines, cocaine, opiates, benzodiazepines, and methadone. Significant difference was found between the positivity rates of cannabinoids in cosmetically treated hair and that in urine indicating urine to be a better alternative to the use of the hair matrix even when cosmetically treated. The opposite was observed for the alcohol consumption marker ethyl glucuronide (EtG) for which the positivity rate in cosmetically treated hair was twice that in urine samples. Particularly for alcohol abstinence monitoring, as for the rehabilitative driving licence re-granting medical and psychological assessment (MPA) programme in Germany, it seems that ethyl glucuronide (EtG) in hair presents a much better alternative than urine testing, even when cosmetically treated hair is analyzed. Moreover, segmentation is an additional advantage of hair testing which can provide additional useful information. Copyright © 2014 John Wiley & Sons, Ltd. Language: en

We numerically investigate the formation of stable clusters of overlapping particles in certain systems interacting via purely repulsive, bounded pair potentials. In close vicinity of a first-order phase transition between a disordered and an ordered structure, clusters are encountered already in the fluid phase which then freeze into crystals with multiply occupied lattice sites. These hyper-crystals are characterized by a number of remarkable features that are in clear contradiction to our experience with harshly repulsive systems: upon compression, the lattice constant remains invariant, leading to a concomitant linear growth in the cluster population with density; further, the freezing and melting lines are to high accuracy linear in the density-temperature plane, and the conventional indicator that announces freezing, that is, the Hansen-Verlet value of the first peak of the structure factor, attains for these soft systems much higher values than for their hard-matter counterparts. Our investigations are based on the generalized exponential model of index 4 (i.e., Phi(r) approximately exp[-(r/sigma)4]). The properties of the phases involved are calculated via liquid state theory and classical density functional theory. Monte Carlo simulations for selected states confirm the theoretical results for the structural and thermodynamic properties of the system. These numerical data, in turn, fully corroborate an approximate theoretical framework that was recently put forward to explain the clustering phenomenon for systems of this kind (Likos, C. N.; Mladek, B. M.; Gottwald, D.; Kahl, G. J. Chem. Phys. 2007, 126, 224502).

We have studied the structural properties of a fluid of linear polymer chains under the external field of a disordered porous matrix, represented by a frozen liquid configuration of star polymer particles. To this end, and using a coarse-grained description of both fluid and matrix, we have applied the replica Ornstein-Zernike formalism to this particular case where soft and ultra-soft interparticle potentials have come into play. Our systematic investigations point out how the matrix influences the structural properties of the confined polymer fluid.

By compressing a (two-dimensional) system of ultrasoft, cluster-forming particles via a combined thermo- and barostat, a cluster crystal is created that is obviously not in its equilibrium: launching from such a configuration expansion and compression runs leads to systems that differ in their density distinctively from the one of the initial state. With our analysis we can discard dynamic lag as the origin of these discrepancies and can confirm that they occur due to the complex interplay between the two mechanisms which accompany any change in volume of a cluster crystal, namely the variation of the lattice spacing and of the number of lattice sites. In our investigations we show that only combined melting and annealing runs are able to transform the system in its equilibrium. This is evidenced by the fact that the density of this state remains unchanged even after the application of combined compression- and expansion-experiments as well as combined melting- and annealing-runs; furthermore, the cluster crystal shows an essentially perfect hexagonal arrangement.

Inverse patchy colloids (IPCs) differ from conventional patchy particles because their patches repel (rather than attract) each other and attract (rather than repel) the part of the colloidal surface that is free of patches. These particular features occur, e.g. in heterogeneously charged colloidal systems. Here we consider overall neutral IPCs carrying two, relatively small, polar patches. Previous studies of the same model under planar confinement have evidenced the formation of branched, disordered aggregates composed of ring-like structures. We investigate here the bulk behavior of the system via molecular dynamics simulations, focusing on both the structure and the dynamics of the fluid phase in a wide region of the phase diagram. Additionally, the simulation results for the static observables are compared to the Associative Percus Yevick solution of an integral equation approach based on the multi-density Ornstein-Zernike theory. A good agreement between theoretical and numerical quantities is observed even in the region of the phase diagram where the slowing down of the dynamics occurs.

The behavior of a fluid with competing interaction ranges adsorbed in a controlled pore size disordered matrix is studied by means of grand canonical Monte Carlo simulations in order to analyze the effects of confinement. The disordered matrix model is constructed from a two-dimensional non-additive hard-sphere fluid (which shows close to its demixing critical point large fluctuations in the concentration), after a subsequent quenching of the particle positions and removal of one of the components. The topology of the porous network is analyzed by means of a Delaunay tessellation procedure. The porous cavities are large enough to allow for cluster formation, which is however somewhat hindered as a result of the confinement, as seen from the comparison of cluster size distributions calculated for the fluid under confinement and in the bulk. The occurrence of lamellar phases is impeded by the disordered nature of the porous network. Analysis of two-dimensional density maps of the adsorbed fluid for given matrix configurations shows that clusters tend to build up in specific locations of the porous matrix, so as to minimize inter-cluster repulsion.

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., 24 pages, 2 figures, Review Article

Integral equation theory calculations within the mean spherical approximation (MSA) and grand canonical Monte Carlo (MC) simulations are employed to study the phase behaviour of a symmetrical binary fluid mixture in the presence of a field arising from unequal chemical potentials of the two particle species. Attention is focused on the case for which, in the absence of a field, the phase diagram exhibits a first order liquid-liquid transition in addition to the liquid-vapor transition. We find that in the presence of a field, two possible subtypes of phase behaviour can occur, these being distinguished by the relationship between the critical lines in the full phase diagram of temperature, density, and concentration. We present the detailed form of the respective phase diagrams as calculated from MSA and compare with results from the MC simulations, finding good overall agreement., 8 pages, 2 figures

In an effort to approach a quantitative determination of the phase diagram of polydisperse mixtures of (possibly size-asymmetric) charged hard spheres (CHS), we propose a concept that explicitly takes into account the association of the ions into dimers, i.e., an effect that is known to be of high relevance close to the liquid-vapor phase boundary in the restricted primitive model (RPM). We use the polymer mean spherical approximation (PMSA) to calculate the properties of this polydisperse mixture of dimers. The concept turns out to be a truncatable free energy model, thus the phase diagram can be calculated by solving a highly non-linear set of equations. We present the full phase diagram (in terms of cloud and shadow curves and binodals) and discuss fractionation effects for a model system.

We present results of a systematic investigation of the properties of the generalized Gaussian core model of index n. The potential of this system interpolates via the index n between the potential of the Gaussian core model and the penetrable sphere system, thereby varying the steepness of the repulsion. We have used both conventional and self-consistent liquid state theories to calculate the structural and thermodynamic properties of the system; reference data are provided by computer simulations. The results indicate that the concept of self-consistency becomes indispensable to guarantee excellent agreement with simulation data; in particular, structural consistency (in our approach taken into account via the zero separation theorem) is obviously a very important requirement. Simulation results for the dimensionless equation of state, βP/%, indicate that for an indexvalue of 4 a clustering transition, possibly into a structurally ordered phase might set in as the system is compressed. Ми представляємо результати системного дослідження властивостей узагальненої моделі гаусового кора з індексом n. Потенціал такої системи з допомогою індексу n, що визначає крутизну притягання, дозволяє інтерполювати взаємодії від моделі гаусового кора до системи проникних сфер. Для розрахунку структурних і термодинамічних властивостей нами використовувалися як традиційні, так і самоузгоджені версії теорії рідкого стану, а система відліку бралася з комп’ютерного експерименту. Результати показують, що для отримання доброго узгодження з даними моделювання концепція самоузгодження стає обов’язковою; зокрема структурне узгодження, яке у нашому підході проводиться через теорему нульового розділення, є надзвичайно важливою вимогою. Результати моделювання для обезрозміреного рівняння стану, βP/%, вказують на те, що для моделі з індексом 4 може виникати структурно впорядкована фаза (перехід кластерування) при стисканні системи.

A review is given on recent studies of charged colloidal suspensions and polyelectrolytes both in static and non-equilibrium situations. As far as static equilibrium situations are concerned, we discuss three different problems: 1) Sedimentation density profiles in charged suspensions are shown to exhibit a stretched non-bariometric wing at large heights and binary suspensions under gravity can exhibit an analog of the brazil-nut effect known from granular matter, i.e. the heavier particles settle on top of the lighter ones. 2) Soft polyelectrolyte systems like polyelectrolyte stars and microgels show an ultra-soft effective interaction and this results into an unusual equilibrium phase diagram including reentrant melting transitions and stable open crystalline lattices. 3) The freezing transition in bilayers of confined charged suspensions is discussed and a reentrant behaviour is obtained. As far as nonequilibrium problems are concerned, we discuss an interface instability in oppositely driven colloidal mixtures and discuss possible approaches to simulate electrokinetic effects in charged suspensions.

We have investigated the phase behaviour of a symmetric binary mixture with particles interacting via hard-core Yukawa potentials. To calculate the thermodynamic properties we have used the mean spherical approximation (MSA), a conventional liquid state theory, and the closely related self-consistent Ornstein–Zernike approximation which is defined via an MSA-type closure relation, requiring, in addition, thermodynamic self-consistency between the compressibility and the energy-route. We investigate on a quantitative level the effect of the self-consistency requirement on the phase diagram and on the critical behaviour and confirm the existence of three archetypes of phase diagram, which originate from the competition between the first order liquid/vapour transition and the second order demixing transition.

The solution of the mean spherical approximation for a polydisperse fluid mixture of particles interacting via a set of (factorizable) hard-sphere Yukawa potentials allows one to represent those thermodynamic quantities that are relevant to determine phase coexistence (i.e., pressure and chemical potential) by a limited number of (generalized) moments. Being thus a member of “truncatable free energy models,” the equilibrium conditions reduce to a set of coupled and highly nonlinear equations; we have solved these relations, we have determined phase diagrams for polydisperse fluid mixtures (i.e., cloud and shadow curves as well as binodals), and have extracted the daughter distribution functions of the coexisting phases.

It is by now well established that the self-consistent Ornstein-Zernike approximation (SCOZA) is able to provide accurate data for the location of the liquid-vapour coexistence curve and for the critical point of simple fluids. However, up to now applications of the SCOZA were restricted only to a rather small variety of interatomic potentials. In this contribution we present an extension of the SCOZA which offers access to a considerably larger class of systems. As an example, we present results for model fullerenes Cn ≥ 60 (characterised by the sphericalized Girifalco potentials) and thus contribute to the still open question on the existence of the liquid phase and on the size of the liquid pocket in the phase diagram of fullerenes.

We present expressions that allow the determination of the structural and thermodynamic properties of a polydisperse liquid mixture in contact with a polydisperse matrix; polydispersity of the fluid and of the matrix are described by distribution functions f f(σ) and f m (σ), where σ and σ characterize the size of the fluid and of the matrix particles. The formalism is based on the replica trick (to describe the properties of a fluid in contact with a matrix) and on the expansion of size-dependent functions in terms of orthogonal polynomials associated with the distribution functions f f(σ) and f m (σ). In our expressions structural and thermodynamic properties are calculated from coefficient functions (of these expansions) which are obtained from a numerical solution of the generalized replica Ornstein-Zernike equations, solved along with a suitable closure relation.

We develop a density functional for hard-sphere mixtures which keeps the structure of Rosenfeld's fundamental measure theory (FMT) whilst inputting the Mansoori–Carnahan–Starling–Leland bulk equation of state. Density profiles for the pure hard-sphere fluid and for some binary mixtures adsorbed at a planar hard wall obtained from the present functional exhibit some improvement over those from the original FMT. The pair direct correlation function c(2) (r) of the pure hard-sphere fluid, obtained from functional differentiation, is also improved. When a tensor weight function is incorporated for the pure system our functional yields a good description of fluid–solid coexistence and of the properties of the solid phase.