All Publications

We present a study of the digitized quantum annealing protocol proposed by R. Barends et al. [Nature (London) 534, 222 (2016)NATUAS0028-083610.1038/nature17658]. Our analysis, performed on the benchmark case of a transverse Ising chain problem, shows that the algorithm has a well-defined optimal working point for the annealing time τPopt, scaling as τPopt∼P, where P is the number of digital Trotter steps, beyond which the residual energy error shoots up toward the value characteristic of the maximally disordered state. We present an analytical analysis for the translationally invariant transverse Ising chain case, but our numerical evidence suggests that this scenario is more general, surviving, for instance, the presence of disorder.

We show that a hierarchy of topological phases in one dimension - a topological Devil's staircase - can emerge at fractional filling fractions in interacting systems, whose single-particle band structure describes a topological or a crystalline topological insulator. Focusing on a specific example in the BDI class, we present a field-theoretical argument based on bosonization that indicates how the system, as a function of the filling fraction, hosts a series of density waves. Subsequently, based on a numerical investigation of the low-lying energy spectrum, Wilczek-Zee phases, and entanglement spectra, we show that they are symmetry protected topological phases. In sharp contrast to the non-interacting limit, these topological density waves do not follow the bulk-edge correspondence, as their edge modes are gapped. We then discuss how these results are immediately applicable to models in the AIII class, and to crystalline topological insulators protected by inversion symmetry. Our findings are immediately relevant to cold atom experiments with alkaline-earth atoms in optical lattices, where the band structure properties we exploit have been recently realized.

We discuss and qualify the connection between two separate phenomena in the physics of nanoscale friction, general in nature and relevant to experiments. The first is thermally assisted lubricity (TAL), i.e., the low-velocity regime where a nanosized dry slider exhibits a viscouslike friction despite corrugations that would otherwise imply hard stick-slip friction. The second is the occurrence of negative dissipated work (NDW) events in sampling the work probability distribution. The abundance, or scarcity due to insufficient sampling, of these NDW events implies experimental fulfillment or violation of the celebrated Jarzynski equality (JE) of nonequilibrium statistical mechanics. We show, both analytically and in simulations of the one-dimensional point slider Prandtl-Tomlinson model, that a general crossover can be individuated as the total frictional work per cycle crosses kBT. Below such crossover, the TAL regime holds, the dissipation is essentially linear, and the numerical validation for the JE is feasible (i.e., does not require an exponentially large sampling size). Above it, the dissipation profile departs from linearity and gains its hard stick-slip features, and the mandatory sampling for the JE becomes exponentially large. In addition, we derive a parameter-free formula expressing the linear velocity coefficient of viscous friction, correcting previous empirically parameterized expressions. With due caution, the connection between friction and work tails can be extended beyond a single degree of freedom to more complex sliders, thus inviting realistic crosscheck experiments. Of importance for experimental nanofriction will be the search for NDW tails in the sliding behavior of trapped cold ions, and alternatively checking for TAL in the sliding pattern of dragged colloid monolayers as well as in forced protein unwinding.

Simulated quantum annealing (SQA) is a classical computational strategy that emulates a quantum annealing (QA) dynamics through a path-integral Monte Carlo whose parameters are changed during the simulation. Here we apply SQA to the one-dimensional transverse field Ising chain, where previous works have shown that, in the presence of disorder, a coherent QA provides a quadratic speedup with respect to classical simulated annealing, with a density of Kibble-Zurek defects decaying as ρKZQA∼(log10τ)-2 as opposed to ρKZSA∼(log10τ)-1, τ being the total annealing time, while for the ordered case both give the same power law ρKZQA≈ρKZSA∼τ-1/2. We show that the dynamics of SQA, while correctly capturing the Kibble-Zurek scaling τ-1/2 for the ordered case, is unable to reproduce the QA dynamics in the disordered case at intermediate τ. We analyze and discuss several issues related to the choice of the Monte Carlo moves (local or global in space), the time-continuum limit needed to eliminate the Trotter-discretization error, and the long autocorrelation times shown by a local-in-space Monte Carlo dynamics for large disordered samples.

We investigate the use of variational wave functions that mimic stochastic recurrent neural networks, specifically, unrestricted Boltzmann machines, as guiding functions in projective quantum Monte Carlo (PQMC) simulations of quantum spin models. As a preliminary step, we investigate the accuracy of such unrestricted neural network states as variational Ansätze for the ground state of the ferromagnetic quantum Ising chain. We find that by optimizing just three variational parameters, independently on the system size, accurate ground-state energies are obtained, comparable to those previously obtained using restricted Boltzmann machines with few variational parameters per spin. Chiefly, we show that if one uses optimized unrestricted neural network states as guiding functions for importance sampling, the efficiency of the PQMC algorithms is greatly enhanced, drastically reducing the most relevant systematic bias, namely, the one due to the finite random-walker population. The scaling of the computational cost with the system size changes from the exponential scaling characteristic of PQMC simulations performed without importance sampling, to a polynomial scaling, apparently even at the ferromagnetic quantum critical point. The important role of the protocol chosen to sample hidden-spin configurations, in particular at the critical point, is analyzed. We discuss the implications of these findings for what concerns the problem of simulating adiabatic quantum optimization using stochastic algorithms on classical computers.

We study the robustness of the quantization of the Hall conductivity in the Harper-Hofstadter model towards the details of the protocol with which a longitudinal uniform driving force Fx(t) is turned on. In the vector potential gauge, through Peierls substitution, this involves the switching on of complex time-dependent hopping amplitudes e-iAx(t) in the ◯ direction such that ∂tAx(t)=Fx(t). The switching on can be sudden, Fx(t)=θ(t)F, where F is the steady driving force, or more generally smooth Fx(t)=f(t/t0)F, where f(t/t0) is such that f(0)=0 and f(1)=1. We investigate how the time-averaged (steady-state) particle current density jy in the ŷ direction deviates from the quantized value jyh/F=n due to the finite value of F and the details of the switching-on protocol. Exploiting the time periodicity of the Hamiltonian Ĥ(t), we use Floquet techniques to study this problem. In this picture the (Kubo) linear response F→0 regime corresponds to the adiabatic limit for Ĥ(t). In the case of a sudden quench jyh/F shows F2 corrections to the perfectly quantized limit. When the switching on is smooth, the result depends on the switch-on time t0: For a fixed t0 we observe a crossover force F∗ between a quadratic regime for FF∗. The crossover F∗ decreases as t0 increases, eventually recovering the topological robustness. These effects are in principle amenable to experimental tests in optical lattice cold atomic systems with synthetic gauge fields.

We study the effect of a thermal environment on the quantum annealing dynamics of a transverse-field Ising chain. The environment is modeled as a single Ohmic bath of quantum harmonic oscillators weakly interacting with the total transverse magnetization of the chain in a translationally invariant manner. We show that the density of defects generated at the end of the annealing process displays a minimum as a function of the annealing time, the so-called optimal working point, only in rather special regions of the bath temperature and coupling strength plane. We discuss the relevance of our results for current and future experimental implementations with quantum annealing hardware.

Figure 6 has a typographical error in the key: The labels corresponding to RWA and without the RWA are wrongly switched. The caption as well as all the rest of the paper are nevertheless correct and unaffected by this error. We show the exact figure here below with its original caption from the paper. (Figure Presented).

The quantum motion of nuclei, generally ignored in the physics of sliding friction, can affect in an important manner the frictional dissipation of a light particle forced to slide in an optical lattice. The density matrix-calculated evolution of the quantum version of the basic Prandtl–Tomlinson model, describing the dragging by an external force of a point particle in a periodic potential, shows that purely classical friction predictions can be very wrong. The strongest quantum effect occurs not for weak but for strong periodic potentials, where barriers are high but energy levels in each well are discrete, and resonant Rabi or Landau–Zener tunneling to states in the nearest well can preempt classical stick–slip with nonnegligible efficiency, depending on the forcing speed. The resulting permeation of otherwise unsurmountable barriers is predicted to cause quantum lubricity, a phenomenon which we expect should be observable in the recently implemented sliding cold ion experiments.

Recent model simulations discovered unexpected nonmonotonic features in the wear-free dry phononic friction as a function of the sliding speed. Here we demonstrate that a rather straightforward application of linear-response theory, appropriate in a regime of weak slider-substrate interaction, predicts frictional one-phonon singularities which imply a nontrivial dependence of the dynamical friction force on the slider speed and/or coupling to the substrate. The explicit formula which we derive reproduces very accurately the classical atomistic simulations when available. By modifying the slider-substrate interaction the analytical understanding obtained provides a practical means to tailor and control the speed dependence of friction with substantial freedom.

We study Thouless pumping out of the adiabatic limit. Our findings show that despite its topological nature, this phenomenon is not generically robust to nonadiabatic effects. Indeed, we find that the Floquet diagonal ensemble value of the pumped charge shows a deviation from the topologically quantized limit which is quadratic in the driving frequency for a sudden switch on of the driving. This is reflected also in the charge pumped in a single period, which shows a nonanalytic behavior on top of an overall quadratic decrease. Exponentially small corrections are recovered only with a careful tailoring of the driving protocol. We also discuss thermal effects and the experimental feasibility of observing such a deviation.

In simple ferromagnetic quantum Ising models characterized by an effective double-well energy landscape the characteristic tunneling time of path-integral Monte Carlo (PIMC) simulations has been shown to scale as the incoherent quantum-tunneling time, i.e., as 1/Δ2, where Δ is the tunneling gap. Since incoherent quantum tunneling is employed by quantum annealers (QAs) to solve optimization problems, this result suggests that there is no quantum advantage in using QAs with respect to quantum Monte Carlo (QMC) simulations. A counterexample is the recently introduced shamrock model (Andriyash and Amin, arXiv:1703.09277), where topological obstructions cause an exponential slowdown of the PIMC tunneling dynamics with respect to incoherent quantum tunneling, leaving open the possibility for potential quantum speedup, even for stoquastic models. In this work we investigate the tunneling time of projective QMC simulations based on the diffusion Monte Carlo (DMC) algorithm without guiding functions, showing that it scales as 1/Δ, i.e., even more favorably than the incoherent quantum-tunneling time, both in a simple ferromagnetic system and in the more challenging shamrock model. However, a careful comparison between the DMC ground-state energies and the exact solution available for the transverse-field Ising chain indicates an exponential scaling of the computational cost required to keep a fixed relative error as the system size increases.

Revealing phase transitions of solids through mechanical anomalies in the friction of nanotips sliding on their surfaces, a successful approach for continuous transitions, is still an unexplored tool for first-order ones. Owing to slow nucleation, first-order structural transformations occur with hysteresis, comprised between two spinodal temperatures where, on both sides of the thermodynamic transition, one or the other metastable free energy branches terminates. The spinodal transformation, a collective one-shot event without heat capacity anomaly, is easy to trigger by a weak external perturbation. Here we show that even the gossamer mechanical action of an AFM-tip can locally act as a trigger, narrowly preempting the spontaneous spinodal transformation, and making it observable as a nanofrictional anomaly. Confirming this expectation, the CCDW-NCCDW first-order transition of the important layer compound 1T-TaS2 is shown to provide a demonstration of this effect.

The realization of the Hofstadter model in a strongly anisotropic ladder geometry has now become possible in one-dimensional optical lattices with a synthetic dimension. In this work, we show how the Hofstadter Hamiltonian in such ladder configurations hosts a topological phase of matter which is radically different from its two-dimensional counterpart. This topological phase stems directly from the hybrid nature of the ladder geometry and is protected by a properly defined inversion symmetry. We start our analysis by considering the paradigmatic case of a three-leg ladder which supports a topological phase exhibiting the typical features of topological states in one dimension: robust fermionic edge modes, a degenerate entanglement spectrum, and a nonzero Zak phase; then, we generalize our findings - addressable in the state-of-the-art cold-atom experiments - to ladders with a higher number of legs.

We introduce and study the adiabatic dynamics of free-fermion models subject to a local Lindblad bath and in the presence of a time-dependent Hamiltonian. The merit of these models is that they can be solved exactly, and will help us to study the interplay between nonadiabatic transitions and dissipation in many-body quantum systems. After the adiabatic evolution, we evaluate the excess energy (the average value of the Hamiltonian) as a measure of the deviation from reaching the final target ground state. We compute the excess energy in a variety of different situations, where the nature of the bath and the Hamiltonian is modified. We find robust evidence of the fact that an optimal working time for the quantum annealing protocol emerges as a result of the competition between the nonadiabatic effects and the dissipative processes. We compare these results with the matrix-product-operator simulations of an Ising system and show that the phenomenology we found also applies for this more realistic case.

We study the approach to the adiabatic limit in periodically driven systems. Specifically focusing on a spin-1/2 in a magnetic field we find that, when the parameters of the Hamiltonian lead to a quasi-degeneracy in the Floquet spectrum, the evolution is not adiabatic even if the frequency of the field is much smaller than the spectral gap of the Hamiltonian. We argue that this is a general phenomenon of periodically driven systems. Although an explanation based on a perturbation theory in cannot be given, because of the singularity of the zero frequency limit, we are able to describe this phenomenon by means of a mapping to an extended Hilbert space, in terms of resonances of an effective two-band Wannier-Stark ladder. Remarkably, the phenomenon survives in presence of dissipation towards an environment and can be therefore easily experimentally observed.

We compare the performance of quantum annealing (QA, through Schrödinger dynamics) and simulated annealing (SA, through a classical master equation) on the p-spin infinite range ferromagnetic Ising model, by slowly driving the system across its equilibrium, quantum or classical, phase transition. When the phase transition is second order (p=2, the familiar two-spin Ising interaction) SA shows a remarkable exponential speed-up over QA. For a first-order phase transition (p≥3, i.e., with multispin Ising interactions), in contrast, the classical annealing dynamics appears to remain stuck in the disordered phase, while we have clear evidence that QA shows a residual energy which decreases towards zero when the total annealing time τ increases, albeit in a rather slow (logarithmic) fashion. This is one of the rare examples where a limited quantum speedup, a speedup by QA over SA, has been shown to exist by direct solutions of the Schrödinger and master equations in combination with a nonequilibrium Landau-Zener analysis. We also analyze the imaginary-time QA dynamics of the model, finding a 1/τ2 behavior for all finite values of p, as predicted by the adiabatic theorem of quantum mechanics. The Grover-search limit p(odd)=∞ is also discussed.

We revisit here the issue of thermally assisted Quantum Annealing by a detailed study of the dissipative Landau-Zener problem in the presence of a Caldeira-Leggett bath of harmonic oscillators, using both a weak-coupling quantum master equation and a quasiadiabatic path-integral approach. Building on the known zero-temperature exact results [Wubs, Phys. Rev. Lett. 97, 200404 (2006)PRLTAO0031-900710.1103/PhysRevLett.97.200404], we show that a finite temperature bath can have a beneficial effect on the ground-state probability only if it couples also to a spin direction that is transverse with respect to the driving field, while no improvement is obtained for the more commonly studied purely longitudinal coupling. In particular, we also highlight that, for a transverse coupling, raising the bath temperature further improves the ground-state probability in the fast-driving regime. We discuss the relevance of these findings for the current quantum-annealing flux qubit chips.

In this work we study the entanglement entropy of a uniform quantum Ising chain in transverse field undergoing a periodic driving of period τ. By means of Floquet theory we show that, for any subchain, the entanglement entropy tends asymptotically to a value τ-periodic in time. We provide a semi-analytical formula for the leading term of this asymptotic regime: It is constant in time and obeys a volume law. The entropy in the asymptotic regime is always smaller than the thermal one: because of integrability the system locally relaxes to a generalized Gibbs ensemble (GGE) density matrix. The leading term of the asymptotic entanglement entropy is completely determined by this GGE density matrix. Remarkably, the asymptotic entropy shows marked features in correspondence to some non-equilibrium quantum phase transitions undergone by a Floquet state analog of the ground state.

We show clear evidence of a quadratic speedup of a quantum annealing (QA) Schrödinger dynamics over a Glauber master equation simulated annealing (SA) for a random Ising model in one dimension, via an equal-footing exact deterministic dynamics of the Jordan-Wigner fermionized problems. This is remarkable, in view of the arguments of H. G. Katzgraber et al. [Phys. Rev. X 4, 021008 (2014)2160-330810.1103/PhysRevX.4.021008], since SA does not encounter any phase transition, while QA does. We also find a second remarkable result: that a "quantum-inspired" imaginary-time Schrödinger QA provides a further exponential speedup, i.e., an asymptotic residual error decreasing as a power law τ-μ of the annealing time τ.