Non-interacting Spin-1/2 Particles in Non-commuting External Magnetic Fields

Similar documents
PHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

ψ ij has the eigenvalue

A how to guide to second quantization method.

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

1 Matrix representations of canonical matrices

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

Lecture 14: Forces and Stresses

5.03, Inorganic Chemistry Prof. Daniel G. Nocera Lecture 2 May 11: Ligand Field Theory

Quantum spin system with on-site exchange in a magnetic field

The non-negativity of probabilities and the collapse of state

Density matrix. c α (t)φ α (q)

GREENBERGER- HORNE- ZEILINGER (GHZ) STATES IN QUANTUM DOT MOLECULE

Advanced Quantum Mechanics

ECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

One-sided finite-difference approximations suitable for use with Richardson extrapolation

2.3 Nilpotent endomorphisms

PHYS 705: Classical Mechanics. Canonical Transformation II

Quantum Mechanics I - Session 4

The exponential map of GL(N)

V.C The Niemeijer van Leeuwen Cumulant Approximation

Formulas for the Determinant

PSEUDO ORBIT EXPANSION FOR THE RESONANCE CONDITION ON QUANTUM GRAPHS AND THE RESONANCE ASYMPTOTICS

Difference Equations

24. Atomic Spectra, Term Symbols and Hund s Rules

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:

ρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

P A = (P P + P )A = P (I P T (P P ))A = P (A P T (P P )A) Hence if we let E = P T (P P A), We have that

Change. Flamenco Chuck Keyser. Updates 11/26/2017, 11/28/2017, 11/29/2017, 12/05/2017. Most Recent Update 12/22/2017

THEOREMS OF QUANTUM MECHANICS

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

A Hybrid Variational Iteration Method for Blasius Equation

Supplemental document

On the correction of the h-index for career length

DETAILED DERIVATION OF THE GENERAL MASTER EQUATION IN QUANTUM OPTICS

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

arxiv:quant-ph/ Jul 2002

Norm Bounds for a Transformed Activity Level. Vector in Sraffian Systems: A Dual Exercise

Note on the Electron EDM

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

Susceptibility and Inverted Hysteresis Loop of Prussian Blue Analogs with Orthorhombic Structure

A property of the elementary symmetric functions

An Introduction to Morita Theory

5 The Rational Canonical Form

Representation theory and quantum mechanics tutorial Representation theory and quantum conservation laws

Structure and Drive Paul A. Jensen Copyright July 20, 2003

14 The Postulates of Quantum mechanics

Foundations of Arithmetic

Linear Approximation with Regularization and Moving Least Squares

Solution 1 for USTC class Physics of Quantum Information

1 Vectors over the complex numbers

12. The Hamilton-Jacobi Equation Michael Fowler

8.6 The Complex Number System

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

CHAPTER 14 GENERAL PERTURBATION THEORY

8. Superfluid to Mott-insulator transition

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

Indeterminate pin-jointed frames (trusses)

arxiv: v2 [quant-ph] 29 Jun 2018

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

arxiv:cond-mat/ v2 [cond-mat.mes-hall] 3 Jan 2006

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

h-analogue of Fibonacci Numbers

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

A FORMULA FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIAGONAL MATRIX

SUPPLEMENTARY INFORMATION

Numerical Heat and Mass Transfer

Open string operator quantization

Dynamics of a Superconducting Qubit Coupled to an LC Resonator

Canonical transformations

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution

LECTURE 9 CANONICAL CORRELATION ANALYSIS

Supplementary Materials for

MMA and GCMMA two methods for nonlinear optimization

EEE 241: Linear Systems

1 Rabi oscillations. Physical Chemistry V Solution II 8 March 2013

Quantum Mechanics for Scientists and Engineers. David Miller

Randić Energy and Randić Estrada Index of a Graph

One-Dimensional Ising Chain with Competing Interactions: Exact Results and Connection with Other Statistical Models

Linear system of the Schrödinger equation Notes on Quantum Mechanics

Section 8.3 Polar Form of Complex Numbers

RG treatment of an Ising chain with long range interactions

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

Army Ants Tunneling for Classical Simulations

Affine transformations and convexity

More metrics on cartesian products

Explicit constructions of all separable two-qubits density matrices and related problems for three-qubits systems

On the set of natural numbers

) is the unite step-function, which signifies that the second term of the right-hand side of the

Rate of Absorption and Stimulated Emission

Large distance asymptotic behavior of the emptiness formation probability of the XXZ spin- 1 2 Heisenberg chain

A new Approach for Solving Linear Ordinary Differential Equations

Transcription:

EJTP 6, No. 0 009) 43 56 Electronc Journal of Theoretcal Physcs Non-nteractng Spn-1/ Partcles n Non-commutng External Magnetc Felds Kunle Adegoke Physcs Department, Obafem Awolowo Unversty, Ile-Ife, Ngera Receved 6 July 008, Accepted 13 January 009, Publshed 0 February 009 Abstract: We obtan, n one dmenson, all the energy levels of a system of non-nteractng spn-1/ partcles n non-commutng external magnetc felds. Examples of how to ncorporate nteractons as perturbatons are gven for the Isng model n two orthogonal felds and for the XZ model n two orthogonal felds. c Electronc Journal of Theoretcal Physcs. All rghts reserved. Keywords: Quantum Spn Systems; Isng Systems; Magnetc Felds; Nonlnear Dynamcs and Nonlnear Dynamcal Systems; Magnetc Ansotropy; Numercal Smulaton Studes PACS 008): 75.10.Hk; 75.30.Ds; 05.45.-a 1. Introducton Hamltonan spn models nvolvng two external non-commutng magnetc felds are beng ncreasngly studed these days [1,, 3]. It s almost always the case that such models cannot be solved exactly, n the presence of spn nteractons. Interacton s usually due to ether nearest neghbour exchange or next nearest neghbour exchange or both. As s common practce, approxmate calculatons can be done by frst startng wth an exactly solvable model and then ntroducng the nteractons as perturbatons. Ths s done n ths paper. We start by consderng the smple model descrbed by the Hamltonan H = h x N =1 S x h z N S z, 1) =1 where S x and S z are spn-1/ operators on the th lattce ste of a one-dmensonal chan of N partcles. h x and h z are the external transverse and longtudnal magnetc felds, respectvely, measured n unts where the splttng factor and Bohr magneton are unty. adegoke@daad-alumn.de and adegoke@rushpost.com

44 Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 It wll often be convenent to wrte H = N H, =1 where H = h x S x h zs z. The model 1) s exactly solvable and n ths work we wll obtan all the energy levels and the correspondng states. The man mport of the result wll be that nteractons can be ncluded as perturbatons on the model 1), ths approach beng effectve n obtanng accurate results for models ncorporatng two external felds.. The One-partcle Model Snce the N spn-1/ partcles descrbed by 1) are non-nteractng, all results can be obtaned from the Hamltonan for a sngle partcle. We drop the ste subscrpts n 1) and wrte H ε for the one spn system and wrte H ε = h x S x h z S z, where S x and S z are smply the spn-1/ operators S x = 1 01,S z = 1 1 0 10 0 1 n a bass wth S z dagonal and wth = 1. We can denote the bass states n the S z bass by the set {, } andnthes x bass by {, }, sothat S z = 1 S z = 1 S x = 1 S x = 1. S x s obtaned from S z by the smlarty transformaton S x = PS z P 1, where P s the untary matrx P = 1 1 1 = P 1. 1 1

Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 45 Thus, the S x bass states and the S z bass states are related by = 1 + ) = 1 ), wth the nverse relaton = 1 + ) = 1 ). The Hamltonan of the one partcle system n the S z bass s The normalzed egenstates of H ε are wth egenenergy h x H ε = 1 h z. h x h z h x + ε 0 = h x + h z + ) h x + h z h z + ), h x + h z and wth egenenergy ε 1 = h ε 0 = x + h z, h x + h x + h z ) h x + h z h z ), h x + h z h ε 1 =+ x + h z. Now that we have obtaned the egenstates and the correspondng energes of the one-partcle Hamltonan, we return to the general Hamltonan 1). 3. Energy Levels and Degeneraces From the results of the prevous secton and the fact that the partcles are non-nteractng, t s clear that there are N + 1 energy levels. We obtan them presently.

46 Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 3.1 Ground State The ground state E 0 of the system of non-nteractng N spn-1/ partcles s the drect product state E 0 = ε 0 ε 0 ε 0 ε 0 = ε 0 ) N. ) The ground state energy of the system s found as follows H E 0 = H E 0 =H 1 ε 0 ) ε 0 ε 0 ε 0 + ε 0 H ε 0 ) ε 0 ε 0 + ε 0 ε 0 H 3 ε 0 ) ε 0 ε 0 + + ε 0 ε 0 ε 0 ε 0 H N ε 0 ). Snce H ε 0 = ε 0 ε 0 and there are N terms n the above sum, we thus fnd that H E 0 = Nε 0 E 0. The ground state energy of the model s therefore E 0 = Nε 0 = N h x + h z. We shall have more to say about the ground state, n the next secton, but presently we obtan the other energy levels. 3. Excted States A frst excted state FES) of the model s E 1 1 = ε 0 ε 0 ε 0 ε 0 ε 1. We see at once that the FES s N fold degenerate, snce ε 1 can occur anywhere n the drect product state. The superscrpt 1 was affxed n antcpaton. The remanng N 1 states that are degenerate wth E 1 1 are E 1 = ε 0 ε 0 ε 0 ε 1 ε 0 E 1 3 = ε 0 ε 0 ε 1 ε 0 ε 0 E 1 N = ε 1 ε 0 ε 0 ε 0 ε 0.

Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 47 The energy of the frst excted state s found from H E 1 1 = H E 1 1 =H 1 ε 0 ) ε 0 ε 0 ε 0 ε 1 + ε 0 H ε 0 ) ε 0 ε 0 ε 1 + + ε 0 ε 0 ε 0 H N ε 1 ) From whch t follows that =N 1)ε 0 E 1 1 + ε 1 E 1 1 =[N 1)ε 0 + ε 1 ] E 1 1. E 1 =N 1)ε 0 + ε 1 =N )ε 0. In general, a kth excted state has k ε 1 factors n the drect product state. Snce there are N,k) ways of arrangng the k ε 1 factors, ths means that the kth excted state has degeneracy g E k ), gven by g E k )= N k = N! k!n k)!. The degenerate energy E k s easly found by the above scheme to be E k =N k)ε 0 + kε 1 =N k)ε 0. We see that only the ground state E 0 and the Nth excted state E N are nondegenerate. 4. An Explct Expresson for the Ground State We recall from equaton ) that the ground state s the drect product state where ε 0 = E 0 = ε 0 ) N, 3) h x + h z + ) h x + h z h x h + z + ). h x + h z

48 Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 Performng the bnomal expanson suggested by equaton 3), we have that E 0 = N h N m x h z + m h x z) + h Sm [ h x + h z + ) ] N/, 4) h x + h z m=0 where S m s the lnear combnaton of the N,m) states wth m spns up n the total S z bass), that s, the lnear combnaton of all states wth total S z = m N/. We can check the trval lmts of 4): 1) h z =0 H = h x ) h x =0 that s Puttng h z = 0 n 4) gves E 0 = N m=0 = 1 ) N S x h N m x h m x S m h x )N/ N S m m=0 = 1 + )N ) N = ) N = H = h z Here the ground state s trvally the all spns up, non-degenerate ferromagnetc state E 0 = S z wth energy E 0 = Nh z / Puttng h x = 0 n n 4), we see that the only non-vanshng term n the sum occurs when m = N and so, E 0 = h z) N h z ) S N N =.

Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 49 5. Example Applcatons In ths secton, we wll consder two models wth nteracton, the nteractons wll be treated as perturbatons. Specfcally we wll obtan the ground state energy of the Isng model n two felds and also the XZ model n non-commutng felds to second order n nearest neghbour exchange nteractons. We wll also compute frst order correctons to the energy of the frst excted state for each model. We frst note that for the one-partcle system and S z ε 0 = 1 h x + h z + ) h x + h z h x h + z + ) h x + h z = α S z ε 1 = 1 h x + h x + = β, so that we have the matrx elements h z ) h x + h z h z ) h x + h z and Furthermore, ε 0 S z ε 0 = ε 0 α = 1 h z h x+h z ε 1 S z ε 0 = ε 1 α = 1 h x = ε 0 β h x+h z ε 1 S z ε 1 = ε 1 β = 1 h z. h x+h z S x ε 0 = 1 h x + h z + ) h x + h z h x h + z + ) h x + h z S x ε 1 = 1 h x + h z ) h x + h z h x h + z ), h x + h z 5) so that we have also the followng marx elements: ε 0 S x ε 0 = ε 0 λ = 1 h x = ε 1 α h x +h z ε 1 S x ε 0 = ε 1 λ = 1 h z = ε 0 S x ε 1 = ε 0 δ h x +h z ε 1 S x ε 1 = ε 1 δ = 1 h x = ε 1 α = ε 0 λ h x +h z

50 Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 We are now ready to compute the energy correctons. 5.1 Perturbaton Expanson of the Ground State Energy of the Isng Model n Non-Commutng Felds The Isng model n non-commutng magnetc felds s descrbed by the Hamltonan where H I = j S z S+1 z h x S x h z S z = V + H, 6) V = j S z S z +1 and H s as gven n 1). The nearest neghbour exchange nteracton j s assumed postve, so that we have a ferromagnetc model. Furthermore, we assume a perodc boundary condton, so that SN+ z = Sz. For weak nteracton, V can be treated as a perturbaton. 5.1.1 Frst Order Energy Correcton The frst order correcton to the ground state energy of the Hamltonan 6) s gven by ΔE 1) 0 I = E 0 V E 0 = j E 0 S z Sz +1 E 0. Now, S z S+1 z E 0 = S1S z z E 0 + SS z 3 z E 0 + + SN z Sz 1 E 0 = α α ε 0 ε 0 ε 0 + ε 0 α α ε 0 ε 0 + ε 0 ε 0 α α ε 0 7) + + α ε 0 ε 0 ε 0 α. Multplyng from the left by E 0,wehave E 0 S z Sz +1 E 0 = N ε 0 α ) = N 4 h z h x +. h z

Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 51 We therefore have that the frst order correcton to the ground state energy of the Isng model n non-commutng magnetc felds s 5.1. Second Order Correcton ΔE 1) 0 I = jn 4 h z h x +. h z The second order correcton to the ground state energy of the Isng model n noncommutng felds has the form ΔE ) 0 I = E E V E 0 E 0 E. 8) From the form of 7) and the fact that ε 0 ε 1 = 0 t s clear that only states wth E = E 1 and those wth E = E can contrbute to the sum n 8). We consder them n turns. When E = E 1,weobservethat E 1 1 = ε 0 ε 0 ε 0 ε 1 gves 1 E 1 S z Sz +1 E 0 = ε 1 ε 0 ε 0 ε 0 ε 0 ε 0 ε 0 ε 0 α α + ε 1 ε 0 ε 0 ε 0 α ε 0 ε 0 ε 0 ε 0 α = ε 1 α ε 0 α. We have smlar results from the remanng N 1 states that are degenerate wth E 1 1.TheE = E 1 states therefore contrbute N ε 1 α ε 0 α ) j) ε 0 = Nj 8ε 0 h x h z h x + h z) to the sum n 8). We have used 5) to substtute for ε 1 α and ε 0 α and we have also used the fact that E 0 E k =kε 0. As for the states wth E = E,onlyN those wth two consecutve ε 1 factors) of the N,) states contrbute to the sum. Ther contrbuton s N ε 1 α ) j) 4ε 0 = Nj h 4 x 64ε 0 h x + h z). The second order correcton to the ground state energy of the Isng model n noncommutng felds s therefore ΔE ) 0 I = Nj 8ε 0 h x h z h x + h z) + Nj 64ε 0 h 4 x h x + h z).

5 Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 5. Perturbaton expanson of the ground state energy of the XZ model n non-commutng felds. The XZ model n one dmenson, n the presence of non-commutng external felds s descrbed by the Hamltonan H xz = j 1 N =1 S x S x +1 j N =1 S z S z +1 h x N =1 S x h z where j 1 > 0andj > 0 are the nearest neghbour exchange nteracton energy. Agan we assume a perodc boundary condton, so that SN+ x = Sx and SN+ z = Sz. For weak nearest neghbour nteractons, after calculatons smlar to that n the prevous secton, we obtan the followng correctons to the ground state energy: ΔE 1) 0 XZ = N 1 4 h x + j1 h h x + j hz) z and ΔE ) 0 XZ = N h x h ) z j 8ε 0 h x + 1 + j h z ) + N 1 j 64ε 0 h x + h z) 1 h 4 x + ) j h4 z. The ground state energy of the XZ model n non-commutng felds, to second order n nteractons j 1 and j s thus E 0XZ = N h x + h z N 1 j1 h 4 h x + h x + j hz) z + N 8ε 0 h x h z h x + h z) j 1 + j + N 64ε 0 1 h x + h z ) j 1 h 4 x + j h4 z 5.3 Frst Order Correcton to the Energy of the Frst Excted State FES) of the Isng Model n Non-Commutng Felds The frst excted state, n the absence of nteractons, s N fold degenerate, wth the N states as gven n 3.). We therefore use degenerate perturbaton theory to fnd the frst order correcton to the frst excted state energy of the Isng model n non-commutng felds. S z S+1 z E 1 1 = α α ε 0 ε 0 ε 1 + ε 0 α α ε 0 ε 0 ε 1 + + ε 0 ε 0 ε 0 ε 0 α β + α ε 0 ε 0 ε 0 ε 0 β. ) ). N =1 S z,

Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 53 We therefore have that 1 E 1 S z Sz +1 E 1 1 =N ) ε 0 α ) + ε 0 α ε 1 β ) = h z N h x + h z 4 1. A smlar calculaton for the remanng N 1 states shows that the dagonal elements of the N N perturbaton matrx V, are equal, beng gven by: V = E 1 j ) S z Sz +1 E 1 ) = jh z N h x + h z 4 1 = f. 9) We can also compute V 1 thus V 1 = E 1 j ) S z Sz +1 E 1 1 = j ε 0 ε 1 ε 0 ε 0 ε 0 α β = j ε 1 α ε 0 β = j h x 4 h x +. h z In fact, for m n, wehave V mn = m E 1 j ) S z Sz +1 E 1 n = j h x δ 4 h x + h m,n±1 z = V nm = g. 10) The N N perturbaton matrx V s therefore fg00 0 0 g gfg0 0 0 0 0 gfg 0 0 0 00gf g 0 0 V = 000g f g.. 000... g f g g 000 0 g f

54 Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 where f and g are as gven n 9) and 10). Clearly f and g are negatve for N>4. The smallest egenvalue of V s g + f and s non-degenerate, so that the correcton to the frst excted state energy, to frst order n nteracton j s ΔE 1) 1 I =g + f = j h x h x + h z jh z h x + h z ) N 4 1. The nearest neghbour nteractons thus lft the degeneracy n the frst excted state of the model. The energy of the frst excted state of the Isng model n non-commutng felds, to frst order n j s therefore E 1I = 1 N ) h x + h z j h x h x + h z jh z h x + h z ) N 4 1. 5.4 Frst Order Correcton to the Energy of the Frst Excted State FES) of the XZ Model n Non-Commutng Felds A completely analogous treatment to the above gves the energy of the frst excted state of the XZ model n non-commutng felds, to frst order n the nteractons j 1 and j as E 1XY ) N h = 1 x + h z 1 j h x + j 1 h z h x + h z N 4 1 ) j1 h x + j h z h x + h z. Concluson We have obtaned all the energy levels of a system of N non-nteractng spn-1/ partcles n non-commutng external magnetc felds, n one dmenson. The energy of the kth level s E k =N k)ε 0,whereε 0 = h x + h z )/ andk =0, 1,,...,N. The kth energy level was found to be N,k) = N!/[k!N k)!]-fold degenerate. An expct expresson for the ground state was also derved equaton 4)). It s not dffcult to gve smlar explct expressons for the excted states. A scheme for obtanng the energes of the excted states was hghlghted n secton 3.. Snce most real systems do nteract, examples of how the nteracton terms of an Hamltonan can be ncluded by Raylegh-Schrödnger perturbaton was demonstrated by fndng energy correctons to the ground state energy of the Isng model n non-commutng felds, as well as to the ground state energy of the XZ model n non-commutng felds. Correctons to the frst excted state energy were also calculated.

Electronc Journal of Theoretcal Physcs 6, No. 0 009) 43 56 55 Acknowledgments The author s grateful to the DAAD for a scholarshp and thanks the Physcs Insttute, Unverstät Bayreuth for hosptalty. References [1] D.V.DmtrevandV.Y.Krvnov,Phys.Rev.B70, 144414 004). [] M. Kenzelmann, R. Coldea, D. A. Tennant, D. Vsser, M. Hofmann, P. Smebdl, and Z. Tylczynsk, Phys. Rev B. 65, 14443 00). [3] A. Ovchnnkov, D. V. Dmtrev, V. Y. Krvnov, and V. O. Cheranovsk, Phys. Rev. B 68, 14406 003).