Simulations of composite carbon films with nanotube inclusions

We study the interfacial structure, stability, and elastic properties of composite carbon films containing nanotubes. Our Monte Carlo simulations show that Van der Waals forces play a vital role in shaping up the interfacial geometry, producing a curved graphitic wall surrounding the tubes. The most stable structures are predicted to have intermediate densities, high anisotropies, and increased elastic moduli compared to pure amorphous carbon films.Comment: 3 pages, 3 figures, to appear in Appl. Phys. Let

Ionic force field optimization based on single-ion and ion-pair solvation properties

Molecular dynamics simulations of ionic solutions depend sensitively on the force fields employed for the ions. To resolve the fine differences between ions of the same valence and roughly similar size and in particular to correctly describe ion-specific effects, it is clear that accurate force fields are necessary. In the past, optimization strategies for ionic force fields either considered single-ion properties (such as the solvation free energy at infinite dilution or the ion-water structure) or ion-pair properties (in the form of ion-ion distribution functions). In this paper we investigate strategies to optimize ionic force fields based on single-ion and ion-pair properties simultaneously. To that end, we simulate five different salt solutions, namely CsCl, KCl, NaI, KF, and CsI, at finite ion concentration. The force fields of these ions are systematically varied under the constraint that the single-ion solvation free energy matches the experimental value, which reduces the two-dimensional $\{\sigma,\epsilon\}$ parameter space of the Lennard Jones interaction to a one dimensional line for each ion. From the finite-concentration simulations, the pair-potential is extracted and the osmotic coefficient is calculated, which is compared to experimental data. We find a strong dependence of the osmotic coefficient on the force field, which is remarkable as the single-ion solvation free energy and the ion-water structure remain invariant under the parameter variation. Optimization of the force field is achieved for the cations Cs$^+$ and K$^+$, while for the anions I$^-$ and F$^-$ the experimental osmotic coefficient cannot be reached. This suggests that in the long run, additional parameters might have to be introduced into the modeling, for example by modified mixing rules.Comment: 29 pages, 6 figures, to appear in J. Chem. Phy

Electrokinetic Lattice Boltzmann solver coupled to Molecular Dynamics: application to polymer translocation

We develop a theoretical and computational approach to deal with systems that involve a disparate range of spatio-temporal scales, such as those comprised of colloidal particles or polymers moving in a fluidic molecular environment. Our approach is based on a multiscale modeling that combines the slow dynamics of the large particles with the fast dynamics of the solvent into a unique framework. The former is numerically solved via Molecular Dynamics and the latter via a multi-component Lattice Boltzmann. The two techniques are coupled together to allow for a seamless exchange of information between the descriptions. Being based on a kinetic multi-component description of the fluid species, the scheme is flexible in modeling charge flow within complex geometries and ranging from large to vanishing salt concentration. The details of the scheme are presented and the method is applied to the problem of translocation of a charged polymer through a nanopores. In the end, we discuss the advantages and complexities of the approach

Going beyond standard mixing rules

Using molecular dynamics (MD) simulations in conjunction with the SPC/E water model, we optimize ionic force-field parameters for seven different halide and alkali ions, considering a total of eight ion-pairs. Our strategy is based on simultaneous optimizing single-ion and ion-pair properties, i.e., we first fix ion water parameters based on single-ion solvation free energies, and in a second step determine the cation-anion interaction parameters (traditionally given by mixing or combination rules) based on the Kirkwood-Buff theory without modification of the ion-water interaction parameters. In doing so, we have introduced scaling factors for the cation-anion Lennard-Jones (LJ) interaction that quantify deviations from the standard mixing rules. For the rather size-symmetric salt solutions involving bromide and chloride ions, the standard mixing rules work fine. On the other hand, for the iodide and fluoride solutions, corresponding to the largest and smallest anion considered in this work, a rescaling of the mixing rules was necessary. For iodide, the experimental activities suggest more tightly bound ion pairing than given by the standard mixing rules, which is achieved in simulations by reducing the scaling factor of the cation-anion LJ energy. For fluoride, the situation is different and the simulations show too large attraction between fluoride and cations when compared with experimental data. For NaF, the situation can be rectified by increasing the cation-anion LJ energy. For KF, it proves necessary to increase the effective cation-anion Lennard-Jones diameter. The optimization strategy outlined in this work can be easily adapted to different kinds of ions

Probing the sp^2 dependence of elastic moduli in ultrahard diamond films

The structural and elastic properties of diamond nanocomposites and ultrananocrystalline diamond films (UNCD) are investigated using both empirical potentials and tight binding schemes. We find that both materials are extremely hard, but their superb diamondlike properties are limited by their sp^2 component. In diamond composites, the sp^2 atoms are found in the matrix and far from the interface with the inclusion, and they are responsible for the softening of the material. In UNCD, the sp^2 atoms are located in the grain boundaries. They offer relaxation mechanisms which relieve the strain but, on the other hand, impose deformations that lead to softening. The higher the sp^2 component the less rigid these materials are.Comment: 10 pages, 3 figures. to appear in Diamond and Relarted Material