6 Symmetry Properties of Electronic and Photonic Band Structures W.Hergert,M.Dane,and D.Kodderitzsch Martin-Luther-University Halle-Wittenberg,Department of Physics, Von-Seckendorff-Platz 1,06120 Halle,Germany Abstract.Group theoretical investigations have a huge potential to simplify calcu- lations in solid state theory.We will discuss the application of group theory to elec- tronic and photonic band structures.The symmetry properties of the Schrodinger equation and Maxwell's equations as well will be investigated.We have developed methods to simplify group theoretical investigations based on the computer algebra system Mathematica. 6.1 Introduction The majority of physical systems exhibit intrinsic symmetries which can be used to simplify the solution of the equations governing these systems.Group theory as a mathematical tool plays an important role to classify the solu- tions within the context of the underlying symmetries.Extensive use of group theory has been made to simplify the study of electronic structure or vibra- tional modes of solids or molecules.There exists a number of excellent books illustrating the use of group theory.[6.1-6.4] Similarities between the solution of the Schrodinger equation,or the Kohn-Sham equations in the framework of density functional theory for a crystal,and the solution of Maxwell's equations have been already pointed out by Joannopoulos et al..6.5 Therefore,it is clear that the same group theoretical concepts like in the theory of electronic band structures including two-dimensional and three-dimensional structures,surfaces as well as defects should be applicable to photonic band structure calculations,if we take into account the vectorial nature of the electromagnetic field. There are several publications about group theoretical investigations of photonic crystals in the literature.Sakoda has extensively studied the sym- metry properties of two-dimensional and three-dimensional photonic crystals (cf.[6.6-6.9),starting from a plane wave representation of the electromag- netic fields.Group theoretical investigations are done also by Ohtaka and Tanabe.[6.10]They investigate the symmetry properties of photonic crystals represented by an array of dielectric spheres.In this case a series expansion in terms of vector spherical harmonics is used. The problem at the end for electronic and the photonic band structure calculations is:How to apply group theory in the actual research work,if one goes away from all the textbook examples.The aim of the paper is to show, W.Hergert,M.Dane,and D.Kodderitzsch,Symmetry Properties of Electronic and Photonic Band Structures,Lect.Notes Phys.642,103-125(2004) http://www.springerlink.com/ C Springer-Verlag Berlin Heidelberg 2004
6 Symmetry Properties of Electronic and Photonic Band Structures W. Hergert, M. D¨ane, and D. K¨odderitzsch Martin-Luther-University Halle-Wittenberg, Department of Physics, Von-Seckendorff-Platz 1, 06120 Halle, Germany Abstract. Group theoretical investigations have a huge potential to simplify calculations in solid state theory. We will discuss the application of group theory to electronic and photonic band structures. The symmetry properties of the Schr¨odinger equation and Maxwell’s equations as well will be investigated. We have developed methods to simplify group theoretical investigations based on the computer algebra system Mathematica. 6.1 Introduction The majority of physical systems exhibit intrinsic symmetries which can be used to simplify the solution of the equations governing these systems. Group theory as a mathematical tool plays an important role to classify the solutions within the context of the underlying symmetries. Extensive use of group theory has been made to simplify the study of electronic structure or vibrational modes of solids or molecules. There exists a number of excellent books illustrating the use of group theory. [6.1–6.4] Similarities between the solution of the Schr¨odinger equation, or the Kohn-Sham equations in the framework of density functional theory for a crystal, and the solution of Maxwell’s equations have been already pointed out by Joannopoulos et al.. [6.5] Therefore, it is clear that the same group theoretical concepts like in the theory of electronic band structures including two-dimensional and three-dimensional structures, surfaces as well as defects should be applicable to photonic band structure calculations, if we take into account the vectorial nature of the electromagnetic field. There are several publications about group theoretical investigations of photonic crystals in the literature. Sakoda has extensively studied the symmetry properties of two-dimensional and three-dimensional photonic crystals (cf. [6.6–6.9]), starting from a plane wave representation of the electromagnetic fields. Group theoretical investigations are done also by Ohtaka and Tanabe. [6.10] They investigate the symmetry properties of photonic crystals represented by an array of dielectric spheres. In this case a series expansion in terms of vector spherical harmonics is used. The problem at the end for electronic and the photonic band structure calculations is: How to apply group theory in the actual research work, if one goes away from all the textbook examples. The aim of the paper is to show, W. Hergert, M. D¨ane, and D. K¨odderitzsch, Symmetry Properties of Electronic and Photonic Band Structures, Lect. Notes Phys. 642, 103–125 (2004) http://www.springerlink.com/ c Springer-Verlag Berlin Heidelberg 2004
104 W.Hergert,M.Dane,and D.Kodderitzsch that computer algebra tools are appropriate to simplify group theoretical discussions connected with the calculation of electronic and photonic band structures. After a short discussion of the usefulness of computer algebra systems, we will introduce basic concepts of group theory.Representation theory will be discussed next,followed by the analysis of the symmetry properties of the Schrodinger equation and of Maxwell's equations.We want to solve the electronic and photonic problem for a lattice periodic situation.In the one case we have a periodic potential V(r)=V(r+R)in the other case a periodic dielectric constant e(r)=e(r+R).The consequences of that periodicity will be considered.A simplification of the solution of both kinds of problems is possible,if symmetry-adapted basis functions are used.This is discussed in more detail for the electronic problem.At the end we apply group theory to the calculation of photonic band structures. 6.2 Group Theory Packages for Computer Algebra Systems Computer algebra(CA)systems like Mathematica or Maple have been devel- oped to allow formal mathematical manipulations.Nowadays those systems are complete in such a sense,that formal manipulations,numerical calcula- tions,as well as graphical representations are possible in an easy and intuitive way with the same software.Apart from such general purpose CA systems, there exist also systems developed for special applications in mathematics. Group theoretical considerations,although conceptionally easy,lead very of- ten to time-consuming algebraic calculations,which are error-prone.There- fore group theory is an excellent field for the application of CA systems.Some systems for abstract group theory are available [6.11],but they are not very helpful for considerations in solid state theory.K.Shirai developed a Math- ematica package for group theory in solid state physics [6.12].We followed similar ideas but tried to make the package more easy to use. We have constructed a package for the CA system Mathematica which allows to do group theoretical manipulations which occur in solid state theory. 6.13]The software allows basic considerations with point groups,contains tight-binding theory for the electronic structure of solids,but also special applications to photonic crystals.The package is accompanied by an on line help,which is integrated in Mathematica's help system.All considerations discussed in this paper can be found in a Mathematica notebook which is part of our package.We will give references to the commands implemented in the package throughout the paper
104 W. Hergert, M. D¨ane, and D. K¨odderitzsch that computer algebra tools are appropriate to simplify group theoretical discussions connected with the calculation of electronic and photonic band structures. After a short discussion of the usefulness of computer algebra systems, we will introduce basic concepts of group theory. Representation theory will be discussed next, followed by the analysis of the symmetry properties of the Schr¨odinger equation and of Maxwell’s equations. We want to solve the electronic and photonic problem for a lattice periodic situation. In the one case we have a periodic potential V (r) = V (r+R) in the other case a periodic dielectric constant (r) = (r + R). The consequences of that periodicity will be considered. A simplification of the solution of both kinds of problems is possible, if symmetry-adapted basis functions are used. This is discussed in more detail for the electronic problem. At the end we apply group theory to the calculation of photonic band structures. 6.2 Group Theory Packages for Computer Algebra Systems Computer algebra (CA) systems like Mathematica or Maple have been developed to allow formal mathematical manipulations. Nowadays those systems are complete in such a sense, that formal manipulations, numerical calculations, as well as graphical representations are possible in an easy and intuitive way with the same software. Apart from such general purpose CA systems, there exist also systems developed for special applications in mathematics. Group theoretical considerations, although conceptionally easy, lead very often to time-consuming algebraic calculations, which are error-prone. Therefore group theory is an excellent field for the application of CA systems. Some systems for abstract group theory are available [6.11], but they are not very helpful for considerations in solid state theory. K. Shirai developed a Mathematica package for group theory in solid state physics [6.12]. We followed similar ideas but tried to make the package more easy to use. We have constructed a package for the CA system Mathematica which allows to do group theoretical manipulations which occur in solid state theory. [6.13] The software allows basic considerations with point groups, contains tight-binding theory for the electronic structure of solids, but also special applications to photonic crystals. The package is accompanied by an on line help, which is integrated in Mathematica’s help system. All considerations discussed in this paper can be found in a Mathematica notebook which is part of our package. We will give references to the commands implemented in the package throughout the paper
6 Symmetry Properties of Electronic and Photonic Band Structures 105 6.3 Basic Concepts in Group Theory In this section we will introduce the basic concepts of group theory.We will focus on definitions which will be necessary for the following discussions.To illustrate the concepts,let us discuss the symmetry group of a square in two dimensions,i.e.we are interested in all transformations which transform the square into itself.(cf.Fig.6.1)This example will be of later use for the discussion of photonic band structures. A set g ofelements A.B,C...is called a group if the following four axioms are fulfilled:i)There exists an operation,often called multiplication,which associates every pair of elements of g with another element of g:A Eg,BE g→A·B=C,C∈g,ii)The associative law is valid:A,B,C∈g→ (A·B).C=A·(B.C)=A·B.C,iii)in the set exists an identity element: A,E∈g→A·E=E·A=A,iv)for all A∈G exists an inverse element A-1∈G with A·A-1=A-1·A=E. The symmetry group of the square consists of rotations of /2 around the z-axis and mirror operations.The normal vectors of the mirror planes are the x and y-axis,Oa and Ob.A mirror symmetry may be expressed as a twofold rotation,followed by an inversion.All the symmetry operations of the group of the square,named C4v are: CAv=[E,Cz,CAz,Cz,ICx,IC2y:IC2a,IC2b} (6.1) If the group theory package is included in a Mathematica notebook by means of the command Needs ["GroupTheory'Master'"]1,the newly defined com- mand c4v=InstallGroup["C4v"]will install the group in terms of the ro- tation matrices of the elements.The matrices are stored in the list c4v.The command c4vs=GetSymbol[c4v]will transform the elements of the group into symbolic form (cf.(6.1)).The symbols will be stored in the list c4vs. Operations on the group can be done in both representations of the group el- ements.The group multiplication is implemented in our Mathematica package by a redefinition of the infix-operator (see 6.25). Fig.6.1.Symmetry of a square. 1 Mathematica commands are marked using this particular font
6 Symmetry Properties of Electronic and Photonic Band Structures 105 6.3 Basic Concepts in Group Theory In this section we will introduce the basic concepts of group theory. We will focus on definitions which will be necessary for the following discussions. To illustrate the concepts, let us discuss the symmetry group of a square in two dimensions, i.e. we are interested in all transformations which transform the square into itself. (cf. Fig. 6.1) This example will be of later use for the discussion of photonic band structures. A set G of elements A,B,C ... is called a group if the following four axioms are fulfilled: i) There exists an operation, often called multiplication, which associates every pair of elements of G with another element of G: A ∈ G, B ∈ G → A · B = C, C ∈ G, ii) The associative law is valid: A, B, C ∈G→ (A · B)· C = A ·(B · C) = A · B · C, iii) in the set exists an identity element: A, E ∈G→ A · E = E · A = A, iv) for all A ∈ G exists an inverse element A−1 ∈ G with A · A−1 = A−1 · A = E. The symmetry group of the square consists of rotations of π/2 around the z-axis and mirror operations. The normal vectors of the mirror planes are the x and y-axis, −→ Oa and −→ Ob. A mirror symmetry may be expressed as a twofold rotation, followed by an inversion. All the symmetry operations of the group of the square, named C4v are: C4v = E,C2z, C4z, C−1 4z ,IC2x,IC2y,IC2a,IC2b (6.1) If the group theory package is included in a Mathematica notebook by means of the command Needs["GroupTheory‘Master‘"]1, the newly defined command c4v=InstallGroup["C4v"] will install the group in terms of the rotation matrices of the elements. The matrices are stored in the list c4v. The command c4vs=GetSymbol[c4v] will transform the elements of the group into symbolic form (cf.(6.1)). The symbols will be stored in the list c4vs. Operations on the group can be done in both representations of the group elements. The group multiplication is implemented in our Mathematica package by a redefinition of the infix-operator ⊕ (see [6.25]). Fig. 6.1. Symmetry of a square. 1 Mathematica commands are marked using this particular font
106 W.Hergert,M.Dane,and D.Kodderitzsch For the graphical user interface of Mathematica we designed two addi- tional palettes.The palette SymmetryElements contains all the symmetry elements of the 32 point groups.The palette PointGroup contains all the operations defined in the package. The inspection of the multiplication table(MultiplicationTable[c4v]) of the group shows,that the table is not symmetric,indicating that the group multiplication is not commutative in this group.Therefore the group is not abelian (AbelianQ [c4vs]returns False.).All elements of the group can be also generated by successive multiplication of only two elements of the group, the so called generators.For our example we find the generators C4z,IC2x. (Generators[group_])2 Any set of elements of a group which obeys all the group postulates is called a subgroup of the group.Examples for subgroups are:[E,CCzC and [E,IC2y}(SubGroupQ[group_,subgroup_]performs a test). Another structure of the group is introduced by the definition of conjugate elements.An element B of the group g is said to be conjugate to A if their exists a group element X such that B=XAX-1.A class(classes [c4vs]) is a collection of mutually conjugate elements of a group.The classes of the group Cav are:(E),(Cz),(Cz,C),(ICx,ICy),(IC2a,ICb).It is easy to prove that:i)E always forms a class on its own,ii)every group element in g is a member of some class of g,iii)no element can be the member of two classes of g,iv)if g is an abelian group,every element forms a class on its own. We have already seen,that the group in the example can be represented by the abstract symbols or the rotation matrices.The relation between groups is formulated more rigorously in terms of homomorphism and isomorphism. Two groups g=[A,B,C,...}and C'=[A',B',C',...}are called homo- morphic,if i)to each element of g corresponds one and only one element of G'ii)for all elements holds A-A',B-B'AB-A'B'.The set of elements which are mapped to E'E g in a homomorphism is called ker- nel of the homomorphism Nk.In contrast to the homomorphism the map- ping in the isomorphism is one-to-one.Whereas the mapping of the symbols on the rotation matrices constitutes an isomorphism,the following mapping E,Cz,Caz,C+1 and IC2x;IC2y:IC2a:IC2b-1 constitutes a ho- momorphism. 6.4 Representation Theory 6.4.1 Matrix Representations of Groups Matrir representations of symmetry groups are the essential tools to inves- tigate symmetry properties of solutions of field equations as Schrodinger's equation or Maxwell's equations. 2 The formal arguments of the commands indicate,which type of information has to be plugged in.For more information use the help system or the example notebook
106 W. Hergert, M. D¨ane, and D. K¨odderitzsch For the graphical user interface of Mathematica we designed two additional palettes. The palette SymmetryElements contains all the symmetry elements of the 32 point groups. The palette PointGroup contains all the operations defined in the package. The inspection of the multiplication table (MultiplicationTable[c4v]) of the group shows, that the table is not symmetric, indicating that the group multiplication is not commutative in this group. Therefore the group is not abelian (AbelianQ[c4vs] returns False.). All elements of the group can be also generated by successive multiplication of only two elements of the group, the so called generators. For our example we find the generators C4z,IC2x. (Generators[group_]) 2 Any set of elements of a group which obeys all the group postulates is called a subgroup of the group. Examples for subgroups are: E,C2z, C4z, C−1 4z and {E,IC2y} (SubGroupQ[group_,subgroup_] performs a test). Another structure of the group is introduced by the definition of conjugate elements. An element B of the group G is said to be conjugate to A if their exists a group element X such that B = XAX−1. A class (Classes[c4vs]) is a collection of mutually conjugate elements of a group. The classes of the group C4v are: (E), (C2z), (C4z, C−1 4z ), (IC2x,IC2y), (IC2a,IC2b). It is easy to prove that: i) E always forms a class on its own, ii) every group element in G is a member of some class of G, iii) no element can be the member of two classes of G, iv) if G is an abelian group, every element forms a class on its own. We have already seen, that the group in the example can be represented by the abstract symbols or the rotation matrices. The relation between groups is formulated more rigorously in terms of homomorphism and isomorphism. Two groups G = {A, B, C, . . . } and G = {A , B , C ,... } are called homomorphic, if i) to each element of G corresponds one and only one element of G ii) for all elements holds A → A , B → B → AB → A B . The set of elements which are mapped to E ∈ G in a homomorphism is called kernel of the homomorphism NK. In contrast to the homomorphism the mapping in the isomorphism is one-to-one. Whereas the mapping of the symbols on the rotation matrices constitutes an isomorphism, the following mapping E,C2z, C4z, C−1 4z ⇒ +1 and IC2x,IC2y,IC2a,IC2b ⇒ −1 constitutes a homomorphism. 6.4 Representation Theory 6.4.1 Matrix Representations of Groups Matrix representations of symmetry groups are the essential tools to investigate symmetry properties of solutions of field equations as Schr¨odinger’s equation or Maxwell’s equations. 2 The formal arguments of the commands indicate, which type of information has to be plugged in. For more information use the help system or the example notebook.
6 Symmetry Properties of Electronic and Photonic Band Structures 107 A group of square matrices,with the matrix multiplication as relation between the elements,which is homomorphic to a group g is called a matrix representation of g.Each element Ag corresponds to a matrix I(A): T(A)T(B)=I(C)V A,B,CEg T(E)=E (identity matrix) (6.2) T(A-1)=T(A)-1 It is obvious that an I-dimensional representation consisting of a set of(Ix l)- matrices can be transformed in another representation by means of a similar- ity transformation I(A)=S-1I(A)S using a non-singular(I x I)-matrix S. The rotation matrices of a finite symmetry group form a three-dimensional matrix representation in that sense.The mapping of the group elements of Cav onto the numbers +1,-1 represents also,in this case a one-dimensional, matrix representation of the group. From two matrix representations T and T2 of g with the dimensions l and l2,a(l1+l2)-dimensional representation I can be built,forming the block matrices T(A (6.3) The representation I is called the direct sum of the representations T and T2:T=IT2.Representations which can be transformed into direct sums by similarity transformations are called reducible.If such a transformation is not possible,the representation is called irreducible.The rotation matrices of the group Cav are given by: 100 -100 010 T(E)= 010 (C2z) 0-10 I(Caz)= 100 001 001 001 0-10 -100 100 r(IC)= 10 0 (IC2x)= 010 -10 001 001 001 0-10 010 T(IC2a)= -100 T(IC2b)= 100 001 001 The rotation matrices are block matrices and therefore reducible into a two- dimensional and a one-dimensional matrix representation of Cav. An important quantity of a matrix representation which will not change under a similarity transformation is the trace of the rotation matrix.This trace is called the character x of the representation matrix
6 Symmetry Properties of Electronic and Photonic Band Structures 107 A group of square matrices, with the matrix multiplication as relation between the elements, which is homomorphic to a group G is called a matrix representation of G. Each element A ∈ G corresponds to a matrix Γ(A): Γ(A)Γ(B) = Γ(C) ∀ A, B, C ∈ G Γ(E) = E (identity matrix) (6.2) Γ(A−1) = Γ(A) −1 It is obvious that an l-dimensional representation consisting of a set of (l×l)- matrices can be transformed in another representation by means of a similarity transformation Γ (A) = S−1Γ(A)S using a non-singular (l × l)-matrix S. The rotation matrices of a finite symmetry group form a three-dimensional matrix representation in that sense. The mapping of the group elements of C4v onto the numbers +1, −1 represents also, in this case a one-dimensional, matrix representation of the group. From two matrix representations Γ1 and Γ2 of G with the dimensions l1 and l2,a(l1 + l2)-dimensional representation Γ can be built, forming the block matrices Γ(A) = Γ1(A) 0 0 Γ2(A) . (6.3) The representation Γ is called the direct sum of the representations Γ1 and Γ2 : Γ = Γ1⊕Γ2. Representations which can be transformed into direct sums by similarity transformations are called reducible. If such a transformation is not possible, the representation is called irreducible. The rotation matrices of the group C4v are given by: Γ(E) = 100 010 001 Γ(C2z) = −100 0 −1 0 0 01 Γ(C4z) = 0 10 −100 0 01 Γ(IC−1 4z ) = 0 −1 0 100 001 Γ(IC2x) = −100 0 10 0 01 Γ(IC2y) = 100 0 −1 0 001 Γ(IC2a) = 0 −1 0 −100 0 01 Γ(IC2b) = 010 100 001 The rotation matrices are block matrices and therefore reducible into a twodimensional and a one-dimensional matrix representation of C4v. An important quantity of a matrix representation which will not change under a similarity transformation is the trace of the rotation matrix. This trace is called the character χ of the representation matrix.