RESUMEN
The molecular structure of lubricant additives controls not only their adsorption and dissociation behaviour at the nanoscale, but also their ability to reduce friction and wear at the macroscale. Here, we show using nonequilibrium molecular dynamics simulations with a reactive force field that tri(s-butyl)phosphate dissociates much faster than tri(n-butyl)phosphate when heated and compressed between sliding iron surfaces. For both molecules, dissociative chemisorption proceeds through cleavage of carbon-oxygen bonds. The dissociation rate increases exponentially with temperature and stress. When the rate-temperature-stress data are fitted with the Bell model, both molecules have similar activation energies and activation volumes and the higher reactivity of tri(s-butyl)phosphate is due to a larger pre-exponential factor. These observations are consistent with experiments using the antiwear additive zinc dialkyldithiophosphate. This study represents a crucial step towards the virtual screening of lubricant additives with different substituents to optimise tribological performance.
RESUMEN
The goal of the eighth industrial fluid properties simulation challenge was to test the ability of molecular simulation methods to predict the adsorption of organic adsorbates in activated carbon materials. In particular, the eighth challenge focused on the adsorption of perfluorohexane in the activated carbon BAM-109. Entrants were challenged to predict the adsorption in the carbon at 273 K and relative pressures of 0.1, 0.3, and 0.6. The predictions were judged by comparison to a benchmark set of experimentally determined values. Overall good agreement and consistency were found between the predictions of most entrants.
RESUMEN
The primary goal of the eighth industrial fluid properties simulation challenge was to test the ability of molecular simulation methods to predict the adsorption of organic adsorbates in activated carbon materials. The challenge focused on the adsorption of perfluorohexane in the activated carbon standard BAM-P109 (Panne and Thünemann 2010). Entrants were challenged to predict the adsorption of perfluorohexane in the activated carbon at a temperature of 273 K and at relative pressures of 0.1, 0.3, and 0.6. The relative pressure (P/Po) is defined as that relative to the bulk saturation pressure predicted by the fluid model at a given temperature (273 K in this case). The predictions were judged by comparison to a set of experimentally determined values, which are published here for the first time and were not disclosed to the entrants prior to the challenge. Benchmark experimental studies, described herein, were also carried out and provided to entrants in order to aid in the development of new force fields and simulation methods to be employed in the challenge. These studies included argon, carbon dioxide, and water adsorption in the BAM-P109 activated carbon as well as X-ray diffraction, X-ray microtomography, photoelectron spectroscopy, and atomic emission spectroscopy studies of BAM-P109. Several concurrent studies were carried out for the BAM-P108 activated carbon (Panne and Thünemann 2010). These are included in the current manuscript for comparison.
RESUMEN
We describe the development of a density-dependent transferable coarse-grain model of crystalline hexahydro-1,3,5-trinitro-s-triazine (RDX) that can be used with the energy conserving dissipative particle dynamics method. The model is an extension of a recently reported one-site model of RDX that was developed by using a force-matching method. The density-dependent forces in that original model are provided through an interpolation scheme that poorly conserves energy. The development of the new model presented in this work first involved a multi-objective procedure to improve the structural and thermodynamic properties of the previous model, followed by the inclusion of the density dependency via a conservative form of the force field that conserves energy. The new model accurately predicts the density, structure, pressure-volume isotherm, bulk modulus, and elastic constants of the RDX crystal at ambient pressure and exhibits transferability to a liquid phase at melt conditions.
RESUMEN
A suite of computational tools is described for particle-based mesoscale simulations of the nonequilibrium dynamics of energetic solids, including mechanical deformation, phase transitions, and chemical reactivity triggered by shock or thermal loading. The method builds upon our recent advances both in generating coarse-grain models under high strains and in developing a variant of dissipative particle dynamics (DPD) that includes chemical reactions. To describe chemical reactivity, a coarse-grain particle equation-of-state was introduced into the constant-energy DPD variant that rigorously treats complex chemical reactions and the associated chemical energy release. As illustration of these developments, we present simulations of shock compression of an RDX crystal and its thermal decomposition under high temperatures. We also discuss our current efforts toward a highly scalable domain-decomposition implementation that extends applicability to micrometer-size simulations. With appropriate parametrization, the method is applicable to other materials whose dynamic response is driven by microstructural heterogeneities.
RESUMEN
Using molecular simulation, four types of nanoporous carbons are examined as adsorbents for the separation of CO(2)/CH(4) mixtures at ambient temperature and pressures up to 10 MPa. First, the adsorption selectivity of CO(2) is investigated in carbon slit pores and single-walled carbon nanotube bundles in order to find the optimal pore dimensions for CO(2) separation. Then, the adsorptive properties of the optimized slit pore and nanotube bundle are compared with two realistic nanoporous carbon models: a carbon replica of zeolite Y and an amorphous carbon. For the four carbon models, adsorption isotherms and isosteric heats of adsorption are presented for both pure components and mixtures. Special attention is given to the calculation of excess isotherms and isosteric heats, which are necessary to assess the performance of model nanoporous materials in the context of experimental measurements. From these results, we discuss the impact that variables such as pore size, pore morphology, pressure and mixture composition have on the performance of nanoporous carbons for CO(2) separation.
RESUMEN
Molecular modeling at the electronic and atomistic levels plays an important and complementary role to experimental studies of confinement effects. Theory and atomistic simulation can provide fundamental understanding, determine the limits of well known macroscopic laws such as Kelvin's equation, provide predictions for systems that are difficult to study via experiment (e.g. adsorption of highly toxic gases), and can be used to gain detailed molecular level information that may not be accessible in the laboratory (e.g. the local structure and composition of confined phases). We describe the most important and useful methods that are based firmly on quantum mechanics and statistical mechanics, including ab intio and classical density functional theories, and Monte Carlo and molecular dynamics simulation. We discuss their strengths and limitations. We then describe examples of applications of these methods to adsorption and equilibrium properties, including testing the Kelvin equation, determination of pore size distributions and capillary phenomena. Applications to self and transport diffusion, including single-file and anomalous diffusion, and viscous flow in nanoporous materials are described. The use of these methods to understand confinement effects on chemical reactions in heterogeneous media is treated, including effects on reaction equilibria, rates and mechanism. Finally we discuss the current status of molecular modeling in this area, and the outlook and future research needs for the next few years. The treatment is suitable for the general technical reader.
RESUMEN
The transition from single-file diffusion to Fickian diffusion in narrow cylindrical pores is investigated for systems of rigid single-walled armchair carbon nanotubes, solvated with binary mixtures of Lennard-Jones fluids (Ar/Ne, Ar/Kr, and Ar/Xe). A range of effects is examined including the mixture concentration, the size ratio of the two components, and the nanotube diameter. The transition from single-file to Fickian diffusion in varying carbon nanotube diameters is analyzed in terms of the Fickian self-diffusivity and the single-file mobility of the mixture components. It is found that the single-file to Fickian carbon nanotube transition diameter is a unique property of the individual molecule's diameter and remains unchanged regardless of the mixture composition. In applications of binary mixtures, each component may crossover from single-file to Fickian diffusion in a different carbon nanotube diameter, giving rise to bimodal diffusion in some nanotubes. This transition allows for one species to diffuse in single-file while the other diffuses by a Fickian mechanism, yielding orders of magnitude difference between the self-diffusional rates of the two molecules. This phenomenon might be further extended to alter the diffusional motion of molecules in nanoporous materials.
RESUMEN
The adsorption and diffusion mechanisms of argon at 120 K were examined in a (25,0) single-walled carbon nanotube (SWCNT) bundle using a combination of Grand Canonical Monte Carlo and microcanonical molecular dynamics simulations. Interstices between the SWCNTs provided the most energetically favorable adsorption sites and filled completely at low relative pressure, followed by adsorption in the SWCNTs. We calculated the self-diffusivities from the average mean squared displacements of argon molecules. In both flexible and rigid bundles, we observed a bimodal diffusion mechanism, with single-file diffusion occurring in the interstitial sites and Fickian diffusion in the SWCNTs. Strong system size effects were observed in our simulations. The largest system sizes showed very little influence of the nanotube flexibility on the diffusion of argon even at the lowest pressures studied.
RESUMEN
This work focuses on the dynamic properties and energetics of the protein/peptide drug during its transport through carbon nanotubes (CNTs). A systematic study was performed on the interaction between the peptide and the CNTs. In the molecular dynamics (MD) simulations, the protein/peptide molecule Zadaxin is observed to be encapsulated inside the nanotube after its spontaneous insertion and oscillates around the center of the tube, where the van der Waals interaction energy is observed to be a minimum. Furthermore, it is found by performing steered MD simulations that the pulling force applied to the peptide reaches a maximum value, which demonstrates the ability of the CNTs to trap protein/peptide drugs. Such effects, attributed to van der Waals interactions, can be influenced by varying the lengths and diameters of the CNTs. Longer nanotubes provide a broader area to trap the peptide, while smaller nanotubes are able to encapsulate the peptide with a deeper interaction energy well. This investigation provides insights into nanoscale pharmaceutical drug delivery devices.