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.

Finite amplitude envelope solitons in pair-ion plasmas

The nonlinear coupling between finite amplitude ion thermal waves (ITWs) and quasistationary density perturbations in a pair-ion plasma is considered. A generalized nonlinear Schrödinger equation is derived for the ITW electric field envelope, accounting for large amplitude quasistationary plasma slow motion describing the ITW ponderomotive force. The present theory accounts for the trapping of ITWs in a large amplitude ion density hole. The small amplitude limit is considered and exact analytical solutions are obtained. Finite amplitude solutions are obtained numerically and their characteristics are discussed.

Nonlinear dynamics and solitons in Debye bi-crystals

The nonlinear dynamics of longitudinal dust lattice waves propagating in a dusty plasma bi-crystal is investigated. A “diatomic”-like one-dimensional dust lattice configuration is considered, consisting of two distinct dust grain species with different charges and masses. Two different frequency dispersion modes are obtained in the linear limit, namely, an optical and an acoustic wave dispersion branch. Nonlinear solitary wave solutions are shown to exist in both branches, by considering the continuum limit for lattice excitations in different nonlinear potential regimes. For this purpose, a generalized Boussinesq and an extended Korteweg de Vries equation is derived, for the acoustic mode excitations, and their exact soliton solutions are provided and compared. For the optic mode, a nonlinear Schrödinger-type equation is obtained, which is shown to possess bright- (dark-) type envelope soliton solutions in the long (short, respectively) wavelength range. Optic-type longitudinal wavepackets are shown to be generally unstable in the continuum limit, though this is shown not to be the rule in the general (discrete) case.

Effect of quantization of vibrations on the structural properties of crystals

We study the structural effects produced by the quantization of vibrational degrees of freedom in periodic crystals at zero temperature. To this end we introduce a methodology based on mapping a suitable subspace of the vibrational manifold and solving the Schrödinger equation in it. A number of increasingly accurate approximations ranging from the quasiharmonic approximation (QHA) to the vibrational self-consistent field (VSCF) method and the exact solution are described. A thorough analysis of the approximations is presented for model monatomic and hydrogen-bonded chains, and results are presented for a linear H-F chain where the potential-energy surface is obtained via first-principles electronic structure calculations. We focus on quantum nuclear effects on the lattice constant and show that the VSCF is an excellent approximation, meaning that correlation between modes is not extremely important. The QHA is excellent for covalently bonded mildly anharmonic systems, but it fails for hydrogen-bonded ones. In the latter, the zero-point energy exhibits a nonanalytic behavior at the lattice constant where the H atoms center, which leads to a spurious secondary minimum in the quantum-corrected energy curve. An inexpensive anharmonic approximation of noninteracting modes appears to produce rather good results for hydrogen-bonded chains for small system sizes. However, it converges to the incorrect QHA results for increasing size. Isotope effects are studied for the first-principles H-F chain. We show how the lattice constant and the H-F distance increase with decreasing mass and how the QHA proves to be insufficient to reproduce this behavior.

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.

From the UV to the static-field limit: rates and scaling laws of intense-field ionization of helium

We present high-accuracy calculations of ionization rates of helium at UV (195 nm) wavelengths. The data are obtained from full-dimensionality integrations of the helium-laser time-dependent Schrödinger equation. Comparison is made with our previously obtained data at 390 nm and 780 nm. We show that scaling laws introduced by Parker et al extend unmodified from the near-infrared limit into the UV limit. Static-field ionization rates of helium are also obtained, again from time-dependent full-dimensionality integrations of the helium Schrödinger equation. We compare the static-field ionization results with those of Scrinzi et al and Themelis et al, who also treat the full-dimensional helium atom, but with time-independent methods. Good agreement is obtained.

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.

Time delay between photoemission from the 2p and 2s subshells of neon

The R-matrix incorporating time (RMT) method is a method developed recently for solving the time-dependent Schrödinger equation for multielectron atomic systems exposed to intense short-pulse laser light. We have employed the RMT method to investigate the time delay in the photoemission of an electron liberated from a 2p orbital in a neon atom with respect to one released from a 2s orbital following absorption of an attosecond xuv pulse. Time delays due to xuv pulses in the range 76-105 eV are presented. For an xuv pulse at the experimentally relevant energy of 105.2 eV, we calculate the time delay to be 10.2±1.3 attoseconds (as), somewhat larger than estimated by other theoretical calculations, but still a factor of 2 smaller than experiment. We repeated the calculation for a photon energy of 89.8 eV with a larger basis set capable of modeling correlated-electron dynamics within the neon atom and the residual Ne ion. A time delay of 14.5±1.5 as was observed, compared to a 16.7±1.5 as result using a single-configuration representation of the residual Ne⁺ ion. 

Modulational instability in asymmetric coupled wave functions

The evolution of the amplitude of two nonlinearly interacting waves is considered, via a set of coupled nonlinear Schrödinger-type equations. The dynamical profile is determined by the wave dispersion laws (i.e. the group velocities and the group velocity dispersion terms) and the nonlinearity and coupling coefficients, on which no assumption is made. A generalized dispersion relation is obtained, relating the frequency and wave-number of a small perturbation around a coupled monochromatic (Stokes') wave solution. Explicitly stability criteria are obtained. The analysis reveals a number of possibilities. Two (individually) stable systems may be destabilized due to coupling. Unstable systems may, when coupled, present an enhanced instability growth rate, for an extended wave number range of values. Distinct unstable wavenumber windows may arise simultaneously.

Nonlinear Modulated Envelope Electrostatic Wavepacket Propagation in Plasmas

A brief review of the occurrence of amplitude modulated structures in space and laboratory plasmas is provided, followed by a theoretical analysis of the mechanism of carrier wave (self-) interaction, with respect to electrostatic plasma modes. A generic collisionless unmagnetized fluid model is employed. Both cold-(zero-temperature) and warm-(finite temperature) fluid descriptions are considered and compared. The weakly nonlinear oscillation regime is investigated by applying a multiple scale (reductive perturbation) technique and a Nonlinear Schrödinger Equation (NLSE) is obtained, describing the evolution of the slowly varying wave amplitude in time and space. The amplitude’s stability profile reveals the possibility of modulational instability to occur under the influence of external perturbations. The NLSE admits exact localized envelope (solitary wave) solutions of bright (pulses) or dark (holes, voids) type, whose characteristics depend on intrinsic plasma parameters. The role of perturbation obliqueness (with respect to the propagation direction), finite temperature and — possibly — defect (dust) concentration is explicitly considered. The relevance of this description with respect to known electron-ion (e-i) as well as dusty (complex) plasma modes is briefly discussed. © 2004 American Institute of Physics

Electromagnetic rogue waves in beam-plasma interactions

The occurrence of rogue waves (freak waves) associated with electromagnetic pulse propagation interacting with a plasma is investigated, from first principles. A multiscale technique is employed to solve the fluid Maxwell equations describing weakly nonlinear circularly polarized electromagnetic pulses in magnetized plasmas. A nonlinear Schrödinger (NLS) type equation is shown to govern the amplitude of the vector potential. A set of non-stationary envelope solutions of the NLS equation are considered as potential candidates for the modeling of rogue waves (freak waves) in beam-plasma interactions, namely in the form of the Peregrine soliton, the Akhmediev breather and the Kuznetsov-Ma breather. The variation of the structural properties of the latter structures with relevant plasma parameters is investigated, in particular focusing on the ratio between the (magnetic field dependent) cyclotron (gyro-)frequency and the plasma frequency. © 2013 IOP Publishing Ltd.

Semiclassical relativistic fluid theory for electrostatic envelope modes in dense electron-positron-ion plasmas:Modulational instability and rogue waves

A fluid model is used to describe the propagation of envelope structures in an ion plasma under the influence of the action of weakly relativistic electrons and positrons. A multiscale perturbative method is used to derive a nonlinear Schrödinger equation for the envelope amplitude. Criteria for modulational instability, which occurs for small values of the carrier wavenumber (long carrier wavelengths), are derived. The occurrence of rogue waves is briefly discussed. © Cambridge University Press 2013.

Freak waves and electrostatic wavepacket modulation in a quantum electron-positron-ion plasma

The occurrence of rogue waves (freak waves) associated with electrostatic wavepacket propagation in a quantum electron-positron-ion plasma is investigated from first principles. Electrons and positrons follow a Fermi-Dirac distribution, while the ions are subject to a quantum (Fermi) pressure. A fluid model is proposed and analyzed via a multiscale technique. The evolution of the wave envelope is shown to be described by a nonlinear Schrödinger equation (NLSE). Criteria for modulational instability are obtained in terms of the intrinsic plasma parameters. Analytical solutions of the NLSE in the form of envelope solitons (of the bright or dark type) and localized breathers are reviewed. The characteristics of exact solutions in the form of the Peregrine soliton, the Akhmediev breather and the Kuznetsov-Ma breather are proposed as candidate functions for rogue waves (freak waves) within the model. The characteristics of the latter and their dependence on relevant parameters (positron concentration and temperature) are investigated. © 2014 IOP Publishing Ltd.

Amplitude modulation of quantum-ion-acoustic wavepackets in electron-positron-ion plasmas: Modulational instability, envelope modes, extreme waves

A semirelativistic fluid model is employed to describe the nonlinear amplitude modulation of low-frequency (ionic scale) electrostatic waves in an unmagnetized electron-positron-ion plasma. Electrons and positrons are assumed to be degenerated and inertialess, whereas ions are warm and classical. A multiscale perturbation method is used to derive a nonlinear Schrödinger equation for the envelope amplitude, based on which the occurrence of modulational instability is investigated in detail. Various types of localized ion acoustic excitations are shown to exist, in the form of either bright type envelope solitons (envelope pulses) or dark-type envelope solitons (voids, holes). The plasma configurational parameters (namely, the relativistic degeneracy parameter, the positron concentration, and the ionic temperature) are shown to affect the conditions for modulational instability significantly, in fact modifying the associated threshold as well as the instability growth rate. In particular, the relativistic degeneracy parameter leads to an enhancement of the modulational instability mechanism. Furthermore, the effect of different relevant plasma parameters on the characteristics (amplitude, width) of these envelope solitary structures is also presented in detail. Finally, the occurrence of extreme amplitude excitation (rogue waves) is also discussed briefly. Our results aim at elucidating the formation and dynamics of nonlinear electrostatic excitations in superdense astrophysical regimes.

Control of ionization and dissociation of H2+ by elliptically polarized ultra-short VUV laser pulses

Resonance-enhanced multiphoton ionization of H₂ ⁺ exposed to elliptically polarized VUV laser pulses is investigated. Differential cross sections for nuclei and electron are obtained using numerical solutions of the time-dependent Schrödinger equation. In this work in progress, we explore the dependence of the dissociative ionization observables with the polarization of the light.