Marivi Fernandez-Serra

Network equilibration and first-principles liquid water

Motivated by the very low diffusivity recently found in ab initio simulations of liquid water, we have studied its dependence with temperature, system size, and duration of the simulations. We use ab initio molecular dynamics (AIMD), following the Born-Oppenheimer forces obtained from density-functional theory (DFT). The linear-scaling capability of our method allows the consideration of larger system sizes (up to 128 molecules in this study), even if the main emphasis of this work is in the time scale. We obtain diffusivities that are substantially lower than the experimental values, in agreement with recent findings using similar methods. A fairly good agreement with D(T) experiments is obtained if the simulation temperature is scaled down by approximately 20%. It is still an open question whether the deviation is due to the limited accuracy of present density functionals or to quantum fluctuations, but neither technical approximations (basis set, localization for linear scaling) nor the system size (down to 32 molecules) deteriorate the DFT description in an appreciable way. We find that the need for long equilibration times is consequence of the slow process of rearranging the H-bond network (at least 20 ps at AIMDs room temperature). The diffusivity is observed to be very directly linked to network imperfection. This link does not appear an artifact of the simulations, but a genuine property of liquid water.

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.

Electrons and Hydrogen-Bond Connectivity in Liquid Water

The network connectivity in liquid water is revised in terms of electronic signatures of hydrogen bonds (HBs) instead of geometric criteria, in view of recent x-ray absorption studies. The analysis is based on ab initio molecular-dynamics simulations at ambient conditions. Even if instantaneous threadlike structures are observed in the electronic network, they continuously reshape in oscillations reminiscent of the and modes in ice (tau approximately 170 fs). However, two water molecules initially joined by a HB remain effectively bound over many periods regardless of its electronic signature.

Direct detection of sub-GeV dark matter with semiconductor targets

A bstractDark matter in the sub-GeV mass range is a theoretically motivated but largely unexplored paradigm. Such light masses are out of reach for conventional nuclear recoil direct detection experiments, but may be detected through the small ionization signals caused by dark matter-electron scattering. Semiconductors are well-studied and are particularly promising target materials because their O

Two exchange-correlation functionals compared for first-principles liquid water

\mathcal{O}

GW and Bethe-Salpeter study of small water clusters

(1 eV) band gaps allow for ionization signals from dark matter particles as light as a few hundred keV. Current direct detection technologies are being adapted for dark matter-electron scattering. In this paper, we provide the theoretical calculations for dark matter-electron scattering rate in semiconductors, overcoming several complications that stem from the many-body nature of the problem. We use density functional theory to numerically calculate the rates for dark matter-electron scattering in silicon and germanium, and estimate the sensitivity for upcoming experiments such as DAMIC and SuperCDMS. We find that the reach for these upcoming experiments has the potential to be orders of magnitude beyond current direct detection constraints and that sub-GeV dark matter has a sizable modulation signal. We also give the first direct detection limits on sub-GeV dark matter from its scattering off electrons in a semiconductor target (silicon) based on published results from DAMIC. We make available publicly our code, QEdark, with which we calculate our results. Our results can be used by experimental collaborations to calculate their own sensitivities based on their specific setup. The searches we propose will probe vast new regions of unexplored dark matter model and parameter space.

Ferroelectric PbTiO3/SrRuO3 Superlattices with Broken Inversion Symmetry

The electron density in the peptide bonds of crambin, a plant-seed hydrophobic globular protein with 46 residues and 642 atoms, is studied both theoretically and experimentally. The results of density functional calculations of crambin in vacuo for the deformation electron density in its peptide bonds are compared with the electronic distribution obtained from ultra-high-resolution X-ray crystallography. The comparison is centered on the average peptide-bond map, where the experimental results are clearest. Theory is then used to ascertain on differences among peptide bonds in different chemical environments.

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

The first-principles description of liquid water using ab initio molecular dynamics (AIMD) based on density functional theory (DFT) has recently been found to require long equilibration times, giving too low diffusivities and a clear over-structuring of the liquid. In the light of these findings we compare here the room-temperature description offered by two different exchange correlation functionals: BLYP, the most popular for liquid water so far, and RPBE, a revision of the widely used PBE. We find for RPBE a less structured liquid with radial distribution functions closer to experiment ones than the ones of BLYP. The diffusivity obtained with RPBE for heavy water is still 20% lower than the corresponding experimental value, but it represents a substantial improvement on the BLYP value, one order of magnitude lower than experiment. These characteristics and the hydrogen-bond (HB) network imperfection point to an effective temperature ∼3% lower than the actual simulation temperature for the RPBE liquid, as compared with BLYPs ∼17% deviation. The too long O-O average nearest-neighbor distance observed points to an excessively weak HB, possibly compensating more fundamental errors in the DFT description.