Claudio Pastorino

Static and dynamic properties of the interface between a polymer brush and a melt of identical chains

Molecular-dynamics simulations of a short-chain polymer melt between two brush-covered surfaces under shear have been performed. The end-grafted polymers which constitute the brush have the same chemical properties as the free chains in the melt and provide a soft deformable substrate. Polymer chains are described by a coarse-grained bead-spring model, which includes excluded volume and backbone connectivity of the chains. The grafting density of the brush layer offers a way of controlling the behavior of the surface without altering the molecular interactions. We perform equilibrium and nonequilibrium molecular-dynamics simulations at constant temperature and volume using the dissipative particle dynamics thermostat. The equilibrium density profiles and the behavior under shear are studied as well as the interdigitation of the melt into the brush, the orientation on different length scales (bond vectors, radius of gyration, and end-to-end vector) of free and grafted chains, and velocity profiles. The obtained boundary conditions and slip length show a rich behavior as a function of grafting density and shear velocity.

Comparison of Ternary Bilayer Mixtures with Asymmetric or Symmetric Unsaturated Phosphatidylcholine Lipids by Coarse Grained Molecular Dynamics Simulations

We studied the phase behavior of various ternary bilayer mixtures composed of cholesterol, an unsaturated lipid, and a fully saturated lipid, by means of molecular dynamics simulations of the MARTINI coarse grain model. We aimed at comparing lateral organization and local properties of bilayers containing phosphatidylcholine (PC) lipids, either with two unsaturated tails (symmetric), or one unsaturated and one saturated tail (asymmetric), as the low-melting component of the mixture. The number of unsaturations per chain was systematically varied in both classes of unsaturated lipids, to account for its consequences in segregation. In the asymmetric unsaturated PCs, the saturated tail was kept identical to the hydrophobic chains of the fully saturated lipid component. Membranes with a symmetric or an asymmetric unsaturated lipid, with the same kind of unsaturated chain, show different phase behavior. Symmetric polyunsaturated PCs set the separation in two phases. Instead, the asymmetric polyunsaturated lipids induced nonideal mixing of components in single-phase bilayers. A significative drop of temperature, within the accessible temperature range, enhances the segregation in mixtures with the more unsaturated asymmetric PC, but still within a single phase. This different phase behavior between membranes with symmetric and asymmetric unsaturated PCs is also observed for lipids with the same total number of unsaturations. On the other hand, the degree of unsaturation per se enhances the segregation, by increasing the composition fluctuations in single-phase membranes with asymmetric PC lipids, and raising the line tension in the two-phase bilayer mixtures with symmetric polyunsaturated PCs. Dynamic clusters of unsaturated asymmetric lipids can be identified. The clusters show no correlation between leaflets, as observed for the phase domains in mixtures with the symmetric polyunsaturated PCs. Interestingly, we found that asymmetric PC lipids have a preferential orientation such that their saturated tails increase their density toward the periphery of the clusters, facing regions enriched in the fully saturated lipids and cholesterol. The degree of unsaturation increases the cluster size and also enhances the anisotropy of the orientation. The surface density of cholesterol follows a gradient that favors its interaction with the saturated tails. Such gradients in composition lead to gradients in order parameters, such as the conformational order and the area of the tails, which increases away from the unsaturated lipid clusters. We compared, in addition, differences in hydrophobic length mismatch between acyl chains of the low-melting and high-melting components, in mixtures containing either symmetric or asymmetric unsaturated lipids.

Polyunsaturated and Saturated Phospholipids in Mixed Bilayers: A Study from the Molecular Scale to the Lateral Lipid Organization

Polyunsaturated lipids are remarkably flexible molecules, with a great influence on the membrane structure and dynamics, affecting from mechanical properties to domain segregation. In this work, we studied phospholipid mixtures of dipalmitoylphosphatidylcholine (DPPC) and diunsaturated phosphatidylcholine lipids (diunsPC) of different lengths, by means of molecular dynamic simulations of a coarse-grained interaction model. These diunsPC:DPPC binary mixtures show nonideal behavior characterized by one mixed phase with composition fluctuations on a length scale of nanometers. Motivated by this observation, we studied comprehensively the characteristics of molecular structure as a function of the compositional gradient. We analyzed orientational order profiles, density distributions, and pair-pair correlation functions between the molecule residues. We observed that, in diunsPC-enriched regions, DPPC tails become expanded and disordered, especially toward the membrane center. On the other hand, in the more condensed DPPC-enriched patches, diunsaturated acyl chains become displaced toward the interface instead of stretching along the membrane normal. From the comparison of the two diunsPC lipids of different tail length, we measured that the presence of a longer terminal saturated segment induces better mixing with DPPC, and most interestingly eliminates the up-down composition correlation measured with the shorter tail-diunsPC. At molecular level, there is a reduced redistribution of densities and changes in the local order as a function of composition. We interpret these results as indicative that the packing incompatibility between polyunsaturated and saturated lipids rules their mixing behavior.

Dissipative Particle Dynamics: A Method to Simulate Soft Matter Systems in Equilibrium and Under Flow

We provide examples and a concise review of the method of Dissipative Particle Dynamics (DPD), as a simulation tool to study soft matter systems and simple liquids in equilibrium and under flow. DPD was initially thought as a simulation method, which in combination with soft potentials, could simulate “fluid particles” with suitable hydrodynamic correlations. Then DPD evolved to a generic “thermostat” to simulate systems in equilibrium and under flow, with arbitrary interaction potential among particles. We describe the application of the method with soft potentials and other coarse-grain models usually used in polymeric and other soft matter systems. We explain the advantages, common problems and limitations of DPD, in comparison with other thermostats widely used in simulations. The implementation of the DPD forces in a working Molecular Dynamics (MD) code is explained, which is a very convenient property of DPD. We present various examples of use, according to our research interests and experiences, and tricks of trade in different situations. The use of DPD in equilibrium simulations in the canonical ensemble, the grand canonical ensemble at constant chemical potential, and stationary Couette and Poiseuille flows is explained. It is also described in detail the use of different interaction models for molecules: soft and hard potentials, electrostatic interactions and bonding interactions to represent polymers. We end this contribution with our personal views and concluding remarks.

Brushes of semiflexible polymers in equilibrium and under flow in a super-hydrophobic regime

We performed molecular dynamics simulations to study the equilibrium and flow properties of a liquid in a nano-channel with confining surfaces coated with a layer of grafted semiflexible polymers. The coverage spans a wide range of grafting densities from essentially isolated chains to dense brushes. The end-grafted polymers were described by a bead spring model with a harmonic potential to include the bond stiffness of the chains. We varied the rigidity of the chains, from fully flexible polymers to rigid rods, in which the configurational entropy of the chains is negligible. The brush-liquid interaction was tuned to obtain a super-hydrophobic channel, in which the liquid did not penetrate the polymer brush, giving rise to a Cassie-Baxter state. Equilibrium properties such as brush height and bending energy were measured, varying the grafting density and the stiffness of the polymers. We also studied the characteristics of the brush-liquid interface and the morphology of the polymer chains supporting the liquid for different bending rigidities. Non-equilibrium simulations were performed, moving the walls of the channel in opposite directions at constant speed, obtaining a Couette velocity profile in the bulk liquid. The molecular degrees of freedom of the polymers were studied as a function of the Weissenberg number. Also, the violation of the no-slip boundary condition and the slip properties were analyzed as a function of the shear rate, grafting density and bending stiffness. At high grafting densities, a finite slip length independent of the shear rate or bending constant was found, while at low grafting densities a very interesting non-monotonic dependence on the bending constant is observed.

Measuring the composition-curvature coupling in binary lipid membranes by computer simulations

The coupling between local composition fluctuations in binary lipid membranes and curvature affects the lateral membrane structure. We propose an efficient method to compute the composition-curvature coupling in molecular simulations and apply it to two coarse-grained membrane models-a minimal, implicit-solvent model and the MARTINI model. Both the weak-curvature behavior that is typical for thermal fluctuations of planar bilayer membranes as well as the strong-curvature regime corresponding to narrow cylindrical membrane tubes are studied by molecular dynamics simulation. The simulation results are analyzed by using a phenomenological model of the thermodynamics of curved, mixed bilayer membranes that accounts for the change of the monolayer area upon bending. Additionally the role of thermodynamic characteristics such as the incompatibility between the two lipid species and asymmetry of composition are investigated.

Mean properties and free energy of a few hard spheres confined in a spherical cavity

We use analytical calculations and event-driven molecular dynamics simulations to study a small number of hard sphere particles in a spherical cavity. The cavity is also taken as the thermal bath so that the system thermalizes by collisions with the wall. In that way, these systems of two, three, and four particles, are considered in the canonical ensemble. We characterize various mean and thermal properties for a wide range of number densities. We study the density profiles, the components of the local pressure tensor, the interface tension, and the adsorption at the wall. This spans from the ideal gas limit at low densities to the high-packing limit in which there are significant regions of the cavity for which the particles have no access, due the conjunction of excluded volume and confinement. The contact density and the pressure on the wall are obtained by simulations and compared to exact analytical results. We also obtain the excess free energy for N = 4, by using a simulated-assisted approach in which we combine simulation results with the knowledge of the exact partition function for two and three particles in a spherical cavity.

Droplet Transport in a Nanochannel Coated by Hydrophobic Semiflexible Polymer Brushes: The Effect of Chain Stiffness

We study the influence of chain stiffness on droplet flow in a nanochannel, coated with semiflexible hydrophobic polymers by means of nonequilibrium molecular dynamics simulations. The studied system is then a moving droplet in the slit channel, coexisting with its vapor and subjected to periodic boundary conditions in the flow direction. The polymer chains, grafted by the terminal bead to the confining walls, are described by a coarse-grained model that accounts for chain connectivity, excluded volume interactions and local chain stiffness. The rheological, frictional and dynamical properties of the brush are explored over a wide range of persistence lengths. We find a rich behavior of polymer conformations and concomitant changes in the friction properties over the wide range of studied polymer stiffnesses. A rapid decrease in the droplet velocity was observed as the rigidity of the chains is increased for polymers whose persistence length is smaller than their contour length. We find a strong relation between the internal dynamics of the brush and the droplet transport properties, which could be used to tailor flow properties by surface functionalization. The monomers of the brush layer, under the droplet, present a collective treadmill belt like dynamics which can only be present due the existence of grafted chains. We describe its changes in spatial extension upon variations of polymer stiffness, with bidimensional velocity and density profiles. The deformation of the polymer brushes due to the presence of the droplet is analyzed in detail. Lastly, the droplet-gas interaction is studied by varying the liquid to gas ratio, observing a 16% speed increase for droplets that flow close to each other, compared to a train of droplets that present a large gap between consecutive droplets.

Molecular conformation of linear alkane molecules: From gas phase to bulk water through the interface

We studied the behavior of long chain alkanes (LCAs) as they were transferred from gas to bulk water, through the liquid-vapor interface. These systems were studied using umbrella sampling molecular dynamics simulation and we have calculated properties like free energy profiles, molecular orientation, and radius of gyration of the LCA molecules. The results show changes in conformation of the solutes along the path. LCAs adopt pronounced molecular orientations and the larger ones extend appreciably when partially immersed in the interface. In bulk water, their conformations up to dodecane are mainly extended. However, larger alkanes like eicosane present a more stable collapsed conformation as they approach bulk water. We have characterized the more probable configurations in all interface and bulk regions. The results obtained are of interest for the study of biomatter processes requiring the transfer of hydrophobic matter, especially chain-like molecules like LCAs, from gas to bulk aqueous systems through the interface.

Structure, thermodynamic properties, and phase diagrams of few colloids confined in a spherical pore

We study a system of few colloids confined in a small spherical cavity with event driven molecular dynamics simulations in the canonical ensemble. The colloidal particles interact through a short range square-well potential that takes into account the basic elements of attraction and excluded-volume repulsion of the interaction among colloids. We analyze the structural and thermodynamic properties of this few-body confined system in the framework of inhomogeneous fluids theory. Pair correlation function and density profile are used to determine the structure and the spatial characteristics of the system. Pressure on the walls, internal energy, and surface quantities such as surface tension and adsorption are also analyzed for a wide range of densities and temperatures. We have characterized systems from 2 to 6 confined particles, identifying distinctive qualitative behavior over the thermodynamic plane T - ρ, in a few-particle equivalent to phase diagrams of macroscopic systems. Applying the extended law of corresponding states, the square well interaction is mapped to the Asakura-Oosawa model for colloid-polymer mixtures. We link explicitly the temperature of the confined square-well fluid to the equivalent packing fraction of polymers in the Asakura-Oosawa model. Using this approach, we study the confined system of few colloids in a colloid-polymer mixture.