Density functional calculations for polymers and clusters - progress and limitations

The applicability of density functional (DF) methods has progressed greatly since the first workshop of this series ten years ago. Applications that show both the successes and the limitations can be found in the fields of: (a) the structures of the isomers of atomic clusters. and (b) the structure of organic molecules and polymers, and their reactions with additional molecules. We shall review some of the results and the lessons to be learned from them. (C) 2001 Elsevier Science B.V. All rights reserved.

Computational challenges in atomic, molecular and optical physics

Six challenges are discussed. These are the laser-driven helium atom; the laser-driven hydrogen molecule and hydrogen molecular ion: electron scattering (with ionization) from one-electron atoms; the vibrational and rotational structure of molecules such as H-3(+) and water at their dissociation limits; laser- heated clusters; and quantum degeneracy and Bose-Einstein condensation. The first four concern fundamental few-body systems where use of high-performance computing (HPC) is currently making possible accurate modelling from first principles. This leads to reliable predictions and support for laboratory experiment as well as true understanding of the dynamics. Important aspects of these challenges addressable only via a terascale facility are set out. Such a facility makes the last two challenges in the above list meaningfully accessible for the first time, and the scientific interest together with the prospective role for HPC in these is emphasized.

Hydrogen bonding and collective proton modes in clusters and periodic layers of squaric acid: A density functional study

Hydrogen bonding in clusters and extended layers of squaric acid molecules has been investigated by density functional computations. Equilibrium geometries, harmonic vibrational frequencies, and energy barriers for proton transfer along hydrogen bonds have been determined using the Car-Parrinello method. The results provide crucial parameters for a first principles modeling of the potential energy surface, and highlight the role of collective modes in the low-energy proton dynamics. The importance of quantum effects in condensed squaric acid systems has been investigated, and shown to be negligible for the lowest-energy collective proton modes. This information provides a quantitative basis for improved atomistic models of the order-disorder and displacive transitions undergone by squaric acid crystals as a function of temperature and pressure. (C) 2001 American Institute of Physics.

An isothermal-isobaric Langevin thermostat for simulating nanoparticles under pressure: Application to Au clusters

We present a method for simulating clusters or, molecules subjected to an external pressure, which is exerted by a pressure-transmitting medium. It is based on the canoninical Langevin thermostat, but extended in such a way that the Brownian forces are allowed to operate only from the region exterior to the cluster. We show that the frictional force of the Langevin thermostat is linked to the pressure of the reservoir in a unique way, and that this property manifests itself when the particle it acts upon is not pointlike but has finite dimensions. By choosing appropriately the strength of the random forces and the friction coefficient, both temperature and pressure can be controlled independently. We illustrate the capabilities of this new method by calculating the compressibility of small gold clusters under pressure.

Double ionization of atomic negative ions in an intense laser field

In recent years there have been many studies of multiple ionization of closed shell rare gas atoms by intense laser fields. Until now no similar work has been done in the study of more diverse targets such as negative ions where low binding energies and strong electron correlations could yield distinctive behaviour. We present the first results of ionization of more than one electron from a range of atomic negative ions by intense laser pulses. Although these pulses are long by modern standards, and tend to produce sequential ionization in atoms, the positive ion yields from the negative ions do not depend predictably on the ionization potentials. This suggests that there may, intriguingly, be an alternative mechanism enhancing double ionization at low intensities.

SrTiO3(001)(2x1) reconstructions: First-principles calculations of surface energy and atomic structure compared with scanning tunneling microscopy images

(1x1) and (2x1) reconstructions of the (001) SrTiO3 surface were studied using the first-principles full-potential linear muffin-tin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The (1x1) unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the (2x1) Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the (1x1) surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the (2x1) reconstructions would be easier to image than the (1x1) surfaces. The possibility that the presence of oxygen vacancies would destabilise the (1x1) surfaces was also investigated. If the (1x1) surfaces are unstable then there exists the further possibility that the (2x1) DL-TiO2 reconstruction [Natasha Erdman Nature (London) 419, 55 (2002)] is stable in a TiO2-rich environment and for p(O2)>10(-18) atm.

Bismuth embrittlement of copper is an atomic size effect

Embrittlement by the segregation of impurity elements to grain boundaries is one of a small number of phenomena that can lead to metallurgical failure by fast fracture(1). Here we settle a question that has been debated for over a hundred years(2): how can minute traces of bismuth in copper cause this ductile metal to fail in a brittle manner? Three hypotheses for Bi embrittlement of Cu exist: two assign an electronic effect to either a strengthening(3) or weakening(4) of bonds, the third postulates a simple atomic size effect(5). Here we report first principles quantum mechanical calculations that allow us to reject the electronic hypotheses, while supporting a size effect. We show that upon segregation to the grain boundary, the large Bi atoms weaken the interatomic bonding by pushing apart the Cu atoms at the interface. The resolution of the mechanism underlying grain boundary weakening should be relevant for all cases of embrittlement by oversize impurities.

The VLT-FLAMES survey of massive stars: Observations in the Galactic clusters NGC 3293, NGC 4755 and NGC 6611

We introduce a new survey of massive stars in the Galaxy and the Magellanic Clouds using the Fibre Large Array Multi- Element Spectrograph ( FLAMES) instrument at the Very Large Telescope ( VLT). Here we present observations of 269 Galactic stars with the FLAMES- Giraffe Spectrograph ( R similar or equal to 25 000), in fields centered on the open clusters NGC3293, NGC4755 and NGC6611. These data are supplemented by a further 50 targets observed with the Fibre- Fed Extended Range Optical Spectrograph ( FEROS, R = 48 000). Following a description of our scientific motivations and target selection criteria, the data reduction methods are described; of critical importance the FLAMES reduction pipeline is found to yield spectra that are in excellent agreement with less automated methods. Spectral classifications and radial velocity measurements are presented for each star, with particular attention paid to morphological peculiarities and evidence of binarity. These observations represent a significant increase in the known spectral content of NGC3293 and NGC4755, and will serve as standards against which our subsequent FLAMES observations in the Magellanic Clouds will be compared.

Positronium scattering by atomic hydrogen with inclusion of target excitation channels

Calculations are reported for positronium (Ps) scattering by atomic hydrogen (H) in the energy range 0-6.5 eV in a coupled- pseudostate approximation in which excitation and ionization channels of both the Ps and the H are taken into account. The approximation contains an accurate representation of the van der Waals coefficient. Results are presented for phase shifts, scattering lengths, effective ranges, and various cross sections including partial wave, total, and ortho-para conversion cross sections. An analysis of the possible spin transitions is provided and the energy of the positronium hydride (PsH) bound state is determined. Substantial differences are found from earlier work within the frozen target approximation, now clearly confirming the importance of target excitation channels. Good agreement is obtained with recent calculations of S-wave phase shifts and scattering lengths using the stabilization method. Convergence to the exact binding energy for PsH appears to be slow. Resonances corresponding to unstable states of the positron orbiting H- are seen in the electronic spin singlet partial waves. The importance of the H- formation channel is discussed.

Atomic collisions involving positrons

In this short review we look at bound states, positron-atom scattering. positronium-atom scattering. positronium-positronium scattering, cold antihydrogen and annihilation. (c) 2005 Published by Elsevier B.V.

Ionization of atomic hydrogen and He+ by slow antiprotons

We study the ionization of H(1s), He+(1s) and He+(2s) by antiprotons in the energy range from 0.1 to 500 keV. We adopt a semiclassical single centre close-coupling approach in which the wavefunction for the electron is expanded in a B-spline basis centred on the nucleus of the atom/ion. Comparison is made with existing theoretical calculations and available experimental data. The results are encouraging.

An investigation into the use of atomic datasets in simulations of the Ni-like gadolinium x-ray laser

The effect of differing the datasets used in the modelling of the Ni-like Gd x-ray laser (XRL) is examined through the 1.50 hydro-atomic code, EHYBRID. Two atomic datasets, including energy levels and radiative and collisional excitation rates, are used as input data for the code. It is found that the behaviour of the XRL is somewhat different than might be expected from superficial examination of the atomic data. The similarities in the gain profiles at low densities are found to have encouraging implications. in our attempts to model XRLs.