SIMLA: Simulating particle dynamics in intense laser and other electromagnetic fields via classical and quantum electrodynamics

We present the Fortran program SIMLA, which is designed for the study of charged particle dynamics in laser and other background fields. The dynamics can be determined classically via the Lorentz force and Landau–Lifshitz equations or, alternatively, via the simulation of photon emission events determined by strong-field quantum-electrodynamics amplitudes and implemented using Monte-Carlo routines. Multiple background fields can be included in the simulation and, where applicable, the propagation direction, field type (plane wave, focussed paraxial, constant crossed, or constant magnetic), and time envelope of each can be independently specified.

Quantum-information processing for a coherent superposition state via a mixed entangled coherent channel

An entangled two-mode coherent state is studied within the framework of 2 x 2-dimensional Hilbert space. An entanglement concentration scheme based on joint Bell-state measurements is worked out. When the entangled coherent state is embedded in vacuum environment, its entanglement is degraded but not totally lost. It is found that the larger the initial coherent amplitude, the faster entanglement decreases. We investigate a scheme to teleport a coherent superposition state while considering a mixed quantum channel. We find that the decohered entangled coherent state may be useless for quantum teleportation as it gives the optimal fidelity of teleportation less than the classical limit 2/3.

Non-local quantum gates: a cavity-quantum-electrodynamics implementation

The problems related to the management of large quantum registers could be handled in the context of distributed quantum computation: unitary non-local transformations among spatially separated local processors are realized performing local unitary transformations and exchanging classical communication. In this paper, a scheme is proposed for the implementation of universal non-local quantum gates such as a controlled NOT (CNOT) and a controlled quantum phase gate (CQPG). The system chosen for their physical implementation is a cavity-quantum-electrodynamics (CQED) system formed by two spatially separated microwave cavities and two trapped Rydberg atoms. The procedures to follow for the realization of each step necessary to perform a specific non-local operation are described.

Quantum manifestations of closed orbits in the photoexcitation scaled spectrum of the hydrogen atom in crossed fields

The scaled photoexcitation spectrum of the hydrogen atom in crossed electric and magnetic fields has been obtained by means of accurate quantum mechanical calculation using a new algorithm. Closed orbits in the corresponding classical system have also been obtained, using a new, efficient and practical searching procedure. Two new classes of closed orbit have been identified. Fourier transforming each photoexcitation quantum spectrum to yield a plot against scaled action has allowed direct comparison between peaks in such plots and the scaled action values of closed orbits, Excellent agreement has been found with all peaks assigned.

Plane wave and Coulomb asymptotics

A simple plane wave solution of the Schrodinger-Helmholtz equation is a quantum eigenfunction obeying both energy and linear momentum correspondence principles. Inclusion of the outgoing wave with scattering amplitude f asymptotic development of the plane wave, we show that there is a problem with angular momentum when we consider forward scattering at the point of closest approach and at large impact parameter given semiclassically by (l + 1/2)/k where l is the azimuthal quantum number and may be large (J. Leech et al., Phys. Rev. Lett. 88. 257901 (2002)). The problem is resolved via non- uniform, non-standard analysis involving the Heaviside step function, unifying classical, semiclassical and quantum mechanics, and the treatment is extended to the case of pure Coulomb scattering.

Non-linear conductance of disordered quantum wires

The self-consistent electron potential in a current-carrying disordered quantum wire is spatially inhomogeneous due to the formation of resistivity dipoles across scattering centres. In this paper it is argued that these inhomogeneities in the potential result in a suppression of the differential conductance of such a wire at finite applied voltage. A semi-classical argument allows this suppression, quadratic in the voltage, to be related directly to the amount of intrinsic defect scattering in the wire. This result is then tested against numerical calculations.

Reducing quantum control for spin-spin entanglement distribution

We present a protocol that sets maximum stationary entanglement between remote spins through scattering of mobile mediators without initialization, post-selection or feedback of the mediators' state. No time-resolved tuning is needed and, counterintuitively, the protocol generates two-qubit singlet states even when classical mediators are used. The mechanism responsible for this effect is resilient against non-optimal coupling strengths and dephasing affecting the spins. The scheme uses itinerant particles and scattering centres and can be implemented in various settings. When quantum dots and photons are used a striking result is found: injection of classical mediators, rather than quantum ones, improves the scheme efficiency.

Noise resilience and entanglement evolution in two nonequivalent classes of quantum algorithms

The speedup provided by quantum algorithms with respect to their classical counterparts is at the origin of scientific interest in quantum computation. However, the fundamental reasons for such a speedup are not yet completely understood and deserve further attention. In this context, the classical simulation of quantum algorithms is a useful tool that can help us in gaining insight. Starting from the study of general conditions for classical simulation, we highlight several important differences between two nonequivalent classes of quantum algorithms. We investigate their performance under realistic conditions by quantitatively studying their resilience with respect to static noise. This latter refers to errors affecting the initial preparation of the register used to run an algorithm. We also compare the evolution of the entanglement involved in the different computational processes.

Quasi-classical model of non-destructive wavepacket manipulation by intense ultrashort nonresonant laser pulses

A quasi-classical model (QCM) of nuclear wavepacket generation, modification and imaging by three intense ultrafast near-infrared laser pulses has been developed. Intensities in excess of 10(13) W cm(-2) are studied, the laser radiation is non-resonant and pulse durations are in the few-cycle regime, hence significantly removed from the conditions typical of coherent control and femtochemistry. The 1s sigma ground state of the D-2 precursor is projected onto the available electronic states in D-2(+) (1s sigma(g) ground and 2p sigma(u) dissociative) and D+ + D+ (Coulomb explosion) by tunnel ionization by an ultrashort 'pump' pulse, and relative populations are found numerically. A generalized non-adiabatic treatment allows the dependence of the initial vibrational population distribution on laser intensity to be calculated. The wavepacket is approximated as a classical ensemble of particles moving on the 1s sigma(g) potential energy surface (PES), and hence follow trajectories of different amplitudes and frequencies depending on the initial vibrational state. The 'control' pulse introduces a time-dependent polarization of the molecular orbital, causing the PES to be modified according to the dynamic Stark effect and the transition dipole. The trajectories adjust in amplitude, frequency and phase-offset as work is done on or by the resulting force; comparing the perturbed and unperturbed trajectories allows the final vibrational state populations and phases to be determined. The action of the 'probe' pulse is represented by a discrete internuclear boundary, such that elements of the ensemble at a larger internuclear separation are assumed to be photodissociated. The vibrational populations predicted by the QCM are compared to recent quantum simulations (Niederhausen and Thumm 2008 Phys. Rev. A 77 013404), and a remarkable agreement has been found. The applicability of this model to femtosecond and attosecond time-scale experiments is discussed and the relation to established femtochemistry and coherent control techniques are explored.

Electromagnetic pulse compression and energy localization in quantum plasmas

The evolution of the intensity of a relativistic laser beam propagating through a dense quantum plasma is investigated, by considering different plasma regimes. A cold quantum fluid plasma and then a thermal quantum description(s) is (are) adopted, in comparison with the classical case of reference. Considering a Gaussian beam cross-section, we investigate both the longitudinal compression and lateral/longitudinal localization of the intensity of a finite-radius electromagnetic pulse. By employing a quantum plasma fluid model in combination with Maxwell's equations, we rely on earlier results on the quantum dielectric response, to model beam-plasma interaction. We present an extensive parametric investigation of the dependence of the longitudinal pulse compression mechanism on the electron density in cold quantum plasmas, and also study the role of the Fermi temperature in thermal quantum plasmas. Our numerical results show pulse localization through a series of successive compression cycles, as the pulse propagates through the plasma. A pulse of 100 fs propagating through cold quantum plasma is compressed to a temporal size of approximate to 1.35 attosecond and a spatial size of approximate to 1.08 10(-3) cm. Incorporating Fermi pressure via a thermal quantum plasma model is shown to enhance localization effects. A 100 fs pulse propagating through quantum plasma with a Fermi temperature of 350 K is compressed to a temporal size of approximate to 0.6 attosecond and a spatial size of approximate to 2.4 10(-3) cm. (c) 2010 Elsevier B.V. All rights reserved.

Thermal and quantum fluctuations in chains of ultracold polar molecules

Ultracold polar molecules, in highly anisotropic traps and interacting via a repulsive dipolar potential, may form one-dimensional chains at high densities. According to classical theory, at low temperatures there exists a critical value of the density at which a second-order phase transition from a linear to a zigzag chain occurs. We study the effect of thermal and quantum fluctuations on these self-organized structures using classical and quantum Monte Carlo methods, by means of which we evaluate the pair correlation function and the static structure factor. Depending on the parameters, these functions exhibit properties typical of a crystalline or of a liquid system. We compare the thermal and the quantum results, identifying analogies and differences. Finally, we discuss experimental parameter regimes where the effects of quantum fluctuations on the linear-zigzag transition can be observed.

Nonlinear electrostatic excitations of charged dust in degenerate ultra-dense quantum dusty plasmas

The linear and nonlinear properties of low-frequency electrostatic excitations of charged dust particles (or defects) in a dense collisionless, unmagnetized Thomas-Fermi plasma are investigated. A fully ionized three-component model plasma consisting of electrons, ions, and negatively charged massive dust grains is considered. Electrons and ions are assumed to be in a degenerate quantum state, obeying the Thomas-Fermi density distribution, whereas the inertial dust component is described by a set of classical fluid equations. Considering large-amplitude stationary profile travelling-waves in a moving reference frame, the fluid evolution equations are reduced to a pseudo-energy-balance equation, involving a Sagdeev-type potential function. The analysis describes the dynamics of supersonic dust-acoustic solitary waves in Thomas-Fermi plasmas, and provides exact predictions for their dynamical characteristics, whose dependence on relevant parameters (namely, the ion-to-electron Fermi temperature ratio, and the dust concentration) is investigated. An alternative route is also adopted, by assuming weakly varying small-amplitude disturbances off equilibrium, and then adopting a multiscale perturbation technique to derive a Korteweg–de Vries equation for the electrostatic potential, and finally solving in terms for electric potential pulses (electrostatic solitons). A critical comparison between the two methods reveals that they agree exactly in the small-amplitude, weakly superacoustic limit. The dust concentration (Havnes) parameter h = Zd0nd0/ne0 affects the propagation characteristics by modifying the phase speed, as well as the electron/ion Fermi temperatures. Our results aim at elucidating the characteristics of electrostatic excitations in dust-contaminated dense plasmas, e.g., in metallic electronic devices, and also arguably in supernova environments, where charged dust defects may occur in the quantum plasma regime.

Statistical-mechanical description of classical test-particle dynamics in the presence of an external force field: modelling noise and damping from first principles

Aiming to establish a rigorous link between macroscopic random motion (described e.g. by Langevin-type theories) and microscopic dynamics, we have undertaken a kinetic-theoretical study of the dynamics of a classical test-particle weakly coupled to a large heat-bath in thermal equilibrium. Both subsystems are subject to an external force field. From the (time-non-local) generalized master equation a Fokker-Planck-type equation follows as a "quasi-Markovian" approximation. The kinetic operator thus defined is shown to be ill-defined; in specific, it does not preserve the positivity of the test-particle distribution function f(x, v; t). Adopting an alternative approach, previously introduced for quantum open systems, is proposed to lead to a correct kinetic operator, which yields all the expected properties. A set of explicit expressions for the diffusion and drift coefficients are obtained, allowing for modelling macroscopic diffusion and dynamical friction phenomena, in terms of an external field and intrinsic physical parameters.

Almost all quantum states have nonclassical correlations

Quantum discord quantifies nonclassical correlations in a quantum system including those not captured by entanglement. Thus, only states with zero discord exhibit strictly classical correlations. We prove that these states are negligible in the whole Hilbert space: typically a state picked out at random has positive discord and, given a state with zero discord, a generic arbitrarily small perturbation drives it to a positive-discord state. These results hold for any Hilbert-space dimension and have direct implications for quantum computation and for the foundations of the theory of open systems. In addition, we provide a simple necessary criterion for zero quantum discord. Finally, we show that, for almost all positive-discord states, an arbitrary Markovian evolution cannot lead to a sudden, permanent vanishing of discord.

Local temperature in quantum thermal states

We consider blocks of quantum spins in a chain at thermal equilibrium, focusing on their properties from a thermodynamical perspective. In a classical system the temperature behaves as an intensive magnitude, above a certain block size, regardless of the actual value of the temperature itself. However, a deviation from this behavior is expected in quantum systems. In particular, we see that under some conditions the description of the blocks as thermal states with the same global temperature as the whole chain fails. We analyze this issue by employing the quantum fidelity as a figure of merit, singling out in detail the departure from the classical behavior. As it may be expected, we see that quantum features are more prominent at low temperatures and are affected by the presence of zero-temperature quantum phase transitions. Interestingly, we show that the blocks can be considered indeed as thermal states with a high fidelity, provided an effective local temperature is properly identified. Such a result may originate from typical properties of reduced subsystems of energy-constrained Hilbert spaces. Finally, the relation between local and global temperatures is analyzed as a function of the size of the blocks and the system parameters.