A microwave realization of the chiral orthogonal, unitary, and symplectic ensembles
A. Rehemanjiang M. Richter U. Kuhl H. -J. Stöckmann



Fachbereich Physik der Philipps-Universität Marburg, D-35032 Marburg, Germany

School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, UK
Université Côte d'Azur, CNRS, Institut de Physique de Nice (InPhyNi), 06108 Nice, France

Fachbereich Physik der Philipps-Universität Marburg, D-35032 Marburg, Germany
Université Côte d'Azur, CNRS, Institut de Physique de Nice (InPhyNi), 06108 Nice, France

Fachbereich Physik der Philipps-Universität Marburg, D-35032 Marburg, Germany



Random matrix theory has proven very successful in the understanding of the spectra of chaotic systems.
Depending on symmetry with respect to time reversal and the presence or absence of a spin 1/2 there are three ensembles, the Gaussian orthogonal (GOE), Gaussian unitary (GUE), and Gaussian symplectic (GSE) one.
With a further particle-antiparticle symmetry the chiral variants of these ensembles, the chiral orthogonal, unitary, and symplectic ensembles (the BDI, AIII, and CII in Cartan's notation) appear.
A microwave study of the chiral ensembles is presented using a linear chain of evanescently coupled dielectric cylindrical resonators.
In all cases the predicted repulsion behavior between positive and negative eigenvalues for energies close to zero could be verified.




05.45.Mt

A microwave realization of the chiral orthogonal, unitary, and symplectic ensembles
    H.-J. Stöckmann
    December 30, 2023
===================================================================================



Random matrix theory originally had been developed by Wigner, Dyson, Mehta <cit.> and others as a tool to describe the spectral properties of chaotic systems.
For time-reversal symmetric systems the Hamiltonian H commutes with the time-reversal operator T, HT=TH, where T^2=1 for systems without spin 1/2, and T^2=-1 in the presence of a spin 1/2 <cit.>.
The three options (T^2=1, no T, T^2=-1) give rise to the three classical random matrix ensembles, the Gaussian orthogonal (GOE), unitary (GUE) and symplectic (GSE) one, respectively.
Further there may be a chiral symmetry, i. e., an operator C anticommuting with H, HC=-CH, again with the two options C^2=1 and C^2=-1.
Such a symmetry exists, e. g., for the Dirac equation<cit.>.
All possible combinations of (T^2=1, no T, T^2=-1) and (C^2=1, no C, C^2=-1) yield a total of nine random matrix ensembles.
Together with the last remaining option (no T, no C, but CT) one finally ends up with the ten-fold way <cit.>.

For the GOE there is an abundant number of realizations, see Sec. 3.2 of Ref. stoe99, but for the GUE the number of experiments is still small <cit.>.
The GSE has been realized only recently by us in a peculiarly designed microwave network mimicking a spin 1/2 <cit.>.

For the new ensembles systematic random matrix studies are still missing though there are a lot of studies of systems showing chiral symmetry <cit.>.
In the present work a microwave study of the chiral relatives of the classical ensembles shall be presented, the chiral orthogonal (chOE), the chiral unitary (chUE), and the chiral symplectic (chSE) ensemble (the BDI, the AIII, and the CII in Cartan's notation).
We omit the `G' in the notation, since the ensembles studied by us are partly not Gaussian.

For a chiral symmetry particles and anti-particles are not really needed.
Sufficient is a system consisting of two subsystems I and II with interactions only between I and II, but no internal interactions within I or II.
The Hamiltonian for such a situation may be written as

    H=([   0   A; A^†   0 ]) ,

where the diagonal blocks belong to the two subsystems, and the off-diagonal blocks describe the interaction.
Hamiltonian (<ref>) is chiral symmetric, since it anticommutes with C=diag( 1_n,- 1_m).
The characteristic polynomial of H is given by

    χ(E)= |
                   [ E· 1_n     -A;   -A^† E· 1_m;        ]|=
        E^n-m|E^2· 1_m-A^† A| ,

where an elementary matrix transformation has been applied.
Equation (<ref>) has a number of important implications:
(i) For n>m there are μ=n-m zero eigenvalues.
They do not depend on the interaction between the two subsystems and can only be destroyed by lifting the chiral symmetry.
(ii) All other eigenvalues appear in pairs E_k and E_-k=-E_k (k=1,2,…).
(iii) For energies far from E=0 the statistical properties of the chiral ensembles approach those of the classical ones, but close to E=0 the eigenvalues E_1, E_2, … feel the proximity of their partners E_-1, E_-2, …, resulting in an oscillatory modulation in the density of states ρ (E) and a possible repulsion of the non-zero eigenvalue pairs close to E=0 <cit.>,

    ρ (E) ∼ |E|^α+μβ ,

where β is the universality index known already from the classical ensembles (β=1, 2, 4 for the GOE, the GUE, and the GSE, respectively), and α=β -1 for their chiral relatives <cit.>.

The main building blocks for our experimental realization of the chiral ensembles are dielectric cylindrical resonators (h=5 mm, r=3.8 mm) with an index of reflection of n=6, see inset in Fig. <ref>(bottom).
They are placed between aluminum bottom and top plates separated by a distance of H=11 mm.

There are two types of resonance modes, the transverse magnetic (TM) and transverse electric (TE) one.
In the experiment we used the lowest TE mode with magnetic field parallel and electric field perpendicular to the cylinder axis.
As the resonators were not perfectly equal, the mean eigenfrequency ν_0=13.6 GHz spread by about 3 MHz corresponding to 20 percent of the line width.
In the experiment ν_0 corresponds to the energy zero and ν-ν_0 to the energy.
As ν_0 is below the cutoff frequency ν_cutoff=c/(2H)=9.369 GHz, where c is the velocity of light, the resonators are evanescently coupled allowing to vary the coupling strength via the distance.
For further experimental details it is referred to <cit.>.
Reflections S_11, S_22 and transmissions S_12, S_21 have been measured in the range 6.7 to 7.0 GHz using a vector network analyzer (Agilent 8720ES).



For the realization of the chOE the set-up shown in
Fig. <ref>(bottom) had been used.
Up to five resonators are placed in a row.
The resonances are excited by a bent antenna placed close to the leftmost resonator.
Figure <ref>(top) shows typical reflection spectra for
N = 2, 3, 4 and 5 resonators obtained by removing one resonator after the other
from the right.
The chiral symmetry is clearly evident from the spectra.
For N>5 it became more and more difficult to extract all resonances.
Furthermore localization effects eventually became important.

The situation can be mapped onto a Hamiltonian with the eigenfrequencies of the resonators in the diagonal, the coupling constants in the secondary diagonals, and zeros everywhere else.
Ordering row and columns with odds sites first and even ones secondly, the Hamiltonian takes the structure (<ref>).
For N=2 the eigenfrequencies are given by
ν_±=ν_0 ±√(|a|^2+Δ^2/4),
where Δ is the difference of the eigenfrequencies of the two unperturbed resonators, and a the coupling constant due to the evanescent coupling.
Figure <ref>(bottom) shows the eigenfrequencies for a two-resonators system in
dependence of their distance.
Such two-resonator measurements have been used for a calibration of the coupling constants in terms of the
distance.

For the chOE some hundred realizations have been taken, with Gaussian distributed coupling constants,


    p_β(a)∼ a^β -1e^-a^2/2σ^2 ,


where β=1, and a cutoff at the maximal available coupling
constant a_max= 102.56 MHz.
With σ= 0.354 a_max this means a truncation of 0.5 percent of the area in the tail of the distribution.
For N=2 and 3 the A matrix is completely filled, i.e., the resulting ensembles are thus Gaussian.
For increasing N A contains more and more zero matrix elements, hence these ensembles are no longer Gaussian.
While this modifies the details of the modulation of the ensemble averaged density of states close to E=0, it does not change the repulsion behavior (<ref>).
Figure <ref>, left panel, shows the resulting ensemble averaged density of states ρ(ν-ν_0) for linear chains of lengths N=2, 3, 4, and 5.





For N≤ 5 the zeros of the characteristic polynomial (<ref>) can be obtained as roots of quadratic equations.
It remains a Gaussian average over the coupling constants to calculate the ensemble averaged density of states <cit.>.
The dashed lines in Fig. <ref>, left panel, show the analytical expressions for ρ(ν-ν_0).
For odd N a delta peak at ν=ν_0 is predicted, and a linear repulsion of the eigenfrequencies, see Eq. (<ref>).
For even N there should be no repulsion.

All these features are found in the experiment.
There is only one small mismatch resulting from the fact that eigenfrequencies of the resonators are not identical, but differ by some MHz.
Furthermore, the eigenfrequency of the leftmost resonator is detuned by the nearby antenna.
This results in a hole in the distribution at ν=ν_0 for N=2, and to a smaller extent also for N=4.

Each non-zero element in the diagonal block of the Hamiltonian
(<ref>) destroys the chiral symmetry.
This may be due to next-nearest neighbor contributions, or by different
site energies as in the present situation.
This shortcoming is unavoidable: 
Even if the eigenfrequencies of the isolated resonators would be identical they are detuned by the presence of the other ones.
A manifestation of this detuning is the shift of the center of gravity of the two-resonator spectrum observed for small distances, see Fig. <ref>(bottom).

Perturbation of the chiral symmetry results in two imperfections: 
(i) Shift of the center of gravity of the spectrum, 
(ii) A left-right asymmetry between the “positron” and the “electron” part of the spectrum.
For the linear chain with an odd number of elements this means in particular that the zero energy peak is smeared out.
In addition the repulsion behavior is distorted as is evident, e. g., from Fig. <ref>, left panel, for N=3.

Therefore we studied in addition another quantity, the two-point correlation function ρ_+1,-1(Δν) giving the probability density to find a frequency distance Δν between the states ν_1 and ν_-1.
Because of the chiral symmetry ρ_+1,-1(Δν) is, up to a factor of two, identical with the ensemble averaged density ρ_1(ν-ν_0) of state 1.
Figure <ref>, right panel, shows the results, together with the random matrix expectation.
To accentuate the repulsion behavior, the integrated pair correlation function I_+1,-1(Δν)=∫_0^Δν dν' ρ_+1,-1(ν') is shown in the inset in a log-log plot.
Because of the integration the repulsion exponents increase by one, i. e., repulsion exponents should be 1 and 2 for odd and even N, respectively.
For all cases the expected behavior is verified.
For odd N the agreement between theory and experiment is nearly perfect.
Only for the lowest shown Δν values small deviations appear.
Here the number of realizations entering is below ten and the results are no longer statistically trustworthy.
For even N the deviations are somewhat larger, again a manifestation of the fact that the eigenfrequencies of the resonators are not identical.



Since the zeros of the characteristic polynomial (<ref>) depend only on the moduli of the coupling constants but not their phases, a convenient way to realize chUE and chSE statistics would be to reuse the same experimental setup and to only take a modified function for the coupling constants, again Eq. (<ref>), but now with β=2, 4 for the chUE and the chSE, respectively.
We preferred a more ambitious approach and studied systems really having a broken time reversal (β=2) or symplectic symmetry (β=4).

In our realization of the chUE the break of time reversal symmetry was achieved by means of circulators, which had been used already for this purpose previously <cit.>.
A circulator is a microwave device with three ports, where waves entering via ports 1, 2, 3 exit through ports 2, 3, 1, respectively.
The inset of Fig. <ref>(bottom) shows the set-up.
Now the resonators are at fixed distances of 85 mm, too large for a direct coupling which is achieved instead by means of circulators.
The resonances had been excited by loop antennas placed on top of the resonators requiring a height H=16 mm.

Figure <ref>(top) shows typical spectra for N=2, 3, and 4 resonators.
Unfortunately, the circulators introduce a considerable broadening making an analysis difficult.
Therefore, the harmonic inversion technique has been used, allowing for an analysis of the spectra also for overlapping resonances <cit.>.
The deviations between the original and the reconstructed spectra are hardly seen, demonstrating the reliability of the technique.
For N>4 a reliable analysis of the spectra was no longer possible.

Since the circulators introduce directionality, the coupling constants are complex.
A Gaussian distribution of the real and imaginary parts results in a distribution p_2(|a|) (see Eq. <ref>) for the modulus, and a uniform distribution for the phase.
The phase could be varied by attaching different terminators to the third port of the circulators.
Again a two-resonator measurement was used to relate the distance between circulators and resonators to the coupling.
The red symbols in Fig. <ref>(bottom) show the result for an open-end terminator attached.
The shift of the center of gravity of the two resonances is large, comparable to the splitting.
For short-end terminators (blue dots) the effect is even more dramatic.
Therefore we refrained from varying the phases by attaching different terminators to the circulators and used only the open-end terminators.
But even here for small circulator distances the shift of the center of gravity turned out to be too large.



Therefore, only distances d≥ 2 mm were used corresponding to a_max= 4.47 MHz, with σ= 0.288 a_max.
The results are shown in the left column of Fig. <ref>.
Though there are significant deviations in detail, the general features of the theoretical predictions are correctly reproduced.
For N=3 and 4 the center of gravity of the spectra had been adjusted individually by less than 1 MHz.
This shift is caused by a detuning of the eigenfrequencies of the inner resonators due to the presence of the second circulators, an effect not accounted for by the calibration.
All these deficiencies drop out for the two-point correlation function ρ_+1,-1(Δν).
Therefore, for this quantity the whole d range depicted in Fig. <ref> could be used.
The results are presented in the right column of Fig. <ref>.
The solid lines correspond to the theoretical expectations for the unperturbed systems.
The inset, showing the integrated two-point correlation in a log-log plot, illustrates that for the chUE, too, the expected repulsion exponents are reproduced.



Recently we realized the Gaussian symplectic ensemble in a microwave graph mimicking a spin 1/2 <cit.>.
The main ingredients had been two subgraphs, complex conjugates of each other, coupled by a pair of bonds with a phase shift of π in one of the bonds and no phase shift in the other one.
These ideas may be taken over one-to-one to the present situation.
The minimal configuration needs four resonators, 1, 1̅, 2, 2̅, see Fig. <ref>(left).
Resonators 1 and 2, and resonators 1̅ and 2̅ are coupled by circulators with opposite sense of propagation thus making the two subsystems complex conjugates of each other.
Resonators 1 and 2̅, and resonators 1̅ and 2 are coupled by cables with a length difference of 1.762 cm corresponding to a phase difference of π at ν=ν_0.

Though four resonators are involved, the configuration corresponds to N=2, where the resonators with and without bars mimic the two spin components.
The corresponding characteristic polynomial is given by

    χ(E)=|E· 1 -H|= (E^2-|a|^2-|b|^2)^2 .

There are hence doubly degenerate eigenvalues at

    E_+1,-1=±√(|a|^2+|b|^2) .

The four-resonator system thus shows both the chiral symmetry, with E_+1 and E_-1=-E_+1 coming in pairs, and the symplectic symmetry, with the characteristic Kramers doublet structure of the spectrum.

Again the harmonic inversion technique was used to analyze the spectra.
Figure <ref>(right) shows the distance distribution of the two Kramers doublets, where the distribution function p_2(|a|) (see Eq. <ref>) was used both for |a| and |b| (with a_max = 7.13 MHz and σ = 0.322 a_max).
Hence, Eq. (<ref>) implies a cubic repulsion of E_+1 and E_-1, in accordance with Eq. (<ref>).
This repulsion is perfectly reproduced in the experiment.

An extension to larger N values seem hardly feasible.
Each additional pair of resonators would mean four more bonds and a corresponding increase of absorption.
Therefore we have to be content with this demonstration for N=2.


H. Schomerus, Lancaster University, is thanked for illuminating discussions.
This work was funded by the Deutsche Forschungsgemeinschaft via the individual grants STO 157/16-2 and KU 1525/3-1 as well as the European Commission through the H2020 programm by the Open Future Emerging Technology “NEMF21” Project (664828).









20
	fxundefined [1]
		ifx#1
	
	fnum [1]
		#1firstoftwo
		secondoftwo
		
	fx [1]
		#1firstoftwo
		secondoftwo
		
	
	
	
	
	
	noop [0]secondoftwo
	
	ref[1]@startlink#1@href
	href[1]#1@endlink
	anitize@url [0]`
12`$12`&12`#12`1̂2`_12`%12
	startlink[1]
	endlink[0]
	
	rl [1]href #1
	
	
	
	
	
	
	
	
	
	
	
	
	@bib@innerbibempty
	
	[Dyson(1962)]dys62c
	author author F. J. Dyson, title title The threefold
				way. Algebraic structure of symmetry groups and ensembles in quantum
				mechanics, 10.1063/1.1703863 journal
		journal J. Math. Phys. volume 3,		pages 1199 (year 1962)NoStop
	[Mehta(1991)]meh91
	author author M. L. Mehta, 10.1016/C2009-0-22297-5 title
			Random Matrices (Revised and Enlarged Second Edition) (publisher Academic Press, address San Diego, year 1991)NoStop
	[Haake(2010)]haa10
	author author F. Haake, 10.1007/978-3-642-05428-0 title Quantum Signatures of Chaos. 3rd edition (publisher
	Springer, address Berlin, year
	2010)NoStop
	[Verbaarschot and Zahed(1993)]ver93
	author author J. J. M.			Verbaarschot and author I. Zahed, title title
			Spectral density of the QCD Dirac operator near zero virtuality,	10.1103/PhysRevLett.70.3852 journal journal Phys. Rev. Lett. volume 70, pages 3852 (year 1993)NoStop
	[Zirnbauer(1996)]zir96
	author author M. R. Zirnbauer, title title Riemannian
				symmetric superspaces and their origin in random-matrix theory, 10.1063/1.531675 journal journal J.
				Math. Phys. volume 37, pages 4986
		(year 1996)NoStop
	[Altland and Zirnbauer(1997)]alt97c
	author author A. Altland and author M. R. Zirnbauer, title title Nonstandard
				symmetry classes in mesoscopic normal-superconducting hybrid structures,	10.1103/PhysRevB.55.1142 journal journal Phys. Rev. B volume 55, pages 1142 (year 1997)NoStop
	[Stöckmann(1999)]stoe99
	author author H.-J. Stöckmann, 10.2277/0521592844 title Quantum Chaos - An Introduction (publisher
	University Press, address Cambridge, year
	1999)NoStop
	[So et al.(1995)So,
		Anlage, Ott, and Oerter]so95
	author author P. So, author S. M. Anlage,
		author E. Ott,  and author R. N. Oerter, title
	title Wave chaos experiments with and without time
				reversal symmetry: GUE and GOE statistics, 10.1103/PhysRevLett.74.2662 journal journal
			Phys. Rev. Lett. volume 74, pages
		2662 (year 1995)NoStop
	[Stoffregen et al.(1995)Stoffregen, Stein, Stöckmann,
		Kuś, and Haake]sto95b
	author author U. Stoffregen, author J. Stein,
		author H.-J. Stöckmann,
		author M. Kuś,  and author F. Haake, title
	title Microwave billiards with broken time reversal
				symmetry, 10.1103/PhysRevLett.74.2666 journal journal Phys. Rev. Lett. volume 74, pages 2666 (year
		1995)NoStop
	[Ławniczak et al.(2010)Ławniczak, Bauch, Hul, and		Sirko]law10
	author author M. Ławniczak, author S. Bauch,
		author O. Hul,  and author L. Sirko, title
	title Experimental investigation of the enhancement
				factor for microwave irregular networks with preserved and broken time
				reversal symmetry in the presence of absorption, 10.1103/PhysRevE.81.046204 journal journal Phys.
				Rev. E volume 81, pages 046204
		(year 2010)NoStop
	[Rehemanjiang et al.(2016)Rehemanjiang, Allgaier, Joyner,
		Müller, Sieber, Kuhl,		and Stöckmann]reh16
	author author A. Rehemanjiang, author M. Allgaier, author C. H. Joyner, author S. Müller,
		author M. Sieber, author U. Kuhl,  and author
		H.-J. Stöckmann, title
	title Microwave realization of the Gaussian
				symplectic ensemble, 10.1103/PhysRevLett.117.064101
	journal journal Phys. Rev. Lett. volume 117, pages 064101 (year
		2016)NoStop
	[Beenakker(2015)]bee15
	author author C. W. J.			Beenakker, title title
			Random-matrix theory of Majorana fermions and topological
				superconductors, 10.1103/RevModPhys.87.1037 journal journal Rev. Mod. Phys. volume
			87, pages 1037 (year 2015)NoStop
	[Ivanov(2002)]iva02
	author author D. A. Ivanov, title title The
				supersymmetric technique for random-matrix ensembles with zero
				eigenvalues, 10.1063/1.1423765 journal
		journal J. Math. Phys. volume 43,		pages 126 (year 2002)NoStop
	[Barkhofen et al.(2013)Barkhofen, Bellec, Kuhl, and		Mortessagne]bar13a
	author author S. Barkhofen, author M. Bellec,
		author U. Kuhl,  and author F. Mortessagne, title title Disordered graphene and boron nitride in
				a microwave tight-binding analog, 10.1103/PhysRevB.87.035101 journal journal Phys.
				Rev. B volume 87, pages 035101
		(year 2013)NoStop
	[Bellec et al.(2013)Bellec,
		Kuhl, Montambaux, and Mortessagne]bel13b
	author author M. Bellec, author U. Kuhl,
		author G. Montambaux,  and		author F. Mortessagne,	title title Tight-binding couplings in
				microwave artificial graphene, 10.1103/PhysRevB.88.115437 journal journal Phys.
				Rev. B volume 88, pages 115437
		(year 2013)NoStop
	[Richter et al.(2019)Richter, Rehemanjiang, Kuhl, and		Stöckmann]arxRic19
	author author M. Richter, author A. Rehemanjiang, author U. Kuhl,
		 and author H.-J. Stöckmann, @noop  title Microwave
				studies of the three chiral ensembles in chains of coupled dielectric
				resonators,  (year 2019), note to be
		publishedNoStop
	[Wall and Neuhauser(1995)]wal95
	author author M. R. Wall and author D. Neuhauser, title title Extraction,
				through filter-diagonalization, of general quantum eigenvalues or classical
				normal mode frequencies from a small number of residues or a short-time
				segment of a signal. I. Theory and application to a quantum-dynamics
				model, 10.1063/1.468999 journal journal J. Chem. Phys. volume 102, pages 8011 (year 1995)NoStop
	[Mandelshtam and Taylor(1997)]man97a
	author author V. A. Mandelshtam and author H. S. Taylor, title title Harmonic
				inversion of time signals and its applications, 10.1063/1.475324 journal journal J. Chem.
				Phys. volume 107, pages 6756
		(year 1997)NoStop
	[Main(1999)]mai99
	author author J. Main, title title Use of harmonic
				inversion techniques in semiclassical quantization and analysis of quantum
				spectra, 10.1016/S0370-1573(98)00131-8 journal journal Phys. Rep. volume
			316, pages 233 (year 1999)NoStop
	[Kuhl et al.(2008)Kuhl,
		Höhmann, Main, and Stöckmann]kuh08b
	author author U. Kuhl, author R. Höhmann,
		author J. Main,  and author H.-J. Stöckmann, title title Resonance widths in open microwave
				cavities studied by harmonic inversion, 10.1103/PhysRevLett.100.254101 journal journal
			Phys. Rev. Lett. volume 100, pages
		254101 (year 2008)NoStop



