Search

We study the shock-induced hot spot formation mechanism of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine-based explosives by molecular dynamics, compare different kinds of desensitizers and different shock velocities. A set of programs is written to calculate the physical picture of shock loading. Based on the simulations and analyses, the hot spots are found at the interface and are heated by plastic work in three ways: the interface intrinsic dissipation, the pore collapse, and the coating layer deformation. The work/heat transition rate is proved to be increasing with a loading speed. [ABSTRACT FROM AUTHOR]

This paper demonstrates a method for calculating the tension of a system with a curved interface from a molecular dynamics simulation. To do so, the pressure of a subset of the system is determined by applying a local (virtual) mechanical deformation, fitting the response to that of a bulk fluid, and then using the Young-Laplace equation to infer the tension of the interface. The accuracy of the method is tested by calculating the local pressure of a series of water simulations at various external pressures. The tension of a simulated curved octane-water interface is computed with the method and compares well with the planar tension (≈ 46.7 dyn/cm). Finally, an ambiguity is resolved between the Harasima and Irving-Kirkwood methods of calculating the local pressure as a means for computing the tension. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations have been carried out to obtain the interfacial and coexistence properties of soft-sphere attractive Yukawa (SAY) fluids with short attraction range, κ = 10, 9, 8, 7, 6, and 5. All our simulation results are new. These data are also compared with the recently reported results in the literature of hard-core attractive Yukawa (HAY) fluids. We show that the interfacial and coexistence properties of both potentials are different. For the surveyed systems, here we show that all coexistence curves collapse into a master curve when we rescale with their respective critical points and the surface tension curves form a single master curve when we plot γ* vs. T/Tc. [ABSTRACT FROM AUTHOR]

Bulk and surface properties of the ionic liquids 1-alkyl-3-methyl-imidazolium iodides ([Cnmim]I) were simulated by classical molecular dynamics using all atom non-polarizable force field (n = 4, butyl; 6, hexyl; 8, octyl). The structure of ionic liquids were initially optimized by density functional theory and atomic charges obtained by CHELPG method. Reduction of partial atomic charges (by 20% for simulation of density and surface tension, and by 10% for viscosity) found to improve the accuracy, while a non-polarizable force field was applied. Additionally, the simulation ensembles approach the equilibrium faster when the charge reduction is applied. By these refined force field parameters, simulated surface tensions in the range of 323-393 k are quite in agreement with the experiments. Simulation of temperature dependent surface tension of [C4mim]I well beyond room temperature (up to 700 K) permits prediction of the critical temperature in agreement with that predicted from experimental surface tension data. Simulated densities in the range of 298-450 K for the three ionic liquids are within 0.8% of the experimental data. Structural properties for [C4mim]I were found to be in agreement with the results of Car-Parrinello molecular dynamics simulation we performed, which indicates a rather well-structured cation-anion interaction and occurs essentially through the imidazolium ring cation. Diffusion coefficient changes with alkyl chain length in the order of [C8mim]I > [C6mim]I > [C4mim]I for the cation and the anion. Formation of a dense domain in subsurface region is quite evident, and progressively becomes denser as the alkyl chain length increases. Bivariate orientational analysis was used to determine the average orientation of molecule in ionic liquids surface, subsurface, and bulk regions. Dynamic bisector-wise and side-wise movement of the imodazolium ring cation in the surface region can be deduced from the bivariate maps. Atom-atom density profile and bivariate analysis indicate that the imidazolium cation takes a spoon like configuration in the surface region and the tilt of alkyl group is a function length of alkyl chain exposing as linear as possible to the vapor phase. [ABSTRACT FROM AUTHOR]

Data for the flow rate of water in carbon nanopores is widely scattered, both in experiments and simulations. In this work, we aim at precisely quantifying the characteristic large slip length and flow rate of water flowing in a planar graphene nanochannel. First, we quantify the slip length using the intrinsic interfacial friction coefficient between water and graphene, which is found from equilibrium molecular dynamics (EMD) simulations. We then calculate the flow rate and the slip length from the streaming velocity profiles obtained using non-equilibrium molecular dynamics (NEMD) simulations and compare with the predictions from the EMD simulations. The slip length calculated from NEMD simulations is found to be extremely sensitive to the curvature of the velocity profile and it possesses large statistical errors. We therefore pose the question: Can a micrometer range slip length be reliably determined using velocity profiles obtained from NEMD simulations? Our answer is 'not practical, if not impossible' based on the analysis given as the results. In the case of high slip systems such as water in carbon nanochannels, the EMD method results are more reliable, accurate, and computationally more efficient compared to the direct NEMD method for predicting the nanofluidic flow rate and hydrodynamic boundary condition. [ABSTRACT FROM AUTHOR]

First-principles molecular dynamics (FPMD) simulations and static quantum chemical (QC) calculations are used to evaluate the tensile strengths, σc, of interfaces consisting of (0001) surfaces of α-Al2O3 separated by small organic species. The evaluation of σc with FPMD was achieved by performing simulations in which the simulation cell was extending in a direction normal to the fracture plane until rupture of the interface occurred. The static QC calculations employed an approach which treated fracture of the interface as a competition between uniform extension of the simulation cell and crack formation at the rupture site, which is analogous to that used in the construction of universal binding energy relationships. The results showed that the static QC calculations accurately reproduced the FPMD simulations with respect to tensile strength and the cell extension at which rupture occurred, provided that the rupture site employed in the static calculations matched the site at which rupture occurred during the FPMD simulations. A simple strategy for identifying the rupture site, even in complex systems containing many potential rupture sites, is proposed. Overall, the work extends the calculation of tensile strengths with static QC methods to highly heterogeneous interfaces, thus providing a computationally efficient alternative to demanding FPMD simulations for this purpose. [ABSTRACT FROM AUTHOR]

A new molecular dynamics simulation method, with coupling to external baths, is used to perform equilibrium simulations on polyamide-6,6 trimers nanoconfined between graphene surfaces, in equilibrium with the bulk polymer. The method is coupled with the reverse nonequilibrium molecular dynamics simulation technique to exchange heat in the direction normal to the surfaces. To be able to study the effect of confinement on the heat conductance in nanoconfined pores, in this work a number of simulations on systems with different pore sizes are done. It is concluded that the coefficient of heat conductivity depends on the degree of polymer layering between the surfaces and on the pore width. Our results further indicate a considerable temperature drop at the interface between the surfaces and polymer. The calculated Kapitza lengths depend on the intersurface distance and on the layering of the polymer nanoconfined between the surfaces. [ABSTRACT FROM AUTHOR]

While the interaction of colloidal particles (sizes in excess of 100 nm) with liquid interfaces may be understood in terms of continuum models, which are grounded in macroscopic properties such as surface and line tensions, the behaviour of nanoparticles at liquid interfaces may be more complex. Recent simulations [D. L. Cheung and S. A. F. Bon, Phys. Rev. Lett. 102, 066103 (2009)] of nanoparticles at an idealised liquid-liquid interface showed that the nanoparticle-interface interaction range was larger than expected due, in part, to the action of thermal capillary waves. In this paper, molecular dynamics simulations of a Lennard-Jones nanoparticle in a binary Lennard-Jones mixture are used to confirm that these previous results hold for more realistic models. Furthermore by including attractive interactions between the nanoparticle and the solvent, it is found that the detachment energy decreases as the nanoparticle-solvent attraction increases. Comparison between the simulation results and recent theoretical predictions [H. Lehle and M. Oettel, J. Phys. Condens. Matter 20, 404224 (2008)] shows that for small particles the incorporation of capillary waves into the predicted effective nanoparticle-interface interaction improves agreement between simulation and theory. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations have been used to investigate the kinetics of spontaneous cavitation and crystallization in a Lennard-Jones liquid at negative pressures in the temperature range where these processes compete with each other. The nucleation rate has been calculated in NVE and NpT ensembles by the method of mean lifetime and the transition interface sampling method with parallel path swapping. The data obtained have been used to determine in the framework of classical nucleation theory the value of the ratio of the solid-liquid and the liquid-void interfacial free energy for critical crystals and cavities and the values of their volumes at points where the cavitation rate of the liquid is equal to the rate of its crystallization. [ABSTRACT FROM AUTHOR]

Though the local dielectric constant at interfaces is an important phenomenological parameter in the analysis of surface spectroscopy, its microscopic definition has been uncertain. Here, we present a full molecular theory on the local field at interfaces with the help of molecular dynamics simulation, and thereby provide microscopic basis for the local dielectric constant so as to be consistent to the phenomenological three-layer model of interface systems. To demonstrate its performance, we applied the theory to the water/vapor interface, and obtained the local field properties near the interface where the simple dielectric model breaks down. Some computational issues pertinent to Ewald calculations of the dielectric properties are also discussed. [ABSTRACT FROM AUTHOR]

Commonly, the confinement effects are studied from the grand canonical Monte Carlo (GCMC) simulations from the computation of the density of liquid in the confined phase. The GCMC modeling and chemical potential (μ) calculations are based on the insertion/deletion of the real and ghost particle, respectively. At high density, i.e., at high pressure or low temperature, the insertions fail from the Widom insertions while the performing methods as expanded method or perturbation approach are not efficient to treat the large and complex molecules. To overcome this problem we use a simple and efficient method to compute the liquid's density in the confined medium. This method does not require the precalculation of μ and is an alternative to the GCMC simulations. From the isothermal-isosurface-isobaric statistical ensemble we consider the explicit framework/liquid external interface to model an explicit liquid's reservoir. In this procedure only the liquid molecules undergo the volume changes while the volume of the framework is kept constant. Therefore, this method is described in the NpnAVfT statistical ensemble, where N is the number of particles, pn is the normal pressure, Vf is the volume of framework, A is the surface of the solid/fluid interface, and T is the temperature. This approach is applied and validated from the computation of the density of the methanol and water confined in the mesoporous cylindrical silica nanopores and the MIL-53(Cr) metal organic framework type, respectively. [ABSTRACT FROM AUTHOR]

Density functional molecular dynamics simulations have been performed in the NVT ensemble (moles (N), volume (V) and temperature (T)) on a system formed by ten acetone molecules at a temperature of 2000 K and density ρ = 1.322 g cm-3. These conditions resemble closely those realized at the interface of an acetone vapor bubble in the early stages of supercompression experiments and result in an average pressure of 5 GPa. Two relevant reactive events occur during the simulation: the condensation of two acetone molecules to give hexane-2,5-dione and dihydrogen and the isomerization to the enolic propen-2-ol form. The mechanisms of these events are discussed in detail. [ABSTRACT FROM AUTHOR]

In this paper we calculate the interfacial resistances to heat and mass transfer through a liquid-vapor interface in a binary mixture. We use two methods, the direct calculation from the actual nonequilibrium solution and integral relations, derived earlier. We verify, that integral relations, being a relatively faster and cheaper method, indeed gives the same results as the direct processing of a nonequilibrium solution. Furthermore we compare the absolute values of the interfacial resistances with the ones obtained from kinetic theory. Matching the diagonal resistances for the binary mixture we find that kinetic theory underestimates the cross coefficients. The heat of transfer is, as a consequence, correspondingly larger. [ABSTRACT FROM AUTHOR]

We consider the computation of the interfacial properties of molecular chains from direct simulation of the vapor-liquid interface. The molecules are modeled as fully flexible chains formed from tangentially bonded monomers with truncated Lennard-Jones interactions. Four different model systems comprising of 4, 8, 12, and 16 monomers per molecule are considered. The simulations are performed in the canonical ensemble, and the vapor-liquid interfacial tension is evaluated using the test area and the wandering interface methods. In addition to the surface tension, we also obtain density profiles, coexistence densities, critical temperature and density, and interfacial thickness as functions of temperature, paying particular attention to the effect of the chain length on these properties. According to our results, the main effect of increasing the chain length (at fixed temperature) is to sharpen the vapor-liquid interface and to increase the width of the biphasic coexistence region. As a result, the interfacial thickness decreases and the surface tension increases as the molecular chains get longer. The interfacial thickness and surface tension appear to exhibit an asymptotic limiting behavior for long chains. A similar behavior is also observed for the coexistence densities and critical properties. Our simulation results indicate that the asymptotic regime is reached for Lennard-Jones chains formed from eight monomer segments. We also include a preliminary study on the effect of the cutoff distance on the interfacial properties. Our results indicate that all of the properties exhibit a dependence with the distance at which the interactions are truncated, though the relative effect varies from one property to the other. The interfacial thickness and, more particularly, the interfacial tension are found to be strongly dependent on the particular choice of cutoff, whereas the density profiles and coexistence densities are, in general, less sensitive to the truncation. [ABSTRACT FROM AUTHOR]

The Wolf method for the calculation of electrostatic interactions is applied in a liquid phase and at the liquid-vapor interface of water and its results are compared with those from the Ewald sums method. Molecular dynamics simulations are performed to calculate the radial distribution functions at room temperature. The interface simulations are used to obtain the coexisting densities and surface tension along the coexistence curve. The water model is a flexible version of the extended simple point charge model. The Wolf method gives good structural results, fair coexistence densities, and poor surface tensions as compared with those obtained using the Ewald sums method. [ABSTRACT FROM AUTHOR]

We report molecular dynamics simulations of three globular proteins: ubiquitin, apo-calbindin D9K, and the C-terminal SH2 domain of phospholipase C-γ1 in explicit water. The proteins differ in their overall charge and fold type and were chosen to represent to some degree the structural variability found in medium-sized proteins. The length of each simulation was at least 15 ns, and larger than usual solvent boxes were used. We computed radial distribution functions, as well as orientational correlation functions about the surface residues. Two solvent shells could be clearly discerned about charged and polar amino acids. Near apolar amino acids the water density near such residues was almost devoid of structure. The mean residence time of water molecules was determined for water shells about the full protein, as well as for water layers about individual amino acids. In the dynamic properties, two solvent shells could be characterized as well. However, by comparison to simulations of pure water it could be shown that the influence of the protein reaches beyond 6 Å, i.e., beyond the first two shells. In the first shell (r≤=3.5 Å), the structural and dynamical properties of solvent waters varied considerably and depended primarily on the physicochemical properties of the closest amino acid side chain, with which the waters interact. By contrast, the solvent properties seem not to depend on the specifics of the protein studied (such as the net charge) or on the secondary structure element in which an amino acid is located. While differing considerably from the neat liquid, the properties of waters in the second solvation shell (3.5≤r≤=6 Å) are rather uniform; a direct influence from surface amino acids are already mostly shielded. [ABSTRACT FROM AUTHOR]

We report molecular dynamics (MD) simulations of three protein-water systems (ubiquitin, apo-calbindin D9K, and the C-terminal SH2 domain of phospholipase C-γ1), from which we compute the dielectric properties of the solutions. Since two of the proteins studied have a net charge, we develop the necessary theory to account for the presence of charged species in a form suitable for computer simulations. In order to ensure convergence of the time correlation functions needed for the analysis, the minimum length of the MD simulations was 20 ns. The system sizes (box length, number of waters) were chosen so that the resulting protein concentrations are comparable to experimental conditions. A dielectric component analysis was carried out to analyze the contributions from protein and water to the frequency-dependent dielectric susceptibility χ(ω) of the solutions. Additionally, an even finer decomposition into protein, two solvation shells, and the remaining water (bulk water) was carried out. The results of these dielectric decompositions were used to study protein solvation at mesoscopic resolution, i.e., in terms of protein, first and second solvation layers, and bulk water. This study, therefore, complements the structural and dynamical analyses at molecular resolution that are presented in the companion paper. The dielectric component contributions from the second shell and bulk water are very similar in all three systems. We find that the proteins influence the dielectric properties of water even beyond the second solvation shell, in agreement with what was observed for the mean residence times of water molecules in protein solutions. By contrast, the protein contributions, as well as the contributions of the first solvation shell, are system specific. Most importantly, the protein and the first water shell around ubiquitin and apo-calbindin are anticorrelated, whereas the first water shell around the SH2 domain is positively correlated. [ABSTRACT FROM AUTHOR]

We describe a novel algorithm that includes the effect of host lattice flexibility into molecular dynamics simulations that use rigid lattices. It uses a Lowe-Andersen thermostat for interface-fluid collisions to take the most important aspects of flexibility into account. The same diffusivities and other properties of the flexible framework system are reproduced at a small fraction of the computational cost of an explicit simulation. We study the influence of flexibility on the self-diffusion of simple gases inside single walled carbon nanotubes. Results are shown for different guest molecules (methane, helium, and sulfur hexafluoride), temperatures, and types of carbon nanotubes. We show, surprisingly, that at low loadings flexibility is always relevant. Notably, it has a crucial influence on the diffusive dynamics of the guest molecules. [ABSTRACT FROM AUTHOR]

Using the molecular dynamics simulation technique, we study the effect of a highly ordered physisorbed alkyl monolayer on the surface induced liquid structure near a smooth solid–liquid interface. The system considered here is 1,12-dodecanediol adsorbed at the benzene–graphite interface. We compare the results on the graphite–1,12-dodecanediol–benzene interface to previous simulations on the bare graphite–benzene interface. We find that the surface induced liquid structure near the adsorbate covered interface differs from the structure induced by the bare interface only within the first solvation shell. Beyond the first solvation shell, there is no discernible difference, even though the surface induced structure is still pronounced. In addition, we study the structure of the hydrogen bond stabilized molecular adsorbate network. Here we find that the hydrogen bonds link the 1,12-dodecanediol molecules in a chain-like fashion across the boundaries of the herringbone structure formed by the adsorbate. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations are used to study the structure of the vapor–liquid interface of three-dimensional fluids. Particles interact via a truncated Lennard-Jones pair potential in the absence of external fields. The effect of the surface area on the surface tension is investigated. It is found that the surface tension increases with the decrease of the surface area. However, this finite-size effect is pronounced only in small surface areas. In addition, our simulation results show that the finite-size correction of the surface tension is directly proportional to the reciprocal of the surface area, in accord with the prediction of the Gaussian model of capillary waves. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We have used molecular dynamics simulations to calculate the interfacial tension of hydroxylated SiO2 nanoparticles under different temperatures and solutions (helium and brine with monovalent and divalent salts). In order to benchmark the atomistic model, quartz SiO2 interfacial tension was measured based on inverse gas chromatography under He atmosphere. The experimental interfacial tension values for quartz were found between 0.512 and 0.617 N/m. Our calculated results for the interfacial tension of silica nanoparticles within helium atmosphere was 0.676 N/m, which is higher than the value found for the system containing He/α-quartz (0.478 N/m), but it is similar to the one found for amorphous silica surface. We have also studied the interfacial tension of the nanoparticles in electrolyte aqueous solution for different types and salts concentrations (NaCl, CaCl2, and MgCl2). Our calculations indicate that adsorption properties and salt solutions greatly influence the interfacial tension in an order of CaCl2 > MgCl2 > NaCl. This effect is due to the difference in distribution of ions in solution, which modifies the hydration and electrostatic potential of those ions near the nanoparticle. [ABSTRACT FROM AUTHOR]

Monte Carlo and molecular dynamics computer simulations have been used to study the structure and dynamics of the interlayer aqueous solution in a colloidal sodium laponite clay at 277 K. The system studied has a clay-clay spacing of 34.06 Å, and contains 1200 interlayer water molecules and 24 sodium counterions. The density profiles for interlayer species show two distinct layers of surface water as one moves away from the clay particles. The innermost of these layers is strongly oriented to form hydrogen bonds to the surface oxygen atoms. Radially averaged pair distributions have been calculated as a function of distance from the clay surfaces, and show that throughout our system the water structure is significantly perturbed from the bulk. In particular, we observe an increase in the second nearest-neighbor oxygen-oxygen distance, similar to that reported for low-density water at 268 K [A. K. Soper and M. A. Ricci, Phys. Rev. Lett. 84, 2881 (2000)]. The majority of the sodium counterions are fully hydrated by six water molecules. These hydrated ions have a strong tendency to remain close to the solid surfaces, as so-called "outer-sphere" complexes. However, we also observe cations further from the clay sheets, in the diffuse layer. Diffusion of water and cations in the plane of the clay sheets is comparable to that in the bulk, but is significantly reduced normal to the clay sheets. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]