author: A. E. Garcı́a Pérez ; W. Aoki ; S. Inoue ; S. G. Ryan ; T.  K. Suzuki ; M. Chiba
bibliography: 12289.bib
date: ;
title: $^6$Li/$^7$Li estimates for metal-poor stars

The presence of the lithium-6 isotope in some metal-poor stars is a matter of surprise because of the high values observed. Non-standard models of Big Bang nucleosynthesis and pre-Galactic cosmic ray fusion and spallation have been proposed to explain these values. However, the observations of this light isotope are challenging which may make some detections disputable. The goal was to determine $^6\mathrm{Li}/^7\mathrm{Li}$ for a sample of metal-poor stars; three of them have been previously studied and the remaining two are new for this type of study. The purpose was to increase, if possible, the number of lithium-6 detections and to confirm previously published results. Spectra of the resonance doublet line of neutral lithium 670.78 nm were taken with the High Dispersion Spectrograph at the Subaru 8.2 m-telescope for a sample of five metal-poor stars (${-3.12\le\mathrm{[Fe/H]}\le-2.19}$). The contribution of lithium-6 to the total observed line profile was estimated from the 1D-LTE analysis of the line asymmetry. Observed asymmetries could be reproduced assuming isotopic abundance ratios $^6\mathrm{Li}/^7\mathrm{Li}$ of the order of: $0.004$ for , $\sim 0.010$ for and , $0.025$ for and $0.047$ for . We found that these results were very sensitive to several of the assumptions made in the analysis, in particular, the treatment of the residual structure in the analysed spectra. Our final estimates for the errors are respectively $\Delta$$^6\mathrm{Li}/^7\mathrm{Li}$ $= \pm 0.028, 0.029, 0.039, 0.025$ and $0.039$. The $^6\mathrm{Li}/^7\mathrm{Li}$ ratios for the sample are comparable to or even lower than these error values, so that detections of lithium-6 can not safely be claimed despite of the high resolving power ($R \sim 95\,000$) and $S/N$ (400-600).

# Introduction

The abundances of light elements in metal-poor stars are a good indicator of what happened early in the Universe and the Galaxy . The origins of the two stable isotopes of the light element lithium, $^6\mathrm{Li}$ and $^7\mathrm{Li}$, are still puzzling. Their derived abundances in metal-poor stars are difficult to reconcile simultaneously with standard predictions. Different production sites in the early Galaxy have been discussed: Big Bang nucleosynthesis , cosmic ray spallation and/or $\alpha$-fusion and neutrino spallation (see e.g. ). The most widely discussed picture is that lithium-6 has a cosmic-ray-induced origin (principally $\alpha$-fusion) and lithium-7 a predominantly primordial one .

The observed abundances of $^6\mathrm{Li}$ are used to discriminate between different production sites: pre-Galactic (e.g. ), the Galactic interstellar medium (e.g. ), superbubbles (e.g. ), fast supernova eject (e.g. ), in-situ stellar flares , etc. The problem with this isotope is that detections for metal-poor stars are very challenging and very rare. reported the first probable detection of $^6\mathrm{Li}$ in HD 84937 ($^6\mathrm{Li}/^7\mathrm{Li} = 0.05\pm0.02$) which was later confirmed by other works . A few years after this first detection, a detection for another star, , ), was reported. Possible detections for two other stars were provided by and : BD$+42^\circ\,2667$ ($^6\mathrm{Li}/^7\mathrm{Li} = 0.10$) and G 271-162 ($^6\mathrm{Li}/^7\mathrm{Li} = 0.02\pm0.01$) respectively. There are results published for other metal-poor stars also, however the values are mostly upper limits. Our group determined $^6\mathrm{Li}/^7\mathrm{Li}< 0.018$ for the metal-poor star HD140283 ; this star may have partly depleted its initial content . The work of has the highest number of published detections in metal-poor stars: nine. The detections for these stars appear to lie in a plateau that, if confirmed, would be difficult to explain within the standard framework of light element production. One would expect an evolution of the lithium-6 abundances with metallicity because of the contribution of Galactic cosmic ray spallation. Analyses of $^6$Li for very metal-poor stars ($\mathrm{[Fe/H]}<-2.5$) are still very limited, so new samples of a few stars have a large impact.

A possible solution to the apparent existence of a plateau and also to the high abundance of that plateau could lie in a pre-Galactic origin of the cosmic rays (see, however ). Non-standard models of Big Bang nucleosynthesis (BBN), generally involving particles of uncertain character, are also considered as a possible solution. The importance of establishing whether or not $^6$Li has been reliably detected in metal-poor stars is thus readily apparent: the previously published detections have been a significant motivator of investigations of the role of supersymmetry in modifying nucleosynthesis during the Big Bang. The detection (or otherwise) of $^6$Li could therefore potentially provide crucial, rare observational constraints on the nature of certain hypothesised supersymmetric particles and dark matter candidates, some of which cannot be studied experimentally with current or near-term experimental facilities (e.g. ).

Although detections for $^7\mathrm{Li}$ are more common than for $^6\mathrm{Li}$, the origin of the heavier isotope is also quite intriguing. observed essentially the same lithium content, the so called Spite value, in a sample of metal-poor dwarf stars with effective temperatures $5800\le T_\mathrm{eff}\le 6250$, despite their metallicities ranging from $-2.2\le\mathrm{[Fe/H]}\le-1.1$. A plausible interpretation for the plateau is that it is close to the primordial value. $^7\mathrm{Li}$ is the only light isotope besides H, D, $^4$He and $^3$He that, according to standard BBN models, is predicted to be produced in significant amounts in the Big Bang. Predictions for the primordial lithium content based on standard Big Bang nucleosynthesis (SBBN) in combination with the baryon density inferred from the *Wilkinson Microwave Anisotropy Probe Mission* (WMAP) give a value two or three times higher than the observed plateau value, $A\mathrm{(Li)}\footnote{$A(\mathrm{Li}) = \log{N_{\mathrm{Li}}/N_\mathrm{H}}$+12} = 2.65$ vs $2.1$-$2.3$ dex .

Lithium is destroyed by nuclear reactions at high temperatures. In the presence of mixing, material from the stellar atmospheres can reach depths of high temperatures where lithium can burn or can mix with material that has burnt it. According to some stellar models, little or no lithium-7 destruction is expected for stars in the Spite plateau (e.g. and references therein) so that these stars should retain their initial content. However atomic diffusion is expected to act in warm main-sequence stars and alter the initial photospheric chemical composition (e.g. ). Stellar models including diffusion and/or stellar rotation have been suggested to reconcile the observations of metal-poor stars with the SBBN predictions, e.g and . Observers have searched for signatures of mixing by measuring the lithium-7 abundance dispersion and by studying the behaviour of the abundances with stellar temperatures and metallicity. There is not consensus yet whether the Spite plateau is a plateau or not and whether lithium abundances are dispersed or not (e.g. ). Lithium-6 is however more fragile than lithium-7 (burning at temperatures $T \sim 1.6\times10^6$ K versus $T \sim 2.0\times10^6$ K), so an observable content of the lighter isotope sets constraints on the degree of depletion of the other one. Stellar models without diffusion or rotational mixing predict part destruction of lithium-6 in the pre-main-sequence phase for most of the stars with published lithium-6 detections . According to neither the rotational mixing models of nor the diffusion models (including turbulent diffusion) from seem to give a solution to both lithium problems simultaneously. A high depletion of lithium-6 would be difficult to reconcile with observations of the other stable isotope because a very high initial content of lithium-6 would generally imply a significant lithium-7 production by cosmic ray nucleosynthesis. Nevertheless, recent results for extremely metal-poor stars show a correlation between lithium-7 abundance and metallicity which may have an interpretation in terms of Galactic evolution by cosmic ray spallation.

Given the importance of the published lithium-6 detections for the understanding of Big Bang and cosmic ray nucleosynthesis, it is paramount to confirm the detections and, if possible, to add more. This is particularly important because of the very weak signal presented by $^6\mathrm{Li}$: a small depression of the red wing of the $^7\mathrm{Li}$ doublet, at a level comparable to possible sources of systematic error. We present here $^6\mathrm{Li}/^7\mathrm{Li}$ estimates based on observations taken with the Subaru 8.2m-telescope for five metal-poor stars, of which three are in common with previous studies. The observations are described in § , while the stellar parameter determinations and the spectral synthesis are in § . In §  we present the stellar and instrumental broadening calculations. In §  we explain the procedure followed to get our best isotopic abundance ratio estimates, and we present the results and their sensitivity to different uncertainties, while in §  we discuss the results in terms of the observational quality and the sensitivity to the spectral analysis. We finish with conclusions in § .

#  Observations

The $^6\mathrm{Li}/^7\mathrm{Li}$ estimates presented here are based on the spectral analysis of the observed resonance doublet line at 670.78 nm. Spectra were taken for eight metal-poor stars with the High Dispersion Spectrograph (HDS, ) mounted on the Subaru 8.2 m-telescope on the nights of May 17-19, 2005 (UT). Results for five stars are presented here; the observations of three other stars were not completed in this run to achieve the data quality required for this study. Values for exposure times are listed in Table , together with the $V$ magnitude, date of the observations, exposure time, the central wavelength of the instrumental setting employed and the measured heliocentric radial velocity of the star. Total exposures were split into several shorter ones for cosmic ray rejection. Three of the stars were observed in two consecutive nights. Observations were taken using a $2048\times4100$ EEV CCD ( $13.5\,\mu$m pixel size), a slit width of 0.4 arcsec (0.20 mm) and two standard instrumental settings: one with central wavelength at 608.11 nm (StdYd) and the other at 620.78 nm (StdYb) (0.002 nm/pixel dispersion). Special care was taken for the calibration; a high number of flat-field frames ($\sim20$) were acquired and ThAr lamp spectra were taken at several times. Fast rotator stars were also observed, in case their spectra were required for flat-fielding correction.

The observational data were reduced under using the echelle package and following the standard procedure: bias subtraction, flat-fielding correction, spectra extraction, wavelength calibration, normalisation and combination of individual spectra (spectral orders were not merged). For the wavelength calibration, an average of 20-40 Th lines per spectral order were identified. The wavelengths of these lines were used to convert pixel values into wavelengths. The wavelength solutions are based on a fifth order polynomial fit to the transformation. Typical values for the rms of these fits are 0.06 pm. If necessary, the fits were interpolated in time before being applied to the observed spectra. The normalisation consisted of dividing the spectra by their cubic spline fit with a wavelength scale of 0.3 nm. Spectra were co-added and weighted according to their exposure times.

The co-added spectra are of very high quality. The signal-to-noise ratio ($S/N$) per pixel was estimated from the measurements of noise on both sides of 670.78 nm line in a window of 0.07 nm, values are in the range 417-611. The resolving power ($R = \lambda/\Delta\lambda_\mathrm{FWHM}$) was estimated from the full width at half maximum (FWHM) of a Gaussian fit to a set of ThAr lines lying in the spectral order that contains the lithium resonance doublet line; results are in the range $R \sim$`<!-- -->`{=html}90 000-100 000.

#  Stellar parameters, model atmospheres and spectral synthesis

The observed spectra were analysed using spectral synthesis. The first step in the analysis was to get an estimate of the fundamental parameters describing the stellar atmospheres: effective temperature $T_\mathrm{eff}$, the surface gravity $\log g$, the metallicity $\mathrm{[Fe}/\mathrm{H}]$ and the microturbulence $\xi_\mathrm{t}$. The values listed in Table  were taken from , with the exception of the metallicity. The choice of the stellar parameters is not very important for the accuracy of the $^6\mathrm{Li}/^7\mathrm{Li}$ final results (e.g. ), however we have tried to have the maximum consistency. Although for some stars there are more recent values than the values assumed here, we have chosen those which were determined using the same indicators and in a similar way: $(b-y)$ and $(V-K)$ colors for $T_\mathrm{eff}$, and Hipparcos parallaxes for $\log{g}$. The uncertainty in the measured parallax of is very high so an estimate coming from the calibration of the Strömgren indices based on Hipparcos parallaxes was used instead. Our metallicity estimates come from the measured equivalent widths of a set of lines: $\lambda\,614.92$, $\lambda\,624.76$, $\lambda\,641.69$, $\lambda\,643.27$, and $\lambda\,645.64$ nm. Iron abundances were derived and taken as a metallicity estimate. The solar values were determined from the measured equivalent widths of these lines as in \[but now using new 3D solar C and O abundance values\]. For , the lines above were not detected, but one at 523.4 nm was, giving $\mathrm{[Fe/H]}=-2.97$. As our metallicity calculation for this star would be based on only one line, we preferred to take the value from that is determined from several lines: $\mathrm{[Fe/H]} = -3.12$.

Equipped with the stellar parameters (see Table ), model atmospheres were computed using the code version in . The code solves the equations of hydrostatic equilibrium and energy equilibrium assuming that the atmospheres have either plane-parallel or spherical symmetry (1D); we assumed plane-parallel geometry in our computations. The radiation and the matter are assumed to be coupled, so the radiation can be described by the local atmospheric conditions; radiation is said to be in local thermodynamical equilibrium (LTE). The same assumptions were made for the spectral synthesis carried out with the code from the Uppsala package.

The validity of these assumptions is a subject of discussion. Atmospheric convective motions can produce asymmetries in spectral line profiles. These asymmetries cannot be reproduced with the use of 1D model atmospheres, so in principle our results for isotopic abundance ratios based on observations of the profile asymmetry may be in error. A more appropriate treatment of the convection, like the one employed in the construction of hydrodynamic model atmospheres (3D), may be required for a better estimate of the line asymmetry due to the isotope lithium-6 . However Asplund et al. investigated the effects of the use of 3D model atmospheres on their $^6\mathrm{Li}/^7\mathrm{Li}$ calculations and found no significant differences with respect to the use of 1D models. The investigation was done assuming LTE conditions. Studies of the formation of the 670.7 line (e.g. ) suggest that lithium abundances derived from this line observed in spectra of metal-poor stars may be affected by departures from LTE by up to 0.03 dex. showed that for the metal-poor star HD 74000, 1D-LTE and 3D-NLTE results for $^6\mathrm{Li}/^7\mathrm{Li}$ disagree (0.024 vs 0.006). The 3D-NLTE analysis of Cayrel et al. suggests there is no $^6\mathrm{Li}$ in the star. One of the purposes of our work is to compare our results with values reported in other works which are basically 1D-LTE results; that is why we have made a 1D-LTE assumption. The 3D-NLTE effects cannot be neglected when discussing implications of our results for the modelling of chemical evolution of the universe.

The 670.7 line is a multiplet and has two isotopic components, $^6\mathrm{Li}$ and $^7\mathrm{Li}$. The assumed values of wavelengths and $\log{gf}$ of these transitions are the same as in . Laboratory measurements suggest they are known with a high accuracy. Atomic data for the lines used for metallicity determinations were taken from and from for lines other than these.

```latex
\begin{table*}
%\begin{minipage}[t]{\columnwidth}
\caption{Observational information}
\label{stepar}
\centering
\renewcommand{\footnoterule}{}  % to avoid a line before footnotes
\begin{tabular}{lccccccccc}
\hline \hline
Star& $V$&Date & Exp. Time  & Inst. Setting $\lambda_0$ & $v_\mathrm{rad\ H}$&$T_{\mathrm{eff}}$ &$\log{g}$ & [Fe/H] & $\xi_t$ \\
 &(mag)&(UT)& (s) & (nm) &(km\,s$^{-1})$&(K)&(cgs)&(dex)&(km\,s$^{-1}$)\\
\hline

BD$-04^\circ\,3208$   & 9.99 & 18-05-05 & $4\times1800$ &620.78 (StdYb)& 58.7 & 6338   & 4.00   & $-2.28$  & 1.5 \\
 & &19-05-05 &$2\times1800$\\
G 64-37                       &11.14&  18-05-05 &$7\times1800$ &620.78 (StdYb) &82.1& 6318   & 4.16   & $-3.12$ & 1.5 \\
&&19-05-05 &$7\times1800$\\
BD$+02^\circ\,3375$  & 9.93&19-05-05 & $4\times1800$ &620.78 (StdYb) &$-397.1$&     5855   & 4.16   &$-2.19$  & 1.5 \\
BD$+20^\circ\,3603$  & 9.69 & 18-05-05 &  $1\times1800$ &620.78 (StdYb) &$-241.6$&   6092   & 4.04   & $-2.22$  & 1.5 \\
&&&$2\times2400$\\
&&19-05-05&$1\times1800$\\
BD$+26^\circ\,3578$  &9.37& 17-05-05  & $3\times1800$  & 608.11 (StdYd) &$-128.2$ &  6239   & 3.87   & $-2.33$  & 1.5 \\
&&&$1\times2400$\\
%081543  &     6514   & 4.23   & -2.72  & 1.5 & 17,18-05-05 & 620.78\\
\hline
\end{tabular}
%\end{minipage}
\end{table*}
```

#  Stellar and instrumental broadening

Our estimates for the isotopic abundance ratios ($^6\mathrm{Li}/^7\mathrm{Li}$) are based on the comparison of the observed and synthesised profile of the resonance 670.78 nm doublet line. Before the comparison, the synthetic spectrum has to be convolved by a function that reproduces the instrumental and the stellar line broadening; for this work we have assumed a Gaussian function. We adopted the approach of including the stellar and instrumental broadening contribution in a single free parameter, the broadening parameter $\zeta_\mathrm{broad}$ ($= \Delta\lambda_\mathrm{FWHM}$ \[km s$^{-1}$\]), that was determined from the analysis of up to 26 observed spectral lines: (588.9 nm), (470.2, 552.8 nm), (558.8, 612.2, 616.2, 643.9 nm), (444.3, 446.8, 450.1, 453.4, 456.3, 457.1 nm), (441.5, 452.8, 491.8, 492.0, 523.3, 526.9, 532.8, 537.1, 539.7, 561.5 nm) and (458.3, 492.3, 501.8 nm).

Different $R$ values for the different lines could compromise the accuracy of our determinations, however $R$ measurements from ThAr lamp lines in the range $\lambda = 550$ to 700 nm did not show significant systematic variations with wavelength ($\Delta R \sim 1\%$). The spread of values at a given wavelength was higher, though still not so significant ($\sim \pm 5\%$).

The other two free parameters of the analysis were the rest wavelength of the observed spectra and the chemical abundance ($\log\epsilon_X$) of the line studied [^1]. The free parameters were varied until the best fit to the observations was achieved, using the $\chi^2$ technique.

Values of the broadening parameter versus the measured equivalent widths of the spectral lines are shown in Fig.  for all the observed stars of our sample. Not all lines were measurable in all stars. The data in Fig.  suggest that lines stronger than $\sim 7$ pm require lower broadening values so mean values from the weaker lines were determined for $\zeta_\mathrm{broad}$ and taken as our best estimate, see Table . The weakest lines observed have equivalent widths in the range 0.4 to 3.8 pm depending on the star. For comparison, the observed Li lines have equivalent widths in the range 1.6 to 3.1 pm (see Table ), values are shown in Fig.  as dashed vertical lines; only has a lithium line that is much weaker than the minimum equivalent width included in its broadening calculations (2.6 vs 3.8 pm). We assumed that for this star the broadening parameter does not change significantly with equivalent width in the weak-line regime.

The mean values of the broadening parameter, listed in Table , were well determined for all the stars with typical standard deviation values of 0.15-0.40 km s$^{-1}$. The observed 612.2 line and its best fit is shown in Figure  where you can appreciate the low sensitivity of the fit to changes in the broadening parameter of the order of $\pm0.5\,\mathrm{km}\,\mathrm{s}^{-1}$, this value is higher than the maximum of the $\sigma_\zeta$ values listed in column 3 in Table .

#  Isotopic abundance ratios

## Calculations

The $^6\mathrm{Li}$ isotope is shifted by $\sim 0.015$ nm relative to $^7\mathrm{Li}$ in the neutral lithium resonance doublet line. This produces an asymmetry in the already asymmetric profile of the red wing of the line. Under the assumptions of 1D model atmospheres, observed asymmetries may only be due to the multiplet configuration of the lithium resonance line and to the isotopic shifts (for 3D effects, see discussion in § ).

The best fit to the observations of the lithium line was determined from a $\chi^2$ analysis ($\chi^2 = \sum{(O_{i}-S_{i})^2/\sigma^{2}_\mathrm{S/N}}$; where $O_{i} =$ normalised observation in pixel $i$, $S_{i} =$ normalised synthesis and $\sigma_\mathrm{S/N} =$ noise based on the $S/N$) with three free parameters. A first, approximate, continuum normalisation was done in . The continuum level of the spectra was redetermined and fixed prior to the $\chi^2$ analysis, not taken as a free parameter as in . The continuum levels were determined by linear interpolation between two windows (see Figure ), one on each side of the lithium line. Different regions were selected for different stars in order to minimise the presence of structure. The presence of ripples that remain after this continuum normalisation compromises the accuracy of the continuum level as well as the determination of the $S/N$. We return to this issue below.

|                      |                        |                |     |     |     |
|:---------------------|:----------------------:|:--------------:|:---:|:---:|:---:|
| Star                 | $\zeta_\mathrm{broad}$ | $\sigma_\zeta$ |     |     |     |
|                      |     (km s$^{-1})$      | (km s$^{-1}$)  |     |     |     |
| BD $-04^\circ\,3208$ |          6.61          |      0.18      |     |     |     |
| G 64-37              |          6.24          |      0.39      |     |     |     |
| BD $+02^\circ\,3375$ |          5.20          |      0.32      |     |     |     |
| BD $+20^\circ\,3603$ |          5.84          |      0.23      |     |     |     |
| BD $+26^\circ\,3578$ |          6.49          |      0.16      |     |     |     |

Broadening parameter values and their standard deviation.

```latex
\begin{table*}
%\begin{minipage}[t]{\columnwidth}
\caption{Equivalent widths and $S/N$ measurements of the \ion{Li}{i} 670.78\,nm together with the best results for the $\chi^2$ three free parameters and the sensitivity of the analysis.}
\label{lithium}
\centering
\renewcommand{\footnoterule}{}  % to avoid a line before footnotes
\begin{tabular}{lccccccccccc}
\hline \hline
Star & $W_\lambda$& $S/N$ &$ \Delta\lambda$ & $A(^7$Li$)$ &\multicolumn{3}{c}{$^6\mathrm{Li}/^7\mathrm{Li}$}&&\multicolumn{3}{c} {$\Delta^6\mathrm{Li}/^7\mathrm{Li}$} \\
\cline{6-8}
\cline{10-12}

%\cline{8-11}
%\rule{0pt}{2.6ex}
&(pm)& &(nm)&(dex)\\

&&&&&ncor &cor & best-esti && cont   &  $\Delta\lambda_\mathrm{cover}$& $\zeta_\mathrm{broad}-\sigma_\zeta$\\
\hline
BD $-04^\circ\,3208$ &  2.37& 602 & 0.01036& 2.27 & 0.047& $0.043$ &$0.047\pm0.039$&&  $-0.037$ & 0.003& 0.009\\
G 64-37 & 1.60 & 496 & $-0.00125$& 2.04  & 0.006 &$0.000$ &$0.006\pm0.039$&& 0.028&0.023&0.013\\
BD $+02^\circ\,3375$ &  3.10 & 417& 0.00745 & 2.04 & 0.007 & $-0.001$ &$0.007\pm0.029$&&0.023 &0.007& 0.013 \\
%BD$+20^\circ\,3603$ &  2.56&492 & 0.00713 & 2.14 & $4.026\pm3.29$& 0.0079 & 0.0261 & $-0.0028$ & 0.0114&$-0.0011$&0.0019\\
BD $+20^\circ\,3603$ &  2.56&437& 0.00626& 2.14 & 0.004&0.025 & $0.025\pm0.025$ & &  $-0.003$ &0.002& 0.011\\

BD $+26^\circ\,3578$ &  2.37& 611&  0.00163 & 2.21 &0.074 &$0.004$ &$0.004\pm0.028$&&  0.017 &  0.014& 0.016\\

\hline
\end{tabular}
%\end{minipage}
\end{table*}
```

With the broadening parameter and continuum level fixed, the free parameters of the $\chi^2$ analysis were the lithium-7 abundance $A(^7\mathrm{Li})$, the isotopic lithium abundance ratio $^6\mathrm{Li}/^7\mathrm{Li}$ and the wavelength zero point of the observations.The three free parameters were allowed to vary until $\chi^2$ was minimised. They were changed in steps of: 0.02 dex in lithium abundance, 0.01 in isotopic abundance ratios and 0.0002 nm in wavelength shift, with typically a $11\times13\times11$ points grid. Wavelength points in a range from 670.743 to 670.823 nm were included in the $\chi^2$ calculations (36-37 wavelength points). Our best $^6\mathrm{Li}/^7\mathrm{Li}$ estimates for all the stars are listed in Table  (column 8) together with the $S/N$ around the resonance neutral lithium line (column 3), the wavelength shift (column 4) and the lithium-7 abundance (column 5). Figures ,  ,  ,   and   show the observed metal-poor spectra and their best synthetic fits for $^6\mathrm{Li}/^7\mathrm{Li}$ = 0.00, 0.05, 0.10 together with their residuals (observations $-$ synthesis in $\%$). The goodness of the final fits can be appreciated from the $\Delta\chi^2$ values ($\Delta\chi^2 = \chi{^2}-$ the minimum value) presented in these figures. Note that the $\chi{^2}$ values were not divided by the degree of freedom so the values are non-reduced $\chi{^2}$ values.

##  Contamination

Possible "contamination" of the lithium line was investigated by comparing the spectra of the observed metal-poor stars with the spectra of hot fast rotator stars taken the same night. This comparison is shown in Figure  where two plots are shown for stars observed in two different nights, each corresponding to a different date of observation. The associated hot fast rotator spectra of the four metal-poor stars , , and show tiny features of peak-to-valley range $\sim 0.003$-$0.006$ (in normalised flux units) that suggest that their observed lithium line suffer from some contamination. At most one or two of such features could be identified with diffuse interstellar bands. We presume the majority is residual fringing or residual flat-field errors. Note that, although the observed ripples would not significantly affect common abundance determinations, they can have a great impact on the isotopic abundance ratios. One approach to quantify the effect of these ripples in the $^6\mathrm{Li}/^7\mathrm{Li}$ results is to divide the metal-poor spectra by the spectra of the hot fast rotator stars. The division has two effects: one is to alter the flux in the Li line directly, and the other is to require a revision of the continuum fit, which again, affects the normalised flux in the line. The correction is not perfect; there are features in the observed metal-poor spectra that do not match the features in the fast-rotator spectra (e.g. ) so we cannot be sure the lithium line is free of ripples at the level of accuracy that we require for our determinations even after this correction. Results for the isotopic abundance ratios coming from non-corrected and the corrected (divided) spectra are presented in Table , columns 6 and 7 respectively.

The observations of were from only one night and, unfortunately, the $S/N$ of the fast rotator spectra for that night was not very high, $S/N \sim 250$. To avoid degrading the quality of the analysed spectra, the division by a fast rotator spectrum should be done only when necessary and so for this star it was not done. Had the division been done for this star, the signal-to-noise of the analysed observations would have decreased ($S/N \sim437$ to 300) but the $^6\mathrm{Li}/^7\mathrm{Li}$ results would have not changed significantly, 0.007 vs $-0.001$ (see values in Table ). For this star we have adopted the $^6\mathrm{Li}/^7\mathrm{Li}$ value based on the non-corrected spectra (higher $S/N$) as our best estimate, while the ripple effect ($^6\mathrm{Li}/^7\mathrm{Li}$$|_{cor}-$$^6\mathrm{Li}/^7\mathrm{Li}$$|_{ncor}$) is included in the error bar.

The problem of the low $S/N$ of the fast rotator spectrum was not such a big issue for the other stars observed that night, , and , as they were observed also on a second night. Therefore, the $S/N$ per night required to achieve the high target $S/N$ was lower than for the stars which were observed on only one night. Still, the low $S/N$ makes it difficult to judge by eye whether there is, or not, contamination.

Observations of suggest no contamination for the first night, whereas it is hard to say for the second night. The $^6\mathrm{Li}/^7\mathrm{Li}$ corrections for are very small, $\Delta$$^6\mathrm{Li}/^7\mathrm{Li}$ $\sim 0.007$, so for this star as well, no correction was done. The comparison of the spectra of and its associated fast rotator does not show any sign of strong contamination. Hence, our isotopic ratio results coming from non-corrected and corrected spectra are quite similar, 0.047 vs 0.043, our best estimate is based on the non-corrected spectra.

Results for were more affected by the division of spectra: $^6\mathrm{Li}/^7\mathrm{Li}$ $\sim 0.004$ (non-corrected) vs $0.025$ (corrected). We consider it more appropriate for our best estimate to adopt a value based on the corrected spectra, given that there are indications of contamination.

This was also the case for ; ripples are observed in the hot spectrum at the wavelengths of the 670.8 nm line, therefore we applied the correction to get our $^6\mathrm{Li}/^7\mathrm{Li}$ best estimate, $0.004$ (corrected) vs $0.074$ (non-corrected). The match between the ripples in the metal-poor spectrum and in the associated fast rotator is quite impressive (see Figure ). The observed residual fringing in the fast rotator spectra recorded on May $17^{\mathrm{th}}$, 2005 was greater than on the other nights when a different instrumental setting was used. was the only programme star observed on that night (see Table ). The observed spectra shown in Figures  and are the divided spectra.

##  Sensitivity

The sensitivities of our results to the different assumptions made in the analysis have been studied. Concerning the normalisation, changes of $0.1\%$ in the continuum level, where the typical line depth is $10\%$, has an effect of 0.02 on the $^6\mathrm{Li}/^7\mathrm{Li}$ results, so the normalisation of the spectra is an important issue.

After the initial continuum fit (see $\S$ ), the spectrum of showed a remaining tilt, so it was re-normalised in using a linear function similar to that employed in the $\chi^2$ analysis. After this was done, the match between the ripples in the spectrum of and in the associated fast rotator spectrum was excellent, hence we have used a linear fit in our analysis. For stars with only minor differences in the measured continuum level on either side of the lithium line, the use of the mean value of the continuum fluxes, rather than a linear interpolation, does not have a big effect on the $^6\mathrm{Li}/^7\mathrm{Li}$ results. However, for and , the isotope ratio changed by $\Delta$$^6\mathrm{Li}/^7\mathrm{Li}$ $= 0.005$ when using the mean-value approach.

The level measurements can change with the spectral window choice, different areas were selected for different stars and they are shown by horizontal bars in Fig. . A change of choice within the range of possibilities would affect the results significantly for all the stars except . Even a shift of the window in the left wing by $<0.1$ nm will change the isotope abundance ratios significantly, the changes ranging from $\Delta$$^6\mathrm{Li}/^7\mathrm{Li}$ $= -0.037$ () to $0.028$ (), see column 9 in Table . and are the stars with the highest uncertainty in the continuum location among our sample. For these two stars, the use of the fast rotator spectrum did not help for discriminating between two different windows in the left wing, either because the small structures observed in the fast rotator spectra did not match the ones in the metal-poor spectra or because the spectra of the comparison star had a low $S/N$. The result of this was that two normalisation levels were possible. Summarising, isotopic results for these stars have a high uncertainty because of the uncertain continuum location.

Another aspect to take into account is the number of wavelength points included in the $\chi^2$ calculation, in particular in the red wings of the lithium line. A decrease of typically four or five wavelength points ($\sim 0.010$ nm) in the red wing affects the isotope ratio results by a few tenths of one percent for three stars (Table , column 10). The stars and were more affected by these changes, $\Delta$$^6\mathrm{Li}/^7\mathrm{Li}$ = 0.023 and 0.014 respectively.

Results also show sensitivities to the assumed values of macroturbulence; we determined new isotope ratio results for broadening parameters values decreased by their $\sigma_\zeta$. We find a mean sensitivity of $\Delta$$^6\mathrm{Li}/^7\mathrm{Li}$ $\sim 0.01$ (see values in column 11 of Table ), although the sensitivity for is even higher, 0.016.

We have investigated the effects of changing the stellar parameters for one particular case, , according to their uncertainty: $\Delta T_\mathrm{eff} = 100\,$K, $\Delta\log{g} = 0.2$, $\Delta[\mathrm{Fe}/\mathrm{H}] = 0.2\,$dex and $\Delta \xi_\mathrm{t} = 0.2\,$km s$^{-1}$. Isotopic ratio determinations are not very sensitive to these changes directly, $\Delta^6\mathrm{Li}/^7\mathrm{Li}<$ 0.004 (see top row of Table ). A second effect of changing the stellar parameters is to alter the inferred broadening parameter (see second row of Table ). For example, an increase in the microturbulence results in lower $\zeta_\mathrm{broad}$ values; the broadening parameter compensates particularly for those parameters most strongly connected with the profile width: $\Delta \xi_\mathrm{t}$ and $\log{g}$. We have shown already (Table , column 11) that the isotopic ratio is quite sensitive to the broadening parameter. The total effect of increasing the microturbulence will be a $^6\mathrm{Li}/^7\mathrm{Li}$ increase of $\sim 0.003$ (fourth row of Table ), $-0.003$ from microturbulence directly, and $+0.006$ from the resulting change in $\zeta_\mathrm{broad}$ (third row of Table ). A total error associated with the uncertainty in the stellar parameters including the broadening parameter is estimated from the quadrature sum of the values in fourth row in Table  and is $\Delta$$^6\mathrm{Li}/^7\mathrm{Li}$ $\sim 0.006$. Note that this excludes the direct effect of the uncertainty in the broadening parameter, which is given in Table , column 11. Summarising, the two first rows of Table  list the partial sensitivities to the stellar parameters, the third row the partial sensitivity of $^6\mathrm{Li}/^7\mathrm{Li}$ results associated to the change in the broadening parameter as a result of changing the stellar parameters, and the last row the total combined effect (stellar parameters + broadening parameter). The quadrature sum of all the uncertainties are listed in the last column.

A summary of all the individual sensitivities are given as $\Delta$$^6\mathrm{Li}/^7\mathrm{Li}$ values in Table  (columns 9 to 11) and the final row of Table . The quadrature sum of these plus the $^6\mathrm{Li}/^7\mathrm{Li}$$|_{cor}-$$^6\mathrm{Li}/^7\mathrm{Li}$$|_{ncor}$ values were taken as our best estimates of the total $^6\mathrm{Li}/^7\mathrm{Li}$ errors (Table , column 8). We considered the division by a hot fast rotator spectrum to be essential, not merely optional, for , so the $^6\mathrm{Li}/^7\mathrm{Li}$ error bar for this star does not include the uncertainty associated with that correction, otherwise the error bar would be overestimated.

Another concern is whether asymmetries are intrinsic to even single lines. Following , we have studied the asymmetry of the observed profiles of two lines at 639.36 nm and 649.50 nm in the spectrum of . For that, the profiles were re-scaled to the same line depth, they were averaged and after that they were convolved with a Gaussian function of width 2.39 km s$^{-1}$ to correct for the different thermal broadening of the Fe and Li lines (see ). The magnitude of the asymmetries shown in the determined profile depends on the assumed normalisation; the presence of fringing can introduce asymmetries. Continuum levels that are acceptable choices can produce an asymmetry whose additional depth is of the same order of magnitude as that produced by a displaced line of equivalent width $4\%$ of the main one. Therefore, despite the quality of our observations, it is difficult to determine whether the asymmetry is real or just a product of the normalisation.

```latex
\begin{table*}
%\begin{minipage}[t]{\columnwidth}
\caption{Sensitivity of the {$^6\mathrm{Li}/^7\mathrm{Li}$} and the broadening parameter results to the stellar parameters uncertainties for \object{BD $-04^\circ\,3208$.}}
\label{senstepar}
\centering
\renewcommand{\footnoterule}{}  % to avoid a line before footnotes
\begin{tabular}{lcccccc}
\hline \hline
Star      &  $\Delta T_{\mathrm{eff}}$ &$\Delta\log{g}$ &$\Delta$ [Fe/H] & $\Delta \xi_t$ & $\Delta \zeta_\mathrm{broad}$ & $\sigma _ {^6\mathrm{Li}/^7\mathrm{Li}}$\\
&(K)&(cgs)&(dex)&(\,km\,s$^{-1}$)&(\,km\,s$^{-1}$)\\
& $+100$& $+0.2$& $+0.2$& $+0.2$ & $+0.2$\\
\hline
$\Delta^6\mathrm{Li}/^7\mathrm{Li}_\mathrm{direct}$ & $-0.004$ & $-0.004$ & $0.001$ & $-0.003$ & $-0.008$\\
$\Delta \zeta_\mathrm{broad}$ \footnote{km\,s$^{-1}$} & 0.03 & $-0.08$ & $-0.02$ & $-0.15$& &\\
$\Delta^6\mathrm{Li}/^7\mathrm{Li}_\mathrm{broad}$ & $-0.001$ & $0.003$ & $0.001$ & $0.006$ & & \\
$\Delta^6\mathrm{Li}/^7\mathrm{Li}_\mathrm{total}$ & $-0.005$ & $-0.001$ & $0.002$ & $0.003$ & & 0.006\\
\hline
$^3$ (\,km\,s$^{-1}$)\\
\end{tabular}
\end{table*}
```

## Summary

Our best estimates of the $^6\mathrm{Li}/^7\mathrm{Li}$ ratios in these stars are set out in Table , column 8 and displayed in Fig. . The tabulated uncertainty is a quadrature sum of the contribution from uncertainties in: whether or not division by a hot fast rotator is performed (§  and Table , columns 6 and 7), $\sim 0.010$, changes resulting from different choices of continuum windows (§  and Table , column 9), $\sim 0.022$; changes resulting from different choices of the wavelength range over which the observed and synthetic spectra are compared via the $\chi^2$ statistics (Table , column 10), $\sim 0.010$; the star's broadening parameter (Table , column 3 and Table , column 11), $\sim 0.012$; and the star's atmospheric parameters (Table , row 1), $\sim 0.006$. We note, but do not include, further possible error sources: whether a mean value of the continuum windows or a linear interpolation between continuum windows is adopted as "the continuum"; and whether even single spectral lines are asymmetric. The tabulated uncertainties range from $\pm 0.025$ to $\pm 0.039$, and are dominated by the uncertainty in the continuum, even though the $S/N$ exceeds 400 in every case. The reason that the continuum placement continues to exert such an influence on the outcome of the analysis lies substantially in the problem of small structures seen in the reduced data, and that although they are small by most subjective measures, they are significant relative to the very weak signal of $^6$Li, the spectral line of which has a central depth of probably only $\leq 0.003$ normalised flux units (i.e. $\leq 4\%$ of the depth of the $^7$Li, which in turn is $\sim 0.08$ flux units). Based on the extreme sensitivity of the isotope ratio results to alternative choices in the spectral analysis, we are certain that we could not infer $^6\mathrm{Li}/^7\mathrm{Li}$ to better than $\pm0.025$, or perhaps more realistically $\pm 0.040$, despite having data of $R = 95\,000$ and $S/N \approx 400$ to 600. Moreover, we note that the formal errors which we quote are not the maximum possible values and in some cases will underestimate the real error, particulary since they result from the quadrature sum of several distinct contributions. The continuum from $\sigma_\zeta$ in particular is only a $1\sigma$ value.

#  Discussion

According to models of big bang nucleosynthesis and cosmic-ray spallation, low lithium-6 abundances are predicted in stars of low metallicity. Detections of this isotope in these stars would have thus a strong impact on production models of light elements. All the tests we have done, e.g. comparing the results based on non-corrected and corrected spectra, on different normalisation levels, on different wavelength ranges covered in the analysis, suggest that our best estimates for $^6\mathrm{Li}/^7\mathrm{Li}$ are uncertain to at least 0.025-0.040. Given than our best estimates of $^6\mathrm{Li}/^7\mathrm{Li}$ are of the same order of magnitude as these values, we do not claim any significant detection.

One of the goals of this paper was to observe targets whose lithium-6 content was already determined from observations taken with other spectrographs, and show whether, or not, the published results were corroborated by our observations. There are three objects in common with other studies: , and (color symbols in Figure , *see the electronic edition of the Journal for a color version of the figure*).

The first was observed in and (red symbols in Fig. ). They reported different isotopic ratios results: $^6\mathrm{Li}/^7\mathrm{Li}=0.05\pm0.03$ ($S/N \sim 420$) and $^6\mathrm{Li}/^7\mathrm{Li}=0.01\pm0.013$ ($S/N \sim 680$) respectively. suggested that the differences between their value and that of could be due to the differences in the determination of the broadening parameter. We have used a similar method to that employed by , the difference being that we did not use a rotational profile (metal-poor stars have low rotational velocities). We obtain $\zeta_\mathrm{broad}=6.50$ vs 5.65 $\mathrm{km}\, \mathrm{s}^{-1}$ . We have data of similar $S/N$ ($\sim 642$) but slightly lower resolving power ($R \sim 95\,000\ \mathrm{vs}\ 115\,000$). Our best estimate for this star ($^6\mathrm{Li}/^7\mathrm{Li}$ $= 0.004\pm0.028$) agrees with the result from . Much higher values $^6\mathrm{Li}/^7\mathrm{Li}$ $= 0.074$ would be determined from the non-corrected spectra, but based on the match between the residual features seen in our spectra of this star and the corresponding fast rotator (see Figure , bottom right), we are confident that division (correction) by the fast rotator is essential for . Our observed bright fast rotator star spectra suggest a residual fringing of the order of $0.4-0.6\%$ (peak-to-valley) whereas report a value much lower for theirs $\le 0.2\%$ ($\pm 0.1\%$). The amplitude of the fringe patterns of the HDS spectra (peak-to-valley difference) at 670.0 nm is typically $5\%$, which results in a pattern of a few tenth percent after flat-fielding. The width of the fringe pattern is much larger than the Li absorption feature, and the effect on the line analysis is usually not severe. We note, however, that the effect is not negligible in the setting adopted on May 17th 2005 (see §5.2). According to , their residual fringing is negligible although they do not quote any magnitude so it is difficult to judge whether fringing can account for the differences between the results.

Our best estimate for , the other star in common with (blue symbols in Fig. , *see the electronic edition of the Journal for a color version of the figure*) is $0.047\pm0.039$ which agrees well with their value $0.042\pm0.019$. Our error bar is higher so we do not claim any detection in contrast to them. Concerning the quality of the spectra, we have slightly higher $S/N$ ($\sim 600$ vs 530).

The other star in common with is (green symbols in Fig. ). We have data of similar quality to theirs ($S/N \sim 430$), though their $R$ was a bit higher ($\sim115\,000$). Our best estimate of $0.025\pm0.025$ is higher than their best estimate of $0.00\pm0.02$ but consistent within the uncertainties. If no fast rotation correction was done to our spectrum, a $^6\mathrm{Li}/^7\mathrm{Li}$ value closer to 0.00 would better reproduce the observations. Differences in the assumed stellar parameters for these stars are not significant so they should not affect the results.

Stellar parameters assumed for the stars and are similar to one another. The stars also have spectra of similiar quality. Their isotope abundance ratio results are different but agree within their error bars, $0.047\pm0.039$ vs. $0.004\pm0.028$; the continuum location for in particular is quite uncertain. A change in the continuum level within the uncertaintity for this star would bring the best $^6\mathrm{Li}/^7\mathrm{Li}$ estimates for the two stars closer, $\sim 0.01$ vs $\sim 0.00$, though given that their temperatures differ by $\sim 100$ K, they may have experienced different degrees of $^6\mathrm{Li}$ depletion . The stellar parameters values assumed for and are slightly different and as in the previous case, their isotopic estimates are different but agree within their error bars. The agreement gets better if no correction to the observed spectra is made for . According to our error estimates, the $^6\mathrm{Li}/^7\mathrm{Li}$ ratios of these two stars are the most accurate.

We are not aware of any other lithium-6 determinations in the literature for the two targets and . Unfortunately, the data quality for these two stars were among the lowest of the observed sample. Our best estimate suggests that we can only determine an upper limit of $\sim 0.03$ (1$\sigma$) for . The temperature assumed for this star is the lowest among our sample and low enough, according to the isochrones in , for the star to have suffered from severe $^6\mathrm{Li}$ depletion ($\sim 95\%$ depletion). In fact, published detections of $^6\mathrm{Li}$ at metallicities similar to that of . reported a $^6\mathrm{Li}$ detection for other star of $[\mathrm{Fe}/\mathrm{H}]<-2.5$, LP 815-43 ($^6\mathrm{Li}/^7\mathrm{Li}$$= 0.046\pm0.022$).

The $^6\mathrm{Li}/^7\mathrm{Li}$ values do not correlate with $S/N$; we determined high isotopic ratio values from both high $S/N$ spectra () and low $S/N$ (), and the same applies to low values of the isotopic ratio.

The principle conclusion of this study concerns not the isotope ratio itself but rather the uncertainties associated with the measurements. Despite having data of high resolving power ($R = 95\,000$) and $S/N \sim$ 400 to 600, a thorough examination of error contributions has led us to conclude that, for our dataset, we could not confidently determine the isotope fraction to better than $2.5\%$ to $4.0\%$. Moreover, the detections achieved by other groups for metal-poor stars are around this limit, as we show with the $^6\mathrm{Li}$ contribution functions in Figure  determined from the linear fit to the $^7\mathrm{Li}$ observed abundances. It is quite possible that other studies have achieved a higher precision than was achieved in our study, despite the similarity of the initial datasets. Although, notice that observations show lower levels of fringing than ours. In contrast, we interpret our results as indicating that extreme caution is required in finally quantifying the uncertainties, in order to establish $^6\mathrm{Li}/^7\mathrm{Li}$ ratios to the level of precision required to be confident of detections of $^6$Li at the $4\%$ level or lower. We are not confident that we have significant detections of $^6$Li in any of our five stars.

#  Conclusions

The observed spectra of the resonance line were analysed for five metal-poor stars. Our best fits to the observations suggest that only the star may have a high lithium-6 content, $^6\mathrm{Li}/^7\mathrm{Li}$ $= 0.047 \pm 0.039$. However, the uncertainties of the calculations are comparable to the calculations themselves, such that positive detections cannot be claimed for any of them. In general the analysis is very sensitive to the assumptions made for the continuum location, the residual fringing treatment and even the wavelength range covered in the analysis. Our error bars ($\Delta$$^6\mathrm{Li}/^7\mathrm{Li}$ $\sim 0.025$-$0.040$) tell us that detections may have been possible if the stars had a higher lithium-6 content ($^6\mathrm{Li}/^7\mathrm{Li}$ $\geq$ 0.050-0.080) following a $2\sigma$ criteria.

Not surprising, we observed asymmetries in the profile of the lines as did. However, we have major concerns because the magnitude of the asymmetries depends on the spectral normalisation carried out, and the presence of residual fringing could also introduce spurious asymmetries.

Our $^6\mathrm{Li}/^7\mathrm{Li}$ results are comparable to recent values in the literature, but our error estimates are higher, the result of a detailed analysis of a range of uncertain aspects of the analysis associated with the significant presence of structure in our reduced data-set (summarised in § ). In contrast, tabulated only the uncertainties due to finite $S/N$ —which they find usually to be their dominant error— but not the wide range of error sources encountered in our work. This puts limitations on the detections so that only upper limits of $^6\mathrm{Li}/^7\mathrm{Li}$ were inferred for the stars observed in this work.

One possible, but unjustified, generalisation of our results would be to question the reliability of all $^6\mathrm{Li}/^7\mathrm{Li}$ measurements at levels of 4% or 5%. Making this generalisation would then call into question the evidence for a possible upper envelope or plateau in the $^6$Li abundance, which has caught the interest of astrophysicists in recent years. We believe there is a value in continuing research into the $^6\mathrm{Li}/^7\mathrm{Li}$ ratios in metal-poor stars, including the fullest treatment of possible error sources, to determine whether such an upper envelope or plateau exists.

[^1]: $\log\epsilon_X = \log{N_{X}/N_\mathrm{H}} + 12$


author: Haroon Asghar; John. G. McInerney
title: Asymmetric dual-loop feedback to suppress spurious tones and reduce timing jitter in self-mode-locked quantum dash lasers emitting at 1.55 $\mu$m

Semiconductor mode-locked diode lasers (MLLs) are compact, rugged and efficient sources of ultra-short, intense, and high repetition frequency optical pulses with many applications such as, all-optical clock recovery, Lidar, optical frequency combs, and telecommunications \[1-3\]. A major limitation of MLLs for most practical applications is their very high timing jitter and phase noise, as spontaneous emission noise and cavity losses make MLLs prone to broad linewidths and therefore substantial phase noise \[4\]. To improve timing jitter, several experimental methods such as single cavity feedback \[5-8\], coupled optoelectronic oscillators (OEOs) \[9\], injection locking \[10-12\] and dual-loop feedback \[13-15\] have been proposed and demonstrated. Of the stabilization techniques demonstrated to date, optical feedback is a promising approach in which an additional reflector creates a compound cavity with a high quality factor, with no need for an external RF or optical source. Due to the existence of the extra mirror, side-bands resonant with the round trip time of the external cavity are generated which affect the overall timing jitter and quality of the RF spectra. To overcome these issues, optoelectronic feedback \[9\] can also be utilized to stabilize timing jitter and to suppress cavity side-modes by conversion of the optical oscillation (using a fast photodetector) to an electrical signal used in a long feedback loop. This technique does not require an RF source, but requires optical-to-electrical conversion. Recently, a simpler dual-loop feedback technique \[13-15\] without optical/electrical conversion has been demonstrated to improve timing jitter of the MLLs and to filter or suppress the unwanted spurious side-bands. Dual-loop configurations proposed to date \[13, 14\] yields sub-kHz linewidth but produces additional noise peaks at frequencies resonant with the inverse of the delay time in the second cavity. This is undesirable in many applications where low noise and flat spectra are required, as in frequency comb generation. Recently the influence of the second feedback delay on side-mode suppression \[16\] and timing jitter \[17\] has been studied numerically. In this letter, we report experimental investigation to eliminate these adverse dynamical effects using asymmetric dual-loop feedback by appropriately choosing the length of second feedback cavity. Best side-mode suppression and lower timing jitter relative to single loop feedback is achieved with the length ratio between the two cavities $\sim$ 8x. It was further observed that RF linewidth and integrated timing jitter were reduced by increasing the length of the second cavity. Our findings suggest that noise stabilization and side-mode suppression depends strongly on additional feedback delay times.

Devices under investigation are two-section InAs/InP quantum dash mode-locked lasers (QDash MLLs) whose active layers have 9 InAs quantum dash monolayers grown by gas source molecular beam epitaxy (GSMBE) embedded within two barriers and separate confinement heterostructure (SCH) layers (dash in a barrier structure). Both barriers and SCHs consisted of $In_{0.8}Ga_{0.2}As_{0.4}P_{0.6}$ quaternary materials with $\lambda_{g}$ = 1.55 $\mu$m \[18\]. Total cavity length was 2030 $\mu$m with absorber lengths 240 $\mu$m (length ratio $\sim$ 11.8%), giving pulse repetition frequency $\sim$ 20.7 GHz ($I_{Gain}$ = 300 mA) and average free-space output powers of a few mW. Mode-locking was obtained without reverse bias applied to the absorber section, and the heat sink temperature was fixed at $19^0$C. This is a two-section device but works similarly to a single section self-mode-locked laser, since the absorber is not biased; in this case the amount of minimal residual absorption does not affect the mode-locking mechanism \[12\]. The absorber and gain sections were isolated by a resistance 9 k$\Omega$. The QDash MLL was mounted p-side up on AlN submounts and copper blocks with active temperature control and electrical contacts formed by wire-bonding.

A schematic for the dual-loop technique is depicted in Fig. 1. For single and dual-loop feedback configurations, a calibrated fraction of light was fed back through port 1 of an optical circulator, then injected into the laser cavity via port 2. Optical coupling loss from port 2 to port 3 was -0.64 dB. The output of the circulator was sent to a semiconductor optical amplifier (SOA) with gain of 9.8 dB, then split into two arms by a 50/50 coupler. 50% went to an RF spectrum analyzer (Keysight E-series, E4407B) via a 21 GHz photodiode and to optical spectrum analyzers (Ando AQ6317B and Advantest Q8384). The other 50% of power was split into two equal parts by a 3-dB splitter. For single loop feedback, all power passed through loop-I. For dual-loop configurations (feedback loops-I and-II) the power was split into two loops at the 3-dB splitter. Feedback strengths in both loops were controlled by variable optical attenuators and monitored using power meters. In this experimental arrangement, the length of loop-I was fixed to 160m while the length of second feedback loop was varied in three chosen lengths: 20, 53 and 80m. Polarization controllers in each loop plus one polarization controller before port 1 of the circulator ensured the light fed back through both loops matched the emitted light polarizations to maximize feedback effectiveness. In this experiment the feedback ratio into gain section was limited to $\sim$ -22 dB.

We calculated RMS timing jitter from single sideband (SSB) phase noise spectra measured at the fundamental RF pulse repitition frequency (20.7 GHz) using: $$\sigma_{RMS}=\frac{1}{2 \pi f_{ML}}\sqrt{2 \int_{f_{d}}^{f_{u}} L(f)\,df}$$where $f_{ML}$ is the pulse repetition rate, $f_{u}$ and $f_{d}$ are the upper and lower integration limits. *L(f)* is the single sideband (SSB) phase noise spectrum, normalized to the carrier power per Hz. To measure RMS timing jitter of the laser in more detail, single-sideband (SSB) noise spectra for the fundamental harmonic repetition frequency were measured. To assess this, RF spectra at several spans around the repetition frequency were measured from small (finest) to large (coarse) resolution bandwidths. The corresponding ranges for frequency offsets were then extracted from each spectrum and superimposed to obtain SSB spectra normalized for power and per unit frequency bandwidth. The higher frequency bound was set to 100 MHz (instrument limited).

To observe the RF spectrum with single-loop feedback, the length of the loop was initially set at 160m; optimally stable resonance occurred when the feedback length was fine tuned using an optical delay line (ODL-I) (which spanned 0 to 84 ps in steps of 1.67 ps). Such optimization provides a resonant condition (at delay setting=13 ps) under which the RF linewidth was reduced from 100 kHz free-running to as low as 4 kHz, with integrated timing jitter to 0.7 ps from 3.9 ps \[integrated from 10 kHz - 100 MHz\]. Measured phase noise traces as functions of frequency offset from fundamental mode-locked frequency and RF spectra for free running laser (green line) and single loop feedback (gray line) are given in Fig. 5 and Fig. 6, respectively. Under similar delay settings, external cavity side-modes appear in the RF spectrum with frequency spacing 1.28 MHz, the inverse of the loop round trip delay. RF spectra are shown in Fig. 2(a) (gray line) and (b) (gray line), using spans of 10 MHz and 100 MHz, respectively. Frequency resonances can be seen in both frequency spans which contribute significantly to timing jitter, particularly for the longer feedback cavities as they are closer to the main peak and are less suppressed \[11\]. To eliminate these fluctuations and to improve the side-mode suppression ratio (SMSR), dual-loop feedback was implemented as described in the next section.

To assess the suppression of these frequency resonances, a shorter feedback cavity corresponding to half the period of noise-induced oscillations of loop-I was introduced. Feedback strengths of both cavities were equalized using variable optical attenuators (Att-I and Att-II) plus polarization controllers (PC-I and PC-II) in both loops. One optical delay (ODL-II) was adjusted to full resonance and the length of the other delay line (ODL-I) was tuned over the maximum range available 0-84 ps. When the optical delay lines ODL-I and ODL-II were fine tuned (ODL-I=15 ps and ODL-II=25 ps) so that every second mode of loop-I coincided precisely with a mode of loop-II, maximum \>30 dB suppression in the first order side-mode was achieved. However, additional fluctuations (modal overlaps) appeared at frequencies resonant with the inverse of the length of second delay time which becomes the carrier signal. These noise fluctuations depends on the ratio of the loop lengths. Here these fluctuations at frequency spacing 2.60 MHz, are consistent with the length of the second feedback loop 80m. RF spectra of the asymmetric dual-loop configuration are shown in Fig. 2(a) (red line) and (b) (red line), using spans 10 MHz and 100 MHz, respectively. In this fully resonant configuration, the RF linewidth narrowed to \< 1 kHz (instrument limited), with timing jitter reduced to 295 fs. Phase noise trace and RF spectra are shown in the Fig. 5 (blue line) and Fig. 6 (blue line), respectively. The RF spectra illustrated in Fig. 2(a) and (b) show that this feedback configuration is not suitable to achieve effective suppression in frequency resonances, as the second delay time will be resonant with the second mode of the first feedback loop which restrict many practical applications where flat and side-band free RF spectra are required. To improve on this situation, a different dual-loop feedback configuration with non-resonant shorter second loop (53m) was investigated, described in the next section.

In the dual-loop configuration presented in this section, the length of loop-I was initially set to 160m while that of loop-II was 53m. Upon fine tuning of both optical delay lines (ODL-I=13 ps and ODL-II=15 ps), when the second delay time was resonant with the third harmonic of the first loop, suppression of the first two frequency resonances occurred, while the third harmonic (modal overlap) was unsuppressed. This harmonic was observed at frequency offset 3.9 MHz, corresponding to the 53m length of the outer feedback loop. RF spectra for the asymmetric dual-loop configuration are shown in Fig. 3(a) (red line) and (b) (red line), using spans 10 MHz and 100 MHz, respectively. In this configuration, when both external cavities are fully resonant, the RF linewidth narrows to 2 kHz with integrated timing jitter 0.45 ps. Phase noise trace and RF spectra are shown in Fig. 5 (black line) and Fig. 6 (black line) respectively. These experimentally measured results show that external cavity side-modes cannot be optimally suppressed by simply choosing the second feedback delay time to be a fraction of the first. To achieve stable and flat RF spectra, we designed an asymmetric dual-loop feedback configuration for effective suppression of external cavity side-modes. Which produced flat RF spectra close to the main peak compared to conventional single- and dual-loop feedback.

In this asymmetric dual-loop feedback configuration the length of loop-I was fixed (160m) and loop-II was set $\sim$ 8x shorter than loop-I. Fine tuning of both cavities (ODL-I=15 ps and ODL-II=21 ps) produced precise coincidence of every eighth mode of loop-I with a mode of loop-II, so that strong side-mode suppression occurred and all feedback-induced side-modes and spectral resonances were eliminated under frequency span 10 MHz. RF spectra for this dual-loop feedback configuration (red line) are shown in Fig. 4(a) and (b) with frequency spans 10 MHz and 100 MHz respectively. Furthermore, when both external cavities are fully resonant, the RF linewidth narrows to 8 kHz with integrated timing jitter 0.6 ps. The phase noise trace and RF spectra are shown in Fig. 5 (red line) and Fig. 6 (red line) respectively. In this configuration, the RF linewidth was higher than for single loop feedback, but measured timing jitter was lower; this is due to suppression of external cavity side-modes. In addition, weak modal overlap (with intensity $\sim$ -6 dBm) in RF spectra of dual-loop feedback was noticed at 10.2 MHz frequency spacing, consistent with our 20m outer loop; this is shown in Fig. 4(b)(red line). This behavior shows that effective suppression of external cavity side-modes and reduced timing jitter can be achieved by appropriately fine-tuning the length of the second feedback loop. It should be noted that length of loop-II ($\sim$ 20m) is only optimal in our specific experimental setup. Further reduction in the length of second feedback loop is not possible, as the combined variable optical attenuator, optical delay line, polarization controller and 3-dB coupler have minimum length of $\sim$ 20m. Better suppression of external cavity side-modes could be achieved in an arrangement not subject to this limitation, such as photonic integrated circuit.

Measured phase noise traces for free-running condition (green line), single-loop (gray line) and dual-loop feedback with loop-II at 20m (red line), 53m (black line) and 80m (blue line) as functions of frequency offset from the fundamental mode-locked frequency are given in Fig. 5

Comparison of RF linewidth and integrated timing jitter under stable resonant conditions, for three chosen lengths of second feedback cavity is shown in inset of Fig. 6. Measured integrated timing jitter in all dual-loop configurations was lower than for single-loop feedback. However, best suppression in external cavity side-modes was achieved with the second delay time $\sim$ 8x shorter than the first. Furthermore, the integrated timing jitter in this case was 16% lower than for single-loop feedback. Reduction in timing jitter occurs due to suppression of external cavity side-modes relative to single-loop feedback.

In the literature \[13, 14\], it was experimentally observed that side-mode suppression was achieved when both feedback delays had a common multiple. This shows that effective suppression in external cavity side-modes is highly dependent on the length of the second loop. Recently, the influence of the second loop delay on suppression of external cavity side-modes \[16\] and timing jitter \[17\] was studied numerically. In this work, experimentally measured suppression of cavity side-modes and integrated timing jitter as a function of the second cavity delay correspond well with published numerical simulations \[16, 17\].

In summary, an asymmetric dual-loop feedback method has been demonstrated which suppress additional noise resonances found in conventional single- and dual-loop feedback schemes. These results show that dual-loop feedback with precise alignment of the second loop delay effectively suppresses external cavity side modes and produces flat RF spectra closer to the main peak. Furthermore, by increasing the length of the second loop, significant reduction in RF linewidth and integrated timing jitter was produced. Our experimental results have validated recently published numerical simulations. Using this method, stable side-band-free integrated photonic oscillators based on mode-locked lasers may be developed which are feasible and attractive for many applications in optical telecommunications, time-domain multiplexing, frequency comb generation and as synchronized pulse sources or multi-wavelength lasers for wavelength-diversity or multiplexing.

**Funding.** The authors acknowledge financial support from Science Foundation Ireland (grant 12/IP/1658) and the European Office of Aerospace Research and Development (grant FA9550-14-1-0204).