Relativistic laser pulse compression in plasmas with a linear axial density gradient

The self-compression of a relativistic Gaussian laser pulse propagating in a non-uniform plasma is investigated. A linear density inhomogeneity (density ramp) is assumed in the axial direction. The nonlinear Schrodinger equation is first solved within a one-dimensional geometry by using the paraxial formalism to demonstrate the occurrence of longitudinal pulse compression and the associated increase in intensity. Both longitudinal and transverse self-compression in plasma is examined for a finite extent Gaussian laser pulse. A pair of appropriate trial functions, for the beam width parameter (in space) and the pulse width parameter (in time) are defined and the corresponding equations of space and time evolution are derived. A numerical investigation shows that inhomogeneity in the plasma can further boost the compression mechanism and localize the pulse intensity, in comparison with a homogeneous plasma. A 100 fs pulse is compressed in an inhomogeneous plasma medium by more than ten times. Our findings indicate the possibility for the generation of particularly intense and short pulses, with relevance to the future development of tabletop high-power ultrashort laser pulse based particle acceleration devices and associated high harmonic generation. An extension of the model is proposed to investigate relativistic laser pulse compression in magnetized plasmas.

Discretization and solution of the inhomogeneous Bogoliubov-de Gennes equations

In the presence of inhomogeneities, defects and currents, the equations describing a Bose-condensed ensemble of alkali atoms have to be solved numerically. By combining both linear and nonlinear equations within a Discrete Variable Representation framework, we describe a computational scheme for the solution of the coupled Bogoliubov-de Gennes (BdG) and nonlinear Schrodinger (NLS) equations for fields in a 3D spheroidal potential. We use the method to calculate the collective excitation spectrum and quasiparticle mode densities for excitations of a Bose condensed gas in a spheroidal trap. The method is compared against finite-difference and spectral methods, and we find the DVR computational scheme to be superior in accuracy and efficiency for the cases we consider. (C) 2004 Elsevier B.V. All rights reserved.

Electrostatic solitary waves in the presence of excess superthermal electrons: modulational instability and envelope soliton modes

The nonlinear dynamics of electrostatic solitary waves in the form of localized modulated wavepackets is investigated from first principles. Electron-acoustic (EA) excitations are considered in a two-electron plasma, via a fluid formulation. The plasma, assumed to be collisionless and uniform (unmagnetized), is composed of two types of electrons (inertial cold electrons and inertialess kappa-distributed superthermal electrons) and stationary ions. By making use of a multiscale perturbation technique, a nonlinear Schrodinger equation is derived for the modulated envelope, relying on which the occurrence of modulational instability (MI) is investigated in detail. Stationary profile localized EA excitations may exist, in the form of bright solitons (envelope pulses) or dark envelopes (voids). The presence of superthermal electrons modifies the conditions for MI to occur, as well as the associated threshold and growth rate. The concentration of superthermal electrons (i.e., the deviation from a Maxwellian electron distribution) may control or even suppress MI. Furthermore, superthermality affects the characteristics of solitary envelope structures, both qualitatively (supporting one or the other type, for different.) and quantitatively, changing their characteristics (width, amplitude). The stability of bright and dark-type nonlinear structures is confirmed by numerical simulations.

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.

Nonlinear wavepackets in pair-ion and electron-positron-ion plasmas

The occurrence of amplitude-modulated electrostatic and electromagnetic
wavepackets in pair plasmas is investigated. A static additional charged background species is considered, accounting for dust defects or for heavy ion
presence in the background. Relying on a two-fluid description, a nonlinear
Schrodinger type evolution equation is obtained and analyzed, in terms of the
slow dynamics of the wave amplitude. Exact envelope excitations are obtained,
modelling envelope pulses or holes, and their characteristics are discussed.

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]

Dynamics of dark hollow Gaussian laser pulses in relativistic plasma

Optical beams with null central intensity have potential applications in the field of atom optics. The spatial and temporal evolution of a central shadow dark hollow Gaussian (DHG) relativistic laser pulse propagating in a plasma is studied in this article for first principles. A nonlinear Schrodinger-type equation is obtained for the beam spot profile and then solved numerically to investigate the pulse propagation characteristics. As series of numerical simulations are employed to trace the profile of the focused and compressed DHG laser pulse as it propagates through the plasma. The theoretical and simulation results predict that higher-order DHG pulses show smaller divergence as they propagate and, thus, lead to enhanced energy transport. © 2013 American Physical Society.

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.

Optimization by quantum annealing: Lessons from simple cases

We investigate the basic behavior and performance of simulated quantum annealing (QA) in comparison with classical annealing (CA). Three simple one-dimensional case study systems are considered: namely, a parabolic well, a double well, and a curved washboard. The time-dependent Schrodinger evolution in either real or imaginary time describing QA is contrasted with the Fokker-Planck evolution of CA. The asymptotic decrease of excess energy with annealing time is studied in each case, and the reasons for differences are examined and discussed. The Huse-Fisher classical power law of double-well CA is replaced with a different power law in QA. The multiwell washboard problem studied in CA by Shinomoto and Kabashima and leading classically to a logarithmic annealing even in the absence of disorder turns to a power-law behavior when annealed with QA. The crucial role of disorder and localization is briefly discussed.