Quantum Annealing for Heisenberg Spin Chains

Size: px
Start display at page:

Download "Quantum Annealing for Heisenberg Spin Chains"

Transcription

1 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.

2 . 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

3 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.

4 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 π.

5 . 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,, ().

6 . 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.

7 . 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.

8 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.

9 . 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.

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

Quantum Theory Assignment 3

Quantum Theory Assignment 3 Quatum Theory Assigmet 3 Assigmet 3.1 1. Cosider a spi-1/ system i a magetic field i the z-directio. The Hamiltoia is give by: ) eb H = S z = ωs z. mc a) Fid the Heiseberg operators S x t), S y t), ad

More information

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

Hydrogen (atoms, molecules) in external fields. Static electric and magnetic fields Oscyllating electromagnetic fields Hydroge (atoms, molecules) i exteral fields Static electric ad magetic fields Oscyllatig electromagetic fields Everythig said up to ow has to be modified more or less strogly if we cosider atoms (ad ios)

More information

The Born-Oppenheimer approximation

The Born-Oppenheimer approximation The Bor-Oppeheimer approximatio 1 Re-writig the Schrödiger equatio We will begi from the full time-idepedet Schrödiger equatio for the eigestates of a molecular system: [ P 2 + ( Pm 2 + e2 1 1 2m 2m m

More information

Math 2784 (or 2794W) University of Connecticut

Math 2784 (or 2794W) University of Connecticut ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really

More information

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka) 7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.

More information

Lecture 25 (Dec. 6, 2017)

Lecture 25 (Dec. 6, 2017) Lecture 5 8.31 Quatum Theory I, Fall 017 106 Lecture 5 (Dec. 6, 017) 5.1 Degeerate Perturbatio Theory Previously, whe discussig perturbatio theory, we restricted ourselves to the case where the uperturbed

More information

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

Chapter 7: The z-transform. Chih-Wei Liu Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability

More information

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

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle Similarity betwee quatum mechaics ad thermodyamics: Etropy, temperature, ad Carot cycle Sumiyoshi Abe 1,,3 ad Shiji Okuyama 1 1 Departmet of Physical Egieerig, Mie Uiversity, Mie 514-8507, Japa Istitut

More information

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

Quantum Simulation: Solving Schrödinger Equation on a Quantum Computer Purdue Uiversity Purdue e-pubs Birc Poster Sessios Birc Naotechology Ceter 4-14-008 Quatum Simulatio: Solvig Schrödiger Equatio o a Quatum Computer Hefeg Wag Purdue Uiversity, wag10@purdue.edu Sabre Kais

More information

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

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001. Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;

More information

PHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions

PHYC - 505: Statistical Mechanics Homework Assignment 4 Solutions PHYC - 55: Statistical Mechaics Homewor Assigmet 4 Solutios Due February 5, 14 1. Cosider a ifiite classical chai of idetical masses coupled by earest eighbor sprigs with idetical sprig costats. a Write

More information

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 } χ.

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 } χ. PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio

More information

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

The Heisenberg versus the Schrödinger picture in quantum field theory. Dan Solomon Rauland-Borg Corporation 3450 W. Oakton Skokie, IL USA 1 The Heiseberg versus the chrödiger picture i quatum field theory by Da olomo Raulad-Borg Corporatio 345 W. Oakto kokie, IL 677 UA Phoe: 847-324-8337 Email: da.solomo@raulad.com PAC 11.1-z March 15, 24

More information

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

C/CS/Phys C191 Deutsch and Deutsch-Josza algorithms 10/20/07 Fall 2007 Lecture 17 C/CS/Phs C9 Deutsch ad Deutsch-Josza algorithms 0/0/07 Fall 007 Lecture 7 Readigs Beeti et al., Ch. 3.9-3.9. Stolze ad Suter, Quatum Computig, Ch. 8. - 8..5) Nielse ad Chuag, Quatum Computatio ad Quatum

More information

Chapter 4. Fourier Series

Chapter 4. Fourier Series Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,

More information

Probability, Expectation Value and Uncertainty

Probability, Expectation Value and Uncertainty Chapter 1 Probability, Expectatio Value ad Ucertaity We have see that the physically observable properties of a quatum system are represeted by Hermitea operators (also referred to as observables ) such

More information

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

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 The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

Math 312 Lecture Notes One Dimensional Maps

Math 312 Lecture Notes One Dimensional Maps Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,

More information

Matsubara-Green s Functions

Matsubara-Green s Functions Matsubara-Gree s Fuctios Time Orderig : Cosider the followig operator If H = H the we ca trivially factorise this as, E(s = e s(h+ E(s = e sh e s I geeral this is ot true. However for practical applicatio

More information

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

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT TR/46 OCTOBER 974 THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION by A. TALBOT .. Itroductio. A problem i approximatio theory o which I have recetly worked [] required for its solutio a proof that the

More information

Chapter 2 The Monte Carlo Method

Chapter 2 The Monte Carlo Method Chapter 2 The Mote Carlo Method The Mote Carlo Method stads for a broad class of computatioal algorithms that rely o radom sampligs. It is ofte used i physical ad mathematical problems ad is most useful

More information

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 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. Adrei Tokmakoff, MIT Departmet of Chemistry, 2/13/2007 1-1 574 TIME-DEPENDENT QUANTUM MECHANICS 1 INTRODUCTION 11 Time-evolutio for time-idepedet Hamiltoias The time evolutio of the state of a quatum system

More information

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

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled 1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how

More information

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

Areas and Distances. We can easily find areas of certain geometric figures using well-known formulas: Areas ad Distaces We ca easily fid areas of certai geometric figures usig well-kow formulas: However, it is t easy to fid the area of a regio with curved sides: METHOD: To evaluate the area of the regio

More information

INTEGRATION BY PARTS (TABLE METHOD)

INTEGRATION BY PARTS (TABLE METHOD) INTEGRATION BY PARTS (TABLE METHOD) Suppose you wat to evaluate cos d usig itegratio by parts. Usig the u dv otatio, we get So, u dv d cos du d v si cos d si si d or si si d We see that it is ecessary

More information

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

Chapter 9 - CD companion 1. A Generic Implementation; The Common-Merge Amplifier. 1 τ is. ω ch. τ io Chapter 9 - CD compaio CHAPTER NINE CD-9.2 CD-9.2. Stages With Voltage ad Curret Gai A Geeric Implemetatio; The Commo-Merge Amplifier The advaced method preseted i the text for approximatig cutoff frequecies

More information

Polynomial Functions and Their Graphs

Polynomial Functions and Their Graphs Polyomial Fuctios ad Their Graphs I this sectio we begi the study of fuctios defied by polyomial expressios. Polyomial ad ratioal fuctios are the most commo fuctios used to model data, ad are used extesively

More information

Section 13.3 Area and the Definite Integral

Section 13.3 Area and the Definite Integral Sectio 3.3 Area ad the Defiite Itegral We ca easily fid areas of certai geometric figures usig well-kow formulas: However, it is t easy to fid the area of a regio with curved sides: METHOD: To evaluate

More information

1 Adiabatic and diabatic representations

1 Adiabatic and diabatic representations 1 Adiabatic ad diabatic represetatios 1.1 Bor-Oppeheimer approximatio The time-idepedet Schrödiger equatio for both electroic ad uclear degrees of freedom is Ĥ Ψ(r, R) = E Ψ(r, R), (1) where the full molecular

More information

Physics 232 Gauge invariance of the magnetic susceptibilty

Physics 232 Gauge invariance of the magnetic susceptibilty Physics 232 Gauge ivariace of the magetic susceptibilty Peter Youg (Dated: Jauary 16, 2006) I. INTRODUCTION We have see i class that the followig additioal terms appear i the Hamiltoia o addig a magetic

More information

Kinetics of Complex Reactions

Kinetics of Complex Reactions Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet

More information

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

3/21/2017. Commuting and Non-commuting Operators Chapter 17. A a Commutig ad No-commutig Operators Chapter 17 Postulate 3. I ay measuremet of the observable associated with a operator A the oly values that will ever be observed are the eige values, a, which satisfy

More information

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 =

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 = Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,

More information

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A)

REGRESSION (Physics 1210 Notes, Partial Modified Appendix A) REGRESSION (Physics 0 Notes, Partial Modified Appedix A) HOW TO PERFORM A LINEAR REGRESSION Cosider the followig data poits ad their graph (Table I ad Figure ): X Y 0 3 5 3 7 4 9 5 Table : Example Data

More information

Introduction to Signals and Systems, Part V: Lecture Summary

Introduction to Signals and Systems, Part V: Lecture Summary EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive

More information

Time-Domain Representations of LTI Systems

Time-Domain Representations of LTI Systems 2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable

More information

Exercises and Problems

Exercises and Problems HW Chapter 4: Oe-Dimesioal Quatum Mechaics Coceptual Questios 4.. Five. 4.4.. is idepedet of. a b c mu ( E). a b m( ev 5 ev) c m(6 ev ev) Exercises ad Problems 4.. Model: Model the electro as a particle

More information

Warped, Chirp Z-Transform: Radar Signal Processing

Warped, Chirp Z-Transform: Radar Signal Processing arped, Chirp Z-Trasform: Radar Sigal Processig by Garimella Ramamurthy Report o: IIIT/TR// Cetre for Commuicatios Iteratioal Istitute of Iformatio Techology Hyderabad - 5 3, IDIA Jauary ARPED, CHIRP Z

More information

Principle Of Superposition

Principle Of Superposition ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give

More information

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

5.76 Lecture #33 5/08/91 Page 1 of 10 pages. Lecture #33: Vibronic Coupling 5.76 Lecture #33 5/8/9 Page of pages Lecture #33: Vibroic Couplig Last time: H CO A A X A Electroically forbidde if A -state is plaar vibroically allowed to alterate v if A -state is plaar iertial defect

More information

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

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as

More information

Introduction to Computational Molecular Biology. Gibbs Sampling

Introduction to Computational Molecular Biology. Gibbs Sampling 18.417 Itroductio to Computatioal Molecular Biology Lecture 19: November 16, 2004 Scribe: Tushara C. Karuarata Lecturer: Ross Lippert Editor: Tushara C. Karuarata Gibbs Samplig Itroductio Let s first recall

More information

Frequency Response of FIR Filters

Frequency Response of FIR Filters EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we itroduce the idea of the frequecy respose of LTI systems, ad focus specifically o the frequecy respose of FIR filters.. Steady-state

More information

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

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer. 6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio

More information

Fastest mixing Markov chain on a path

Fastest mixing Markov chain on a path Fastest mixig Markov chai o a path Stephe Boyd Persi Diacois Ju Su Li Xiao Revised July 2004 Abstract We ider the problem of assigig trasitio probabilities to the edges of a path, so the resultig Markov

More information

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

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,

More information

Lecture III-2: Light propagation in nonmagnetic

Lecture III-2: Light propagation in nonmagnetic A. La Rosa Lecture Notes ALIED OTIC Lecture III2: Light propagatio i omagetic materials 2.1 urface ( ), volume ( ), ad curret ( j ) desities produced by arizatio charges The objective i this sectio is

More information

Chimica Inorganica 3

Chimica Inorganica 3 himica Iorgaica Irreducible Represetatios ad haracter Tables Rather tha usig geometrical operatios, it is ofte much more coveiet to employ a ew set of group elemets which are matrices ad to make the rule

More information

6.003 Homework #3 Solutions

6.003 Homework #3 Solutions 6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the

More information

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

For example suppose we divide the interval [0,2] into 5 equal subintervals of length Math 120c Calculus Sec 1: Estimatig with Fiite Sums I Area A Cosider the problem of fidig the area uder the curve o the fuctio y!x 2 + over the domai [0,2] We ca approximate this area by usig a familiar

More information

MDIV. Multiple divisor functions

MDIV. Multiple divisor functions MDIV. Multiple divisor fuctios The fuctios τ k For k, defie τ k ( to be the umber of (ordered factorisatios of ito k factors, i other words, the umber of ordered k-tuples (j, j 2,..., j k with j j 2...

More information

A Block Cipher Using Linear Congruences

A Block Cipher Using Linear Congruences Joural of Computer Sciece 3 (7): 556-560, 2007 ISSN 1549-3636 2007 Sciece Publicatios A Block Cipher Usig Liear Cogrueces 1 V.U.K. Sastry ad 2 V. Jaaki 1 Academic Affairs, Sreeidhi Istitute of Sciece &

More information

Some properties of Boubaker polynomials and applications

Some properties of Boubaker polynomials and applications Some properties of Boubaker polyomials ad applicatios Gradimir V. Milovaović ad Duša Joksimović Citatio: AIP Cof. Proc. 179, 1050 (2012); doi: 10.1063/1.756326 View olie: http://dx.doi.org/10.1063/1.756326

More information

Exponents. Learning Objectives. Pre-Activity

Exponents. Learning Objectives. Pre-Activity Sectio. Pre-Activity Preparatio Epoets A Chai Letter Chai letters are geerated every day. If you sed a chai letter to three frieds ad they each sed it o to three frieds, who each sed it o to three frieds,

More information

FIR Filter Design: Part II

FIR Filter Design: Part II EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we cosider how we might go about desigig FIR filters with arbitrary frequecy resposes, through compositio of multiple sigle-peak

More information

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

The Z-Transform. (t-t 0 ) Figure 1: Simplified graph of an impulse function. For an impulse, it can be shown that (1) The Z-Trasform Sampled Data The geeralied fuctio (t) (also kow as the impulse fuctio) is useful i the defiitio ad aalysis of sampled-data sigals. Figure below shows a simplified graph of a impulse. (t-t

More information

SECTION 2 Electrostatics

SECTION 2 Electrostatics SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial

More information

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

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet

More information

PHY4905: Nearly-Free Electron Model (NFE)

PHY4905: Nearly-Free Electron Model (NFE) PHY4905: Nearly-Free Electro Model (NFE) D. L. Maslov Departmet of Physics, Uiversity of Florida (Dated: Jauary 12, 2011) 1 I. REMINDER: QUANTUM MECHANICAL PERTURBATION THEORY A. No-degeerate eigestates

More information

The Random Walk For Dummies

The Random Walk For Dummies The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oe-dimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli

More information

Name Solutions to Test 2 October 14, 2015

Name Solutions to Test 2 October 14, 2015 Name Solutios to Test October 4, 05 This test cosists of three parts. Please ote that i parts II ad III, you ca skip oe questio of those offered. The equatios below may be helpful with some problems. Costats

More information

Solution of Quantum Anharmonic Oscillator with Quartic Perturbation

Solution of Quantum Anharmonic Oscillator with Quartic Perturbation ISS -79X (Paper) ISS 5-0638 (Olie) Vol.7, 0 Abstract Solutio of Quatum Aharmoic Oscillator with Quartic Perturbatio Adelaku A.O. Departmet of Physics, Wesley Uiversity of Sciece ad Techology, Odo, Odo

More information

True Nature of Potential Energy of a Hydrogen Atom

True Nature of Potential Energy of a Hydrogen Atom True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial

More information

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

Things you should know when you leave Discussion today for one-electron atoms: E = -R Thigs ou should kow whe ou leave Discussio toda for oe-electro atoms: = -.79 0-8 J = -.6eV ΔEmatter=E-Em ; Ioizatio Eerg=E E(iitial) ΔΕlight=hνlight= IE +KE. Cosider the followig eerg levels of

More information

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

PRELIMINARY EXAMINATION Department of Physics University of Florida Part A, January, 2016, 09:00 12:00. Instructions Studet ID Number: PRELIMINRY EXMINTION Part, Jauary, 6, 9: : Istructios. You may use a calculator ad CRC Math tables or equivalet. No other tables or aids are allowed or required. You may NOT use programmable

More information

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

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) = AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,

More information

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

For example suppose we divide the interval [0,2] into 5 equal subintervals of length Math 1206 Calculus Sec 1: Estimatig with Fiite Sums Abbreviatios: wrt with respect to! for all! there exists! therefore Def defiitio Th m Theorem sol solutio! perpedicular iff or! if ad oly if pt poit

More information

THE KALMAN FILTER RAUL ROJAS

THE KALMAN FILTER RAUL ROJAS THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider

More information

Lanczos-Haydock Recursion

Lanczos-Haydock Recursion Laczos-Haydock Recursio Bor: Feb 893 i Székesehérvár, Hugary- Died: 5 Jue 974 i Budapest, Hugary Corelius Laczos From a abstract mathematical viewpoit, the method for puttig a symmetric matrix i three-diagoal

More information

Formation of A Supergain Array and Its Application in Radar

Formation of A Supergain Array and Its Application in Radar Formatio of A Supergai Array ad ts Applicatio i Radar Tra Cao Quye, Do Trug Kie ad Bach Gia Duog. Research Ceter for Electroic ad Telecommuicatios, College of Techology (Coltech, Vietam atioal Uiversity,

More information

Lecture 2: Monte Carlo Simulation

Lecture 2: Monte Carlo Simulation STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?

More information

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

x a x a Lecture 2 Series (See Chapter 1 in Boas) Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio

More information

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

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

COMPUTING SUMS AND THE AVERAGE VALUE OF THE DIVISOR FUNCTION (x 1) + x = n = n. COMPUTING SUMS AND THE AVERAGE VALUE OF THE DIVISOR FUNCTION Abstract. We itroduce a method for computig sums of the form f( where f( is ice. We apply this method to study the average value of d(, where

More information

PC5215 Numerical Recipes with Applications - Review Problems

PC5215 Numerical Recipes with Applications - Review Problems PC55 Numerical Recipes with Applicatios - Review Problems Give the IEEE 754 sigle precisio bit patter (biary or he format) of the followig umbers: 0 0 05 00 0 00 Note that it has 8 bits for the epoet,

More information

The Wave Function and Quantum Reality

The Wave Function and Quantum Reality The Wave Fuctio ad Quatum Reality Sha Gao Uit for History ad Philosophy of Sciece & Cetre for Time, SOPHI Uiversity of Sydey, Sydey, NSW 006, Australia Abstract. We ivestigate the meaig of the wave fuctio

More information

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.

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. Solutio set III.. Szabo & Ostlud:.,.,.,.5,.7. These problems are fairly straightforward ad I will ot discuss them here.. N! N! i= k= N! N! N! N! p p i j pi+ pj i j i j i= j= i= j= AA ˆˆ= ( ) Pˆ ( ) Pˆ

More information

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

arxiv: v1 [cond-mat.stat-mech] 25 Jan 2017 Iteractig Multi-particle Classical Szilard Egie P. S. Pal a,b ad A. M. Jayaavar a,b a Istitute of Physics, Sachivalaya Marg, Bhubaeswar-75005, Idia b Homi Bhabha atioal Istitute, Traiig School Comple,

More information

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

COURSE INTRODUCTION: WHAT HAPPENS TO A QUANTUM PARTICLE ON A PENDULUM π 2 SECONDS AFTER IT IS TOSSED IN? COURSE INTRODUCTION: WHAT HAPPENS TO A QUANTUM PARTICLE ON A PENDULUM π SECONDS AFTER IT IS TOSSED IN? DROR BAR-NATAN Follows a lecture give by the author i the trivial otios semiar i Harvard o April 9,

More information

arxiv: v1 [math-ph] 5 Jul 2017

arxiv: v1 [math-ph] 5 Jul 2017 O eigestates for some sl 2 related Hamiltoia arxiv:1707.01193v1 [math-ph] 5 Jul 2017 Fahad M. Alamrai Faculty of Educatio Sciece Techology & Mathematics, Uiversity of Caberra, Bruce ACT 2601, Australia.,

More information

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

MIT Department of Chemistry 5.74, Spring 2005: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff MIT Departmet of Chemistry 5.74, Sprig 5: Itroductory Quatum Mechaics II Istructor: Professor Adrei Tomaoff p. 97 ABSORPTION SPECTRA OF MOLECULAR AGGREGATES The absorptio spectra of periodic arrays of

More information

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

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 Problem ( poits) Evaluate the itegrals Z p x 9 x We ca draw a right triagle labeled this way x p x 9 From this we ca read off x = sec, so = sec ta, ad p x 9 = R ta. Puttig those pieces ito the itegralrwe

More information

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

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu

More information

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

Filter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and Filter bas Separately, the lowpass ad highpass filters are ot ivertible T removes the highest frequecy / ad removes the lowest frequecy Together these filters separate the sigal ito low-frequecy ad high-frequecy

More information

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

mx bx kx F t. dt IR I LI V t, Q LQ RQ V t, Lecture 5 omplex Variables II (Applicatios i Physics) (See hapter i Boas) To see why complex variables are so useful cosider first the (liear) mechaics of a sigle particle described by Newto s equatio

More information

Chapter 9: Numerical Differentiation

Chapter 9: Numerical Differentiation 178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical

More information

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.

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. 0_0905.qxd //0 :7 PM Page SECTION 9.5 Alteratig Series Sectio 9.5 Alteratig Series Use the Alteratig Series Test to determie whether a ifiite series coverges. Use the Alteratig Series Remaider to approximate

More information

TIME-CORRELATION FUNCTIONS

TIME-CORRELATION FUNCTIONS p. 8 TIME-CORRELATION FUNCTIONS Time-correlatio fuctios are a effective way of represetig the dyamics of a system. They provide a statistical descriptio of the time-evolutio of a variable for a esemble

More information

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

L = n i, i=1. dp p n 1 Exchageable sequeces ad probabilities for probabilities 1996; modified 98 5 21 to add material o mutual iformatio; modified 98 7 21 to add Heath-Sudderth proof of de Fietti represetatio; modified 99 11

More information

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

Quantum Mechanics I. 21 April, x=0. , α = A + B = C. ik 1 A ik 1 B = αc. Quatum Mechaics I 1 April, 14 Assigmet 5: Solutio 1 For a particle icidet o a potetial step with E < V, show that the magitudes of the amplitudes of the icidet ad reflected waves fuctios are the same Fid

More information

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

LECTURE 14. Non-linear transverse motion. Non-linear transverse motion LETURE 4 No-liear trasverse motio Floquet trasformatio Harmoic aalysis-oe dimesioal resoaces Two-dimesioal resoaces No-liear trasverse motio No-liear field terms i the trajectory equatio: Trajectory equatio

More information

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

HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples: 5.6 4 Lecture #3-4 page HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT Do t restrict the wavefuctio to a sigle term! Could be a liear combiatio of several wavefuctios e.g. two terms:

More information

The Minimum Distance Energy for Polygonal Unknots

The Minimum Distance Energy for Polygonal Unknots The Miimum Distace Eergy for Polygoal Ukots By:Johaa Tam Advisor: Rollad Trapp Abstract This paper ivestigates the eergy U MD of polygoal ukots It provides equatios for fidig the eergy for ay plaar regular

More information

Sequences I. Chapter Introduction

Sequences I. Chapter Introduction Chapter 2 Sequeces I 2. Itroductio A sequece is a list of umbers i a defiite order so that we kow which umber is i the first place, which umber is i the secod place ad, for ay atural umber, we kow which

More information

Chapter 2 Feedback Control Theory Continued

Chapter 2 Feedback Control Theory Continued Chapter Feedback Cotrol Theor Cotiued. Itroductio I the previous chapter, the respose characteristic of simple first ad secod order trasfer fuctios were studied. It was show that first order trasfer fuctio,

More information

CALCULATION OF FIBONACCI VECTORS

CALCULATION OF FIBONACCI VECTORS CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College

More information

a b c d e f g h Supplementary Information

a b c d e f g h Supplementary Information Supplemetary Iformatio a b c d e f g h Supplemetary Figure S STM images show that Dark patters are frequetly preset ad ted to accumulate. (a) mv, pa, m ; (b) mv, pa, m ; (c) mv, pa, m ; (d) mv, pa, m ;

More information

Marcinkiwiecz-Zygmund Type Inequalities for all Arcs of the Circle

Marcinkiwiecz-Zygmund Type Inequalities for all Arcs of the Circle Marcikiwiecz-ygmud Type Iequalities for all Arcs of the Circle C.K. Kobidarajah ad D. S. Lubisky Mathematics Departmet, Easter Uiversity, Chekalady, Sri Laka; Mathematics Departmet, Georgia Istitute of

More information

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

Diffusivity and Mobility Quantization. in Quantum Electrical Semi-Ballistic. Quasi-One-Dimensional Conductors Advaces i Applied Physics, Vol., 014, o. 1, 9-13 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/aap.014.3110 Diffusivity ad Mobility Quatizatio i Quatum Electrical Semi-Ballistic Quasi-Oe-Dimesioal

More information