\section{Introduction}
\label{sec:intro}
\change{Laser spectroscopy is a powerful tool to gain insight into the structure of atoms and molecules. Spectroscopy of optical transitions in highly charged ions (HCIs) has been proposed for tests of fundamental physics and for superior optical clocks. For example, due to high relativistic contributions to the binding energy of their electrons, some transitions in HCIs are exceptionally sensitive to a change of the fine structure constant \cite{berengut2010enhanced}. Also, small polarizabilities and high internal fields render HCIs insensitive to external perturbations, promising small systematic shifts in atomic clocks based on optical transitions in HCIs \cite{derevianko_highly_2012,dzuba_optical_2015}. However, laser spectroscopy of HCIs is still in its infancy. The highest resolution in laser spectroscopy of HCIs has been demonstrated on the dipole forbidden $^{2}P_{1/2}$ to $^{2}P_{3/2}$ transition in Ar$^{13+}$ with a pulsed dye laser, reaching a resolution of $\sim$\,400\,MHz of this line with an excited state lifetime of 10\,ms \cite{mackel2011laser}. The measurement was limited by the large kinetic energy of the HCIs in the electron beam ion trap. Sympathetic cooling to the mK regime by embedding the Ar$^{13+}$ ions in a laser cooled crystal of Be$^{+}$ ions has been demonstrated in the Cryogenic Paul Trap Experiment (CryPTEx) \cite{schmoger2015coulomb}. The seven orders of magnitude reduction in the temperature of the ions enables Doppler-limited spectroscopy with a resolution on the order of 100\,kHz. Laser spectroscopy of this transition and similar transitions in other HCIs requires a continuous-wave spectroscopy laser that can be adjusted to operate over a large wavelength range and exhibits a linewidth and long-term drift below 100\,kHz as described in this work.}

\change{There are several ways to achieve this goal.} The spectroscopy laser can be locked to a passive cavity \cite{young_visible_1999,kessler_sub-40-mhz-linewidth_2012,hafner_8_2015}, to a stabilized optical frequency comb \cite{benkler_robust_2013,inaba_spectroscopy_2013,rohde_phase-predictable_2014,nicolodi2014spectral,scharnhorst2015high}, or to a broad atomic or molecular transition either directly \cite{wieman_doppler-free_1976,corwin_frequency-stabilized_1998,ye_ultrasensitive_1998,shirley1982modulation,bjorklund_frequency_1983,mccarron2008modulation} or via a transfer cavity \cite{riedle1994stabilization}. In our specific case, the laser will be transported from PTB in Braunschweig, where it is built and calibrated, to the MPIK in Heidelberg, where the HCI spectroscopy will be performed. While stabilization to a passive cavity can in principle provide the required performance goals, a significant engineering effort is required to make it transportable \cite{vogt2011demonstration,chen2014compact,leibrandt2011field,leibrandt2011spherical} \change{ and to eliminate unexpected shifts and non-predictable long-term drifts of its resonance frequencies \cite{dube2009narrow}}. For decades, lasers have been stabilized to atomic references that can in principle provide light with very small long-term drift.
However, typically two issues have to be overcome: often there are no suitable atomic transitions close to the desired spectroscopy laser frequency, and the linewidth of the transition may be broader than requir ed for the target laser linewidth.
Both issues can be resolved using a frequency comb referenced to a stable cavity or atomic reference that allows to transfer the stability of the reference to the spectroscopy laser on all relevant time scales \cite{nicolodi2014spectral,scharnhorst2015high}. However, a comb and a suitable reference are required for such an approach, both of which are resource intensive and therefore not generally available in all laboratories.
\change{For bridging frequency ranges of up to a few THz, an EOM inside a Fabry-Perot cavity can be used to produce a narrow band frequency comb \cite{matthey2015compact}.} As an alternative, a cavity can be used to transfer the stability of a reference laser to the spectroscopy laser \cite{riedle1994stabilization,bohlouli-zanjani_optical_2006,uetake2009frequency,rohde2010diode,albrecht_laser_2012,yin_narrow-linewidth_2015,Bridge:16}. To achieve simultaneous resonance of both lasers with the transfer cavity, at least one of them must be tunable by one free spectral range (FSR). For common cavity lengths the FSR is on the order of 1\,GHz.

Here, we report on a setup to stabilize a spectroscopy laser with a wavelength of 882\,nm via a transfer cavity to a reference laser operating at the Rb D$_2$-line wavelength of 780\,nm. The reference laser is stabilized to the F=3 to F'=4 transition of the $^{85}$Rb D$_2$-line using the modulation transfer spectroscopy technique \cite{shirley1982modulation,mccarron2008modulation}. Frequency drifts of the atomic reference are suppressed by an active temperature stabilization of the Rb vapor cell and magnetic shielding. We use an evacuated cavity to transfer the \change{stability of the reference laser} to the spectroscopy laser, thus avoiding environmental influences that affect the \change{performance of the transfer lock}. The length of the cavity is actively stabilized to the reference laser frequency using a \change{Pound-Drever-Hall (PDH) \cite{black_introduction_2001} sideband locking technique to provide gap-less tuning.} The spectroscopy laser is stabilized to the transfer cavity in a standard PDH configuration.

\change{Using an optical frequency comb that is locked to a long-term-stable H-maser,} we show that the resulting frequency stability of the reference and spectroscopy laser is below 50\,kHz on a time scale of several days. The short term linewidth is measured via a comparison with an ultra-stable laser. Our setup enables us to narrow the linewidth of the spectroscopy laser to well below that of the reference laser, which reduces the requirements for the latter. \change{We operate the laser system at 882\,nm to produce the spectroscopy wavelength of 441\,nm for Ar$^{13+}$ after frequency doubling. However, the Ti:Sapphire laser can be broadly tuned to access optical transitions in other HCIs. The current set of transfer cavity mirrors supports a wide range of wavelengths spanning from $775\dots 890$\,nm. This range can be further extended to other spectroscopy laser wavelengths by adapting the mirror coating.}

\section{Experimental setup and results} \label{sec:setup}

The experimental setup is schematically shown in Fig.\,\ref{setup}. It is divided into three parts: the reference laser, the cavity transfer lock, and the spectroscopy laser, each of which will be discussed in detail in the following.

\begin{figure}[htb]
\centering\includegraphics[width=0.45\textwidth]{Zeichnung_TiSa_bereinigt}
\caption{Simplified experimental setup showing the reference laser, transfer lock setup and spectroscopy laser. For details about the individual parts of the setup refer to the corresponding sections. PBS: polarizing beam splitter, EOM: electro-optic modulator, PD: photo diode, DM: dichroic mirror, PZT: piezo actuator, PI: proportional-integral controller, VCO: voltage controlled oscillator, PLL: phase-locked loop, HV amp: high-voltage amplifier, PM: phase modulation, $\lambda/2$: half-wave plate, $\lambda/4$: quarter-wave plate, V$_{\textnormal{tune}}$: mixer output/feedback voltage\label{setup}}
\end{figure}

\subsection{Reference laser} \label{sec:ref}

The reference laser consists of a distributed feedback laser diode (Eagleyard EYP-DFB-0780-00080) mounted in a temperature-stabilized housing (Thorlabs TCLDM9). \change{The low noise laser diode driver (Thorlabs LDC201CU) allows modulating the diode pump current with a bandwidth of 100~kHz for frequency stabilisation.} The free running linewidth of the laser diode is specified to be around 2\,MHz. The output power of $\sim$\,50\,mW is divided into two parts: 5\,mW are used to stabilize the frequency to an atomic resonance by modulation transfer spectroscopy (MTS) \cite{shirley1982modulation,mccarron2008modulation}, the remainder is sent \change{through a fiberized electro-optic modulator to provide the modulation for locking the length of} the transfer cavity.

The optical setup for the frequency stabilization of the reference laser is shown in Fig.\,\ref{setup} in the upper left box. A linearly polarized pump beam is phase modulated at $\subsc{\omega}{MTS}\simeq 9.4$\,MHz. It is aligned to counterpropagate along an unmodulated probe beam in a temperature-stabilized rubidium vapor cell.
\change{The laser is subject to dispersion in the rubidium vapor which, around resonance, converts the phase modulation of the pump beam to amplitude modulation. Thus, the saturation of the medium becomes} modulated with $\subsc{\omega}{MTS}$. The counterpropagating probe beam experiences a modulated absorption, which is subsequently detected as amplitude modulation on a photo diode. Homodyne mixing of the photo diode signal with $\subsc{\omega}{MTS}$ yields a dispersion-shaped output signal around the atomic resonance. \change{This MTS signal is tailored by a PI controller and fed back as a modulation of the current via the laser diode driver with a typically achieved feedback bandwidth of 7~kHz.}\par

\begin{figure}[htb]
\centering\includegraphics[width=0.45\textwidth]{MTS}
\caption{Saturated absorption (black) and modulation transfer spectroscopy signal (red). Two Doppler broadened transitions contribute to the D$_2$-line for each isotope. Each consists of three hyperfine transitions. Here the transitions for $^{85}$Rb F=3 to F'=2,3,4 and $^{87}$Rb F=2 to F'=1,2,3 are shown. Only the cycling transitions F=3 to F'=4 ($^{85}$Rb) and F=2 to F'=3 ($^{87}$Rb) produce a strong dispersive signal in the modulation transfer spectroscopy. The inset shows \change{a higher resolution scan across} the $^{85}$Rb transition used to stabilize the reference laser. \label{MTS}}
\end{figure}
In Fig.\,\ref{MTS} the saturated absorption and MTS signals are shown when scanning the laser frequency \change{across} the D2-line\change{s} originating from the F=3 hyperfine ground state in $^{85}$Rb and the F=2 ground state in $^{87}$Rb \cite{steckrubidium85,steckrubidium87}. The advantages of the MTS technique for absolute frequency stabilization can be seen in the figure. Only the closed-cycle transitions F=3 to F'=4 in $^{85}$Rb and F=2 to F'=3 in $^{87}$Rb create a \change{significant} error signal, resulting in a simple identification and reproducibility of the lock point. Furthermore, \change{in contrast to frequency modulation or polarization spectroscopy \cite{mccarron2008modulation} MTS is free of a Doppler background}.
It has been shown that an absolute frequency accuracy below 1\,kHz can in principle be achieved using MTS \cite{MartinezdeEscobar:15}.\par
The vapor cell is temperature stabilized to $\sim$\,10\,mK root mean square (RMS) deviation from the set-point. It is magnetically shielded with two insulated layers of mu-metal. A housing around the cell mount with small openings for the laser beam reduces air flow.\par
The reference laser frequency is locked to the F=3 to F'=4 hyperfine transition in $^{85}$Rb which exhibits the error signal with the steepest slope in our setup.

\subsection*{Performance of the reference laser} \label{sec:perf_ref}
\change{The reference laser performance and systematic shifts are studied by comparing the laser frequency to a hydrogen maser. For this, part of the laser power is overlapped with the output of a maser-stabilized optical frequency comb and a beat note is observed on a high bandwidth photo diode and recorded by frequency counters. We use two independent counters to validate the frequency data.}

Since the rubidium vapor is enclosed in a fixed volume, any temperature change translates directly into a pressure change. \change{This} not only broadens the transition lines, but also shifts them in energy \cite{PhysRev.124.800}. By changing the set-point of the temperature stabilization while observing the beat frequency over time on an electronic counter, a temperature sensitivity of $\sim$\,100\,kHz/K is measured for the transition \change{used here. The temperature fluctuations of the stabilized cell are $\sim$\,10\,mK,} leading to a frequency instability on the few-kHz level.

The differential magnetic field shift of the hyperfine levels used for laser stabilization is 2.3\,MHz/mT \cite{steckrubidium85}. With typical magnetic field fluctuations on the order of $10^{-2}$\,mT and shielding as described above, magnetic field shifts are negligible on the kHz level.\par
\change{The error signal of MTS consists of a signal with odd symmetry with respect to the resonance frequency (originating from phase modulation) and a signal with even symmetry (originating from amplitude modulation). Residual amplitude modulation (RAM) at the EOM resonance frequency thus introduces an asymmetry in the error signal \cite{jaatinen1998possible}. The effect of varying laser power on the error signal is a linear scaling of the voltage response. For a symmetric error signal this only affects the slope. For an asymmetric signal, however, the position of the zero-crossing on the frequency axis depends on laser power. RAM-induced asymmetry thus introduces a laser power dependence of the lock point.
}

In order to study the influence of laser power fluctuations on the laser frequency we vary the optical power going into the MTS setup while recording both laser power and beat frequency. For relative changes below 10\,\% the laser frequency responds to a good approximation linearly to a change in laser power (see Fig.~\ref{pow_sens}). \change{Choosing the laser power level for the MTS is a tradeoff between a larger error signal for high power and a narrower linewidth for lower power. The optimum power level is defined by a maximum in the error signal slope. After determining this} optimum power level, a sensitivity for absolute laser frequency shift per unit optical power shift \change{was} extracted. \change{With the given laser and rf power} we obtain a sensitivity of 91\,kHz/mW. We observe laser power fluctuations up to $\pm 8\%$ or 0.35 mW.  The resulting fluctuations in the zero-crossing are the main source of residual frequency fluctuations in the reference laser. \change{For more demanding applications, this effect could be further suppressed by active laser power stabilization.}
\begin{figure}[htb]
\centering\includegraphics[width=0.45\textwidth]{power_sensitivity}
\caption{Frequency shift of the reference laser relative to optimal laser power recorded versus laser power in the modulation transfer spectroscopy setup. \change{The error bar for each data point is well within the symbol representing it.} For small laser power deviations a linear dependence of 91 kHz per mW is extracted (dashed line). \label{pow_sens}}
\end{figure}

\change{The linewidth of the reference laser was estimated through comparison with a stable optical reference. For this, we stabilized the optical frequency comb to the optical reference \cite{scharnhorst2015high} and recorded the spectral distribution of the beat note between the reference laser and the frequency comb with a spectrum analyzer.} This signal is a convolution of the $\sim$\,100\,kHz broad comb tooth with the reference laser line. \change{From a Gaussian fit to the signal we derive an upper bound for the linewidth of the reference laser of 400\,kHz.}

\begin{figure}[htb]
\centering\includegraphics[width=0.45\textwidth]{drift}
\caption{Reference laser frequency measurements relative to the absolute frequency $\nu_0=384.229\,240\,061$\,THz of the first data point. Each data point is averaged over 30\,seconds. \change{The error bars for each measurement are well within the size of the symbol representing it.} The offset of 1.6\,MHz from the literature value \cite{steckrubidium85} arises from a finite magnetic field, a pressure shift and the asymmetric error signal. \label{rb_drift}}
\end{figure}

The long-term frequency stability is evaluated by several comparisons against an active hydrogen maser frequency standard with the help of the frequency comb over the course of 4 days. Fig.\,\ref{rb_drift} shows the frequency deviation of the reference laser for each measurement relative to the first measurement. In a time interval of 4~days the change in frequency is well within a 100\,kHz band, the RMS deviation being 20\,kHz. At a wavelength of 780\,nm this corresponds to a fractional instability of the absolute laser frequency of $5\times 10^{-11}$.

\subsection{Transfer lock} \label{sec:cav}
A schematic overview of the cavity transfer lock is shown in the upper right-hand box in Fig.~\ref{setup}. To sustain a transfer lock bridging 100\,nm in the optical/infrared range, care has to be taken to suppress differential dispersion inside the cavity arising from temperature, pressure and relative humidity fluctuations in the air. To that end, evacuated cavities are commonly employed \cite{hough1984dye,helmcke1987new,zhao1998computer}.

\change{We employ a cavity spacer made from a 25\,mm diameter Invar rod with a through hole. Together with a ring piezo attached to one end of the rod, the total cavity length amounts to 100\,mm.} Two curved mirrors with $R=250$\,mm are glued to both ends of the spacer, yielding a symmetric cavity with a free spectral range of about 1500\,MHz. The length of the cavity was chosen such that higher order modes do not coincide with PDH modulation signals. \change{This approach avoids distortion of the main signal from beam pointing fluctuations that change the optical power in these modes \cite{amairi2013reducing}.}

The broadband high-reflective coating of the cavity mirrors reflects more than 99\,\% of the light in the wavelength range between $775\dots 890$\,nm. For the employed wavelengths of 780 and 882\,nm, we measure a finesse of about 300 and 400, respectively. The spacer is placed in a small vacuum chamber with vibration insulated mounting. An attached 2\,l/s ion pump maintains a vacuum of $\sim$\,$10^{-8}$\,mbar.

The electronic sideband locking technique \cite{thorpe2008laser} based on PDH stabilization is used to stabilize the length of the cavity to the 780\,nm reference laser. The basic idea is to create a sideband with a tunable offset to the carrier laser frequency. This sideband itself exhibits phase modulation sidebands for PDH stabilization. If the variable offset spans a range of half an FSR, the cavity can be locked to one or the other of the offset sidebands. \change{This allows the cavity to be locked to the reference laser for any desired spectroscopy laser frequency.}
The electronic signal necessary to create the above mentioned sidebands is a phase modulated rf signal with variable carrier frequency up to FSR/2, in our case 750\,MHz.
To create the phase modulated rf signal (see upper right box in Fig.\,\ref{setup}), a VCO (Mini-Circuits ZX95-800C) at $798\dots 803$\,MHz is locked to a fixed-frequency PLL at 800\,MHz, referenced to a maser-stabilized 10\,MHz source. \change{The control input of the VCO used to stabilize the carrier frequency is modulated via a bias tee with $\subsc{\omega}{PDH} \simeq 9.3$\,MHz, thus phase modulating the VCO output at the PDH frequency.} Mixing the VCO output with a signal generator at $810\dots 1600$\,MHz creates sum and difference frequencies. A low-pass filter with a cutoff at 800\,MHz \change{is used} to remove the sum frequency. \change{The difference frequency $\subsc{\omega}{var}$ can be varied between $10\dots 800$\,MHz, while $\subsc{\omega}{PDH}$ remains fixed to $\sim$\,9.3\,MHz.} This signal is applied to a fiber-coupled EOM (Photline NIR-MPX800-LN-05) with an rf bandwidth of 5\,GHz.\par
Fig.\,\ref{SBL} shows a trace of the error signal obtained in our setup together with the corresponding cavity transmission \change{recorded by a photo diode}. \change{Both sidebands at $\pm \subsc{\omega}{var}$ from the carrier can be used for locking the cavity, thus covering one free spectral range of 1500\,MHz.}

\begin{figure}[htb]
\centering\includegraphics[width=0.45\textwidth]{sideband_lock}
\caption{Photo diode signal of the transfer cavity transmission around the transition frequency $\nu_0$ for F=3 to F'=4 in $^{85}$Rb  (red) and corresponding PDH error signal (black).
\change{The inset shows a zoom of the sideband at $\subsc{\omega}{var}$, showing the PDH sidebands at $\subsc{\omega}{var}\pm\subsc{\omega}{PDH}$.} The slopes of the error signals at $\pm\subsc{\omega}{var}$ have opposite signs. \label{SBL}}
\end{figure}

The feedback loop stabilizing the cavity length to the reference laser is electronically bandwidth limited to 150\,Hz. Therefore, frequency noise on the reference laser with Fourier frequencies above 150\,Hz does not affect the cavity length. This allows us to transfer the mean frequency stability of the reference laser without being limited by its high-frequency noise.

\subsection{Spectroscopy laser} \label{sec:spec}

A Ti:Sapphire ring laser (TekhnoScan TIS-SF-07) is used as a light source at 882 nm. It is optically pumped with 10\,watts at 532\,nm by a diode pumped solid state laser (Verdi G10). \change{On the order of 20\,mW of the total output power of $\sim$\,500\,mW is used for locking the laser to the transfer cavity. The remaining power is frequency doubled in an external second harmonic generation stage using a temperature-tuned periodically-poled KTP crystal. An output power of 200 mW at 441 nm is routinely achieved.}

The beam towards the transfer cavity is phase modulated at 91\,MHz in a resonant-circuit EOM (TimeBase EOM-320IR). A PDH setup generates an error signal from the cavity-reflected beam. The signal is modified by a PI controller which feeds back to two feedback channels acting upon different piezo actuators inside the \change{laser} resonator. \change{The fast channel drives a fast piezo actuator with low stroke. A small mirror mounted on this actuator compensates} high-frequency fluctuations with low modulation index. The resonance frequency of this feedback channel is above 50\,kHz. \change{The slow feedback channel drives three large piezo ring actuators with mounted cavity mirrors.} This channel compensates low-frequency fluctuations up to several 100\,Hz with high amplitude.

\subsection*{Performance of the spectroscopy laser} \label{sec:perf_spec}
The linewidth of the locked spectroscopy laser is characterized by comparison with an ultra-stable laser at 1542\,nm, locked to a high-finesse passive cavity resulting in a linewidth of $\sim$\,1\,Hz. \change{For both lasers a beat note with the frequency comb and additionally the carrier envelope offset is recorded using a 200 MS/s digital oscilloscope which digitizes $10^7$ data points within 50 ms. A numerical evaluation according to \cite{telle2002kerr} yields the frequency difference over time between our spectroscopy laser and the ultra-stable laser. As the frequency variations of the ultra-stable laser can be neglected on that time scale, fluctuations in the frequency difference can be assigned to absolute frequency fluctuations of the spectroscopy laser.} The inset of Fig.~\ref{linewidth} shows the distribution \change{of the frequency fluctuations of the spectroscopy laser} around the mean value. \change{As the spectral distribution is not Gaussian}, an estimate for the linewidth is inferred from the frequency noise spectral density. In Ref.~\cite{di2010simple} a method is presented to extract a laser linewidth from arbitrary spectral noise distributions.
\begin{figure}[htb]
\centering\includegraphics[width=0.45\textwidth]{linewidth}
\caption{Frequency noise spectral density $S_\nu$ of the spectroscopy laser in transfer locked state (solid line). The $\beta$-separation line (dashed line) separates the areas of high and low modulation index. The inset shows a 50\,ms snapshot of the laser frequency compared to another laser with $\sim$\,1\,Hz linewidth. \label{linewidth}}
\end{figure}
The frequency noise spectrum can be divided into two regions. One with a high noise level compared to the Fourier frequency, thus exerting a modulation index $\beta > 1$, and one with a low noise level compared to the Fourier frequency, yielding a low modulation index. The first region contributes to the linewidth, while the latter mostly affects the wings of the laser line. These regions are separated by the so-called $\beta$-separation line. \change{The full width at half maximum (FWHM) of the laser linewidth is extracted by integrating the frequency noise spectral density, $S_\nu$, up to the Fourier frequency where it intersects the $\beta$-separation line. In our case this is at about 2.2 kHz (see Fig.\,\ref{linewidth}). The lower cut-off frequency is given by the observation time of 50 ms and corresponds to 40 Hz. For the integration we most likely overestimate $S_\nu$ by assuming a constant value between 0-40 Hz.} Applying this method we derive a FWHM of 55\,kHz. The RMS deviation of the laser frequency distribution in the inset of Fig.\,\ref{linewidth} is 23\,kHz. \change{We expect that frequency doubling in an external cavity doubles the linewidth without introducing significant additional noise \cite{liu2007narrow}. Frequency doubling is performed with a PPKTP crystal in a bow-tie external doubling cavity. We routinely achieve an output power of 200\,mW at 441\,nm with 500\,mW power at the fundamental frequency.}

\FloatBarrier
We characterize \change{the laser frequency} stabilization by comparing the power spectral densities (PSD) of the in-loop error signal for two cases: a tight lock with optimized parameters and a loose lock where the stabilization is just sufficient to keep the laser frequency on the central slope of the PDH error signal. That way, we obtain a nearly linear voltage response for frequency fluctuations. Fig.\,\ref{noise} shows PSD spectra up to 100 kHz recorded with a fast Fourier transform device. It can be seen that the feedback loop suppresses noise within a bandwidth of about 18\,kHz.

\begin{figure}[htb]
\centering\includegraphics[width=0.45\textwidth]{tight-vs-loose}
\caption{Phase noise power spectral density of the Ti:Sapphire laser in tight lock (black) and loose lock (red). A noise suppression of up to 50\,dB is established by the cavity lock. We infer a feedback bandwidth of around 18\,kHz for the fast feedback channel from the crossing point of the two curves. The difference in PSD above 18\,kHz between both curves \change{arises from excursions of the laser frequency beyond the linear region of the PDH error signal, leading to high-frequency harmonics of
low-frequency noise components.} \label{noise}}
\end{figure}

With the transfer lock engaged, the frequency instability of the spectroscopy laser is measured the same way as before for the reference laser. The wavelength of the spectroscopy laser is tuned to 871\,nm where an output port of a frequency comb is available. The beat frequency between spectroscopy laser and frequency comb is monitored over several days.
\begin{figure}[htb]
\centering\includegraphics[width=0.45\textwidth]{stability_TiSa}
\caption{(a) Beat frequency of the transfer locked spectroscopy laser with a frequency comb at 871\,nm. (b) Absolute and relative overlapping Allan deviation of the transfer locked spectroscopy laser up to averaging times of $6.5\times 10^{4}$\,s. \label{drift_tisa}}
\end{figure}
Fig.\,\ref{drift_tisa} (a) shows a frequency trace taken over 90~hours. A linear fit to the data suggests a frequency drift of $(0.28 \pm 0.30)$\,kHz/h. The RMS frequency deviation is 44\,kHz. More information on the stability can be extracted from the Allan deviation in Fig.\,\ref{drift_tisa} (b). We achieve a frequency instability below $10^{-10}$ for averaging times up to 18\,h. At that point the Allan deviation flattens out, indicating no linear drift on this time scale. \change{The oscillations in the laser frequency visible in Fig.\,\ref{drift_tisa} (a) (also visible as local maxima in the Allan deviation in Fig.\,\ref{drift_tisa} (b)), are most likely of thermal origin. The two characteristic times of about 4 s and 2000 s can be traced back to the control loop for laser diode temperature and room temperature, respectively. Oscillations in temperature affect the laser power and any spurious etalon fringes between optical elements. This leads to frequency fluctuations on these time scales which are, however, within the 50\,kHz tolerance required for the intended application.}\par
We expect that the lower limit in Allan deviation is given by the instability of the reference laser. The stability of the spectroscopy laser is further deteriorated by offset fluctuations in the PDH error signals due to residual amplitude modulation.

\section{Conclusion}
\change{We have presented a laser setup that is capable of providing a gapless output range spanning over 100\,nm in the infrared spectrum. We have demonstrated that the fast linewidth and long-term frequency drifts are below 50\,kHz over several days.}
These values were verified by measurements with a referenced optical frequency comb. We show that the spectral impurity of the reference laser does not limit the spectroscopy laser linewidth due to low-pass filtering of frequency noise of the transfer cavity. The transfer locking technique demonstrated here can be extended to cover a broad range of wavelengths in the optical and infrared, limited only by the cavity mirror coating. \change{To our knowledge this is the best combined performance in terms of linewidth and long-term stability for a transfer lock using a broad alkali D-line as reference transition. Better results can be achieved on narrower atomic transitions but require an intrinsically narrower linewidth of the reference laser \cite{Bridge:16}.}\par
Current limitations in our frequency stabilization setup arise mainly from the reference laser. A simple improvement would be a laser power stabilization of the reference laser and a temperature stabilization of the EOM used for the modulation transfer spectroscopy to reduce time-varying RAM. Furthermore, \change{better thermal isolation} of the rubidium vapor cell would further improve the stability of the reference. To narrow the short term linewidth of the spectroscopy laser, a higher feedback bandwidth is necessary. An intracavity EOM or external AOM \cite{haubrich1996modified,boyd1991frequency} could be used for this purpose. Additionally, the cavity mirror coating should then be tailored in such a way, that the finesse is substantially higher for the spectroscopy laser while keeping a value on the order of 100 for the reference laser.\par
With this setup we anticipate to resolve the line of the fine structure transition $1s^{2}2s^{2}2p\ ^{2}P_{1/2} -\ ^{2}P_{3/2}$ at 441\,nm in Coulomb-crystallized Ar$^{13+}$ ions in a Paul trap to within 100\,kHz.

\section*{Acknowledgments}
We acknowledge support from DFG through QUEST. I.D.L. acknowledges a
fellowship from the Alexander von Humboldt Foundation. This work was
funded by PTB.


\section{Introduction}\label{sect:intro}
The action of a dynamo in stars with convective envelopes leads to the generation of magnetic field \citep[cf.][]{dynamoLR}. The magnetic field emerges on the stellar surface, forming structures such as dark spots and bright faculae, and then rises  further to fill the outer stellar atmosphere \citep[cf.][]{Sami_B}. The emergent magnetic field can be traced via sophisticated techniques such as Doppler or Zeeman-Doppler imaging \citep[see reviews by][and references therein]{berdyugina2005a, LR_Reiners},  or through the non-thermal UV \citep[][]{Hall_LR} and  X-ray  \citep[][]{Gudeletal2004} emission.

An alternative way of studying stellar magnetic fields is provided by brightness variations caused by surface magnetic features. The transit of bright and dark magnetic features across the visible stellar disc as the star rotates, as well as the evolution of magnetic features  cause  changes of stellar brightness on the timescale of the stellar rotation. The overall modulation of the stellar disc coverage by magnetic features over the activity cycle leads to brightness changes on timescales of the stellar activity cycle.

The brightness variations of the Sun and Sun-like stars, dominated by magnetism except at the shortest and longest timescales, has been studied for the last few decades. The variations of solar brightness were detected with the start of regular spaceborne  Total Solar Irradiance, TSI (total radiative flux from the Sun, normalized to one AU) monitoring \citep[][]{hickeyetal1980, wilson1981}. In particular, it was established that TSI varies on the 27-day rotational and 11-year activity timescales \citep[see reviews by][] {Claus_rev,koppetal2014}. The same has been noted of the spectrally-resolved solar irradiance, the so-termed Spectral Solar Irradiance, SSI \citep{Floydetal2003,harderetal2009, delandandcebula2012}.

Synoptic observations at the Lowell and Fairborn observatories led to a detection of photometric variability in a few dozens of Sun-like stars \citep{lockwoodetal1997, radicketal1998, lockwoodetal2007, halletal2009}. While there have been some attempts to retrieve short-term  (i.e. night-to-night) stellar photometric variability from the Lowell data \citep[see, e.g.,][]{radicketal1998}, the large uncertainties in the individual measurements make these data more suitable for studying the long-term variability on timescales of stellar activity cycles.

The advent of the  CoRoT \citep{COROT2,COROT} and Kepler \citep{KEPLER}  missions initiated a new era in the study of stellar photometric variabilities. CoRoT and Kepler allowed measuring brightness variations on timescales up to about three months with unprecedented precision and cadence.
The uncertainties in the long-term calibration and relatively short period of CoRoT and Kepler measurements makes them more suitable for studying  stellar brightness variations on the timescale of stellar rotation making them complementary to the Lowell and Fairborn data.
There have been a number of studies \citep[e.g.][]{basrietal2010, basrietal2011,basrietal2013, McQuillanetal2012, McQuillanetal2014, Garciaetal2014}, performed with Kepler data,  aimed at understanding whether solar variability is typical or rather weak compared to the majority of Sun-like stars and to assess the fraction of stars which are more variable than the Sun. These studies have important applications for solar physics and climate research as they allow an assessment of the dynamic range of solar variability and consequently of its influence on the Earth's climate.

The important question in studies of stellar photometric variability is which kind of magnetic features dominate the photometric variability. In particular,  whether photometric variability is attributed to large concentrations of magnetic field  which form dark spots or to smaller magnetic elements, the ensemble of which forms bright faculae and network. It is known that {\it cyclic} (i.e. over the activity cycle) photometric variability of the Sun and Sun-like stars less active than the Sun is faculae-dominated, while that of more active stars is spot-dominated \citep{Radicketal1990, lockwoodetal2007,halletal2009,Lockwood2013_stars}.
This observation, together with the general patterns of variation (e.g. dependence of the brightness variability on mean chromospheric activity) established by the Lowell and Fairborn programs
for the long-term {\it cyclic} photometric  variability of Sun-like stars  is often extrapolated to shorter timescales of stellar rotations \citep[cf.][]{Gilliland2011,Harrison2012,Lewis2013,Dumusque1,Dumusque2}. At the same time the
applicability of such an extrapolation  to  the rotational timescale has not been directly addressed in the literature so far. One of the aims of the present paper is to determine if this assumption is correct.

There have been several studies aimed at connecting photometric variability of Sun-like stars with magnetic features on their surfaces. For example, \cite{Lanzaetal2003,Lanzaetal2006} developed a model for simulating brightness variations of the Sun and Sun-like stars. They assumed that stellar magnetic features consist of a mixture of spot and facular areas in a fixed proportion. \cite{Lanzaetal2009a,Lanzaetal2009b}  applied this model to light-curves of  CoRoT-4a  and CoRoT-Exo-2a measured by CoRoT, while \cite{Gondoin2008} applied it to light-curves of $\epsilon$ Eri and $\kappa$ Ceti measured by the MOST microsatellite. Their analysis suggested that the facular to spot area ratio is somewhat lower for more active stars \citep[cf.][]{foukal1998,solankiandunruh2012,Shapiro2014_stars}.

Facular signatures have been found in the frequency spectrum of the green photometry recorded by the Sun PhotoMeter on the SoHO/VIRGO experiment \citep{Karoff2012} as well as in the frequency spectra of the high-cadence Keplar light curves of several Sun-like stars \citep{Karoffetal2013}. Efforts have also been made to estimate the effect of faculae on stellar radial velocity signatures \citep[see][and references therein]{Meunieretal2010,Jeffersetal2014,Dumusque1,Dumusque2}. However, the effect of faculae on the rotational stellar brightness variability measured by Kepler and CoRoT has not been addressed so far.

In this paper we employ the SATIRE \citep[Spectral And Total Irradiance Reconstruction, see][]{fliggeetal2000, krivovaetal2003} model to calculate the relative contributions of faculae and spots to magnetically-driven solar brightness variability. We study how these contributions depend on the wavelength, timescale, and position of the observer relative to the ecliptic plane. In particular, we investigate whether facular signatures can be seen in the Kepler  light curves of a star identical to the Sun.

In Sect.~\ref{sect:model} we briefly describe the SATIRE approach.
In Sect.~\ref{sect:Earth}, we calculate the contribution by faculae and spots to solar brightness variability, {\it as apparent to an Earth-bound observer}, and the wavelength- and timescale-dependence. In Sect.~\ref{sect:out} we investigate how the relative role of facular and spot contributions varies with the position of the observer relative to the ecliptic plane. We also examine if the change in apparent solar brightness variability with observer perspective is similar at rotational and cyclical timescales, as is often assumed in the literature.
In Sect.~\ref{sect:K-C-out} we calculate the contributions of faculae and spots to the Kepler light curves of a star identical to the Sun.  Finally, we summarize our main results in Sect.~\ref{sect:conc}.

\section{The model}\label{sect:model}
\subsection{SATIRE-S}\label{subsect:SAT}
The calculations presented in this paper are based on the SATIRE model. SATIRE attributes the solar brightness variability on timescales of a day or longer to the time-dependent contributions from dark  and bright surface magnetic features.
It models the visible solar disc as comprising of quiet Sun, and three classes of magnetic features: sunspot umbrae, sunspot penumbrae, and faculae (which encompasses all non-sunspot magnetic activity, including network). The spectra of the quiet Sun and magnetic features  at different disc positions have been calculated by \cite{sat_spectra} with the ATLAS9 code \citep{kurucz1992,ATLAS9_CK}.  Depending on the timescale, different approaches are used to derive the solar disc coverage by magnetic features.
\begin{figure*}
\resizebox{\hsize}{!}{\includegraphics{Fig_TSI.eps}}
\caption{The TSI from the SATIRE-S reconstruction (red curves) and from the PMOD composite (black) for the entire period of observations (panel a), and for the annual intervals representing high (panels b and c) and low (panel d) levels of solar activity. Three black rectangles in panel a constrain the periods and TSI ranges shown in the other panels. {  The presented SATIRE-S reconstruction (panel a) is based on full-disc observations of intensity and magnetic flux from the Kitt Peak Vacuum Telescope, the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory, and the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory. The parts of the reconstruction displayed in panels b, c and d are based on the MDI data.}}
\label{fig:KL}
\end{figure*}

There are several implementations of SATIRE, differing by the solar observations used to infer faculae and sunspot coverage \citep{SATIRE}. Here, we utilize the version based on full-disc magnetograms and intensity images, SATIRE-S, the suffix denoting "satellite-era"  \citep{balletal2014, yeoetal2014}.
This makes SATIRE-S the only model returning TSI and SSI simultaneously and based on direct measurements of solar surface magnetic field rather than on indirect proxies of solar magnetic activity \citep[see also discussion in][]{SATIRE}.

In SATIRE-S, the Spectral Solar Irradiance $S(t,\lambda)$ is calculated by summing up the contributions from image pixels over the solar disc, i.e.
\begin{equation}
S(t,\lambda)= S^{\rm QS}(\lambda)  + \sum\limits_{mn} \sum\limits_k \left ( I_{mn}^k(\lambda) - I_{mn}^{\rm QS}(\lambda)  \right ) \, \alpha_{mn}^k(t) \, \Delta \Omega.
\label{eq:sum}
\end{equation}
The first summation on the right side of Eq.~(\ref{eq:sum}) is done over the magnetogram pixels (with $m$ and $n$ being the  abscissa and ordinate of the pixel, respectively). The second summation is done over the magnetic components of the solar atmosphere (i.e. umbrae, penumbrae, and faculae).  $I_{mn}^k(\lambda)$ and $I_{mn}^{\rm QS}$ are the intensities at wavelength $\lambda$ emergent at the position of the pixel \{m,n\} from the magnetic component $k$ and the quiet Sun, respectively. $\alpha_{mn}^k(t)$ is the coverage of the  pixel \{m,n\} by the component $k$, and $\Delta \Omega$ is the solid angle of the solar disc segment corresponding to one magnetogram pixel when magnetogram is obtained at 1 au from the Sun. $ S^{\rm QS}(\lambda) $ is the solar irradiance as it would be observed if the visible part of the solar disc were free from spots and faculae (i.e. completely covered by the quiet Sun). It is given by
\begin{equation}
S^{\rm QS}(\lambda)= \sum\limits_{mn}  I_{mn}^{\rm QS}(\lambda)  \, \Delta \Omega.
\label{eq:total}
\end{equation}

SATIRE-S has been demonstrated to successfully reproduce the apparent variability in TSI and SSI measurements from various monitoring missions
 \citep[see][and references therein]{balletal2012,balletal2014, yeoetal2014, Sat_Mn}. Figure~\ref{fig:KL} compares daily values returned by SATIRE-S with the observed TSI \citep[according to the PMOD composite;][] {pmod_comp} on the solar cycle timescale (Fig.~\ref{fig:KL}a) and on the rotational timescale (Fig.~\ref{fig:KL}b--d). Conspicuous dips in the observations  (e.g. in September 2000 and in August 2002) are associated with transits of large sunspot groups. The agreement between SATIRE-S and observations appears to worsen somewhat in 2008 but this is associated with the increased noise contribution to the observations (note the difference in vertical scales between panels b or c and d of Fig.~\ref{fig:KL}). Overall, SATIRE-S replicates 92\% of the TSI variability in PMOD composite and 96\% of the TSI variability when only the most reliable 1996--2013 magnetograms and TSI measurements are considered \citep[both numbers are according to][]{yeoetal2014}.

\subsection{The solar brightness variability observed from the ecliptic}\label{subsect:ecl}
Using daily observations from the Michelson Doppler Imager onboard the Solar and Heliospheric Observatory \citep[SOHO/MDI;][]{scherreretal1995} and from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory \citep[SDO/HMI;][]{HMI}, we reconstructed SSI over the period of 2 February 1999 and 1 August 2014. Although MDI observations extend back to 1996, we disregard the pre-1999 data
\citep[following][]{balletal2012,yeoetal2014}. \cite{balletal2012} reported evidence that the response of the MDI magnetograph might have changed over the SoHO vacation. In addition to the $S(\lambda, t)$ time series we calculate ``{\it facular}'' $S_{\rm fac}(\lambda, t)$ and ``{\it spot}'' $S_{\rm spot}(\lambda, t)$ time series by omitting the spot term in Eq.~(\ref{eq:sum}) in the first case and the facular term in the second case.
The $\alpha_{mn}^k(t)$ values have been utilised from \cite{yeoetal2014}, who employed the harmonising procedure which allows avoiding any inconsistencies between the SOHO/MDI and SDO/HMI segments.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{fig_im.eps}}
\caption{The umbrae (plotted in black), penumbrae (plotted in blue), and faculae and network (both plotted in red)   as they would be seen on the solar disc on October 31, 2003 about midnight UTC  at four different inclinations: $i=90^{\circ}$ (i.e. from the ecliptic plane, panel a), $i=60^{\circ}$ (panel b), $i=30^{\circ}$ (panel c), and $i=0^{\circ}$ (panel d).}
\label{fig:images}
\end{figure}

We define the amplitude of the solar brightness variability on the activity timescale as
\begin{equation}
\Delta {\rm SSI} (\lambda) = \frac{1}{2} \, \frac{ <{\rm SSI}(\lambda, t)>_{\rm 2000} -  <{\rm SSI}(\lambda, t)>_{\rm 2008}}{  <{\rm SSI}(\lambda, t)>_{\rm 2000}  +  <{\rm SSI}(\lambda, t)>_{\rm 2008} },
\label{eq:11}
\end{equation}
where the annual averaging is performed over the years representing periods of high (year 2000) and low (year 2008) solar activity.

The amplitude of the solar brightness variability on the timescale of solar rotation is defined as
\begin{equation}
\sigma {\rm SSI} (\lambda) = {\rm RMS} \left (     \frac {   {\rm SSI}(\lambda, t) - <{\rm SSI}(\lambda, t)>_{81}}{<{\rm SSI}(\lambda, t)>}    \right ),
\label{eq:rot}
\end{equation}
where  ``RMS'' stands for the root mean square. The angle brackets in the numerator and denominator of  Eq.~(\ref{eq:rot})  denote the 81-day running mean and the average over the entire time series, respectively.

We feed Eqs.~(\ref{eq:11})--(\ref{eq:rot}) with $S(\lambda, t)$, $S_{\rm fac}(\lambda, t)$, and $S_{\rm spot}(\lambda, t)$ time series to calculate the {total}  amplitude of the solar brightness variability, as well as its {\it facular} and {\it spot} components, respectively.

\subsection{The solar brightness variability observed out-of-ecliptic}\label{subsect:out}
While the Earth-bound observer always have a near-equatorial perspective of the Sun (hereafter we will neglect  the angle  between solar equator and ecliptic which is very small, being equal to $\sim 7.25{^{\circ}}$), the same is not true of other stars.
To calculate the solar brightness variability as it would be seen from an out-of-ecliptic view-point we project the apparent coverages of magnetic features obtained by \cite{yeoetal2014}  from SOHO/MDI and SDO/HMI measurements to out-of-ecliptic viewing (see Fig.~\ref{fig:images}).

The projected distributions of magnetic features and corresponding  projected values of pixel coverages  $\alpha_{mn}^{k\,i}(t)$  are used to calculate $S^i(\lambda, t)$, $S^i_{\rm fac}(\lambda, t)$, and $S^i_{\rm spot}(\lambda, t)$ time series, where $i$ is the angle between the solar rotational axis and direction to the observer (hereafter referred to as the inclination). We then substitute $S^i(\lambda, t)$, $S^i_{\rm fac}(\lambda, t)$, and $S^i_{\rm spot}(\lambda, t)$ time series into Eq.~(\ref{eq:11}--\ref{eq:rot}) to recalculate the amplitude of total solar brightness variability and its {\it facular} and {\it spot} components for arbitrary inclination $i$.

The magnetic field measurements, and consequently the solar surface  coverage by faculae and network, are rather uncertain towards the poles due to foreshortening.
At the same time the contribution of polar faculae and network to solar brightness variability is very small independently of the inclination \citep[see detailed discussion in][]{luis2012} so that we do not expect that the uncertainty in their disc coverages here will affect our results.

An out-of-ecliptic perspective reveals magnetic features on the far side of the Sun that would otherwise be obscured from an ecliptic viewpoint. Since we do not know the true longitudinal distribution of magnetic features on the far side of the Sun, in our calculations the distribution of magnetic features on the far-side of the Sun {   is obtained by projecting the distribution of magnetic features on the near-side of the Sun on to the far side assuming that the two sides are point symmetric with respect to each other through the centre of the Sun.  In other words the magnetic features on the far-side of the Sun are assumed to be $180^{\circ}$ in longitude apart and with opposite latitudes from the features observed on the near-side of the Sun.}

We do not expect this assumption to introduce any unphysical effects on the calculated cyclic variability since the distribution of magnetic flux should, on average, not differ between the near and far sides of the Sun. The shortcoming of this algorithm for calculating the brightness variability on the rotational timescale is that it does not allow proper accounting for the magnetic features rotating into and out of view. For example, magnetic features will abruptly appear and disappear on the solar disc even when the Sun is observed from its rotational axis. However, the main contribution to the solar brightness variability on the rotational timescale comes not from the disappearance and  appearance of magnetic features but rather from their evolution and the centre-to-limb variations (CLV) of their contrasts. In Appendix \ref{app:dev} we show that our omission of the proper accounting for the disappearance and  appearance of magnetic features has only minor impact on the calculated amplitude of solar brightness variability on the rotational timescale.

\section{Facular and spot components of solar brightness variability observed from the ecliptic}\label{sect:Earth}
In this section we utilize the model presented in Sect.~\ref{sect:model} to assess the contributions of faculae and spots to solar brightness variability observed from the ecliptic plane on the 11-year activity cycle (Sect.~\ref{sub:11}) and on rotational (Sect.~\ref{sub:rot}) timescales.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{Fig2.eps}}
\caption{The amplitude of the 11-year solar brightness variability (black curves) and its facular (red) and spot (blue) components expressed in mmag.  We use axes with different scales shortward and longward of 400 nm (left and right parts  and Y-axes of the figure, respectively separated by the vertical black line). The green shaded contours indicate the transmission curves of the Str{\"o}mgren filters $b$ and $y$ (centred at 467 and 547 nm, respectively). The scale for the transmission curves is chosen so that transmission values 0\% and 100\%  correspond to 0 and 5 mmag, respectively (according to the right axis).
\label{fig:prof_act}
\end{figure}

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{Fig1abc.eps}}
\caption{The amplitude of the rotational solar brightness variability (black curves) and its facular (red) and spot (blue) components calculated for the 1999-2014 period (panel a), as well as averaged over the years 2000 (panel b) and 2008 (panel c). The green shaded contours indicate the transmission curves of the Str{\"o}mgren filters $b$ and $y$. The dark (light) shaded areas show Kepler (CoRoT) total spectral efficiency.}
\label{fig:prof_rot}
\end{figure}

\subsection{Activity cycle timescale}\label{sub:11}
In Fig.~\ref{fig:prof_act} we plot the amplitudes of the 11-year solar brightness variability and its  facular and spot components (see Sect.~\ref{subsect:ecl}) as functions of wavelength.  One can see that the facular component computed with SATIRE-S overweights that of the spots in the UV, visible, and almost entire near-IR spectral domains.

The amplitude of the facular component strongly decreases with wavelength. This comes in part from the diminishing sensitivity of the Planck function to the temperature contrast beween facular and quiet regions.
Due to higher temperature contrast between spot and quiet Sun, the spot component decreases with wavelength not as fast as the facular one so that the regime of the solar brightness variability on the activity cycle timescale switches  from faculae- to spot-dominated about 1200 nm.

Figure~\ref{fig:prof_act} demonstrates that the spectral dependences of the spot and facular components are remarkably different. The facular component has a very complex profile due to the strong influence of  Fraunhofer lines \citep{sat_spectra, yvonne2008,profiles}. The effect of Fraunhofer lines on the spot component is much more subtle  and exposes itself as barely visible spikes superposed on otherwise undulating dependence \citep[see discussion in][]{yvonne2008,profiles}.

The facular component drops abruptly longward of the CH G-band at 430 nm, triggering a similar drop in overall solar brightness variability. As a result, solar brightness variability almost reaches zero (i.e. facular and spot components nearly compensate each other) in the Str{\"o}mgren $b$  filter and is relatively small in the Str{\"o}mgren  $y$ filter (see Fig.~\ref{fig:prof_act}). These filters are widely used in ground-based synoptic programs aimed at the brightness monitoring of Sun-like stars \citep[see e.g.][]{radicketal1998,lockwoodetal2007,halletal2009,Lockwood2013_stars}.

To compare the solar brightness variations to those of Sun-like stars the variability of the solar Str{\"o}mgren  $(b+y)/2$ flux  is often estimated by scaling the TSI variability, since no direct measurements of solar brightness variability in Str{\"o}mgren $b$  and $y$ filters exist.
This is usually done assuming that the solar brightness  variation is caused by a change of the solar effective temperature \citep[][and references therein]{lockwoodetal2007}.  Such an approach yields 1.1 mmag for the amplitude of the 11-year solar cycle variability, identical to that of what is arguably the best solar twin, 18 Sco \citep[cf.][]{petitetal2008}, 1.1 mmag \citep[see detailed discussion in][]{Lockwood2013_stars}. Our calculations show that after a proper radiative transfer computation of spectra of solar features, the solar variability in Str{\"o}mgren  $(b+y)/2$ photometry is significantly lower at 0.37 mmag (with facular and spot components being equal to 1.2 and -0.83 mmag, respectively).

\subsection{Timescale of solar rotation}\label{sub:rot}
In Fig.~\ref{fig:prof_rot}a we plot the amplitude of the rotational brightness variability of the Sun and its facular and spot components calculated employing the entire SATIRE-S time series considered in this study (2 February 1999 -- 1 August 2014, see Sect.~\ref{sect:model}). This time series covers more than a half of solar cycle 23 (including the maximum in 2000) and the ascending part of cycle 24.

The amplitude of the solar brightness variability closely follows its facular component shortward of a rather sharp threshold between faculae- and spot-dominated regimes and its spot component longward of this threshold
\citep[see also low panels of Fig.~8 from][who calculated the RMS variability of the SATIRE-S SSI time series between May and August 2004 as well as the spot and facular component of the RMS variability]{yvonne2008}. As a consequence, the rotational solar brightness  variability depends very strongly on wavelength in the UV, but is
monotonous and smooth in the visible and IR spectral domains. In contrast to the 11-year activity cycle timescale, the transition from the faculae- to spot-dominated regimes  does not lead to a cancellation of the facular and spot contributions to solar brightness variability and to a subsequent decrease of the variability's amplitude.

The comparison of Fig.~\ref{fig:prof_act} with  Fig.~\ref{fig:prof_rot}a reveals one more crucial difference between brightness variability  on the rotational  and 11-year activity timescales:  the transition from faculae- to spot-dominated regime in the case of the rotational variability happens  at substantially shorter wavelengths than for the 11-year variability. This has to do with different spatial distributions of bright (i.e. faculae and network) and dark (i.e. spots) magnetic features on the solar surface. Faculae and the network are distributed more homogeneously than the spots and,
consequently, the relative contribution of faculae and network to the brightness variability on the rotational timescale is smaller than on the 11-year timescale.
For example, a magnetic component of the solar atmosphere homogeneously distributed on the solar surface will contribute to the 11-year brightness variability (defined by Eq.~\ref{eq:11}) if the disc area coverage of this feature slowly changes with the solar cycle, but will contribute only negligibly to the rotational brightness variability (defined by Eq.~\ref{eq:rot}). As a result, the transition  from the faculae- to the spot-dominated regime of rotational variability happens already around 400 nm,  immediately longward of the CN violet system \citep[which amplifies the facular contrast and thus ensures the faculae-dominated regime of variability, see discussion in][]{profiles}.

Figures~\ref{fig:prof_rot}b~and~\ref{fig:prof_rot}c present the amplitude of the rotational brightness variability of the Sun and its  facular and spot components calculated employing two annual segments of the SATIRE-S reconstruction, representing conditions of the activity maximum (2000, Fig.~\ref{fig:prof_rot}b) and minimum (2008, Fig.~\ref{fig:prof_rot}c).
Comparison of the amplitudes calculated over the entire period (Fig.~\ref{fig:prof_rot}a) and those over the 2000 segment (Fig.~\ref{fig:prof_rot}b) shows that the relative spot contribution to the SSI variability slightly increases towards activity maximum.

Solar activity was exceptionally low in 2008 (with annually averaged sunspot number being equal to 2.9 and 265 days being spotless; both numbers according to WDC-SILSO, Royal Observatory of Belgium, Brussels). This renders the amplitude of the rotational brightness variability in 2008 to be almost an order of magnitude less than in 2000.

Both CoRoT and Kepler sensitivity functions peaks in the visible part of the solar spectrum (see Fig.~\ref{fig:prof_rot}), where the amplitudes of the brightness variability averaged over 1999-2014 period and 2000 segment are spot-dominated. In 2008 spot and faculae contribution are, however, comparable when integrated over the entire spectral domains of  CoRoT and Kepler sensitivity.

\section{Facular and spot components of solar brightness variability observed out-of-ecliptic}\label{sect:out}
In this Section we investigate how the regimes of  solar brightness variability depend on the position of the observer relative to the ecliptic plane.

The effect of the out-of-ecliptic observations on the 11-year solar brightness variability has been studied in detail by \cite{knaacketal2001} but with a more rudimentary model than the one presented in this study. While  \cite{knaacketal2001} uniformly distributed faculae and sunspots in two latitude bands (one in each solar hemisphere) assuming a fixed ratio between sunspot umbra and penumbra areas, we utilize the observed positions and disc area coverages of faculae, spots umbrae and penumbrae.

In Sect.~\ref{sub:11_out} we briefly discuss the effect of the out-of-ecliptic observations on the 11-year solar brightness variability and compare our results with those of \cite{knaacketal2001}. In Sect.~\ref{sub:rot_out} we show how the observations from an out-of-ecliptic vantage point affect the rotational solar brightness variability. We note that the calculations of the rotational solar brightness variability require knowledge of the spot and facular positions as well as  the day-to-day evolution of their disc area coverages and thus could not be performed with the \cite{knaacketal2001} model.

\subsection{Activity cycle timescale}\label{sub:11_out}
\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{Fig2_angle.eps}}
\caption{The amplitude of the 11-year solar brightness variability (black curves in panel a) and its facular  and spot  components (red and blue curves in panel b, respectively) expressed in mmag. As in Fig.~\ref{fig:prof_rot} the green shaded contours indicate the transmission curves of the Str{\"o}mgren filters $b$ and $y$.  The amplitude values are plotted for four values of solar inclination: 90$^{\circ}$ (thick solid), 60$^{\circ}$ (thin solid), 30$^{\circ}$ (dotted), and 0$^{\circ}$ (dashed). }
\label{fig:prof_act_incl}
\end{figure}

\begin{figure*}
\resizebox{\hsize}{!}{\includegraphics{Fig3.eps}}
\caption{The amplitude of the rotational solar brightness variability (black curves) and its facular (red) and spot (blue) components calculated for the 1999-2014 period at three values of the inclination: 60$^{\circ}$ (panel a),  30$^{\circ}$ (panel b),  and 0$^{\circ}$ (panel c). The dark (light) shaded areas show Kepler (CoRoT) total spectral efficiency.}
\label{fig:prof_rot_incl}
\end{figure*}

In Fig.~\ref{fig:prof_act_incl} we plot the amplitude of the 11-year solar brightness variability (Fig.~\ref{fig:prof_act_incl}a) and its facular and spot components (Fig.~\ref{fig:prof_act_incl}b) calculated at four different inclinations.

When the observer moves from the ecliptic plane, the apparent distribution of magnetic features generally shifts towards the visible solar limb (see Fig.~\ref{fig:images}). The absolute value of the spot brightness contrast decreases monotonically  from the solar disc centre towards the limb \citep[see detailed discussion in][]{knaacketal2001}. Furthermore, the shift of the spot towards the solar limb decreases its apparent disc coverage due to the foreshortening effect. As a result, the absolute value of the spot component decreases monotonically  when the observer moves away from the ecliptic plane (i.e. inclination decreases) independently of the wavelength.

The facular brightness contrast  increases towards the limb in the visible spectral domain, so that the CLV of the facular contrast almost exactly compensates the foreshortening effect. As a result,  the amplitude of the facular component only marginally depends on the inclination in the visible spectral domain. At shorter wavelengths,
the CLV of the facular contrast gets weaker and even changes sign in the UV \citep[i.e. faculae get dimer towards the limb, see discussion in][and references therein]{yeoetal2013} so that the CLV and foreshortening effects do not compensate each other. Consequently, the amplitude of the facular variability in the UV noticeably decreases when the observer moves from the ecliptic plane.

As a result the 11-year solar brightness variability decreases with inclination longward of about 400 nm and increases with inclination at shorter wavelengths (see  Fig.~\ref{fig:prof_act_incl}a). Such a behaviour of the 11-year solar brightness variability was pointed out by  \cite{knaacketal2001} (see their Fig.~7).  \cite{knaacketal2001} estimated that the 11-year variability of solar Str{\"o}mgren $(b+y)/2$ flux increases by a factor of 1.3 when observed at inclination  $i=57^{\circ}$ and by a factor of 2.2 when observed at $i=0^{\circ}$, i.e. along the rotational axis (both numbers are relative to the brightness variability observed from the ecliptic plane). Our calculations yield similar factors of 1.35 and 2.7, respectively. The small difference between our results and those of \cite{knaacketal2001} indicate that the simplification of the distribution of magnetic features on the solar surface employed  by \cite{knaacketal2001} for the solar case, and later by \cite{Shapiro2014_stars} for more active stars, has only a minor effect on the calculations of the 11-year brightness variability.

\begin{figure*}
\resizebox{\hsize}{!}{\includegraphics{Fig4.eps}}
\caption{Solar total (black curves),  facular (red), and spot (blue) light curves as they would appear to Kepler. The light curves are plotted for four values of solar inclination:  90$^{\circ}$ (panel a), 60$^{\circ}$ (panel b), 30$^{\circ}$ (panel c), and 0$^{\circ}$ (panel d). The mean levels of all plotted light curves (indicated by the dashed lines) have been offset relative to each other for clarity,  so that only the variations of the brightness, but not the offsets, have a physical meaning.}
\label{fig:LC}
\end{figure*}

Our calculations indicate  that the 11-year variability of the solar brightness is faculae-dominated in the UV and visible spectral domains independently of the inclination.

\subsection{Timescale of solar rotation}\label{sub:rot_out}
In Fig.~\ref{fig:prof_rot_incl} we present  the amplitude of the rotational solar brightness variability and its facular and spot components as a function of wavelength calculated at three out-of-ecliptic positions of the observer. The amplitudes plotted in Fig.~\ref{fig:prof_rot_incl} refer to the entire period considered in this study (2 February 1999 - 1 August 2014).

Comparing Figs.~\ref{fig:prof_act_incl}~and~\ref{fig:prof_rot_incl} shows that the shift of the observer from the ecliptic plane affects facular and spot components of the solar brightness variability on the rotational timescale in a very similar way as it does on the 11-year timescale. In particular, while the facular component  only marginally depends on the inclination in the visible domain of the spectrum, the strength of the spot component noticeably decreases when the observer moves away from the ecliptic plane.

Despite such a similarity in the variability of the individual components  on the 11-year and rotational timescales, the behavior of the total brightness variability on these two timescales differs. The 11-year variability in the visible  is faculae-dominated so that the decrease of the spot component for the out-of-ecliptic observer causes an increase of the brightness variability (see Fig.~\ref{fig:prof_act_incl}). In contrast, the rotational variability  in the visible spectral domain is spot-dominated, so that the decrease of the spot component leads to weaker overall variability (see Fig.~\ref{fig:prof_rot_incl}).
Thus the amplitudes of solar brightness variability on the rotational and activity cycle timescales show opposite trends with the inclination $i$. Consequently special care is needed when extrapolating the dependence of the brightness variability on the inclination from the 11-year to the rotational timescale  \citep[cf. e.g.][]{Gilliland2011,Harrison2012,Lewis2013} and vice versa.

\section{The solar brightness variability as it would be observed in broadband fllters}\label{sect:K-C-out}
In this Section we simulate the light curves of the Sun as they would be observed by Kepler at various inclinations.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{corr.eps}}
\caption{{  The dependence of the Pearson correlation coefficients between spot and total light curves ($\rho_{\rm tot, spot}$, blue symbols) and between facular and total light curves ($\rho_{\rm tot, fac}$, red symbols) on inclination.
The correlation coefficients are calculated for solar brightness variations in 2002 as they would be observed in the Kepler passband (asterisks) and in the Str{\"o}mgren $(b+y)/2$ photometry (triangles).}}
\label{fig:corr}
\end{figure}

To calculate total, facular, and spot fluxes as they would be measured by Kepler ($F^i(t)$, $F^i_{\rm fac}(t)$, and $F^i_{\rm spot}(t)$, respectively), we convolve $S^i(\lambda, t)$, $S^i_{\rm fac}(\lambda, t)$, and $S^i_{\rm spot}(\lambda, t)$  time series  (see Sect.~\ref{subsect:out}) with the Kepler spectral efficiency profile  $\cal{P} (\lambda)$. These time series are then used to calculate the light curves of solar brightness with daily cadence,  $F^i(t)/F^i(t_0)$ (hereafter, total light curves),  as well as the facular $F^i_{\rm fac}(t)/F^i(t_0)$ and spot $F^i_{\rm spot}(t)/F^i(t_0)$ light curves (i.e. light curves, which only account for the brightness variations caused by faculae and spots, respectively).
Here  $t_0$ is the time of a reference measurement.  The choice of $t_0$  is not important, since only the variation of the solar brightness, but not its absolute value is essential for the present study.

In Fig.~\ref{fig:LC} we show solar total, facular, and spot light curves calculated for 2002, a year of relatively high solar activity (with annual sunspot number being equal to 104) at four different inclinations. In line with the discussion in Sect.~\ref{sub:rot} solar brightness variations as they would be measured by Kepler from the ecliptic plane ($i=90^{\circ}$) are essentially driven by spots, i.e. total and spot light curves change in phase (the Pearson correlation coefficient between total and spot light curves for 2002 is $\rho_{\rm tot, spot}^{90}=0.82$, while it is only  $\rho_{\rm tot, fac}^{90}=0.35$ between the total and facular light curves).

The shift of the observer $30^{\circ}$ out-of-ecliptic towards the South ($i=60^{\circ}$) has only marginal effect on the facular light curve. Interestingly some of the dips due to sunspots remain unchanged (e.g. in April 2002) while others are strongly reduced at $i=60^{\circ}$ (e.g. in July and August 2002). The former are caused by spot groups in the Southern hemisphere, while the latter are due to Northern hemisphere spots (note that the observer is now looking at the Sun from the south).

It also has a small effect on the dips in the spot light curve caused by sunspot groups in the South solar hemisphere (e.g. in April 2002). At the same time such a shift of the observer noticeably affects the dips caused  by sunspot groups in the North solar hemisphere (e.g. in July and August 2002). Overall brightness variations at $i=60^{\circ}$  are mainly associated with spots, although the relative contribution of faculae is slightly higher compared to observations from the ecliptic plane ($\rho_{\rm tot, spot}^{60}=0.7$ and  $\rho_{\rm tot, fac}^{60}=0.49$). For $i=30^{\circ}$ and
 $i=0^{\circ}$, the spot contribution is significantly reduced while facular contribution is similar to the $i=90^{\circ}$ case. Consequently, solar brightness variations for $i=30^{\circ}$ and $i=0^{\circ}$ are mainly due to faculae  ($\rho_{\rm tot, spot}^{30}=0.41$, $\rho_{\rm tot, spot}^{0}=0.24$ and $\rho_{\rm tot, fac}^{30}=0.7$, $\rho_{\rm tot, fac}^{0}=0.83$) as explained in Sect.~\ref{sub:rot_out}.

{  In Fig.~\ref{fig:corr} we plot the dependence of $\rho_{\rm tot, fac}$ and  $\rho_{\rm tot, spot}$ values on inclination. For high inclinations $\rho_{\rm tot, spot} > \rho_{\rm tot, fac}$, i.e. for the observer located close to the ecliptic, the solar brightness variations are mainly associated with spots. The situation reverses at about $i=45^{\circ}$ and the day-to-day solar brightness variability is mainly brought about by faculae at low inclinations. Contribution of spots to solar brightness variations is slightly larger when it is observed in the Str{\"o}mgren $(b+y)/2$ photometry than in the Kepler passband. We note that Kepler and CoRoT have very similar passbands (see the shaded areas in  Fig.~\ref{fig:prof_rot} and Fig.~\ref{fig:prof_rot_incl}) so that the Kepler and CoRoT curves will be hardly distinguishable in Fig.~\ref{fig:corr}.
 }

{  In Fig.~\ref{fig:RMS} we plot the RMS values of solar brightness variations in the Str{\"o}mgren $(b+y)/2$ photometry and in the Kepler passband as a function of inclination.  One can see that the amplitude of solar rotational brightness variability observed in Kepler passband decreases by a factor of about 2.5 when the observer moves from the ecliptic plane ($i=90^{\circ}$) to the solar rotation axis ($i=0^{\circ}$). A larger value of about 6 was recently found by \cite{Borgniet2015} for the decrease of the TSI variability.
The difference between our result and that of \cite{Borgniet2015}  might be due to the difference in utilised distributions of magnetic features and their contrasts. We note that effect of the inclination on the 11-year TSI variability found by \cite{Borgniet2015}  is also different from that of \cite{knaacketal2001} \citep[see discussion in ][]{Borgniet2015}.

Interestingly, our value of the RMS variability of solar brightness observed from the ecliptic in the Str{\"o}mgren $(b+y)/2$ photometry (about 0.33 mmag) is close to the value given by \cite{radicketal1998} for the solar short-term photometric variations (about 0.5 mmag, see their Fig.~7). A small difference between two values might be attributed to different reference periods. Our value refers to the 1999--2014 period, when the solar brightness variations were generally lower than during the 1980--1989 period (see Fig.~\ref{fig:KL}), the value that \cite{radicketal1998} refer to. We note that scaling the TSI variability to calculate the variability of the solar Str{\"o}mgren $(b+y)/2$ flux, broadly employed in the literature \citep[and utilised by][]{radicketal1998}, works reasonably on the rotational timescale, but leads to a substantial overestimate of the variability on the 11-year activity timescale (see Sect.~\ref{sub:11}).

}

\cite{basrietal2011,basrietal2013} introduced a useful metric of stellar photometric variability  $R_{\rm var} (t_{\rm len})$, which describes the range of photometric changes in a light curve over a given period in time $t_{\rm len}$. It is calculated by sorting all flux points for the light curve
and taking the difference between the 95-th and 5-th percentile. \cite{basrietal2013} utilised  $R_{\rm var} (30\,\, {\rm days})$ values. Here we follow up on their approach and  divide the $F^i(t)$, $F_{\rm fac}^i(t)$, and $F_{\rm spot}^i(t)$ time series calculated for 2 February 1999 -- 1 August 2014 period into  30-day segments. For each segment we calculate $R_{\rm var} (30\,\, {\rm days})$  value by subtracting the second lowest brightness value from the second highest. We note that since our model is free from the spurious outliers, which plague measured lights curves, subtracting the weakest value from the strongest would give a very similar result.
Individual $R_{\rm var} (30\,\, {\rm days})$ values correspond to different phases of the solar activity cycle  and thus are well-suited for a comparison with Sun-like stars which are observed at arbitrary phases of their cycles  \citep[see detailed discussion in][and, in particular, their Fig.~5]{basrietal2013}.

Figure~\ref{fig:tail} shows values of the tail distribution functions for $R_{\rm var} (30\,\, {\rm days})$ calculated for total, facular, and spot solar light curves. The tail distribution functions give the fraction of $R_{\rm var} (30\,\, {\rm days})$ values which are larger than a certain value. In the case of the Earth-based observer (Fig.~\ref{fig:tail}a) the tail distribution function of the total solar brightness variations closely follows the facular component at small amplitudes but abruptly switches to follow the spot component at an amplitude of about 0.15 mmag. This suggests that the  variations of solar brightness with small amplitude are associated with faculae, while variations with larger amplitudes are brought about by the passage of spots.

The behaviour of the tail distribution functions at $i=60^{\circ}$ (Fig.~\ref{fig:tail}b) is similar to the case of the Earth-based observer, but the transition from faculae-driven to spot-driven variations happens at slightly higher amplitude of about 0.2 mmag. In contrast, the solar brightness variations at inclination  $i=30^{\circ}$  (Fig.~\ref{fig:tail}c)  are mainly faculae-dominated and only the strongest changes (stronger  than $\approx$0.4--0.5 mmag) are due to spots. At $i=0^{\circ}$ (see Fig.~\ref{fig:tail}d)  the solar brightness variations are caused by faculae independently of the amplitude of variations. As in the case of the 11-year variability,  spots counteract the facular contribution, slightly decreasing the total variations. We note that at  $i=0^{\circ}$  the variations are solely caused by the evolution of magnetic features.

\begin{figure}
\resizebox{\hsize}{!}{\includegraphics{RMS.eps}}
\caption{{  The RMS variability of total (black symbols), facular (red symbols), and spot (blue symbols) solar light curves as they would be observed in the Kepler passband (asterisks) and in Str{\"o}mgren $(b+y)/2$ photometry (triangles). The RMS values are calculated for the 2 February 1999 -- 1 August 2014  period.}}
\label{fig:RMS}
\end{figure}

\begin{figure*}
\resizebox{\hsize}{!}{\includegraphics{Fig5.eps}}
\caption{The tail distribution functions (see the definition in the main text) of solar brightness  variations as they would be measured by {  Kepler (solid curves), CoRoT (dotted curves in panel b), and in Str{\"o}mgren $(b+y)/2$ photometry (dotted curves in panel d)}. Black, red, and blue curves correspond to total, facular, and spot light curves, respectively. The variations are measured by $R_{\rm var} (30\,\, {\rm days})$ and are given in parts per thousand (ppt, 1 mmag $\approx$ 0.92 ppt). Plotted are the tail distribution functions calculated at four values of solar inclination: $i=90^{\circ}$ (panel a), $i=60^{\circ}$ (panel b), $i=30^{\circ}$ (panel c), and $i=0^{\circ}$ (panel d). }
\label{fig:tail}
\end{figure*}

We show the tail distribution functions of solar brightness  variations  as they would be measured by CoRoT at  $i=60^{\circ}$ in Fig.~\ref{fig:tail}b. The tail distribution functions of the spot component are basically identical for Kepler and CoRoT, while the facular component is slightly larger for CoRoT. The somewhat stronger effect of the different passband  on the facular component is probably due to the non-monotonic dependence  of the facular component of brightness variability on wavelength, in particular,  the two peaks about 390 nm and 430 nm (associated with the CN violet system and CH G-band, respectively).

\section{Conclusions}\label{sect:conc}
We employed the SATIRE-S model to assess the relative contributions of spots and faculae to solar brightness variations depending on the wavelength, inclination, and timescale of the variability. This allowed us to identify the main drivers of brightness variations of a star identical to the Sun as observed by ground-based or spaceborne telescopes.

On the 11-year activity cycle timescale, the solar brightness variations are faculae-dominated in the UV, visible, and near-infrared spectral domains. In particular, if the Sun were observed in Str{\"o}mgren filters $b$ and $y$ within the Lowell or Fairborn synoptic programs, the variations of its seasonally-averaged brightness would be faculae-dominated. The Str{\"o}mgren filters $b$ and $y$ turn out to be somewhat unfortunate choices for studying the variability of Sun-like stars in that the facular contrast is very small there, resulting in an almost full compensation of the facular and spot contributions to the 11-year solar brightness variations in these passbands. Consequently, the amplitude of the 11-year variability of the solar Str{\"o}mgren $(b+y)/2$ flux  observed from the ecliptic plane is only about 0.4 mmag, which is significantly smaller than the amplitude often assumed in the literature \citep[about 1 mmag, see][]{lockwoodetal1997, radicketal1998, lockwoodetal2007}.
Thus the activity cycle brightness variability of a true solar twin would hardly be detectable with ground-based observations in the Str{\"o}mgren filters $b$ and $y$.

Interestingly, \cite{Lockwood2013_stars} found that a few dozens of stars with chromospheric activity similar to that of the Sun have significantly higher photometric variability than the Sun.  {  It has been shown that this discrepancy cannot be explained by the special position of the Earth-bound observer who always sees the Sun from the near-equator plane \citep{knaacketal2001,Shapiro2014_stars,Borgniet2015}.} This, however, does not necessarily imply that the Sun is anomalous with respect to its stellar cohort. One can expect that the facular contrast in Sun-like stars depends on the stellar brightness temperature \citep[cf.][]{Beecketal2015} and on the stellar metallicity \citep[since the facular contrast is strongly affected by the weak atomic and molecular lines, see][]{profiles}. A small change of the facular contrast might break a delicate balance between facular and spot contributions to Str{\"o}mgren $b$ and $y$ photometry and significantly increase the stellar brightness variability. Hence it is possible that stars with slightly different parameters show brightness variability significantly stronger than that of the Sun.

On the rotational timescale the situation is different.
The RMS solar brightness variations on the rotational timescale averaged over the 1999--2014 period is dominated by sunspots longward of 400 nm and by faculae at shorter wavelengths. Thus solar brightness variations  observed by Kepler or CoRoT from the ecliptic plane would to a large extend be due to spots. The only exception is periods of very low solar activity when the spot contribution is weak and becomes comparable to that of the faculae.

In agreement with previous studies \citep[e.g.][]{knaacketal2001}, a relocation of the observer away from the ecliptic plane  does not change the regime of solar Str{\"o}mgren $(b+y)/2$  flux variations on the 11-year timescale  but {\it increases} the amplitude of such variations (by a factor of up to 2.7 times, for observations along the solar rotational axis).
However, on the rotational timescale a relocation of the observer from the ecliptic plane has an opposite effect on the brightness variations. In general, it {\it decreases} the amplitude of solar brightness variations and, in particular, the amplitude of solar brightness variations as they appear to Kepler and CoRoT  (by a factor of up to 2.5  for observations along the solar rotation axis). The facular contribution to solar brightness variations on the rotational timescale increases with decreasing inclination, so that solar brightness variations on this timescale observed by Kepler and CoRoT at  $i=30^{\circ}$ would be, to a large extent, associated with faculae (with the exception of the strongest variations, with amplitude over 0.4--0.5 mmag, which are produced by large sunspot groups).

Our results indicate that the intercomparison of the observed photometric trends obtained for the cyclic (e.g. with the Lowell and Fairborn ground-based surveys) and for the rotational (e.g. from the Kepler and CoRoT data) timescales must be carried out with caution.
In particular, the observation that the Sun exhibits the same level of the {\it rotational} variability as most  Kepler stars \citep{basrietal2013} does not contradict its relatively low variability on the {\it 11-year activity} timescale \citep{lockwoodetal2007, Lockwood2013_stars}.

While Kepler and CoRoT measure stellar brightness variations with a cadence of 30 minutes and 512 sec (in the nominal regimes), respectively, our model allows calculating the solar brightness values only with  daily cadence. Hence it cannot be applied for simulating solar brightness variations on periods shorter than a day and, in particular, on timescales of several hours  typical for planetary transits. As a next step, we plan to extend our model to simulate the variations of solar brightness on periods shorter than a days as well as to Sun-like stars with various effective temperatures and metallicities.

\begin{acknowledgements}
The research leading to this paper has received funding from the People Programme (Marie Curie Actions) of the European Union's Seventh Framework Programme (FP7/2007-2013) under REA grant agreement No. 624817. It also got financial support  from the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea.
\end{acknowledgements}

\bibliographystyle{aa}