Magnetically quantized continuum distorted-wave theory in atomic and molecular collisions

The continuum distorted-wave eikonal initial-state (CDW-EIS) theory of Crothers and McCann (J Phys B 1983, 16, 3229) used to describe ionization in ion-atom collisions is generalized (G) to GCDW-EIS to incorporate the azimuthal angle dependence of each CDW in the final-state wave function. This is accomplished by the analytic continuation of hydrogenic-like wave functions from below to above threshold, using parabolic coordinates and quantum numbers including magnetic quantum numbers, thus providing a more complete set of states. At impact energies lower than 25 keVu(-1), the total ionization cross-section falls off, with decreasing energy, too quickly in comparison with experimental data. The idea behind and motivation for the GCDW-EIS model is to improve the theory with respect to experiment by including contributions from nonzero magnetic quantum numbers. We also therefore incidentally provide a new derivation of the theory of continuum distorted waves for zero magnetic quantum numbers while simultaneously generalizing it. (C) 2004 Wiley Periodicals, Inc.

Quantum annealing of an ising spin-glass by Green's function Monte Carlo

We present an implementation of quantum annealing (QA) via lattice Green's function Monte Carlo (GFMC), focusing on its application to the Ising spin glass in transverse field. In particular, we study whether or not such a method is more effective than the path-integral Monte Carlo- (PIMC) based QA, as well as classical simulated annealing (CA), previously tested on the same optimization problem. We identify the issue of importance sampling, i.e., the necessity of possessing reasonably good (variational) trial wave functions, as the key point of the algorithm. We performed GFMC-QA runs using such a Boltzmann-type trial wave function, finding results for the residual energies that are qualitatively similar to those of CA (but at a much larger computational cost), and definitely worse than PIMC-QA. We conclude that, at present, without a serious effort in constructing reliable importance sampling variational wave functions for a quantum glass, GFMC-QA is not a true competitor of PIMC-QA.

Single ionization of helium by antiproton impact

We study the ionization of helium (fie) in collision with antiprotons (p) in the energy range from 10 keV to 1000 keV. We adopt a semiclassical single center close coupling approach in which the wave function for the electron is expanded in a B-spline basis centered on the nucleus of the He atom, The calculations are performed using two different models: the independent particle (IP) model and the one-electron (OE) approximation. The interaction between the active electron and the rest of the atom, i.e. passive electron and nucleus, is represented by a model potential. The results obtained are compared with experimental data as well as with existing theoretical calculations. (c) 2005 Published by Elsevier B.V.

Positronium-atom collisions

New results are presented for Ps(1s) scattering by H(1s), He(1(1)S) and Li(2s). Calculations have been performed in a coupled state framework, usually employing pseudostates, and allowing for excitation of both the Ps and the atom. In the Ps(1s)-H(1s) calculations the H- formation channel has also been included using a highly accurate H- wave function. Resonances resulting from unstable states in which the positron orbits H- have been calculated and analysed. The new Ps(1s)-He(1(1)S) calculations still fail to resolve existing discrepancies between theory and experiment at very low energies. The possible importance of the Ps(-) formation channel in all three collision systems is discussed. (C) 2004 Elsevier B.V. All rights reserved.

Time-dependent R -matrix theory for ultrafast atomic processes

We describe an ab initio nonperturbative time-dependent R-matrix theory for ultrafast atomic processes. This theory enables investigations of the interaction of few-femtosecond and -attosecond pulse lasers with complex multielectron atoms and atomic ions. A derivation and analysis of the basic equations are given, which propagate the atomic wave function in the presence of the laser field forward in time in the internal and external R-matrix regions. To verify the accuracy of the approach, we investigate two-photon ionization of Ne irradiated by an intense laser pulse and compare current results with those obtained using the R-matrix Floquet method and an alternative time-dependent method. We also verify the capability of the current approach by applying it to the study of two-dimensional momentum distributions of electrons ejected from Ne due to irradiation by a sequence of 2 as light pulses in the presence of a 780 nm laser field.

Electron-impact excitation of Fe II: Effective collision strengths for optically allowed fine-structure transitions

In this paper. we present collision strengths and Maxwellian averaged effective collision strengths for the electron-impact excitation of Fe II. We consider specifically the optically allowed lines for transitions from the 3d(6)4s and 3d(7) even parity configuration states to the 3d(6)4p odd parity configuration levels. The parallel suite of Breit-Pauli codes are utilized to compute the collision cross-sections where relativistic effects are included explicitly in both the target and the scattering approximation. A total of 100 LS or 262-jj levels formed from the basis configurations 3d(6)4s, 3d(7) and 3d(6)4p were included in the wave-function representation of the target, including all doublet. quartet and sextet terms. The Maxwellian averaged effective collision strengths are computed across a wide range of electron temperatures from 100 to 100,000 K, temperatures of importance in astrophysical and plasma applications. A detailed comparison is made with previous works and significant differences were found to occur for some of the transitions considered. We conclude that in order to obtain converged collision strengths and effective collision strengths for these allowed transitions it is necessary to include contributions from partial waves up to L = 50 explicitly in the calculation, and in addition, account for contributions from even higher partial waves through a "top up" procedure.

Effective collision strengths for electron-impact excitation of S X

Effective collision strengths computed by the R-matrix method are presented for the electron-impact excitation of nitrogen-like S X. The total wave function used in the expansion includes the lowest 11 eigenstates of S X which arise from the 2s(2)2p(3), 2s2p(4), 2p(5) and 2s(2)2p(2)3s configurations. These 11 LS target states correspond to 22 fine-structure levels, giving 231 possible transitions. All the effective collision strengths for these transitions are tabulated in the range log T(K) = 4.6 to log T(K) = 6.7. The energy level values and oscillator strengths for allowed transitions are also tabulated. The effective collision strengths were calculated by averaging the electron collision strengths over a Maxwellian distribution of velocities. The present effective collision strengths are the only results currently available for these fine-structure transition rates. (C) 2000 Academic Press.

Density-potential mapping in time-dependent density-functional theory

The key questions of uniqueness and existence in time-dependent density-functional theory are usually formulated only for potentials and densities that are analytic in time. Simple examples, standard in quantum mechanics, lead, however, to nonanalyticities. We reformulate these questions in terms of a nonlinear Schroedinger equation with a potential that depends nonlocally on the wave function.

Combined R-matrix eigenstate basis set and finite-difference propagation method for the time-dependent Schrodinger equation: The one-electron case

In this work we present the theoretical framework for the solution of the time-dependent Schrödinger equation (TDSE) of atomic and molecular systems under strong electromagnetic fields with the configuration space of the electron’s coordinates separated over two regions; that is, regions I and II. In region I the solution of the TDSE is obtained by an R-matrix basis set representation of the time-dependent wave function. In region II a grid representation of the wave function is considered and propagation in space and time is obtained through the finite-difference method. With this, a combination of basis set and grid methods is put forward for tackling multiregion time-dependent problems. In both regions, a high-order explicit scheme is employed for the time propagation. While, in a purely hydrogenic system no approximation is involved due to this separation, in multielectron systems the validity and the usefulness of the present method relies on the basic assumption of R-matrix theory, namely, that beyond a certain distance (encompassing region I) a single ejected electron is distinguishable from the other electrons of the multielectron system and evolves there (region II) effectively as a one-electron system. The method is developed in detail for single active electron systems and applied to the exemplar case of the hydrogen atom in an intense laser field.

Low-energy excitations of a boson pair in a double-well trap

The states of a boson pair in a one-dimensional double-well potential are investigated. Properties of the ground and lowest excited states of this system are studied, including the two-particle wave function, momentum pair distribution, and entanglement. The effects of varying both the barrier height and the effective interaction strength are investigated.

Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

Positron scattering and annihilation in hydrogenlike ions

Diagrammatic many-body theory is used to calculate the scattering phase shifts, normalized annihilation rates Zeff, and annihilation ? spectra for positron collisions with the hydrogenlike ions He+, Li2+, B4+, and F8+. Short-range electron-positron correlations and longer-range positron-ion correlations are accounted for by evaluating nonlocal corrections to the annihilation vertex and the exact positron self-energy. The numerical calculation of the many-body theory diagrams is performed using B-spline basis sets. To elucidate the role of the positron-ion repulsion, the annihilation rate is also estimated analytically in the Coulomb-Born approximation. It is found that the energy dependence and magnitude of Zeff are governed by the Gamow factor that characterizes the suppression of the positron wave function near the ion. For all of the H-like ions, the correlation enhancement of the annihilation rate is found to be predominantly due to corrections to the annihilation vertex, while the corrections to the positron wave function play only a minor role. Results of the calculations for s-, p-, and d-wave incident positrons of energies up to the positronium-formation threshold are presented. Where comparison is possible, our values are in excellent agreement with the results obtained using other, e.g., variational, methods. The annihilation-vertex enhancement factors obtained in the present calculations are found to scale approximately as 1+(1.6+0.46l)/Zi, where Zi is the net charge of the ion and l is the positron orbital angular momentum. Our results for positron annihilation in H-like ions provide insights into the problem of positron annihilation with core electrons in atoms and condensed matter systems, which have similar binding energies.

Liouville-space R-matrix-Floquet description of atomic radiative processes involving autoionizing states in the presence of intense electromagnetic fields

A reduced-density-operator description is developed for coherent optical phenomena in many-electron atomic systems, utilizing a Liouville-space, multiple-mode Floquet–Fourier representation. The Liouville-space formulation provides a natural generalization of the ordinary Hilbert-space (Hamiltonian) R-matrix-Floquet method, which has been developed for multi-photon transitions and laser-assisted electron–atom collision processes. In these applications, the R-matrix-Floquet method has been demonstrated to be capable of providing an accurate representation of the complex, multi-level structure of many-electron atomic systems in bound, continuum, and autoionizing states. The ordinary Hilbert-space (Hamiltonian) formulation of the R-matrix-Floquet method has been implemented in highly developed computer programs, which can provide a non-perturbative treatment of the interaction of a classical, multiple-mode electromagnetic field with a quantum system. This quantum system may correspond to a many-electron, bound atomic system and a single continuum electron. However, including pseudo-states in the expansion of the many-electron atomic wave function can provide a representation of multiple continuum electrons. The 'dressed' many-electron atomic states thereby obtained can be used in a realistic non-perturbative evaluation of the transition probabilities for an extensive class of atomic collision and radiation processes in the presence of intense electromagnetic fields. In order to incorporate environmental relaxation and decoherence phenomena, we propose to utilize the ordinary Hilbert-space (Hamiltonian) R-matrix-Floquet method as a starting-point for a Liouville-space (reduced-density-operator) formulation. To illustrate how the Liouville-space R-matrix-Floquet formulation can be implemented for coherent atomic radiative processes, we discuss applications to electromagnetically induced transparency, as well as to related pump–probe optical phenomena, and also to the unified description of radiative and dielectronic recombination in electron–ion beam interactions and high-temperature plasmas.

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.