scholarly journals Molecular Simulation of the Adsorption and Diffusion in Cylindrical Nanopores: Effect of Shape and Fluid–Solid Interactions

Molecules ◽  
2019 ◽  
Vol 24 (3) ◽  
pp. 608 ◽  
Author(s):  
Harry Cárdenas ◽  
Erich Müller

We report on molecular simulations of model fluids composed of three tangentially bonded Lennard-Jones interaction sites with three distinct morphologies: a flexible “pearl-necklace” chain, a rigid “stiff” linear configuration, and an equilateral rigid triangular ring. The adsorption of these three models in cylindrical pores of diameters 1, 2, and 3 nm and with varying solid–fluid strength was determined by direct molecular dynamics simulations, where a sample pore was placed in contact with a bulk fluid. Adsorption isotherms of Type I, V, and H1 were obtained depending on the choice of pore size and solid–fluid strength. Additionally, the bulk-phase equilibria, the nematic order parameter of the adsorbed phase, and the self-diffusion coefficient in the direction of the pore axis were examined. It was found that both the molecular shape and the surface attractions play a decisive role in the shape of the adsorption isotherm. In general, the ring molecules showed a larger adsorption, while the fully flexible model showed the smallest adsorption. Morphology and surface strength were found to have a lesser effect on the diffusion of the molecules. An exceptional high adsorption and diffusion, suggesting an enhanced permeability, was observed for the linear stiff molecules in ultraconfinement, which was ascribed to a phase transition of the adsorbed fluid into a nematic liquid crystal.

2017 ◽  
Vol 82 (6) ◽  
pp. 681-694 ◽  
Author(s):  
Saeid Yeganegi ◽  
Vahid Moeini ◽  
Zohreh Doroodi

Molecular dynamics simulations were performed to study thermodynamics and structural properties of expanded caesium fluid. Internal pressure, radial distribution functions (RDFs), coordination numbers and diffusion coefficients have been calculated at temperature range 700?1600 K and pressure range 100?800 bar. We used the internal pressure to predict the metal?non-metal transition occurrence region. RDFs were calculated at wide ranges of temperature and pressure. The coordination numbers decrease and positions of the first peak of RDFs slightly increase as the temperature increases and pressure decreases. The calculated self-diffusion coefficients at various temperatures and pressures show no distinct boundary between Cs metallic fluid and its expanded fluid where it continuously increases with temperature.


Author(s):  
Nicola Molinari ◽  
Jonathan P. Mailoa ◽  
Boris Kozinsky

We show that strong cation-anion interactions in a wide range of lithium-salt/ionic liquid mixtures result in a negative lithium transference number, using molecular dynamics simulations and rigorous concentrated solution theory. This behavior fundamentally deviates from the one obtained using self-diffusion coefficient analysis and agrees well with experimental electrophoretic NMR measurements, which accounts for ion correlations. We extend these findings to several ionic liquid compositions. We investigate the degree of spatial ionic coordination employing single-linkage cluster analysis, unveiling asymmetrical anion-cation clusters. Additionally, we formulate a way to compute the effective lithium charge that corresponds to and agrees well with electrophoretic measurements and show that lithium effectively carries a negative charge in a remarkably wide range of chemistries and concentrations. The generality of our observation has significant implications for the energy storage community, emphasizing the need to reconsider the potential of these systems as next generation battery electrolytes.<br>


2007 ◽  
Vol 130 (1) ◽  
pp. 111-116 ◽  
Author(s):  
Artem B. Mamonov ◽  
Rob D. Coalson ◽  
Mark L. Zeidel ◽  
John C. Mathai

Determining the mechanisms of flux through protein channels requires a combination of structural data, permeability measurement, and molecular dynamics (MD) simulations. To further clarify the mechanism of flux through aquaporin 1 (AQP1), osmotic pf (cm3/s/pore) and diffusion pd (cm3/s/pore) permeability coefficients per pore of H2O and D2O in AQP1 were calculated using MD simulations. We then compared the simulation results with experimental measurements of the osmotic AQP1 permeabilities of H2O and D2O. In this manner we evaluated the ability of MD simulations to predict actual flux results. For the MD simulations, the force field parameters of the D2O model were reparameterized from the TIP3P water model to reproduce the experimentally observed difference in the bulk self diffusion constants of H2O vs. D2O. Two MD systems (one for each solvent) were constructed, each containing explicit palmitoyl-oleoyl-phosphatidyl-ethanolamine (POPE) phospholipid molecules, solvent, and AQP1. It was found that the calculated value of pf for D2O is ∼15% smaller than for H2O. Bovine AQP1 was reconstituted into palmitoyl-oleoyl-phosphatidylcholine (POPC) liposomes, and it was found that the measured macroscopic osmotic permeability coefficient Pf (cm/s) of D2O is ∼21% lower than for H2O. The combined computational and experimental results suggest that deuterium oxide permeability through AQP1 is similar to that of water. The slightly lower observed osmotic permeability of D2O compared to H2O in AQP1 is most likely due to the lower self diffusion constant of D2O.


2021 ◽  
Vol 29 ◽  
pp. 95-115
Author(s):  
Rafal Kozubski ◽  
Graeme E. Murch ◽  
Irina V. Belova

We review the results of our Monte Carlo simulation studies carried out within the past two decades in the area of atomic-migration-controlled phenomena in intermetallic compounds. The review aims at showing the high potential of Monte Carlo methods in modelling both the equilibrium states of the systems and the kinetics of the running processes. We focus on three particular problems: (i) the atomistic origin of the complexity of the ‘order-order’ relaxations in γ’-Ni3Al; (ii) surface-induced ordering phenomena in γ-FePt and (iii) ‘order—order’ kinetics and self-diffusion in the ‘triple-defect’ β-NiAl. The latter investigation demonstrated how diverse Monte Carlo techniques may be used to model the phenomena where equilibrium thermodynamics interplays and competes with kinetic effects.


2019 ◽  
Author(s):  
Thiago José Pinheiro dos Santos ◽  
Charlles Abreu ◽  
Bruno Horta ◽  
Frederico W. Tavares

Mass transport coefficients play an important role in process design and in compositional grading of oil reservoirs. As experimental measurements of these properties can be costly and hazardous, Molecular Dynamics simulations emerge as an alternative approach. In this work, we used Molecular Dynamics to calculate the self-diffusion coefficients of methane/n-hexane mixtures at different conditions, in both liquid and supercritical phases. We evaluated how the finite box size and the choice of the force field affect the calculated properties at high pressures. Results show a strong dependency between self-diffusion and the simulation box size. The Yeh-Hummer analytical correction [J. Phys. Chem. B, 108, 15873 (2004)] can attenuate this effect, but sometimes makes the results depart from experimental data due to issues concerning the force fields. We have also found that different all-atom and united-atom models can produce biased results due to caging effects and to different dihedral configurations of the n-alkane.


Sign in / Sign up

Export Citation Format

Share Document