The effect of alumina nanoparticles on the thermal properties of PMMA: a molecular dynamics simulation

Metal oxides, as one of the most promising flame retardant additives, improve the fire retardant and the thermal stability properties of polymers. In the present study, molecular dynamics (MD) simulations based on the united atom model were applied to study the effect of alumina nanoparticles on the density, thermal conductivity, heat capacity, and thermal diffusivity of isotactic poly(methyl methacrylate) (is-PMMA). Thermal diffusivity of PMMA and PMMA/alumina nanocomposite were investigated through calculating thermal conductivity, density and heat capacity in the range of 300–700?K. Heat capacity can be calculated using fluctuations properties of energy. Thermal conductivity was calculated through the nonequilibrium molecular dynamics (NEMD) simulation by Fourier’s law approach. Our results show that the addition of alumina nanoparticles decreases the heat capacity and increases the glass transition temperature (T_g), thermal conductivity and thermal diffusivity of the PMMA. Therefore, the addition of alumina nanoparticles to PMMA improves the fire retardancy of the polymer. In addition, we illustrate the links between the intermolecular and bulk properties of PMMA in the presence of the alumina nanoparticles.     

Differences in the adsorption and diffusion behaviour of water and non-polar gases in nanoporous carbon: role of cooperative effects of pore confinement and hydrogen bonding

We investigate the effect of pore confinement and molecular geometry on the adsorption and self-diffusion of H₂O, CO₂, Ar, CH₄, C₃H₆, SF₆ and C₅H₁₂, in a realistic model of nanoporous silicon carbide derived carbon (SiC-DC), constructed using hybrid reverse Monte Carlo simulation. Adsorption isotherms, adsorbate–adsorbate and adsorbate–adsorbent contributions to the isosteric heat of adsorption are determined to study the effect of pore confinement, microporosity and molecular geometry on adsorption of these gases. We describe the cooperative effect of pore confinement and hydrogen bonding on the formation of water clusters and anomalous adsorption behaviour of water compared with non-polar gases. We find that, in contrast to literature results based on the slit-pore model, pore-filling does not occur below the saturation pressure in hydrophobic amorphous carbon materials such as SiC-DC and activated carbon fibre. We also compare self-diffusivities and activation energy barriers of water and non-polar gases in the microporous structure of SiC-DC to identify underlying correlations with molecular properties. We demonstrate that the self-diffusivity of water deviates considerably from the correlation between diffusivity and molecular kinetic diameter observed for non-polar gases. This is attributed to the reduced diffusivity of water, and its relatively large energy barrier at high loadings despite its small kinetic diameter, which is due to the blocking effect of water clusters at pore entries.     

N-octane diffusivity enhancement via carbon dioxide in silica slit-shaped nanopores – a molecular dynamics simulation

Equilibrium molecular dynamics simulations were conducted to study the competitive adsorption and diffusion of mixtures containing n-octane and carbon dioxide confined in slit-shaped silica pores of width 1.9 nm. Atomic density profiles substantiate strong interactions between CO₂ molecules and the protonated pore walls. Non-monotonic change in n-octane self-diffusion coefficients as a function of CO₂ loading was observed. CO₂ preferential adsorption to the pore surface is likely to attenuate the surface adsorption of n-octane, lower the activation energy for n-octane diffusivity, and consequently enhance n-octane mobility at low CO₂ loading. This observation was confirmed by conducting test simulations for pure n-octane confined in narrower pores. At high CO₂ loading, n-octane diffusivity is hindered by molecular crowding. Thus, n-octane diffusivity displays a maximum. In contrast, within the concentration range considered here, the self-diffusion coefficient predicted for CO₂ exhibits a monotonic increase with loading, which is attributed to a combination of effects including the saturation of the adsorption capacity of the silica surface. Test simulations suggest that the results are strongly dependent on the pore morphology, and in particular on the presence of edges that can preferentially adsorb CO₂ molecules and therefore affect the distribution of these molecules equally on the pore surface, which appears to be required to provide the effective enhancement of n-octane diffusivity.     

Molecular dynamics simulations of a bacterial autotransporter: NalP from Neisseria meningitidis

NalP is an autotransporter secretory protein found in the outer membrane of Neisseria meningitidis. The crystal structure of the NalP translocator domain revealed a transmembrane β-barrel containing a central α-helix. The role of this α-helix, and of the conformational dynamics of the β-barrel pore have been studied via atomistic molecular dynamics simulations. Three simulations, each of 10 ns duration, of NalP embedded within a solvated DMPC bilayer were performed. The helix was removed from the barrel interior in one simulation. The conformational stability of the protein is similar to that of other outer membrane proteins, e.g., OmpA, in comparable simulations. The transmembrane β-barrel is stable even in the absence of the α-helix. Removal of the helix results in an influx of water into the pore region, suggesting the helix acts as a ‘plug’. Water molecules entering the resultant pore form hydrogen bonds with the barrel lining that compensate for the loss of helix-barrel hydrogen bonds. The dimensions of the pore fluctuate over the course of the simulation revealing it to be flexible, but only wide enough to allow transport of the passenger domain in an unfolded or extended conformation. The simulations help us to understand the role of the central helix in plugging the pore and in maintaining the width of the barrel, and show that the NalP monomer is sufficient for the transport of the passenger domain in an unfolded or extended conformation.     

Dimensionality reduction of free energy profiles of benzene in silicalite-1: calculation of diffusion coefficients using transition state theory

This study describes the development and application of a new computational methodology for calculating the self-diffusivity of sorbate molecules strongly confined within shape-selective nanoporous materials. An umbrella sampling strategy, employing repulsive walls to confine the sorbate within specific regions of the pore space, is invoked to extract free energy profiles with respect to the sorbate degrees of freedom. Based on these profiles, it is shown how the multidimensional problem of translational diffusion of benzene in flexible silicalite-1 can be reduced first to a six-dimensional problem, then to a three-dimensional (3D) problem and finally, to a 1D problem. A 3D free energy distribution is accumulated as a function of the benzene centre of mass position and ultimately reduced to a set of 1D profiles for the benzene centre of mass along the pore axes. From these profiles, the rate constants for jumps executed by benzene between sorption sites are calculated using transition state theory; from the latter rate constants, the low-occupancy self-diffusivity is obtained using the MESoRReD method [Kolokathis PD, Theodorou DN. On solving the master equation in spatially periodic systems. J Chem Phys. 2012;137:034112]. The activation energy for diffusion and preferred orientations in the various sorption states in silicalite are in very favourable agreement with available experimental measurements.     

Computer Simulation of Isothermal Mass Transport in Graphite Slit Pores

Abstract

Results are presented from a simulation study of the mass transport of oxygen and nitrogen through graphite slit pores. The work is motivated by an attempt to understand the molecular origins of the kinetic selectivity displayed when air is separated into its major components using pressure swing adsorption. A combination of non-equilibrium molecular dynamics (NEMD), equilibrium molecular dynamics (EMD) and grand canonical Monte Carlo methods has been employed in our study to extract the maximum information. Transport diffusivities, self-diffusivities, permeabilities and Darken thermodynamic factors have been calculated as a function of pore width and temperature for pure component oxygen and nitrogen. In addition, new EMD simulation data for an 80:20 mixture of nitrogen and oxygen is reported, including a direct calculation of the Stefan-Maxwell coefficients. The results are discussed in terms of the oxygen selectivity and the possible mechanisms, which increase or decrease this quantity.

We find that the pore width behaviour of the diffusion coefficients consists of three distinct regimes: a regime at larger pore widths in which single component diffusion coefficients are largely independent of pore width, an optimum pore width at which both diffusivities increase substantially but the slit pore is selective towards nitrogen, and a regime at very low pore widths at which the diffusivities decrease sharply, but the slits are selective towards oxygen. The mechanism behind each of these regimes is discussed in terms of “entropic” effects and potential barrier heights.

We have also found that permeability selectivity is substantially reduced in a mixture of the two gases with a composition similar to that of air. Cross diffusion coefficients in the mixture have been calculated and shown to be non-negligible.     

Self diffusion of argon in flexible,single wall,carbon nanotubes

The high-throughput Condor environment now allows many simulations to be performed on related systems, whether the focus is on improving the statistics or on broadening the range of conditions under which these simulations run. We illustrate the scope of the approach by using equilibrium molecular dynamics (EMD) to calculate self-diffusivities of argon atoms diffusing through single wall carbon nanotubes (SWNT). The diameters of the tubes and their helicities were varied and different argon loadings were studied. We also considered the effect of the rigidity/flexibility of the tube on the diffusivity. We found that the helicity and flexibility of the tubes have almost no noticeable influences. The size of the pore had a small effect, but the diffusivity depended essentially on the fluid loading.     

Adsorption and diffusion of CO2 and CH4 in covalent organic frameworks: an MC/MD simulation study

Covalent organic frameworks (COFs) are a promising gas separation material which have been developed recently. In this work, we have used grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) simulations to investigate the adsorption and diffusion properties of CO₂ and CH₄ in five recent synthesised COF materials. We have also considered the properties of amino-modified COFs by adding –NH₂ group to the five COFs. The adsorption isotherm, adsorption/diffusion selectivity, self/transport diffusion coefficients have been examined and discussed. All of the five COFs exhibit promising adsorption selectivity which is higher than common nanoporous materials. An S-shaped adsorption isotherm can be found for CO₂ instead of CH₄ adsorption. The introduction of –NH₂ group is effective at low pressure region (<200?kPa). The diffusion coefficients are similar for TS-COFs but increase with the pore size for PI-COFs, and the diffusion coefficients seem less dependent on the –NH₂ groups.     

Transport diffusivities of fluids in nanopores by non-equilibrium molecular dynamics simulation

We present a method to study fluid transport through nanoporous materials using highly efficient non-equilibrium molecular dynamics simulations. A steady flow is induced by applying an external field to the fluid particles within a small slab of the simulation cell. The external field generates a density gradient between both sides of the porous material, which in turn triggers a convective flux through the porous medium. The heat dissipated by the fluid flow is released by a Gaussian thermostat applied to the wall particles. This method is effective for studying diffusivities in a slit pore as well as more natural, complex wall geometries. The dependence of the diffusive flux on the external field sheds light on the transport diffusivities and allows a direct calculation of effective diffusivities. Both pore and fluid particle interactions are represented by coarse-grained molecular models in order to present a proof-of-concept and to retain computational efficiency in the simulations. The application of the method is demonstrated in two different scenarios, namely the effective mass transport through a slit pore and the calculation of the effective self-diffusion through this system. The method allows for a distinction between diffusive and convective contributions of the mass transport.     

Intestinal transport of TRH analogs through PepT1: the role of in silico and in vitro modeling

The present study involves molecular docking, molecular dynamics (MD) simulation studies, and Caco‐2 cell monolayer permeability assay to investigate the effect of structural modifications on PepT1‐mediated transport of thyrotropin releasing hormone (TRH) analogs. Molecular docking of four TRH analogs was performed using a homology model of human PepT1 followed by subsequent MD simulation studies. Caco‐2 cell monolayer permeability studies of four TRH analogs were performed at apical to basolateral and basolateral to apical directions. Inhibition experiments were carried out using Gly‐Sar, a typical PepT1 substrate, to confirm the PepT1‐mediated transport mechanism of TRH analogs. P_app of the four analogs follows the order: NP‐1894 < NP‐2378 < NP‐1896 < NP‐1895. Higher absorptive transport was observed in the case of TRH analogs, indicating the possibility of a carrier‐mediated transport mechanism. Further, the significant inhibition of the uptake of Gly‐Sar by TRH analogs confirmed the PepT1‐mediated transport mechanism. Glide docking scores of all the four analogues were in good agreement with their transport rates, suggesting the role of substrate binding affinity in the PepT1‐mediated transport of TRH analogs. MD simulation studies revealed that the polar interactions with amino acid residues present in the active site are primarily responsible for substrate binding, and a downward trend was observed with the increase in bulkiness at the N‐histidyl moiety of TRH analogs. Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling of Platinum Clusters in H-Mordenite

The size, location and structure of Pt clusters in H-mordenite have been investigated by molecular mechanics energy minimization and molecular dynamics simulation techniques using the Catalysis software of Molecular Simulations (MSI). Lattice energy minimizations are performed to study the effects of the specific framework aluminum positions on the location and stability of monoatomic Pt sites in H-mordenite. The lattice energies relative to the siliceous platinum-aluminosilicate structure reveal that the stability of a single Pt atom in H-mordenite is remarkably influenced by the specific location of the Al atoms in the lattice. At the studied Si/Al ratio of two Al ions per unit cell, a stabilization of the H-mordenite lattice upon Pt deposition is obtained. Moreover, lattice energy calculations on Pt/aluminosilicate mordenites of different metal contents per unit cell have been performed. An optimum size for the aggregate confined to the 12-ring main channel that is almost independent of the Pt content per mordenite unit cell has been found. The structural features of the resulting clusters at the end of molecular dynamics simulations on Pt/alumina-mordenites reflect a strong metal-zeolite interaction. The present results are consistent with a previous molecular dynamics simulation study on the structure of platinum deposited on SiO₂ surfaces.     

Combined Diffusive and Viscous Transport of Methane in a Carbon Slit Pore

Abstract

The transport of mass through porous materials can occur by essentially two different mechanisms: (1) diffusion and (2) viscous flow. The former occurs when there is a gradient in chemical potential of the pore fluid, while the latter occurs in the presence of a pressure gradient. In general, fluid transport occurs by both of these mechanisms and their respective contributions to the total intra-pore flux are approximately additive. Experimentally, there is no unambiguous way of determining the individual contributions to the total flux of these two modes of transport. Fortunately, molecular simulations does provide a solution.

We present a novel simulation method in which the separate contributions to the total flux are determined. The method involves the use of two non-equilibrium molecular dynamics techniques: dual control volume grand canonical molecular dynamics (DCV GCMD) and an algorithm for simulating planar Poiseuille flow. We apply this technique to study the combined (viscous and diffusive) transport of methane through single slit-shaped graphite pores of width 2.5, 5.0 and 10.0 methane diameters. We find that the viscous contribution to the total intrapore flux through each of these pores is 10%, 15% and 34%, respectively.     

Comparative simulations of aquaporin family: AQP1, AQPZ, AQP0 and GlpF

Molecular dynamics simulations were performed for four members of the aquaporin family (AQP1, AQPZ, AQP0, and GlpF) in the explicit membrane environment. The single-channel water permeability, pf, was evaluated to be GlpF approximately AQPZ > AQP1 > AQP0, while their relative pore sizes were GlpF > AQP1 > AQPZ > AQP0. This relation between pf and pore size indicates that water permeability was determined not only by the channel radius, but also another competing factor. Analysis of water dynamics revealed that this factor was the single-file nature of water transport.     

The influences of protonation state of histidine on aromatic/arginine region of aquaporin-1 protein

Aquaporin-1 (AQP1) is widely distributed in the epithelial tissue for water absorption and secretion. The histidine (His182) residue in the aromatic/arginine (ar/R) constriction region plays an important role in transporting water through the membrane. In this study, we have performed a total of 46 ns equilibrium molecular dynamics (MD) simulations, and obtained the influence of His182 in two protonation states (Hsd is the proton at N_δ and Hse is the proton at N_?) on the ar/R region. Water permeation rate shows that it is easier for water molecules to permeate the ar/R region of the AQP1 with residue in the Hsd state than in the Hse state. The minimum radii of the pore in the ar/R region were calculated during the last 10 ns MD simulation. We have analysed the correlation among the state of the pore (open or close), the minimum radius of the ar/R region and the dihedral angles < C_β-C_γ-C_δ-N_? of Arg197. The results show that the minimum radius can be used to mark the state of the pore.     

Amino acid modified carbon nanotubes with optimal pore size for chiral separation

Precisely controlling pore size of porous materials is of great importance for chiral separation, but a great challenge in practical applications. In contrast, the molecular dynamics (MD) simulation can be quite a convenient way to determine the effect of the pore dimension on the chiral resolution performances and thus to define the optimal pore size. In this work, inner-wall functionalised carbon nanotubes (CNTs) were used as porous materials and D- and L-phenylalanine were selected as chiral probes. The enantioseparation behaviour was investigated via varying the pore diameter of CNTs, controlling the grafting amount of chiral selectors and tuning the spacer length. Results show that varying the pore size has a significant effect on the enantioselectivity. Additionally, the effect of the introduction of varying the grafting ratio and tuning the spacer length on the chiral separation performance was also examined in this work. It was found that varying the grafting ratio, especially the spacer length between substrates and selectors, could also be one of the most effective alternatives to improving enantioselectivity. Our findings can provide a guidance for the practical applications in the chiral separation.     

On the Approximation of Solvent Effects on the Conformation and Dynamics of Cyclosporin A by Stochastic Dynamics Simulation Techniques

Abstract

The molecular simulation technique of stochastic dynamics (SD) is tested by application to the immunosuppressive drug cyclosporin A (CPA). Two stochastic dynamics simulations are performed, one (SD_CCl4 ) with atomic friction coefficients proportional to the viscosity of the nonpolar solvent CCl₄, and one (SD_H2O) with atomic friction coefficients corresponding to an aqueous solution. The atomic friction coefficients are also taken proportional to an approximate expression for the atomic accessible surface area. The properties of both stochastic dynamics simulations are compared to those of two full molecular dynamics (MD) simulations of cyclosporin A, one in a box with 591 CCl₄ molecules, and one in a box with 632 H₂O molecules.

The properties of cyclosporin A as found in the molecular dynamics simulation in CCl₄ are well reproduced by the SD_CCl4 simulation. This indicates that the neglect of a mean force reresenting the average solvent effects on the solute is justified in the case of nonpolar solvents. For polar solvents, like water, this mean force may not be neglected. The SD_H2O simulation of cyclosporin A clearly fails to reproduce the amount of hydrogen bonding found in the molecular dynamics stimulation of cyclosporin A in water.

A comparison with a molecular dynamics simulation of cyclosporin A in vacuo shows that both the SD_CCl4 and the SD_H2O simulation come closer to the properties of the molecular dynamics simulations in CCl₄ and in H₂O than a molecular dynamics simulation in vacuo.     

Molecular dynamics simulation of O2 diffusion in polydimethylsiloxane (PDMS) and end-linked PDMS networks

Molecular dynamics simulations are used to compute diffusion coefficients for O₂ molecules in polydimethylsiloxane (PDMS) and end-linked PDMS networks. The PDMS chains and penetrants are modelled using a hybrid interatomic potential which treats the Si and O atoms along the chain backbone explicitly while coarse-graining the methyl side groups and penetrants. In PDMS models with different molecular weights, diffusivity of the O₂ penetrants is found to modestly decrease with an increase in chain length. To match typical experimental conditions, the end-linked PDMS networks are constructed with a PDMS to crosslinking (CL) molecule mass ratio of 5:1 or 10:1, demanding that the number of CL molecules exceeds the number of PDMS chains in each model. Despite end-linking, the presence of non-bonded CL molecules promotes increased O₂ diffusivity in comparison with uncrosslinked PDMS. Temperature dependence is captured using the Williams–Landel–Ferry equation.     

Diffusion of sucrose and dextran through agar gel membranes

Mass transfer limitations severely impede the performance of bioreactions involving large molecules by gel-entrapped microorganisms. This paper describes a quantitative investigation of such diffusional limitations in agar gel membranes. Sucrose and commercial dextran fractions with (weight-average) molecular weights ranging from 10,000 to 2,000,000 Da were used as standard diffusants. For all tested solutes but sucrose, the values of the agar/water partition coefficients highlighted steric hindrance at the entrance of the membrane pores. The effective diffusivity of sucrose in agar was similar to that in water. All dextran fractions, however, displayed restricted diffusion in the agar membranes. Their effective diffusivities were a decreasing function of the agar content of the gel membrane (0.5, 1.0, or 1.5% w/v). The effective diffusivity in a given membrane decreased as the molecular weight of the diffusing molecule increased. T500 (ucbar|M_w = 470,000 Da) and ucbar|M_w = 1,950,000 Da) fractions were unable to diffuse through 1.0 or 1.5% agar membranes. The diffusion data did not agree with the classical (Renkin) model for a hard sphere diffusing through a cylindrical pore. These results are discussed in terms of gel and diffusant characteristics.     

On the weak dependence of water diffusivity on the degree of hydrophobicity of acetylated hydroxypropyl xylan

We used molecular dynamics (MD) simulation to study the diffusion of water at low concentrations in a series of chemically modified xylans, a major hemicellulose, including hydroxypropyl xylan (HPX) and acetoxypropyl xylan (APX) which is essentially acetylated HPX, with different degrees of acetylation (i.e., different degrees of hydrophobicity) at 400 K, a temperature well above the glass transition temperatures of the materials. We used one HPX and three APX models. The three APX models were constructed using the HPX model by substituting 1, 2 or 3 hydroxyl moieties on its repeating units, respectively.     

Multiscale simulation studies of geometrical effects on solution transport through nanopores

Due to intrinsic properties, solid-state nanopores are widely used in nanopore technology. Different geometries (cylindrical (CY), hourglass (HG) and conical (CO)) of artificial nanopores have been fabricated and studied. Each was found to promote different transport abilities experimentally. To explore such pore effects, the combination of finite element (FE) and molecular dynamics (MD) simulations with applied electric filed (150 mV) were performed. The dimension of anion-selective protein pore was used as a nanopore template. Different pore geometries with a narrowest diameter ranging from 1.8 to 1.8 μm were studied here. Firstly, we found that the narrowest regions at a pore orifice in CO and constriction site in HG maximise water velocity and consequently control a water flow rate. Secondly, CY triggers the highest water flux, but low ion selectivity, whilst the funnel-like geometries (HG and CO) enhance the ion selectivity significantly. Both HG and CO show similar degrees of permeant flux and selectivity. The orifice and constriction site in CO and HG are the main player for selectivity and permeation control. Thirdly, the transport properties are tuneable by changing the flow direction in asymmetric CO pore. The tip-to-base flow in CO obviously promotes stronger anion selectivity than the base-to-tip one.