Fully differential cross sections for 3.6 MeV u(-1) AuZp++He collisions

Recent experimental data for fully differential cross sections have been compared to various continuum-distorted-wave eikonal-initial-state models without much success, despite good agreement with double-differential cross sections. A four-body model is formulated here and results are presented both when the internuclear potential is omitted and when it is included. They are compared with recent experimental data for fully differential cross sections for 3.6 MeV/u Au-P(Z)++He collisions, Z(P)=24,53.

Ionization from the outer shell of Ar by proton impact

Recent results for proton-argon total ionization cross sections [Kirchner Phys. Rev. Lett. 79, 1658 (1997)] show large disagreement between theory and experiment for energies below 80 keV. To address this problem we have employed a recently developed theoretical method with a more pragmatic approach to the charge screening both in the initial and final channels. The target is considered as a one-electron atom and the interactions between this active electron and remaining target electrons are treated by a model potential including both short- and long-range effects. In the final channel the usual product of two continuum distorted wave functions each associated with a distinct electron-nucleus interaction is used. New results in the present calculation show good agreement in total cross sections for the energy range 10-300 keV with the measurement of Rudd [Rev. Mod. Phys. 57, 965 (1985)].

On the semiclassical approach to cold atomic collisions

We extend the semiclassical description of two-state atomic collisions to low energies for which the impact parameter treatment fails. The problem reduces to solving a system of first-order differential equations with coefficients whose semiclassical asymptotes experience the Stokes phenomenon in the complex coordinate plane. Primitive semiclassical and uniform Airy approximations are discussed.

Multi-ionization of helium and lithium using the independent electron and independent event models with intrinsic CDW

Cross sections for the multi-ionization of He and Li are presented for impact energies in the range of 50 to 1000 keV/amu. These are calculated using the eikonal initial state approximation to represent the input and exit channels of the active electrons. The ionization process is simulated in a variety of ways, most notably an attempt to account for the effects of electron correlation via the inclusion of a continuum density of states (CDS) term. Inadequacies, of the CDW formulation at small impact parameters, and of the models themselves, are discussed and conclusions are drawn on what repercussions this has for the cross sections calculated.

Total and single differential cross sections for simple resonant collisions using a fully orthonormal continuum-distorted-wave basis

Simple electron capture processes are studied using an orthonormal two state continuum-distorted-wave (CDW) basis. The suitability of the basis set is tested by comparing predictions for total and differential cross sections with available experimental data. Overall good agreement is obtained and the authors conclude that a relatively small CDW basis set may be suitable to model a wide variety of low-energy collisions if the members of this extended set are astutely chosen.

Other uses of the program 'ARGON.f90'

The purpose of this communication is to show that the program 'ARGON.f90' can be simply extended to model ionization from the excited states of atoms where the active electron has a principal quantum number less than or equal to 3. This fact is illustrated by considering a relatively simple collision involving a proton and a neutral hydrogen atom with principal quantum number n = 2. (C) 2005 Published by Elsevier B.V.

Green function for the diffusion limit of one-dimensional continuous time random walks

It is shown how the fractional probability density diffusion equation for the diffusion limit of one-dimensional continuous time random walks may be derived from a generalized Markovian Chapman-Kolmogorov equation. The non-Markovian behaviour is incorporated into the Markovian Chapman-Kolmogorov equation by postulating a Levy like distribution of waiting times as a kernel. The Chapman-Kolmogorov equation so generalised then takes on the form of a convolution integral. The dependence on the initial conditions typical of a non-Markovian process is treated by adding a time dependent term involving the survival probability to the convolution integral. In the diffusion limit these two assumptions about the past history of the process are sufficient to reproduce anomalous diffusion and relaxation behaviour of the Cole-Cole type. The Green function in the diffusion limit is calculated using the fact that the characteristic function is the Mittag-Leffler function. Fourier inversion of the characteristic function yields the Green function in terms of a Wright function. The moments of the distribution function are evaluated from the Mittag-Leffler function using the properties of characteristic functions and a relation between the powers of the second moment and higher order even moments is derived. (C) 2004 Elsevier B.V. All rights reserved.

Escape times for rigid Brownian rotators in a bistable potential from the time evolution of the Green function and the characteristic time of the probability evolution

The greatest relaxation time for an assembly of three- dimensional rigid rotators in an axially symmetric bistable potential is obtained exactly in terms of continued fractions as a sum of the zero frequency decay functions (averages of the Legendre polynomials) of the system. This is accomplished by studying the entire time evolution of the Green function (transition probability) by expanding the time dependent distribution as a Fourier series and proceeding to the zero frequency limit of the Laplace transform of that distribution. The procedure is entirely analogous to the calculation of the characteristic time of the probability evolution (the integral of the configuration space probability density function with respect to the position co-ordinate) for a particle undergoing translational diffusion in a potential; a concept originally used by Malakhov and Pankratov (Physica A 229 (1996) 109). This procedure allowed them to obtain exact solutions of the Kramers one-dimensional translational escape rate problem for piecewise parabolic potentials. The solution was accomplished by posing the problem in terms of the appropriate Sturm-Liouville equation which could be solved in terms of the parabolic cylinder functions. The method (as applied to rotational problems and posed in terms of recurrence relations for the decay functions, i.e., the Brinkman approach c.f. Blomberg, Physica A 86 (1977) 49, as opposed to the Sturm-Liouville one) demonstrates clearly that the greatest relaxation time unlike the integral relaxation time which is governed by a single decay function (albeit coupled to all the others in non-linear fashion via the underlying recurrence relation) is governed by a sum of decay functions. The method is easily generalized to multidimensional state spaces by matrix continued fraction methods allowing one to treat non-axially symmetric potentials, where the distribution function is governed by two state variables. (C) 2001 Elsevier Science B.V. All rights reserved.

Nonlinear response of permanent dipoles in a uniaxial potential to alternating fields

It is shown how the existing theory of the dynamic Kerr effect and nonlinear dielectric relaxation based on the noninertial Brownian rotation of noninteracting rigid dipolar particles may be generalized to take into account interparticle interactions using the Maier-Saupe mean field potential. The results (available in simple closed form) suggest that the frequency dependent nonlinear response provides a method of measuring the Kramers escape rate (or in the analogous problem of magnetic relaxation of fine single domain ferromagnetic particles, the superparamagnetic relaxation time).

A variational approach to the relaxation of single domain magnetic particles based on Brown's model

Brown's model for the relaxation of the magnetization of a single domain ferromagnetic particle is considered. This model results in the Fokker-Planck equation of the process. The solution of this equation in the cases of most interest is non- trivial. The probability density of orientations of the magnetization in the Fokker-Planck equation can be expanded in terms of an infinite set of eigenfunctions and their corresponding eigenvalues where these obey a Sturm-Liouville type equation. A variational principle is applied to the solution of this equation in the case of an axially symmetric potential. The first (non-zero) eigenvalue, corresponding to the largest time constant, is considered. From this we obtain two new results. Firstly, an approximate minimising trial function is obtained which allows calculation of a rigorous upper bound. Secondly, a new upper bound formula is derived based on the Euler-Lagrange condition. This leads to very accurate calculation of the eigenvalue but also, interestingly, from this, use of the simplest trial function yields an equivalent result to the correlation time of Coffey et at. and the integral relaxation time of Garanin. (C) 2004 Elsevier B.V. All rights reserved.

The role of Mittag-Leffler functions in anomalous relaxation

The inertia-corrected Debye model of rotational Brownian motion of polar molecules was generalized by Coffey et al. [Phys. Rev. E, 65, 32 102 (2002)] to describe fractional dynamics and anomalous rotational diffusion. The linear- response theory of the normalized complex susceptibility was given in terms of a Laplace transform and as a function of frequency. The angular-velocity correlation function was parametrized via fractal Mittag-Leffler functions. Here we apply the latter method and complex-contour integral- representation methods to determine the original time-dependent amplitude as an inverse Laplace transform using both analytical and numerical approaches, as appropriate. (C) 2004 Elsevier B.V. All rights reserved.

An exactly solvable four transition point model

Heavy particle collisions, in particular low-energy ion-atom collisions, are amenable to semiclassical JWKB phase integral analysis in the complex plane of the internuclear separation. Analytic continuation in this plane requires due attention to the Stokes phenomenon which parametrizes the physical mechanisms of curve crossing, non-crossing, the hybrid Nikitin model, rotational coupling and predissociation. Complex transition points represent adiabatic degeneracies. In the case of two or more such points, the Stokes constants may only be completely determined by resort to the so-called comparison- equation method involving, in particular, parabolic cylinder functions or Whittaker functions and their strong-coupling asymptotics. In particular, the Nikitin model is a two transition-point one-double-pole problem in each half-plane corresponding to either ingoing or outgoing waves. When the four transition points are closely clustered, new techniques are required to determine Stokes constants. However, such investigations remain incomplete, A model problem is therefore solved exactly for scattering along a one-dimensional z-axis. The energy eigenvalue is b(2)-a(2) and the potential comprises -z(2)/2 (parabolic) and -a(2) + b(2)/2z(2) (centrifugal/centripetal) components. The square of the wavenumber has in the complex z-plane, four zeros each a transition point at z = +/-a +/- ib and has a double pole at z = 0. In cases (a) and (b), a and b are real and unitarity obtains. In case (a) the reflection and transition coefficients are parametrized by exponentials when a(2) + b(2) > 1/2. In case (b) they are parametrized by trigonometrics when a(2) + b(2) <1/2 and total reflection is achievable. In case (c) a and b are complex and in general unitarity is not achieved due to loss of flux to a continuum (O'Rourke and Crothers, 1992 Proc. R. Sec. 438 1). Nevertheless, case (c) coefficients reduce to (a) or (b) under appropriate limiting conditions. Setting z = ht, with h a real constant, an attempt is made to model a two-state collision problem modelled by a pair of coupled first-order impact parameter equations and an appropriate (T) over tilde-tau relation, where (T) over tilde is the Stueckelberg variable and tau is the reduced or scaled time. The attempt fails because (T) over tilde is an odd function of tau, which is unphysical in a real collision problem. However, it is pointed out that by applying the Kummer exponential model to each half-plane (O'Rourke and Crothers 1994 J. Phys. B: At. Mel. Opt. Phys. 27 2497) the current model is in effect extended to a collision problem with four transition points and a double pole in each half-plane. Moreover, the attempt in itself is not a complete failure since it is shown that the result is a perfect diabatic inelastic collision for a traceless Hamiltonian matrix, or at least when both diagonal elements are odd and the off-diagonal elements equal and even.

Electron-impact ionization of atomic hydrogen close to threshold

A systematic study of the ionization of atomic hydrogen by electron impact from 0.3 eV to a few eV above the ionization threshold has been carried out using a semiclassical-quantal calculation. Differential and integrated cross sections are presented at 0.3 eV above the energy threshold. Triple- differential cross sections (TDCS) are presented at constant theta(12) geometry where theta(12)=180degrees and 150degrees. Good agreement is achieved with the measurement [Roder, Phys. Rev. Lett. 79, 1666 (1997)] and calculations based on exterior complex scaling at 2 eV and 4 eV above threshold. Results of triple-differential cross sections are also presented at 0.3, 0.5, and 1.0 eV above threshold at both theta(12)=180degrees and 150degrees. At theta(12)=180degrees the small local maximum in the TDCS around theta(1)=90degrees reported by Pan and Starace [Phys. Rev. A 45, 4588 (1992)] at 0.5 eV above threshold is not observed in our calculation at energies down to 0.3 eV above threshold. The shape of our double differential cross sections seems to disagree qualitatively with the available calculations as we found two local maxima around 15degrees and 165degrees in our calculation. Single differential cross sections in our formulation appear naturally as a function of total excess energy E and, therefore, constant for all combinations of individual electron energies E-1 and E- 2 with E=E-1+E-2. Total ionization cross sections are also compared with measurement and available theoretical calculations and found to be in reasonably good agreement up to 10 eV above ionization threshold.