Understanding the Mechanical Properties of Ultradeformable Liposomes Using Molecular Dynamics Simulations

Improving drug delivery efficiency to solid tumor sites is a central challenge in anticancer therapeutic research. Our previous experimental study (Guo et al., Nat. Commun. 2018, 9, 130) showed that soft, elastic liposomes had increased uptake and accumulation in cancer cells and tumors in vitro and in vivo respectively, relative to rigid particles. As a first step toward understanding how liposomes' molecular structure and composition modulates their elasticity, we performed all-atom and coarse-grained classical molecular dynamics (MD) simulations of lipid bilayers formed by mixing a long-tailed unsaturated phospholipid with a short-tailed saturated lipid with the same headgroup. The former types of phospholipids considered were 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoleoyl-sn-glycero-3-phosphocholine (termed here DPMPC). The shorter saturated lipids examined were 1,2-diheptanoyl-sn-glycero-3-phosphocholine (DHPC), 1,2-didecanoyl-sn-glycero-3-phosphocholine (DDPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Several lipid concentrations and surface tensions were considered. Our results show that DOPC or DPMPC systems having 25-35 mol % of the shortest lipids DHPC or DDPC are the least rigid, having area compressibility moduli KA that are ∼10% smaller than the values observed in pure DOPC or DPMPC bilayers. These results agree with experimental measurements of the stretching modulus and lysis tension in liposomes with the same compositions. These mixed systems also have lower areas per lipid and form more uneven x-y interfaces with water, the tails of both primary and secondary lipids are more disordered, and the terminal methyl groups in the tails of the long lipid DOPC or DPMPC wriggle more in the vertical direction, compared to pure DOPC or DPMPC bilayers or their mixtures with the longer saturated lipid DLPC or DMPC. These observations confirm our hypothesis that adding increasing concentrations of the short unsaturated lipid DHPC or DDPC to DOPC or DPMPC bilayers alters lipid packing and thus makes the resulting liposomes more elastic and less rigid. No formation of lipid nanodomains was noted in our simulations, and no clear trends were observed in the lateral diffusivities of the lipids as the concentration, type of secondary lipid, and surface tension were varied.

Integration of an LPAR1 Antagonist into Liposomes Enhances Their Internalization and Tumor Accumulation in an Animal Model of Human Metastatic Breast Cancer

Lysophosphatidic acid receptor 1 (LPAR1) is elevated in breast cancer. The deregulation of LPAR1, including the function and level of expression, is linked to cancer initiation, progression, and metastasis. LPAR1 antagonists, AM095 or Ki16425, may be effective therapeutic molecules, yet their limited water solubility hinders in vivo delivery. In this study, we report on the synthesis of two liposomal formulations incorporating AM095 or Ki16425, embedded within the lipid bilayer, as targeted nanocarriers for metastatic breast cancer (MBC). The data show that the Ki16425 liposomal formulation exhibited a 50% increase in internalization by MBC mouse epithelial cells (4T1) and a 100% increase in tumor accumulation in a mouse model of MBC compared with that of a blank liposomal formulation (control). At the same time, normal mouse epithelial cells (EpH-4Ev) internalized the Ki16425 liposomal formulation 25% lesser than the control formulation. Molecular dynamics simulations show that the integration of AM095 or Ki16425 modified the physical and mechanical properties of the lipid bilayer, making it more flexible in these liposomal formulations compared with liposomes without drug. The incorporation of an LPAR1 antagonist within a liposomal drug delivery system represents a viable therapeutic approach for targeting the LPA-LPAR1 axis, which may hinder the progression of MBC.

Molecular Dynamics Simulations of Nanostructures Formed by Hydrophobins and Oil in Seawater

Classical molecular dynamics simulations using the Martini coarse-grained force field were performed to study oil nanodroplets surrounded by fungal hydrophobin (HP) proteins in seawater. The class I EAS and the class II HFBII HPs were studied along with two model oils, namely, benzene and n-decane. Both HPs exhibit free energy minima at the oil-seawater interface, which is deeper in benzene compared to the n-decane systems. Larger constraint forces are required to keep both HPs within the n-decane phase compared to inside benzene, with HFBII being more affine to benzene compared to EAS. Smaller surface tensions are observed at benzene-seawater interfaces coated with HPs compared to their n-decane counterparts. In the latter the surface tension remains unchanged upon increases in the concentration of HPs, whereas in benzene systems adding more HPs lead to decreases in surface tension. EAS has a larger tendency to cluster together in the interface compared to HFBII, with both HPs having larger coordination numbers when surrounding benzene droplets compared to when they are around n-decane nanoblobs. The HP-oil nanostructures in seawater examined have radii of gyration ranging between 2 and 12 nm, where the n-decane structures are larger and have more irregular shapes compared to the benzene systems. The n-decane molecules within the nanostructures form a compact spherical core, with the HPs partially covering its surface and clustering together, conferring irregular shapes to the nanostructures. The EAS with n-decane structures are larger and have more irregular shapes compared to their HFBII counterparts. In contrast, in the HP-benzene structures both HPs tend to penetrate the oil part of the droplet. The HFBII-benzene structures having the larger oil/HP ratios examined tend to be more compact and spherical compared to their EAS counterparts; however, some of the HFBII-benzene systems that have smaller oil/HP ratios have a more elongated structure compared to their EAS counterparts. This simulation study provides insights into HP-oil nanostructures that are smaller than the oil droplets and gas bubbles recently studied in experiments and, thus, might be challenging to examine with experimental techniques.

Deep Eutectic Solvents Mixed with Fluorinated Refrigerants for Absorption Refrigeration: A Molecular Simulation Study

Molecular simulations were performed to evaluate mixtures of fluorinated refrigerants with deep eutectic solvents (DESs), for potential use in single-effect absorption refrigeration cycles that use low quality waste heat sources at temperatures of ∼80 °C. The refrigerants considered were the hydrofluorocarbon R245fa and the hydrofluoroolefins R1234zeE and HFO1336mzzE, whereas the DESs evaluated were 1:2 molar mixtures of choline chloride with either ethylene glycol (ethaline) or levulinic acid (levuline) as hydrogen bond donors (HBDs). Assuming the same cycle operating conditions, the waste heat cycle efficiency η was computed for all working fluid mixtures from molecular simulation results of the mixture densities and Henry's law constants of the refrigerants in the DESs, coupled with phase equilibrium calculations and the enthalpies of the pure refrigerants. The largest efficiency was obtained for the mixture R245fa-ethaline (η = 6.82), followed by R245fa-levuline (η = 4.64) and HFO1336mzzE-levuline (η = 2.10). These modest efficiency values could be further increased by tailoring the cycle operating conditions to each particular refrigerant-DES system, as well as optimizing our choice of working fluid mixtures, neither of which we attempted in this study. Strong interactions were observed between the chlorine anions and some of the hydrogen atoms of the refrigerants, but in general the cation-refrigerant and HBD-refrigerant interactions are weaker compared to the refrigerant-refrigerant interactions. Refrigerant molecules have the largest diffusivities and make the cations, anions and HBD to move faster compared to systems of DESs without refrigerant; in general, species in refrigerant-ethaline mixtures have larger diffusivities compared to those for refrigerant-levuline mixtures. We also computed waste heat cycle efficiencies for the same R134a-DES mixtures studied in our previous work, finding significant differences between the efficiencies determined from molecular simulation data and those determined before using the COSMO-RS approach using two standard parametrizations. This observation suggests that further work is needed to improve the accuracy of the COSMO-RS predictions for these systems.

On the Performance of Confined Deep Eutectic Solvents and Ionic Liquids for Separations of Carbon Dioxide from Methane: Molecular Dynamics Simulations

Classical molecular dynamics simulations were used to investigate the performance of slit graphite and titania (rutile) pores of 5.2 nm in width, partially and completely filled with deep eutectic solvents (DESs) or ionic liquids (ILs), for gas separations of a carbon dioxide-methane mixture of 5:95 molar ratio and temperatures and pressures on the order of 318 K and 100 bar, respectively. The DESs studied were ethaline and levuline (1:2 molar mixtures of choline chloride with ethylene glycol or levulinic acid), and the IL considered was 1- n-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [bmim⁺][NTf₂^-]. The performance of these systems in terms of solubility selectivity, diffusion selectivity, and permselectivity was compared against the performance of the bulk solvents (which could also be viewed as a model system for the micrometer-sized pores of a supported IL or DES membrane) and against carbon and rutile pores without preadsorbed solvent. The best performance in terms of permselectivity is obtained for bulk levuline and by rutile pores fully filled by ethaline, followed by graphite pores filled by ethaline and the IL. Empty rutile pores have the largest value of solubility selectivity, followed by bulk ethaline and rutile pores completely filled by the IL. The largest values of diffusivity selectivity were observed for bulk levuline, followed by ethaline completely filling a rutile nanopore and a graphite nanopore completely filled with the IL. These observations are rationalized by examining local density profiles and interaction energies among the different entities in our systems. In general, systems of nanopores fully filled by solvents, as well as the bulk solvents, have larger permselectivities than pores partially filled by the IL or the DESs. Drops of 2-3 orders of magnitude are observed in the gas diffusivity in pores filled with solvents with respect to systems of empty pores, which may be problematic if gas permeation is mainly controlled by diffusion. However, if adsorption dominates the gas permeation within the membrane, our results suggest that systems of levuline in the micrometer-sized pores of a supported DES membrane or ethaline confined in the rutile nanopores of a supported DES phase material might represent promising systems for gas separation.

Computational Evaluation of Mixtures of Hydrofluorocarbons and Deep Eutectic Solvents for Absorption Refrigeration Systems

We used computational tools to evaluate three working fluid mixtures for single-effect absorption refrigeration systems, where the generator (desorber) is powered by waste or solar heat. The mixtures studied here resulted from combining a widely used hydrofluorocarbon (HFC) refrigerant, R134a, with three common deep eutectic solvents (DESs) formed by mixing choline chloride (hydrogen bond acceptor, HBA) with urea, glycerol, or ethylene glycol as the hydrogen bond donor (HBD) species. The COSMOtherm/TmoleX software package was used in combination with refrigerant data from NIST/REFPROP, to perform a thermodynamic evaluation of absorption refrigeration cycles using the proposed working fluid mixtures. Afterward, classical MD simulations of the three mixtures were performed to gain insight on these systems at the molecular level. Larger cycle efficiencies are obtained when R134a is combined with choline chloride and ethylene glycol, followed by the system where glycerol is the HBD, and finally that where the HBD is urea. MD simulations indicate that the local density profiles of all species exhibit very sharp variations in systems containing glycerol or urea; furthermore, the Henry's law constants of R134a in these two systems are larger than those observed for the HFC in choline chloride and ethylene glycol, indicating that R134a is more soluble in the latter DES. Interaction energies indicate that the R134a-R134a interactions are weaker in the system where ethylene glycol is the HBD, as compared to in the other DES. Radial distribution functions confirm that in all systems, the DES species do not form strong directional interactions (e.g., hydrogen bonds) with the R134a molecules. Relatively strong interactions are observed between the Cl anions and the hydrogen atoms in R134a; however, the atom-atom interactions between R134a and the cation and HBD species are weaker and do not play a significant role in the solvation of the refrigerant. In all systems, R134a has the largest diffusion coefficients, followed by the HBD, the anion and the cation; the diffusion coefficients are the largest in the systems containing ethylene glycol, followed by those having glycerol and urea. This work is our first step toward our long-term goal of designing and demonstrating optimal working fluid mixtures for use in absorption refrigeration systems. Our results suggest that COSMO-RS can be used to perform a rapid screening of a large number of working fluid mixtures, and select a few candidates for further exploration using molecular simulations and experiments. These latter approaches can be used to refine the accuracy of the COSMO-RS predictions, and to optimize the selection of optimal working fluid mixtures for demonstration in absorption refrigeration systems powered by solar or waste heat sources.

Molecular simulation of homogeneous nucleation of crystals of an ionic liquid from the melt

The homogeneous nucleation of crystals of the ionic liquid [dmim(+)][Cl(-)] from its supercooled liquid phase in the bulk (P = 1 bar, T = 340 K, representing a supercooling of 58 K) was studied using molecular simulations. The string method in collective variables [Maragliano et al., J. Chem. Phys. 125, 024106 (2006)] was used in combination with Markovian milestoning with Voronoi tessellations [Maragliano et al., J. Chem. Theory Comput. 5, 2589-2594 (2009)] and order parameters for molecular crystals [E. E. Santiso and B. L. Trout, J. Chem. Phys. 134, 064109 (2011)] to sketch a minimum free energy path connecting the supercooled liquid and the monoclinic crystal phases, and to determine the free energy and the rates involved in the homogeneous nucleation process. The physical significance of the configurations found along this minimum free energy path is discussed with the help of calculations based on classical nucleation theory and with additional simulation results obtained for a larger system. Our results indicate that, at a supercooling of 58 K, the liquid has to overcome a free energy barrier of the order of 60 kcal/mol and to form a critical nucleus with an average size of about 3.6 nm, before it reaches the thermodynamically stable crystal phase. A simulated homogeneous nucleation rate of 5.0 × 10(10) cm(-3) s(-1) was obtained for our system, which is in reasonable agreement with experimental and simulation rates for homogeneous nucleation of ice at similar degrees of supercooling. This study represents our first step in a series of studies aimed at understanding the nucleation and growth of crystals of organic salts near surfaces and inside nanopores.

Bubble bursting as an aerosol generation mechanism during an oil spill in the deep-sea environment: molecular dynamics simulations of oil alkanes and dispersants in atmospheric air/salt water interfaces

Potential of mean force (PMF) calculations and molecular dynamics (MD) simulations were performed to investigate the properties of oil n-alkanes [i.e., n-pentadecane (C15), n-icosane (C20) and n-triacontane (C30)], as well as several surfactant species [i.e., the standard anionic surfactant sodium dodecyl sulfate (SDS), and three model dispersants similar to the Tween and Span species present in Corexit 9500A] at air/salt water interfaces. This study was motivated by the 2010 Deepwater Horizon (DWH) oil spill, and our simulation results show that, from the thermodynamic point of view, the n-alkanes and the model dispersants have a strong preference to remain at the air/salt water interface, as indicated by the presence of deep free energy minima at these interfaces. The free energy minimum of these n-alkanes becomes deeper as their chain length increases, and as the concentration of surfactant species at the interface increases. The n-alkanes tend to adopt a flat orientation and form aggregates at the bare air/salt water interface. When this interface is coated with surfactants, the n-alkanes tend to adopt more tilted orientations with respect to the vector normal to the interface. These simulation results are consistent with the experimental findings reported in the accompanying paper [Ehrenhauser et al., Environ. Sci.: Processes Impacts 2013, in press, (DOI: 10.1039/c3em00390f)]. The fact that these long-chain n-alkanes show a strong thermodynamic preference to remain at the air/salt water interfaces, especially if these interfaces are coated with surfactants, makes these species very likely to adsorb at the surface of bubbles or droplets and be ejected to the atmosphere by sea surface processes such as whitecaps (breaking waves) and bubble bursting. Finally, the experimental finding that more oil hydrocarbons are ejected when Corexit 9500A is present in the system is consistent with the deeper free energy minima observed for the n-alkanes at the air/salt water interface at increasing concentrations of surfactant species.

Bubble bursting as an aerosol generation mechanism during an oil spill in the deep-sea environment: laboratory experimental demonstration of the transport pathway

Oil spills in the deep-sea environment such as the 2010 Deep Water Horizon oil spill in the Gulf of Mexico release vast quantities of crude oil into the sea-surface environment. Various investigators have discussed the marine transport and fate of the oil into different environmental compartments (air, water, sediment, and biota). The transport of the oil into the atmosphere in these previous investigations has been limited to only evaporation, a volatility dependent pathway. In this work, we studied the aerosolization of oil spill matter via bursting bubbles as they occur during whitecaps in a laboratory aerosolization reactor. By evaluating the alkane content in oil mousse, crude oil, the gas phase, and particulate matter we clearly demonstrate that aerosolization via bursting bubbles is a solubility and volatility independent transport pathway for alkanes. The signature of alkane fractions in the native oil and aerosolized matter matched well especially for the less volatile alkanes (C20-C29). Scanning electron microscope interfaced with energy dispersive X-ray images identified the carbon fractions associated with salt particles of aerosols. Theoretical molecular dynamics simulations in the accompanying paper lend support to the observed propensity for alkanes at air-salt water interfaces of breaking bubbles and the produced droplets. The presence of a dispersant in the aqueous phase increased the oil ejection rate at the surface especially for the C20-C29 alkanes. The information presented here emphasizes the need to further study sea-spray aerosols as a possible transport vector for spilled oil in the sea surface environment.

Molecular modeling of the green leaf volatile methyl salicylate on atmospheric air/water interfaces

Methyl salicylate (MeSA) is a green leaf volatile (GLV) compound that is emitted in significant amounts by plants, especially when they are under stress conditions. GLVs can then undergo chemical reactions with atmospheric oxidants, yielding compounds that contribute to the formation of secondary organic aerosols (SOAs). We investigated the adsorption of MeSA on atmospheric air/water interfaces at 298 K using thermodynamic integration (TI), potential of mean force (PMF) calculations, and classical molecular dynamics (MD) simulations. Our molecular models can reproduce experimental results of the 1-octanol/water partition coefficient of MeSA. A deep free energy minimum was found for MeSA at the air/water interface, which is mainly driven by energetic interactions between MeSA and water. At the interface, the oxygenated groups in MeSA tend to point toward the water side of the interface, with the aromatic group of MeSA lying farther away from water. Increases in the concentrations of MeSA lead to reductions in the height of the peaks in the MeSA-MeSA g(r) functions, a slowing down of the dynamics of both MeSA and water at the interface, and a reduction in the interfacial surface tension. Our results indicate that MeSA has a strong thermodynamic preference to remain at the air/water interface, and thus chemical reactions with atmospheric oxidants are more likely to take place at this interface, rather than in the water phase of atmospheric water droplets or in the gas phase.

Nontemplated approach to tuning the spectral properties of cyanine-based fluorescent nanoGUMBOS

Template-free controlled aggregation and spectral properties in fluorescent organic nanoparticles (FONs) is highly desirable for various applications. Herein, we report a nontemplated method for controlling the aggregation in near-infrared (NIR) cyanine-based nanoparticles derived from a group of uniform materials based on organic salts (GUMBOS). Cationic heptamethine cyanine dye 1,1',3,3,3',3'-hexamethylindotricarbocyanine (HMT) was coupled with five different anions, viz., [NTf(2)(-)], [BETI(-)], [TFPB(-)], [AOT(-)], and [TFP4B(-)], by an ion-exchange method to obtain the respective GUMBOS. The nanoGUMBOS obtained via a reprecipitation method were primarily amorphous and spherical (30-100 nm) as suggested by selected area electron diffraction (SAED) and transmission electron microscopy (TEM). The formation of tunable self-assemblies within the nanoGUMBOS was characterized using absorption and fluorescence spectroscopy in conjunction with molecular dynamics simulations. Counterion-controlled spectral properties observed in the nanoGUMBOS were attributed to variations in J/H ratios with different anions. Association with the [AOT(-)] anion afforded predominant J aggregation enabling the highest fluorescence intensity, whereas [TFP4B(-)] disabled the fluorescence due to predominant H aggregation in the nanoparticles. Analyses of the stacking angle of the cations based on molecular dynamic simulation results in [HMT][NTf(2)], [HMT][BETI], and [HMT][AOT] dispersed in water and a visual analysis of the representative simulation snapshots also imply that the type of aggregation was controlled through the counterion associated with the dye cation.

Modeling micelle-templated mesoporous material SBA-15: atomistic model and gas adsorption studies

We report the development of a realistic molecular model for mesoporous silica SBA-15, which includes both the large cylindrical mesopores and the smaller micropores in the pore walls. The methodology for modeling the SBA-15 structure involves molecular and mesoscale simulations combined with geometrical interpolation techniques. First, a mesoscale model is prepared by mimicking the synthesis process using lattice Monte Carlo simulations. The main physical features of this mesoscale pore model are then carved out of an atomistic silica block; both the mesopores and the micropores are incorporated from the mimetic simulations. The calculated pore size distribution, surface area, and simulated TEM images of the model structure are in good agreement with those obtained from experimental samples of SBA-15. We then investigate the adsorption of argon in this structure using Grand Canonical Monte Carlo (GCMC) simulations. The adsorption results for our SBA-15 model are compared with those for a similar model that does not include the micropores; we also compare with results obtained in a regular cylindrical pore. The simulated adsorption isotherm for the SBA-15 model shows semiquantitative agreement with the experimental isotherm for a SBA-15 sample having a similar pore size. We observe that the presence of the micropores leads to increased adsorption at low pressure compared to the case of a model without micropores in the pore walls. At higher pressures, for all models, the filling proceeds via the monolayer-multilayer adsorption on the mesopore surface followed by capillary condensation, which is mainly controlled by the mesopore diameter and is not influenced by the presence of the micropores.

Quadrupolar particles in a nematic liquid crystal: effects of particle size and shape

We investigate the effects of particle size and shape on the quadrupolar (Saturn-ring-like) defect structures formed by a nematic liquid crystal around nm-sized and mum -sized particles with spherical and spherocylindrical shapes. We also report results for the potentials of mean force in our systems, calculated using a mesoscale theory for the tensor order parameter Q of the nematic. Our results indicate that for pairs of nm-sized particles in close proximity, the nematic forms "entangled hyperbolic" defect structures regardless of the shape of the nanoparticles. In our calculations with nanoparticles we did not observe any other entangled or unentangled defect structures, in contrast to what was reported for pairs of mum -sized spherical particles. Such a finding suggests that the "entangled hyperbolic" defect structures are the most stable for pairs of nanoparticles in close proximity. For pairs of mum -sized particles, our results indicate that the nematic forms entangled "figure-of-eight" defect structures around pairs of spheres and spherocylinders. Our results suggest that the transition between "entangled hyperbolic" and figure-of-eight defect structures takes place when the diameter of the particle is between D=100 nm and 1 microm . We have also calculated the torques that develop when pairs of spherocylindrical nanoparticles in a nematic approach each other. Our calculations suggest that the nematic-mediated interactions between the nm-sized particles are fairly strong, up to 5700 k{B}T for the case of pairs of spherocylindrical nanoparticles arranged with their long axis parallel to each other. Furthermore, these interactions can make the particles to bind together at specific locations, and thus could be used to assemble the particles into ordered structures with different morphologies.

Nanoparticles in nematic liquid crystals: Interactions with nanochannels

A mesoscale theory for the tensor order parameter Q is used to investigate the structures that arise when spherical nanoparticles are suspended in confined nematic liquid crystals (NLCs). The NLC is "sandwiched" between a wall and a small channel. The potential of mean force is determined between particles and the bottom of the channels or between several particles. Our results suggest that strong NLC-mediated interactions between the particles and the sidewalls of the channels, on the order of hundreds of k(B)T, arise when the colloids are inside the channels. The magnitude of the channel-particle interactions is dictated by a combination of two factors, namely, the type of defect structures that develop when a nanoparticle is inside a channel, and the degree of ordering of the nematic in the region between the colloid and the nanochannel. The channel-particle interactions become stronger as the nanoparticle diameter becomes commensurate with the nanochannel width. Nanochannel geometry also affects the channel-particle interactions. Among the different geometries considered, a cylindrical channel seems to provide the strongest interactions. Our calculations suggest that small variations in geometry, such as removing the sharp edges of the channels, can lead to important reductions in channel-particle interactions. Our calculations for systems of several nanoparticles indicate that linear arrays of colloids with Saturn ring defects, which for some physical conditions are not stable in a bulk system, can be stabilized inside the nanochannels. These results suggest that nanochannels and NLCs could be used to direct the assembly of nanoparticles into ordered arrays with unusual morphologies.

Anisotropic nanoparticles immersed in a nematic liquid crystal: defect structures and potentials of mean force

We report results for the potential of mean force (PMF) and the defect structures that arise when spherocylindrical nanoparticles are immersed in a nematic liquid crystal. Using a dynamic field theory for the tensor order parameter Q of the liquid crystal, we analyzed configurations, including one, two, and three elongated particles, with strong homeotropic anchoring at their surfaces. For systems with one nanoparticle, the most stable configuration is achieved when the spherocylinder is placed with its long axis perpendicular to the far-field director, for which the defect structure consists of an elongated Saturn ring. For systems with two or three nanoparticles with their long axes placed perpendicular to the far-field director, at small separations the defect structures consist of incomplete Saturn rings fused with new disclination rings orthogonal to the original ones, in analogy to results previously observed for spherical nanoparticles. The shape of these orthogonal rings depends on the nanoparticles' configuration, i.e., triangular, linear, or parallel with respect to their long axis. A comparison of the PMFs indicates that the latter configuration is the most stable. The stability of the different arrays depends on whether orthogonal disclination rings form or not, their size, and the curvature effects in the interparticle regions. Our results suggest that the one-elastic-constant approximation is valid for the considered systems; similar results were obtained when a three-constant expression is used to represent the elastic free energy. The attractive interactions between the elongated particles were compared to those observed for spheres of similar diameters. Similar interparticle energies were observed for linear arrays; in contrast, parallel and triangular arrays of spherocylinders yielded interactions that were up to 3.4 times stronger than those observed for spherical particles.

Adsorption of simple gases in MCM-41 materials: the role of surface roughness

This paper reports the development and testing of atomistic models of silica MCM-41 pores. Model A is a regular cylindrical pore having a constant section. Model B has a surface disorder that reproduces the morphological features of a pore obtained from an on-lattice simulation that mimics the synthesis process of MCM-41 materials. Both models are generated using a similar procedure, which consists of carving the pore out of an atomistic silica block. The differences between the two models are analyzed in terms of small angle neutron scattering spectra as well as adsorption isotherms and isosteric heat curves for Ar at 87 K and Xe at 195 K. As expected for capillary condensation in regular nanopores, the Ar and Xe adsorption/desorption cycles for model A exhibit a large hysteresis loop having a symmetrical shape, i.e., with parallel adsorption and desorption branches. The features of the adsorption isotherms for model B strongly depart from those observed for model A. Both the Ar and Xe adsorption branches for model B correspond to a quasicontinuous pore filling that involves coexistence within the pore of liquid bridges and gas nanobubbles. As in the case of model A, the Ar adsorption isotherm for model B exhibits a significant hysteresis loop; however, the shape of the loop is asymmetrical with a desorption branch much steeper than the adsorption branch. In contrast, the adsorption/desorption cycle for Xe in model B is quasicontinuous and quasireversible. Comparison with adsorption and neutron scattering experiments suggests that model B is too rough at the molecular scale but reproduces reasonably the surface disorder of real MCM-41 at larger length scales. In contrast, model A is smooth at small length scales in agreement with experiments but seems to be too ordered at larger length scales.

Molecular modeling of freezing of simple fluids confined within carbon nanotubes

We report Monte Carlo simulation results for freezing of Lennard-Jones carbon tetrachloride confined within model multiwalled carbon nanotubes of different diameters. The structure and thermodynamic stability of the confined phases, as well as the transition temperatures, were determined from parallel tempering grand canonical Monte Carlo simulations and free-energy calculations. The simulations show that the adsorbate forms concentric molecular layers that solidify into defective quasi-two-dimensional hexagonal crystals. Freezing in such concentric layers occurs via intermediate phases that show remnants of hexatic behavior, similar to the freezing mechanism observed for slit pores in previous works. The adsorbate molecules in the inner regions of the pore also exhibit changes in their properties upon reduction of temperature. The structural changes in the different regions of adsorbate occur at temperatures above or below the bulk freezing point, depending on pore diameter and distance of the adsorbate molecules from the pore wall. The simulations show evidence of a rich phase behavior in confinement; a number of phases, some of them inhomogeneous, were observed for the pore sizes considered. The multiple transition temperatures obtained from the simulations were found to be in good agreement with recent dielectric relaxation spectroscopy experiments for CCl(4) confined within multiwalled carbon nanotubes.

An experimental study of melting of CCl4 in carbon nanotubes

We report dielectric relaxation spectroscopy measurements of the melting point of carbon tetrachloride confined within open-tip multi-walled carbon nanotubes with two different pore diameters, 4.0 and 2.8 nm. In both cases, a single transition temperature well above the bulk melting point was obtained for confined CCl4. These results contrast with what was obtained in our previous measurements using carbon nanotubes with a pore diameter of 5.0 nm, where multiple transition temperatures both above and below the bulk melting point of CCl4 were observed. Our experimental measurements are consistent with our recent molecular simulation results (F. R. Hung, B. Coasne, E. E. Santiso, K. E. Gubbins, F. R. Siperstein and M. Sliwinska-Bartkowiak, J. Chem. Phys., 2005, 122, 144706). Although the simulations overestimate the temperatures in which melting upon confinement occurs, both simulations and experiments suggest that all regions of adsorbate freeze at the same temperature, and that freezing occurs at higher temperatures upon reduction of the pore diameter.