Dielectronic recombination in He-like, Li-like, and Be-like highly charged ions in the KLL and KLM manifolds

This paper reports a systematic study of the dependence on atomic number of the dielectronic recombination resonance strengths for He-like, Li-like and Be-like ions. Recent measurements of dielectronic recombination resonance strengths for the KLL and KLM manifolds for iron, yttrium, iodine, holmium, and bismuth are also described. The resonance strengths were normalized to calculated electron impact ionization cross sections. The measured resonance strengths generally agree well with theoretical calculations using the distorted wave approximation. However, KLM resonance strength measurements on high atomic number open-shell ions gave higher values than those suggested by calculations. Using recently measured data, along with existing results, scaling laws have been generated as a function of atomic number for He-like, Li-like, and Be-like ions in the KLL and KLM manifolds.

Electron shell contributions to gamma-ray spectra of positron annihilation in noble gases

Gamma-ray positron annihilation spectra of the noble gases are simulated using computational chemistry tools for the bound electron wavefunctions and plane-wave approximation for the low-energy positron. The present annihilation line shapes, i.e. the full width at half maximum, Delta epsilon, of the gamma-ray annihilation spectra for He and Ar (valence) agree well with available independent atomic calculations using a different algorithm. For other noble gases they achieve moderate agreement with the experimental measurements. It is found that the contributions of various atomic electron shells to the spectra depend significantly on their principal quantum number n and orbital angular momentum quantum number l. The present study further reveals that the outermost ns electrons of the noble gases exhibit spectral line shapes in close agreement with those measured, indicating (as expected) that the measurements are not due to a simple sum over the momentum densities for all atomic electrons. The robust nature of the present approach makes it possible for us to proceed to more complex molecular systems using the tools of modern computational chemistry.

Modelling the mechanical response of the reed-mouthpiece-lip system of a clarinet. Part II: A lumped model approximation

A non-linear lumped model of the reed-mouthpiece-lip system of a clarinet is formulated, in which the lumped parameters are derived from numerical experiments with a finite-difference simulation based on a distributed reed model. The effective stiffness per unit area is formulated as a function of the pressure signal driving the reed, in order to simulate the effects of the reed bending against the lay, and mass and damping terms are added as a first approximation to the dynamic behaviour of the reed. A discrete-time formulation is presented, and its response is compared to that of the distributed model. In addition, the lumped model is applied in the simulation of clarinet tones, enabling the analysis of the effects of using a pressure-dependent stiffness per unit area on sustained oscillations. The analysed effects and features are in qualitative agreement with players' experiences and experimental results obtained in prior studies.

MEASUREMENT AND CALCULATION OF ABSOLUTE SINGLE AND MULTIPLE CHARGE EXCHANGE CROSS SECTIONS FOR Feq+ IONS IMPACTING H2O

Charge exchange (CE) plays a fundamental role in the collisions of solar- and stellar-wind ions with lunar and planetary exospheres, comets, and circumstellar clouds. Reported herein are absolute cross sections for single, double, triple, and quadruple CE of Feq+ (q = 5-13) ions with H2O at a collision energy of 7q keV. One measured value of the pentuple CE is also given for Fe9+ ions. An electron cyclotron resonance ion source is used to provide currents of the highly charged Fe ions. Absolute data are derived from knowledge of the target gas pressure, target path length, and incident and charge-exchanged ion currents. Experimental cross sections are compared with new results of the n-electron classical trajectory Monte Carlo approximation. The radiative and non-radiative cascades following electron transfers are approximated using scaled hydrogenic transition probabilities and scaled Auger rates. Also given are estimates of cross sections for single capture, and multiple capture followed by autoionization, as derived from the extended overbarrier model. These estimates are based on new theoretical calculations of the vertical ionization potentials of H2O up to H2O10+.

Many-body theory calculations of positron binding to negative ions

A many-body theory approach developed by the authors [Phys. Rev. A 70, 032720 (2004)] is applied to positron bound states and annihilation rates in atomic systems. Within the formalism, full account of virtual positronium (Ps) formation is made by summing the electron-positron ladder diagram series, thus enabling the theory to include all important many-body correlation effects in the positron problem. Numerical calculations have been performed for positron bound states with the hydrogen and halogen negative ions, also known as Ps hydride and Ps halides. The Ps binding energies of 1.118, 2.718, 2.245, 1.873 and 1.393 eV and annihilation rates of 2.544, 2.482, 1.984, 1.913 and 1.809 ns^{-1}, have been obtained for PsH, PsF, PsCl, PsBr and PsI, respectively.

Large-scale density functional calculations of the surface properties of the Wigner crystal

The surface properties of the jellium model have been investigated by large supercell computations in the density functional theory-local spin-density (DFT-LSD) approach for planar slabs with up to 1000 electrons. A wide interval of densities has been explored, extending into the stability range of the Wigner crystal. Most computations have been carried out on nominally paramagnetic samples with an equal number of spin-up and spin-down electrons. The results show that within DFT-LSD spontaneous spin polarization and charge localization start nearly simultaneously at the surface for r(s) similar to 20, then, with decreasing density, they progress toward the center of the slab. Electrons are fully localized and spin polarized at r(s) = 30. At this density the charge distribution is the superposition of disjoint charge blobs, each corresponding to one electron. The distribution of blobs displays both regularities and disorder, the first being represented by well-defined planes and simple in-plane geometries, and the latter by a variety of surface defects. The surface energy, surface dipole, electric polarisability, and magnetization pattern have been determined as a function of density. All these quantities display characteristic anomalies at the density of the localization transition. The analysis of the low-frequency electric conductivity shows that in the fluid paramagnetic regime the in-plane current preferentially flows in the central region of the slab and the two spin channels are equally conducting. In the charge localized, spin-polarized regime, conductivity is primarily a surface effect, and an apparent asymmetry is observed in the two spin currents.

Strong-field approximation for Coulomb explosion of H-2(+) and D-2(+) by short intense laser pulses

We present a simple quantum mechanical model to describe Coulomb explosion of H-2(+) and D-2(+) by short, intense infrared laser pulses. The model is based on the length gauge version of the molecular strong-field approximation and is valid when the process of dissociation prior to ionization is negligible. The results are compared with recent experimental data for the proton kinetic energy spectrum [Th. Ergler , Phys. Rev. Lett. 95, 093001 (2005); D. S. Murphy , J. Phys. B 40, S359 (2007)]. Using a Franck-Condon distribution over initial vibrational states, the theory reproduces the overall shape of the spectrum with only a small overestimation of slow protons. The agreement between theory and experiment can be made perfect by using a non-Frank-Condon initial distribution characteristic for H-2(+) (D-2(+)) targets produced by strong-field ionization of H-2 (D-2). For comparison, we also present results obtained by two different tunneling models for this process.

An understanding of chemoselective hydrogenation on crotonaldehyde over Pt(111) in the free energy landscape: The microkinetics study based on first-principles calculations

Microkinetic model is developed in the free energy landscape based on density functional theory (DFT) to quantitatively investigate the reaction mechanism of chemoselective partial hydrogenation of crotonaldehyde to crotyl alcohol over Pt(1 1 1) at the temperature of 353 K. Three different methods (mobile, immobile and collision theory models) were carried out to obtain free energy barrier of adsorption/desorption processes. The results from mobile and collision theory models are similar. The calculated TOFs from both models are close to the experiment value. However, for the immobile model, in which the free energy barrier of desorption approaches the energy barrier, the calculated TOF is 2 orders of magnitude lower than the other models. The difficulty of adsorption/ desorption may be overestimated in the immobile model. In addition, detailed analyses show that for the surface hydrogenation elementary steps, the entropy and internal energy effects are small under the reaction condition, while the zero-point-energy (ZPE) correction is significant, especially for the multi-step hydrogenation reaction. The total energy with the ZPE correction approaches to the full free energy calculation for the surface reaction under the reaction condition. (c) 2011 Elsevier B.V. All rights reserved.

Approximating the Set of Local Minima in Partial RNA Folding Landscapes

Motivation: We study a stochastic method for approximating the set of local minima in partial RNA folding landscapes associated with a bounded-distance neighbourhood of folding conformations. The conformations are limited to RNA secondary structures without pseudoknots. The method aims at exploring partial energy landscapes pL induced by folding simulations and their underlying neighbourhood relations. It combines an approximation of the number of local optima devised by Garnier and Kallel (2002) with a run-time estimation for identifying sets of local optima established by Reeves and Eremeev (2004).

Results: The method is tested on nine sequences of length between 50 nt and 400 nt, which allows us to compare the results with data generated by RNAsubopt and subsequent barrier tree calculations. On the nine sequences, the method captures on average 92% of local minima with settings designed for a target of 95%. The run-time of the heuristic can be estimated by O(n2D?ln?), where n is the sequence length, ? is the number of local minima in the partial landscape pL under consideration and D is the maximum number of steepest descent steps in attraction basins associated with pL.