Picosecond metrology of laser-driven proton bursts.

Tracking primary radiation-induced processes in matter requires ultrafast sources and high precision timing. While compact laser-driven ion accelerators are seeding the development of novel high instantaneous flux applications, combining the ultrashort ion and laser pulse durations with their inherent synchronicity to trace the real-time evolution of initial damage events has yet to be realized. Here we report on the absolute measurement of proton bursts as short as 3.5±0.7 ps from laser solid target interactions for this purpose. Our results verify that laser-driven ion acceleration can deliver interaction times over a factor of hundred shorter than those of state-of-the-art accelerators optimized for high instantaneous flux. Furthermore, these observations draw ion interaction physics into the field of ultrafast science, opening the opportunity for quantitative comparison with both numerical modelling and the adjacent fields of ultrafast electron and photon interactions in matter.

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules.

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water.

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice--leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point.

Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water.

We demonstrate a model for stoichiometric and reduced titanium dioxide intended for use in molecular dynamics and other atomistic simulations and based in the polarizable ion tight binding theory. This extends the model introduced in two previous papers from molecular and liquid applications into the solid state, thus completing the task of providing a comprehensive and unified scheme for studying chemical reactions, particularly aimed at problems in catalysis and electrochemistry. As before, experimental results are given priority over theoretical ones in selecting targets for model fitting, for which we used crystal parameters and band gaps of titania bulk polymorphs, rutile and anatase. The model is applied to six low index titania surfaces, with and without oxygen vacancies and adsorbed water molecules, both in dissociated and non-dissociated states. Finally, we present the results of molecular dynamics simulation of an anatase cluster with a number of adsorbed water molecules and discuss the role of edge and corner atoms of the cluster.

Irradiation of water ice by C(+) ions in the cosmic environment.

We present a first-principles MD (FPMD) study of the interaction of low-energy, positively charged carbon (C(+)) projectiles with amorphous solid water clusters at 30 K. Reactions involving the carbon ion at an initial energy of 11 and 1.7 eV with a 30-molecule cluster have been investigated. Simulations indicate that the neutral isoformyl radical, COH(•), and carbon monoxide, CO, are the dominant products of these reactions. All of these reactions are accompanied by the transfer of a proton from the reacting water molecule to the ice, where it forms a hydronium ion. We find that COH(•) is formed either via a direct, "knock-out", mechanism following the impact of the C(+) projectile upon a water molecule or by creation of a COH2(+) intermediate. The direct mechanism is more prominent at higher energies. CO is generally produced following the dissociation of COH(•). More frequent production of the formyl radical, HCO(•), is observed here than in gas-phase calculations. A less commonly occurring product is the dihydroxymethyl, CH(OH)2(•), radical. Although a minor result, its existence gives an indication of the increasing chemical complexity that is possible in such heterogeneous environments.

Electron-induced hydrogen loss in uracil in a water cluster environment.

Low-energy electron-impact hydrogen loss due to dissociative electron attachment (DEA) to the uracil and thymine molecules in a water cluster environment is investigated theoretically. Only the A(')-resonance contribution, describing the near-threshold behavior of DEA, is incorporated. Calculations are based on the nonlocal complex potential theory and the multiple scattering theory, and are performed for a model target with basic properties of uracil and thymine, surrounded by five water molecules. The DEA cross section is strongly enhanced when the attaching molecule is embedded in a water cluster. This growth is due to two effects: the increase of the resonance lifetime and the negative shift in the resonance position due to interaction of the intermediate negative ion with the surrounding water molecules. A similar effect was earlier found in DEA to chlorofluorocarbons.

Organic synthesis in the interstellar medium by low-energy carbon irradiation.

We present a first principles molecular dynamics (FPMD) study of the interaction of low-energy neutral carbon projectiles with amorphous solid water clusters at 30 K. Reactions involving the carbon atom at an initial energy of 11 and 1.7 eV with 30-molecule clusters have been investigated. Simulations indicate that the formation of hydroxymethylene, an intermediate in formaldehyde production, dominates at the higher energy. The reaction proceeds by fragmenting a water molecule, binding the carbon to the OH radical, and saturating the C valence with a hydrogen atom that can arise from the originally dissociated water molecule, or through a chain of proton transfer events. We identified several possible pathways for the formation of HCOH. When the initial collision occurs at the periphery of the cluster, we observe the formation of CO and the evaporation of water molecules. At the lower energy water fragmentation is not favorable, thus leading to the formation of weakly bound carbon-water complexes.

Electronic stopping power in gold: the role of d electrons and the H/He anomaly.

The electronic stopping power of H and He moving through gold is obtained to high accuracy using time-evolving density-functional theory, thereby bringing usual first principles accuracies into this kind of strongly coupled, continuum nonadiabatic processes in condensed matter. The two key unexplained features of what observed experimentally have been reproduced and understood: (i) The nonlinear behavior of stopping power versus velocity is a gradual crossover as excitations tail into the d-electron spectrum; and (ii) the low-velocity H/He anomaly (the relative stopping powers are contrary to established theory) is explained by the substantial involvement of the d electrons in the screening of the projectile even at the lowest velocities where the energy loss is generated by s-like electron-hole pair formation only.

A tight binding model for water.

We demonstrate for the first time a tight binding model for water incorporating polarizable oxygen atoms. A novel aspect is that we adopt a "ground up" approach in that properties of the monomer and dimer only are fitted. Subsequently we make predictions of the structure and properties of hexamer clusters, ice-XI and liquid water. A particular feature, missing in current tight binding and semiempirical hamiltonians, is that we reproduce the almost two-fold increase in molecular dipole moment as clusters are built up toward the limit of bulk liquid. We concentrate on properties of liquid water, particularly dielectric constant and self diffusion coefficient, which are very well rendered in comparison with experiment. Finally we comment on the question of the contrasting densities of water and ice which is central to an understanding of the subtleties of the hydrogen bond.

Evidence of short-range screening in shock-compressed aluminum plasma.

We have investigated the angular variation in elastic x-ray scattering from a dense, laser-shock-compressed aluminum foil. A comparison of the experiment with simulations using an embedded atom potential in a molecular dynamics simulation shows a significantly better agreement than simulations based on an unscreened one-component plasma model. These data illustrate, experimentally, the importance of screening for the dense plasma static structure factor.

Compton scatter profiles for warm dense matter.

In this paper, we discuss the possibility of using x-ray Compton scattering as a probe of the outer electronic structure of ions immersed in warm dense matter. It is proposed that the x-ray free-electron lasers currently under construction will provide an ideal tool for this, with the main pulse being used to create a uniform well-defined sample and the third harmonic providing a clean monochromatic probe. We model the plasma photon scatter spectrum by combining self-consistent finite-temperature electronic structure calculations with molecular dynamics simulations of the ion-ion structure factor. In particular, we present bound-free Compton profiles that are more accurate that those obtained using form factor or impulse approximations.

Competition of charge-density waves and superconductivity in sulfur.

A one-dimensional charge-density wave (CDW) instability is shown to be responsible for the formation of the incommensurate modulation of the atomic lattice in the high-pressure phase of sulfur. The coexistence of, and competition between, the CDW and the superconducting state leads to the previously observed increase of T{c} up to 17 K, which we attribute to the suppression of the CDW instability, the same phenomenology found in doped layered dichalcogenides.

Neutral and charged 1-butyl-3-methylimidazolium triflate clusters: equilibrium concentration in the vapor phase and thermal properties of nanometric droplets.

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P approximately 0.01 Pa) and temperatures (T >or= 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0

Simulation of interfaces between room temperature ionic liquids and other liquids.

The structure and properties of the interfaces between the room temperature ionic liquid dimethylimidazolium chloride ([dmim]Cl) and different Lennard-Jones fluids and between ionic liquid and water have been studied by molecular dynamics simulations, and compared to the ionic liquid-vapour interface. Two contrasting types of interface were investigated, thermodynamically stable interfaces between ionic liquid and vapour and between ionic liquid and Lennard-Jones fluids, and diffusing interfaces between miscible phases of different compositions involving water. The density profiles of different species through the interface are presented. The cations and water molecules near the former type of interface are aligned relative to the surface, but no orientational preference was found near or in the broad diffusing interface. The ionic liquid has a negative electrostatic potential relative to vapour or Lennard-Jones fluid, but is more positive than pure water. This contrast is explained in terms of the relative importance of orientation and concentration differences in the two types of interface.

Ferroelectricity and isotope effects in hydrogen-bonded KDP crystals.

Based on an accurate first principles description of the energetics in H-bonded potassium-dihydrogen-phosphate crystals, we conduct a first study of nuclear quantum effects and of the changes brought about by deuteration. Tunneling is allowed only for clusters involving correlated protons and heavy ion displacements, the main effect of deuteration being a depletion of the proton probability density at the O-H-O bridge center, which in turn weakens its proton-mediated covalent bonding. The ensuing lattice expansion couples self-consistently with the proton off-centering, thus explaining both the giant isotope effect and its close connection with geometrical effects.