Quantum Annealing for Heisenberg Spin Chains

Similar documents
Quantum Theory Assignment 3

Hydrogen (atoms, molecules) in external fields. Static electric and magnetic fields Oscyllating electromagnetic fields

The Born-Oppenheimer approximation

Math 2784 (or 2794W) University of Connecticut

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

Lecture 25 (Dec. 6, 2017)

Chapter 7: The z-transform. Chih-Wei Liu

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle

Quantum Simulation: Solving Schrödinger Equation on a Quantum Computer

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

PHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions

Assignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.

The Heisenberg versus the Schrödinger picture in quantum field theory. Dan Solomon Rauland-Borg Corporation 3450 W. Oakton Skokie, IL USA

C/CS/Phys C191 Deutsch and Deutsch-Josza algorithms 10/20/07 Fall 2007 Lecture 17

Chapter 4. Fourier Series

Probability, Expectation Value and Uncertainty

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Math 312 Lecture Notes One Dimensional Maps

Matsubara-Green s Functions

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT

Chapter 2 The Monte Carlo Method

The time evolution of the state of a quantum system is described by the time-dependent Schrödinger equation (TDSE): ( ) ( ) 2m "2 + V ( r,t) (1.

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled

Areas and Distances. We can easily find areas of certain geometric figures using well-known formulas:

INTEGRATION BY PARTS (TABLE METHOD)

Chapter 9 - CD companion 1. A Generic Implementation; The Common-Merge Amplifier. 1 τ is. ω ch. τ io

Polynomial Functions and Their Graphs

Section 13.3 Area and the Definite Integral

1 Adiabatic and diabatic representations

Physics 232 Gauge invariance of the magnetic susceptibilty

Kinetics of Complex Reactions

3/21/2017. Commuting and Non-commuting Operators Chapter 17. A a

Problem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

Introduction to Signals and Systems, Part V: Lecture Summary

Time-Domain Representations of LTI Systems

Exercises and Problems

Warped, Chirp Z-Transform: Radar Signal Processing

Principle Of Superposition

5.76 Lecture #33 5/08/91 Page 1 of 10 pages. Lecture #33: Vibronic Coupling

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

Introduction to Computational Molecular Biology. Gibbs Sampling

Frequency Response of FIR Filters

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.

Fastest mixing Markov chain on a path

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Lecture III-2: Light propagation in nonmagnetic

Chimica Inorganica 3

6.003 Homework #3 Solutions

For example suppose we divide the interval [0,2] into 5 equal subintervals of length

MDIV. Multiple divisor functions

A Block Cipher Using Linear Congruences

Some properties of Boubaker polynomials and applications

Exponents. Learning Objectives. Pre-Activity

FIR Filter Design: Part II

The Z-Transform. (t-t 0 ) Figure 1: Simplified graph of an impulse function. For an impulse, it can be shown that (1)

SECTION 2 Electrostatics

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

PHY4905: Nearly-Free Electron Model (NFE)

The Random Walk For Dummies

Name Solutions to Test 2 October 14, 2015

Solution of Quantum Anharmonic Oscillator with Quartic Perturbation

True Nature of Potential Energy of a Hydrogen Atom

Things you should know when you leave Discussion today for one-electron atoms:

PRELIMINARY EXAMINATION Department of Physics University of Florida Part A, January, 2016, 09:00 12:00. Instructions

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =

For example suppose we divide the interval [0,2] into 5 equal subintervals of length

THE KALMAN FILTER RAUL ROJAS

Lanczos-Haydock Recursion

Formation of A Supergain Array and Its Application in Radar

Lecture 2: Monte Carlo Simulation

x a x a Lecture 2 Series (See Chapter 1 in Boas)

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

6.3 Testing Series With Positive Terms

COMPUTING SUMS AND THE AVERAGE VALUE OF THE DIVISOR FUNCTION (x 1) + x = n = n.

PC5215 Numerical Recipes with Applications - Review Problems

The Wave Function and Quantum Reality

1. Szabo & Ostlund: 2.1, 2.2, 2.4, 2.5, 2.7. These problems are fairly straightforward and I will not discuss them here.

arxiv: v1 [cond-mat.stat-mech] 25 Jan 2017

COURSE INTRODUCTION: WHAT HAPPENS TO A QUANTUM PARTICLE ON A PENDULUM π 2 SECONDS AFTER IT IS TOSSED IN?

arxiv: v1 [math-ph] 5 Jul 2017

MIT Department of Chemistry 5.74, Spring 2005: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff

Z ß cos x + si x R du We start with the substitutio u = si(x), so du = cos(x). The itegral becomes but +u we should chage the limits to go with the ew

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

Filter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and

mx bx kx F t. dt IR I LI V t, Q LQ RQ V t,

Chapter 9: Numerical Differentiation

Alternating Series. 1 n 0 2 n n THEOREM 9.14 Alternating Series Test Let a n > 0. The alternating series. 1 n a n.

TIME-CORRELATION FUNCTIONS

L = n i, i=1. dp p n 1

Quantum Mechanics I. 21 April, x=0. , α = A + B = C. ik 1 A ik 1 B = αc.

LECTURE 14. Non-linear transverse motion. Non-linear transverse motion

HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples:

The Minimum Distance Energy for Polygonal Unknots

Sequences I. Chapter Introduction

Chapter 2 Feedback Control Theory Continued

CALCULATION OF FIBONACCI VECTORS

a b c d e f g h Supplementary Information

Marcinkiwiecz-Zygmund Type Inequalities for all Arcs of the Circle

Diffusivity and Mobility Quantization. in Quantum Electrical Semi-Ballistic. Quasi-One-Dimensional Conductors

Transcription:

LA-UR # - Quatum Aealig for Heiseberg Spi Chais G.P. Berma, V.N. Gorshkov,, ad V.I.Tsifriovich Theoretical Divisio, Los Alamos Natioal Laboratory, Los Alamos, NM Istitute of Physics, Natioal Academy of Scieces, Kiev, Ukraie Physics Departmet, Polytechic Uiversity, Brookly, New York, Abstract We suggest usig the method of quatum aealig for computig the groud state of the Heiseberg spi chais. Our iitial Hamiltoia describes a spi system i a highly ouiform magetic field. The iitial Hamiltoia gradually trasforms ito the Heiseberg Hamiltoia, which describes the exchage iteractio ad the Zeema iteractio with the uiform magetic field. We demostrate quatum aealig for -, -, ad -spi systems. I particular, we cosider the alterate (ferro- ad atiferromagetic) exchage itegral. By chagig the sig of the exchage itegrals we switch from the frustrated to the correspodig o-frustrated spi system. We discuss the similarity ad differece betwee the frustrated ad o-frustrated spi systems.. Itroductio Recetly, the quatum aealig method has bee widely used to compute the groud state of Isig spi chais []. Quatum aealig assumes that the iitial Hamiltoia of the system has a simple groud state ad does ot commute with the Isig Hamiltoia. The iitial Hamiltoia smoothly (adiabatically) trasforms ito the Isig Hamiltoia. I this case a simple iitial groud state smoothly trasforms ito the complicated groud state of the actual Isig Hamiltoia. Typically, the iitial Hamiltoia describes the iteractio of the spi system with the trasverse magetic field. The Isig spi chai has umerous applicatios to mathematically itractable problems. From the other side the Heiseberg spi chais may be used to create atomic-scale magetic structures, i particular, for the data storage devices. Heiseberg spi chais are actively studied ow for applicatios i aotechology []. Cosequetly, we cosider the applicatio of the quatum aealig to the Heiseberg spi chai. The uiform trasverse magetic field caot be used ow because the isotropic Heiseberg Hamiltoia iteracts with ay compoet of the total spi. Thus, we have to cosider iteractio with ouiform magetic field, which does ot commute with the Heiseberg Hamiltoia.

. Two-spi Heiseberg system I order to uderstad the properties of the quatum aealig computatio we start with a trivial two-spi Heiseberg system. We cosider the followig timedepedet Hamiltoia of the system: ˆ () ( z z t ) ( x x H t = JSS + γ B S + S + γ B S S )( t / τ ), () τ where < t < τ, J is the exchage iteractio costat, γ is the absolute value of the electro gyromagetic ratio, BB is the uiform magetic field which poits i the positive z - directio, B is the magitude of the o-uiform magetic field which poits i the positive x - directio at the locatio of the first spi ad poits i the egative x - directio at the locatio of the secod spi. A o-uiform magetic field is itroduced i order to fid the groud state of the Heiseberg system usig the quatum aealig approach. At time t = the Hamiltoia () describes two idepedet spis iteractig with the o-uiform magetic field. Thus, the eigestate of the Hamiltoia () at t = represets the first spi poitig i the egative x -directio ad the secod spi poitig i the positive x -directio. We choose this state as the iitial state for the quatum aealig. At t = τ the Hamiltoia () describes the coupled Heiseberg spis iteractig with the uiform magetic field BB, which poits i the positive z - directio. If τ is large eough oe ca expect that the groud state of the iitial Hamiltoia Ĥ(t = ) will trasfer smoothly to the groud state of the fial Hamiltoia Ĥ(t = τ ). Below we describe our computer simulatios for the twospi system. We use BB as the uit of the magetic field, γb B as the uit of frequecy, γħbb as the uit of eergy, ad /(γħb B ) as the uit of time. I these uits we have chose the followig parameters of our Hamiltoia: J =, B =, ad τ =. () Positive value of J correspods to the simplest ferromagetic groud state. We take z the eigestates of the operators, S k, as the basis states ad represet the wave fuctio as a superpositio of all basis states with the amplitudes C. The iitial values of the amplitudes are C = - ½, C = - ½, C = ½, C = ½. () Here ad below we assume that for a sigle spi state represets the dow directio, ad state represets the up directio. For two or more spis we use the same rules to trasfer from the biary to the decimal otatio. Thus, for the two

spis state represets, state represets, state represets, ad state represets. Note that we cout spis from the left to the right. For example, the state idicates that the first spi poits dow while the secod spi poits up. Fig. a shows evolutio of the probabilities p = C. Oe ca see that the probability p gradually approaches to uity, while all the other probabilities approach to zero. Cosequetly, the spi system approaches to the ferromagetic groud state with spi poitig i the egative z - directio. Fig. a shows probabilities p for the atiferomagetic exchage iteractio J =. I this case probabilities p ad p approach the value. while the two other probabilities vaish. This reflects evolutio to the atiferromagetic state S =, which is a superpositio of the two states ad with equal probabilities ad phase differece π. (Phase is ot show i Figs. a ad a.) I Figs. b ad b, we show also the chage of the eergy levels for J = ad J =, correspodigly. These graphs do ot have a direct coectio to the quatum aealig. They are preseted to better uderstad the spi systems uder cosideratio.. The Heiseberg system of three ad more spis For three or more spis we use a Hamiltoia similar to () ˆ () k t ( / ) ( ) k+ Ht = JkmSS k m + γ B Sz + γ B tτ S x k, () km, k τ k At time t = the spis poit i the directio of their local magetic field (i.e. i the positive x - directio for eve k, ad the egative x - directio for odd k ). We assume a closed chai of spis. If the product of the exchage itegrals is a odd umber the the spi system is frustrated. Fig. a demostrates quatum aealig to the ferromagetic groud state for the - spi system with J km =. Fig. b shows the evolutio of the eergy levels for the same system. Note that the probability p does ot icrease mootoically as it did for the -spi system. Fig. a shows quatum aealig for the frustrated -spi system: J = J =, J = -, ad Fig. b shows the chage of the eergy levels. I this case, the groud state is a superpositio of the states = ad = with equal probabilities ad phase shift π. Oe ca see that the probabilities p ad p mootoically icrease durig the process of quatum aealig. Surprisigly, quatum aealig for the frustrated system appears simpler tha for the ferromagetic system.

Next, we cosider a more complicated system of spis. The iitial eergy levels for t = are show i Fig.. All eergy levels except for the lowest ad the highest oes are degeerate. Fig. shows the fial eergy levels at t = τ for the case of the ferromagetic exchage iteractio, J km =. I the scale of the picture we observe a almost cotiuous eergy spectrum except for the highest levels. The groud state of the system is the ferromagetic state with the total spi S = / ad the z- z compoet of the total spi S = /. The first excited state is the same z ferromagetic state with S = / ad S = /. Correspodigly, the spacig betwee the groud ad the first excited state is exactly oe uit, which is ot visible i the scale of Fig.. Quatum aealig drives the spi system to the ferromagetic groud state. Now, we cosider a system of spis with alterate exchage itegral. This spi system is show i Fig.. The five double lies i the Fig. correspod to positive values of the exchage itegral J =, ad the four solid lies correspod to egative values, J =. The fial eergy spectrum is show i Fig.. The two lowest eergy levels have a Zeema spacig of oe uit, the same as for the ferromagetic case. I the scale of Fig. these two levels are represeted by a sigle dot. These two levels are separated by a gap of approximately. uits from the followig quasi-cotiuous spectrum. The secod gap appears ear the highest levels. The groud state obtaied by the quatum aealig is a complicated superpositio of the basis states. (See Fig..) The maximum probability correspods to the basis state : p.. Fig. demostrates the structure of the state : the circles at locatios,,, ad show the spis, which poit up. Fially we discuss the same spi system with the opposite sig of the exchage itegrals J km. This spi system is frustrated because the umber of egative exchage itegrals is odd. The eergy spectrum of the frustrated spi chai is very similar to that of the o-frustrated oe. (See Fig..) At the same time, the structure of the groud state, which is obtaied usig quatum aealig, is much more complicated. Fig. shows the values of the amplitudes C. Oe ca see that the groud state of the frustrated spi system has four domiat states istead of a sigle oe for the o-frustrated system. These domiat states are show i Fig.. The correspodig values of the amplitudes are: C C =., C = -., C =., C = -.. () Thus, we have the two pairs of the domiat states. The two states i a pair have approximately equal probabilities ad a phase shift of π.

. Coclusio We have demostrated the successful applicatio of the quatum aealig method to Heiseberg spi chais. I order to use quatum aealig we apply a highly o-uiform magetic field. The correspodig Hamiltoia does ot commute with both the Zeema ad the exchage terms of the basic Hamiltoia. First, we demostrated fast quatum aealig for the simple two- ad three-spi systems. The we cosidered a more complicated -spi system. We have show that the frustrated spi system has almost the same eergy spectrum as the correspodig o-frustrated system but the structure of the groud state for the two systems is sigificatly differet. The o-frustrated system has a sigle domiat state while the correspodig frustrated system has the two pairs of the domiat states. The states i each pair have approximately equal probability ad a phase shift of π. This work was carried out uder the auspices of the Natioal Nuclear Security Admiistratio of the U.S. Departmet of Eergy at Los Alamos Natioal Laboratory uder Cotract No. DE-AC-NA. Refereces. Arab Das, Bikas Chakrabarti (Eds), Quatum Aealig ad Other Optimizatio Methods, Lecture Notes i Physics, Publisher: Spriger Berli / Heidelberg, V. /.. Mark Wilso, Physics Today,, ().

. p a. = E (t) b.........., - -...... Fig.. a-probabilities of the basis states, which cotribute to the groud state of the -spi system with the ferromagetic exchage iteractio, as a fuctio of time i the process of the quatum aealig; b- the eergy levels of the spi system as a fuctio of time.. p a, E (t) b.... = -....... -...... Fig.. The same as i Fig. for the -spi system with the atiferromagetic exchage iteractio.

. p a.. = E (t) b.. - -....... -...... Fig.. The same as i Fig. - for the three spi system with the ferromagetic exchage iteractio.. p a........... = - - E (t) b -...... Fig.. The same as i Figs. - for the frustrated -spi system.

- - - - E (t=) Fig.. The eergy levels of the -spi system at the begiig of the quatum aealig. - - E (t =τ) Fig.. The eergy of the statioary states of the - spi system with the ferromagetic exchage iteractio E (t =τ) Fig.. The -spi chai with the alterate ferro- atiferromagetic couplig (solid lie - J =, dоuble lie - J = ). - - - - Fig.. The eergy of the statioary states of the ofrustrated -spi system with the alterate exchage iteractio.

. p... Fig.. Probabilities of the basis states, which cotribute to the groud state of the o-frustrated -spi system with the alterate exchage iteractio. Fig.. The structure of the basis state =, which makes the mai cotributio to the groud state. Circles idicate spis poitig up. E (t =τ). C (t =τ)... - - - - Fig.. The eergy of the statioary states of the frustrated -spi system with the alterate exchage iteractio. -. -. -. Fig.. Amplitudes of the basis states, which cotribute to the groud state of the frustrated -spi system with the alterate exchage iteractio.

= = = = Fig.. The structure of the basis states, which make the mai cotributio to the groud state. Circles idicate spis poitig up.