On the Application of Widom's Test Particle Method to Homogeneous and Inhomogeneous Fluids

Abstract

Chemical potentials of a homogeneous and an inhomogeneous Lennard-Jones fluid have been determined by molecular dynamics simulations on the vector computer CYBER 205 by applying essentially the fictitious test particle method of Widom. For the homogeneous fluid we find, contrary to the previous result of Guillot and Guissani, that the simulated chemical potential is independent of the particle number. The crucial point, however, is a sufficiently large cut-off radius in the evaluation of the Boltzmann factor. Comparing with our WCA-type perturbation theory, we get agreement in the chemical potentials within 0.1 kT up to the density n[sgrave]³ = 0.80 and a difference of 0.2 kT at n[sgrave]³ = 0.85. For the inhomogeneous case we consider a fluid in a cylindrical pore and integrate Widom's equation over a certain probe volume as suggested earlier by us. Chemical potentials are then calculated independently in five different probe volumes, which are cylindrical shells. The results agree well from the second to the fourth shell. Inaccuracies in the innermost cylinder can be easily explained by bad statistics. In the shell close to the wall the extremely high local density is responsible for the inaccuracies. Extending the probe volume over all cylindrical shells besides the one closest to the wall is thought to yield rather reliable results for the chemical potential. As a by-product of the simulations we also obtained diffusion coefficients, which are given in an appendix.     

Dissipative particle dynamics of non-spherical particles using a Gaussian density model

The Gaussian density molecular model has been adapted for dissipative particle dynamics. The model, when combined with a soft potential, is shown to be a very flexible mesoscale model exhibiting a wide range of phase behaviour. The soft potential allows relatively large time steps to be used and hence a more rapid equilibration. In addition, the model can be used to study both uniaxial and biaxial systems. We have undertaken a number of pilot studies and have demonstrated that the Gaussian model is able to identify nematic–isotropic phase transitions in liquid crystals and the formation of ordered discotic phases.     

Dissipative particle dynamics simulation on the meso-scale structure of diblock copolymer under cylindrical confinement

The meso-scale structure of symmetric diblock copolymer under cylindrical confinement is studied by dissipative particle dynamics (DPD). The simulation results show that coiled cylindrical geometry is favored in the presence of larger cylinder radius (R/L ₀>～1.5), and the number of rings depends on the cylinder radius. Because of the cylinder wall's selectivity, each block can form the central core, but only the preferential block forms the outmost layer. An approximately linear relationship exists between structure transition point, which is approximately in proportion to the 3/5 exponential of chain length of copolymer and number of layers. As the cylinder radius is decreased, a helical morphology is found. Lamellae parallel to the underside of the cylinder appear when the cylinder radius is made smaller (R/L ₀ < ～1.1).     

Heat transfer simulation using energy conservative dissipative particle dynamics

Dissipative particle dynamics with energy conservation (eDPD) was used to investigate conduction heat transfer in two dimensions under steady-state condition. Various types of boundary condition were implemented to the conduction domain. Besides, 2D conduction with internal heat generation was studied and the heat generation term was used to measure the thermal conductivity and diffusivity of the eDPD system. The boundary conditions used include both the Neumann and Dirichlet boundary conditions. The Neumann boundary condition was applied via adiabatic surfaces and surfaces exposed to convection heat transfer. The DPD simulations were compared to analytical solutions and finite-difference techniques. It was found that DPD appropriately predicts the temperature distribution in the conduction regime. Details of boundary condition implementation and thermal diffusivity measurement are also described in this paper.     

An Extended Gibbs Ensemble

Abstract

A general extended Gibbs ensemble, obtained by augmenting the standard Gibbs ensemble by intermediate states in the spirit of the scaled particle method of Nezbeda and Kolafa [Molec. Simul., 5, 391 (1991)], is introduced. The intermediate states span the states with different number of particles in the simulation boxes and facilitate the transfer of particles even in such complex systems as e.g., mixtures of very different components, systems of flexible polymeric molecules, or systems at very high densities. A general formulation of the ensemble is given and two implementations are considered in detail. The method is very general and is exemplified by studying the fluid-fluid coexistence in a dense binary mixture of the hard-sphere and square-well fluids. It is found that its efficiency is about by factor three greater in comparison with the standard Gibbs ensemble simulations.     

Manufacture of liposomes by isopropanol injection: characterization of the method

Unilamellar liposomes are conventionally prepared by rapid injection of an ethanolic solution of lipids into an aqueous medium. The aim of the present study was to control, more efficiently, vesicle diameter by using an alternative solvent. The results show that isopropanol injection is a good alternative to ethanol injection for the manufacture of liposomes. Particle size can be controlled by the variation of process parameters, such as stirring speed of the aqueous phase and injection flow rate of lipid-isopropanol solution. Diameter of vesicles obtained by this method is less affected by the nature of phospholipid, as well as lipid concentration, than in the ethanol-injection process. In addition, the vesicles are generally smaller (approximately 40–210?nm). Accurate characterization of the particles, by fluorescence, ³¹P-NMR, and cryo–transmission electron microscopy, showed that particles are formed of a single lipid bilayer around an aqueous cavity. We thus provide the scientific community with a fully characterized alternative method to produce unilamellar vesicles.     

Gene Transfer by Means of Lipo- and Polyplexes: Role of Clathrin and Caveolae-Mediated Endocytosis

In this paper we address the contribution of different endocytic pathways to the intracellular uptake and processing of differently sized latex particles and of plasmid DNA complexes by means of fluorescence microscopy and FACS analysis. By using a number of specific inhibitors of either clathrin-dependent or caveolae-dependent endocytosis we were able to discriminate between these two pathways. Latex particles smaller than 200 nm were internalized exclusively by clathrin-mediated endocytosis, whereas larger particles entered the cells via a caveolae-dependent pathway.

The route of uptake of plasmid DNA complexes appears strongly dependent on the nature of the complexes. Thus, lipoplexes containing the cationic lipid DOTAP, were exclusively internalized by a clathrin-dependent mechanism, while polyplexes prepared from the cationic polymer polyethyleneimine (PEI) were internalized in roughly equal proportions by both pathways. Upon incubation of cells with lipoplexes containing the luciferase gene abundant luciferase expression was observed, which was effectively blocked by inhibitors of clathrin-dependent endocytosis but not by inhibitors of the caveolae-dependent uptake mechanism. By contrast, luciferase transfection of the cells with polyplexes was unaffected by inhibition of clathrin-mediated endocytosis, but was nearly completely blocked by inhibitors interfering with the caveolae pathway. The results are discussed with respect to possible differences in the mechanism by which plasmid DNA is released from lipoplexes and polyplexes into the cytosol and to the role of size in the uptake and processing of the complexes. Our data suggest that improvement of non-viral gene transfection could very much benefit from controlling particle size, which would allow targeting of particle internalization via a non-degradative pathway, involving caveolae-mediated endocytosis.     

Ultrahigh frequency lensless ultrasonic transducers for acoustic tweezers application

Similar to optical tweezers, a tightly focused ultrasound microbeam is needed to manipulate microparticles in acoustic tweezers. The development of highly sensitive ultrahigh frequency ultrasonic transducers is crucial for trapping particles or cells with a size of a few microns. As an extra lens would cause excessive attenuation at ultrahigh frequencies, two types of 200‐MHz lensless transducer design were developed as an ultrasound microbeam device for acoustic tweezers application. Lithium niobate single crystal press‐focused (PF) transducer and zinc oxide self‐focused transducer were designed, fabricated and characterized. Tightly focused acoustic beams produced by these transducers were shown to be capable of manipulating single microspheres as small as 5 µm two‐dimensionally within a range of hundreds of micrometers in distilled water. The size of the trapped microspheres is the smallest ever reported in the literature of acoustic PF devices. These results suggest that these lensless ultrahigh frequency ultrasonic transducers are capable of manipulating particles at the cellular level and that acoustic tweezers may be a useful tool to manipulate a single cell or molecule for a wide range of biomedical applications. Biotechnol. Bioeng. 2013; 110: 881–886. © 2012 Wiley Periodicals, Inc.