• Photonics Research
  • Vol. 10, Issue 2, 574 (2022)
Flore Hentinger1, Melissa Hedir1, Bruno Garbin1, Mathias Marconi2, Li Ge3、4, Fabrice Raineri1、5, Juan A. Levenson1, and Alejandro M. Yacomotti1、*
Author Affiliations
  • 1Centre de Nanosciences et de Nanotechnologies, CNRS, Université Paris-Sud, Université Paris-Saclay, 91120 Palaiseau, France
  • 2Université Côte d’Azur, Institut de Physique de Nice, CNRS-UMR 7010, Sophia Antipolis, France
  • 3Department of Physics and Astronomy, College of Staten Island, CUNY, Staten Island, New York 10314, USA
  • 4Graduate Center, CUNY, New York, New York 10016, USA
  • 5Université de Paris, 75205 Paris Cedex 13, France
  • show less
    DOI: 10.1364/PRJ.440050 Cite this Article Set citation alerts
    Flore Hentinger, Melissa Hedir, Bruno Garbin, Mathias Marconi, Li Ge, Fabrice Raineri, Juan A. Levenson, Alejandro M. Yacomotti. Direct observation of zero modes in a non-Hermitian optical nanocavity array[J]. Photonics Research, 2022, 10(2): 574 Copy Citation Text show less
    NH zero mode in a three coupled cavities array: CMT model. (a) NH Hamiltonian. (b) Black symbols and lines: eigenvalues in the CMT model with τ0=0.2 and g=12.3 (inset: intensity and phase spatial distribution at M2 laser threshold). Red symbols and lines: eigenvalues in the CD-CMT model, for increasing carrier density N with a Gaussian profile N(x;x1)=exp[−(x−x1)2/σ2]N centered at cavity 1 (τc=7 ps, x1=60, σ=50, αH=2; other parameters can be found in Appendix A). (c) Logarithmic spectral intensity as a function of the center of a Gaussian pump spot, computed from Eq. (2) in the text. The pump spot profile is P(x;X)=exp[−(x−X)2/σ2], and γn(X)=P(xn;X), with σ=50 and γ=6, i.e., below laser threshold. Solid lines correspond to Re(εj) for j=1 (M1, blue), j=2 (M2, burgundy), and j=3 (M3, orange).
    Fig. 1. NH zero mode in a three coupled cavities array: CMT model. (a) NH Hamiltonian. (b) Black symbols and lines: eigenvalues in the CMT model with τ0=0.2 and g=12.3 (inset: intensity and phase spatial distribution at M2 laser threshold). Red symbols and lines: eigenvalues in the CD-CMT model, for increasing carrier density N with a Gaussian profile N(x;x1)=exp[(xx1)2/σ2]N centered at cavity 1 (τc=7  ps, x1=60, σ=50, αH=2; other parameters can be found in Appendix A). (c) Logarithmic spectral intensity as a function of the center of a Gaussian pump spot, computed from Eq. (2) in the text. The pump spot profile is P(x;X)=exp[(xX)2/σ2], and γn(X)=P(xn;X), with σ=50 and γ=6, i.e., below laser threshold. Solid lines correspond to Re(εj) for j=1 (M1, blue), j=2 (M2, burgundy), and j=3 (M3, orange).
    Three coupled PhC cavities. (a) Artist view of the system, featuring controllable coupling by means of the two barriers (highlighted with dashed boxes) in which holes are modified. Two sublattices A and B can be defined, where couplings only take place between cavities belonging to different sublattices. Bottom: schematic representation showing how the presence of the barriers modifies the cavity detuning. The sublattice detuning is then Δω. (b) The 3D-FDTD simulation results of three coupled PhC cavities showing the evolution of the mode resonant wavelengths (symbols, left axis) as a function of the barrier parameter. Solid lines (right axis) are CMT predictions using polynomial approximations of g(h) and Δ(h) obtained from two coupled cavity FDTD simulations (see Appendix B). Two regions can be distinguished: the low (light green, |g|≳|Δω|) and high sublattice detuning regions (light blue, |g|<|Δω|); the near-zero coupling region is highlighted with a striped background. (c) SEM image of a sample of three coupled L3 InP-based PhC cavities. Inset: QW PL. (d) Normalized reflectivity spectra as a function of the barrier parameter from resonant scattering experiments. These results can be interpreted as the optical response of the coupled cavity system in the Hermitian limit.
    Fig. 2. Three coupled PhC cavities. (a) Artist view of the system, featuring controllable coupling by means of the two barriers (highlighted with dashed boxes) in which holes are modified. Two sublattices A and B can be defined, where couplings only take place between cavities belonging to different sublattices. Bottom: schematic representation showing how the presence of the barriers modifies the cavity detuning. The sublattice detuning is then Δω. (b) The 3D-FDTD simulation results of three coupled PhC cavities showing the evolution of the mode resonant wavelengths (symbols, left axis) as a function of the barrier parameter. Solid lines (right axis) are CMT predictions using polynomial approximations of g(h) and Δ(h) obtained from two coupled cavity FDTD simulations (see Appendix B). Two regions can be distinguished: the low (light green, |g||Δω|) and high sublattice detuning regions (light blue, |g|<|Δω|); the near-zero coupling region is highlighted with a striped background. (c) SEM image of a sample of three coupled L3 InP-based PhC cavities. Inset: QW PL. (d) Normalized reflectivity spectra as a function of the barrier parameter from resonant scattering experiments. These results can be interpreted as the optical response of the coupled cavity system in the Hermitian limit.
    Spatially resolved PL measurements in the large detuning regime. (a) Experimental results showing spectral intensity maps upon spatial scanning of a pump spot for a=416 nm and h=−20%, Ppump=0.8 μW. The position of the sample is changed by means of the piezo-electric voltage. (b) CD-CMT predictions (σ=40, N=0.22, αH=3). Near-field [(a)-i to (a)-v, (b)-i to (b)-iii] and far-field [(a)-vi to (a)-x, (b)-iv to (b)-vi] images are displayed at the selected pump spot positions marked with horizontal dashed lines. Vertical white dashed lines in the near field indicate the position of the cavities, and kx=0 in the far field.
    Fig. 3. Spatially resolved PL measurements in the large detuning regime. (a) Experimental results showing spectral intensity maps upon spatial scanning of a pump spot for a=416  nm and h=20%, Ppump=0.8  μW. The position of the sample is changed by means of the piezo-electric voltage. (b) CD-CMT predictions (σ=40, N=0.22, αH=3). Near-field [(a)-i to (a)-v, (b)-i to (b)-iii] and far-field [(a)-vi to (a)-x, (b)-iv to (b)-vi] images are displayed at the selected pump spot positions marked with horizontal dashed lines. Vertical white dashed lines in the near field indicate the position of the cavities, and kx=0 in the far field.
    Observation of the zero mode in the low detuning regime. (a) Experimental results showing spectral intensity maps upon spatial scanning of a pump spot for a=408 nm and h=0%, Ppump=0.8 μW. The position of the sample is changed by means of the piezo-electric voltage. (b) CD-CMT predictions (σ=40, N=0.27, αH=3). Near-field [(a)-i to (a)-iv, (b)-i to (b)-iii] and far-field [(a)-v to (a)-viii, (b)-iv to (b)-vi] images are displayed at the selected pump spot positions marked with horizontal dashed lines. Spectral filters are used in order to remove contributions from other modes, the spectral bandwidth being represented by the horizontal extension of dashed boxes. Vertical white dashed lines in the near field indicate the position of the cavities, and kx=0 in the far field.
    Fig. 4. Observation of the zero mode in the low detuning regime. (a) Experimental results showing spectral intensity maps upon spatial scanning of a pump spot for a=408nm and h=0%, Ppump=0.8  μW. The position of the sample is changed by means of the piezo-electric voltage. (b) CD-CMT predictions (σ=40, N=0.27, αH=3). Near-field [(a)-i to (a)-iv, (b)-i to (b)-iii] and far-field [(a)-v to (a)-viii, (b)-iv to (b)-vi] images are displayed at the selected pump spot positions marked with horizontal dashed lines. Spectral filters are used in order to remove contributions from other modes, the spectral bandwidth being represented by the horizontal extension of dashed boxes. Vertical white dashed lines in the near field indicate the position of the cavities, and kx=0 in the far field.
    Two coupled PhC cavities’ 3D FDTD numerical simulations as a function of the barrier perturbation parameter h (see main text). (a) Hybrid mode frequencies. (b) Coupling and cavity detuning extracted from the hybrid mode frequency splitting and average, respectively. Solid lines are third-order polynomial fittings, which are subsequently used for CMT calculation with three coupled cavities in the main text.
    Fig. 5. Two coupled PhC cavities’ 3D FDTD numerical simulations as a function of the barrier perturbation parameter h (see main text). (a) Hybrid mode frequencies. (b) Coupling and cavity detuning extracted from the hybrid mode frequency splitting and average, respectively. Solid lines are third-order polynomial fittings, which are subsequently used for CMT calculation with three coupled cavities in the main text.
    Single PhC cavity 3D FDTD numerical simulations in the presence of one or two barriers at the sides.
    Fig. 6. Single PhC cavity 3D FDTD numerical simulations in the presence of one or two barriers at the sides.
    Phase diagram underlying the transition from sublattice delocalization and zero modes to sublattice localization and mode coalescence. (a)–(e) Experimental PL intensity maps under pump spot position scanning across the coupled cavity system. (a)–(c) a=416 nm, Ppump=0.8 μW; (d), (e) a=408 nm, Ppump=1.1 μW. (f)–(j) CD-CMT and (k)–(o) linear CMT calculations with parameters αH=3, σ=40. (f)–(h) N=0.23, (i) N=0.24, and (j) N=0.27; (k)–(m) γ=1, (n) γ=1.1, and (o) γ=1.3.
    Fig. 7. Phase diagram underlying the transition from sublattice delocalization and zero modes to sublattice localization and mode coalescence. (a)–(e) Experimental PL intensity maps under pump spot position scanning across the coupled cavity system. (a)–(c) a=416  nm, Ppump=0.8  μW; (d), (e) a=408  nm, Ppump=1.1  μW. (f)–(j) CD-CMT and (k)–(o) linear CMT calculations with parameters αH=3, σ=40. (f)–(h) N=0.23, (i) N=0.24, and (j) N=0.27; (k)–(m) γ=1, (n) γ=1.1, and (o) γ=1.3.
    Zero-mode extinction due to full decoupling. (a) Experimental result on a sample with h=−20% featuring fabrication imperfections. A detuning of about 1 nm is observed between the extreme cavities. Notice that although three modes are observed around a pump position of 12 V, the three modes are localized into each cavity featuring different frequencies. The middle cavity is decoupled from the extreme one because of large negative h (localized phase). (b) CD-CMT simulations with h=−20%, N=0.225, αH=3, σ=40, and the cavity centered at x=140 is red-detuned by −6τ−1, equivalent to Δλ∼1.08 nm.
    Fig. 8. Zero-mode extinction due to full decoupling. (a) Experimental result on a sample with h=20% featuring fabrication imperfections. A detuning of about 1 nm is observed between the extreme cavities. Notice that although three modes are observed around a pump position of 12 V, the three modes are localized into each cavity featuring different frequencies. The middle cavity is decoupled from the extreme one because of large negative h (localized phase). (b) CD-CMT simulations with h=20%, N=0.225, αH=3, σ=40, and the cavity centered at x=140 is red-detuned by 6τ1, equivalent to Δλ1.08  nm.
    HermitianNon-Hermitian
    Definitionεj=0Re(εj)=0
    Symmetry of the entire spectrumChiral (εk=εj)Non-Hermitian particle-hole (εk=εj*)
    Dimension of existence in the complex plane0D (at the origin)1D (along the imaginary axis)
    WavefunctionDark on one sublatticeπ/2 phase shift between sublattices
    Excitation methodResonant excitation (fluorescence)Resonant excitation or incoherent pumping (lasing zero mode)
    Robustness againstCoupling disorderCoupling disorder and/or gain/loss disorder
    Table 1. Comparative Table with the Salient Properties of Hermitian and Non-Hermitian Photonic Zero Modes
    Flore Hentinger, Melissa Hedir, Bruno Garbin, Mathias Marconi, Li Ge, Fabrice Raineri, Juan A. Levenson, Alejandro M. Yacomotti. Direct observation of zero modes in a non-Hermitian optical nanocavity array[J]. Photonics Research, 2022, 10(2): 574
    Download Citation