# 99% beta factor and directional coupling of quantum dots to fast light in photonic crystal waveguides determined by hyperspectral imaging

###### Abstract

Spontaneous emission from excitonic transitions in InAs/GaAs quantum dots embedded in photonic crystal waveguides at 5 K into non-guided and guided modes is determined by direct hyperspectral imaging. This enables measurement of the absolute coupling efficiency into the guided modes, the beta factor, directly, without assumptions on decay rates used previously. Notably, we found beta factors above 90% over a wide spectral range of 40 meV in the fast light regime, reaching a maximum of ()%. We measure the directional emission of the circularly polarized transitions in a magnetic field into counter-propagating guided modes, to deduce the mode circularity at the quantum dot sites. We find that points of high directionality, up to 97%, correlate with a reduced beta factor consistent with their positions away from the mode field antinode. By comparison with calibrated finite-difference time-domain simulations, we use the emission energy, mode circularity and beta factor to estimate the quantum dot position inside the photonic crystal waveguide unit cell.

## I Introduction

Quantum dots (QDs) embedded in photonic crystal waveguides (PCWGs) are a promising system to implement quantum technologies. Due to the broadband coupling Rao and Hughes (2007); Lecamp et al. (2007), the system can be used for high-efficiency on-chip single photon sources. Furthermore, the strong lateral confinement of light results in a significant longitudinal component of the electromagnetic mode field, which allows for local circular polarization, and therefore selective coupling of circularly polarized dipoles into a single mode Rodríguez-Fortuo et al. (2013); Junge et al. (2013); Mitsch et al. (2014); Söllner et al. (2015); Young et al. (2015); Lodahl et al. (2017). Recent experiments have shown QD spin-photon path conversion and photon path-dependent QD spin initialization using this mechanism Coles et al. (2016, 2017), which is robust against disorder Lang et al. (2015). In conjunction with the recently demonstrated spin-controlled photon switching Javadi et al. (2018a) and super-radiant emission from two coupled QDs in a PCWG Kim et al. (2018), these results show the potential of such a system for the implementation of scalable quantum technologies on chip Lodahl et al. (2015); Mahmoodian et al. (2016). A fundamental requisite of quantum technology based on QD on PCWGs is that the spontaneous emission (SE) from the emitter couples exclusively to the designed channels of the system, which are typically waveguide (WG) modes, and not to other background channels creating losses. For a given WG mode, the probability of a QD exciton (QDE) to emit into the mode is called the beta factor, defined as

(1) |

where and are the decay rates to the selected WG mode, and other, typically non-guided, modes, respectively, and is the non-radiative decay rate. In previous works, the beta factor was estimated using lifetime measurements Lund-Hansen et al. (2008); Thyrrestrup et al. (2010); Dewhurst et al. (2010); Laucht et al. (2012a); Hoang et al. (2012); Arcari et al. (2014). The challenge in these measurements is to determine all three decay rates. A common approach assumed to be given by the difference between the decay rate of a QDE coupled to the WG mode, , and the decay rate of similar QDEs not coupled to the WG, , so that , and . This accuracy of this analysis depends on the accuracy of the underlying assumption that all QDEs have the same decay rate into non-guided modes, and the same non-radiative decay. The non-radiative decay is dependent on QD charging and local defects and is thus determined by properties beyond the photonic environment. Assuming no influence of local defects or tunneling, it is negligible for neutral excitons in InAs/GaAs QDs. Even in charged excitons, where Auger processes provide a non-radiative decay, this rate in the order of s Kurzmann et al. (2016) is three orders of magnitude below typical radiative decay rates of 1 ns Langbein et al. (2004a). The radiative decay into other modes , however, is likely to be significantly modified by the local dielectric environment of the QDs in PCWG structures, as recent calculations Javadi et al. (2018b) have highlighted. Therefore, the analysis reported in previous work is expected to exhibit significant systematic errors in the determined beta factor, as recently pointed out Rigal et al. (2018).

In the present work we use direct hyperspectral imaging to determine the emitted power, avoiding assumptions on decay rates altogether. The SE from QDs embedded in the PCWG along the WG is imaged onto the input slit of an imaging spectrometer, and is measured spatially and spectrally resolved. In this way, the SE guided by the WG and coupled at the ends of the WG by grating couplers into free space is measured together with the SE emerging from the QD directly into free space. The beta factor is determined using the relative emission powers, after correcting for the propagation losses and relative efficiencies of the couplers. Using an external magnetic field in Faraday (out of plane) direction, the QDE states are split into two spectrally resolvable transitions with opposite circular polarisation. Depending on the QD location within the unit cell of the PCWG, the two transitions couple differently to WG modes of opposite propagation direction. Using the spectrally resolved emission from the two couplers, the emission into the two counterpointing WG modes is measured. Using the powers of the two transitions emerging from the two couplers, we quantify the directionality Söllner et al. (2015); Coles et al. (2016) of the emission, and deduce the WG mode directionality at the QD site. The statistical distribution of the determined beta factor versus over a large ensemble of QDs, which are expected to randomly sample the in-plane area of the PCWG, shows that above 90% are mostly found for directionalities below 80%, and vice versa. The experimental results are analysed using detailed electromagnetic simulations of the PCWG structure. Specifically, we calculate and versus position inside the PCWG, and comparing with experiments, we estimate the QD position within the PCWG unit cell on a 10 nm length scale. This position is expected to affect the exciton dephasing, both via the interaction with surface states, and through a modification of the local phonon density of states Tighineanu et al. (2018).

## Ii Sample and Methods

The investigated sample is a GaAs photonic crystal membrane of 125 nm nominal thickness, with a single layer of InAs/GaAs QDs at the center. The QD area density was about cm, corresponding to an average distance of about 300 nm, or about one QD per unit cell of the PCWG. The QDs are n-doped using a Si delta-doping layer with a density of about cm, 10 nm below the QDs. Given the inhomogeneous size and spatial distribution of the QDs, this leads to a distribution of QD electron charging. Multiply charged QDs provide broader SE multiplets due to final state damping and spin-splitting. The emission lines analyzed were sharp lines which did not show a resolvable fine-structure splitting. Since typical values of fine-structure splittings for these QDs are in the few 10 μeV range Langbein et al. (2004b), and considering the doping, we attribute these lines to negatively charged exciton transitions. Even though the charging is relevant for applications using the electron spin Young et al. (2015); Javadi et al. (2018a), it is not important for the experimental results shown in this work. A fine-structure splitting much smaller than the magnetic-field induced Zeeman splitting of around 50 μeV, however, is important for the purity of the circular polarization of the transitions in a magnetic field Bayer et al. (2002).

The PCWG is created by a line of missing holes in the periodic hexagonal pattern of round air holes of separation , giving rise to guided modes within the two-dimensional photonic band gap for modes with dominating in-plane electric fields Johnson et al. (1999). An illustration of the investigated PCWGs is shown in Fig. 1a. The PCWGs investigated here are 100 long, with nm, and have 6 rows of holes on each side. The PCWGs are terminated by Bragg reflector couplers, which couple light propagating in the PCWG to free space Faraon et al. (2008). Two PCWGs (called and ) have been analyzed, which differ by their hole radius to period ratio , being 0.24 and 0.26, respectively.

The sample is mounted in a low vibration closed cycle cryostat (Montana Cryostation) on a XYZ piezo stage (Attocube) with a spatial resolution around 0.1 μm, allowing focusing and lateral alignment. For the measurements shown in this work, the sample temperature was 5 K. A microscope objective (MO) with a numerical aperture NA of 0.85 is mounted inside the cryostat on the cold shield, having a temperature of about 30 K. This avoids sample heating by thermal radiation, which can be an issue, specifically with membrane structures due to their reduced thermal contact to the substrate, when using the more common geometry with a room temperature window close to the sample and an external MO. A permanent magnet can be mounted below the sample to provide a magnetic field of =0.45 T, where labels the direction normal to the sample plane, as indicated in Fig. 1a.

The QD SE was excited by a laser at a wavelength of nm focused onto the sample to a sub-μm spot by the MO (see also appendix F.1). The SE is filtered with a colour filter (Schott RG680) transmitting wavelengths above 680 nm. The SE from the QDs can be visualized by two imaging cameras, one for the near field (NF), imaging the real space at the sample, and one for the far field (FF), imaging the reciprocal space at the sample. For spectral imaging, the real space is imaged onto the input slit of an imaging spectrometer with a focal length of 1.9 m, a 1200 l/mm holographic grating of () mm size, 900 nm blaze wavelength, and detected by a CCD (Roper Pixis) of square pixels of 20 m size. For all the measurements performed in this work, the input slit aperture was 20 m, corresponding to 639 nm at the sample plane. The corresponding spectral resolution (full width at half maximum (FWHM)) is 8 μeV at 880 nm (1.41 eV). A white light Köhler illumination is integrated with the main PL setup, to simultaneously visualize the PCWG sample and the PL emission from the QDs. The origin of the -axis along the WG is chosen at the center of the WG.

A reflection image of the sample is shown in Fig. 1b. The red spot indicates the excitation laser at a specific position along the WG. The excitation at a photon energy of 1.96 eV (633 nm wavelength) creates electron-hole pairs in the GaAs membrane, which subsequently relax by phonon emission towards the GaAs band gap around 1.52 eV, before being captured into the highly strained InGaAs wetting layer, where they further relax to the wetting layer band gap around 1.42 eV, and finally into the QDs, which emit in the energy range between 1.37 eV and 1.41 eV. The direct imaging allows measurement of the diffusion length in the planar region of the sample, which was found to be about 3-4 μm on the membrane, as shown in appendix F.1. The described carrier relaxation process is complex, with the formation of excitons in GaAs and the wetting layer also playing a role. Once captured into a QD, the carriers relax to the ground state within tens of picoseconds, from which they radiatively recombine, emitting a photon into the local photonic mode structure, consisting of the WG modes, confining light, and non-guided modes, rapidly escaping to free space on both sides of the slab. The emission is imaged onto the input slit of the spectrometer, indicated as black dashed line in Fig. 1b. Importantly, the WG has been aligned along the slit, in order to collect and image the emission from the whole WG and the couplers. The SE spectrum from the QDs is therefore detected spatially resolved along the WG, as exemplified in Fig. 1c, covering an energy range of about 15 meV for a given spectrometer center position. Exciting at μm, we observe most of the free space emission close to the excitation, while signals around μm arise from the WG couplers, representing the QD SE into the WG modes. The spectrum from the bottom WG coupler is indicated as a white line, integrated over μm. The center signal, close to the excitation position, is attributed to QD SE into the non-guided modes, excited by the carrier relaxation and diffusion processes described before. However, we also observed QD free-space emission from positions far away from the excitation spot, which we attribute to indirect excitation, where the wetting layer emission, which is coupled into the WG mode, is absorbed by QD excited states (see appendix B).

## Iii Photonic band structure

Determining the photonic band structure of the investigated sample is crucial for a quantitative comparison with numerical simulations. In the literature, near-field scanning optical spectroscopy Sukhorukov et al. (2009) and interferometric techniques Gersen et al. (2005) have been used to determine the guided mode dispersion in PCWGs. In other cases, the simulation parameters were adjusted to approximately reproduce the measured transmission window Laucht et al. (2012b). Alternatively, Fabry-Pérot fringes were used to calculate the group index of the guided mode Arcari et al. (2014).

Here, we use Fourier imaging Langbein et al. (1999); Langbein and Hvam (2002) to directly measure the band dispersion within the light cone Thomas et al. (2007). We use pulsed laser excitation, with center energies corresponding to band wavevectors within the light cone, coupled to the WG via the bottom coupler. The laser polarization is set orthogonal to the WG direction, in order to select the fundamental mode (see appendix G.3). To suppress the reflection of the excitation laser, we use a rectangular aperture in an intermediate image plane, corresponding to (8.7,19.0) μm size at the PCWG. We measure both the NF and the FF of the emission along the WG as a function of energy, from which we can determine the propagation losses and the wavevector of the corresponding Bloch mode. An example of a NF measurement for =1.4565 eV is shown in Fig. 2b, with the corresponding FF measurement shown in Fig. 2c. The NF profile along the WG exhibits significant fluctuations, reflecting the coherent nature of the emission. After background subtraction, it can be fit with an exponential decay, as shown in the right panel of Fig. 2b, from which we obtain the WG loss , as further detailed in appendix F. The resulting are given in the inset of Fig. 2a as function of , showing an increase with , which is attributed to both the radiative losses of the WG mode in the light cone and the absorption by the wetting layer, as discussed in more detail in appendix D.

Turning to the FF measurements, we note that the accessible range of the in-plane wavevector is limited by the NA of the MO to , with the free-space wavevector . The resulting cut-off is visible in Fig. 2c and from this we can calibrate the -space of our measurements (see appendix F.3). The measured FF pattern shows a stripe at around μm, elongated in direction, which is the WG radiation loss. The WG mode wavevector is given by , and the large extension in is due to the small extension of the WG mode in . Knowing that the excited WG mode is propagating in positive direction, it is interesting to note that it exhibits a negative wavevector – clear evidence of the negative group velocity of the WG mode, so that the phase velocity, given by , is opposite to the group velocity, which is along the propagation direction.

The finite width in is due to the exponential decay of the field along the direction, the finite size of the imaged region, and the finite bandwidth of the excitation laser, as discussed in appendix F.2. We determine the propagation wavevector by fitting the FF profile with a Lorentzian, as shown in the right panel of Fig. 2c. The FF profile is obtained by averaging over , indicated by the white dashed lines. The resulting are shown in Fig. 2a, with the blue circles, noting that the reflection symmetry of the WG allows to use the absolute value.

We now use this measurement of the WG mode wavevectors to fine-tune the sample parameters used in simulations. The membrane thickness and hole diameters have been measured by scanning electron microscopy (SEM), and the photonic crystal (PC) period has a negligible fabrication error. We note that the refractive index of GaAs has previously been used as adjustable parameter, taking low temperature values of Rao and Hughes (2007), 3.5 Laucht et al. (2012a), and 3.45 Arcari et al. (2014). However, the GaAs refractive index is well known Gehrsitz et al. (2000) and somewhat larger than these values. We therefore use in our simulations the known refractive index at low temperatures, including its dispersion. The QD layer embedded in the slab corresponds to less than 1% of its thickness and is made of a similar material as GaAs, so that we neglect its effect on the refractive index.

It is known that fabrication by selective etching, and oxidation of GaAs over time Lukeš (1972); DeSalvo et al. (1996), can remove a surface layer of GaAs. We therefore use an effective thickness of a removed surface layer as a parameter, to match the simulations to the measurements, as detailed in appendix G.1. The calculated band structure for nm is shown in Fig. 2a, together with a typical free space SE of the QDs embedded in the PCWG. The measured spectrum consists of sharp lines at energies below 1.41 eV, and a broad emission at higher energies, which we attribute to the wetting layer. Sharp lines superimposed to the wetting layer emission are attributed to localized excitons.

## Iv Results and discussion

Using the spectral imaging as shown in Fig. 1, we identify the SE of individual QDs from the top and bottom couplers and from the QD position into non-guided modes. The excitation position can be adjusted along the WG to maximize the QD emission. The measurements were performed at a low excitation power to avoid multiexciton emission.

Examples of the SE detected for three different QDs are shown in Fig. 3a-c, over a small spectral range covering the splitting of μeV created by the magnetic field. As the exciton has spin projection , the Zeeman splitting between these states is given by , with the Bohr magneton. Using T, we calculate an exciton g-factor , comparable to the value of 1.2 reported in Coles et al. (2016). The beta factor can be calculated as the ratio between the power emitted from the couplers and the total emitted power as

(2) |

where , , are the detected SE powers of a single QD from the bottom and top couplers and at the QD site, respectively, corrected for losses and relative coupler efficiencies (see appendix D), and labels the helicity of the QD transition. We note that we have not corrected the values for the simulated collection efficiency of free-space and coupler emission discussed in appendix G.5, which would increase the beta factor, since the free space emission is collected typically twice as efficiently, depending on the QD position. Furthermore we determine the QD transition energy and its position along the WG, , as detailed in appendix A.

For the QD SE shown in Fig. 3a, we find . The QD is located close to the top coupler, as highlighted in the figure by the white dashed rectangle. Interestingly, the QD is visible for excitation at , which we attribute to reabsorption of the wetting layer emission propagating along the PCWG. We verified that is independent of the excitation position (see appendix B), as expected given that is determined by the local photon density of states, which is independent of the QD excitation pathway. Fig. 3b,c show the SE of QDs with and , respectively. The directional emission is evident in the asymmetry of the power from top and bottom coupler. The degree of directionality depends on the WG mode circularity at the QD site, which we calculate here from the measured powers as

(3) |

which by using the ratio of power ratios is independent of the efficiencies of the couplers. is related to the directionality used in Ref. [Coles et al., 2016] by , making it equal for small circularities but avoiding saturation for high circularities (see green line in Fig. 3d). At a circularly polarized point (C-point), diverges, while at a linearly polarized point (L-point), is zero. Residual reflections from the couplers might limit the maximum value of . From the maximum measured , we can deduce that these reflections are below 5%, see appendix E.

Notably, the of the QDs shown in Fig. 3a-c decreases with increasing . To investigate this further, we determined circularity and beta factor for many QDs, in and , as shown in Fig. 3d. We find that low are more likely at higher circularity. This is consistent with the fact that C-points require both transversal and longitudinal fields, and thus do not occur at the mode field anti-node. We find positive and negative circularities (red and blue symbols, respectively), randomly distributed for both investigated PCWGs. In the same figure, we show that the corresponding directionality is above 0.9 for some of the QDs, while in Ref. [Coles et al., 2016] the reported are below 0.8. Furthermore, we find above 0.9 over an emission energy range of about 40 meV (see Fig. 3e). In particular, in we find QDs with above 0.9 for energies between 1.372 eV and 1.408 eV. We note that in most of this spectral range both fundamental and the higher order mode are predicted to be present, see Fig. 2. The coupling of a QD to these two modes depends on its energy and its position within the PCWG unit cell, as discussed in appendix G.6. However, no obvious feature of the presence of the two modes is visible in the dependence of the beta factor on energy (Fig. 3e). Importantly, these results show that efficient coupling occurs even in region of low group index within the light cone. In Fig. 3f we plot as a function of the total SE power from the QD for excitation at , using for all QDs an excitation power of 0.5 μW at the sample. We find above 0.90 over three order of magnitudes in SE power. This indicates that the emission power is not dictated by the beta factor, but is governed mostly by the local carrier capture dynamics, which depends on the disorder landscape in the wetting layer close to the QDs. The circularity and QD energy also do not show a correlation with the coarse QD position along the WG, as shown in appendix C.

The WG mode circularity at the QD site, together with the QD beta factor and energy, identifies possible positions of the QD inside the PCWG unit cell. We therefore performed finite-difference time-domain (FDTD) simulations for with nm, extracting the simulated beta factor and circularity at positions across the PCWG unit cell, and at the energies of the measured QDEs, as detailed in appendix G. The simulated positions are organised on a checkerboard grid with a step of nm. For a comparison with the measured quantities and , the simulated beta factor , and circularity , have been corrected to take into account the projections of the QD emission onto the fundamental and higher modes and the corresponding polarisation resolved collection efficiencies, which include the couplers, the NA of the MO, the input slit of the spectrometer and the efficiency of the spectrometer, as detailed in appendix G, yielding and . To identify the position of the QD, we compare simulations and experiment taking into account the experimental errors, by evaluating the normalized residual for each QD, given by

(4) |

with

(5) |

In these expressions, and are the experimental errors (standard deviation) for the measured and , respectively. and are the simulated positions in the unit cell, according to the reference system given in Fig. 1a. Note that agreement between simulation and experiment within the estimated errors is obtained for values of of the order of unity. In general, the minimum of the residual indicates the most likely QD position inside the unit cell. In Fig. 4 the results are shown for QDs in , ordered with increasing energy. Generally, for most QDs we find positions with residuals below 5. Examples of QDs showing residuals below 2, labelled (a,b,d,e), enable to extract a well-defined position of the QD, with a precision of a few 10 nm considering the region over which is increasing by unity (i.e. one standard deviation) from its minimum. These positions are distributed over the unit cell, as expected. One example in (c), shows for all simulated positions – possibly, this QD is situated beyond the range covered in the simulations.

## V Conclusions

In conclusion, we have shown that direct spectral imaging allows measurement of the beta factor and the directional emission, without assumptions on radiative decay rates, by using only the relative powers emitted into the WG mode and to free space. We found a maximum beta factor of % in the fast light regime. Beta factors above 90% are mainly found for quantum dots located at sites with small WG mode circularity, consistent with the fact that circular points occur away from the field antinodes. Using Fourier imaging to measure the band dispersion of the WG mode within the light cone, we calibrate FDTD simulations, allowing us to locate the QD positions inside the PCWG unit cell with a few 10 nm precision, from their measured beta factor and circularity. These results are promising for the suitability of the system for photon blockade and more advanced quantum technology, and the methods presented can be used to identify suited QDs. Furthermore, the position determination can be used to determine the precision of QD site-control techniques, important for the development of useful and scalable devices.

## Acknowledgments

This work was supported by the EPSRC under grant EP/M020479/1, and partially supported by the GW4 Accelerator ”Southwest Quantum Dot Quantum Technologies GW4-AF3-004. The simulations were performed in the Advanced Computing Research Centre, University of Bristol. B.L. acknowledges support by an EPSRC studentship DTA 1407622. R.O. acknowledges her EPSRC fellowship EP/N003381/1. Assistance in sample growth and fabrication by Monika Emmerlin, Johannes Beetz and Sebastian Maier is gratefully acknowledged.

## Contribution statement

L.S. performed the measurements, analyzed the data, and drafted the manuscript, with support by F.M., W.L., and B.L.. B.L. performed the simulations, with support from D.B., E.A.M., L.S., and W.L.. A.B.Y., R.O., S.H., and C.S. provided the sample and participated in discussions and manuscript revision. W.L. guided the work.

## Appendix A Analysis of spectral images

Here we discuss the analysis of spectral imaging data to retrieve the emission powers used in Eq. (2) and Eq. (3). After subtracting the dark background, the data are divided by the integration time and multiplied by the CCD gain of 2 electrons per count, resulting in data as given in Fig. 5a in photoelectrons per time and pixel. All data were taken in a magnetic field of T, providing a Zeeman splitting into circularly polarized transitions, as visible in Fig. 5a. In order to determine the total detected power emitted from the QD directly in free space, the emission is integrated along the waveguide over the number of pixels contained in the region of interest. This region is highlighted in Fig. 5a with dashed green lines. In order to find the peak positions and areas, we first fit the QD free space emission using a sum of two Voigt profiles, with equal linewidth parameters, as shown in Fig. 5c. Using the obtained linewidths and energies, the emission from the couplers is fitted varying the amplitudes only. The areas of the fitted peaks, having units of photoelectrons per time, are taken as the measured powers , with and . These powers are corrected for WG losses and relative coupler efficiencies, as described in Sec. D, to obtain the powers used in Eq. (2).

In order to obtain the QD position along the WG, , we first separate the QD signal from the excitation background by averaging over two regions, one including the QD signal and the other one spectrally shifted up by about 100 μeV containing the spectrally broad background from other emitters, as exemplified by the white and orange dashed lines in Fig. 5a, respectively. Phonon-assisted transitions of the QD are spectrally broad and provide a background about two orders of magnitude below the zero-phonon line emission, which can be neglected in this analysis. Subtracting the excitation background, we isolate the QD signal, as shown in Fig. 5b, which is then fitted using three Gaussian peaks, to determine the position of the two couplers and of the QD free space emission . The axis is calibrated using the known distance between the couplers. The position error resulting from the fit is typically around 100 nm (standard deviation). Similarly, by fitting the background from other emitters, we determine the excitation position along the WG.

## Appendix B beta factor versus excitation position

In Fig. 5d we show the total SE and the beta factor as a function of the excitation position along the PCWG for one specific QD. The total SE is the sum of the SE emitted from the bottom and top couplers and from the QD position into free space, as in Fig. 3f. The total SE varies over three order of magnitude for different excitation positions along the WG, showing a maximum at the QD location determined from the free space SE as described in appendix A. The beta factor determined versus exctation position, is also shown in Fig. 5d, is constant within 1% standard deviation. This shows that the determination of the beta factor is robust against changes of excitation conditions. This variation is taken as a lower limit for the error in the beta factor. We note that out of the 31 analyzed QDs, only a few show significant non-local excitation and all of these emit between 1.393 eV and 1.377 eV. Non-local excitation is attributed to re-absorption of wetting layer emission into a WG mode. We can speculate that the QD absorption occurs into the trion p-shell (see Supplementary Material of Ref. [Coles et al., 2017] and Ref. [Benny et al., 2012]), which is about 40 meV above the s-shell, consistent with the energy separation between the wetting layer emission around 1.42 eV and the QD emission.

## Appendix C QD emission energy and WG mode circularity versus position

One would expect that there is no correlation of energy and WG mode circularity with the QD position along the PCWG. The experimental data are given in Fig. 6, and exhibit a random distribution of energy and circularity along the QD position, consistent with this expectation.

## Appendix D Relative efficiencies of couplers

To determine the relative efficiencies of the two couplers, we perform the following analysis on data taken with magnetic field. For a given excitation position along the PCWG, we define the distances and between the excitation position and the top and bottom couplers, respectively, with the length m of the PCWG, as sketched in the inset of Fig. 7a. Accounting for propagation losses , the SE powers emerging from the couplers are proportional to , with the efficiency of the couplers, where . Here we assume that most of the emission occurs at the excitation position, consistent with the experimental findings of the imaging along the WG, as discussed in appendix B. To apply the model, we spectrally integrate the SE from the bottom and top couplers over the range covered by the CCD camera ( meV), thereby averaging over the directionality of the individual QDs, resulting in the powers . The power ratio is fitted by

(6) |

with , as shown in Fig. 7a for the spectral range centered at =1.425 eV for . A good agreement between data and fit is observed. We repeated the analysis for the 4 spectral ranges considered in the experiment, and show in Fig. 7b-c the resulting and versus . We find that has no significant dependence on , and we use the average value of to correct the emission intensities of all the data from this WG. The analysis has been repeated with , giving a similar result . The loss increases with , as expected from the increasing absorption by the wetting layer and the increasing radiation losses (see appendix G.2). Notably, the values are consistent with the loss measured by NF imaging, shown as black circles in Fig. 7c.

The determined and are used to calculate the powers and unaffected by loss and relative coupling efficiency, used in Eq. (2), as

(7) | ||||

## Appendix E Back-reflections at PCWG termination

Back-reflections at the couplers can affect the observed circularity, since light emitted in a given direction is detected after the reflection at the opposite coupler. To estimate the influence of this effect for the investigated sample, we consider the measured powers of the QD with the highest mode circularity we found. To first order in the reflection coefficient , the reflection of the measured power from the bottom coupler is adding to the measured power from top coupler, and vice versa. The powers unaffected by the reflection are then given by , and . Requiring that these powers are positive, we find from the measured powers of all investigated QDs shown in Fig. 3 an upper limit for of about 3-5 %. The reflection leads to a systematic underestimation of the absolute circularity and directionality, and for 5% is limiting the measured directionality to and the measured circularity to . The dynamic range of our experimental data is sufficient to measure above 4, and above 99.9%, and is thus not limiting the results presented.

## Appendix F NF and FF analysis

In Fig. 2b in the main text, we show an example of the NF of the radiation losses along the waveguide. The data was obtained using the sum of the 100 subsequent frames of a video taken by a Sony DCR-TRV620E digital 8 camcorder in nightshot mode. The dark background was subtracted for all data shown. The NF emission has been derived as the difference between the emission profile along the waveguide and a background profile next to the waveguide. Each profile is averaged over 1.7 μm in the -direction, as indicated in Fig. 2b with white and red dashed lines respectively. The axis is calibrated using the known distance between the couplers. We repeated the analysis for different excitation energies, and calculated the corresponding loss coefficients. The result is shown in the inset of Fig. 2a. The origin of these losses is discussed in appendix G.2.

### f.1 Carrier diffusion length

The emission imaging allows to measure the carrier diffusion length on the photonic crystal membrane. In particular, we image the PL emission from the QDs for excitation about 3 μm sideways offset from WG, on the unstructured free-standing membrane, as shown in Fig. 8. The spatially resolved emission is extended compared to the excitation, due to carrier diffusion between excitation and recombination. In order to evaluate the diffusion length, we average over about 2.5 μm in the -direction (see white dashed lines). The corresponding profile, given as solid white line, is fitted with a Gaussian profile, showing a standard deviation of about 1.4 μm, which represents the carrier diffusion length along the direction. Some of the carriers diffuse into the WG region, exciting QDs that couple to the guided mode, as is evident from the weak emission observed from the couplers.

### f.2 Far-field width

The far field profile (see Fig. 2c) has a finite FWHM, , in the direction. There are two contributions to the observed width: (i) the spectral width of the excitation laser, which translates via the PCWG dispersion into a FWHM , and (ii) the spatial profile of the NF in the direction, which is given by an exponential decay of the field amplitude due to losses, delimited by the aperture size . The observed field amplitude is thus modelled, neglecting constant factors, as

(8) |

with the Heaviside function and the aperture length . Fourier-transforming and taking the absolute square, we find the corresponding intensity in -space is given by

(9) |

For large losses across the aperture, , the constant cosh term dominates and is a Lorentzian, while for small losses, , the cos term dominates, resulting in a sinc function. We fit Eq. (9) with a Lorentzian to determine the equivalent FWHM . The laser spectrum, converted into using the linear dispersion of about μm eV in the relevant range (see Fig. 2), and , together with the corresponding Lorentzian fits to determine , and , respectively, are shown in Fig. 9a for =1.4565 eV, as used in Fig. 2b-c. The FF profile is then given by the convolution of (i) and (ii), which for Lorentzians keeps a Lorentzian, having a width given by the sum of the widths, here . The result of this analysis is shown in Fig. 9b for all measured . Generally a good quantitative agreement is found. The remaining difference between and of about 0.1 μm could be related to a slight defocus of the FF imaging.

### f.3 Calibration of FF imaging

In order to calibrate the in-plane wavevector in the FF imaging, we measured the FF of the emission of the top coupler when exciting into the bottom coupler at =1.4102 eV, as shown in Fig. 10a. To use the known NA of the objective for calibration, we determine the cut-off radius, using the polar coordinate representation of the data, as shown in Fig. 10b. We adjusted the center of the coordinate system to obtain a constant maximum radius, shown by the vertical dashed line. We average the range between the two horizontal dashed lines in Fig. 10b, obtaining the profile shown in Fig. 10c. We take the NA radius to be the value of at half step height, and calibrate this radius to μm. A relative error of a few % rms of this calibration is estimated.

## Appendix G Simulation methods

### g.1 Geometry and Parameters

FDTD calculations using the package Meep Oskooi et al. (2010) were carried out to model the properties of the PCWGs under study. All calculations used a resolution of 24 points per lattice constant, and a cubic Yee-lattice. The permittivity of GaAs at low temperature K including its dispersion Gehrsitz et al. (2000) is given by

(10) | ||||

with the photon energy in units of eV. This permittivity was implemented in the FDTD dispersion model as shown in Fig. 11 over the relevant energy range.

Single unit cell simulations were carried out with periodic boundary conditions along and a current source at and a given position in . Initially these calculations made use of the nominal sample parameters: lattice constant nm, hole radius and slab height nm. Perfectly matched layers (PMLs) were placed on the out-of-plane simulation facets at and adiabatic absorbers Oskooi et al. (2008) on the in-plane facets at to absorb light tunneling through the cladding. Six layers of air holes were placed either side of the waveguide, matching the sample. A sketch of the geometry covering one unit cell of the PCWG is shown in Fig. 12.

To calculate the WG mode dispersion, we implemented periodic boundary conditions with a phase factor over a single PCWG unit cell, and extracted the mode resonances from the time-domain Fourier transforms as function of the phase.

Simulations with the nominal structural parameters showed significant deviations from the measured dispersion. As discussed in the main text, a single parameter was introduced, describing the thickness of material removed from all surfaces of the structure, for example by etching during fabrication or subsequent oxidation. Simulations were carried out with the radius of the holes expanded by , and the height of the slab reduced by . We found that around 7-8 nm produces a good match with experiment, see Fig. 13, and used nm for the remaining calculations, for both investigated WGs.

### g.2 PCWG loss

A significant number of QDs measured are in resonance with WG modes with propagation wavevectors within the light cone, with the mode eigenfrequency and the speed of light . Modes in this part of the dispersion emit into free space and thus experience propagation loss, even in an ideal structure. In addition to this out-of-plane loss, all modes will experience an in-plane loss due to the finite extension of the lateral PC, allowing light to tunnel through the PC cladding either side of the WG.

The loss was simulated with the method described in appendix G.1, choosing the symmetry of the simulation and sources either even or odd to excite only one of the two WG modes (see appendix G.3). The simulation geometry is as sketched in Fig. 12, with the height extended to and the PML thickness to , in order to allow for better separation of the out-of-plane loss from the in-plane propagation. We use a source with a Gaussian time dependence of 29 fs standard deviation, truncated at 5 standard deviations either side of the peak. To maximize excitation of the WG mode of interest, the source center frequency was set to the WG mode frequency at the simulated , as calculated in appendix G.1. The simulation was run for 1000 time-steps, ending about fs after the source is switched off, at which point the remaining field of non-guided modes is negligible, as is evidenced by a stable field distribution, apart from a global oscillating and decaying pre-factor with the complex mode frequency . The Poynting vector of this field was then evaluated, providing the power flux density of the selected WG mode. This distribution was used to calculate the flux through the plane, representing the propagating flux along the PCWG, and the flux through the planes representing the out-of-plane loss, and through the planes, representing the in-plane loss. The transmission coefficient per unit cell is then given by , from which the loss coefficient is determined, as shown in Fig. 14 for the two WG modes as a function of energy. The loss coefficient considering in-plane loss only (using ) is also given. The even (higher order) mode has a dispersion showing a maximum, and thus presents two -vectors for a given energy over a significant range. One branch corresponds to the part of the even mode inside the light cone, showing high loss, dominated by out-of-plane loss, while the other is outside the light cone and shows only the small remaining in-plane loss due to the tunnelling through the cladding. The odd (fundamental) mode shows increasing loss with increasing energy. Below 1.38 eV, the mode is outside the light cone, so that the loss is strongly reduced to the in-plane loss only. The measured loss is indicated by the circular data points, and show a good agreement with the calculated loss of the odd mode. The additional loss in the measurements above 1.43 eV is attributed to absorption in the wetting layer, not taken into account in the simulations.

From the same simulations we also obtain the mode field amplitude and polarisation. In Fig. 15, the absolute value of the electric field amplitude , the ellipticity and the orientation of the polarisation ellipse are shown in the plane , for . Note that for the field component is zero. Orientation and ellipticity are defined as

(11) | ||||

with the Stokes parameters , and and . We can clearly see the spatial extension of the modes in their amplitude, and their circular and linear points in their ellipticity. The fundamental mode is linear in the center (yellow in ), and circular of opposite helicity (black and white) in the center of the unit cell along , close to the first holes. Other linear and circular points are present, but at much reduced mode field amplitudes. The higher order mode at a small is dominated by linear polarization, and at large has circular points at the edge of the unit cell, at lower field strengths. This is a consequence of the temporal symmetry that ensures that no component of circular polarisation can exist at the band edge points and in this type of waveguide Lang et al. (2016). Thus in the proximity of these one finds either a reduced component of circular polarisation and/or a reduction in the electric field strength in the regions with circular polarisation.

To give some insight into the radiative loss of the modes, 2D spatial Fourier transforms were carried out on the fields of the guided modes. The fields in the =0 plane of the slab were used. In this plane, symmetry ensures that the three non-zero field components are , and . The profile of appears nearly identical to , and we choose to use for the following discussion , the Fourier transform with wavevector . Components at , inside the light cone , can couple to free space modes, having an intensity fraction , with , which is indicative of the loss rate per unit time Akahane et al. (2005). For this calculation, the mode fields were extracted from an eigensolver Johnson and Joannopoulos (2001), using , and nm. The fraction is given in Fig. 16a for both WG modes as function of the Bloch wavevector , and the distribution of is shown in Fig. 16b for selected .

The fundamental mode for small (see Fig. 16b panels I and II) is in the continuum above the bandgap of the PC cladding, and thus able to propagate in the directions, leading to a small extension in . This concentration reduces the overlap with the light cone, hence the low in Fig. 16a in this region. For higher (see panel III) the mode enters the bandgap, confining it in and thus broadening it in . This results in an increased , until leaves the light cone (panel IV). In contrast, the higher order mode (see bottom row of Fig. 16b) remains inside the bandgap for all , so that is decreasing monotonically with increasing , reaching zero at the edge of the light cone, .

In order to estimate the propagation loss coefficient per unit distance, the fraction is scaled as

(12) |

In this expression, the factor is the attempt rate with which the light is emitted, with the emission fraction per attempt. The factor 1/, with the group velocity , converts the resulting loss per unit time into the loss per unit propagation length. The resulting loss is given in Fig. 17. This approach only accounts for out-of-plane losses and assumes an infinite WG without absorption and disorder. The loss thus becomes exactly zero outside the light cone, beyond the lower end of the curves in Fig. 17. Notably, this estimate of loss is qualitatively reproducing the one calculated via FDTD (see Fig. 14) for both modes, and is quantitatively about a factor of 3 lower.

### g.3 Coupler efficiency

FDTD simulations were used to determine the efficiency of the couplers. A sketch of the simulation volume used is shown in Fig. 18. In order to calculate the efficiency and reflectivity of the coupler, two simulations are used. The coupler simulation contains the actual structure, and a calibration simulation is identical except that the coupler is removed and replaced by a continuation of the PCWG. The power in the PCWG mode travelling towards the coupler, called the input power , is determined by the calibration simulation, given by the power through the flux plane indicated in blue in Fig. 18. After determining the power through the same flux plane in the coupler simulation, the coupler reflectivity is calculated as . Using the power radiated out of the coupler, , through the flux plane indicated in red in Fig. 18, located 13 nm above the slab, the coupler efficiency is calculated as . Note that due to the inversion symmetry in , the same power is emitted towards both sides of the slab WG.

In the simulations, an electrical dipole source (yellow star in Fig. 18) with a time-varying current given by is used, with the time , the central frequency , the duration , and the dipole vector . Note that MEEP implements such a source using discrete time derivatives of Gaussians for better performance. The source frequency was set to , at the centre of the measured QD distribution ( eV). The source has a duration of fs standard deviation, which results in a frequency standard deviation of meV, covering the range of QD energies measured experimentally. To convert the simulated fields from time-domain to frequency-domain, we use MEEPs flux plane function Oskooi et al. (2010).

The source and simulation symmetry were chosen to select for either even or odd modes. For the odd mode, the source was polarised along the -direction (), and placed at coordinates relative to the centre of the unit cell (see Fig. 15). For the even mode, two in-phase -polarised sources are used () placed at . The symmetries were exploited to gain a factor of 4 reduction in simulation time and memory, using an even mirror plane at , and a mirror plane at with a parity matching the simulated mode Oskooi et al. (2010).

The simulation size is in , in , and in direction, as shown in Fig. 18. Thick adiabatic absorbers covering the reflecting simulation boundaries in direction were found to be necessary, since PCWG modes, specifically at low group velocity, are easily reflected from the spatially varying absorption in the absorbers. The absorbers are implemented via an electric and magnetic conductivity, and we used a scaling of this conductivity proportional to the sixth power of the depth into the absorber, which provided lower reflections compared to using second and forth power. For the odd mode, in the calibration simulation, the absorber thickness in the positive -direction was set to . The negative direction simulation boundary is less critical, as any reflections from this boundary will match between the calibration and coupler simulations, so that only the weaker reflection of the coupler reflection in the coupler simulation needs to be sufficiently suppressed. We used a thickness of . The other simulation facets do not have photonic structures intersecting them, and thus can be treated with shorter absorbers. We used in the directions, and PMLs of thickness in . In the coupler simulation, there is no photonic lattice at the boundary, and an absorber of thickness was used. The frequency window of interest includes slow group velocity regions of the even mode, making them even more sensitive to being reflected by absorbers. In the even-mode simulations, the absorber thickness on the boundaries were increased by . The absorber strength was set to 1/2 (1/18) of MEEP’s default Oskooi et al. (2010), for the odd (even) mode, respectively. We found that varied by about 0.7% in calibration simulations changing the absorber thickness on the facets to . Furthermore, displacing the source by 1 unit cell along the waveguide, altered the fluxes by about 0.7%. All relevant calculations were repeated with less absorber depth on the facets to ensure the accuracy of the results, and a change the measured fluxes of a few % was observed.

The simulated efficiency and reflectivity is shown in Fig. 19a as a function of energy for the fundamental WG mode. Note that the coupler efficiency of about 40% is referring to a single-sided emission, and that double sided efficiencies are twice as high considering the reflection symmetry in .

We note that we find a reflectivity of 10-20% for the odd mode in the relevant QD energy range (see Fig. 3), while the measured circularity provides an upper limit (see appendix E) of about 5%. There are two main aspects contributing to this difference: (i) the singly reflected light propagates up to further. Using the typical propagation losses of μm (see Fig. 7) over a distance of we find a transmission factor of 0.6. (ii) is measured close to the coupler, and thus contains contributions from non-guided modes. For the even mode, the light reflected is furthermore distributed between the two branches inside and outside of the light cone, with the former having strong propagation losses. Both aspects reduce the reflectivity relevant for the measured circularity.

In order to simulate the efficiency of coupling into the spectrometer, the far field of the couplers was extracted from the simulation. The MEEP near-field to far-field transform is used to calculate the full 6-components electromagnetic field vector at points homogeneously sampling a hemisphere of radius , much larger than any simulation feature. The resulting far-field intensity distribution in is indicated in Fig. 19b for three frequencies. The far-field was then limited to the objective collection range, and transformed back into the near field. In detail, for each far-field point, having coordinates on the hemisphere, we calculate the in-plane wavevector . We then rotate the 6-component electromagnetic field at that point from the radial to the propagation direction, resulting in , and simulating the transformation of the emitted field by the aplanatic objective into its back focal plane. The imaging of the coupler onto the spectrometer slit is then calculated by multiplying the field at each point with a 2D plane wave , and summing the vector-field plane waves from every far-field point within the microscope NA (). This procedure provides the near-field in a small angle approximation, valid for the small NA of about 0.03 of the image at the spectrometer slit, which is magnified by a factor of 31.3. This near field (see Fig. 19c) is then transmitted through the spectrometer slit (dashed lines), providing a near-field collection efficiency. The coupler efficiency, and the near-field collection efficiency, are multiplied to produce the spectrometer efficiency (see Fig. 19a), which averages to 14.4% in the simulated range. This is the efficiency with which photons emitted by a QD into the PCWG mode will enter the spectrometer, neglecting propagation loss in the WG. Losses occurring for all detected light, such as reflection loss of the optics, and detector quantum efficiencies, are not considered here as they are not influencing the measured beta factor . These are estimated to provide an additional factor of around 30% in the setup used, resulting in an estimated overall detection efficiency of QD emission around 4%.

The analysis above used the sum of the intensities of the two polarization components of the propagating field. To take into account the polarization dependence of the detection efficiency, due to the diffraction grating of the spectrometer, we separated the intensity transmitted through the spectrometer slit into and polarized components, and find 18% of the power in and 82% in polarization, at eV.

We performed an equivalent analysis for the even mode, as shown in Fig. 20. This mode shows two minima and one maximum within the PC bandgap. Close to these extrema, the mode group velocity tends to zero, thus producing a slow light regime, which gives rise to difficulties in the numerical treatment, leading to unreliable results above 1.42 eV and below 1.385 eV in our simulations. In the remaining range, the calculated coupler efficiency is slightly lower and the reflectivity is significantly higher than for the odd mode. The far-field pattern (see Fig. 20b) is more structured, and the near-field (see Fig. 20c) is wider than for the odd mode. Overall this results in a collection efficiency of the higher order mode of 10% at =1.403 eV, of which 73% is polarized.

To determine the corrected for Eq. (5), we have used a collection efficiency of 14% for the fundamental mode, calculated as the average over the simulated range, and 10% for the higher order mode.

### g.4 Beta factor and directionality

The beta factor and directionality were calculated using further FDTD simulations. We use of a simulation volume of size in , in , and in direction, as sketched in Fig. 21. A circular dipole source, , is placed in the central unit cell at , at a given position in . The source has the same time-dependence as used for the coupler simulation in appendix G.3. The beta factor and directionality of the source is calculated from the power the source radiates into the different channels. To obtain accurate results, we have to suppress reflections from the boundaries of the simulation volume. Similar to appendix G.3, we use adiabatic absorbers of thickness covering the boundaries of the simulation Oskooi et al. (2008), leaving the central 9 unit cells unperturbed. In direction, the photonic band gap results in a strong reduction of the field, and the boundaries are in the region of unstructured slab, so that absorbers with thicknesses of were found to suffice. In -direction, extending in free-space, PMLs of thickness were used. In the directions the absorption profile is quartic while in the others it is quadratic. The mirror symmetry at is exploited in the calculation.

Flux planes in the WG forward and backward direction record the frequency-resolved transmitted powers and , respectively. They are placed away from the source in , and are wide in , covering the entire photonic cladding but not the absorbers. In they extend beyond either side of the slab, to contain the mode in the WG and most of its evanescent tail, while limiting the overlap with non-guided light. A third flux plane records the power radiated out of the slab () in direction, which is doubled to account for the symmetry. These three quantities determine the beta factor and the directionality by:

(13) |

Note that refers to the total fraction of the QD emission into all guided modes and at each frequency we consider there can be up to 6 such modes.

In a single mode waveguide hosting a circular dipole source, at each location and frequency is equal to the third Stokes parameter of the mode field. In multimode waveguides is given by an average of the mode ’s weighted by their densities of states Lang et al. (2016). For comparison with experiment can be converted into a circularity value by:

(14) |

In each simulation these quantities are extracted as a function of frequency for a given source location. To probe the spatial dependence of the beta factor and circularities, simulations with the source dipole at different locations were performed. As the dipole source is moved in the direction it is moved closer to one waveguide flux plane and further from the other. For modes with propagation loss this will impact the measured directionality as the light travels further in one direction than the other. To compensate for this, all calculations were run twice, once with the dipole at and again at . The lattice symmetry ensures that the directionality must be identical at these two locations. Averaging the powers of the two simulations accounts for the effect of being closer to one flux plane or the other, removing the impact of losses on the directionality, as long as they are small over a unit cell. The effect of the losses on the beta factor remain, due to the propagation loss between the source location and the forwards/backwards waveguide flux planes. We can see in Fig. 14 that the propagation loss over this distance of about 1 μm is negligible (%) for most of the modes we consider. The only exception is the lossy branch of the even mode, which experiences significant attenuation over this distance. The energy lost by this mode is radiated out of the slab and collected by the flux plane measuring the free space emission. The mode is so lossy that it is part of the free space emission, and this is the case also in the measurements.

The beta factors inferred from these calculations are plotted in Fig. 22 versus the directionality. Each data point corresponds to a location in the WG unit cell and an emission energy. The locations are colour coded according to the map shown. The hue varies along the axis and the saturation along . Note that the breaking of symmetry in originates from the display of only half the unit cell, with the other half being mirror symmetric, mapping to and to .

To assess the sensitivity of these results on the choice of the etch parameter , simulations were performed at 5 random locations with nm instead of nm. After accounting for an overall frequency shift of (6.7 meV) the mean absolute change of was 0.024, and the change in directionality was 0.014. These small changes reflects the fact that the main effect of altering the etch parameter is a shift in frequency while hardly changing the mode profiles.

### g.5 Free-space emission

The far field intensity distributions in space of the free space emission of circular dipoles at various locations in the PCWG unit cell, taken from the simulations discussed in appendix G.4, are shown in Fig. 23 for eV. They show a significant variation with source location. Note that the missing mirror symmetry about arises from the circularly polarized dipole excitation which lacks this symmetry. Using these profiles, we determine the corresponding free space collection efficiency into the spectrometer, as shown in Fig. 24b, considering the objective NA and the spectrometer input slit as described in Sec. G.3. We can observe that the regions of low values correspond to highly asymmetric far-field profiles in Fig. 23. The average efficiency is about 29%, and varying from 18% to 33% with the location of the QD. We can consider these collection efficiencies to improve the accuracy of the beta factors determined in the experiment. We note that free space emission is collected with up to a factor of 2 higher efficiency than the WG mode emission, suggesting that the real beta factors are up to a factor of two closer to unity than reported in the main text.

To investigate the validity of the assumption of free-space emission strength independent of position, used in previous works reporting on the beta factor, we show in Fig. 24a the free space emission power as function of position, relative to bulk GaAs, for 8 energies covering the simulation range. We observe variations over two orders of magnitude, from 0.05 to 5, implying that the assumption is not justified. This points towards large systematic errors in the beta factors retrieved in previous works Lund-Hansen et al. (2008); Thyrrestrup et al. (2010); Dewhurst et al. (2010); Laucht et al. (2012a); Hoang et al. (2012); Arcari et al. (2014). Large spatial variations in the free space emission are also indicated by the simulations in Javadi et al. (2018b), which also predicts horizontal bands of reduced free space coupling as seen in Fig. 24.

For =1.4208 eV, we see large increases in the loss rate. This is attributed to the branch of the even mode in the light cone, as the energy is close to the frequency maximum in the even mode, where its group velocity tends to zero. For a low group velocity the emission is enhanced by the slow light, and the propagation length is reduced, resulting in an emission into free space close to the emitter. Indeed, the pattern resembles the electric field intensity pattern of the even (higher order) seen in Fig. 15.

### g.6 Mode Separation

The FDTD calculations provide the beta factor and directionality. The values observed in the experiment will be slightly modified from these values by the different collection efficiencies of coupler and free space emission. As the collection efficiencies of the even and odd mode are different, adjusting for these efficiencies requires separation of the power between the two modes. To achieve this separation, the electrical field profiles ) and group indices extracted from the loss FDTD simulations (see appendix G.2) are used. The group indices are calculated using the spatially integrated Poynting vector Chen et al. (2010). These are used to calculate the density of states (DOS) into the forward and backward directions at each spatial location using the Purcell factor expressions presented in Rao and Hughes (2007), adapted to separate the intensity in a given direction. As we are only concerned with ratios between modes at a common frequency, these expressions can be simplified to

(15) |

Using the fields associated with each mode, respectively, allows us to calculate the fraction of the power emitted into the forward even mode, for a given location and frequency. We neglect the branch of the even mode in the light cone, as the power in it is lost while propagating to the coupler. Using the resulting fractions the collection efficiencies of the couplers can be corrected while accounting for the different efficiencies for the two modes.

## Appendix H Correction of simulations

In this section we describe how the raw data presented in Sec. G.4 have been corrected in order to be comparable to the experimental data. The forward and the backward propagating fluxes are firstly separated into the fundamental and higher order modes projections as calculated in Sec. G.6. The fluxes of each mode are then multiplied by the corresponding collection efficiencies due to the couplers, the NA of the MO and the input slit of the imaging spectrometer. We recall that the collection efficiencies are averaged over the spectral range covered by the simulations. The resulting fluxes from each mode are then separated into transverse and longitudinal polarisation components using the analysis presented in Sec. G.3. Finally, we include the polarisation dependent efficiency and of the spectrometer grating for transverse and longitudinal polarisations respectively. The polarisation resolved fluxes are multiplied by the coefficients and , and the resulting fluxes are then summed up, separately for each mode, to obtain the final corrected fluxes. and are determined by measuring the grating polarisation efficiency ratio , with the normalisation , such that for polarization independent detection we have .

For the free space emission we first multiply the out-of-plane flux by the position dependent collection efficiency, as described in Sec. G.5. In this case the collection efficiency is calculated at eV, which is approximately the central energy of the measured QDs. We calculated the energy dependence at one location and we found it to be negligible. In order to include the polarisation efficiency of the grating spectrometer, we consider the free space emission to be circularly polarised.

## References

- Rao and Hughes (2007) V. S. C. M. Rao and S. Hughes, Phys. Rev. B 75, 205437 (2007).
- Lecamp et al. (2007) G. Lecamp, P. Lalanne, and J. P. Hugonin, Phys. Rev. Lett. 99, 023902 (2007).
- Rodríguez-Fortuo et al. (2013) F. J. Rodríguez-Fortuo, G. Marino, P. Ginzburg, D. O’Connor, A. Martínez, G. A. Wurtz, and A. V. Zayats, Science 340, 328 (2013).
- Junge et al. (2013) C. Junge, D. O’Shea, J. Volz, and A. Rauschenbeutel, Phys. Rev. Lett. 110, 213604 (2013).
- Mitsch et al. (2014) R. Mitsch, C. Sayrin, B. Albrecht, P. Schneeweiss, and A. Rauschenbeutel, Nat. Commun. 5, 5713 (2014).
- Söllner et al. (2015) I. Söllner, S. Mahmoodian, S. L. Hansen, L. Midolo, A. Javadi, G. Kiršanske, T. Pregnolato, H. El-Ella, E. H. Lee, J. D. Song, S. Stobbe, and P. Lodahl, Nat. Nano. 10, 775 (2015).
- Young et al. (2015) A. Young, A. Thijssen, D. Beggs, P.Androvitsaneas, L. Kuipers, J. Rarity, S. Hughes, and R. Oulton, Phys. Rev. Lett. 115, 153901 (2015).
- Lodahl et al. (2017) P. Lodahl, S. Mahmoodian, S. Stobbe, A. Rauschenbeutel, P. Schneeweiss, J. Volz, H. Pichler, and P. Zoller, Nature 541, 473 (2017).
- Coles et al. (2016) R. Coles, D. Price, J. Dixon, B. Royall, E. Clarke, P. Kok, M. Skolnick, A. Fox, and M. Makhonin, Nat. Commun. 7, 11183 (2016).
- Coles et al. (2017) R. J. Coles, D. M. Price, B. Royall, E. Clarke, M. S. Skolnick, A. M. Fox, and M. N. Makhonin, Phys. Rev. B 95, 121401 (2017).
- Lang et al. (2015) B. Lang, D. M. Beggs, A. B. Young, J. G. Rarity, and R. Oulton, Phys. Rev. A 92, 063819 (2015).
- Javadi et al. (2018a) A. Javadi, D. Ding, M. H. Appel, S. M. M. C. Löbl, I. Söllner, R. Schott, C. Papon, T. Pregnolato, S. Stobbe, L. Midolo, T. Schröder, A. D. Wieck, A. Ludwig, R. J. Warburton, and P. Lodahl, Nature Nanotechnology 13, 398 (2018a).
- Kim et al. (2018) J.-H. Kim, S. Aghaeimeibodi, C. J. K. Richardson, R. P. Leavitt, and E. Waks, Arxiv (2018).
- Lodahl et al. (2015) P. Lodahl, S. Mahmoodian, and S. Stobbe, Review of Modern Physics 87, 347 (2015).
- Mahmoodian et al. (2016) S. Mahmoodian, P. Lodahl, and A. S. Sørensen, Phys. Rev. Lett. 117, 240501 (2016).
- Lund-Hansen et al. (2008) T. Lund-Hansen, S. Stobbe, B. Julsgaard, H. Thyrrestrup, T. Sünner, M. Kamp, A. Forchel, and P. Lodahl, Phys. Rev. Lett. 101, 113903 (2008).
- Thyrrestrup et al. (2010) H. Thyrrestrup, L. Sapienza, and P. Lodahl, Appl. Phys. Lett. 96, 231106 (2010).
- Dewhurst et al. (2010) S. J. Dewhurst, D. Granados, D. J. P. Ellis, A. J. Bennett, R. B. Patel, I. Farrer, D. Anderson, G. A. C. Jones, D. A. Ritchie, and A. J. Shields, Appl. Phys. Lett. 96, 031109 (2010).
- Laucht et al. (2012a) A. Laucht, S. Pütz, T. Günthner, N. Hauke, R. Saive, S. Frédérick, M. Bichler, M.-C. Amann, A. Holleitner, M. Kaniber, and J. J. Finley, Phys. Rev. X 2, 011014 (2012a).
- Hoang et al. (2012) T. B. Hoang, J. Beetz, L. Midolo, M. Skacel, M. Lermer, M. Kamp, S. Höfling, L. Balet, N. Chauvin, and A. Fiore, Appl. Phys. Lett. 100, 061122 (2012).
- Arcari et al. (2014) M. Arcari, I. Söllner, A. Javadi, S. L. Hansen, S. Mahmoodian, J. Liu, H. Thyrrestrup, E. Lee, J. Song, S. Stobbe, and P. Lodahl, Phys. Rev. Lett. 113, 093603 (2014).
- Kurzmann et al. (2016) A. Kurzmann, A. Ludwig, A. D. Wieck, A. Lorke, and M. Geller, Nano Lett. 16, 3367 (2016).
- Langbein et al. (2004a) W. Langbein, P. Borri, U. Woggon, V. Stavarache, D. Reuter, and A. D. Wieck, Phys. Rev. B 70, 033301 (2004a).
- Javadi et al. (2018b) A. Javadi, S. Mahmoodian, I. Söllner, and P. Lodahl, Journal of the Optical Society of America B: Optical Physics 35, 514 (2018b).
- Rigal et al. (2018) B. Rigal, B. Dwir, A. Rudra, I. Kulkova, A. Lyasota, and E. Kapon, Appl. Phys. Lett. 112, 051101 (2018).
- Tighineanu et al. (2018) P. Tighineanu, C. L. Dreeßen, C. Flindt, P. Lodahl, and A. S. Sørensen, Phys. Rev. Lett. 120, 257401 (2018).
- Langbein et al. (2004b) W. Langbein, P. Borri, U. Woggon, V. Stavarache, D. Reuter, and A. D. Wieck, Phys. Rev. B 69, 161301(R) (2004b).
- Johnson et al. (1999) S. Johnson, S.Fan, P. Villeneuve, J. D. Joannopoulos, and L. Kolodziejski, Phys. Rev. B 60, 5751 (1999).
- Faraon et al. (2008) A. Faraon, I. Fushman, D. Englund, N. Stoltz, P. Petroff, and J. Vučković, Opt. Express 16, 12154 (2008).
- Sukhorukov et al. (2009) A. A. Sukhorukov, S. Ha, I. V. Shadrivov, D. A. Powell, and Y. S. Kivshar, Optics Express 17, 3716 (2009).
- Gersen et al. (2005) H. Gersen, T. Karle, R. Engelen, W. Bogaerts, J. Korterik, N. van Hulst, T. Krauss, and L. Kuipers, Physical Review Letters 94, 073903 (2005).
- Laucht et al. (2012b) A. Laucht, T. Günthner, S. Pütz, R. Saive, S. Frédérick, N. Hauke, M. Bichler, M.-C. Amann, A. W. Holleitner, M. Kaniber, and J. J. Finley, Journal of Applied Physics 112, 093520 (2012b).
- Langbein et al. (1999) W. Langbein, J. M. Hvam, and R. Zimmermann, Phys. Rev. Lett. 82, 1040 (1999).
- Langbein and Hvam (2002) W. Langbein and J. M. Hvam, Phys. Rev. Lett. 88, 047401 (2002).
- Thomas et al. (2007) N. L. Thomas, R. Houdré, M. V. Kotlyar, D. O’Brien, and T. F. Krauss, J. Opt. Soc. Am. B 24, 2964 (2007).
- Gehrsitz et al. (2000) S. Gehrsitz, F. K. Reinhart, C. Gourgon, N. Herres, A. Vonlanthen, and H. Sigga, Journal of Applied Physics 87, 7825 (2000).
- Lukeš (1972) F. Lukeš, Surf. Sci. 30, 91 (1972).
- DeSalvo et al. (1996) G. C. DeSalvo, C. A. Bozada, J. L. Ebel, D. C. Look, J. P. Barrette, C. L. A. Cerny, R. W. Dettmer, J. K. Gillespie, C. K. Havasy, T. J. Jenkins, K. Nakano, C. I. Pettiford, T. K. Quach, J. S. Sewell, and G. D. Via, Journal of The Electrochemical Society 143, 3652 (1996).
- Bayer et al. (2002) M. Bayer, G. Ortner, O. Stern, A. Kuther, A. A. Gorbunov, A. Forchel, P. Hawrylak, S. Fafard, K. Hinzer, T. L. Reinecke, S. N. Walck, J. P. Reithmaier, F. Klopf, and F. Schäfer, Phys. Rev. B 65, 195315 (2002).
- Benny et al. (2012) Y. Benny, Y. Kodriano, E. Poem, D. Gershoni, T. A. Truong, and P. M. Petroff, Phys. Rev. B 86, 085306 (2012).
- Oskooi et al. (2010) A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. Joannopoulos, and S. G. Johnson, Comput. Phys. Commun. 181, 687 (2010).
- Oskooi et al. (2008) A. Oskooi, L. Zhang, Y. Avniel, and S. Johnson, Opt. Express 16, 11376 (2008).
- Lang et al. (2016) B. Lang, D. M. Beggs, and R. Oulton, Phil. Trans. R. Soc. A 374, 2075 (2016).
- Akahane et al. (2005) Y. Akahane, T. Asano, B.-S. Song, and S. Noda, Opt. Express 13, 1202 (2005).
- Johnson and Joannopoulos (2001) S. G. Johnson and J. D. Joannopoulos, Opt. Express 8, 173 (2001).
- Chen et al. (2010) Y. Chen, T. R. Nielsen, N. Gregersen, P. Lodahl, and J. Mørk, Phys. Rev. B 81, 125431 (2010).