Entanglement Replication in Driven Dissipative Many-Body systems

We study the dissipative dynamics of two independent arrays of many-body systems, locally driven by a common entangled field. We showthat in the steady state the entanglement of the driving field is reproduced in an arbitrarily large series of inter-array entangled pairs over all distances. Local nonclassical driving thus realizes a scale-free entanglement replication and long-distance entanglement distribution mechanism that has immediate bearing on the implementation of quantum communication networks. 

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Nonequilibrium processes from generalized Langevin equations: Realistic nanoscale systems connected to two thermal baths

We extend the generalized Langevin equation (GLE) method [L. Stella, C. D. Lorenz, and L. Kantorovich, Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the dissipative central system and its surrounding baths. Following the original GLE approach, the extended Langevin dynamics scheme is modified to take into account two sets of auxiliary degrees of freedom corresponding to the mapping of the vibrational properties of each bath. These auxiliary variables are then used to solve the non-Markovian dissipative dynamics of the central region. The resulting algorithm is used to study a model of a short Al nanowire connected to two baths. The results of the simulations using the GLE-2B approach are compared to the results of other simulations that were carried out using standard thermostatting approaches (based on Markovian Langevin and Nosé-Hoover thermostats). We concentrate on the steady-state regime and study the establishment of a local temperature profile within the system. The conditions for obtaining a flat profile or a temperature gradient are examined in detail, in agreement with earlier studies. The results show that the GLE-2B approach is able to treat, within a single scheme, two widely different thermal transport regimes, i.e., ballistic systems, with no temperature gradient, and diffusive systems with a temperature gradient.

Dissipative Optomechanics in a Michelson--Sagnac Interferometer

Dissipative optomechanics studies the coupling of the motion of an optical element to the decay rate of a cavity. We propose and theoretically explore a realization of this system in the optical domain, using a combined Michelson-Sagnac interferometer, which enables a strong and tunable dissipative coupling. Quantum interference in such a setup results in the suppression of the lower motional sideband, leading to strongly enhanced cooling in the non-sideband-resolved regime. With state-of-the-art parameters, ground-state cooling and low-power quantum-limited position transduction are both possible. The possibility of a strong, tunable dissipative coupling opens up a new route towards observation of such fundamental optomechanical effects as nonlinear dynamics. Beyond optomechanics, the suggested method can be readily transferred to other setups involving nonlinear media, atomic ensembles, or single atoms.

Electrostatic shock dynamics in superthermal plasmas

The propagation of ion acoustic shocks in nonthermal plasmas is investigated, both analytically and numerically. An unmagnetized collisionless electron-ion plasma is considered, featuring a superthermal (non-Maxwellian) electron distribution, which is modeled by a ?-(kappa) distribution function. Adopting a multiscale approach, it is shown that the dynamics of low-amplitude shocks is modeled by a hybrid Korteweg-de Vries-Burgers (KdVB) equation, in which the nonlinear and dispersion coefficients are functions of the ? parameter, while the dissipative coefficient is a linear function of the ion viscosity. All relevant shock parameters are shown to depend on ?: higher deviations from a pure Maxwellian behavior induce shocks which are narrower, faster, and of larger amplitude. The stability profile of the kink-shaped solutions of the KdVB equation against external perturbations is investigated. The spatial profile of the shocks is found to depend upon the dispersion and the dissipation term, and the role of the interplay between dispersion and dissipation is elucidated.

Current-induced atomic dynamics, instabilities, and Raman signals: Quasiclassical Langevin equation approach

We derive and employ a semiclassical Langevin equation obtained from path integrals to describe the ionic dynamics of a molecular junction in the presence of electrical current. The electronic environment serves as an effective nonequilibrium bath. The bath results in random forces describing Joule heating, current-induced forces including the nonconservative wind force, dissipative frictional forces, and an effective Lorentz-type force due to the Berry phase of the nonequilibrium electrons. Using a generic two-level molecular model, we highlight the importance of both current-induced forces and Joule heating for the stability of the system. We compare the impact of the different forces, and the wide-band approximation for the electronic structure on our result. We examine the current-induced instabilities (excitation of runaway "waterwheel" modes) and investigate the signature of these in the Raman signals.

Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

Electron-scale dissipative electrostatic solitons in multi-species plasmas

The linear and nonlinear properties of small-amplitude electron-acoustic solitary waves are investigated via the fluid dynamical approach. A three-component plasma is considered, composed of hot electrons, cold electrons, and ions (considered stationary at the scale of interest). A dissipative (wave damping) effect is assumed due to electron-neutral collisions. The background (hot) electrons are characterized by an energetic (excessively superthermal) population and are thus modeled via a κ-type nonthermal distribution. The linear characteristics of electron-acoustic excitations are discussed, for different values of the plasma parameters (superthermality index κ and cold versus hot electron population concentration β). Large wavelengths (beyond a threshold value) are shown to be overdamped. The reductive perturbation technique is used to derive a dissipative Korteweg de-Vries (KdV) equation for small-amplitude electrostatic potential disturbances. These are expressed by exact solutions in the form of dissipative solitary waves, whose dynamics is investigated analytically and numerically. Our results should be useful in elucidating the behavior of space and experimental plasmas characterized by a coexistence of electron populations at different temperatures, where electron-neutral collisions are of relevance.