Source author record

Emil A. Yuzbashyan

Emil A. Yuzbashyan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Breakdown of the Migdal-Eliashberg theory and a theory of lattice-fermionic superfluidity

We show that the Migdal-Eliashberg theory loses validity at a finite value $λ_c$ of the electron-phonon coupling $λ$ regardless of the underlying model Hamiltonian. The value of $λ_c$ is approximately between 3.0 and 3.7. The new phase that emerges at $λ>λ_c$ breaks the lattice translational symmetry. Depending on the filling fraction and crystal symmetry, it is an insulator or a Fermi liquid. Its characteristic feature is a gap or a pronounced depression of the fermionic density of states near the Fermi level. We establish the breakdown from within the Migdal-Eliashberg theory by demonstrating that the normal state specific heat is negative for $λ\ge 3.7$ and the quasiparticle lifetime vanishes in the strong coupling limit. At fixed $λ>λ_c$, the transition to the new phase occurs at a critical temperature higher than the superconducting transition temperature. In addition, there is a first order phase transition between the new phase and the superconducting state as we vary $λ$ across $λ_c$ at fixed temperature. We put forward a new theory - lattice-fermionic theory of superfluidity - that bridges the gap between the Migdal-Eliashberg approach and the physics at stronger coupling. At small $λ$, our theory reduces to the Migdal-Eliashberg theory and, past $λ_c$, it describes the new phase and a range of other phenomena.

preprint2022arXiv

Migdal-Eliashberg theory as a classical spin chain

We formulate the Migdal-Eliashberg theory of electron-phonon interactions in terms of classical spins by mapping the free energy to a Heisenberg spin chain in a Zeeman magnetic field. Spin components are energy-integrated normal and anomalous Green's functions and sites of the chain are fermionic Matsubara frequencies. The Zeeman field grows linearly with the spin coordinate and competes with ferromagnetic spin-spin interaction that falls off as the square of the inverse distance. The spin-chain representation makes a range of previously unknown properties plain to see. In particular, infinitely many new solutions of the Eliashberg equations both in the normal and superconducting states emerge at strong coupling. These saddle points of the free-energy functional correspond to spin flips. We argue that they are also fixed points of kinetic equations and play an essential role in far from equilibrium dynamics of strongly coupled superconductors. Up to an overall phase, the frequency-dependent gap function that minimizes the free energy must be non-negative. There are strong parallels between our Eliashberg spins and Anderson pseudospins, though the two sets of spins never coincide.

preprint2022arXiv

Simulating dynamical phases of chiral $p+ i p$ superconductors with a trapped ion magnet

Two-dimensional $p+ i p$ superconductors and superfluids are systems that feature chiral behavior emerging from the Cooper pairing of electrons or neutral fermionic atoms with non-zero angular momentum. Their realization has been a longstanding goal because they offer great potential utility for quantum computation and memory. However, they have so far eluded experimental observation both in solid state systems as well as in ultracold quantum gases. Here, we propose to leverage the tremendous control offered by rotating two-dimensional trapped-ion crystals in a Penning trap to simulate the dynamical phases of two-dimensional $p+ip$ superfluids. This is accomplished by mapping the presence or absence of a Cooper pair into an effective spin-1/2 system encoded in the ions' electronic levels. We show how to infer the topological properties of the dynamical phases, and discuss the role of beyond mean-field corrections. More broadly, our work opens the door to use trapped ion systems to explore exotic models of topological superconductivity and also paves the way to generate and manipulate skyrmionic spin textures in these platforms.

preprint2022arXiv

Singularities in the Loschmidt echo of quenched topological superconductors

We study the Loschmidt echo in the quenched two-dimensional $p$-wave topological superconductor. We find that if this superconductor is quenched out of the critical point separating its topological and non-topological phases into either of the two gapful phases, its Loschmidt echo features singularities occurring periodically in time where the second derivative of the Loschmidt echo over time diverges logarithmically. Conversely, we give arguments towards $s$-wave superconductors not having singularities in their Loschmidt echo regardless of the quench. We also demonstrate that the conventional mean field theory calculates classical echo instead of its quantum counterpart, and show how it should be modified to capture the full quantum Loschmidt echo.

preprint2022arXiv

Superconductivity near a quantum critical point in the extreme retardation regime

We study fermions at quantum criticality with extremely retarded interactions of the form $V(ω_l)=(g/|ω_l|)^γ$, where $ω_l$ is the transferred Matsubara frequency. This system undergoes a normal-superconductor phase transition at a critical temperature $T=T_c$. The order parameter is the frequency-dependent gap function $Δ(ω_n)$ as in the Eliashberg theory. In general, the interaction is extremely retarded for $γ\gg 1$, except at low temperatures $γ>3$ is sufficient. We evaluate the normal state specific heat, $T_c$, the jump in the specific heat, $Δ(ω_n)$ near $T_c$, and the Landau free energy. Our answers are asymptotically exact in the limit $γ\to\infty$. At low temperatures, we prove that the global minimum of the free energy is nondegenerate and determine the order parameter, the free energy, and the specific heat. These answers are exact for $T\to0$ and $γ>3$. We also uncover and investigate an instability of the $γ$ model: negative specific heat at $T\to0$ and just above $T_c$.

preprint2016arXiv

Integrable matrix theory: Level statistics

We study level statistics in ensembles of integrable $N\times N$ matrices linear in a real parameter $x$. The matrix $H(x)$ is considered integrable if it has a prescribed number $n>1$ of linearly independent commuting partners $H^i(x)$ (integrals of motion) $\left[H(x),H^i(x)\right] = 0$, $\left[H^i(x), H^j(x)\right]$ = 0, for all $x$. In a recent work, we developed a basis-independent construction of $H(x)$ for any $n$ from which we derived the probability density function, thereby determining how to choose a typical integrable matrix from the ensemble. Here, we find that typical integrable matrices have Poisson statistics in the $N\to\infty$ limit provided $n$ scales at least as $\log{N}$; otherwise, they exhibit level repulsion. Exceptions to the Poisson case occur at isolated coupling values $x=x_0$ or when correlations are introduced between typically independent matrix parameters. However, level statistics cross over to Poisson at $ \mathcal{O}(N^{-0.5})$ deviations from these exceptions, indicating that non-Poissonian statistics characterize only subsets of measure zero in the parameter space. Furthermore, we present strong numerical evidence that ensembles of integrable matrices are stationary and ergodic with respect to nearest neighbor level statistics.

preprint2016arXiv

Integrals of motion for one-dimensional Anderson localized systems

Anderson localization is known to be inevitable in one dimension for generic disordered models. Since localization leads to Poissonian energy level statistics, we ask if localized systems possess "additional" integrals of motion as well, so as to enhance the analogy with quantum integrable systems. We answer this in the affirmative in the present work. We construct a set of nontrivial integrals of motion for Anderson localized models, in terms of the original creation and annihilation operators. These are found as a power series in the hopping parameter. The recently found Type-1 Hamiltonians, which are known to be quantum integrable in a precise sense, motivate our construction. We note that these models can be viewed as disordered electron models with infinite-range hopping, where a similar series truncates at the linear order. We show that despite the infinite range hopping, all states but one are localized. We also study the conservation laws for the disorder free Aubry-Andre model, where the states are either localized or extended, depending on the strength of a coupling constant. We formulate a specific procedure for averaging over disorder, in order to examine the convergence of the power series. Using this procedure in the Aubry-Andre model, we show that integrals of motion given by our construction are well-defined in localized phase, but not so in the extended phase. Finally, we also obtain the integrals of motion for a model with interactions to lowest order in the interaction.

preprint2016arXiv

Rotationally invariant ensembles of integrable matrices

We construct ensembles of random integrable matrices with any prescribed number of nontrivial integrals and formulate integrable matrix theory (IMT) -- a counterpart of random matrix theory (RMT) for quantum integrable models. A type-M family of integrable matrices consists of exactly N-M independent commuting N-by-N matrices linear in a real parameter. We first develop a rotationally invariant parametrization of such matrices, previously only constructed in a preferred basis. For example, an arbitrary choice of a vector and two commuting Hermitian matrices defines a type-1 family and vice versa. Higher types similarly involve a random vector and two matrices. The basis-independent formulation allows us to derive the joint probability density for integrable matrices, in a manner similar to the construction of Gaussian ensembles in the RMT.

preprint2016arXiv

Superconductivity near a quantum-critical point: The special role of the first Matsubara frequency

Near a quantum-critical point in a metal strong fermion-fermion interaction mediated by a soft collective boson gives rise to incoherent, non-Fermi liquid behavior. It also often gives rise to superconductivity which masks the non-Fermi liquid behavior. We analyze the interplay between superconductivity and fermionic incoherence for a set of quantum-critical models with effective dynamical interaction between low-energy fermions. We argue that superconductivity always develops above a quantum-critical point, no matter how strong the fermionic self-energy is. We argue that superconductivity should not be viewed as a pairing of incoherent fermions, as previously thought. Rather, $T_c$ is non-zero due to the fact that the self-energy is suppressed at the two lowest fermionic Matsubara frequencies $ω_m = \pm πT$. We obtain the analytic formula for $T_c$ which reproduces well numerical results for the electron-phonon model at vanishing Debye frequency.

preprint2015arXiv

Generalized microcanonical and Gibbs ensembles in classical and quantum integrable dynamics

We prove two statements about the long time dynamics of integrable Hamiltonian systems. In classical mechanics, we prove the microcanonical version of the Generalized Gibbs Ensemble (GGE) by mapping it to a known theorem and then extend it to the limit of infinite number of degrees of freedom. In quantum mechanics, we prove GGE for maximal Hamiltonians - a class of models stemming from a rigorous notion of quantum integrability understood as the existence of conserved charges with prescribed dependence on a system parameter, e.g. Hubbard U, anisotropy in the XXZ model etc. In analogy with classical integrability, the defining property of these models is that they have the maximum number of independent integrals. We contrast their dynamics induced by quenching the parameter to that of random matrix Hamiltonians.

preprint2014arXiv

Quantum integrability in the multistate Landau-Zener problem

We analyze Hamiltonians linear in the time variable for which the multistate Landau-Zener problem is known to have an exact solution. We show that they either belong to families of mutually commuting Hamiltonians polynomial in time or reduce to the 2 x 2 Landau-Zener problem, which is considered trivially integrable. The former category includes the equal slope, bow-tie, and generalized bow-tie models. For each of these models we explicitly construct the corresponding families of commuting matrices. The equal slope model is a member of an integrable family that consists of the maximum possible number (for a given matrix size) of commuting matrices linear in time. The bow-tie model belongs to a previously unknown, similarly maximal family of quadratic commuting matrices. We thus conjecture that quantum integrability understood as the existence of nontrivial parameter-dependent commuting partners is a necessary condition for the Landau-Zener solvability. Descendants of the 2 x 2 Landau-Zener Hamiltonian are e.g. general SU(2) and SU(1,1) Hamiltonians, time-dependent linear chain, linear, nonlinear, and double oscillators. We explicitly obtain solutions to all these Landau-Zener problems from the 2 x 2 case.

preprint2014arXiv

Quantum quench phase diagrams of an s-wave BCS-BEC condensate

We study the dynamic response of an s-wave BCS-BEC (atomic-molecular) condensate to detuning quenches within the two channel model beyond the weak coupling BCS limit. At long times after the quench, the condensate ends up in one of three main asymptotic states (nonequilibrium phases), which are qualitatively similar to those in other fermionic condensates defined by a global complex order parameter. In phase I the amplitude of the order parameter vanishes as a power law, in phase II it goes to a nonzero constant, and in phase III it oscillates persistently. We construct exact quench phase diagrams that predict the asymptotic state (including the many-body wavefunction) depending on the initial and final detunings and on the Feshbach resonance width. Outside of the weak coupling regime, both the mechanism and the time dependence of the relaxation of the amplitude of the order parameter in phases I and II are modified. Also, quenches from arbitrarily weak initial to sufficiently strong final coupling do not produce persistent oscillations in contrast to the behavior in the BCS regime. The most remarkable feature of coherent condensate dynamics in various fermion superfluids is an effective reduction in the number of dynamic degrees of freedom as the evolution time goes to infinity. As a result, the long time dynamics can be fully described in terms of just a few new collective dynamical variables governed by the same Hamiltonian only with "renormalized" parameters. Combining this feature with the integrability of the underlying (e.g. the two channel) model, we develop and consistently present a general method that explicitly obtains the exact asymptotic state of the system.

preprint2014arXiv

Quench-induced Floquet topological p-wave superfluids

Ultracold atomic gases in two dimensions tuned close to a p-wave Feshbach resonance were expected to exhibit topological superfluidity, but these were found to be experimentally unstable. We show that one can induce a topological Floquet superfluid if weakly interacting atoms are brought suddenly close ("quenched") to such a resonance, in the time before the instability kicks in. The resulting superfluid possesses Majorana edge modes, yet differs from a conventional Floquet system as it is not driven externally. Instead, the periodic modulation is self-generated by the dynamics.

preprint2013arXiv

Quantum integrability in systems with finite number of levels

We consider the problem of defining quantum integrability in systems with finite number of energy levels starting from commuting matrices and construct new general classes of such matrix models with a given number of commuting partners. We argue that if the matrices depend on a (real) parameter, one can define quantum integrability from this feature alone, leading to specific results such as exact solvability, Poissonian energy level statistics and to level crossings.

preprint2013arXiv

Quantum quench in a p + i p superfluid: Winding numbers and topological states far from equilibrium

We study the non-adiabatic dynamics of a 2D p+ip superfluid following a quantum quench of the BCS coupling constant. The model describes a topological superconductor with a non-trivial BCS (trivial BEC) phase appearing at weak (strong) coupling strengths. We extract the exact long-time asymptotics of the order parameter Δ(t) by exploiting the integrability of the classical p-wave Hamiltonian, which we establish via a Lax construction. Three different types of behavior can occur depending upon the strength and direction of the quench. In phase I, the order parameter asymptotes to zero. In phase II, Δ(t) goes to a non-zero constant. Phase III is characterized by persistent oscillations of Δ(t). For quenches within I and II, we determine the topological character of the asymptotic states. We show that two different formulations of the bulk topological winding number, although equivalent in the ground state, must be regarded as independent out of equilibrium. The first number Q characterizes the Anderson pseudospin texture of the initial state; we show that it is conserved. For non-zero Q, this leads to the prediction of a "gapless topological" state when Δ(t) goes to zero. The presence or absence of Majorana edge modes in a sample with a boundary is encoded in the second winding number W, formulated in terms of the retarded Green's function. We show that W can change following a quench across the quantum critical point. We discuss the implications for the (dis)appearance of Majorana edge modes. Finally, we show that the parity of zeros in the bulk out-of-equilibrium Cooper pair distribution function constitutes a Z2-valued quantum number, which is non-zero whenever W differs from Q. The pair distribution can in principle be measured using RF spectroscopy in an ultracold atom realization, allowing direct experimental detection of the bulk Z2 number.

preprint2012arXiv

Interaction-mediated surface state instability in disordered three-dimensional topological superconductors with spin SU(2) symmetry

We show that arbitrarily weak interparticle interactions destabilize the surface states of 3D topological superconductors with spin SU(2) invariance (symmetry class CI), in the presence of non-magnetic disorder. The conduit for the instability is disorder-induced wavefunction multifractality. We argue that time-reversal symmetry breaks spontaneously at the surface, so that topologically-protected states do not exist for this class. The interaction-stabilized surface phase is expected to exhibit ferromagnetic order, or to reside in an insulating plateau of the spin quantum Hall effect.

preprint2011arXiv

BCS superconductivity in metallic nanograins: Finite-size corrections, low energy excitations, and robustness of shell effects

We combine the BCS self-consistency condition, a semiclassical expansion for the spectral density and interaction matrix elements to describe analytically how the superconducting gap depends on the size and shape of a 2d and 3d superconducting grain. In chaotic grains mesoscopic fluctuations of the matrix elements lead to a smooth dependence of the order parameter on the excitation energy. In the integrable case we find shell effects i. e. for certain values of the electron number $N$ a small change in N leads to large changes in the energy gap. With regard to possible experimental tests we provide a detailed analysis of the dependence of the gap on the coherence length and the robustness of shell effects under small geometrical deformations.

preprint2011arXiv

Classification of Parameter-Dependent Quantum Integrable Models, Their Parameterization, Exact Solution, and Other Properties

We study general quantum integrable Hamiltonians linear in a coupling constant and represented by finite NxN real symmetric matrices. The restriction on the coupling dependence leads to a natural notion of nontrivial integrals of motion and classification of integrable families into Types according to the number of such integrals. A Type M family in our definition is formed by N-M nontrivial mutually commuting operators linear in the coupling. Working from this definition alone, we parameterize Type M operators, i.e. resolve the commutation relations, and obtain an exact solution for their eigenvalues and eigenvectors. We show that our parameterization covers all Type 1, 2, and 3 integrable models and discuss the extent to which it is complete for other types. We also present robust numerical observation on the number of energy level crossings in Type M integrable systems and analyze the taxonomy of types in the 1d Hubbard model.

preprint2011arXiv

Quantum quench spectroscopy of a Luttinger liquid: Ultrarelativistic density wave dynamics due to fractionalization in an XXZ chain

We compute the dynamics of localized excitations produced by a quantum quench in the spin 1/2 XXZ chain. Using numerics combining the density matrix renormalization group and exact time evolution, as well as analytical arguments, we show that fractionalization due to interactions in the pre-quench state gives rise to "ultrarelativistic" density waves that travel at the maximum band velocity. The system is initially prepared in the ground state of the chain within the gapless XY phase, which admits a Luttinger liquid (LL) description at low energies and long wavelengths. The Hamiltonian is then suddenly quenched to a band insulator, after which the chain evolves unitarily. Through the gapped dispersion of the insulator spectrum, the post-quench dynamics serve as a "velocity microscope," revealing initial state particle correlations via space time density propagation. We show that the ultrarelativistic wave production is tied to the particular way in which fractionalization evades Pauli-blocking in the zero-temperature initial LL state.

preprint2010arXiv

Quantum quench in 1D: Coherent inhomogeneity amplification and 'supersolitons'

We study a quantum quench in a 1D system possessing Luttinger liquid (LL) and Mott insulating ground states before and after the quench, respectively. We show that the quench induces power law amplification in time of any particle density inhomogeneity in the initial LL ground state. The scaling exponent is set by the fractionalization of the LL quasiparticle number relative to the insulator. As an illustration, we consider the traveling density waves launched from an initial localized density bump. While these waves exhibit a particular rigid shape, their amplitudes grow without bound.

preprint2008arXiv

BCS theory for finite size superconductors

We study finite size effects in superconducting metallic grains and determine the BCS order parameter and the low energy excitation spectrum in terms of size, and shape of the grain. Our approach combines the BCS self-consistency condition, a semiclassical expansion for the spectral density and interaction matrix elements, and corrections to the BCS mean-field. In chaotic grains mesoscopic fluctuations of the matrix elements lead to a smooth dependence of the order parameter on the excitation energy. In the integrable case we observe shell effects when e.g. a small change in the electron number leads to large changes in the energy gap.

preprint2005arXiv

Integrable dynamics of coupled Fermi-Bose condensates

We study the mean-field dynamics of a fermionic condensate interacting with a single bosonic mode (a generalized Dicke model). This problem is integrable and can be mapped onto a corresponding BCS problem. We derive the general solution and a full set of integrals of motion for the time evolution of coupled Fermi-Bose condensates. The present paper complements our earlier study of the dynamics of the BCS model (cond-mat/0407501 and cond-mat/0505493). Here we provide a self-contained introduction to the variable separation method, which enables a complete analytical description of the evolution of the generalized Dicke, BCS, and other similar models.