A comparison of implicit- and explicit-solvent simulations of self-assembly in block copolymer and solute systems

We have developed explicit- and implicit-solvent models for the flash nanoprecipitation process, which involves rapid coprecipitation of block copolymers and solutes by changing solvent quality. The explicit-solvent model uses the dissipative particle dynamics (DPD) method and the implicit-solvent model uses the Brownian dynamics (BD) method. Each of the two models was parameterized to match key properties of the diblock copolymer (specifically, critical micelle concentration, diffusion coefficient, polystyrene melt density, and polyethylene glycol radius of gyration) and the hydrophobic solute (aqueous solubility, diffusion coefficient, and solid density). The models were simulated in the limit of instantaneous mixing of solvent with antisolvent. Despite the significant differences in the potentials employed in the implicit- and explicit-solvent models, the polymer-stabilized nanoparticles formed in both sets of simulations are similar in size and structure; however, the dynamic evolution of the two simulations is quite different. Nanoparticles in the BD simulations have diffusion coefficients that follow Rouse behavior (D ∝ M(-1)), whereas those in the DPD simulations have diffusion coefficients that are close to the values predicted by the Stokes-Einstein relation (D ∝ R(-1)). As the nanoparticles become larger, the discrepancy between diffusion coefficients grows. As a consequence, BD simulations produce increasingly slower aggregation dynamics with respect to real time and result in an unphysical evolution of the nanoparticle size distribution. Surface area per polymer of the stable explicit-solvent nanoparticles agrees well with experimental values, whereas the implicit-solvent nanoparticles are stable when the surface area per particle is roughly two to four times larger. We conclude that implicit-solvent models may produce questionable results when simulating nonequilibrium processes in which hydrodynamics play a critical role.     

Strong correlations and Fickian water diffusion in narrow carbon nanotubes

The authors have used atomistic molecular dynamics (MD) simulations to study the structure and dynamics of water molecules inside an open ended carbon nanotube placed in a bath of water molecules. The size of the nanotube allows only a single file of water molecules inside the nanotube. The water molecules inside the nanotube show solidlike ordering at room temperature, which they quantify by calculating the pair correlation function. It is shown that even for the longest observation times, the mode of diffusion of the water molecules inside the nanotube is Fickian and not subdiffusive. They also propose a one-dimensional random walk model for the diffusion of the water molecules inside the nanotube. They find good agreement between the mean-square displacements calculated from the random walk model and from MD simulations, thereby confirming that the water molecules undergo normal mode diffusion inside the nanotube. They attribute this behavior to strong positional correlations that cause all the water molecules inside the nanotube to move collectively as a single object. The average residence time of the water molecules inside the nanotube is shown to scale quadratically with the nanotube length.     

Nonideality in diffusion of ionic and hydrophobic solutes and pair dynamics in water-acetone mixtures of varying composition

We have performed a series of molecular dynamics simulations of water-acetone mixtures containing either an ionic solute or a neutral hydrophobic solute to study the extent of nonideality in the dynamics of these solutes with variation of composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to change nonmonotonically with the composition of the mixtures showing strong nonideality of their dynamics. Also, the extent of nonideality in the diffusion of these charged solutes is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and acetone molecules in these mixtures which show a somewhat similar changes in the solvation characteristics of charged and dipolar solutes with changes of composition of water-acetone mixtures. The diffusion of the hydrophobic solute, however, shows a monotonic increase with increase of acetone concentration showing its different solvation characteristics as compared to the charged and dipolar solutes. The links between the nonideality in diffusion and solvation structures are further confirmed through calculations of the relevant solute-solvent and solvent-solvent radial distribution functions for both ionic and hydrophobic solutes. We have also calculated various pair dynamical properties such as the relaxation of water-water and acetone-water hydrogen bonds and residence dynamics of water molecules in water and acetone hydration shells. The lifetimes of both water-water and acetone-water hydrogen bonds and also the residence times of water molecules are found to increase steadily with increase in acetone concentration. No maximum or minimum was found in the composition dependence of these pair dynamical quantities. The lifetimes of water-water hydrogen bonds are always found to be longer than that of acetone-water hydrogen bonds in these mixtures. The residence times of water molecules are also found to follow a similar trend.     

Nonideality in diffusion of ionic and neutral solutes and hydrogen bond dynamics in dimethyl sulfoxide-chloroform mixtures of varying composition

Molecular dynamics simulations of charged and neutral solutes in dimethyl sulfoxide (DMSO)-chloroform mixtures reveal pronounced nonideality in the solute diffusion with changes of composition of the mixtures. The diffusion coefficient of the anionic solute first decreases, passes through a minimum at DMSO mole fraction of about 0.50, and then increases to reach its value for pure DMSO. The diffusion coefficients of the cationic and neutral solutes are found to decrease with increase in DMSO content of the solvent mixture. The extent of nonideality in the diffusion and orientational relaxation of solvent molecules is found to be somewhat stronger than that in diffusion of the anionic solute in these mixtures. We have also calculated the relaxation of hydrogen bonds formed between DMSO and chloroform molecules. The lifetimes of DMSO-chloroform hydrogen bonds are found to increase monotonically with increase in DMSO concentration. The average number of hydrogen bonds and their average energies are also computed. It is found that an increase in DMSO concentration causes a decrease in the number of DMSO-chloroform hydrogen bonds per DMSO or chloroform molecules but increases the strength of these hydrogen bonds.     

Impact of the solvent on the conformational isomerism of calix[4]arenes: a study based on continuum solvation models

The influence of solvation on the conformational isomerism of calix[4]arene and p-tert-butylcalix[4]arene has been investigated by using the continuum model reported by Miertus, Scrocco, and Tomasi (MST). The quantum mechanical (QM) and semiclassical (SC) formalisms of the MST model have been considered for two different solvents (chloroform and water). The suitability of the QM-MST and SC-MST methods has been examined by comparison with previous results derived from classical molecular dynamics (MD) simulations with explicit solvent molecules. The application of the continuum model to the solute configurations generated by using in vacuo classical MD simulations provides a fast strategy to evaluate the effects of the solvent on the conformational preferences of calixarenes. These encouraging results allow us to propose the use of continuum models to solutes with complex molecular structures, which are traditionally studied by MD simulations.     

N2,O2水溶液光谱的分子动力学模拟

用分子动力学模拟方法研究了N2和O2水溶液的光谱性质．给出了能描述分子内部运动的溶质－溶剂相互作用势．对溶质和溶剂原子的速度自相关函数（VACF）作了计算．讨论了所得VACF的性质并计算了其谱密度．溶质分子振动谱出现的红移，与液态N2，O2的Raman实验结果相吻合．模拟得出的转动谱表明了溶剂分子对溶质转动运动的阻滞，模拟结果也表明VACF计算对溶液和液体光谱的研究十分有效．     

BROMOC suite: Monte Carlo/Brownian dynamics suite for studies of ion permeation and DNA transport in biological and artificial pores with effective potentials

The transport of ions and solutes by biological pores is central for cellular processes and has a variety of applications in modern biotechnology. The time scale involved in the polymer transport across a nanopore is beyond the accessibility of conventional MD simulations. Moreover, experimental studies lack sufficient resolution to provide details on the molecular underpinning of the transport mechanisms. BROMOC, the code presented herein, performs Brownian dynamics simulations, both serial and parallel, up to several milliseconds long. BROMOC can be used to model large biological systems. IMC‐MACRO software allows for the development of effective potentials for solute–ion interactions based on radial distribution function from all‐atom MD. BROMOC Suite also provides a versatile set of tools to do a wide variety of preprocessing and postsimulation analysis. We illustrate a potential application with ion and ssDNA transport in MspA nanopore. © 2014 Wiley Periodicals, Inc.     

Breakdown of the Stokes-Einstein relationship: role of interactions in the size dependence of self-diffusivity

Einstein and others derived the reciprocal dependence of the self-diffusivity D on the solute radius r(u) for large solutes based on kinetic theory. We examine here (a) the range of r(u) over which Stokes-Einstein (SE) dependence is valid and (b) the precise dependence for small solutes outside of the SE regime. We show through molecular dynamics simulations that there are two distinct regimes for smaller solutes: (i) the interaction or Levitation effect (LE) regime for solutes of intermediate size and (ii) the D proportional, variant 1/r(u)(2) for still smaller solutes. We show that as the solute-solvent size ratio decreases, the breakdown in the Stokes-Einstein relationship leading to the LE regime has its origin in dispersion interaction between the solute and the solvent. These results explain reports of enhanced solute diffusion in solvents existing in the literature seen for small solutes for which no explanation exists.     

Single particle and pair dynamics in water-formic acid mixtures containing ionic and neutral solutes: nonideality in dynamical properties

A series of molecular dynamics simulations of water-formic acid mixtures containing either an ionic solute or a neutral hydrophobic solute has been performed to study the extent of nonideality in the dynamics of these solutes for varying composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to show nonideal behavior with variation of composition, and similar nonideality is also observed for the diffusion and orientational relaxation of solvent molecules in these mixtures. The diffusion coefficient of a neutral hydrophobic solute, however, decreases monotonically with increase in water concentration. We have also investigated some of the pair dynamical properties such as water-water and water-formic acid hydrogen bond relaxation and residence dynamics of water molecules in water and formic acid hydration shells. The lifetimes of water-water hydrogen bonds are found to be longer than those between formic acid carbonyl oxygen-water hydrogen bonds, whereas the lifetimes of formic acid hydroxyl hydrogen-water hydrogen bonds are longer than those of water-water hydrogen bonds. In general, the hydrogen bond lifetimes for both water-water and water-formic acid hydrogen bonds are found to decrease with increase in water concentration. Residence times of water molecules also show the same trend with increase in formic acid concentration. Interestingly, these pair dynamical properties show a monotonic dependence on composition without any maximum or minimum and behave almost ideally with respect to changes in the composition of the mixtures. The present calculations are performed with fixed-charge nonpolarizable models of the solvent and solute molecules without taking into account many-body polarization effects in an explicit manner.     

Molecular diffusion in liquid crystals and chiral discrimination. II. Model calculations

We perform model calculations of the macroscopic diffusion coefficient for a solute moving in a chiral nematic (cholesteric) liquid crystal (LC) phase applying the methodology developed by Frezzato et al. [J. Chem. Phys. 122, 164904 (2005)]. Three types of solutes with different features are studied: ellipsoid [roto-translational coupling (RTC) absent], bent rod (RTC present), and two-blade propeller (with RTC and chiral shape). For each prototype molecule we estimate the effect of cholesteric helix pitch and local order on the diffusion along the helix axis. For the ellipsoidal particle we find that translational diffusion is slowed down by rotation around the short axis. For the chiral solute we show that the enantiomer with shape chirality opposite to that of the LC phase is slowed down more than the other. This provides a proof of principle of the possibility of separating the two enantiomers via transport in a suitable chiral medium.     

Atoms-to-microns model for small solute transport through sticky nanochannels

Modeling the transport of solutes through fluidic systems that have adsorbing surfaces is challenging due to the range of length and time scales involved. The components of such systems typically have dimensions from hundreds of nanometres to microns, whereas adsorption of solutes is sensitive to the atomic-scale structure of the solutes and surfaces. Here, we describe an atomic-resolution Brownian dynamics method for modeling the transport of solutes through sticky nanofluidic channels. Our method can fully recreate the results of all-atom molecular dynamics simulations at a fraction of the computational cost of the latter, which makes simulations of micron-size channels at a millisecond time scale possible without losing information about the atomic-scale features of the system. We demonstrate the capability of our method by simulating the rise and fall of solute concentration in sub-micron-long sticky nanochannels, showing that the atomic-scale features of the channels' surfaces have a dramatic effect on the kinetics of solute transport in and out of the channels. We expect our method to find applications in design and optimization of micro and nanofluidic systems for solute-specific transport and to complement existing approaches to modeling lab-on-a-chip devices by providing atomic scale information at a low computational cost.     

Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules

The aqueous solvation free energies of ionized molecules were computed using a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1, MNDO, and PM3 semiempirical molecular orbital methods for the solute molecule and the TIP3P molecular mechanics model for liquid water. The present work is an extension of our model for neutral solutes where we assumed that the total free energy is the sum of components derived from the electrostatic/polarization terms in the Hamiltonian plus an empirical “nonpolar” term. The electrostatic/polarization contributions to the solvation free energies were computed using molecular dynamics (MD) simulation and thermodynamic integration techniques, while the nonpolar contributions were taken from the literature. The contribution to the electrostatic/polarization component of the free energy due to nonbonded interactions outside the cutoff radii used in the MD simulations was approximated by a Born solvation term. The experimental free energies were reproduced satisfactorily using variational parameters from the vdW terms as in the original model, in addition to a parameter from the one-electron integral terms. The new one-electron parameter was required to account for the short-range effects of overlapping atomic charge densities. The radial distribution functions obtained from the MD simulations showed the expected H-bonded structures between the ionized solute molecule and solvent molecules. We also obtained satisfactory results by neglecting both the empirical nonpolar term and the electronic polarization of the solute, i.e., by implementing a nonpolarization model. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1028–1038, 1999     

Effects of solute structure on local solvation and solvent interactions: results from UV/vis spectroscopy and molecular dynamics simulations

Solvation of heterocyclic amines in CO(2)-expanded methanol (MeOH) has been explored with UV/vis spectroscopy and molecular dynamics (MD) simulations. A synergistic study of experiments and simulations allows exploration of solute and solvent effects on solvation and the molecular interactions that affect absorption. MeOH-nitrogen hydrogen bonds hinder the n-pi* transition; however, CO(2) addition causes a blue shift relative to MeOH because of Lewis acid/base interactions with nitrogen. Effects of solute structure are considered, and very different absorption spectra are obtained as nitrogen positions change. MD simulations provide detailed solvent clustering behavior around the solute molecules and show that the local solvent environment and ultimately the spectra are sensitive to the solute structure. This work demonstrates the importance of atomic-level information in determining the structure-property relationships between solute structure, local salvation, and solvatochromism.     

A cellular automata model of water structuring by a chiral solute

The organization of water around a solute molecule with surface features of varying hydropathic states is studied. A stationary solute molecule and mobile water solvent molecules are modeled using cellular automata dynamics. It is shown that varying hydropathic states of solute molecule surface features influence the relative affinities of water for these features. In the case of a simulated chiral solute, a chiral pattern of associated water molecules binding to the surface is produced. This finding is in agreement with published simulations and circular dichroism measurements. A pattern of water molecules at locations beyond the surface of the solute molecules is detected, evidence of an emergent property in this solvent-solute system.     

A molecular dynamics study of chirality transfer: The impact of a chiral solute on an achiral solvent

A solvation shell may adapt to the presence of a chiral solute by becoming chiral. The extent of this chirality transfer and its dependence on the solute and solvent characteristics are explored in this article. Molecular dynamics simulations of solvated chiral analytes form the basis of the analysis. The chirality induced in the solvent is assessed based on a series of related chirality indexes originally proposed by Osipov [M. A. Osipov et al., Mol. Phys. 84, 1193 (1995)]. Two solvents are considered: Ethanol and benzyl alcohol. Ethanol provides insight into chirality transfer when the solvent interacts with the solute primarily by a hydrogen bond. Several ethanol models have been considered starting with a nonpolarizable model, progressing to a fluctuating charge model, and finally, to a fully polarizable model. This progression provides some insights into the importance of solvent polarizability in the transfer of chirality. Benzyl alcohol, by virtue of the aromatic ring, increases the number of potential solvent-solute interactions. Thus, with these two solvents, the issue of compatibility between the solvent and solute is also considered. The solvation of three chiral solutes is examined: Styrene oxide, acenaphthenol, and n-(1-(4-bromophenyl)ethyl)pivalamide (PAMD). All three solutes have the possibility of hydrogen bonding with the solvent, the last two may also form ring-ring interactions, and the last also has multiple hydrogen bonding sites. For PAMD, the impact of conformational averaging is examined by comparing the chirality transfer about rigid and flexible solutes.     

Dynamics of ionic and hydrophobic solutes in water-methanol mixtures of varying composition

We have carried out a series of molecular-dynamics simulations of water-methanol mixtures containing either an ionic or a neutral atomic solute to investigate the effects of composition of the mixture on the diffusion of these solutes. Altogether, we have considered 17 different systems of varying composition ranging from pure water to pure methanol. The diffusion coefficients of ionic solutes are found to show nonideal behavior with variation of composition of the solvent mixture. The extent of nonideality of the solute diffusion is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and methanol molecules in these mixtures and is attributed to the enhanced stability of the hydrogen bonds and formation of interspecies complexes in the mixtures. The neutral solute shows characteristics of hydrophobic solvation and its diffusion decreases monotonically with increase of methanol concentration. The present simulation results are compared with those of experiments wherever available.     

Effect of Size and Temperature on Water Dynamics inside Carbon Nano-Tubes Studied by Molecular Dynamics Simulation

Water transport inside carbon nano-tubes (CNTs) has attracted considerable attention due to its nano-fluidic properties, its importance in nonporous systems, and the wide range of applications in membrane desalination and biological medicine. Recent studies show an enhancement of water diffusion inside nano-channels depending on the size of the nano-confinement. However, the underlying mechanism of this enhancement is not well understood yet. In this study, we performed Molecular Dynamics (MD) simulations to study water flow inside CNT systems. The length of CNTs considered in this study is 20 nm, but their diameters vary from 1 to 10 nm. The simulations are conducted at temperatures ranging from 260 K to 320 K. We observe that water molecules are arranged into coaxial water tubular sheets. The number of these tubular sheets depends on the CNT size. Further analysis reveals that the diffusion of water molecules along the CNT axis deviates from the Arrhenius temperature dependence. The non-Arrhenius relationship results from a fragile liquid-like water component persisting at low temperatures with fragility higher than that of the bulk water.     

Molecular diffusion in liquid crystals and chiral discrimination. I. Theory

The possibility of using cholesteric phases for discriminating enantiomers of a chiral solute on the basis of their different transport properties, motivates the investigation of the translational diffusion by taking fully into account the roto-translational coupling. In this article a detailed theoretical analysis is presented for the transport properties evaluated according to the asymptotic limit of the mean-squared displacement. A general relation is derived for the transport coefficients, having as main ingredients the mean-field potential due to the mesophase, and the diffusion tensor with its purely translational and rotational components, and with the blocks describing the roto-translational coupling. The application of the theory to nematic phases shows that the roto-translational coupling generates a dynamical contribution reducing the transport coefficients evaluated by taking into account only the translational diffusion components in the center of diffusion. The theory is also specialized to a cholesteric phase with a given helical pitch for the director arrangement, in a form which is suitable for calculations of model systems of chiral solutes to be presented in a forthcoming paper.     

Mutual diffusion coefficients of heptane isomers in nitrogen: a molecular dynamics study

The accurate knowledge of transport properties of pure and mixture fluids is essential for the design of various chemical and mechanical systems that include fluxes of mass, momentum, and energy. In this study we determine the mutual diffusion coefficients of mixtures composed of heptane isomers and nitrogen using molecular dynamics (MD) simulations with fully atomistic intermolecular potential parameters, in conjunction with the Green-Kubo formula. The computed results were compared with the values obtained using the Chapman-Enskog (C-E) equation with Lennard-Jones (LJ) potential parameters derived from the correlations of state values: MD simulations predict a maximum difference of 6% among isomers while the C-E equation presents that of 3% in the mutual diffusion coefficients in the temperature range 500-1000 K. The comparison of two approaches implies that the corresponding state principle can be applied to the models, which are only weakly affected by the anisotropy of the interaction potentials and the large uncertainty will be included in its application for complex polyatomic molecules. The MD simulations successfully address the pure effects of molecular structure among isomers on mutual diffusion coefficients by revealing that the differences of the total mutual diffusion coefficients for the six mixtures are caused mainly by heptane isomers. The cross interaction potential parameters, collision diameter σ(12), and potential energy well depth ?(12) of heptane isomers and nitrogen mixtures were also computed from the mutual diffusion coefficients.