\section{An FU-Orionis-like burst from a high-mass protostar}
\label{sect:introduction}


Accretion-driven luminosity outbursts are a vivid manifestation of variable mass 
accretion onto protostars. They are known as the so-called FU Orionis phenomenon 
in the context of low-mass protostars (Vorobyov \& Basu 2006). More recently, 
this process has been found in models of primordial star formation (Hosokawa et 
al. 2016). In Meyer et al. (2017), using numerical radiation hydrodynamics simulations of a 
collapsing $100\, \rm M_{\odot}$ pre-stellar cores rotating with a ratio of kinetic by 
gravitational energy of 4\% that produces a central massive protostar (cf. Kuiper et al. 2011), we stress 
that present-day forming massive stars also experience variable accretion (Fig.~\ref{fig:disk_density_plots}) and 
show that this process is accompanied by luminous outbursts induced by the 
episodic accretion of gaseous clumps migrating from the circumstellar disk onto 
the protostar (Fig.~\ref{fig:zoom_burst}). 








\section{Observational implications}


We conjecture that luminous flares from regions hosting
forming high-mass star may be an observational
implication of the fragmentation of their accretion disks,
i.e. that those flares constitute a possible tracer of the
fragmentation of their accretion disks. This may apply to
the young star S255IR-NIRS3 that has recently been
associated to a 6.7 GHz methanol maser outburst
(Fujisawa et al. 2015, Stecklum et al. 2016) but also to the
other regions of high-mass star formation from which
originated similar flares (Menten et al. 1991) and which
are showing evidences of accretion flow associated to
massive protostars, see e.g. in W3(OH), W51 and W75.


\begin{figure}
        \centering              
        \begin{minipage}[b]{ 0.57\textwidth}
                        \includegraphics[width=0.98\textwidth]{./zoom.pdf}
        \end{minipage}         
        \caption{ 
                 Accretion rate onto the protostar (in $\rm M_{\odot}\, \rm yr^{-1}$) and total
luminosity of the protostar (in $\rm L_{\odot}$) during the burst. Magenta
dots mark the times of Fig.~\ref{fig:disk_density_plots}b-d.  
                 }      
        \label{fig:zoom_burst}  
\end{figure}


\section*{Acknowledgements}


This study was conducted within the Emmy Noether research group on "Accretion Flows 
and Feedback in Realistic Models of Massive Star Formation" funded by the German Research 
Foundation under grant no. KU 2849/3-1. E.I.V. acknowledges support from the 
Austrian Science Fund (FWF) under research grant I2549-N27 and  RFBR grant 14-02-00719. 


