Optical properties of honeycomb photonic structures Artem D. Sinelnik1, Mikhail V. Rybin2;3,∗ Stanislav Y. Lukashenko1, Mikhail F. Limonov2;3, and Kirill B. Samusev2;3 1Department of Nanophotonics and Metamaterials, ITMO University, St. Petersburg 197101, Russia 2Ioffe Institute, St. Petersburg 194021, Russia and 3Department of Dielectric and Semiconductor Photonics, ITMO University, St.Petersburg 197101, Russia (Dated: September 5, 2018) We study, theoretically and experimentally, optical properties of different types of honeycomb photonic struc- tures, known also as ‘photonic graphene’. First, we employ the two-photon polymerization method to fabricate the honeycomb structures. In experiment, we observe a strong diffraction from a finite number of elements, thus providing a unique tool to define the exact number of scattering elements in the structure by a naked eye. Then, we study theoretically the transmission spectra of both honeycomb single layer and 2D structures of parallel dielectric circular rods. When the dielectric constant of the rod materials " is increasing, we reveal that a two- dimensional photonic graphene structure transforms into a metamaterial when the lowest TE01 Mie gap opens up below the lowest Bragg bandgap. We also observe two Dirac points in the band structure of 2D photonic graphene at the K point of the Brillouin zone and demonstrate a manifestation of the Dirac lensing for the TM polarization. The performance of the Dirac lens is that the 2D photonic graphene layer converts a wave from point source into a beam with flat phase surfaces at the Dirac frequency for the TM polarization. I. INTRODUCTION is one of the physical phenomena that can be observed in the photonic as well as in the electronic realization of graphene. The study of honeycomb structures being known for thou- Photonic properties of a 2D photonic crystal with honey- sands of years is increased in the last few years due to the comb lattice fabricated in silicon using electron-beam lithog- unique physical properties of graphene [1]. Graphene, a two- raphy were calculated numerically and measured experimen- dimensional honeycomb lattice of carbon atoms, has attracted tally [12]. Lattice parameters of the fabricated structures made enormous attention since its first discovery by isolation from it possible to achieve and confirm experimentally the com- bulk graphite using adhesive tape. Graphene and graphene- plete photonic band gap in vicinity of the optical communi- µ based materials demonstrated a very fast development of both cations wavelengths around 1.5 m. Note that the photonic fundamental and practical aspects in optics and electronics band structure of a 2D triangular lattice of air rods in a di- " [2]. A new research hotspot becomes a Dirac point that was electric medium with high can be found in textbooks [13]. r=a = 0:48 first investigated in the electronic energy band structure of For the particular radius and dielectric constant " = 13 graphene [3]. Importantly, this concept penetrated into the , TM and TE band gaps overlap, resulting in an 18.6% field of optics where the so-called photonic graphene, a two- complete photonic band gap. Also the effect of a Dirac point dimensional photonic crystal structure that is analogous to on the transmission spectra was calculated for a 2D photonic graphene, has been studied theoretically and extensively [4– crystal of air rods formed a triangular lattice [5]. Recently it 8] is reported that a photonic crystal with honeycomb-like lattice supports a Dirac mode with a different algebraic-decay field The conical diffraction and the dynamics of optical waves profile [14]. in photonic graphene was studied experimentally [4, 9]. In Nevertheless, the optical studies of photonic graphene are a particular, it was demonstrated that an incident narrow light long way still from completeness. Neither the photonic band beam at 488 nm wavelength, with momentum at the vicin- structure and transmission spectra of 2D crystal with honey- ity of a diabolic point, diffracts in the honeycomb lattice in comb lattice at high dielectric constant " nor optical diffrac- a characteristic conical form, obtaining the shape of a ring tion from photonic graphene were studied theoretically or ex- whose thickness does not broaden, whereas its radius grows perimentally. To improve the situation, in this paper, we cal- linearly with distance [4]. An unconventional edge state at the culate the transmission spectra of photonic graphene, inves- zigzag edge of an optically induced honeycomb lattice resid- tigate optical Laue diffraction both theoretically and exper- ing on the bearded edge was demonstrated experimentally [8]. imentally and, finally, perform a comprehensive studies of Very recently, a novel kind of dispersive edge state was found photonic properties of 2D honeycomb lattice of parallel di- in a photonic honeycomb lattice experimentally that emerges electric circular rods. arXiv:1707.09780v1 [cond-mat.mes-hall] 31 Jul 2017 not only in zigzag and bearded terminations, but also in arm- chair edges [10]. At a microwave frequency, the measured transmission through a structure of photonic graphene enters a pseudodiffusive regime where the transmission scales in- II. FABRICATION OF PHOTONICS LATTICES versely with the thickness of the crystal [11]. The observation of a pseudodiffusive transport regime in a photonic graphene Synthesis of graphene-based nanomaterials includes me- chanical exfoliation from graphite, epitaxial growth from ther- mal deposition of SiC, chemical vapor deposition, solution- based exfoliation from graphite oxide and graphite, bottom-up ∗Electronic address: [email protected] organic synthesis methods [2]. In this study, to fabricate pho- 2 n = 0 n = -1 n = +1 (a) (b) (c) (d) FIG. 1: (a) Schematic of the zero-order (n = 0) and first-order (n = ±1) Laue diffraction from a horizontally oriented chain of scatterers. (b) Schematic of the zero- and first-order Laue diffraction from a hexagonal structure. (c) Fragment of the diffraction pattern of honeycomb structure with its characteristic features discussed in text. (d) Photograph of experimental setup. tonic counterpart of graphene, a honeycomb structure, we em- scatterers varied from tens to tens of thousands. ploy two-photon polymerization method which is also called direct laser writing [15, 16]. This method is based on the non- linear two-photon polymerization of a photosensitive material III. LAUE DIFFRACTION FROM HONEYCOMB in the focus of a femtosecond laser beam and makes it possible STRUCTURE to form a photonic structure with a transverse resolution below 100 nm [17]. In addition, we fabricate 2D photonic structures The samples produced by the two-photon polymerization of hexagonal symmetry (constituted by dense-packed trian- method are characterized by a low dielectric permittivity con- gles). As this was done in our preceding studies, we use the trast. In this case, the optical diffraction is described in the installation and software package from Laser Zentrum Han- Born approximation [19] when the diffraction intensity is de- nover (Germany) [18]. The structures are fabricated from a termined by a product of the squared structural factor S(q), hybrid organic-inorganic material based on zirconium propox- scattering form factor F (q) and polarization factor [20]. Our ide and an Irgacure 369 photo-initiator (Ciba Specialty Chem- calculation in which only the structural factor is taken into ac- icals Inc., Switzerland). The polymerization is performed count allows us to adequately describe the experimental data, with a train of femtosecond pulses (wavelength was 780 nm) and, therefore, we disregard the contributions from the form at a repetition frequency of 80 MHz (12.5 ns between con- factor and the polarization factor. In the 2D case, the diffrac- secutive pulses) from a 50-fs TiF-100F laser (Avesta-Project, tion is determined by S(q) = P exp(iqr ), where r = Russia). The sample is mounted on a two-coordinate mo- i i i a n + a n , q ≡ k − k is the scattering vector, k and k torized linear air-bearing translator (Aerotech Inc., USA), by 1 1 2 2 i s i s are the wave vectors of the incident and scattering waves. The which the laser focus is moved along a 2D scanning path structural factor has maxima in the scattering directions k de- in the XY-plane. The laser radiation is focused in the pho- s termined by the system of Laue equations: (k − k ) · a = toresist volume through the glass substrate with a 100x oil- s i 1;2 2πn . For a 1D linear chain of point-like scatterers r = immersion microscope objective with numerical aperture NA 1;2 n na jS(q)j2 = sin2(Nqa=2)=(qa=2) = 1.4 (Carl Zeiss MicroImaging GmbH, Germany). The cor- 1, one has [21]. When sin(qa=2) ! 0 respondence of the resulting materials to the designed struc- the denominator becomes zero, , the structural jS(q)j2 tures is confirmed by scanning electron microscopy (SEM). factor has strong maxima named “principal” maxima In the structures, the lattice parameters are within the range in what follows. Figure 1c depicts the diffraction pattern of a 0:5µm a 2:0µm in different samples and the number of calculated honeycomb structure, which shows seven principal 6 6 diffraction peaks, one of which is marked by a red square. 3 2 p When the numerator goes to zero, sin (Nqa=2) = 0, the where ζ = jksj sin ' 3a=4, and ' is the angle between ks 2 function jS(q)j becomes zero, and this means that it will and ki. have N − 1 zeroes in the interval between the principal max- The triangular lattice can be represented as a superposition ima and, accordingly, N − 2 more maxima, named “addi- of three systems of linear chains, r1 = a1n1, r2 = a2n2, and tional” in what follows. These additional maxima are marked r3 = a3n3, with the vectors a1, a2, and a3 turned relative to in Fig. 1c by a red rectangle.