Kinetic structures in multi-phase liquid-crystalline composite materials

A mesoscopic kinetic model for phase separation in the presence of liquid crystalline order has been formulated and solved using high performance numerical methods. The thermodynamic phase diagram on temperature–polymer concentration plane indicates the presence of coexistence regions between isotropic and liquid crystalline phases. These regions are partitioned by the phase-separation spinodal and the phase-ordering spinodal. We characterize the morphologies following temperature quenches in the phase diagram. The scenario is completely different from isotropic mixing since the continuous phase exhibits liquid crystalline ordering. Microdomains of the dispersed phase induce long- and short-range forces affecting the kinetics of the phase separation and the emerging structures. Presence of topological defects and elastic distortions around the microdomains formed during the phase separation dominate the morphology. The free energy of the system establishes dynamics and correlations of the morphological structures.     

A molecular dynamics simulation on the effect of different parameters on thermal resistance of graphene-argon interface

Molecular dynamics simulations of argon molecules confined between two parallel graphene sheets are carried out to investigate the parameters affecting heat transfer and thermal properties. These parameters include wall–fluid interaction strength, fluid density and wall temperature. For constant wall temperature simulations, we show that the first two parameters have influence on near-wall fluid density. As a result, the heat transfer at wall–fluid interfaces and thus through argon molecules across the domain will change. Also, we demonstrate that variations in wall temperature rarely affects the density profiles of argon molecules next to the walls. Therefore, in these cases, the variations in thermal resistance at the interface is most dominantly due to wall temperature itself. To analyse the results, the density and temperature profiles and also other parameters including heat flux and temperature gradient of bulk of argon molecules, Kapitza length and argon thermal conductivity are considered. The Kapitza length describes thermal resistance at liquid–solid interface. According to the results, increasing wall–fluid interaction strength leads to greater molecular aggregation of argon molecules near the walls and, consequently, decreasing the Kapitza length. Furthermore, higher fluid density leads to greater thermal resistance at wall–fluid interactions and therefore greater temperature jumps are observed in temperature profiles.     

Packing characteristics of a model system mimicking cytoplasmic bacterial membranes.

The phase diagram of fully hydrated mixtures of dipalmitoylphosphatidylethanolamine and -phosphatidylglycerol was constructed and the coexistence lines of the solidus and liquidus curve calculated based on regular solution theory using two nonideality parameters for each of the phase to account for nonideal and nonsymmetric mixing. Both lipids show nonideal miscibility in the liquid-crystalline phase, while a region of immiscibility exists in the lamellar-gel phase between the mole fraction x(DPPE)=0.05-0.4. Two lines of three-phase coexistence around 35 and 40 degrees C reflects the presence of lipid domains predominantly composed of phosphatidylglycerol as well as of the mixed lipid system. This is reflected in the positive nonideality parameters of the gel phase obtained from the simulation of the phase diagram. Moreover, segregation of pure phosphatidylethanolamine domains was detected in mixtures x(DPPE)>0.9, which formed multilamellar liposomes, while unilamellarity was observed for the mixed lipid systems owing to the presence of the negatively charged phosphatidylglycerol. The packing constraints of these phospholipids, major components of cytoplasmic bacterial membranes, may be of importance in the interaction with various solutes like antimicrobial peptides, and were explained based on the nature of the headgroups and the molecular geometry of the phospholipids.     

Effect of surface-screening parameter of the Yukawa potential model on vapour–liquid phase coexistence and critical-point properties of confined Yukawa fluid

We report the effect of surface-screening parameter of Yukawa potential model on vapour–liquid phase coexistence and critical-point properties of slit–pore-confined Yukawa fluid, using grand canonical transition-matrix Monte Carlo along with the histogram reweighting method. The effect of surface-screening parameter on the vapour–phase coexistence density is insignificant for the studied system. On the other hand, significant effect of surface-screening parameter is observed on liquid phase coexistence density. With increasing surface-screening parameter, liquid phase coexistence density decreases. Critical-point properties have shown monotonic decreasing trends with increase in surface-screening parameter. Moreover, the effect of change of surface-screening parameter is least on critical temperature changes as compared to critical density and critical pressure changes for the studied Yukawa system in this work.     

The effect of ergosterol on dipalmitoylphosphatidylcholine bilayers: a deuterium NMR and calorimetric study

We have studied the effect of ergosterol, an important component of fungal plasma membranes, on the physical properties of dipalmitoylphosphatidylcholine (DPPC) multibilayers using deuterium nuclear magnetic resonance ((2)H NMR) and differential scanning calorimetry (DSC). For the (2)H NMR experiments the sn-1 chain of DPPC was perdeuterated and NMR spectra were taken as a function of temperature and ergosterol concentration. The phase diagram, constructed from the NMR spectra and the DSC thermograms, exhibits both solid-ordered (so) + liquid-ordered (lo) and liquid-disordered (ld) + lo phase coexistence regions with a clear three-phase line. This is the first demonstration that lo domains exist in liquid crystalline membranes containing ergosterol. The domain sizes in the ld+lo phase coexistence region were estimated by analyzing the exchange of labeled DPPC between the two regions, and depend on ergosterol concentration. The DPPC-ergosterol phase diagram is similar to that of the DPPC-cholesterol multibilayer system except that the so+lo and ld+lo phase coexistence regions are considerably broader.     

Temperature-pressure phase diagram of a heterogeneous anionic model biomembrane system: results from a combined calorimetry, spectroscopy and microscopy study

By using Fourier transform infrared (FT-IR) spectroscopy in combination with differential scanning calorimetry (DSC) coupled with pressure perturbation calorimetry (PPC), ultrasound velocimetry, Laurdan fluorescence spectroscopy, fluorescence microscopy and atomic force microscopy (AFM), the temperature and pressure dependent phase behavior of the five-component anionic model raft lipid mixture DOPC/DOPG/DPPC/DPPG/cholesterol (20:5:45:5:25 mol%) was investigated. A temperature range from 5 to 65 °C and a pressure range up to 16 kbar were covered to establish the temperature-pressure phase diagram of this heterogeneous model biomembrane system. Incorporation of 10-20 mol% PG still leads to liquid-ordered (l(o))-liquid-disordered (l(d)) phase coexistence regions over a wide range of temperatures and pressures. Compared to the corresponding neutral model raft mixture (DOPC/DPPC/Chol 25:50:25 mol%), the p,T-phase diagram is - as expected and in accordance with the Gibbs phase rule - more complex, the phase sequence as a function of temperature and pressure is largely similar, however. This anionic heterogeneous model membrane system will serve as a more realistic model biomembrane system to study protein interactions with anionic lipid bilayers displaying liquid-disordered/liquid-ordered domain coexistence over a wide range of the temperature-pressure plane, thus allowing also studies of biologically relevant systems encountered under extreme environmental conditions.     

Richards-like two species population dynamics model

The two-species population dynamics model is the simplest paradigm of inter- and intra-species interaction. Here, we present a generalized Lotka–Volterra model with intraspecific competition, which retrieves as particular cases, some well-known models. The generalization parameter is related to the species habitat dimensionality and their interaction range. Contrary to standard models, the species coupling parameters are general, not restricted to non-negative values. Therefore, they may represent different ecological regimes, which are derived from the asymptotic solution stability analysis and are represented in a phase diagram. In this diagram, we have identified a forbidden region in the mutualism regime, and a survival/extinction transition with dependence on initial conditions for the competition regime. Also, we shed light on two types of predation and competition: weak, if there are species coexistence, or strong, if at least one species is extinguished.     

Structure of Lennard–Jones fluid near large spherical particles: further test of the “universality” of adjustable parameter in perturbation density functional theory

Grand canonical Monte Carlo simulation is used to study the density profiles of Lennard–Jones (LJ) fluid next to a large hard sphere (mimicking a colloidal particle) of various sizes. The LJ fluid in the inhomogeneous system thus maintains equilibrium with the bulk LJ fluid. The chosen density and potential parameters for the bulk fluid correspond to the conditions situated at “dangerous” regions of the phase diagram, i.e. near the critical temperature or close to the gas–liquid coexistence curve. The aim of present extensive simulations is to provide exact data for the broad range of the bulk parameters against which the “universality” of adjustable parameter associated with a perturbation density functional approximation (DFA) can be tested. Here the term “universality” means independence of this parameter on the particular external field responsible for the generation of a non-uniform density profile of the fluid. It is shown that the “universality” of this parameter associated with a third order+second order perturbation DFA holds also in the present case of a large spherical particle as a source of external potential, similarly as established in previous studies dealing with other interaction potential and other external fields [J. Chem. Phys., 122, 064503 (2005); J. Chem. Phys., 123 124708 (2005)]. This DFA can be used as input into a recently proposed framework for the calculation of interparticle potential of mean force (PMF).     

Molecular dynamics simulation of thermal conductivity of an argon liquid layer confined in nanospace

Molecular dynamics simulations were used to study the thermal conductivity of liquid argon ultra thin films confined between two plates spaced several nanometres apart. The research focused on the dependence of the liquid argon thermal conductivity on the liquid layer thickness and the interaction between liquid and solid. The results show that the thermal conductivity of liquid argon ultra thin films confined between two plates depends on the distance between the two plates and the existence of solid-like liquid layering at the liquid–solid interface and the average migration frequency of all liquid molecules. Stronger interactions between the liquid and the solid resulted in a larger number of atoms in the solid-like liquid layer along the surface and hence smaller thermal resistance between the liquid and the solid. However, as the strength of the interaction with the solid increased, the thermal conductivity was reduced due to fewer atoms near the hot solid boundary and less molecular migration.     

Vapor–liquid equilibria and thermophysical behavior of the SPC-HW model for heavy water

The vapor–liquid coexistence curve of the simple point charge heavy-water model (SPC-HW), [J. Chem. Phys., 114, 8064–8067 (2001)] is determined by Gibbs Ensemble Monte-Carlo (GEMC) simulation. The estimated critical conditions of the model based on the Wegner-type expansion for the order parameters and the rectilinear diameter are ρ_c = 0.300 g/cc, T _c = 661 K and P _c = 156 bars. The dielectric constant determined by isothermal–isochoric molecular dynamics is underpredicted along the coexistence curve by 29–44% in comparison with the experimental values. The analysis of the orthobaric temperature dependence of the system microstructure, in terms of the three site–site radial distribution functions, indicates that the first coordination numbers for the oxygen–oxygen and the oxygen–deuterium interactions are ～4.3 ± 0.1 and ～1.9 ± 0.1 at T = 300 K, and decrease by 15 and 55%, respectively, at criticality. The dipole–dipole correlation functions show that the orientational order in heavy water is quickly lost beyond the first oxygen–oxygen coordination shell. The model's second virial coefficient is determined by Monte-Carlo integration and used to aid the interpretation of the predicted phase equilibrium results.     

Sphingomyelin/phosphatidylcholine/cholesterol phase diagram: boundaries and composition of lipid rafts

The ternary system palmitoylsphingomyelin (PSM)/palmitoyloleoylphosphatidylcholine (POPC)/cholesterol is used to model lipid rafts. The phase behavior of the three binary systems PSM/POPC, PSM/cholesterol, and POPC/cholesterol is first experimentally determined. Phase coexistence boundaries are then determined for ternary mixtures at room temperature (23 degrees C) and the ternary phase diagram at that temperature is obtained. From the diagram at 23 degrees C and the binary phase diagrams, a reasonable expectation is drawn for the ternary phase diagram at 37 degrees C. Several photophysical methodologies are employed that do not involve detergent extraction, in addition to literature data (e.g., differential scanning calorimetry) and thermodynamic rules. For the ternary phase diagrams, some tie-lines are calculated, including the one that contains the PSM/POPC/ cholesterol 1:1:1 mixture, which is often used in model raft studies. The diagrams here described are used to rationalize literature results, some of them apparently discrepant, and to discuss lipid rafts within the framework of liquid-ordered/liquid-disordered phase coexistence.     

Molecular simulation of volume of mixing,Helmholtz free energy of mixing and entropy of mixing in bulk fluid mixtures

We describe a kinetic Monte Carlo molecular simulation procedure to calculate the Helmholtz free energy, the entropy and the chemical potentials of all components in a bulk fluid mixture. This allows us to derive the excess properties (volume, free energy and entropy) resulting from the mixing of homogeneous fluids of pure components at constant temperature and pressure. We have chosen neon–xenon mixtures to illustrate our method because of the large difference in collision diameter and well-depth of the interaction energy. When xenon is predominant in the mixture, the volume of mixing is larger. The excess entropy of mixing correlates with the volume of mixing, since a positive excess volume enables more configurations (more possible molecular distributions). The excess thermodynamic quantities as functions of the total density were found to be insensitive to temperature. To investigate the effects of the molecular parameters, we also studied argon–nitrogen and argon–krypton mixtures. The effect of the difference in molecular parameters is in the order: argon–nitrogen < argon–krypton < neon–xenon. A large difference in the well-depth of the interaction energies results in an increase in the excess thermodynamic variables, which is in agreement with the literature McDonald IR. NpT-ensemble Monte Carlo calculations for binary liquid mixtures. Mol Phys. 1972;23(1):41–58; Singer JVL, Singer K. Monte Carlo calculation of thermodynamic properties of binary mixtures of Lennard-Jones (12-6) liquids. Mol Phys. 1972;24(2):357–390.     

Modeling of crude oil biodegradation using two phase partitioning bioreactor

In this work, crude oil biodegradation has been optimized in a solid‐liquid two phase partitioning bioreactor (TPPB) by applying a response surface methodology based d ‐optimal design. Three key factors including phase ratio, substrate concentration in solid organic phase, and sodium chloride concentration in aqueous phase were taken as independent variables, while the efficiency of the biodegradation of absorbed crude oil on polymer beads was considered to be the dependent variable. Commercial thermoplastic polyurethane (Desmopan®) was used as the solid phase in the TPPB. The designed experiments were carried out batch wise using a mixed acclimatized bacterial consortium. Optimum combinations of key factors with a statistically significant cubic model were used to maximize biodegradation in the TPPB. The validity of the model was successfully verified by the good agreement between the model‐predicted and experimental results. When applying the optimum parameters, gas chromatography‐mass spectrometry showed a significant reduction in n‐alkanes and low molecular weight polycyclic aromatic hydrocarbons. This consequently highlights the practical applicability of TPPB in crude oil biodegradation. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 30:797–805, 2014     

An effective way to determine the melting curve

A new ‘two-phase’ simulation method with which to accurately predict the melting curve is proposed. The method requires, as a starting configuration, generating a two-phase coexistence state by employing a suitable ensemble. Examining a change in volume ratio of the two phases upon varying temperature (pressure) under a fixed pressure (temperature) allows us to determine the phase transition point. The Clausius–Clapeyron relationship can then be implemented as a guide to predict the nearby phase transition point. The method was applied to determine the solid–liquid phase boundary of the modified Lennard–Jones system. A better accuracy, as that achieved by the non-equilibrium relaxation method (Asano Y, Fuchizaki K. J Phys Soc Jpn. 2017;86:025001), was obtained but with much less computational cost.     

Corresponding states vapour–liquid phase equilibria of confined square–well fluid

Corresponding states vapour–liquid phase equilibria of confined square-well fluid are studied by means of grand-canonical transition-matrix Monte Carlo simulation and histogram reweighting method. In this study, square-well fluid is considered under hard and attractive slit pore confinements ranging from 1.5 to 40 molecular diameters. Corresponding states vapour–liquid phase coexistence envelopes display insignificant effect of wall?fluid interaction for slit pore confinements ranging from 1.5 to 3 molecular diameters. On the other hand, significant effect of wall?fluid interaction on the corresponding state coexistence envelope is observed for slit pore confinements ranging from 4 to 40 molecular diameters. Moreover, at a given slit width, shrinking in corresponding state coexistence envelope is observed with increase in the wall?fluid interaction. However, in the larger slit pore width of 30 to 40 molecular diameters, shrinking in the corresponding state vapour–liquid coexistence envelopes become indifferent with the stronger wall?fluid interactions studied in this work. Structural behaviour of coexisting phases in slit pores are also investigated through local density profiles, to understand the overall behaviour of corresponding states coexistence envelopes. Fluctuating positive and negative deviations in the corresponding state spreading pressure with respect to corresponding bulk value is observed for studied wall?fluid interactions and slit pore confinements.     

Evidence for a regular distribution of cholesterol in phospholipid bilayers from diphenylhexatriene fluorescence.

Cholesterol/dimyristoylphosphatidylcholine (DMPC) multilamellar vesicles were studied by steady-state fluorescence using diphenylhexatriene (DPH) as a probe. A series of dips were found in the plot of DPH fluorescence intensity versus cholesterol concentration at certain specific cholesterol concentrations. This observation indicates that there are dominant domains in which cholesterol molecules are regularly distributed on a hexagonal superlattice in the acyl chain matrix of DMPC at critical cholesterol concentrations. These concentrations can be predicted by an equation or a mathematical series, except the one at 33 mol %. These dips of DPH fluorescence intensity are temperature dependent. The excellent agreement between experimental data and calculated values as well as similar previous findings of dips and/or kinks in the excimer-over-monomer fluorescence in pyrenephosphatidylcholine/phospholipid mixtures confirm our conclusion about lateral organizations of cholesterol and acyl lipid chains in cholesterol/phospholipid multilamellar vesicles. The regular distribution model at critical concentration is consistent with the phase diagram of cholesterol/DMPC. Using the model of regular distribution, the physical origin of the liquid-disordered (Ld) phase, liquid-ordered phase (Lo), and coexistence of liquid-disordered phase and Lo phase (Lo + Ld) is discussed on the molecular level.     

On the phase transition in a monolayer adsorbed on graphite at temperatures below the 2D-critical temperature

Monte Carlo simulations in the grand ensemble and meso-canonical ensemble in which the adsorbent is connected to a finite reservoir have been used to study adsorption isotherms for monolayer argon adsorption on graphite at temperatures below the 2D-critical temperature in order to elucidate the microscopic details of the 2D-transitions: vapour–solid, vapour–liquid and liquid–solid. An S-shaped van der Waals (vdW) loop was found when a small square surface was used; however, for large square surfaces and rectangular surfaces the isotherms exhibit a vdW-type loop with a vertical segment which indicates the coexistence of two phases separated by a boundary that changes its shape with the loading. This coexistence occurs at the same chemical potential as determined by the mid-density scheme, developed by Do and co-workers (Z. Liu, L. Herrera, V.T. Nguyen, D.D. Do, and D. Nicholson, A Monte Carlo scheme based on mid-density in a hysteresis loop to determine equilibrium phase transition. Mol Simul. 37(11):932–939, 2011; Z. Liu, D.D. Do, and D. Nicholson, A thermodynamic study of the mid-density scheme to determine the equilibrium phase transition in cylindrical pores. Mol Simul. 38(3):189–199, 2011).     

Quantitation of liquid-crystalline ordering in F-actin solutions.

Actin filaments (F-actin) are important determinants of cellular shape and motility. These functions depend on the collective organization of numerous filaments with respect to both position and orientation in the cytoplasm. Much of the orientational organization arises spontaneously through liquid crystal formation in concentrated F-actin solutions. In studying this phenomenon, we found that solutions of purified F-actin undergo a continuous phase transition, from the isotropic state to a liquid crystalline state, when either the mean filament length or the actin concentration is increased above its respective threshold value. The phase diagram representing the threshold filament lengths and concentrations at which the phase transition occurs is consistent with that predicted by Flory's theory on solutions of noninteracting, rigid cylinders (Flory, 1956b). However, in contrast to other predictions based on this model, we found no evidence for the coexistence of isotropic and anisotropic phases. Furthermore, the phase transition proved to be temperature dependent, which suggests the existence of orientation-dependent interfilament interactions or of a temperature-dependent filament flexibility. We developed a simple method for growing undistorted fluorescent acrylodan-labeled F-actin liquid crystals; and we derived a simple theoretical treatment by which polarization-of-fluorescence measurements could be used to quantitate, for the first time, the degree of spontaneous filament ordering (nematic order parameter) in these F-actin liquid crystals. This order parameter was found to increase monotonically with both filament length and concentration. Actin liquid crystals can readily become distorted by a process known as "texturing." Zigzaging and helicoidal liquid crystalline textures which persisted in the absence of ATP were observed through the polarizing microscope. Possible texturing mechanisms are discussed.