\section{Introduction} \label{sec:intro}
UAVs designs can be divided into two principal types: fixed-wing and the rotatory aircraft, both have their respective pros and cons. Nonetheless, it is possible to obtain certain capabilities of both of them by using the so-called convertible aircraft. There exist two main configurations of convertible aircraft: tilt-rotors \cite{1,7040348, GFlores:thesis, flores_cuav_jirs, Flores:IFAC14_convertible} and tail-sitter \cite{4,9,10}. The tail-sitter transition maneuver is considerable more complex than tilt-rotors, since it is required that the complete UAV body tilt around one axis for achieving the \textit{transition maneuver}, an ordered movement that carries the UAV from hover to cruise mode and vice versa, see Fig. \ref{fig:cambio}. This paper is focused to develop and implement a controller especially dedicated to achieve such transition maneuver.

Most of the papers related with this topic tackle the problem from a simplified model without taking into account aerodynamics, bounded limits of actuators, or even smoothness of the controller. In \cite{1,5,14} the control strategy consists in increasing the motor angles until they are horizontally aligned, then a switching strategy between flight modes is performed. In \cite{19} authors develop a robust algorithm to control airplane and hover flight modes considering aerodynamic terms, however it does not present a transition control. In \cite{16} the transition is performed manually and the control relies in modifying the flight modes control weights according to the percentage of the transition progress. In \cite{15,8444168,8327537,LI2018348} it is implemented a controller in for the transition maneuver of a tail-sitter UAV, the controller is designed in the 6-DOF, on the other hand, however \cite{8431535} lacks of a theoretical proof which demonstrate the effectiveness of the transition maneuver. Also, several unified controllers which avoids control switching methods are investigated in \cite{8206359,8463158}.
\begin{figure}[ht!]
	\centering
	\includegraphics[width=250pt,page=1]{graficos.pdf}
	\caption{Transition maneuver, graphical description.}
	\label{fig:cambio}
\end{figure}

In this paper we present a mathematical UAV model which considers aerodynamic effects. The proposed controller takes into account the usual available states in a real scenario. Also, the control design avoids switching dynamics between modes, having a sweet and bounded controllers, which are desired properties in real implementations. The controller designed is based on the time-scale separation property presented in attitude and position UAV dynamics \cite{20,21}. The control design is used saturation functions and Lyapunov stability.

The remainder of this paper is organized as follows. In Section \ref{sec:problem} system model and problem statement are presented. Section \ref{sec:result} presents the control system proposed and the followed strategy to ensure a successful flight mode transition. In Section \ref{sec:sims} simulations experiments demonstrate the effectiveness of the theoretical results. Finally, some conclusions and further work are described in Section \ref{sec:conclusions}.
\section{Problem Formulation} \label{sec:problem}
The tail-sitter UAV have three operation modes: a) hover; b) cruise; and c) transition mode. The last mode consists in the transition phase between cruise (hover) to hover (cruise); this flying mode will be investigated in this paper. In hover to cruise mode the attitude changes from $90$ to approximately $6$ degrees around between the $y$ axis, which is the range for cruise flight mode. It is the same transition from cruise to hover mode but with the angle reversed.
\begin{figure}[ht!]
	\centering
	\includegraphics[width=180pt,page=2]{graficos.pdf}
	\caption{Tail-sitter aircraft; body and world coordinate frames.}
	\label{frames}
\end{figure}
The drone is composed by six actuators: four rotors and two servos which manipulate attitude dynamics. These six actuators results in an over-actuated roll and pitch subsystems. Figs. \ref{fig:cambio} and \ref{frames} depict the tail-sitter drone used in this study.

\subsection{System Model}
The analysis and control design is implemented on the $(x,z)$ plane from the body frame. This is the plane where the most interesting dynamics for the transition maneuver appears \cite{6,8}. Forces affecting aircraft behavior are shown in Fig. \ref{forces}. The system can be described as follows
\begin{figure}[ht!]
	\centering
	\includegraphics[width=180pt]{forces1.pdf}
	\caption{System forces taken into account for mathematical modeling and control design.}
	\label{forces}
\end{figure}
\begin{eqnarray}
\label{sum1}
&\Sigma_{1}&\begin{cases}
        	\dot{u} = -D \cos\alpha + L\sin\alpha + T - g\sin\theta - qw\\
        	\dot{w} = -D \sin\alpha - L\cos\alpha + g\cos\theta + qu
	\end{cases}\\
\label{sum2}
    &\Sigma_{2}& \begin{cases}
        	\dot{\theta} = q\\
        	\dot{q} = \tau.
	\end{cases}
\end{eqnarray}
Model variables and parameters are defined in Table \ref{tab:variables}.
\begin{table}[!hbt]
\begin{center}
\begin{tabular}{|| l | l ||}
\hline
\hline
$u$ & $x$ body velocity expressed in body frame \\
\hline
$w$ & $y$ body velocity expressed in body frame \\
\hline
$\theta$ & Aircraft pitch angle \\
\hline
$q$ & Aircraft angular velocity \\
$D$ & Drag force \\
\hline
$L$ & Lift force \\
\hline
$\alpha$ & Angle of Attack (AoA) \\
\hline
$g$ & Gravity \\
\hline
$T$ & Thrust force (Input control)\\
\hline
$\tau$ & Torque (Input control)\\
\hline
\end{tabular}
\caption{System model variables and parameters.}
\label{tab:variables}
\end{center}
\end{table}
The angle of attack (AoA) is defined as
\begin{equation}
\label{eqn:alpha}
\alpha = \tan^{-1}\frac{w}{u}.
\end{equation}
Control inputs for system (\ref{sum1}-\ref{sum2}) are the thrust $T$ and the torque $\tau$ generated by the difference of thrust between rotors and the aileron angle. These terms can be described as
\begin{eqnarray}
	\label{eqn:thrust}
	T &=& T_{1} + T_{2} \\
    \label{eqn:tau}
    \tau &=& \tau_{\delta} + \tau_{L} + (T_{1} - T_{2})
\end{eqnarray}
where $\tau_{\delta}$ is the part of torque generated by the propellers wind stream in the aileron; $\tau_{L}$ is the torque induced by wind stream shocking in the wing area, and the rest of the torque is originated by the difference in force between the two pairs of motors. The total thrust force is calculated by the sum of each motor's thrust.

System (\ref{sum1})-(\ref{sum2}) can be divided in two subsystems, where $\sum{}{}_{1}$ and $\sum{}{}_{2}$ represent translational and attitude dynamics, respectively. From \cite{20,21}, it is known that $\sum{}{}_{2}$ is faster than $\sum{}{}_{1}$, hence $\theta$ can be used as a virtual controller for $\sum{}{}_{1}$. Subsystem $\sum{}{}_{2}$ is a double integrator and can be fully controlled by $\tau$.

\subsubsection{Aerodynamic forces}
The proposed model (\ref{sum1})-(\ref{sum2}) considers aerodynamic terms explained next. The \textit{lift} $L$ and \textit{drag} $D$ forces are nonlinear functions of several variables. The most important of these variables are: airspeed, AoA and lift and drag coefficients \cite{7}. Lift and drag are defined as follows
\begin{eqnarray}
	\label{eqn:Lforce}
	L &=& KC_{L}V^{2} \\
    \label{eqn:Dforce}
	D &=& KC_{D}V^{2}
\end{eqnarray}
where $K = \frac{\rho S}{2}$; $S$ is the wing hitting area and $\rho$ is air pressure; $V=\sqrt{u^{2} + w^{2}}$ is air speed magnitude. The drag and lift coefficients $(C_{L}, C_{D})$ for the airfoil used in this experiment are described at Fig. \ref{wing_forces} according to the angle of attack $\alpha$. Such coefficients correspond to a symmetrical airfoil NACA-0020 shown at Fig. \ref{airfoil}. They were obtained by numerical simulations made in the software XFR5 with many different Reynold numbers. Several aileron angles $\delta$ are considered to compute different lift and drag coefficients.
\begin{figure}[ht!]
	\centering
	\includegraphics[width=240pt, height=300pt]{cl_cd.pdf}
	\caption{Lift and Drag coefficients w.r.t $\alpha$ angle obtained from XFR5 software and studies presented in \cite{18}. According to the maximum lift to drag ratio $(L/D)$ considering $C_{l}$ and $C_{d}$, it is obtained an AoA $\alpha=6$ $degrees$ tha corresponds to the optimal AoA for cruise flight mode.}
	\label{wing_forces}
\end{figure}
\begin{figure}[ht!]
	\centering
	\includegraphics[width=220pt]{airfoil.pdf}
	\caption{Normalized airfoil profile used in the aircraft (NACA-0020).}
	\label{airfoil}
\end{figure}

\subsection{Problem Statement}
The problem can be summarized in the following
\begin{problem}
Given system (\ref{sum1})-(\ref{sum2}), design a continuous and easy-to implement control with bounded inputs that stabilizes the transition maneuver from hover to cruise mode and vice-versa.
\end{problem}
\section{Main Result} \label{sec:result}
The key idea behind the proposed controller is explained next. Based on time-scale separation principle presented in \cite{20,21}, it is possible to use the variable $\theta$ as a virtual controller for subsystem \eqref{sum1}. Once this virtual controller is designed, its value must be tracked by the subsystem \eqref{sum2}.

Let's start with the design of the virtual controller for $\sum_{1}$. In this case the states $u$ and $w$ should converge to predefined values according to the appropriate characteristics for the actual flight mode. Subsystem \eqref{sum1} can be rewritten as
\begin{eqnarray}
	\label{eqn:udot3}
    \dot{u} &=& f_{1}(u,w,q) + T - g\sqrt[2]{1-\epsilon^{2}}\\
    \label{eqn:wdot3}
    \dot{w} &=& f_{2}(u,w,q) + g\epsilon
\end{eqnarray}
where
\begin{eqnarray*}
	\epsilon &=& \cos\theta\\
    f_{1}(u,w,q) &=& -D \cos\alpha + L\sin\alpha - qw\\
	f_{2}(u,w,q) &=& -D \sin\alpha - L\cos\alpha + qu.
\end{eqnarray*}
From the previous discussion $\epsilon$ is considered as a controller, since it is a function of $theta$ angle. Before proposing the control algorithm we should define some natural assumptions and remarks:
\begin{assumption}
\label{ass:1}
In cruise flight mode it is natural to consider that $u \gg w$, i.e. the horizontal velocity is much greater than vertical velocity.
\end{assumption}
\begin{assumption}
\label{ass:2}
In hover mode $u\approx 0$ and $w\approx 0$.
\end{assumption}
\begin{remark}
\label{rk:alpha}
From \eqref{eqn:alpha}, $u \neq 0$ .
\end{remark}
\begin{remark}
\label{rk:epsilon}
Since \eqref{eqn:udot3} holds, then $\mid \epsilon \mid \leq 1$.
\end{remark}

One important aspect of the control strategy is to design desired values of aircraft velocities $(u_{d}, w_{d})$ according to the speed characteristics of each flight mode. If the transition maneuver is from hover to cruise, $u$ must increase to a value in which the lift force is enough to maintain the aircraft at a constant altitude. On the other hand, $w$ must be designed w.r.t. $u$ in such a way that the aircraft AoA be at the optimal value considering the lift to drag ratio $(L/D)$, which for the lift and drag coefficients shown at Fig. (\ref{wing_forces}) results in $6$ degrees. Equations for $(u_{d}, w_{d})$ are computed as follows
\begin{eqnarray}
\label{eq:ud}
&u_{d}=&\begin{cases}
        	\frac{\arctan \big(a_{u}(\frac{t}{5}-L_{u})\big) }{a_{u}}+L_{u} \hspace{12pt} $if$ \hspace{10pt} \frac{t}{5} > L_{u}\\
            \frac{t}{5}  \hspace{88pt} $if$ \hspace{10pt} \frac{t}{5} \leq L_{u}
	\end{cases}\\
    \label{eq:al}
&\alpha_{d}=&\begin{cases}
        	\frac{\arctan \big(a_{\alpha}(t-L_{\alpha})\big) }{a_{\alpha}}+L_{\alpha}  & $if$ \hspace{10pt} t > L_{\alpha}\\
            t  & $if$ \hspace{10pt} t\leq L_{\alpha}
	\end{cases}\\
    \hspace{10pt}
\label{eq:wd}
    &w_{d}=& u_{d}\tan\alpha_{d}
\end{eqnarray}
where $a_{u} = \frac{\pi}{2(M_{u}-L_{u})}$, $a_{\alpha} = \frac{\pi}{2(M_{\alpha}-L_{\alpha})}$; $(M_u, L_u)$, $(M_{\alpha}, L_{\alpha})$ are positive constants which $L_u \leq M_u$, $L_{\alpha} \leq M_{\alpha}$ and $M_u$, $M_{\alpha}$ are the maximum value the functions could have and $L_u$, $L_{\alpha}$ are the linear region of the functions. For determining the desired AoA from hover to cruise mode, it is supposed that the system starts from $\alpha=0$ and ends in $\alpha=6$ degrees which can be calculated with \eqref{eq:al}. It is important to mention that $(u_d,w_d,\alpha_d)$ are continuously differentiable functions. In the same way, when the transition occurs from cruise to hover, both speeds must be reduced in such a way that the AoA decreases to the point of being equal to zero.

From Remark \ref{rk:epsilon}, the virtual control $epsilon$ is designed by using a continuous and non-decreasing saturation function $\sigma(s)$ satisfying \cite{23}:
\begin{enumerate}[(a)]
\item $s\sigma(s) > 0$	\quad\quad $\forall s \neq 0$.
\item $\mid\sigma(s)\mid \leq M_s$ \quad\space $\forall s \in \mathbb{R}$.
\item $\sigma(s) = s$ when $\mid s \mid \leq L_s$.\\
\end{enumerate}
where $L_s$ and $M_s$ are positive constants which describes the linear part of the function and the saturation value, respectively, hence $L_s \leq M_s$. It is now possible to propose a preliminary control system $(\epsilon, T)$ as follows
\begin{eqnarray}
	\label{eqn:epsilon}
    \epsilon &=& -\sigma_{2}\big(f_{2}(u,w,q) + \sigma_{1}(w - w_{d}) - \dot{w}_{d}\big)\\
    \label{eqn:T}
    T &=& -\sigma_{3}(u - u_{d}) + g\sqrt[2]{1-\epsilon^{2}} - f_{1}(u,w,q) + \dot{u}_{d}.
\end{eqnarray}
From (\ref{eqn:epsilon}) it is possible to obtain $\theta_{d}$. Then it is possible to define controller
\begin{equation}
        \label{control:tau}
      \tau = -k_{\theta}(\theta - \theta_{d}) - k_{q}(q - q_{d}) + \dot{q}_d
\end{equation}
The main result can be summarized in the following
\begin{theorem}
Let system \eqref{sum1}-\eqref{sum2} and consider Assumptions (\ref{ass:1},\ref{ass:2}), then controllers \eqref{eqn:T}, \eqref{control:tau} with $\theta_{d}$ obtained from \eqref{eqn:epsilon} results in a local asymptotically stable closed-loop system.
\end{theorem}
\begin{proof}
Consider a coordinate change over the original system with
\begin{eqnarray*}
	x_{1} = u - u_{d} \hspace{15pt}  x_{2} = w - w_{d} \hspace{15pt} x_{3} = \theta - \theta_{d} \hspace{15pt} x_{4} = q - q_{d}.
\end{eqnarray*}
Since $x_{3}$ and $x_{4}$ states are independent of $x_{1}$ and $x_{2}$, and because it is well known that the structure of the subsystem $\sum_{2}$ is at least GAS, then the problem reduces to investigate stability of states $(x_1,x_2)$. For that, observe that the closed-loop system has a structure similar to that of \cite{23}, and then with the Lyapunov function
\begin{equation}
	\label{eqn:lyapunov2}
    V = \frac{1}{2}x_{1}^{2} + \frac{1}{2}x_{2}^{2}.
\end{equation}
and following the same procedure of \cite{23} one can obtain asymptotic stability of the closed-loop system.
\end{proof}
\section{Simulation Experiments} \label{sec:sims}
In this section a MATLAB simulation is performed, where two experiments are conducted: transition from hover to cruise mode and vice versa. To define the $u_{d}$ and $w_{d}$, the constants used in \eqref{eq:ud}, \eqref{eq:al} and \eqref{eq:wd} were $L_{u}=0.7$, $M_{u}=1$, $L_{\alpha}=4$, $M_{\alpha}=6$ and the saturation and gain constants for the controls were $k_{\theta}=5$, $k_{q}=3$, $L_{\sigma}=0.9$ and $M_{\sigma}=1$. Results are explained next.

\subsection{Hover to cruise flight mode}
During hover, the initial conditions were $u_{0}=0.18$, $w_{0}=0.03$ and $\theta_{0}=80$ degrees. Results are shown in Fig. \ref{fig4} where the evolution over time of $\theta$, $\alpha$ and $\tau$ are presented. Fig. \ref{fig5} depicts the thrust control as $u$ and $w$ evolve over time. According to the conditions established for the hover mode, at the end of the transition it is achieved that $u \gg w$. This is due to the conditions of cruise flight mode. Since the velocity at $u$ is much higher than $w$, it leads to a relatively small AoA as shown in the middle plot of Fig. \ref{fig4}, generating in the same way a lift force in the drone. This allows that the thrust needed to hold the aircraft in the air can be smaller as shown in the lower plot of Fig. \ref{fig5}.

\subsection{Cruise to hover flight mode}
The simulation results for the opposite flight transition are presented with $u_{0}=0.7$, $w_{0}=0.2$ and $\theta_{0}=35$ degrees as the initial conditions. Fig. \ref{angulos2} shows the system evolution w.r.t. the desired angle and its respective AoA. In Fig. \ref{thrust2} it can be seen that the velocities decrease at $\approx0$ which corresponds to the hover condition. Since these speeds decrease, the lift and drag forces decrease in the same way as the angle of attack approaches to zero, then, to keep the drone flying, thrust force will be increased again as can be seen in the third plot of Fig. \ref{thrust2}.

\begin{figure}[ht!]
	\includegraphics[scale=0.65]{hover-cruise.pdf}
	\centering
	\caption{$\theta$ and $\alpha$ evolution from hover to cruise flight mode according to $\theta_{d}$; also the control input $\tau$ is presented.}
	\label{fig4}
\end{figure}
\begin{figure}[ht!]
	\includegraphics[scale=0.65]{hover-cruise1.pdf}
	\centering
	\caption{Velocities $u$ and $w$ during the transition from hover to cruise flight mode, also the thrust control is presented.}
	\label{fig5}
\end{figure}
\begin{figure}[ht!]
	\includegraphics[scale=0.65]{cruise-hover.pdf}
	\centering
	\caption{$\theta$ and $\alpha$ angle evolution from cruise to hover flight mode according to $\theta_{d}$; also the control input $\tau$ is presented.}
	\label{angulos2}
\end{figure}
\begin{figure}[ht!]
	\includegraphics[scale=0.65]{cruise-hover1.pdf}
	\centering
	\caption{$u$ and $w$ velocities during the transition from cruise to hover flight mode, with the respective thrust control.}
	\label{thrust2}
\end{figure}
\section{Conclusions} \label{sec:conclusions}
In this paper, a simple control strategy for the transition maneuver of the tail-sitter UAV is proposed. Such controller is based on the time-scale separation and saturation functions. The design is based on Lyapunov approach. Simulations demonstrate the effectiveness of the controller for achieving transition from hover to cruise mode and vice versa. It is important to mention that the controller have the peculiarity that it does not present any switching, it is smooth and it takes into account the saturation limits imposed by the actuators. Such characteristics are useful for implementation in a real UAV.

In a further work, we can execute a real experiment with the physical prototype of the UAV and investigate the behavior of the real system. Also it could be possible to design a robust control that consider external disturbance as part of the modeling, such as wind gusts terms. Also the control can be designed by analyzing the 6-DOF model.

\bibliographystyle{IEEEtran}


\section{Introduction}
\label{sec:intro}

 Nowadays, optical atomic clocks reach fractional uncertainties at the \SI{2 e-18}{} level \cite{Nicholson2015} and relative instabilities at the \SI{e-17}{} level at one second of averaging time \cite{Sanner2018}. Due to their high accuracy and stability, optical clocks are used in various applications in fundamental research, pushing the limits in tests of the time variation of fundamental constants \cite{Safronova2018b}, tests of general relativity \cite{Dzuba2017}, search for dark matter \cite{Roberts2017,Kalaydzhyan2017} and are envisioned as detectors for gravitational waves \cite{Kolkowitz2016}.  Their high precision makes optical clocks also  important tools for applications like navigation \cite{Schiller2012} and geodesy \cite{Mehlstaubler2018,Flury2016}.

Currently the best performing optical clocks are based on neutral atoms trapped in a magic wavelength optical lattice \cite{Ludlow2015}.  In particular, clocks based on neutral strontium have become widely used over the world in the past, operating uninterrupted over multiple weeks in laboratories and recently started to even perform outside the lab \cite{Grotti2018}.

\begin{figure}[b]
\centering
\includegraphics[width=0.95\linewidth]{Strontium_levelscheme.png}
\caption{The level scheme for the bosonic strontium isotope $^{88}$Sr with  transitions used for laser cooling, repumping, clock operation and detection and their vacuum wavelengths (linewidth given for the cyan transitions).}
\label{fig:Sr_scheme}
\end{figure}

However, the general scheme for laser cooling and trapping of neutral strontium and other alkaline earth atoms is more complex compared to the one used for alkali atoms \cite{Metcalf2003}. The cooling process starts with a first cooling stage utilizing the broad $^1S_0 \rightarrow$  $^1P_1$ transition at \SI{460.9}{nm}, from where atoms can decay towards the meta-stable $^3P_2$ state with a branching of roughly 1 in 50,000 (see Figure \ref{fig:Sr_scheme}),  which limits the MOT lifetime to a few tens of milliseconds \cite{Stellmer2014}. To increase the number of atoms in the MOT several repumping schemes have been employed. Due to the availability of laser sources the majority of strontium lattice clocks use two repump lasers addressing the $^3P_2 \rightarrow$  $^3S_1$ transition at \SI{707}{nm} and the $^3P_0 \rightarrow$  $^3S_1$ transition at \SI{679}{nm}.
Alternatively, only a single repump laser addressing the $^3P_2 \rightarrow$  $^3D_2$ transition at \SI{496.93}{nm} may be used \cite{Moriya2018,Poli2005}. Until now the generation of light at this wavelength relied on second harmonic generation (SHG) from an infrared laser due to the lack of laser diodes directly operating in this range. This results in complex setups with low optical output power but high power consumption that are undesirable for mobile applications.

Newly available gallium-nitride (GaN) laser diodes operating around \SI{495}{nm} pave the way to build extended-cavity diode lasers (ECDL) at \SI{497}{nm}. This allows for more compact and energy efficient laser sources for laser cooling and trapping of strontium.  Optical powers of up to \SI{100}{mW} open up the possibility of sub-Doppler cooling on the almost closed cycling transition $^3P_2 \rightarrow$  $^3D_3$ \cite{Stellmer2014,Hayakawa2018} as already successfully demonstrated with the alkaline earth metal calcium \cite{Gruunert2002}. This could lead to temperatures in the $\mu$K regime, sufficient for effectively loading atoms into the magic wavelength lattice and therefore allows to replace the usual laser system at \SI{689}{nm} that requires prestabilization to an optical cavity.

Here we present the first laser source at \SI{497}{nm} based on an extended-cavity diode laser in Littrow configuration. After presenting the characteristics of the used laser diodes in Chapter 2.1 we present the ECDL design in Chapter 2.2, followed by the characterization of its tuning and spectral characteristics in Chapter 2.3 and 2.4, respectively. We conclude with an outlook on future applications in Chapter 3.

\section{Extended-cavity diode laser at \SI{497}{nm}}
\label{sec:setup}

\subsection{Laser diodes}
\label{subsec:diodes}

In the past years GaN diodes covering the spectrum from \SI{370}{nm} to \SI{530}{nm} have become commercially available. For applications using neutral strontium atoms, diodes operating around the \SI{460}{nm} transition reach output powers of up to \SI{100}{mW} (Nichia corporation, \linebreak  NDB4216),  and with amplifier chips multiple thereof should be available in the near future. Diodes based on InAlGaN multiple quantum wells for the wavelength region around \SI{495}{nm} and \SI{500}{nm}, tunable to the $^3P_2  \rightarrow$ $^3D_2 $ and $^3P_2 \rightarrow$ $^3D_3 $ transition in neutral strontium at \SI{496.36}{nm} and \SI{496.93}{nm}, respectively, were recently introduced from Sharp (GH04955A2G) \cite{SharpGH04955A2G}. According to the data sheet the emission wavelength of these non AR-coated diodes spans from \SI{480}{nm} to \SI{495}{nm}. For our purposes we use diodes from a selected batch with a peak wavelength of \SI{495}{nm}.

\subsection{ECDL design}

The design of our diode laser is shown in Figure \ref{fig:design_ECDL}. It consists of an extended-cavity diode laser in Littrow configuration \cite{Ricci1995} and a mirror for angle compensation \cite{Hawthorn2001}.  Light emitted from the laser diode is collimated using an aspheric lens with a focal length $f = \SI{4.02}{mm}$. The collimated light from the diode illuminates a holographic reflective grating (Thorlabs, GH13-24V) and the $-1^{\text{st}}$ order is being reflected back into the laser diode providing the optical feedback. The reflected $0^{\text{th}}$ order is used as the output. The grating period $d$ of \SI{2400}{l/mm} results in a diffraction angle of \SI{36.5}{\degree}. The diffraction efficiency in the $-1^{\text{st}}$ order is about 38 \% and about 44 \% are reflected into the $0^{\text{th}}$ order.

\begin{figure}[htb]
\centering
\includegraphics[width=0.95\linewidth]{simple_ECDL.png}
\caption{Top-view schematic of the extended-cavity diode laser.
The Peltier element underneath the base plate (outer frame in the picture) is not shown.}
\label{fig:design_ECDL}
\end{figure}

The laser diode is mounted in a press fit aluminum mount inspired by the laser pointer community. This ensures good thermal contact between the diode and the rest of the assembly. The grating is embedded into a mirror mount (Radiant Dyes, MNI-H-U-2-3000) to reduce the cavity length and thus to increase the stability. Both the mirror mount and the diode mount are integrated onto a small aluminium base plate. The temperature is stabilized using a digital TEC controller (Meerstetter Engineering, TEC-1091) by means of a Peltier element and an NTC thermistor that is positioned in the base plate between the diode and the grating, inspired by \cite{Saliba2009b}. The diode mount and the mirror mount are fixed with nylon screws and finally glued to the base plate. A mirror is mounted on a small lever parallel to the grating to compensate for the angle change of the beam due to wavelength tuning with the grating angle \cite{Hawthorn2001}. Our diode laser current source is based on the Libbrecht-Hall design \cite{Libbrecht1993}. \SI{86}{\percent} of laser emission passes an optical isolator (LINOS, FI-488-3SC) and is coupled into a polarization-maintaining optical fiber with an efficiency of $\approx$\,\SI{55}{\percent} for optical spectrum analysis.

\subsection{Power and tuning characteristics}

We measured the output power of the bare laser diode and of the assembled ECDL both at a temperature of \SI{25}{\celsius}. The power was measured before the optical isolator. The results are shown in Figure \ref{fig:PI}. According to the data sheet, the absolute maximum injecting current is \SI{135}{mA} with a typical output power of \SI{55}{mW} at \SI{105}{mA}. In our measurement, we operated the laser diode at injection currents up to \SI{200}{mA}. As can be seen from Figure \ref{fig:PI}, the output power of the free running diode reaches \SI{127}{mW} with only a slight deviation from a linear current/power relationship of \SI{5.9}{\percent}. Therefore, an output power of almost 2.5 times as high as the maximum stated rating is achieved. No thermal roll-over could be observed for even higher injection currents up to \SI{500}{mA}, which we address to a good thermal contact of the laser diode to the heatsink, resulting in an enhanced temperature stabilization of our compact setup. At \SI{500}{mA}, we measured an optical power of \SI{242}{mW}.

\label{subsec:pituning}
\begin{figure}[htb]
\centering
\includegraphics[width=0.95\linewidth]{ECDL_PI_paper.png}
\caption{The output power of the bare laser diode (blue) and of the ECDL with optical feedback from the grating (orange) as a function of the injection current at \SI{25}{\celsius}. The inset shows the reduced laser threshold of the ECDL. }
\label{fig:PI}
\end{figure}

In the ECDL configuration the output power reaches \SI{65}{mW} at an injection current of \SI{200}{mA} and the lasing threshold is reduced from \SI{44}{mA} to \SI{36}{mA} compared to the free running configuration. We are not concerned about the back reflected power from the grating since it is not a major source of diode damage for GaN diodes, as in the case for GaAs diodes \cite{Tomm2017}.
The high output power of this ECDL is useful for a wide range of applications such as laser cooling and spectroscopy at cyan wavelengths around \SI{495}{nm}, where only SHG sources could deliver comparable power levels so far. The power could be further increased by injection locking a second diode in a master slave configuration. To further investigate the suitability of our laser source for atomic physics, wavelength tuning range and the laser linewidth are investigated.

The wavelength of the ECDL is tuned by changing the grating angle. Figure \ref{fig:tuning} shows the emission spectra of our ECDL for different wavelengths along with the emission spectrum of the bare diode (grey shadowed), measured with an optical spectrum analyzer (OSA). All measurements were performed at an injection current of \SI{105}{mA} and at a diode temperature of \SI{25}{\celsius}. A tuning range of \mbox{$\approx$\,\SI{8}{nm}} around the center wavelength, covering the range from \SI{491}{nm} to \SI{499}{nm} is achieved. This is comparable to commercial diode laser systems in the visible part of the spectrum and, to our knowledge, the first tunable diode laser covering this part of the spectrum.

\begin{figure}[htb]
\centering
\includegraphics[width=0.95\linewidth]{ECDL_tuning_paper.png}
\caption{Optical spectra of the ECDL for different wavelengths accessible with our configuration (colored) and for laser diode without feedback (grey) measured at \SI{25}{\celsius}. The red data points show the lasing threshold for each ECDL spectrum.}
\label{fig:tuning}
\end{figure}

The resolution of the measurement is limited by our OSA to \SI{0.1}{nm}. The lasing threshold for each of the spectra is also shown in Figure \ref{fig:tuning}. It is strongly reduced compared to the bare diode and is the smallest for the wavelength of \SI{495}{nm}, the peak wavelength of the bare diode.

\subsection{Emission linewidth}
\label{subsec:Linewidth}

The power output of the ECDL is sufficient for laser spectroscopy and cooling applications with neutral strontium as has been shown above. Another criteria for the suitability of a laser in atomic physics is the linewidth. It should be narrower than the linewidth of the transition we want to address, in our case \SI{2.04}{MHz} for the $^3P_2 \rightarrow$  $^3D_2$ transition at \SI{496.93}{nm}.

\begin{figure}[htb]
\centering
\includegraphics[width=0.95\linewidth]{Beatnote_paper.png}
\caption{Analysis of the beat note between two lasers for a sweep time of \SI{400}{ms} averaged for 35 sweeps. For a single laser this results in a Lorentzian linewidth of \SI{88}{kHz} and a Gaussian linewidth of \SI{1.22}{MHz}. See text for more details.}
\label{fig:beat}
\end{figure}

To investigate the linewidth of the diode laser we assembled a second laser and measured the free running beat note between both of them. \mbox{Figure \ref{fig:beat}} shows the analysis of our beatnote measurement where 35 sweeps were averaged for better statistics. Fitting a Voigt model to the data yields a Lorentzian linewidth of \SI{175 \pm 1}{kHz} and a Gaussian linewidth of \SI{1.72 \pm 0.01}{MHz} at a sweep time of \SI{400}{ms}. The Lorentzian linewidth is divided by 2 and the Gaussian linewidth is divided by $\sqrt{2}$ to obtain the respective linewidths
of a single laser. The resulting Lorentzian linewidth of \SI{88}{kHz} is comparable to commercial diode laser systems in the visible part of the spectrum. The Gaussian linewidth of \SI{1.22}{MHz} stems from external influences and from possible noise in our current driver and temperature controller. The linewidth could therefore be further reduced by low noise electronics and an improved laser housing that is less susceptible to mechanical and acoustic vibrations. Even without further improvement, the linewidth is already lower than the linewidths of the transitions in neutral strontium around \SI{497}{nm}.

\section{Conclusion}
\label{sec:conclusion}
We built and characterized an extended-cavity diode laser operating in the cyan wavelength range around \SI{497}{nm}. The laser provides an optical output power of \SI{65}{mW} at \SI{200}{mA} with only a slight deviation from a linear current/power relationship. A tuning range of $\approx$\,\SI{8}{nm} around the centre wavelength of \SI{495}{nm} was achieved. The linewidth of the laser was determined to be \SI{88}{kHz} for the Lorentzian part and \SI{1.22}{MHz} for the Gaussian part.\\
This compact, simple and low cost laser source can be used for several applications. Laser systems for cooling of strontium atoms can be considerably simplified in two different ways. First, the laser source can be utilized as a repump laser at \SI{496.93}{nm} to replace previous more complex systems and to reduce the number of needed laser sources. Second, sub-Doppler cooling at \SI{496.36}{nm} can be employed to replace the second cooling stage at \SI{689}{nm}. Furthermore these laser source facilitates quantum computing schemes using the $^3P_2 \rightarrow$ $^3D_3$ transition at \SI{496.36}{nm} for population detection \cite{Daley2011}, can be used for experimental investigation of the highly forbidden $5s\rightarrow6s$  transition in rubidium and potentially for compact optical frequency standards taking advantage of cyan transitions in molecular iodine \cite{Cheng2002}. Lastly, experiments employing cold barium ions can profit from convenient diode lasers for the cooling transition \mbox{$6s^2S_{1/2} \rightarrow 6p^2P_{1/2}$} at \SI{493.5}{nm} \cite{Hucul2017}.

\begin{acknowledgements}

This work is supported by the DLR Space Administration with funds provided by the Federal Ministry for Economic Affairs and Energy (BMWi) under grant numbers DLR 50WM1753 and 1857.

The Authors would like to thank Achim Peters for useful discussions and providing resources, Klaus D{\"{o}}ringshoff for useful discussions and proofreading the manuscript.
\end{acknowledgements}

\bibliographystyle{ieeetr}