Scattering theory of cooling and heating in optomechanical systems

We present a one-dimensional scattering theory which enables us to describe a wealth of effects arising from the coupling of the motional degree of freedom of scatterers to the electromagnetic field. Multiple scattering to all orders is taken into account. The theory is applied to describe the scheme of a Fabry-Perot resonator with one of its mirrors moving. The friction force, as well as the diffusion, acting on the moving mirror is derived. In the limit of a small reflection coefficient, the same model provides for the description of the mechanical effect of light on an atom moving in front of a mirror.

Deployment on GPUs of an application in computational atomic physics

This paper describes the deployment on GPUs of PROP, a program of the 2DRMP suite which models electron collisions with H-like atoms and ions. Because performance on GPUs is better in single precision than in double precision, the numerical stability of the PROP program in single precision has been studied. The numerical quality of PROP results computed in single precision and their impact on the next program of the 2DRMP suite has been analyzed. Successive versions of the PROP program on GPUs have been developed in order to improve its performance. Particular attention has been paid to the optimization of data transfers and of linear algebra operations. Performance obtained on several architectures (including NVIDIA Fermi) are presented.

Time-dependent R-matrix theory applied to two-photon double ionization of He

We introduce a time-dependent R-matrix theory generalized to describe double-ionization processes. The method is used to investigate two-photon double ionization of He by intense XUV laser radiation. We combine a detailed B-spline-based wave-function description in an extended inner region with a single-electron outer region containing channels representing both single ionization and double ionization. A comparison of wave-function densities for different box sizes demonstrates that the flow between the two regions is described with excellent accuracy. The obtained two-photon double-ionization cross sections are in excellent agreement with other cross sections available. Compared to calculations fully contained within a finite inner region, the present calculations can be propagated over the time it takes the slowest electron to reach the boundary.

Theory of phonon dissipation in the conduction of stressed Au nanowires

Recent experiments on Au break junctions [Phys. Rev. Lett. 88 (2002) 216803] have characterized the nonlinear conductance of stretched short Au nanowires. They reveal in the voltage range 10-20 meV the signatures of dissipation effects, likely due to phonons in the nanowire, reducing the conductance below the quantized value of 2e(2)/h. We present here a theory, based on a model tight-binding Hamiltonian and on non-equilibrium Green's function techniques, which accounts for the main features of the experiment. The theory helps in revealing details of the experiment which need to be addressed with a more realistic, less idealized, theoretical framework. (C) 2004 Elsevier B.V. All rights reserved.

Harmonic generation by noble-gas atoms in the near-IR regime using ab initio time-dependent R-matrix theory

We demonstrate the capability of ab initio time-dependent R-matrix theory to obtain accurate harmonic generation spectra of noble-gas atoms at near-IR wavelengths between 1200 and 1800 nm and peak intensities up to 1.8 × 10^(14) W/cm^(2). To accommodate the excursion length of the ejected electron, we use an angular-momentum expansion up to Lmax=279. The harmonic spectra show evidence of atomic structure through the presence of a Cooper minimum in harmonic generation for Kr, and of multielectron interaction through the giant resonance for Xe. The theoretical spectra agree well with those obtained experimentally.

Identifying the distinct features of geometric structures for hole trapping to generate radicals on rutile TiO2(110) in photooxidation using density functional theory calculations with hybrid functional

Using density functional theory calculations with HSE 06 functional, we obtained the structures of spin-polarized radicals on rutile TiO2(110), which is crucial to understand the photooxidation at the atomic level, and further calculate the thermodynamic stabilities of these radicals. By analyzing the results, we identify the structural features for hole trapping in the system, and reveal the mutual effects among the geometric structures, the energy levels of trapped hole states and their hole trapping capacities. Furthermore, the results from HSE 06 functional are compared to those from DFT + U and the stability trend of radicals against the number of slabs is tested. The effect of trapped holes on two important steps of the oxygen evolution reaction, i.e. water dissociation and the oxygen removal, is investigated and discussed.

Mn ion dissolution from MnS:a density functional theory study

The dissolution of MnS inclusions could induce pitting corrosion in stainless steels, but its dissolution mechanism is poorly understood at the atomic scale. With the help of ab initio molecular dynamics calculations, one inevitable step in the dissolution of MnS is studied by simulating the process of one Mn ion leaving the surface. The reaction mechanism is determined to contain three steps with two large barriers and a small one, leading to two slow steps in the Mn ion dissolution. Comparing to the Na ion dissolution from NaCl, the barriers of the Mn ion dissolution are much larger, which is a reflection of their different electronic structures.

A density functional theory study of CH2 and H adsorption on Ni(111)

Ab initio total energy calculations within the density functional theory framework have been used to study the adsorption of CH2 and H as well as the coadsorption of CH2 and H on Ni(111). H binds strongly at threefold hollow sites with calculated adsorption energies of 2.60 and 2.54 eV at the face-centered-cubic (fcc) and hexagonal-close-packed (hcp) hollow sites, respectively. Adsorption energies and H-Ni distances are found to agree well with both experimental and theoretical results. CH2 adsorbs strongly at all high symmetry sites with calculated adsorption energies of 3.26, 3.22, 3.14 and 2.36 eV at the fcc, hcp, bridge and top sites, respectively. Optimized structures are reported at all sites, and, in the most stable hollow sites there is considerable internal reorganization of the CH2 fragment. The CH2 molecule is tilted, the hydrogens are inequivalent and the C-H bonds are lengthened relative to the gas phase. In the CH2-H coadsorption systems the adsorbates have a tendency to move toward bridge sites. The bonding of all adsorbates to the surface is analyzed in detail. (C) 2000 American Institute of Physics. [S0021-9606(00)71213-X].

Double ionization in R-matrix theory using a two-electron outer region

We have developed a two-electron outer region for use within R-matrix theory to describe double ionisation processes. The capability of this method is demonstrated for single-photon double ionisation of He in the photon energy region between 80 eV to 180 eV. The cross sections are in agreement with established data. The extended RMT method also provides information on higher-order processes, as demonstrated by the identification of signatures for sequential double ionisation processes involving an intermediate He⁺ state with n=2.

Electron-impact ionization and recombination of M -shell atomic ions in the argon isonuclear sequence

Electron-impact ionization and recombination cross sections and rate coefficients are calculated for M-shell Ar atomic ions using a configuration-average distorted-wave method. The electron-impact ionization calcula- tions are for all atomic ions in the Ar isonuclear sequence. Ionization contributions include both direct ioniza- tion and excitation-autoionization processes. Good agreement is found between theory and experimental crossed-beam measurements for moderately charged ion stages. Comparisons are made with previous theoret- ical calculations where possible.We also generate rate coefficients for neutral argon ionization, based on recent R-matrix with pseudostates calculations. Electron-impact dielectronic recombination is calculated for all M-shell ions of argon. For Ar6+ and Ar7+ the current theoretical results agree well with previous level-resolved distorted-wave calculations. In order to compare with published ionization balance results our dielectronic recombination data are combined with literature values for the higher ion stages and with recent radiative recombination data for all the ion stages. We find significant differences in our equilibrium fractional abun- dances for the M-shell ions, compared with literature values. We relate these differences to the underlying atomic data.

Single-photon single ionization of W+ ions: experiment and theory

Experimental and theoretical results are reported for photoionization of Ta-like (W+) tungsten ions. Absolute cross sections were measured in the energy range 16–245 eV employing the photon–ion merged-beam setup at the advanced light source in Berkeley. Detailed photon-energy scans at 100 meV bandwidth were performed in the 16–108 eV range. In addition, the cross section was scanned at 50 meV resolution in regions where fine resonance structures could be observed. Theoretical results were obtained from a Dirac–Coulomb R-matrix approach. Photoionization cross section calculations were performed for singly ionized atomic tungsten ions in their 5s²5p⁶5d⁴(⁵D)6s ⁶Dj.  J = 1/2, ground level and the associated excited metastable levels with J = 3/2, 5/2, 7/2 and 9/2. Since the ion beams used in the experiments must be expected to contain long-lived excited states also from excited configurations, additional cross-section calculations were performed for the second-lowest term, 5d⁵⁶Sj, J = 5/2, and for the 4F term, 5d³6s^{2 4}Fj, with J = 3/2, 5/2, 7/2 and 9/2. Given the complexity of the electronic structure of W+ the calculations reproduce the main features of the experimental cross section quite well.

Application of [Lambda] to the fourth perturbation theory in calculating the equation of state of rare gas solids and fcc metals

We have calculated the thermodynamic properties of monatomic fcc crystals from the high temperature limit of the Helmholtz free energy. This equation of state included the static and vibrational energy components. The latter contribution was calculated to order A4 of perturbation theory, for a range of crystal volumes, in which a nearest neighbour central force model was used. We have calculated the lattice constant, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the adiabatic and the isothermal bulk modulus, and the Gruneisen parameter, for two of the rare gas solids, Xe and Kr, and for the fcc metals Cu, Ag, Au, Al, and Pb. The LennardJones and the Morse potential were each used to represent the atomic interactions for the rare gas solids, and only the Morse potential was used for the fcc metals. The thermodynamic properties obtained from the A4 equation of state with the Lennard-Jones potential, seem to be in reasonable agreement with experiment for temperatures up to about threequarters of the melting temperature. However, for the higher temperatures, the results are less than satisfactory. For Xe and Kr, the thermodynamic properties calculated from the A2 equation of state with the Morse potential, are qualitatively similar to the A 2 results obtained with the Lennard-Jones potential, however, the properties obtained from the A4 equation of state are in good agreement with experiment, since the contribution from the A4 terms seem to be small. The lattice contribution to the thermal properties of the fcc metals was calculated from the A4 equation of state, and these results produced a slight improvement over the properties calculated from the A2 equation of state. In order to compare the calculated specific heats and bulk moduli results with experiment~ the electronic contribution to thermal properties was taken into account~ by using the free electron model. We found that the results varied significantly with the value chosen for the number of free electrons per atom.

On the equation of state and atomic mean square displacement of crystals

We have presented a Green's function method for the calculation of the atomic mean square displacement (MSD) for an anharmonic Hamil toni an . This method effectively sums a whole class of anharmonic contributions to MSD in the perturbation expansion in the high temperature limit. Using this formalism we have calculated the MSD for a nearest neighbour fcc Lennard Jones solid. The results show an improvement over the lowest order perturbation theory results, the difference with Monte Carlo calculations at temperatures close to melting is reduced from 11% to 3%. We also calculated the MSD for the Alkali metals Nat K/ Cs where a sixth neighbour interaction potential derived from the pseudopotential theory was employed in the calculations. The MSD by this method increases by 2.5% to 3.5% over the respective perturbation theory results. The MSD was calculated for Aluminum where different pseudopotential functions and a phenomenological Morse potential were used. The results show that the pseudopotentials provide better agreement with experimental data than the Morse potential. An excellent agreement with experiment over the whole temperature range is achieved with the Harrison modified point-ion pseudopotential with Hubbard-Sham screening function. We have calculated the thermodynamic properties of solid Kr by minimizing the total energy consisting of static and vibrational components, employing different schemes: The quasiharmonic theory (QH), ).2 and).4 perturbation theory, all terms up to 0 ().4) of the improved self consistent phonon theory (ISC), the ring diagrams up to o ().4) (RING), the iteration scheme (ITER) derived from the Greens's function method and a scheme consisting of ITER plus the remaining contributions of 0 ().4) which are not included in ITER which we call E(FULL). We have calculated the lattice constant, the volume expansion, the isothermal and adiabatic bulk modulus, the specific heat at constant volume and at constant pressure, and the Gruneisen parameter from two different potential functions: Lennard-Jones and Aziz. The Aziz potential gives generally a better agreement with experimental data than the LJ potential for the QH, ).2, ).4 and E(FULL) schemes. When only a partial sum of the).4 diagrams is used in the calculations (e.g. RING and ISC) the LJ results are in better agreement with experiment. The iteration scheme brings a definitive improvement over the).2 PT for both potentials.

Phonon dispersion curves and atomic mean square displacement for several fcc and bcc materials

The atomic mean square displacement (MSD) and the phonon dispersion curves (PDC's) of a number of face-centred cubic (fcc) and body-centred cubic (bcc) materials have been calclllated from the quasiharmonic (QH) theory, the lowest order (A2 ) perturbation theory (PT) and a recently proposed Green's function (GF) method by Shukla and Hiibschle. The latter method includes certain anharmonic effects to all orders of anharmonicity. In order to determine the effect of the range of the interatomic interaction upon the anharmonic contributions to the MSD we have carried out our calculations for a Lennard-Jones (L-J) solid in the nearest-neighbour (NN) and next-nearest neighbour (NNN) approximations. These results can be presented in dimensionless units but if the NN and NNN results are to be compared with each other they must be converted to that of a real solid. When this is done for Xe, the QH MSD for the NN and NNN approximations are found to differ from each other by about 2%. For the A2 and GF results this difference amounts to 8% and 7% respectively. For the NN case we have also compared our PT results, which have been calculated exactly, with PT results calculated using a frequency-shift approximation. We conclude that this frequency-shift approximation is a poor approximation. We have calculated the MSD of five alkali metals, five bcc transition metals and seven fcc transition metals. The model potentials we have used include the Morse, modified Morse, and Rydberg potentials. In general the results obtained from the Green's function method are in the best agreement with experiment. However, this improvement is mostly qualitative and the values of MSD calculated from the Green's function method are not in much better agreement with the experimental data than those calculated from the QH theory. We have calculated the phonon dispersion curves (PDC's) of Na and Cu, using the 4 parameter modified Morse potential. In the case of Na, our results for the PDC's are in poor agreement with experiment. In the case of eu, the agreement between the tlleory and experiment is much better and in addition the results for the PDC's calclliated from the GF method are in better agreement with experiment that those obtained from the QH theory.

Systematic construction of models for the exchange hole of density functional theory.

Dans ce travail, nous étendons le nombre de conditions physiques actuellement con- nues du trou d’échange exact avec la dérivation de l’expansion de quatrième ordre du trou d’échange sphérique moyenne exacte. Nous comparons les expansions de deux- ième et de quatrième ordre avec le trou d’échange exact pour des systèmes atomiques et moléculaires. Nous avons constaté que, en général, l’expansion du quatrième ordre reproduit plus fidèlement le trou d’échange exact pour les petites valeurs de la distance interélectronique. Nous démontrons que les ensembles de base de type gaussiennes ont une influence significative sur les termes de cette nouvelle condition, en étudiant com- ment les oscillations causées par ces ensembles de bases affectent son premier terme. Aussi, nous proposons quatre modèles de trous d’échange analytiques auxquels nous imposons toutes les conditions actuellement connues du trou d’échange exact et la nou- velle présentée dans ce travail. Nous évaluons la performance des modèles en calculant des énergies d’échange et ses contributions à des énergies d’atomisation. On constate que les oscillations causeés par les bases de type gaussiennes peuvent compromettre la précision et la solution des modèles.