Efficient optical pumping using hyperfine levels in {}^{145}Nd{}^{3+}:Y{}_{2}SiO{}_{5} and its application to optical storage

Efficient optical pumping using hyperfine levels in Nd:YSiO and its application to optical storage

Emmanuel Zambrini Cruzeiro    Alexey Tiranov    Jonathan Lavoie Present address: Department of Physics and Oregon Center for Optical Molecular & Quantum Science, University of Oregon, Eugene, OR 97403, USA    Alban Ferrier    Philippe Goldner    Nicolas Gisin    Mikael Afzelius mikael.afzelius@unige.ch Groupe de Physique Appliquée, Université de Genève, CH-1211 Genève, Switzerland PSL Research University, Chimie ParisTech, CNRS, Institut de Recherche de Chimie Paris, 75005 Paris, France Sorbonne Universités, UPMC Univ Paris 06, Paris 75005, France
July 7, 2019

Efficient optical pumping is an important tool for state initialization in quantum technologies, such as optical quantum memories. In crystals doped with Kramers rare-earth ions, such as erbium and neodymium, efficient optical pumping is challenging due to the relatively short population lifetimes of the electronic Zeeman levels, of the order of 100 ms at around 4 K. In this article we show that optical pumping of the hyperfine levels in isotopically enriched Nd:YSiO crystals is more efficient, owing to the longer population relaxation times of hyperfine levels. By optically cycling the population many times through the excited state a nuclear-spin flip can be forced in the ground-state hyperfine manifold, in which case the population is trapped for several seconds before relaxing back to the pumped hyperfine level. To demonstrate the effectiveness of this approach in applications we perform an atomic frequency comb memory experiment with 33% storage efficiency in Nd:YSiO, which is on a par with results obtained in non-Kramers ions, e.g. europium and praseodymium, where optical pumping is generally efficient due to the quenched electronic spin.

I Introduction

In rare-earth-ion doped crystals, one often uses frequency-resolved optical pumping (i.e. spectral hole burning) as a necessary preparation step in various applications, such as optical quantum memories Tittel et al. (2010); Afzelius et al. (2009); Hedges et al. (2010); Jobez et al. (2016), spectrum radio-frequency analysers Babbitt et al. (2007); Linget et al. (2015), and narrow spectral filters for biological imaging Zhang et al. (2012); McAuslan et al. (2012). Common to all these techniques is the need for creating narrow and deep spectral holes in the inhomogeneously broadened optical transition. The most efficient way for realizing deep spectral holes is to optically pump ions into a long-lived ground state, which could be electronic spin or nuclear spin states. Generally it is observed that the efficiency of the optical pumping, i.e. the depth of the spectral hole, strongly depends on the ground-state population lifetime with respect to the radiative lifetime Lauritzen et al. (2008); Afzelius et al. (2010a); Saglamyurek et al. (2015); Zhong et al. (2017).

In non-Kramers ions, such as europium, praseodymium or thulium, the ground state is an electronic singlet in low-symmetry doping sites, such that only nuclear interactions (eg. Zeeman and quadrupole) exist in the ground state Macfarlane and Shelby (1987); Macfarlane (2002). As a result both spin-spin interaction and spin-lattice relaxation rates are low, with lifetimes of many seconds or hours Macfarlane (2002), and efficient optical pumping can be achieved in a wide range of experimental parameters (in terms of doping concentration, temperature and magnetic field).

In Kramers ions such as erbium, neodymium or ytterbium, the ground state is an electronic Zeeman doublet in low-symmetry doping sites, with a magnetic moment in the range of 1–15 , where  GHz/T is the Bohr magneton. The large moment results in strong dipole-dipole interactions between Kramers ions and faster spin-lattice relaxation rates. In terms of applications, however, Kramers ions are interesting as the large Zeeman and hyperfine splittings allow large bandwidth quantum memories Saglamyurek et al. (2011); Clausen et al. (2011); Zhong et al. (2017) and could be interfaced with superconducting qubits working in the 1-10 GHz regime Staudt et al. (2012); Probst et al. (2013).

Recently we presented a detailed study of the spectral hole lifetime in naturally doped Nd:YSiO Cruzeiro et al. (2017), as a function of magnetic field strength and direction, temperature (between 3 and 5.5 K) and Nd doping concentration. It was found that the spectral hole dynamics was dominated by the population relaxation between the electronic Zeeman states of even Nd isotopes having zero nuclear spin . For doping concentrations required for practical applications (a few tens of ppm) the longest measured lifetime was around 160 ms. This relatively short spectral hole lifetime results in residual absorption backgrounds of about 6–7% of the peak absorption, which has been the main factor limiting the efficiency in several quantum storage experiments we have carried out in naturally doped Nd:YSiO Clausen et al. (2011); Usmani et al. (2012); Fröwis et al. (2017). Such background absorptions would not allow very efficient (50% or above) quantum memories based on cavity enhancement Moiseev et al. (2010); Afzelius and Simon (2010), as demonstrated in praseodymium and europium doped materials Sabooni et al. (2013); Jobez et al. (2014). Also, there is an interest in using the 883 nm resonance in Nd:YSiO crystal for biological imaging Walther et al. (2017), but spectral filtering with high dynamic range requires a lower residual absorption. Hence, further application of Nd-doped materials in these research areas would greatly benefit from longer spectral hole lifetimes.

In naturally doped Er:YSiO similar spectral hole lifetimes of up to 130 ms have been achieved using isotopes with , at very low magnetic fields and around 3 K Hastings-Simon et al. (2008). Longer lifetimes can be achieved both in Er:YSiO crystals Probst et al. (2013) and Er-doped fibers Saglamyurek et al. (2015), but at sub-K temperatures. At temperatures above 1 K, Rančić and co-workers recently demonstrated long spectral hole lifetimes of up to 60 seconds in a isotopically enriched Er:YSiO crystal Rančić et al. (2017), where Er has a nuclear spin . In their approach spin-lattice relaxation between the electronic Zeeman states is suppressed by applying a high magnetic field of 3 T or more, in which case the hole lifetime is dominated by relaxation between nuclear states in the lowest electronic Zeeman state. A similar suppression of the electronic spin-lattice relaxation in Nd-doped samples would require much higher magnetic fields, or much lower temperatures, due to four times lower magnetic moment of Nd ions.

As an alternative approach we explore spectral hole lifetimes in moderate magnetic fields and a temperature of about 3 K in an isotopically enriched Nd:YSiO crystal, detailed in Sec. II, where Nd also has a nuclear spin =7/2. In Sec. III we show that the spectral hole lifetime depends strongly on the duration of the hole burning pulse. By burning for durations of around one second, the spectral hole is dominated by a slow decay process with a lifetime that can reach 4 seconds. In Sec. IV we develop a simple model to explain our data, where the basic idea is that optical pumping for long durations increases the probability of flipping the nuclear spin projection in the ground state. Once the nuclear spin has flipped through the optical excitation, the spin flip-flop and spin-lattice relaxation rates are reduced, as these rates depend on nuclear spin mixing induced by the non-secular part of the hyperfine interaction . In Sec. V it is shown that the efficiency of the optical pumping is enhanced as a result of the longer hole lifetimes, with an absorption background of 1-2% of the peak absorption. By applying these results to an optical storage experiment based on an atomic frequency comb memory, we reach storage efficiencies of up to 33%, which is mainly limited by the peak absorption rather than the residual background absorption.

Ii Experimental Details

Figure 1: The optical inhomogeneous absorption profile on the IF transition, at a temperature of 3 K and with no applied magnetic field. The measured peak optical depth is 1.26 (0.44) for the light polarization along the () axis, resulting in an absorption coefficient = 1.05 cm ( = 0.37 cm). The inhomogeneous lineshape is Gaussian with a FWHM linewidth of 7.7 GHz. In the inset we show the holeburning spectrum showing the central hole and the associated side holes due to the hyperfine splitting in the optically excited state F.

Yttrium orthosilicate YSiO crystal belongs to the crystallographic group . It is known as a host material that provides excellent coherence properties due to its low nuclear spin density. Neodymium ions can substitute yttrium ions in two crystallographic sites, both of low symmetry Maksimov et al. (1970). The crystal used in this work was doped with 0.001 at.% concentration of Nd ions, with an isotopic purity of about 90% Maier-Flaig (2013). We note that crystals from the same boule have been used in electron spin measurements presented in Refs Maier-Flaig (2013); Wolfowicz et al. (2015). Our sample was cut with faces perpendicular to the , , and optical extinction axes Li et al. (1992), where light was propagating along the axis (of length 12 mm). The states of interest in this work are the ground I and excited F states of site 1. In the site symmetry, the ground and excited states split into 5 and 2 Kramers doublets, respectively. All measurements presented here were carried out on the transition between the Kramers doublets of lowest energy, at a wavelength of 883.0 nm (in vacuum) for site 1 Beach et al. (1990).

The optical measurements were performed with an external cavity diode laser. Large bandwidth frequency scans (in the GHz range) were performed by scanning a piezo element in the cavity. An acousto-optical modulator (AOM) in double pass configuration provided fine control of the amplitude and frequency of the light (120 MHz bandwidth). The AOM was used to create the optical pulses for the spectral holeburning measurements, as well as for preparing the atomic frequency comb in the light storage experiment.

In Figure 1 we show the optical absorption spectrum recorded at 3 K and zero applied magnetic field. The strongest absorption is obtained with the linear light polarization along the axis. The absorption profile follows a smooth distribution, without any visible substructure due to the hyperfine splittings of the ground and excited states. To reveal transitions hidden within a large inhomogeneous absorption profile one can employ spectral hole burning (SHB) techniques Macfarlane and Shelby (1987). A typical SHB measurement consists of a burn pulse at a fixed frequency and a time-delayed probe pulse whose frequency is scanned around the burn pulse frequency. The spectral hole appears at the burn frequency, while side holes appear at frequencies given by any splittings in the optically excited state Macfarlane and Shelby (1987). In the inset of Figure 1 we show the SHB spectrum in Nd:YSiO with zero applied magnetic field. A set of periodic side holes can indeed be observed, due to the hyperfine levels in the excited state. Unfortunately the large number of expected side holes due to the hyperfine levels does not allow any quantitative analysis of the hole structure. However, the exact SHB spectrum is not of interest here, as we will principally investigate the time-resolved dynamics of the central hole in the SHB spectrum.

Iii Spectral hole decay measurements

The SHB mechanics of interest here is due to trapping of population in ground state hyperfine levels that are not excited by the burn pulse in the SHB sequence. In this way the observed hole decay curve contains information about different relaxation rates within the ground state hyperfine levels. This regime is achieved by optically burning the hole for a duration much longer than the excited state lifetime, which for F is 225 Beach et al. (1990). In the experiments presented here the burn pulse duration varied between 10 ms and 950 ms. To measure the decay of the hole, we varied the delay between the burn and probe pulses and recorded the area of the central spectral hole for each delay. By using the hole area the measurement is independent of spectral diffusion. A range of SHB decay measurements were made with different burn pulse durations and applied magnetic fields. For all measurements, however, the magnetic field orientation was along the axis and the temperature was 3 K.

All SHB decay measurements featured the same overall trend of a fast decay in the beginning, followed by a much slower decay of the hole. A striking observation was that the relative amplitudes of the fast and slow decays depended on the duration of the burn pulse , while the time constants of the fast and slow decays were independent of the burn pulse duration. In Figure 2(a) we show two typical SHB decay curves for a shorter (100 ms) and longer (950 ms) burn pulse, at a magnetic field of 900 mT. In this particular case the fast decay constant was 75 ms, while the slow decay constant was 1.72 s. It is also clearly seen that the fast decay dominates the decay for a short burn pulse, while the slow decay dominates the decay for the long burn pulse.

Figure 2: (a) The spectral hole amplitude as a function of the delay between the burn and probe pulses, for two burn durations of = 100 ms (circles) and 950 ms (squares), respectively. The vertical axis shows the logarithm of the optical depth in the spectral hole. Each curve has been shifted to zero at the shortest delay. The solid lines are fitted curves using the simple model described in the text. (b) Relative amplitudes of the fast (circles) and slow (squares) decay constants, as a function of the burn duration (see text for details).

To investigate the dependence of the relative amplitudes of the fast and slow decays on the burn duration, we recorded and analysed an entire set of decay curves for in the range of 10 to 950 ms. Each set was analysed using only two time constants that were common to all curves, while the relative amplitudes were fitted individually. Specifically, the decay of the spectral hole depth was modelled as . Here and are the amplitudes, and and the time constants, of the fast and slow decays, respectively. The parameters and were fitted for each burn duration, while and were fitted globally for the entire set of burn durations for a given magnetic field strength. We emphasize that all sets of decay curves could be analysed with only two time constants, for all magnetic fields.

In Figure 2(b) we show the relative amplitudes of the fast and slow decays, which we define as and , respectively, as a function of the burn duration for a field of 900 mT. By pumping for sufficiently long, one is able to force the hole to decay predominantly through a slower relaxation mechanism. The crossing point in this particular case is for a burn duration of about 90 ms. In the following section we will present a simple relaxation model to explain the observed SHB curves.

Iv Hyperfine relaxation model

The spectral hole is created by a re-distribution of population among the hyperfine levels in the ground state of Nd, with respect to the thermal distribution, due to the cycling of population through the optically excited state. At most magnetic fields used here, the ground state Zeeman splitting between the branches is larger than the optical broadening, such that the corresponding states can be resolved. In our recent study of Nd ions having no nuclear spin Cruzeiro et al. (2017), the spectral hole decay did not depend on the initial state for fields in this range, so in the following we assume that we pump from for the sake of simplicity, as shown see Fig. 3. The split of the hyperfine states within a given branch, on the other hand, is not resolved optically, as exemplified by the smooth zero-field absorption spectrum shown in Fig. 1. The spectral hole is thus an average of optical pumping of all possible initial states . In addition, in this regime of magnetic fields the hyperfine interaction is a weak perturbation to the Zeeman interaction. Then the separable states describe, to a good approximation, the hyperfine levels in the electronic ground state (see Fig. 3). However, as we will see, a small mixing of states due to the hyperfine interaction is essential in order to understand the different spin relaxation rates.

As population is optically cycled through the excited state, ions can decay to different hyperfine levels. If the ions end up in the opposite state having the same , then it can decay with a relaxation process that does not involve any change in the nuclear spin projection . We denote this relaxation rate , as shown in Fig. 3. For the field of 900 mT used in the experiment shown in Fig. 2, the dominant relaxation process is the direct spin lattice relaxation Cruzeiro et al. (2017). For the spin lattice relaxation (SLR) process this relaxation path is expected to have the highest rate (shortest lifetime)  Larson and Jeffries (1966a), due to the low spin mixing. The hyperfine Hamiltonian , does however open up other relaxation paths, owing to non-secular mixing terms such as and . We denote these relaxation rates by and , which are expected to be significantly slower than  Larson and Jeffries (1966a). Other, pure nuclear relaxation paths, for which , are assumed to be even slower Rančić et al. (2017) and are not considered here.

Using this simple model we can make sense of the spectral hole decay curves shown in Fig. 2, by additionally assuming that the optical decay preferably preserves the nuclear spin projection , as shown by the solid line in Fig. 3. A shorter burn pulse would then preferentially optically pump ions into the opposite state with identical projection, yielding a fast hole relaxation given by . By increasing the burn duration we can increase the chance of optically decaying into neighbouring states, as shown by the dashed lines in Fig. 3, which would yield slower hole decay rates of the order of and/or . In principle, we should then observe two additional rates, while experimentally one slow decay constant is observed. However, in theory the rates and also have dependence, which is also different for the two rates, as discussed by Larson and Jeffries Larson and Jeffries (1966a). In our case we believe that either we measure an effective average rate, which would correspond to the slow decay constant , or one of the rates dominate over the other for the particular field orientation used in this experiment. Further experimental measurements and a detailed theoretical modelling would be necessary to decide upon this question. Yet, the overall behaviour is consistent with this simple model.

Figure 3: Simplified energy level diagram of Nd  for understanding the simple relaxation model presented in Sec. IV. The optical excitation is shown as the red, solid line. The most probable decay path from the excited state is shown as a black solid line, while the less probable paths that change nuclear spin projection are shown as black dashed lines. The fast and slow spin relaxation rates are shown as blue solid and dashed lines, respectively. Three different spin relaxation rates , and are defined. For most of the magnetic fields studied in this article, the ground state splits into two electronic branches (), with eight possible nuclear spin projections for each branch (). For simplicity only three nuclear spin projections are shown. The hyperfine structure of the excited state is not discussed here, hence it is shown as a fuzzy region.

To better understand the relaxation mechanisms governing the fast and slow decays, we measured the hole decay curve for fields in the range of 10 mT to 1.6 Tesla. But before discussing these results, we briefly recall the main findings of Ref. Cruzeiro et al. (2017) where we studied the decay of spectral holes due to relaxation between states for ions with no nuclear spin . At low fields (around 0.5 Tesla) the decay was given by spin-spin relaxation, with a linear increase of the spectral hole lifetime as a function of applied magnetic field . At higher fields the direct SLR process caused a decrease of the lifetime with the well-known scaling  Orbach (1961); Scott and Jeffries (1962); Larson and Jeffries (1966b, a).

In Figure 4 the fast and slow lifetimes are shown as a function of the applied field. As seen the lifetimes generally increase at low fields, reach a maximum around 0.4 Tesla, and then start to decrease for higher fields, similarly to the case of ions with . In Ref. Cruzeiro et al. (2017) it was found that a SLR model including the direct, Raman and Orbach processes explained well the high-field data for ions with . In Figure 4 the blue solid line shows the lifetime predicted by this model, using the fitted parameters Cruzeiro et al. (2017) for this particular magnetic field angle ( axis). The model agrees rather well with the high-field data of the fast decay component. This is supporting our hypothesis that the fast decay is related to the relaxation process that only changes the electronic spin project ( and ), denoted by the rate in Fig. 3. The small difference we believe is due to the fact that the should depend slightly on the value, as discussed by Larson and Jeffries Larson and Jeffries (1966a), while our rate is an average over all nuclear projections . The scaling with magnetic field appears to be somewhat slower than the expected , possibly also due to the averaging over the states. The experimental measurement uncertainty for these short lifetimes does not allow to be certain, however, that there is a real change in the scaling.

Figure 4: The fast (circles) and slow (squares) lifetimes as a function of applied magnetic field along the axis. The blue, solid line shows the SLR model of Ref. Cruzeiro et al. (2017), using the model parameters fitted in that work. The blue, dashed line shows the complete model of Ref. Cruzeiro et al. (2017), consisting of the SLR and flip-flop models, where the latter has been scaled down by a factor of 24 as explained in Sec. IV. The red, solid line is a fitted scaling law .

The scaling of the slower , rates, on the other hand, is expected to be significantly different Larson and Jeffries (1966a); Abragam and Bleaney (1970). Indeed, as the Zeeman energy increases with the magnetic field, the spin mixing due to the non-secular terms of the hyperfine Hamiltonian will decrease. As discussed by Abragam and Bleaney Abragam and Bleaney (1970), this causes the relaxation rates , to decrease, with respect to the relaxation rate, with a factor . Hence, qualitatively one would expect a scaling of the , rates. In Fig. 4, we show a fit of the high-field data of the slow lifetime component to the formula , which yields an exponent . As before, the averaging over the states might be responsible for the difference in scaling. Qualitatively, however, the observed scaling supports our hypothesis that the slow decay is related to the SLR rates , involving nuclear spin flips.

Finally, we comment on the lifetimes measured below 0.5 Tesla, where spin flip-flops limit the spectral hole lifetime for ions with  Cruzeiro et al. (2017). We first note that for fields down to about 0.15 Tesla, the hyperfine Hamiltonian is still a perturbation to the Zeeman Hamiltonian, in which case one would expect the spin flip-flop process to mainly drive the transition shown in Fig. 3. This is because the spin flip-flop process is driven by dipole-dipole interactions between the electronic spins, and changes in the nuclear spin projection are due to the hyperfine Hamiltonian, similarly as for the SLR process. Generally we then expect a similar difference in the and , rates, respectively. This explains the difference in the measured fast and slow lifetimes also in the low-field region, see Fig. 4. In the following we will concentrate on the lifetime of the fast decay.

As compared to the spin flip-flop model that was fitted to the data in Ref. Cruzeiro et al. (2017), here we expect a much slower flip-flop rate as there are much less spins in the opposite state with the same projection (a reduction of 1/8 due to the nuclear states). In addition our sample is less doped (10 ppm compared to 30 ppm), which reduces the dipole-dipole interaction. In Ref. Cruzeiro et al. (2017) we found approximately a linear dependence of the rate on the Nd concentration. So, based on these simple arguments we expect the flip-flop rate to be reduced by factor of 24, qualitatively. This is also supported by the experimental data, as the maximum fast lifetime is 850 ms in this isotopically enriched 10 ppm Nd:YSiO crystal, as opposed to 80 ms for the 30 ppm Nd:YSiO crystal Cruzeiro et al. (2017).

In Fig. 4 we overlay the complete lifetime model of Cruzeiro et al. (2017), shown as a blue dashed line, including the SLR and flip-flop model, where the flip-flop rate has been scaled down by a factor of 24. The excellent agreement with the scaled-down flip-flop model is rather surprising, given the simplicity of the scaling. Admittedly, however, a more detailed study of the angular dependence of the slow and fast lifetimes would be necessary for a more comprehensive understanding of the flip-flop rate between hyperfine levels. It appears clear, however, that the flip-flop process is a limiting factor also for Nd, even at concentrations as low as 10 ppm.

V Application to an optical memory

In this section we discuss a light storage experiment performed in Nd:YSiO, based on the atomic frequency comb (AFC) scheme Afzelius et al. (2009). The AFC scheme is based on the creation of a periodic comb structure in the absorption profile, using a periodic SHB procedure. When an optical input pulse is absorbed by the AFC, it will generate an output pulse after a time (the AFC echo), where is the frequency spacing in the comb. The AFC echo is a delay-line optical memory, which can be used to store quantum states of light Saglamyurek et al. (2011); Clausen et al. (2011); Usmani et al. (2012); Zhou et al. (2015). An on-demand memory can be implemented using a three-level system and optical control pulses Afzelius et al. (2009, 2010b); Gündoğan et al. (2013), but here we restrict ourselves to the basic AFC echo.

The efficiency of the AFC echo depends on the optical depth of the comb peaks, the finesse of the comb and any residual optical depth in between the comb peaks. The efficiency can be calculated using the formula Afzelius et al. (2009) , where is the average optical depth and is an intrinsic dephasing factor due to the comb finesse. The optimal peak shape is squarish Bonarota et al. (2010), in which case and . For a given optical depth , there is an optimal comb finesse . To achieve high efficiencies, the challenge is to make deep spectral holes that preserve the peak optical depth , while reducing such that it does not limit the obtained efficiency significantly.

Figure 5: Optical pumping efficiency measured by the ratio , as a function of the burning time . For the longest burn time the is . The inset shows the AFC structure using a 900 ms burn. The details on the AFC preparation method are described in Ref. Clausen et al. (2011).
Figure 6: The AFC echo efficiency (circles) for different optical depths . The solid (red) line shows the model prediction based on the measured ratio, with the shaded region showing the uncertainty due to the error estimation in . The model predictions for the memory reported in Clausen et al. (2011) () and the perfect memory () are shown as dot-dashed and dashed lines, respectively.

As shown in Fig. 4, the maximum lifetimes are reached around a field of 0.4-0.5 T, both for the fast and slow decays. However, we must also consider that the inhomogeneous optical line, as shown in Fig. 1, splits into two lines for some field due to the ground state splitting, thereby reducing the absorption. For a field along the axis, the ground state -factor is , such that the split corresponds to the optical inhomogeneous linewidth of 7.7 GHz for a field of 0.37 T. To maximize the absorption we set the field to 0.15 T, where we still reach a spectral hole lifetime of s for the slow decay component (see Fig. 4). The light also passes through the crystal 4 times, in order to reach a maximum optical depth of .

In Figure 5 the ratio of the residual to peak optical depth is plotted, as a function of the burning time. The comb was created by scanning the laser while modulating the intensity using the double-pass AOM, the same technique used in Clausen et al. (2011). As expected, the ratio decreases with increasing burn time, while also preserving the peak optical depth, which confirms that the optical pumping efficiency is increasing as the spectral hole decays with a slower rate. For a comb burn time of 900 ms, the ratio is , significantly lower than the 6-7% obtained in the reference experiment for Nd:YSiO crystals in terms of AFC storage efficiency Clausen et al. (2011).

The AFC echo efficiency was then measured as a function of the peak optical depth, using a comb burn time of 1000 ms. The overall optical depth was varied by turning the linear light polarization before the crystal, and using the strong dependence of the optical depth on the polarization Clausen et al. (2012). As shown in Fig. 6, the efficiency reaches a maximum of for the maximum optical depth. The measured efficiency is also compared to the theoretical model using the ratio , showing an excellent agreement. For a longer crystal the model predicts a maximum efficiency of about 43%. The maximum theoretical efficiency for a ratio of is 26%, on the other hand, which shows the improvement achieved in the Nd:YSiO crystal.

The maximum efficiency is comparable to the efficiencies reached in praseodymium and europium doped YSiO crystals without using cavities, which are 35% Amari et al. (2010) and 32% Jobez et al. (2016), respectively. Higher efficiencies (50-55%) have been reached in both materials using cavities Jobez et al. (2014); Sabooni et al. (2013). An analysis using the equations in Ref. Jobez et al. (2014) shows that the loss due to the residual absorption would allow cavity-enhanced efficiencies in Nd:YSiO in the range of 71-89%, for ratios within the experimental error . The uncertainty in this estimation is quite large due to the difficulty to accurately measure the ratio. Also, other losses might also reduce the efficiency Jobez et al. (2014); Sabooni et al. (2013). The main point is that with spectral hole lifetimes of several seconds, using the hyperfine structure in Nd, efficient optical pumping is possible and light storage efficiencies can reach those reached in non-Kramers ions, where optical pumping is known to be efficient due to long-lived nuclear spin states.

Vi Conclusions and Outlook

We have presented spectral hole burning decay measurements in an isotopically enriched Nd:YSiO crystal. It was observed that the spectral hole decays with two characteristic time scales, a fast and a slow one, which can be explained by a model invoking relaxation paths involving and transitions, respectively. The magnetic field dependence of the fast and slow decay times was found to follow the same trend as in naturally doped Nd:YSiO crystals, but with a reduced contribution from spin flip-flops due to the many hyperfine levels. The maximum lifetime for the slow decay reached about 4 seconds, for a magnetic field of 0.5 Tesla. Furthermore, it was shown that the optical pumping efficiency was improved with respect to naturally doped Nd:YSiO crystals, which resulted in AFC echo efficiencies of up to 33%, similar to those obtained previously in europium and praseodymium crystals. An analysis showed that future experiments could reach efficiencies in the range 71-89%, by using cavity enhancement. Another interesting avenue of research would be to investigate coherent optical-spin manipulation, which would open up the possibility to perform on-demand quantum memory operation in Nd:YSiO. In this context we emphasize that a nuclear spin coherence times of 9 ms has already been observed in this material Wolfowicz et al. (2015).


We acknowledge useful discussions and some preliminary measurements by Félix Bussières and Thierry Armici. This work was supported by the European Research Council (ERC-AG MEC), the Swiss programme National Centres of Competence in Research (NCCR) project Quantum Science Technology (QSIT), and the Agence Nationale de la Recherche under grant agreements no. 145-CE26-0037-01 (DISCRYS).


  • Tittel et al. (2010) W. Tittel, M. Afzelius, T. Chanelière, R. Cone, S. Kröll, S. Moiseev,  and M. Sellars, Laser & Photonics Reviews 4, 244 (2010).
  • Afzelius et al. (2009) M. Afzelius, C. Simon, H. de Riedmatten,  and N. Gisin, Phys. Rev. A 79, 052329 (2009).
  • Hedges et al. (2010) M. P. Hedges, J. J. Longdell, Y. Li,  and M. J. Sellars, Nature 465, 1052 (2010).
  • Jobez et al. (2016) P. Jobez, N. Timoney, C. Laplane, J. Etesse, A. Ferrier, P. Goldner, N. Gisin,  and M. Afzelius, Phys. Rev. A 93, 032327 (2016).
  • Babbitt et al. (2007) W. R. Babbitt, M. A. Neifeld,  and K. D. Merkel, Journal of Luminescence 127, 152 (2007), proceedings of the Ninth International Meeting on Hole Burning, Single Molecule, and Related Spectroscopies: Science and Applications.
  • Linget et al. (2015) H. Linget, T. Chanelière, J.-L. Le Gouët, P. Berger, L. Morvan,  and A. Louchet-Chauvet, Phys. Rev. A 91, 023804 (2015).
  • Zhang et al. (2012) H. Zhang, M. Sabooni, L. Rippe, C. Kim, S. Kröll, L. V. Wang,  and P. R. Hemmer, Applied Physics Letters 100 (2012).
  • McAuslan et al. (2012) D. L. McAuslan, L. R. Taylor,  and J. J. Longdell, Applied Physics Letters 101 (2012).
  • Lauritzen et al. (2008) B. Lauritzen, S. R. Hastings-Simon, H. de Riedmatten, M. Afzelius,  and N. Gisin, Phys. Rev. A 78, 043402 (2008).
  • Afzelius et al. (2010a) M. Afzelius, M. U. Staudt, H. de Riedmatten, N. Gisin, O. Guillot-Noël, P. Goldner, R. Marino, P. Porcher, E. Cavalli,  and M. Bettinelli, Journal of Luminescence 130, 1566 (2010a).
  • Saglamyurek et al. (2015) E. Saglamyurek, T. Lutz, L. Veissier, M. P. Hedges, C. W. Thiel, R. L. Cone,  and W. Tittel, Phys. Rev. B 92, 241111 (2015).
  • Zhong et al. (2017) T. Zhong, J. M. Kindem, J. G. Bartholomew, J. Rochman, I. Craiciu, E. Miyazono, M. Bettinelli, E. Cavalli, V. Verma, S. W. Nam, F. Marsili, M. D. Shaw, A. D. Beyer,  and A. Faraon, Science 357, 1392 (2017).
  • Macfarlane and Shelby (1987) R. Macfarlane and R. Shelby, in Spectroscopy of Crystals Containing Rare Earth Ions, Modern Problems in Condensed Matter Sciences, Vol. 21, edited by A. Kaplyankii and R. Macfarlane (Elsevier Science Publishers, North Holland, Amsterdam, 1987) Chap. 4, pp. 51–184.
  • Macfarlane (2002) R. M. Macfarlane, J. Lumin. 100, 1 (2002).
  • Saglamyurek et al. (2011) E. Saglamyurek, N. Sinclair, J. Jin, J. A. Slater, D. Oblak, F. Bussières, M. George, R. Ricken, W. Sohler,  and W. Tittel, Nature 469, 512 (2011).
  • Clausen et al. (2011) C. Clausen, I. Usmani, F. Bussières, N. Sangouard, M. Afzelius, H. de Riedmatten,  and N. Gisin, Nature 469, 508 (2011).
  • Staudt et al. (2012) M. U. Staudt, I.-C. Hoi, P. Krantz, M. Sandberg, M. Simoen, P. Bushev, N. Sangouard, M. Afzelius, V. S. Shumeiko, G. Johansson, P. Delsing,  and C. M. Wilson, Journal of Physics B: Atomic, Molecular and Optical Physics 45, 124019 (2012).
  • Probst et al. (2013) S. Probst, H. Rotzinger, S. Wünsch, P. Jung, M. Jerger, M. Siegel, A. V. Ustinov,  and P. A. Bushev, Phys. Rev. Lett. 110, 157001 (2013).
  • Cruzeiro et al. (2017) E. Z. Cruzeiro, A. Tiranov, I. Usmani, C. Laplane, J. Lavoie, A. Ferrier, P. Goldner, N. Gisin,  and M. Afzelius, Phys. Rev. B 95, 205119 (2017).
  • Usmani et al. (2012) I. Usmani, C. Clausen, F. Bussières, N. Sangouard, M. Afzelius,  and N. Gisin, Nat Photon 6, 234 (2012).
  • Fröwis et al. (2017) F. Fröwis, P. C. Strassmann, A. Tiranov, C. Gut, J. Lavoie, N. Brunner, F. Bussières, M. Afzelius,  and N. Gisin, Nature Communications 8, 907 (2017).
  • Moiseev et al. (2010) S. A. Moiseev, S. N. Andrianov,  and F. F. Gubaidullin, Phys. Rev. A 82, 022311 (2010).
  • Afzelius and Simon (2010) M. Afzelius and C. Simon, Phys. Rev. A 82, 022310 (2010).
  • Sabooni et al. (2013) M. Sabooni, Q. Li, S. Kröll,  and L. Rippe, Phys. Rev. Lett. 110, 133604 (2013).
  • Jobez et al. (2014) P. Jobez, I. Usmani, N. Timoney, C. Laplane, N. Gisin,  and M. Afzelius, New Journal of Physics 16, 083005 (2014).
  • Walther et al. (2017) A. Walther, L. Rippe, L. V. Wang, S. Andersson-Engels,  and S. Kröll, Biomed. Opt. Express 8, 4523 (2017).
  • Hastings-Simon et al. (2008) S. R. Hastings-Simon, B. Lauritzen, M. U. Staudt, J. L. M. van Mechelen, C. Simon, H. de Riedmatten, M. Afzelius,  and N. Gisin, Phys. Rev. B 78, 085410 (2008).
  • Rančić et al. (2017) M. Rančić, M. P. Hedges, R. L. Ahlefeldt,  and M. J. Sellars, Nat Phys online publication, (2017).
  • Maksimov et al. (1970) B. A. Maksimov, V. V. Ilyukhin, U. A. Kharitonov,  and B. N. V, Kristallografiya 15, 926 (1970).
  • Maier-Flaig (2013) H. Maier-Flaig, Electron and nuclear spin properties of doped , PhD thesis, Karlsruhe Institute of Technology (2013).
  • Wolfowicz et al. (2015) G. Wolfowicz, H. Maier-Flaig, R. Marino, A. Ferrier, H. Vezin, J. J. L. Morton,  and P. Goldner, Phys. Rev. Lett. 114, 170503 (2015).
  • Li et al. (1992) C. Li, C. Wyon,  and R. Moncorge, IEEE Journal of Quantum Electronics 28, 1209 (1992).
  • Beach et al. (1990) R. Beach, M. Shinn, L. Davis, R. Solarz,  and W. Krupke, IEEE J. Quantum Electron 26, 1405 (1990).
  • Larson and Jeffries (1966a) G. H. Larson and C. D. Jeffries, Phys. Rev. 145, 311 (1966a).
  • Orbach (1961) R. Orbach, Proc. R. Soc. Lond. Ser. A 264, 458 (1961).
  • Scott and Jeffries (1962) P. L. Scott and C. D. Jeffries, Phys. Rev. 127, 32 (1962).
  • Larson and Jeffries (1966b) G. H. Larson and C. D. Jeffries, Phys. Rev. 141, 461 (1966b).
  • Abragam and Bleaney (1970) A. Abragam and B. Bleaney, Electronic Paramagnetic Resonance of Transition Ions (Clarendon Press, Oxford UK, 1970) chapter 1.12, page 78.
  • Zhou et al. (2015) Z.-Q. Zhou, Y.-L. Hua, X. Liu, G. Chen, J.-S. Xu, Y.-J. Han, C.-F. Li,  and G.-C. Guo, Phys. Rev. Lett. 115, 070502 (2015).
  • Afzelius et al. (2010b) M. Afzelius, I. Usmani, A. Amari, B. Lauritzen, A. Walther, C. Simon, N. Sangouard, J. Minář, H. de Riedmatten, N. Gisin,  and S. Kröll, Phys. Rev. Lett. 104, 040503 (2010b).
  • Gündoğan et al. (2013) M. Gündoğan, M. Mazzera, P. M. Ledingham, M. Cristiani,  and H. de Riedmatten, New Journal of Physics 15, 045012 (2013).
  • Bonarota et al. (2010) M. Bonarota, J. Ruggiero, J. L. L. Gouët,  and T. Chanelière, Phys. Rev. A 81, 033803 (2010).
  • Clausen et al. (2012) C. Clausen, F. Bussières, M. Afzelius,  and N. Gisin, Phys. Rev. Lett. 108, 190503 (2012).
  • Amari et al. (2010) A. Amari, A. Walther, M. Sabooni, M. Huang, S. Kroll, M. Afzelius, I. Usmani, B. Lauritzen, N. Sangouard, H. de Riedmatten,  and N. Gisin, Journal of Luminescence 130, 1579 (2010).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description