The RMT method for many-electron atomic systems in intense short-pulse laser light

We describe a new ab initio method for solving the time-dependent Schrödinger equation for multi-electron atomic systems exposed to intense short-pulse laser light. We call the method the R-matrix with time-dependence (RMT) method. Our starting point is a finite-difference numerical integrator (HELIUM), which has proved successful at describing few-electron atoms and atomic ions in strong laser fields with high accuracy. By exploiting the R-matrix division-of-space concept, we bring together a numerical method most appropriate to the multi-electron finite inner region (R-matrix basis set) and a different numerical method most appropriate to the one-electron outer region (finite difference). In order to exploit massively parallel supercomputers efficiently, we time-propagate the wavefunction in both regions by employing Arnoldi methods, originally developed for HELIUM.

Partial and total multiphoton detachment rates for H- in a perturbative B-spline-based approach

An ab initio approach has been applied to study multiphoton detachment rates for the negative hydrogen ion in the lowest nonvanishing order of perturbation theory. The approach is based on the use of B splines allowing an accurate treatment of the electronic repulsion. Total detachment rates have been determined for two- to six-photon processes as well as partial rates for detachment into the different final symmetries. It is shown that B-spline expansions can yield accurate continuum and bound-state wave functions in a very simple manner. The calculated total rates for two- and three-photon detachment are in good agreement with other perturbative calculations. For more than three-photon detachment little information has been available before now. While the total cross sections show little structure, a fair amount of structure is predicted in the partial cross sections. In the two-photon process, it is shown that the detached electrons mainly have s character. For four- and six-photon processes, the contribution from the d channel is the most important. For three- and five-photon processes p electrons dominate the electron emission spectrum. Detachment rates for s and p electrons show minima as a function of photon energy. © 1994 The American Physical Society.

Positron scattering and annihilation on noble-gas atoms

Positron scattering and annihilation on noble-gas atoms is studied ab initio using many-body theory methods for positron energies below the positronium formation threshold. We show that in this energy range, the many-body theory yields accurate numerical results and provides a near-complete understanding of the positron–noble-gas atom system. It accounts for positron-atom and electron-positron correlations, including the polarization of the atom by the positron and the nonperturbative effect of virtual positronium formation. These correlations have a large influence on the scattering dynamics and result in a strong enhancement of the annihilation rates compared to the independent-particle mean-field description. Computed elastic scattering cross sections are found to be in good agreement with recent experimental results and Kohn variational and convergent close-coupling calculations. The calculated values of the annihilation rate parameter Zeff (effective number of electrons participating in annihilation) rise steeply along the sequence of noble-gas atoms due to the increasing strength of the correlation effects, and agree well with experimental data.

Electron rescattering in strong-field photodetachment of F-

We present ab initio studies of photoelectron spectra for above threshold detachment (ATD) of F- anions in short, 1300 nm and 1800 nm laser pulses. We identify and assess the importance of electron rescattering in strong-field photodetachment of a negative ion through comparison with an analytic, Keldysh-type approach, demonstrating the capability of ab-initio computation in the challenging near-IR regime. We further assess the influence of the strong electron correlation on the photodetachment.

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 CO and atomic oxygen chemisorption on Pt(111)

Ab initio total energy calculations within a density functional theory framework have been performed for CO and atomic oxygen chemisorbed on the Pt(111) surface. Optimised geometries and chemisorption energies for CO and O on four high-symmetry sites, namely the top, bridge, fee hollow and hcp hollow sites, are presented, the coverage in all cases being 0.25 ML. The differences in CO adsorption energies between these sites are found to be small, suggesting that the potential energy surface for CO diffusion across Pt(111) is relatively flat. The 5 sigma and 2 pi molecular orbitals of CO are found to contribute to bonding with the metal. Some mixing of the 4 sigma and 1 pi molecular orbitals with metal states is also observed. For atomic oxygen, the most stable adsorption site is found to be the fee hollow site, followed in decreasing order of stability by the hcp hollow and bridge sites, with the top site being the least stable. The differences in chemisorption energies between sites for oxygen are larger than in the case of CO, suggesting a higher barrier to diffusion for atomic oxygen. The co-adsorption of CO and O has also been investigated. Calculated chemisorption energies for CO on an O/fcc-precovered surface show that of the available chemisorption sites, the top site at the oxygen atom's next-nearest neighbour surface metal atom is the most stable, with the other four sites calculated bring at least 0.29 eV less stable. The trend of CO site stability in the coadsorption system is explained in terms of a 'bonding competition' model. (C) 2000 Elsevier Science B.V. All rights reserved.

A density functional theory study of CO oxidation on Ru(0001) at low coverage

We have performed ab initio density functional theory calculations with the generalized gradient approximation to investigate CO oxidation on Ru(0001). Several reaction pathways and transition states are identified. A much higher reaction barrier compared to that on Pt(111) is determined, confirming that the Ru is very inactive for CO oxidation under UHV conditions. The origin of the reaction barrier was analyzed. It is found that in the transition state the chemisorbed O atom sits in an unfavorable bonding site and a significant competition for bonding with the same substrate atoms occurs between the CO and the chemisorbed O, resulting in the high barrier. Ab initio molecular dynamics calculations show that the activation of the chemisorbed O atom from the initial hcp hollow site (the most stable site) to the bridge site is the crucial step for the reaction. The CO oxidation on Ru(0001) via the Eley-Rideal mechanism has also been investigated. A comparison with previous theoretical work has been made. (C) 2000 American Institute of Physics. [S0021-9606(00)31223-5].

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].

Density functional theory study of the interaction between CO and on a Pt surface:CO/Pt(111), O/Pt(111), and CO/O/Pt(111)

Ab initio total energy calculations within the Density Functional Theory framework were carried out for Pt(111), Pt(111)-p(2x2)-CO, Pt(111)-p(2x2)-O, and Pt(111)-p(2x2)-(CO+O) to provide an insight into the interaction between CO and O on metal surfaces, an important issue in CO oxidation, and also in promotion and poisoning effects of catalysis. The geometrical structures of these systems were optimized with respect to the total energy, the results of which agree with existing experimental values very well. It is found that (i) the local structures of Pt(111)-p(2x2)-(CO+O), such as the bond lengths of C-O, C-Pt, and O-Pt (chemisorbed O atom with Pt), are almost the same as that in Pt(111)-p(2x2)-CO and Pt(111)-p(2x2)-O, respectively, (ii) the total valence charge density distributions in Pt(111)-p(2x2)-(CO+O) are very similar to that in Pt(111)-p(2x2)-CO, except in the region of the chemisorbed oxygen atom, and also nearly identical to that in Pt(111)-p(2x2)-O, apart from in the region of the chemisorbed CO, and (iii) the chemisorption energy of CO on a precovered Pt(111)-p(2x2)-O and the chemisorption energy of O on a precovered Pt(111)-p(2x2)CO are almost equal to that in Pt(111)-p(2x2)-CO and Pt(111)-p(2x2)-O, respectively. These results indicate that the interaction between CO and chemisorbed oxygen on a metal surface is mainly shore range in nature. The discussions of Pt-CO and Pt-O bonding and the interaction between CO and the chemisorbed oxygen atom on Pt(111) are augmented by local densities of states and real space distributions of quantum states.

Insight into CO Activation over Cu(100) under Electrochemical Conditions

The reduction of CO₂ on copper electrodes has attracted great attentions in the last decades, since it provides a sustainable approach for energy restore. During the CO₂ reduction process, the electron transfer to CO_ads is experimentally suggested to be the crucial step. In this work, we examine two possible pathways in CO activation, i.e. to generate COH_ads and CHO_ads, respectively, by performing the state-of-the-art constrained ab initio molecular dynamics simulations on the charged Cu(100) electrode under aqueous conditions, which is close to the realistic electrochemical condition. The free energy profile in the formation of COH_ads via the coupled proton and electron transfer is plotted. Furthermore, by Bader charge analyses, a linear relationship between C-O bond distance and the negative charge in CO fragment is unveiled. The formation of CHO_ads is identified to be a surface catalytic reaction, which requires the adsorption of H atom on the surface first. By comparing these two pathways, we demonstrate that kinetically the formation of COH_ads is more favored than that of CHO_ads, while CHO_ads is thermodynamically more stable. This work reveals that CO activation via COH_ads intermediate is an important pathway in electrocatalysis, which could provide some insights into CO₂ electroreduction over Cu electrodes.

Vibrational excitation of the N2+ first negative (0,0), (1,0) and (2,0) bands by electron impact: a theoretical study using the R-matrix approach

Ab initio cross section calculations for vibronic excitation using the R -matrix approach have been performed on the N 2 + molecular ion complex. A three-state close-coupling expansion is used where the electronic target states; X 2 g + , A 2 u and B 2 u + of the molecular cation are represented by a valence configuration-interaction approximation. A non-adiabatic approximation is invoked to study vibronic excitation for the first three negative bands, (0,0), (1,0) and (2,0) of the X-B transition (B 2 u + v ´ X 2 g + v ´´ ) of N 2 + . Fixed-nuclei and non-adiabatic cross section results are compared with the available experimental data for the (0,0) band and the breakdown of the adiabatic fixed-nuclei approximation is clearly evident for the vibronic excitation of the (1,0) and (2,0) bands in this molecular ion complex.

Correlation Effects in the Compton Profile of Silicon

Ab initio nonlocal pseudopotential variational quantum Monte Carlo techniques are used to compute the correlation effects on the valence momentum density and Compton profile of silicon. Our results for this case are in excellent agreement with the Lam-Platzman correction computed within the local density approximation. Within the approximations used, we rule out valence electron correlations as the dominant source of discrepancies between calculated and measured Compton profiles of silicon.

Compton profiles of Si: Pseudopotential calculation and reconstruction effects

Using an ab initio pseudopotential calculation, we compute Compton profiles of silicon along the (100), (110), and (111) directions, and then reconstruct the pseudo-wave-functions to regain the oscillatory behavior of the all-electron valence wave functions near the atomic cores. We study the effect that this reconstruction has on the Compton profiles and their anisotropies. We find a decrease in the magnitude of the profiles at small wave vectors and in their anisotropies. These changes bring the theoretical predictions closer to experimental results.

Molecular electrostatic properties of ions in an ionic liquid

We have analysed the electronic wave functions from an ab initio simulation of the ionic liquid (room temperature molten salt) dimethyl imidazolium chloride ([dmim][Cl] or [C1mim][Cl]) using localized Wannier orbitals. This allows us to assign electron density to individual ions. The probability distributions of the ionic dipole moments for an isolated ion and for ions in solution are compared. The liquid environment is found to polarize the cation by about 0.7 D and to increase the amplitude of the fluctuations in the dipole moments of both cation and anion. The relative changes in nuclear and electronic contributions are shown. The implications for classical force fields are discussed.

Probing warm dense lithium by inelastic X-ray scattering

One of the grand challenges of contemporary physics is understanding strongly interacting quantum systems comprising such diverse examples as ultracold atoms in traps, electrons in high-temperature superconductors and nuclear matter. Warm dense matter, defined by temperatures of a few electron volts and densities comparable with solids, is a complex state of such interacting matter. Moreover, the study of warm dense matter states has practical applications for controlled thermonuclear fusion, where it is encountered during the implosion phase, and it also represents laboratory analogues of astrophysical environments found in the core of planets and the crusts of old stars, Here we demonstrate how warm dense matter states can be diagnosed and structural properties can be obtained by inelastic X-ray scattering measurements on a compressed lithium sample. Combining experiments and ab initio simulations enables us to determine its microscopic state and to evaluate more approximate theoretical models for the ionic structure.