Local order of liquid water at metallic electrode surfaces.

We study the structure and dynamics of liquid water in contact with Pd and Au (111) surfaces using ab initio molecular dynamics simulations with and without van der Waals interactions. Our results show that the structure of water at the interface of these two metals is very different. For Pd, we observe the formation of two different domains of preferred orientations, with opposite net interfacial dipoles. One of these two domains has a large degree of in-plane hexagonal order. For Au, a single domain exists with no in-plane order. For both metals, the structure of liquid water at the interface is strongly dependent on the use of dispersion forces. The origin of the structural domains observed in Pd is associated to the interplay between water/water and water/metal interactions. This effect is strongly dependent on the charge transfer that occurs at the interface and which is not modeled by current state of the art semi-empirical force fields.

Polar nanoregions in water: a study of the dielectric properties of TIP4P/2005, TIP4P/2005f and TTM3F.

We present a critical comparison of the dielectric properties of three models of water-TIP4P/2005, TIP4P/2005f, and TTM3F. Dipole spatial correlation is measured using the distance dependent Kirkwood function along with one-dimensional and two-dimensional dipole correlation functions. We find that the introduction of flexibility alone does not significantly affect dipole correlation and only affects ɛ(ω) at high frequencies. By contrast the introduction of polarizability increases dipole correlation and yields a more accurate ɛ(ω). Additionally, the introduction of polarizability creates temperature dependence in the dipole moment even at fixed density, yielding a more accurate value for dɛ/dT compared to non-polarizable models. To better understand the physical origin of the dielectric properties of water we make analogies to the physics of polar nanoregions in relaxor ferroelectric materials. We show that ɛ(ω, T) and τD(T) for water have striking similarities with relaxor ferroelectrics, a class of materials characterized by large frequency dispersion in ɛ(ω, T), Vogel-Fulcher-Tammann behaviour in τD(T), and the existence of polar nanoregions.

Room temperature compressibility and diffusivity of liquid water from first principles.

The isothermal compressibility of water is essential to understand its anomalous properties. We compute it by ab initio molecular dynamics simulations of 200 molecules at five densities, using two different van der Waals density functionals. While both functionals predict compressibilities within ~30% of experiment, only one of them accurately reproduces, within the uncertainty of the simulation, the density dependence of the self-diffusion coefficient in the anomalous region. The discrepancies between the two functionals are explained in terms of the low- and high-density structures of the liquid.

Ferroelectric PbTiO3/SrRuO3 superlattices with broken inversion symmetry.

We have fabricated PbTiO3/SrRuO3 superlattices with ultrathin SrRuO3 layers. Because of the superlattice geometry, the samples show a large anisotropy in their electrical resistivity, which can be controlled by changing the thickness of the PbTiO3 layers. Therefore, along the ferroelectric direction, SrRuO3 layers can act as dielectric, rather than metallic, elements. We show that, by reducing the concentration of PbTiO3, an increasingly important effect of polarization asymmetry due to compositional inversion symmetry breaking occurs. The results are significant as they represent a new class of ferroelectric superlattices, with a rich and complex phase diagram. By expanding our set of materials we are able to introduce new behaviors that can only occur when one of the materials is not a perovskite titanate. Here, compositional inversion symmetry breaking in bicolor superlattices, due to the combined variation of A and B site ions within the superlattice, is demonstrated using a combination of experimental measurements and first principles density functional theory.

Anomalous nuclear quantum effects in ice.

One striking anomaly of water ice has been largely neglected and never explained. Replacing hydrogen (1H) by deuterium (2H) causes ice to expand, whereas the normal isotope effect is volume contraction with increased mass. Furthermore, the anomaly increases with temperature T, even though a normal isotope shift should decrease with T and vanish when T is high enough to use classical nuclear motions. In this study, we show that these effects are very well described by ab initio density-functional theory. Our theoretical modeling explains these anomalies, and allows us to predict and to experimentally confirm a counter effect, namely, that replacement of 16O by 18O causes a normal lattice contraction.

Quasi-harmonic approximation of thermodynamic properties of ice Ih, II, and III.

Several thermodynamic properties of ice Ih, II, and III are studied by a quasi-harmonic approximation and compared to results of quantum path integral and classical simulations. This approximation allows to obtain thermodynamic information at a fraction of the computational cost of standard simulation methods, and at the same time permits studying quantum effects related to zero-point vibrations of the atoms. Specifically, we have studied the crystal volume, bulk modulus, kinetic energy, enthalpy, and heat capacity of the three ice phases as a function of temperature and pressure. The flexible q-TIP4P/F model of water was employed for this study, although the results concerning the capability of the quasi-harmonic approximation are expected to be valid independently of the employed water model. The quasi-harmonic approximation reproduces with reasonable accuracy the results of quantum and classical simulations showing an improved agreement at low temperatures (T< 100 K). This agreement does not deteriorate as a function of pressure as long as it is not too close to the limit of mechanical stability of the ice phases.

Anthropomorphic Polydimethylsiloxane silicone-based phantom for Diffuse Optical Imaging.

This work presents a method for constructing phantoms suitable for diffuse optical mammography. They are based on Polydimethylsiloxane silicones, with the characteristic of being anthropomorphic, and having similar mechanical and optical properties as a real breast. These phantoms are useful for testing the performance of diffuse optical imaging devices in the near infrared, both in transmittance and reflectance geometries, since they can be constructed containing inclusions, to simulate breast tumors. An alternative component to be used as scattering agent, that is easier to handle than traditional scattering agents, is also studied. The optical properties of the phantoms were tested varying the concentration of scattering and absorbing agents, while their mechanical properties were modified by adding a silicone fluid to the basic mixture. Finally, the phantoms were tested by Diffuse Optical Imaging experiments, and these images were compared to the ones obtained by conventional ultrasound techniques. Results show that the constructed anthropomorphic phantoms properly reproduce the optical and mechanical characteristics of human breasts, and are suitable to be used in Diffuse Optical Imaging.

The role of hydrogen bonding in water-metal interactions.

The hydrogen bond interaction between water molecules adsorbed on a Pd <111> surface, a nucleator of two dimensional ordered water arrays at low temperatures, is studied using density functional theory calculations. The role of the exchange and correlation density functional in the characterization of both the hydrogen bond and the water-metal interaction is analyzed in detail. The effect of non local correlations using the van der Waals density functional proposed by Dion et al. [M. Dion, H. Rydberg, E. Schröder, D. C. Langreth and B. I. Lundqvist, Phys. Rev. Lett., 2004, 92, 246401] is also studied. We conclude that the choice of this potential is critical in determining the cohesive energy of water-metal complexes. We show that the interaction between water molecules and the metal surface is as sensitive to the density functional choice as hydrogen bonds between water molecules are. The reason for this is that the two interactions are very similar in nature. We make a detailed analogy between the water-water bond in the water dimer and the water-Pd bond at the Pd <111> surface. Our results show a strong similarity between these two interactions and based on this we describe the water-Pd bond as a hydrogen bond type interaction. These results demonstrate the need to obtain an accurate and reliable representation of the hydrogen bond interaction in density functional theory.

Density, structure, and dynamics of water: the effect of van der Waals interactions.

It is known that ab initio molecular dynamics (AIMD) simulations of liquid water at ambient conditions, based on the generalized gradient approximation (GGA) to density functional theory (DFT), with commonly used functionals fail to produce structural and diffusive properties in reasonable agreement with experiment. This is true for canonical, constant temperature simulations where the density of the liquid is fixed to the experimental density. The equilibrium density, at ambient conditions, of DFT water has recently been shown by Schmidt et al. [J. Phys. Chem. B, 113, 11959 (2009)] to be underestimated by different GGA functionals for exchange and correlation, and corrected by the addition of interatomic pair potentials to describe van der Waals (vdW) interactions. In this contribution we present a DFT-AIMD study of liquid water using several GGA functionals as well as the van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)]. As expected, we find that the density of water is grossly underestimated by GGA functionals. When a vdW-DF is used, the density improves drastically and the experimental diffusivity is reproduced without the need of thermal corrections. We analyze the origin of the density differences between all the functionals. We show that the vdW-DF increases the population of non-H-bonded interstitial sites, at distances between the first and second coordination shells. However, it excessively weakens the H-bond network, collapsing the second coordination shell. This structural problem is partially associated to the choice of GGA exchange in the vdW-DF. We show that a different choice for the exchange functional is enough to achieve an overall improvement both in structure and diffusivity.

Insights into the Structure of Liquid Water from Nuclear Quantum Effects on the Density and Compressibility of Ice Polymorphs.

Nuclear quantum effects lead to an anomalous shift of the volume of hexagonal ice; heavy ice has a larger volume than light ice. Furthermore, this anomaly in ice increases with temperature and persists in liquid water up to the boiling point. To gain more insight, we study nuclear quantum effects on the density and compressibility of several ice-like structures and crystalline ice phases. By calculating the anisotropic contributions to the stain tensor, we analyze how the compressibility changes along different directions in hexagonal ice, and find that hexagonal ice is softer along the x- y plane than the z-direction. Furthermore, by performing ab initio density functional theory calculations with a van der Waals functional and with the quasiharmonic approximation, we find an anomalous isotope effect in the bulk modulus of hexagonal ice: heavy ice has a smaller bulk modulus than light ice. In agreement with the experiments, we also obtain an anomalous isotope effect for clathrate hydrate structure I. For the rest of the ice polymorphs, the isotope effect is (i) anomalous for ice IX, Ih, Ic, clathrate, and low density liquid-like (LDL-like) amorphous ice; (ii) normal at T = 0 K and becomes anomalous with increasing temperature for ice IX, II, high density liquid-like (HDL-like) amorphous ices, and ice XV; and (iii) normal for ice VIII up to the melting point. There is a transition from an anomalous isotope effect to a normal isotope effect for both the volume and bulk modulus, as the density (compressibility) of the structures increases (decreases). This result can explain the anomalous isotope effect in liquid water: as the compressibility decreases from the melting point to the compressibility minimum temperature, the difference between the volumes of the heavy and light water rapidly decreases, but the effect stays anomalous up to the boiling temperature as the hydrogen bond network is never completely broken by fully filling all of the interstitial sites.

Bias-dependent local structure of water molecules at a metallic interface.

Understanding the local structure of water at the interfaces of metallic electrodes is a key issue in aqueous-based electrochemistry. Nevertheless a realistic simulation of such a setup is challenging, particularly when the electrodes are maintained at different potentials. To correctly compute the effect of an external bias potential applied to truly semi-infinite surfaces, we combine Density Functional Theory (DFT) and Non-Equilibrium Green's Function (NEGF) methods. This framework allows for the out-of-equilibrium calculation of forces and dynamics, and directly correlates to the chemical potential of the electrodes, which is introduced experimentally. In this work, we apply this methodology to study the electronic properties and atomic forces of a water molecule at the interface of a gold surface. We find that the water molecule tends to align its dipole moment with the electric field, and it is either repelled or attracted to the metal depending on the sign and magnitude of the applied bias, in an asymmetric fashion.

Continuous melting through a hexatic phase in confined bilayer water.

Liquid water is not only of obvious importance but also extremely intriguing, displaying many anomalies that still challenge our understanding of such an a priori simple system. The same is true when looking at nanoconfined water: The liquid between constituents in a cell is confined to such dimensions, and there is already evidence that such water can behave very differently from its bulk counterpart. A striking finding has been reported from computer simulations for two-dimensionally confined water: The liquid displays continuous or discontinuous melting depending on its density. In order to understand this behavior, we have analyzed the melting exhibited by a bilayer of nanoconfined water by means of molecular dynamics simulations. At high density we observe the continuous melting to be related to the phase change of the oxygens only, with the hydrogens remaining liquidlike throughout. Moreover, we find an intermediate hexatic phase for the oxygens between the liquid and a triangular solid ice phase, following the Kosterlitz-Thouless-Halperin-Nelson-Young theory for two-dimensional melting. The liquid itself tends to maintain the local structure of the triangular ice, with its two layers being strongly correlated yet with very slow exchange of matter. The decoupling in the behavior of the oxygens and hydrogens gives rise to a regime in which the complexity of water seems to disappear, resulting in what resembles a simple monoatomic liquid. This intrinsic tendency of our simulated water may be useful for understanding novel behaviors in other confined and interfacial water systems.

Pathways for flip-flop of mono- and di-anionic phospholipids in the erythrocyte membrane.

The inward translocations (flip), from the outer to the inner membrane leaflet of human erythrocytes, of di-anionic NBD-labeled phospholipids containing as a head group phosphate esters of glycolate, butyrate and hydroxyethanesulfonate are slow processes (k = 0.005-0.008 h-1, 37 degrees C) at pH 7.4. A decrease of pH highly stimulates the flip. A major role of the anion exchanger (AE1), band 3, in this flip is indicated by (a) the strong inhibition of the flip (55-85%) by stilbene disulfonates and other inhibitors of anion transport, (b) the stimulation and loss of pH dependence of the flip after modification of band 3 by Woodward's reagent K and NaBH4, and (c) the stimulation of the flip after proteolytic cleavage of band 3 by papain. The flip of mono-anionic NBD-phospholipids with phosphate esters of glycerol, glycol, methanol, butanol and benzyl alcohol is much faster than that of their dianionic analogs (k = 0.04 to > 3.0 h-1, 37 degrees C). It is inhibited by stilbene disulfonates to a decreasing extent (35 to 0%) and is not affected by several reversible inhibitors of anion exchange. This indicates a minor component of band-3-mediated flip and a major component of nonmediated flip. The outward translocations (flop), from the inner to outer membrane leaflet, of both mono- and di-anionic phospholipids are very fast (1.0-5.9 h-1), ATP-dependent and inhibitable by vanadate, fluoride, SH-reagents or Mg(2+)-depletion of cells and thereby likely to be largely mediated by a 'floppase'. The stationary distributions of the NBD-labeled anionic phospholipids are asymmetric to an extent (outer to inner leaflet ratio 2-9) correlating with the ratio of the rates of the outward and the inward translocation. Thus, asymmetry is largely abolished by blockage of the floppase-mediated translocation.

In vitro interaction of mutagenic chromium (VI) with red blood cells.

The interaction of mutagenic Cr(VI) with red blood cells has been studied by ESR spectroscopy. Signals of two Cr(V) species are observed almost immediately after contacting red cells with chromate(VI) aqueous solution at pH 7.4. The signal at go = 1.985, which decays within one hour, is attributed to a Cr(V) complex formed by glutathione due its reducing and chelating ability. The other signal at go = 1.979, which is distinctly more persistent, may indicate that some immobilization of the formed Cr(V) ions takes place on the macromolecular cell components, e.g. glycoproteins.

Enhanced IgG- and complement-independent phagocytosis of sulfatide-enriched human erythrocytes by human monocytes.

Phagocytosis by adherent human monocytes of human erythrocytes (RBC), sulfatide-enriched by incubation with 10(-12) to 10(-9) M cerebroside sulfate, was enhanced approx. 6-fold. Increased phagocytosis was observed only in RBC opsonized with fresh plasma, and not in non-opsonized or serum-opsonized RBC. Increased phagocytosis was immunoglobulin- and complement independent. Thrombospondin and von Willebrand factor, present in plasma but not in serum, and binding selectively to sulfatides, are likely mediators of the enhanced phagocytosis.

ESR investigation of the structural transition of spectrin in ghosts and in solution.

The effect of temperature on the human erythrocyte ghosts and on isolated spectrin, both labelled with a protein spin label (2, 2, 6, 6 tetramethylpiperidin - 1 - oxyl - 4 - maleimide), was studied by the electron spin resonance spin label method. The results show that the spectra associated with the membrane undergo significant changes with temperature and a well-defined transition is detected at 52 degrees C. The observed transition is removed from the membrane by the extraction of spectrin and in the supernatant the transition is unchanged from that exhibited by the native ghosts. The results show that it is possible to observe the heat dependent transition of the spectrin in whole ghosts avoiding the extraction procedure.

Diffuse reflectance optical topography: location of inclusions in 3D and detectability limits.

In the present contribution we investigate the images of CW diffusely reflected light for a point-like source, registered by a CCD camera imaging a turbid medium containing an absorbing lesion. We show that detection of µa variations (absorption anomalies) is achieved if images are normalized to background intensity. A theoretical analysis based on the diffusion approximation is presented to investigate the sensitivity and the limitations of our proposal and a novel procedure to find the location of the inclusions in 3D is given and tested. An analysis of the noise and its influence on the detection capabilities of our proposal is provided. Experimental results on phantoms are also given, supporting the proposed approach.

Optimal finite-range atomic basis sets for liquid water and ice.

Finite-range numerical atomic orbitals are the basis functions of choice for several first principles methods, due to their flexibility and scalability. Generating and testing such basis sets, however, remains a significant challenge for the end user. We discuss these issues and present a new scheme for generating improved polarization orbitals of finite range. We then develop a series of high-accuracy basis sets for the water molecule, and report on their performance in describing the monomer and dimer, two phases of ice, and liquid water at ambient and high density. The tests are performed by comparison with plane-wave calculations, and show the atomic orbital basis sets to exhibit an excellent level of transferability and consistency. The highest-order bases (quadruple-ζ) are shown to give accuracies comparable to a plane-wave kinetic energy cutoff of at least ~1000 eV for quantities such as energy differences and ionic forces, as well as achieving significantly greater accuracies for total energies and absolute pressures.

A first-principles study of the effect of charge doping on the 1D polymerization of C(60).

We study the interplay between charge doping and intermolecular distance in the polymerization of C(60) fullerene chains by means of density functional theory-based first-principles calculations. The potential energy surface analysis shows that both the equilibrium intermolecular distance of the unpolymerized system and the polymerization energy barrier are inversely proportional to the electron doping of the system. We analyze the origin of this charge-induced polymerization effect by studying the behavior of the system's wavefunctions around the Fermi level and the structural modifications of the molecules as a function of two variables: the distance between the centers of the molecules and the number of electrons added to the system.