Search

We study the nonlinear optical properties of heterojunctions made of graphene nanoribbons (GNRs) consisting of two segments with either the same or different topological properties. By utilizing a quantum mechanical approach that incorporates distant-neighbor interactions, we demonstrate that the presence of topological interface states significantly enhances the second- and third-order nonlinear optical response of GNR heterojunctions that are created by merging two topologically inequivalent GNRs. Specifically, GNR heterojunctions with topological interface states display third-order harmonic hyperpolarizabilities that are more than two orders of magnitude larger than those of their similarly sized counterparts without topological interface states, whereas the secondorder harmonic hyperpolarizabilities exhibit a more than ten-fold contrast between heterojunctions with and without topological interface states. Additionally, we find that the topological state at the interface between two topologically distinct GNRs can induce a noticeable red-shift of the quantum plasmon frequency of the heterojunctions. Our results reveal a general and profound connection between the existence of topological states and an enhanced nonlinear optical response of graphene nanostructures and possible other photonic systems., Comment: 16 pages,5 figures

We demonstrate that the influence of Kerr effect on valley-Hall topological transport in graphene metasurfaces can be used to implement an all-optical switch. In particular, by taking advantage of the large Kerr coefficient of graphene, the index of refraction of a topologically-protected graphene metasurface can be tuned via a pump beam, which results in an optically controllable frequency shift of the photonic bands of the metasurface. This spectral variation can in turn be readily employed to control and switch the propagation of an optical signal in certain waveguide modes of the graphene metasurface. Importantly, our theoretical and computational analysis reveals that the threshold pump power needed to optically switch ON/OFF the signal is strongly dependent on the group velocity of the pump mode, especially when the device is operated in the slow-light regime. This study could open up new routes towards active photonic nanodevices whose underlying functionality stems from their topological characteristics., Comment: 7 pages, 6 figures

We derive the $\mathcal{T}$-matrix formalism tailored for the numerical analysis of second-harmonic (SH) generation from arbitrarily-shaped particles made of centrosymmetric optical materials. First, the transfer matrix of a single particle is computed \textit{via} the extended boundary condition method, in which the electromagnetic fields both at fundamental frequency and SH are expanded in vector spherical wave functions and the integral formulation is satisfied away from the surface of the scatterer. We allow for the accurate physical description of the SH sources by taking into account both local surface and nonlocal bulk polarization contributions to the nonlinear polarization density source responsible for the generation of the SH signal in a particle. This single-particle formalism is then extended to arbitrary distributions of particles by incorporating into the formalism linear and nonlinear electromagnetic wave scattering from the particles in the cluster. Importantly from a practical point of view, our method can be applied to particles of arbitrary shape made of optical materials characterized by general frequency-dispersion relations, so that it can describe the linear and nonlinear optical response of clusters of metallic, semiconductor, or polaritonic particles, as well as mixtures of such particles. The approach proposed here is faster and more memory-efficient than well-established numerical techniques, especially in the analysis of spheroidal particles, due to the favourable symmetries of spherical wave basis functions used in the wave scattering analysis., Comment: 19 pages, 12 figures

Nonlinear topological photonics has been attracting increasing research interest, as it provides an exciting photonic platform that combines the advantages of active all-optical control offered by nonlinear optics with the unique features of topological photonic systems, such as topologically-protected defect-immune light propagation. In this paper, we demonstrate that topological interface modes and trivial edge modes of a specially designed graphene metasurface can be coupled in a tunable and optically controllable manner, thus providing an efficient approach to transfer optical power to topologically protected states. This is achieved in a pump-signal configuration, in which an optical pump propagating in a bulk mode of the metasurface is employed to tune the band structure of the photonic system and, consequently, the coupling coefficient and wave-vector mismatch between edge and topological interface modes. This tunable coupling mechanism is particularly efficient due to the large Kerr coefficient of graphene. Importantly, we demonstrate that the required pump power can be significantly reduced if the optical device is operated in the slow-light regime. We perform our analysis using both \textit{ab initio} full-wave simulations and a coupled-mode theory that captures the main physics of this active coupler and observe a good agreement between the two approaches. This work may lead to the design of active topological photonic devices with new or improved functionality., Comment: 12 pages, 8 figures

In this article, we present a $T$-matrix method for numerical computation of second-harmonic generation from clusters of arbitrarily distributed spherical particles made of centrosymmetric optical materials. The electromagnetic fields at the fundamental and second-harmonic (SH) frequencies are expanded in series of vector spherical wave functions, and the single sphere $T$-matrix entries are computed by imposing field boundary conditions at the surface of the particles. Different from previous approaches, we compute the SH fields by taking into account both local surface and nonlocal bulk polarization sources, which allows one to accurately describe the generation of SH in arbitrary clusters of spherical particles. Our numerical method can be used to efficiently analyze clusters of spherical particles made of various optical materials, including metallic, dielectric, semiconductor, and polaritonic materials., Comment: 14 pages, 11 figures

Topological photonics has revolutionized our understanding of light propagation, but most of current studies are focused on designing a static photonic structure. Developing a dynamic photonic topological platform to switch multiple topological functionalities at ultrafast speed is still a great challenge. Here we demonstrate an ultrafast reprogrammable plasmonic topological insulator, where the topological propagation route can be dynamically steered at nanosecond-level switching time, namely more than 10^7 times faster than the current state-of-the-art. This orders-of-magnitude improvement is achieved by using ultrafast electronic switches in an innovative way to implement the programmability. Due to the flexible programmability, many existing photonic topological functionalities can be integrated into this agile topological platform. Our work brings the current studies of photonic topological insulators to a digital and intelligent era, which could boost the development of intelligent and ultrafast photoelectric devices with built-in topological protection.

A novel polarized reflectarray is designed, fabricated, and experimentally characterized to show its flexibility and efficiency to control wave generation and focusing of orbital angular momentum (OAM) vortices with desirable OAM modes in the microwave frequency regime. In order to rigorously study the generation and focusing of OAM, a versatile analytical theory is proposed to theoretically study the compensation phase of reflectarray. Two prototypes of microwave reflectarrays are fabricated and experimentally characterized at 12 GHz, one for generation and one for focusing of OAM-carrying beams. Compared with the OAM-generating reflectarray, the reflectarray for focusing OAM vortex can significantly reduce the beam diameter, and this can further improve the transmission efficiency of the OAM vortex beams. We also show that the numerical and experimental results agree very well. The proposed design method and reflectarrays may spur the development of new efficient approaches to generate and focus OAM vortex waves for applications to microwave wireless communications., Comment: 10 pages, 10 figures

Nonlinear topological photonics, which explores topics common to the fields of topological phases and nonlinear optics, is expected to open up a new paradigm in topological photonics. Here, we demonstrate second-harmonic generation (SHG) via nonlinear interaction of double topological valley-Hall kink modes in all-dielectric photonic crystals (PhCs). We first show that two topological frequency bandgaps can be created around a pair of frequencies, $\omega_0$ and $2\omega_0$, by gapping out the corresponding Dirac points in two-dimensional honeycomb PhCs. Valley-Hall kink modes along a kink-type domain wall interface between two PhCs placed together in a mirror-symmetric manner are generated within the two frequency bandgaps. Importantly, through full-wave simulations and mode dispersion analysis, we demonstrate that tunable, bi-directional phase-matched SHG via nonlinear interaction of the valley-Hall kink modes inside the two bandgaps can be achieved. In particular, by using Stokes parameters associated to the magnetic part of the valley-Hall kink modes, we introduce a new concept, SHG directional dichroism, which is employed to characterize optical probes for sensing chiral molecules. Our work opens up new avenues towards topologically protected nonlinear frequency mixing and active photonic devices implemented in all-dielectric material platforms., Comment: 6 pages, 5 figures

Topologically protected plasmonic modes located inside topological bandgaps are attracting increasing attention, chiefly due to their robustness against disorder-induced backscattering. Here, we introduce a bilayer graphene metasurface that possesses plasmonic topological valley interface modes when the mirror symmetry of the metasurface is broken by horizontally shifting the lattice of holes of the top layer of the two freestanding graphene layers in opposite directions. In this configuration, light propagation along the domain-wall interface of the bilayer graphene metasurface shows unidirectional features. Moreover, we have designed a molecular sensor based on the topological properties of this metasurface using the fact that the Fermi energy of graphene varies upon chemical doping. This effect induces strong variation of the transmission of the topological guided modes, which can be employed as the underlying working principle of gas sensing devices. Our work opens up new ways of developing robust integrated plasmonic devices for molecular sensing., Comment: 4 pages, 5 figures

We demonstrate a novel scheme to dramatically enhance both the second- and third-harmonic generation in a graphene-insulator-graphene metasurface. The key underlying feature of our approach is the existence of a double-resonance phenomenon, namely the metasurface is designed to possess fundamental plasmon resonances at both the fundamental frequency and the higher harmonic. In particular, this dual resonant field enhancement at the two optical frequencies, combined with a favorable spatial overlap of the optical near-fields, lead to the increase of the generated higher harmonic by several orders of magnitude. Remarkably, we demonstrate that by tuning the Fermi energy of the graphene gratings the dual-resonance property can be locked-in over a broad spectral range of ~20 THz, and equally important, the enhanced nonlinear frequency generation process can be readily switched in the same device between the second and third harmonic. This new type of graphene metasurface could open up new avenues towards the development of novel ultra-compact and multi-frequency active photonic nanodevices., Comment: 6 pages, 5 figures

We study topologically-protected four-wave mixing (FWM) interactions in a plasmonic metasurface consisting of a periodic array of nanoholes in a graphene sheet, which exhibits a wide topological bandgap at terahertz frequencies upon the breaking of time-reversal symmetry by a static magnetic field. We demonstrate that due to the significant nonlinearity enhancement and large lifetime of graphene plasmons in specific configurations, a net gain of FWM interaction of plasmonic edge states within the topological bandgap can be achieved with pump power of less than 10 nW. In particular, we find that the effective waveguide nonlinearity coefficient is about 1.1x10^13 1/(Wm), i.e., more than ten orders of magnitude larger than that of commonly used, highly nonlinear silicon photonic nanowires. These findings could pave a new way for developing ultra-low-power-consumption, highly-integrated and robust active photonic systems at deep-subwavelength scale for applications in quantum communications and information processing., Comment: 20 pages, 6 figures

Topological photonics aims to utilize topological photonic bands and corresponding edge modes to implement robust light manipulation, which can be readily achieved in the linear regime of light-matter interaction. Importantly, unlike solid state physics, the common test bed for new ideas in topological physics, topological photonics provide an ideal platform to study wave mixing and other nonlinear interactions. These are well-known topics in classical nonlinear optics but largely unexplored in the context of topological photonics. Here, we investigate nonlinear interactions of one-way edge-modes in frequency mixing processes in topological photonic crystals. We present a detailed analysis of the band topology of two-dimensional photonic crystals with hexagonal symmetry and demonstrate that nonlinear optical processes, such as second- and third-harmonic generation can be conveniently implemented via one-way edge modes of this setup. Moreover, we demonstrate that more exotic phenomena, such as slow-light enhancement of nonlinear interactions and harmonic generation upon interaction of backward-propagating (left-handed) edge modes can also be realized. Our work opens up new avenues towards topology-protected frequency mixing processes in photonics., Comment: 16 pages, 11 figures

We demonstrate that the effective third-order nonlinear susceptibility of a graphene sheet can be enhanced by more than two orders of magnitude by patterning it into a graphene metasurface. In addition, in order to gain deeper physical insights into this phenomenon, we introduce a novel homogenization method, which is subsequently used to characterize quantitatively this nonlinearity enhancement effect by calculating the effective linear and nonlinear susceptibility of graphene metasurfaces. The accuracy of the proposed homogenization method is demonstrated by comparing its predictions with those obtained from the Kramers-Kronig relations. This work may open up new opportunities to explore novel physics pertaining to nonlinear optical interactions in graphene metasurfaces., Comment: 5 pages, 5 figures

A powerful hybrid FDTD--TDDFT method is used to study the interaction between classical plasmons of a gold bowtie nanoantenna and quantum plasmons of graphene nanoflakes (GNFs) placed in the narrow gap of the nanoantenna. Due to the hot-spot plasmon of the bowtie nanoantenna, the local-field intensity in the gap increases significantly, so that the optical response of the GNF is dramatically enhanced. To study this interaction between classical and quantum plasmons, we decompose this multiscale and multiphysics system into two computational regions, a classical and a quantum one. In the quantum region, the quantum plasmons of the GNF are studied using the TDDFT method, whereas the FDTD method is used to investigate the classical plasmons of the bowtie nanoantenna. Our analysis shows that in this hybrid system the quantum plasmon response of a molecular-scale GNF can be enhanced by more than two orders of magnitude, when the frequencies of the quantum and classical plasmons are the same. This finding can be particularly useful for applications to molecular sensors and quantum optics., Comment: 8 pages, 7 figures

In this paper, we propose a novel numerical method for modeling nanostructures containing dispersive and nonlinear two-dimensional (2D) materials, by incorporating a nonlinear generalized source (GS) into the finite-difference time-domain (FDTD) method. Starting from the expressions of nonlinear currents characterizing nonlinear processes in 2D materials, such as second- and third-harmonic generation, we prove that the nonlinear response of such nanostructures can be rigorously determined using two linear simulations. In the first simulation, one computes the linear response of the system upon its excitation by a pulsed incoming wave, whereas in the second one the system is excited by a nonlinear generalized source, which is determined by the linear near-field calculated in the first linear simulation. This new method is particularly suitable for the analysis of dispersive and nonlinear 2D materials, such as graphene and transition-metal dichalcogenides, chiefly because, unlike the case of most alternative approaches, it does not require the thickness of the 2D material. In order to investigate the accuracy of the proposed GS-FDTD method and illustrate its versatility, the linear and nonlinear response of graphene gratings have been calculated and compared to results obtained using alternative methods. Importantly, the proposed GS-FDTD can be extended to 3D bulk nonlinearities, rendering it a powerful tool for the design and analysis of more complicated nanodevices., Comment: 11 pages, 6 figures

We investigate surface modes in plasmonic Bragg fibers composed of nanostructured coaxial cylindrical metal-dielectric multilayers. We demonstrate that the existence of surface modes is determined by the sign of the spatially averaged permittivity of the plasmonic Bragg fiber, $\bar{% \varepsilon}$. Specifically, localized surface modes occur at the interface between the cylindrical core with $\bar{\varepsilon}<0$ and the outermost uniform dielectric medium, which is similar to the topologically protected plasmonic surface modes at the interface between two different one-dimensional planar metal-dielectric lattices with opposite signs of the averaged permittivity. Moreover, when increasing the number of dielectric-metal rings, the propagation constant of surface modes with different azimuthal mode numbers is approaching that of surface plasmon polaritons formed at the corresponding planar metal/dielectric interface. Robustness of such surface modes of plasmonic Bragg fibers is demonstrated too., Comment: 10 pages, 5 figures, Optics Express, to be published

We propose a distant-neighbor quantum-mechanical (DNQM) approach to study the linear and nonlinear optical properties of graphene nanoflakes (GNFs). In contrast to the widely used tight-binding description of the electronic states that considers only the nearest-neighbor coupling between the atoms, our approach is more accurate and general, as it captures the electron-core interactions between all atoms in the structure. Therefore, as we demonstrate, the DNQM approach enables the investigation of the optical coupling between two closely separated, but chemically unbound GNFs. We also find that the optical response of GNFs depends crucially on their shape, size, and symmetry properties. Specifically, increasing the size of nanoflakes is found to shift their accommodated quantum plasmon oscillations to lower frequency. Importantly, we show that by embedding a cavity into GNFs, one can change their symmetry properties, tune their optical properties, or enable otherwise forbidden second-harmonic generation processes., Comment: 22 pages and 7 figures, to appear in APL Photonics

We present a theoretical study of the characteristics of the nonlinear spin-orbital angular momentum coupling induced by second-harmonic generation in plasmonic and dielectric nanostructures made of centrosymmetric materials. In particular, the connection between the phase singularities and polarization helicities in the longitudinal components of the fundamental and second-harmonic optical fields and the scatterer symmetry properties are discussed. By in-depth comparison between the interaction of structured optical beams with plasmonic and dielectric nanostructures, we have found that all-dielectric and plasmonic nanostructures that exhibit magnetic and electric resonances have comparable second-harmonic conversion efficiency. In addition, mechanisms for second-harmonic enhancement for single and chiral clusters of scatterers are unveiled and the relationships between the content of optical angular momentum of the incident optical beams and the enhancement of nonlinear light scattering is discussed. In particular, we formulate a general angular momenta conservation law for the nonlinear spin-orbital angular momentum interaction, which includes the quasi-angular-momentum of chiral structures with different-order rotational symmetry. As a key conclusion of our study relevant to nanophotonics, we argue that all-dielectric nanostructures provide a more suitable platform to investigate experimentally the nonlinear interaction between spin and orbital angular momenta, as compared to plasmonic ones, chiefly due to their narrower resonance peaks, lower intrinsic losses, and higher sustainable optical power.

Remarkable optical and electrical properties of two-dimensional (2D) materials, such as graphene and transition-metal dichalcogenide (TMDC) monolayers, offer vast technological potential for novel and improved optoelectronic nanodevices, many of which relying on nonlinear optical effects in these 2D materials. This article introduces a highly effective numerical method for efficient and accurate description of linear and nonlinear optical effects in nanostructured 2D materials embedded in periodic photonic structures containing regular three-dimensional (3D) optical materials, such as diffraction gratings and periodic metamaterials. The proposed method builds upon the rigorous coupled-wave analysis and incorporates the nonlinear optical response of 2D materials by means of modified electromagnetic boundary conditions. This allows one to reduce the mathematical framework of the numerical method to an inhomogeneous scattering matrix formalism, which makes it more accurate and efficient than previously used approaches. An overview of linear and nonlinear optical properties of graphene and TMDC monolayers is given and the various features of the corresponding optical spectra are explored numerically and discussed. To illustrate the versatility of our numerical method, we use it to investigate the linear and nonlinear multiresonant optical response of 2D-3D heteromaterials for enhanced and tunable second- and third-harmonic generation. In particular, by employing a structured 2D material optically coupled to a patterned slab waveguide, we study the interplay between geometric resonances associated to guiding modes of periodically patterned slab waveguides and plasmon or exciton resonances of 2D materials., Comment: 28 pages, 21 figures

We perform a comparative study of the Anderson localization of light beams in disordered layered photonic nanostructures that, in the limit of periodic layer distribution, possess either a Dirac point or a Bragg gap in the spectrum of the wavevectors. In particular, we demonstrate that the localization length of the Anderson modes increases when the width of the Bragg gap decreases, such that in the vanishingly small bandgap limit, namely when a Dirac point is formed, even extremely high levels of disorders are unable to localize the optical modes located near the Dirac point. A comparative analysis of the key features of the propagation of Anderson modes formed in the Bragg gap or near the Dirac point is also presented. Our findings could provide valuable guidelines in assessing the influence of structural disorder on the functionality of a broad array of optical nanodevices., Comment: 9 pages, 5 figures, to appear in Scientific Reports

The rigorous coupled-wave analysis (RCWA) is one of the most successful and widely used methods for modeling periodic optical structures. It yields fast convergence of the electromagnetic far-field and has been adapted to model various optical devices and wave configurations. In this article, we investigate the accuracy with which the electromagnetic near-field can be calculated by using RCWA and explain the observed slow convergence and numerical artifacts from which it suffers, namely unphysical oscillations at material boundaries due to the Gibb's phenomenon. In order to alleviate these shortcomings, we also introduce a mathematical formulation for accurate near-field calculation in RCWA, for one- and two-dimensional straight and slanted diffraction gratings. This accurate near-field computational approach is tested and evaluated for several representative test-structures and configurations in order to illustrate the advantages provided by the proposed modified formulation of the RCWA., Comment: 13 pages, 12 figures

We demonstrate that plasmonic helical gratings consisting of metallic nanowires imprinted with helical grooves or ridges can be used efficiently to generate plasmonic vortices with radius much smaller than the operating wavelength. In our proposed approach, these helical surface gratings are designed so that plasmon modes with different azimuthal quantum numbers (topological charge) are phase-matched, thus allowing one to generate optical plasmonic vortices with arbitrary topological charge. The general principles for designing plasmonic helical gratings that facilitate efficient generation of such plasmonic vortices are derived and their applicability to the conversion of plasmonic vortices with zero angular momentum into plasmonic vortices with arbitrary angular momentum is illustrated in several particular cases. Our analysis, based both on the exact solutions for the electromagnetic field propagating in the helical plasmonic grating and a coupled-mode theory, suggests that even in the presence of metal losses the fundamental mode with topological charge m=0 can be converted to plasmon vortex modes with topological charge m=1 and m=2 with a conversion efficiency as large as 60%. The plasmonic nanovortices introduced in this study open new avenues for exciting applications of orbital angular momentum in the nanoworld., Comment: 10 pages, 7 figures, to appear in Scientific Reports

We present a comprehensive theoretical analysis and computational study of four-wave mixing (FWM) of optical pulses co-propagating in one-dimensional silicon photonic crystal waveguides (Si-PhCWGs). Our theoretical analysis describes a very general set-up of the interacting optical pulses, namely we consider nondegenerate FWM in a configuration in which at each frequency there exists a superposition of guiding modes. We incorporate in our theoretical model all relevant linear optical effects, including waveguide loss, free-carrier (FC) dispersion and FC absorption, nonlinear optical effects such as self- and cross-phase modulation (SPM, XPM), two-photon absorption (TPA), and cross-absorption modulation (XAM), as well as the coupled dynamics of FCs and optical field. In particular, our theoretical analysis based on the coupled-mode theory provides rigorously derived formulae for linear dispersion coefficients of the guiding modes, linear coupling coefficients between these modes, as well as the nonlinear waveguide coefficients describing SPM, XPM, TPA, XAM, and FWM. In addition, our theoretical analysis and numerical simulations reveal key differences between the characteristics of FWM in the slow- and fast-light regimes, which could potentially have important implications to the design of ultra-compact active photonic devices., Comment: 19 pages, 14 figures

We demonstrate that graphene-based photonic superlattices provide a versatile platform for electrical and all-optical control of photonic beams with deep-subwavelength accuracy. Specifically, by inserting graphene sheets into periodic metallo-dielectric structures one can design optical superlattices that posses photonic Dirac points (DPs) at frequencies at which the spatial average of the permittivity of the superlattice, $\bar{ \varepsilon}$, vanishes. Similar to the well-known zero-$\bar{n}$ bandgaps, we show that these zero-$\bar{\varepsilon}$ DPs are highly robust against structural disorder. We also show that, by tuning the graphene permittivity via the optical Kerr effect or electrical doping, one can induce a spectral variation of the DP exceeding \SI{30}{\nano\meter}, at mid-IR and THz frequencies. The implications of this wide tunability for the photonic Zitterbewegung effect in a vicinity of the DP are explored too., Comment: 5 pages, 5 figures, to appear in Phys. Rev. B as a Rapid Communication

We introduce a versatile numerical method for modeling light diffraction in periodically patterned photonic structures containing quadratically nonlinear non-centrosymmetric optical materials. Our approach extends the generalized source method to nonlinear optical interactions by incorporating the contribution of nonlinear polarization sources to the diffracted field in the algorithm. We derive the mathematical formalism underlying the numerical method and introduce the Fourier-factorization suitable for nonlinear calculations. The numerical efficiency and runtime characteristics of the method are investigated in a set of benchmark calculations: the results corresponding to the fundamental frequency are compared to those obtained from a reference method and the beneficial effects of the modified Fourier-factorization rule on the accuracy of the nonlinear computations is demonstrated. In order to illustrate the capabilities of our method, we employ it to demonstrate strong enhancement of second-harmonic generated in one- and two-dimensional optical gratings resonantly coupled to a slab waveguide. Our method can be easily extended to other types of nonlinear optical interactions by simply incorporating the corresponding nonlinear polarization sources in the algorithm., Comment: 12 pages, 8 figures