Search

We study the non-equilibrium dynamics of the superconducting order parameter in the Hatsugai-Kohmoto (HK) model. In the absence of superconductivity, its ground state is a non-Fermi liquid, whose properties are controlled by the HK interaction. Our protocol involves quantum quenching the HK interaction but leaving the interaction responsible for superconductivity unchanged. We map out the non-equilibrium dynamical phase diagram of the interacting model which contains three phases where, at long times, the order parameter amplitude vanishes, approaches a constant value or persistently oscillates. We also investigate the Loschmidt echo in searching for dynamical quantum phase transition, and find no significant correlation between the three phases and the non-analytic temporal behaviour of the Loschmidt echo. The momentum space entanglement entropy between positive and negative momentum modes, relevant for Coopair pairing, is calculated. Counterintuitively, this momentum space entanglement does not change significantly during the quench dynamics and its value remains reasonably large even for vanishing superconducting order parameter. Nevertheless, its derivative with respect to the HK interaction signals the dynamical phase transition associated to the late time vanishing of superconductivity., Comment: 10 pages, 8 figures

An electron tunneling across a junction integrated into an electric circuit can generate an excitation in the photonic field (electromagnetic environment) and lose energy in the process. Such inelastic tunneling of particles is commonly described using the $P(E)$ theory. In the conventional approach to this theory, the tunneling rate and the electric current through the junction are derived using Fermi's golden rule and by averaging over the environmental photonic degrees of freedom. In this work, we address the same problem of inelastic tunneling due to photonic environment in Lindbladian formalism and we present how the photonic degrees of freedom are traced out in quantum master equation approach. The resulting quantum master equation is parameterized by the same $P(E)$ function and enables us to obtain not only the electric current but various other quantities, for instance the heat current, in a systematic and convenient way. We also demonstrate that the Lindbladian formalism provides a comprehensive description of Bogoliubov quasiparticle tunneling through superconducting junctions and that it properly accounts for the coherence factors. The coherence factors become important if the normal-state density of states is particle-hole asymmetric., Comment: 12 pages, 3 figures

We present a theoretical study of a system consisting of a superconducting island and two quantum dots, a possible platform for building qubits and Cooper pair splitters, in the regime where each dot hosts a single electron and, hence, carries a magnetic moment. We focus on the case where the dots are coupled to overlapping superconductor states and we study whether the spins are ferromagnetically or antiferromagnetically aligned. We show that if the total number of particles is even, the spins align antiferromagnetically in the flatband limit, i.e., when the bandwidth of the superconductor is negligibly small, but ferromagnetically if the bandwidth is finite and above some value. If the total number of particles is odd, the alignment is ferromagnetic independently from the bandwidth. This implies that the results of the flatband limit are applicable only within restricted parameter regime for realistic superconducting qubit systems and that some care is required in applying simplified models to devices such as Cooper pair splitters., Comment: 15 pages, 7 figures

We study the non-equilibrium dynamics of a Luttinger liquid after a simultaneous quantum quench of the interaction and dissipative quench to the environment within the realm of the Lindblad equation. When the couplings to environment satisfy detailed balance, the system is destined to thermalize, which we follow using bosonization. The thermodynamic entropy of the system, which also encodes information about the entanglement with the environment, either exhibits a maximum and minimum or grows monotonically before reaching its thermal value in a universal fashion. The single-particle density matrix reveals rich behaviour in the spatio-temporal "phase diagram", including thermal, prethermal and dissipation enhanced sudden quench Luttinger liquid behaviour as well as dissipation enhanced Fermi liquid response., Comment: 8 pages, 7 figures

We study the autonomous systems of quadratic differential equations of the form $\dot{x}_i(t)=\mathbf{x}(t)^T \mathbf{A}_i \mathbf{x}(t) + \mathbf{v}_i^T \mathbf{x}(t)$ with $\mathbf{x}(t) = (x_1(t),x_2(t),\dots,x_i(t),\dots)$ which, in general, cannot be solved exactly. In the present paper, we present a subclass of analytically solvable quadratic systems, whose solution is realized through a multi-dimensional generalization of the inversion which transforms a quadratic system into a linear system. We provide a constructive algorithm which, on one hand, decides whether the system of differential equations is analytically solvable with the inversion transformation and, on the other hand, provides the solution. The presented results apply for arbitrary, finite number of variables., Comment: 11 pages, 2 figures

We investigate a quantum quench from a critical to an exceptional point. The initial state, prepared in the ground state of a critical hermitian system, is time evolved with a non-hermitian SSH model, tuned to its exceptional point. The single particle density matrix exhibits supersonic modes and multiple light cones, characteristic to non-hermitian time evolution. These propagate with integer multiples of the original Fermi velocity. In the long time limit, the fermionic Green's function decays spatially as $1/x^2$, in sharp contrast to the usual $1/x$ decay of non-interacting fermions. The entanglement entropy is understood as if all these supersonic modes arise from independent quasiparticles (though they do not), traveling with the corresponding supersonic light cone velocity. The entropy production rate decreases with time and develops plateaus during the time evolution, signaling the distinct velocities in the propagation of non-local quantum correlations. At late times, the entanglement entropy saturates to a finite value, satisfying a volume law., Comment: 10 pages, 5 figures

We investigate the stability of a Luttinger liquid, upon suddenly coupling it to a dissipative environment. Within the Lindblad equation, the environment couples to local currents and heats the quantum liquid up to infinite temperatures. The single particle density matrix reveals the fractionalization of fermionic excitations in the spatial correlations by retaining the initial non-integer power law exponents, accompanied by an exponential decay in time with interaction dependent rate. The spectrum of the time evolved density matrix is gapped, which collapses gradually as $-\ln(t)$. The von Neumann entropy crosses over from the early time $-t\ln(t)$ behaviour to $\ln(t)$ growth for late times. The early time dynamics is captured numerically by performing simulations on spinless interacting fermions, using several numerically exact methods. Our results could be tested experimentally in bosonic Luttinger liquids., Comment: 8 pages, 3 figures

We study a one-dimensional Fermi gas in the presence of dissipative coupling to environment through the Lindblad equation. The dissipation involves energy exchange with the environment and favours the relaxation of electrons to excitations. After switching on the dissipation, the system approaches a steady state, which is described by a generalized Gibbs ensemble. The fermionic single particle density matrix resembles deceivingly to that in a hermitian interaction quench. It decays inversely with the distance for short times due to the fermionic correlations in the initial state, which changes into a non-integer power law decay for late times, representing dissipation induced Luttinger liquid behaviour. However, the crossover between the two regions occurs due to dissipation induced damping, and is unrelated to the propagation of excitations. The velocity of information spreading is set by the dissipative coupling, and differs significantly from the original sound velocity. The thermodynamic entropy grows as $-t\ln t$ initially, and saturates to an extensive value. Our results can be tested experimentally in one-dimensional Dirac systems., Comment: 5 pages, 2 figures + supplementary material

Reaching a target quantum state from an initial state within a finite temporal window is a challenging problem due to non-adiabaticity. We study the optimal protocol for swithcing on interactions to reach the ground state of a weakly interacting Luttinger liquid within a finite time $\tau$, starting from the non-interacting ground state. The protocol is optimized by minimizing the excess energy at the end of the quench, or by maximizing the overlap with the interacting ground state. We find that the optimal protocol is symmetric with respect to $\tau/2$, and can be expressed as a functional of the occupation numbers of the bosonic modes in the final state. For short quench durations, the optimal protocol exhibits fast oscillation and excites high energy modes. In the limit of large $\tau$, minimizing energy requires a smooth protocol while maximizing overlap requires a linear quench protocol. In this limit, the minimal energy and maximal overlap are both universal functions of the system size and the duration of the protocol., Comment: 8 pages, 5 figures

We study the non-equilibrium dynamics of the Luttinger model after a quantum quench, when the initial state is a finite temperature thermal equilibrium state. The diagonal elements of the density matrix in the steady state show thermal features for high temperature initial states only, otherwise retain highly non-thermal character. The time evolution of Uhlmann fidelity, which measures the distance between the time evolved and initial states, is evaluated for arbitrary initial temperatures and quench protocols. In the long time limit, the overlap between the time evolved and initial system decreases exponentially with the temperature with a universal prefactor. Within perturbation theory, the statistics of final total energy and work are numerically evaluated in the case of a sudden quench, which yield identical distributions at zero temperature. In both statistics, temperature effects are more significant in small systems. The Dirac-delta peak at the adiabatic ground state energy remains present in the probability distribution of the total energy, but disappears from the work distribution at non-zero initial temperatures., Comment: 11 pages, 4 figures, accepted in PRB

We investigate optical transitions of non-interacting electron systems consisting of two symmetric energy bands touching each other at the Fermi energy (e.g. graphene at half filling). Optical conductivity is obtained using Kubo formula at zero temperature. We show that for particles whose pseudospin direction is determined solely by the direction of their momentum, the optical conductivity has power law frequency dependence with the exponent $(d-2)/z$ where $d$ is the dimension of the system and $z$ is the dynamical exponent. According to our result two-dimensional systems with the above pseudospin characteristics always exhibit frequency-independent optical conductivity.

We analyze the probability distribution function (PDF) of work done on a Luttinger liquid for an arbitrary finite duration interaction quench and show that it can be described in terms a generalized Gibbs ensemble. We construct the corresponding density matrix with explicit intermode correlations, and determine the duration and interaction dependence of the probability of an adiabatic transition and the PDF of non-adiabatic processes. In the thermodynamic limit, the PDF of work exhibits a non-Gaussian maximum around the excess heat, carrying almost all spectral weight. In contrast, in the small system limit most spectral weight is carried by a delta peak at the energy of the adiabatic process, and an oscillating PDF with dips at energies commensurate to the quench duration and with an exponential envelope develops. Relevance to cold atom experiments is also discussed., Comment: 5+3 pages, 1+2 figures, new results on the generalized Gibbs ensemble and on the details of the calculations are added

We investigate the local density of states and Friedel oscillation in graphene around a well localized impurity in Born approximation. In our analytical calculations Green's function technique has been used taking into account both the localized atomic wavefunctions in a tight-binding scheme and the corresponding symmetries of the lattice. As a result we obtained long wavelength oscillations in the density of electrons with long range behavior proportional to the inverse square of the distance from the impurity. These leading oscillations are out of phase on nearby lattice sites (in fact for an extended defect they cancel each other within one unit cell), therefore a probe with resolution worse than a few unit cells will experience only the next to leading inverse cube decay of density oscillations even for a short range scatterer., Comment: 4 pages, 1 figure

We study the quantum critical properties of antiferromagnetism in graphene at T=0 within mean-field (MF) theory. The resulting exponents differ from the conventional MF exponents, describing finite temperature transitions. Motivated by this, we have developed the MF theory of general gapless phases with density of states rho(E) |E|^r, r>-1, with the interaction as control parameter. For r>2, the conventional MF exponents \'a la Landau are recovered, while for -10, therefore no weak-coupling solution exists in this range. This generalizes the results on quantum criticality of the gapless Kondo systems to bulk correlated phases., Comment: 5 pages, 1 figure, 2 tables