Mechanism of enhancement of ferroelectricity of croconic acid with temperature.

A detailed study of the thermal behaviour of atomic motions in the organic ferroelectric croconic acid is presented in the temperature range 5-300 K. Using high-resolution inelastic neutron scattering and first-principles electronic-structure calculations within the framework of density functional theory and a quasiharmonic phonon description of the material, we find that the frequencies of the well defined doublet in inelastic neutron scattering spectra associated with out-of-plane motions of hydrogen-bonded protons decrease monotonically with temperature indicating weakening of these bonding motifs and enhancement of proton motions. Theoretical mean-square displacements for these proton motions are within 5% of experimental values. A detailed analysis of this observable shows that it is unlikely that there is a facile proton transfer along the direction of ferroelectric polarization in the absence of an applied electric field. Calculations predict constrained thermal motion of proton along crystallographic lattice direction c retaining the hydrogen bond motif of the crystal at high temperature. Using the Berry-phase method, we have also calculated the spontaneous polarization of temperature dependent cell structures, and find that our computational model provides a satisfactory description of the anomalous and so far unexplained rise in bulk electric polarization with temperature. Correlating the thermal motion induced lattice strain with temperature dependent spontaneous polarizations, we conclude that increasing thermal strain with temperatures combined with constrained thermal motion along the hydrogen bond motif are responsible of this increase in ferroelectricity at high temperature.

Requirements of first-principles calculations of X-ray absorption spectra of liquid water.

A computational benchmark study on X-ray absorption spectra of water has been performed by means of transition-potential density functional theory (TP-DFT), damped time-dependent density functional theory (TDDFT), and damped coupled cluster (CC) linear response theory. For liquid water, using TDDFT with a tailored CAM-B3LYP functional and a polarizable embedding, we find that an embedding with over 2000 water molecules is required to fully converge spectral features for individual molecules, but a substantially smaller embedding can be used within averaging schemes. TP-DFT and TDDFT calculations on 100 MD structures demonstrate that TDDFT produces a spectrum with spectral features in good agreement with experiment, while it is more difficult to fully resolve the spectral features in the TP-DFT spectrum. Similar trends were also observed for calculations of bulk ice. In order to further establish the performance of these methods, small water clusters have been considered also at the CC2 and CCSD levels of theory. Issues regarding the basis set requirements for spectrum simulations of liquid water and the determination of gas-phase ionization potentials are also discussed.

Molecular reordering processes on ice (0001) surfaces from long timescale simulations.

We report results of long timescale adaptive kinetic Monte Carlo simulations aimed at identifying possible molecular reordering processes on both proton-disordered and ordered (Fletcher) basal plane (0001) surfaces of hexagonal ice. The simulations are based on a force field for flexible molecules and span a time interval of up to 50 µs at a temperature of 100 K, which represents a lower bound to the temperature range of earth's atmosphere. Additional calculations using both density functional theory and an ab initio based polarizable potential function are performed to test and refine the force field predictions. Several distinct processes are found to occur readily even at this low temperature, including concerted reorientation (flipping) of neighboring surface molecules, which changes the pattern of dangling H-atoms, and the formation of interstitial defects by the downwards motion of upper-bilayer molecules. On the proton-disordered surface, one major surface roughening process is observed that significantly disrupts the crystalline structure. Despite much longer simulation time, such roughening processes are not observed on the highly ordered Fletcher surface which is energetically more stable because of smaller repulsive interaction between neighboring dangling H-atoms. However, a more localized process takes place on the Fletcher surface involving a surface molecule transiently leaving its lattice site. The flipping process provides a facile pathway of increasing proton-order and stabilizing the surface, supporting a predominantly Fletcher-like ordering of low-temperature ice surfaces. Our simulations also show that eventual proton-disordered patches on the surface may induce significant local reconstructions. Further, a subset of the molecules on the Fletcher surface are susceptible to forming interstitial defects which might provide active sites for various chemical reactions in the atmosphere.

Complementarity between high-energy photoelectron and L-edge spectroscopy for probing the electronic structure of 5d transition metal catalysts.

We demonstrate the successful use of hard X-ray photoelectron spectroscopy (HAXPES) for selectively probing the platinum partial d-density of states (DOS) in a Pt-Cu nanoparticle catalyst which shows activity superior to pure Pt towards the oxygen-reduction reaction (ORR). The information about occupied Pt d-band states was complemented by Pt L(2)-edge X-ray absorption near-edge spectroscopy (XANES), which probes unoccupied valence states. We found a significant electronic perturbation of the Pt projected d-DOS which was narrowed and shifted to higher binding energy compared to pure platinum. The effect of this electronic structure perturbation on the chemical properties of the nanoparticle surface is discussed in terms of the d-band model. We have thereby demonstrated that the combination of L-edge spectroscopy and HAXPES allows for an experimental derivation of the valence electronic structure in an element-specific way for 5d metal catalysts.

On the range of water structure models compatible with X-ray and neutron diffraction data.

We use the reverse Monte Carlo (RMC) method to critically evaluate the structural information content of diffraction data on bulk water by fitting simultaneously or separately to X-ray and neutron data; the O-H and H-H, but not the O-O, pair-correlation functions (PCFs) are well-described by the neutron data alone. Enforcing at the same time different H-bonding constraints, we generate four topologically different structure models of liquid water, including a simple mixture model, that all equally well reproduce the diffraction data. Although earlier work [Leetmaa, M.; et al. J. Chem. Phys. 2008, 129, 084502] has focused on tetrahedrality in the H-bond network in liquid water, we show here that, even for the O-O-O three-body correlation, tetrahedrality is not strictly defined by the data. We analyze how well two popular MD models (TIP4P-pol2 and SPC/E) reproduce the neutron data in q-space and find differences in important aspects from the experiment. From the RMC fits, we obtain pair-correlation functions (PCFs) that are in optimal agreement with the diffraction data but still show a surprisingly strong variability both in position and height of the first intermolecular (H-bonding) O-H peak. We conclude that, although diffraction data impose important constraints on the range of possible water structures, additional data are needed to narrow the range of possible structure models.