Christopher D. Daub

Water-mediated ordering of nanoparticles in an electric field

Interfacial polar molecules feature a strongly anisotropic response to applied electric field, favoring dipole orientations parallel to the interface. In water, in particular, this effect combines with generic orientational preferences induced by spatial asymmetry of water hydrogen bonding under confined geometry, which may give rise to a Janus interface. The two effects manifest themselves in considerable dependence of water polarization on both the field direction relative to the interface and the polarity (sign) of the field. Using molecular simulations, we demonstrate strong field-induced orientational forces acting on apolar surfaces through water mediation. At a field strength comparable to electric fields around a DNA polyion, the torques we predict to act on an adjacent nanoparticle are sufficient to overcome thermal fluctuations. These torques can align a particle with surface as small as 1 nm2. The mechanism can support electrically controlled ordering of suspended nanoparticles as a means of tuning their properties and can find application in electro-nanomechanical devices.

The influence of molecular-scale roughness on the surface spreading of an aqueous nanodrop

We examine the effect of nanoscale roughness on spreading and surface mobility of water nanodroplets. Using molecular dynamics, we consider model surfaces with sub-nanoscale asperities at varied surface coverage and with different distribution patterns. We test materials that are hydrophobic, and those that are hydrophilic in the absence of surface corrugations. Interestingly, on both types of surfaces, the introduction of surface asperities gives rise to a sharp increase in the apparent contact angle. The Cassie-Baxter equation is obeyed approximately on hydrophobic substrates, however, the increase in the contact angle on a hydrophilic surface differs qualitatively from the behavior on macroscopically rough surfaces described by the Wenzel equation. On the hydrophobic substrate, the superhydrophobic state with the maximal contact angle of 180 degrees is reached when the asperity coverage falls below 25%, suggesting that superhydrophobicity can also be achieved by the nanoscale roughness of a macroscopically smooth material. We further examine the effect of surface roughness on droplet mobility on the substrate. The apparent diffusion constant shows a dramatic slow down of the nanodroplet translation even for asperity coverage in the range of 1% for a hydrophilic surface, while droplets on corrugated hydrophobic surfaces retain the ability to flow around the asperities. In contrast, for smooth surfaces we find that the drop mobility on the hydrophilic surface exceeds that on the hydrophobic one.

Nanoscale Wetting Under Electric Field from Molecular Simulations

Applying an electric field is a well-established experimental method to tune surface wettability. As accessible experimental length scales become shorter, the modification of interfacial properties of water using electric field must come to grips with novel effects existing at the nanoscale. We survey recent progress in understanding these effects on water interfacial tension and on water-mediated interactions using molecular simulations. We highlight the key role of external conditions in determining the systems response to applied electric field. We further discuss the role of appropriate boundary conditions in modeling polar fluids subject to collective polarization. The work reviewed here broadens the basic understanding of applied and internal field effects that can operate in condensed phase systems, from modulating local hydrophilicity/hydrophobicity of engineered and biological surfaces, to surface manipulation in nanofluidic devices.

The role of electron correlation on calculated XH-stretching vibrational band intensities

XH-stretching vibrational band intensities have been calculated for seven small molecules with OH, NH and CH bonds, respectively, using the simple harmonically coupled anharmonic oscillator local m...

Structure of Aqueous Solutions of Monosodium Glutamate

We studied monosodium glutamate (MSG) in aqueous solution using molecular dynamics simulations and compared the results with recent neutron diffraction with isotope contrast variation/empirical potential structure refinement (EPSR) data obtained on the same system (McLain et al. J. Phys. Chem. B 2006, 110, 21251-21258). We used classical simulations with empirical force fields to study both dilute and more concentrated (1.40 M) solutions. To gauge the importance of polarization and other quantum effects, we carried out first-principles molecular dynamics in the dilute case. The glutamate structure was well reproduced by the OPLS/AA and SPC/E force fields: we found a reasonable agreement between the simulations and the experimental data with respect to the hydration numbers for glutamate carboxylate and amine groups and the observation of significant sodium ion-carboxylate binding. However, none of our simulations could reproduce the dramatic reduction in water-water correlations observed experimentally. Simulations showed a large amount of carboxylate-amine binding, as well as segregation of water and glutamate, at moderately high concentrations of MSG. We attribute this result to the breakdown of currently available classical force fields when applied to concentrated ionic solutions, especially large polyatomic ions. We also did not observe the sharing of a water proton by two carboxylate oxygens simultaneously, and we argue against this interpretation of EPSR data on a variety of physical grounds. We offer several suggestions to resolve these discrepancies between simulation and the current interpretation of neutron diffraction data, which should advance the understanding of aqueous ionic solutions in general.

Note: How does the treatment of electrostatic interactions influence the magnitude of thermal polarization of water? The SPC/E model

We investigate how the treatment of electrostatic interactions influences the magnitude of the thermal polarization of water. We performed non-equilibrium molecular dynamics simulations of the extended simple point charge model of water under a thermal gradient, using two different systems: a water droplet confined in a spherical wall where the interactions are computed exactly using the Coulombic potential and a periodic prismatic box using the Wolf and 3D Ewald methods. All the methods reproduce the thermal polarization (TP) of water as well as the direction of the TP field, but the standard implementation of the Wolf method overestimates the strength of the TP field by one order of magnitude, showing that this method might be problematic in simulations involving temperature and/or density gradients.

Lithium ion-water clusters in strong electric fields: a quantum chemical study.

We use density functional theory to investigate the impact that strong electric fields have on the structure and energetics of small lithium ion-water clusters, Li(+)·nH2O, with n = 4 or 6. We find that electric field strengths of ∼0.5 V/Å are sufficient to break the symmetry of the n = 4 tetrahedral energy minimum structure, which undergoes a transformation to an asymmetric cluster consisting of three water molecules bound to lithium and one additional molecule in the second solvation shell. Interestingly, this cluster remains the global minimum configuration at field strengths ≳0.15 V/Å. The 6-coordinated cluster, Li(+)·6H2O, features a similar transition to 5- and 4-coordinated clusters at field strengths of ∼0.2 and ∼0.3 V/Å, respectively, with the tetra-coordinated structure being the global minimum even in the absence of the field. Our findings are relevant to understanding the behavior of the Li(+) ion in aqueous environments under strong electric fields and in interfacial regions where field gradients are significant.

Local Field Factors and Dielectric Properties of Liquid Benzene.

Local electric field factors are calculated for liquid benzene by combining molecular dynamic simulations with a subsequent force-field model based on a combined charge-transfer and point-dipole interaction model for the local field factor. The local field factor is obtained as a linear response of the local field to an external electric field, and the response is calculated at frequencies through the first absorption maximum. It is found that the largest static local field factor is around 2.4, while it is around 6.4 at the absorption frequency. The linear susceptibility, the dielectric constant, and the first absorption maximum of liquid benzene are also studied. The electronic contribution to the dielectric constant is around 2.3 at zero frequency, in good agreement with the experimental value around 2.2, while it increases to 6.3 at the absorption frequency. The π → π* excitation energy is around 6.0 eV, as compared to the gas-phase value of around 6.3 eV, while the experimental values are 6.5 and 6.9 eV for the liquid and gas phase, respectively, demonstrating that the gas-to-liquid shift is well-described.

Thermo-molecular orientation effects in fluids of dipolar dumbbells

We use molecular dynamics simulations in applied thermal gradients to study thermomolecular orientation (TMO) of size-asymmetric dipolar dumbbells with different molecular dipole moments. We find that the direction of the TMO is the same as in apolar dumbbells of the same size, i.e. the smaller atom in the dumbbell tends to orient towards the colder temperature. The ratio of the electrical polarization to the magnitude of the thermal gradient does not vary much with the magnitude of the molecular dipole moment. We also investigate a novel second order TMO that persists even in size-symmetric dipolar dumbbells where molecules have a slight tendency to orient perpendicular to the gradient except very close to the hot region, where (anti-)parallel orientations are preferred. Finally, we investigate rotational correlation functions and characteristic rotational times in these systems in an attempt to model possible spectroscopic signatures of TMO in experiments. Although we cannot detect any difference in integrated rotational times between equilibrium simulations and simulations in a thermal gradient, more careful modelling of the anisotropic rotational dynamics in the thermal gradient may be more successful.

Molecular dynamics simulations to examine structure, energetics, and evaporation/condensation dynamics in small charged clusters of water or methanol containing a single monatomic ion.

We study small clusters of water or methanol containing a single Ca(2+), Na(+), or Cl(-) ion with classical molecular dynamics simulations, using models that incorporate polarizability via the Drude oscillator framework. Evaporation and condensation of solvent from these clusters is examined in two systems, (1) for isolated clusters initially prepared at different temperatures and (2) those with a surrounding inert (Ar) gas of varying temperature. We examine these clusters over a range of sizes, from almost bare ions up to 40 solvent molecules. We report data on the evaporation and condensation of solvent from the clusters and argue that the observed temperature dependence of evaporation in the smallest clusters demonstrates that the presence of heated gas alone cannot, in most cases, solely account for bare ion production in electrospray ionization (ESI), neglecting the key contribution of the electric field. We also present our findings on the structure and energetics of the clusters as a function of size. Our data agree well with the abundant literature on hydrated ion clusters and offer some novel insight into the structure of methanol and ion clusters, especially those with a Cl(-) anion, where we observe the presence of chain-like structures of methanol molecules. Finally, we provide some data on the reparameterizations necessary to simulate ions in methanol using the separately developed Drude oscillator models for methanol and for ions in water.