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.

A discrete time-dependent method for metastable atoms and molecules in intense fields

The full-dimensional time-dependent Schrodinger equation for the electronic dynamics of single-electron systems in intense external fields is solved directly using a discrete method. Our approach combines the finite-difference and Lagrange mesh methods. The method is applied to calculate the quasienergies and ionization probabilities of atomic and molecular systems in intense static and dynamic electric fields. The gauge invariance and accuracy of the method is established. Applications to multiphoton ionization of positronium, the hydrogen atom and the hydrogen molecular ion are presented. At very high laser intensity, above the saturation threshold, we extend the method using a scaling technique to estimate the quasienergies of metastable states of the hydrogen molecular ion. The results are in good agreement with recent experiments. (C) 2004 American Institute of Physics.

A multiconfigurational time-dependent Hartree-Fock method for excited electronic states. I. General formalism and application to open-shell states

The solution of the time-dependent Schrodinger equation for systems of interacting electrons is generally a prohibitive task, for which approximate methods are necessary. Popular approaches, such as the time-dependent Hartree-Fock (TDHF) approximation and time-dependent density functional theory (TDDFT), are essentially single-configurational schemes. TDHF is by construction incapable of fully accounting for the excited character of the electronic states involved in many physical processes of interest; TDDFT, although exact in principle, is limited by the currently available exchange-correlation functionals. On the other hand, multiconfigurational methods, such as the multiconfigurational time-dependent Hartree-Fock (MCTDHF) approach, provide an accurate description of the excited states and can be systematically improved. However, the computational cost becomes prohibitive as the number of degrees of freedom increases, and thus, at present, the MCTDHF method is only practical for few-electron systems. In this work, we propose an alternative approach which effectively establishes a compromise between efficiency and accuracy, by retaining the smallest possible number of configurations that catches the essential features of the electronic wavefunction. Based on a time-dependent variational principle, we derive the MCTDHF working equation for a multiconfigurational expansion with fixed coefficients and specialise to the case of general open-shell states, which are relevant for many physical processes of interest. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3600397]

SOLVING THE INHOMOGENEOUS SCHRODINGER-EQUATION ON A BIDIRECTIONAL PIPELINE OF TRANSPUTERS

In this paper we describe the design of a parallel solution of the inhomogeneous Schrodinger equation, which arises in the construction of continuum orbitals in the R-matrix theory of atomic continuum processes. A prototype system is described which has been programmed in occam2 and implemented on a bi-directional pipeline of transputers. Some timing results for the prototype system are presented, and the development of a full production system is discussed.

Modeling relativistic soliton interactions in overdense plasmas: A perturbed nonlinear Schrodinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrodinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrodinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude

Dissociative ionization of molecules in intense laser fields

Accurate and efficient grid based techniques for the solution of the time-dependent Schrodinger equation for few-electron diatomic molecules irradiated by intense, ultrashort laser pulses are described. These are based on hybrid finite-difference, Lagrange mesh techniques. The methods are applied in three scenarios, namely H-2(+) with fixed internuclear separation, H-2(+) with vibrating nuclei and H-2 with fixed internuclear separation and illustrative results presented.

Dynamic tunnelling ionization of H-+(2) in intense fields

Intense-field ionization of the hydrogen molecular ion by linearly polarized light is modelled by direct solution of the fixed-nuclei time-dependent Schrodinger equation and compared with recent experiments. Parallel transitions are calculated using algorithms which exploit massively parallel computers. We identify and calculate dynamic tunnelling ionization resonances that depend on laser wavelength and intensity, and molecular bond length. Results for lambda similar to 1064 nm are consistent with static tunnelling ionization. At shorter wavelengths lambda similar to 790 nm large dynamic corrections are observed. The results agree very well with recent experimental measurements of the ion spectra. Our results reproduce the single peak resonance and provide accurate ionization rate estimates at high intensities. At lower intensities our results confirm a double peak in the ionization rate as the bond length varies.

Ionization dynamics of laser-driven H-2(+)

We set out aspects of a numerical algorithm used in solving the full-dimensionality time-dependent Schrodinger equation describing the electronic motion of the hydrogen molecular ion driven by an intense, linearly polarized laser pulse aligned along the molecular axis. This algorithm has been implemented within the fixed inter-nuclear separation approximation in a parallel computer code, a brief summary of which is given. Ionization rates are calculated and compared with results from other methods, notably the time-independent Floquet method. Our results compare very favourably with the precise predictions of the Floquet method, although there is some disagreement with other wavepacket calculations. Visualizations of the electron dynamics are also presented in which electron rescattering is observed.

Double-Electron Above-Threshold Ionization Resonances as Interference Phenomena

We report calculations of double ionization energy spectra and momentum distributions of laser-driven helium due to few-cycle pulses of wavelength 195 nm. The results are obtained from full-dimensional numerical integration of the two electron time-dependent Schr¨odinger equation. A momentum-space analysis of doubly ionizing wavepackets shows that the concentric-ring structure of above-threshold double ionization, together with the associated structure of peaks in the total kinetic energy spectrum, may be attributed to wavepacket interference effects, where at least two doubly-ionizing wavepackets from different recollision events populate the same spatial hemisphere.

Robust nonadiabatic molecular dynamics for metals and insulators

We present a new formulation of the correlated electron-ion dynamics (CEID) scheme, which systematically improves Ehrenfest dynamics by including quantum fluctuations around the mean-field atomic trajectories. We show that the method can simulate models of nonadiabatic electronic transitions and test it against exact integration of the time-dependent Schrodinger equation. Unlike previous formulations of CEID, the accuracy of this scheme depends on a single tunable parameter which sets the level of atomic fluctuations included. The convergence to the exact dynamics by increasing the tunable parameter is demonstrated for a model two level system. This algorithm provides a smooth description of the nonadiabatic electronic transitions which satisfies the kinematic constraints (energy and momentum conservation) and preserves quantum coherence. The applicability of this algorithm to more complex atomic systems is discussed.

High-energy cutoff in the spectrum of strong-field nonsequential double ionization

Electron energy distributions of singly and doubly ionized helium in an intense 390 nm laser field have been measured at two intensities (0.8 PW/cm(2) and 1.1 PW/cm(2), where PW equivalent to 10(15) W/cm(2)). Numerical solutions of the full-dimensional time-dependent helium Schrodinger equation show excellent agreement with the experimental measurements. The high-energy portion of the two-electron energy distributions reveals an unexpected 5U(p) cutoff for the double ionization (DI) process and leads to a proposed model for DI below the quasiclassical threshold.

Collapse versus growth for a Bose-Einstein condensate with attractive interactions

Using a model potential approach, we study the time-dependent behavior of a Bose-Einstein condensate with negative scattering length during its collapse in the zero-temperature limit. The condensate is modeled through an effective potential, which linearizes the Schrodinger equation, in order to obtain an intuitive visualization of the dynamics of the condensate. We find that a substantial fraction of the condensate survives the collapse. The origin for this survival is the reappearance of a barrier in the effective potential during the collapse. In contrast to previous calculations, the present calculations indicate that the size of the residual condensate strongly depends on the growth rate of the condensate. The present results are compared to other theoretical calculations and to experimental work.

Tools for analysing configuration interaction wavefunctions

The configuration interaction (CI) approach to quantum chemical calculations is a well-established means of calculating accurately the solution to the Schrodinger equation for many-electron systems. It represents the many-body electron wavefunction as a sum of spin-projected Slater determinants of orthogonal one-body spin-orbitals. The CI wavefunction becomes the exact solution of the Schrodinger equation as the length of the expansion becomes infinite, however, it is a difficult quantity to visualise and analyse for many-electron problems. We describe a method for efficiently calculating the spin-averaged one- and two-body reduced density matrices rho(psi)((r) over bar; (r) over bar' ) and Gamma(psi)((r) over bar (1), (r) over bar (2); (r) over bar'(1), (r) over bar'(2)) of an arbitrary CI wavefunction Psi. These low-dimensional functions are helpful tools for analysing many-body wavefunctions; we illustrate this for the case of the electron-electron cusp. From rho and Gamma one can calculate the matrix elements of any one- or two-body spin-free operator (O) over cap. For example, if (O) over cap is an applied electric field, this field can be included into the CI Hamiltonian and polarisation or gating effects may be studied for finite electron systems. (C) 2003 Elsevier B.V. All rights reserved.

Geometrical dependence in photoionization of H-2(+) in high-intensity, high-frequency, ultrashort laser pulses

The ionization dynamics of H2 + exposed to high-intensity, high-frequency, ultrashort laser pulses is investigated with two theoretical approaches. The time-dependent Schrödinger equation is solved by a direct numerical method, and a simple two-center interference-diffraction model is studied. The energy and angular distributions of the photoelectron for various internuclear distances and relative orientations between the internuclear axis of the molecule and the polarization of the field are calculated. The main features of the photoelectron spectrum pattern are described well by the interference-diffraction model, and excellent quantitative agreement between the two methods is found. The effect of quantal vibration on the photoelectron spectrum is also calculated. We find that vibrational average produces some broadening of the main features, but that the patterns remain clearly distinguishable.

Dissociation of H2+ from a short, intense, infrared laser pulse: proton emission spectra and pulse calibration

The proton energy spectrum from photodissociation of the hydrogen molecular ion by short intense pulses of infrared light is calculated. The time-dependent Schrödinger equation is discretized and integrated. For few-cycle pulses one can resolve vibrational structure, arising from the experimental preparation of the molecular ion. We calculate the corresponding energy spectrum and analyse the dependence on the pulse time delay, pulse length and intensity of the laser for ? ~ 790 nm. We conclude that the proton spectrum is a sensitive probe of both the vibrational populations and phases, and allows us to distinguish between adiabatic and nonadiabatic dissociation. Furthermore, the sensitivity of the proton spectrum from H2+ is a practical means of calibrating the pulse. Our results are compared with recent measurements of the proton spectrum for 65 fs pulses using a Ti:Sapphire laser (? ~ 790 nm) including molecular orientation and focal-volume averaging. Integrating over the laser focal volume, for the intensity I ~ 3 × 1015 W cm-2, we find our results are in excellent agreement with these experiments.