\section{Background}

\subsection{Children's Perception of Privacy}

Privacy has been described as a ``concept in disarray''\cite{solove2008understanding} reflecting the complex and often overlapping conceptualisations. O'Hara proposed ``seven veils'' of privacy, spanning the conceptual, factual, phenomenological, preferential, societal, legal, and moral~\cite{ohara2016seven}. Practical research in privacy has typically spanned multiple such veils; Livingstone suggests that privacy is related to the keeping of information out of the public domain or the effective determination of \emph{what} personal information is made available to \emph{whom}~\cite{livingstone2006children}. Nissenbaum's theory of Contextual Integrity (CI) frames privacy as  ``the right to appropriate flow of personal information"~\cite{nissenbaum2004privacy} where what is ``appropriate'' is based on particular contexts and relationships.  CI has been found to be a useful practical framework to interpret people's perceptions of privacy~\cite{barkhuus2012mismeasurement,kummer2016private}.

Managing privacy is becoming increasingly challenging, and children are particularly susceptible in part due to them having little sense of the risks posed by the accumulation of personal data over time~\cite{Pangrazio:2017:MDM:3097286.3097338}. \comm{The combination of such a lack of understanding, with incentives put in place by apps and platforms to get people to share their data, have yielded the so-called `privacy paradox', a behaviour studied in teenagers~\cite{brandimarte2013misplaced} and adults~\cite{barth2017privacy} alike, in which individuals act in ways contrary to their stated privacy preferences and concerns~\cite{jones2009us,Jia:2015:RLP:2675133.2675287,Pangrazio:2017:MDM:3097286.3097338}.} Recent studies showed that although teenagers from 14 to 18 were typically concerned about being personally identified by their personal data, they failed to perceive the potential threat of re-identification via the particular fragments they shared, e.g., images or geo-location~\cite{Pangrazio:2017:MDM:3097286.3097338}.

Whilst there has been extensive research on \comm{teenagers'} perceptions of privacy, relatively little has been done for children aged below 11. They have often been perceived as being too young to understand privacy or exercise digital independence online~\cite{livingstone2018parentsc}.  However, developmental research has shown that children can develop a ``theory of mind" by the age of 4 \cite{colwell2016secret}, which enables children to recognise the differences between concepts in their minds and those in others, and hence to grasp complicated concepts like ``secrecy''~\cite{colwell2016secret} or ``deception''~\cite{chandler1989small}. Contrary to common expectations, children value their privacy because this enables them to enjoy their experience online \cite{Silva:2017:PCT:3091478.3091479, kumar2018cscw}, such as socialising with their friends and families or experimenting with new games. Existing literature provided us with us with useful insights into children's perception of risks \cite{zhang2016nosy, kumar2018cscw}, however, they they do not tell us how children make these judgements - or whether children can in fact identify a threat and its origination.

\subsection{Managing Children's Privacy Online}
Tools developed for managing children's privacy on mobile platforms \comm{usually comprise} features that monitor and restrict children's online activities~\cite{wisniewski2017parental,livingstone2018parentsa,ofcom2017}. Children often find these tools overly restrictive and invasive of their personal privacy \cite{Ghosh:2018:SVS:3173574.3173698,wisniewski2017parental}. Co-design studies with children under 12 show that they prefer technologies that facilitate their risk coping skill development, promote communication and interaction with their parents, and emphasise restrictions around monitoring~\cite{mcnally2018co}. However, such tools are still scarce.

Several recent studies have recognised the importance of supporting children's learning through approaches like interactive storytelling~\cite{Zhang-Kennedy:2016:TIE:2930674.2935984}, game playing~\cite{maqsood2018day} or co-learning with parents~\cite{hashish2014involving}. They have shown the effectiveness of increasing children's awareness of related online safety issues, such as online personal identity or content appropriateness. This demonstrates that a knowledge scaffolding approach to  a child's means of dealing with privacy risk coping could provide a useful addition to safeguarding their cybersafety.

\subsection{Developmental Stages of Privacy Concepts}

Children's ability to recognise privacy risks may be intrinsically limited by their developmental stage. Children aged 3-5 start to build up friendships; however, at this stage, family interactions are dominant, and many of their online activities are still parent-guided \cite{digital2017}, peer pressure is less of an influence~\cite{ann2000developmental}. Children from 6 to 9, \comm{however}, start to learn about the complexity of relationships and feel the need to fit in to peer social groups~\cite{ann2000developmental}. They are also more involved in online activities and enjoy playing games with their friends~\cite{digital2017}. This social interaction with peers makes them more aware of interpersonal privacy risks, but less of other privacy risks \cite{digital2017}. So while they often care deeply about their personal information being shared with their peers, parents, and others online, they remain unaware of other actors, including platforms, app designers, malicious actors, and others operating in digital ecosystems~\cite{livingstone2018childrenc,Pangrazio:2017:MDM:3097286.3097338}.

Vygotsky's Zone of Proximal Development (ZPD) is a theory that relates the difference between what learners can do independently to what can be achieved by through guidance by a skilled partner~\cite{chaiklin2003zone}. It has been applied to assess the effectiveness of teaching and learning~\cite{miller2010theories}, to identify key barriers in learning by understanding knowledge gaps~\cite{daniels2005introduction}, and in the design of better intelligent tutoring systems (ITS)~\cite{quintana2004scaffolding}. \comm{Inspired by ZPD, we focus on identifying children's misconceptions of privacy risks and highlighting areas of knowledge that could be developed through future active scaffolding. }
\section{Key Findings}

Through careful data analysis of 20-hour interview recordings, we have achieved some positive findings with respect to children's level of awareness. However, at the same time, we have also identified critical gaps in children's online privacy knowledge. Children sometimes may be able to capture certain risks with the right language, however, they may not have fully comprehended the risks. These are critical insights that should draw further attention to families and schools when applying risk mediation to children. The terms used by the participant children and their abilities to apply these terms to different types of risk context has provided valuable inputs to our future design of technologies to facilitate children's risk coping abilities.

\subsection{Children care about online privacy risks}
Our results reinforced that, like teenagers and younger children, our participants valued the positive experiences of being online and keeping their personal space online~\cite{lwin2008protecting}. They enjoyed going online so that they could keep in touch with friends, learn new things or just have fun.

At the same time, children felt `annoyed', `surprised' or `angry' when they felt coerced or not in control. As a result, children cared about, and were sensitive to, who might access their sensitive information (e.g. real names, age, location etc), and applied a range of techniques to safeguard this space, such as by verifying identities through face-to-face interactions or avoiding using real names as usernames.

\subsection{Children still need help with certain privacy risks}

Children demonstrated a strong consciousness of their online identity and the importance of avoiding sharing their real identity or over-sharing their personal information. In these cases, children applied various effective strategies to protect their sensitive personal information, such as `making up a new name', and asking their friend to verify the online invitation in person. However, in contrast to explicit privacy risks, children struggled to associated ``online promotions'' with losing control of personal information.

Several children correctly discussed their interpretation of how YouTubers may try to `persuade' them to watch their videos in order to gain `money' or `more subscribers'. However, none of these children expressed resistance to these video promotions. They treated it as if this is how the Internet works.

Only a few children (3/12 focus groups) recognised how new videos might be related to personalised recommendations. However, they struggled to understand who was performing these recommendations, and as a result, they were less sure about the consequences of their privacy. As a result, the children usually applied the play-and-see strategy to assess the content or apps they came across, without any more complex reasoning.

\subsection{Children's ability to describe risks depends on their actual understanding of risks}

Our study is partially inspired by Vygotsky's Zone of Proximal Development (ZPD), and therefore, we really focused on understanding children's current knowledge by observing the language they used to describe risks. We observed that some terms were repeatedly used by children across different situations --- children were able to describe risks accurately when they could recognise the actual risks. However, when they had only a vague understanding of the risks, they struggled to describe things consistently or to provide a good explanation of what they meant.

When children only made a partial sense of certain types of risks, their use of terms can be inconsistent. For example, although some children were able to describe the scenario of new videos being presented to Bertie using words like `people trying to make money' or `them trying to make you watch more', they struggled to explain who these `people' are and how information might be transmitted to these `[YouTube] channel people' in this context.

Another example is the term `hacking', which has been used by children across different focus groups, but in fact with very different meanings. For example, when trying to explain why a new video was shown following a previous video, children used `hacking' to mean `someone stole my data', `take your account', or `steal from house'. `Hacking' has been used by the same group of children to make sense of other scenarios, including in-app pop-ups or data tracking.

\subsection{Children struggle to comprehend the implications of risks}

Children could miss risks due to a lack of knowledge, or due to their past experiences, which did not lead to any direct consequences related to the risks. For those children who interacted with certain technologies and experienced no implications before, they would tend to be more (over-) confident with technologies.

Children from 7 focus groups treated online video promotions as part of an `autoplay' function of the platform, without questioning how the new content might be presented to them. 12/29 children demonstrated trust in the content provided by their familiar YouTubers. As a result, children reported having been exposed to unexpected content and online baiting. These same children reported that they often saw upsetting content online (e.g. `sometimes in autoplay, it comes up with these really freaky ones like pictures of dead people').

When children not fully comprehending the risks, we also observed that a child's or their peers' personal experiences had a strong influence upon their decision making, even though they didn't always understand what may pose threats to them.

\newpage

\frame{
\textbf{Take-Home Messages For Parents \& Caretakers}
\begin{enumerate}
    \item Continue to talk to children about being careful online, because they are indeed facing challenges in their use of digital devices every day and they really care about protecting their personal space online!

    \item Pay attention to how children describe risks or things that made them uncomfortable online. This is a good start to get to know what they actually mean, and it's always good to follow this up with a question like ``what do you mean?''.

    \item Do not be afraid of talking about digital technologies with children. They would enjoy learning with you and sharing their knowledge with you.

    \item If you are indeed struggling, please take a look at the resources we recommend below. It's never too late to get to know what your child knows about online privacy and how you may help them or how you may help each other!
\end{enumerate}
}

\section{Resources}
\begin{itemize}
    \item \href{commonsensemedia.org}{commonsensemedia.org} would be a good place for finding out additional information about which mobile app would be good for your child, such as Minecraft, Roblox, etc.

    \item \href{https://www.net-aware.org.uk}{NSPCC Net Aware} provides good evidence-based advice on apps that are most popular among children.

    \item \href{childnet.com}{childnet.com} has good resources for talking about online privacy with younger children, such as the story of \href{https://www.childnet.com/resources/smartie-the-penguin}{Smartie the Penguin} or the story book of  \href{https://www.childnet.com/resources/digiduck-stories}{Digiduck}.

    \item Finally, our KOALA Hero is always here to help, either through the \href{https://sites.google.com/view/koala-project-ox/toolkits/cards}{story cards} or regularly updated \href{https://sites.google.com/view/koala-project-ox/reports-and-publications}{public reports}.
\end{itemize}

\section{Introduction}
Today, children are spending more time online than with other media sources, such as televisions or offline video games~\cite{ofcom2017,livingstone2017children}. Among the many kinds of devices now connected to the Internet, mobile devices (such as tablet computers or smartphones) have become the primary means by which children go online~\cite{livingstone2017children}. In the UK, 35\% of children aged five to seven and 52\% of children aged eight to eleven have been reported to have their own tablets  (see Figure 1)~\cite{ofcom2017}; while in the US, ownership of tablets by children in this age group grew five-fold between 2011 and 2013~\cite{commonsense2013}. Children under five are also using smartphones and tablets more often, as the category of apps designed for younger kids continues to expand rapidly~\cite{ofcom2017, lifewire}.

\begin{figure}[h!]
\centering
\begin{bchart}[step=20,max=100,unit=\%]
        \bcbar[label=3-4]{21}
            \smallskip
        \bcbar[label=5-7]{35}
            \smallskip
        \bcbar[label=8-11]{52}
            \smallskip
        \bcbar[label=5-11]{49}
            \smallskip
        \bcbar[label=12-15]{55}

\end{bchart}
\caption{Percentage of Children Owning Their Own Tablets In Proportion to the Total Number of Tablet Accessible to the Household by age~\cite{ofcom2017}}
\end{figure}\label{fig:devices}

However, for children aged between 6 and 10, their awareness of online privacy has been seldom studied in existing literature. It's probably due to the expectation that children of this age cannot comprehend the concept of privacy, and therefore their online safety and privacy should be safeguarded by the trusted grown-ups. However, our previous reports~\cite{DBLP:journals/corr/abs-1809-10841,DBLP:journals/corr/abs-1809-10944} have shown that although parents largely take a protective approach to restrict and monitor what their children can or cannot access on their children's devices, parents do not always fully comprehend all risks. Furthermore, a more recent study led by the LSE has shown that parents struggle to seek for help to mediate their children's use of digital devices, and they largely rely on themselves to figure things out~\cite{livingstone2018parentsb}.

Therefore, the goal of the KOALA project is to, through a series of studies directly interacting with children and their families, to achieve a deeper understanding of the following questions:
\begin{itemize}
    \item How well can children aged 6-10 understand online privacy risks related to their use of tablet computers?
    \item What are the barriers for their parents to mediate children's use of technologies and recognition of risks?
\end{itemize}

Between June and August 2018 we conducted 12 focus group interviews with 29 children, aged 6-10, from UK schools. We found that children in our study had a good understanding of risks related to \textit{inappropriate content, the approach of strangers}, and \textit{oversharing of personal information online}. However, they struggled to fully understand and describe risks related to \textit{online game/video promotions and personal data tracking}. Moreover, children's risk coping strategies depended on their understanding of the risks and their previous experiences: effective risk strategies were applied only if children recognised certain risks or when they \textit{felt} something untoward. These findings demonstrate the importance of learning about potential risks through a multitude of channels, such as school, parents, friends, siblings, and personal experiences~\cite{zhang2017cyberheroes,hashish2014involving}.
\section{Acknowledgement}

We thank all the schools, families and children who have contributed their time and knowledge to our study. Without their support, this study would not have been successful!

The report is based on a full academic publication \href{https://arxiv.org/abs/1901.10245}{`I make up a silly name': Understanding Children's Perception of Privacy Risks Online}, as part of the project of KOALA (https://sites.google.com/view/koala-project-ox/): Kids Online Anonymity \& Lifelong Autonomy, funded by EPSRC Impact Acceleration Account Award, under the grant number of EP/R511742/1.

\section{Contact}
\begin{itemize}
    \item Jun Zhao: jun.zhao@cs.ox.ac.uk
    \item Nigel Shadbolt: nigel.shadbolt@cs.ox.ac.uk
\end{itemize}

\section{Methodology}
We chose the focus group method to elicit children's responses to a collection of hypothetical scenarios that reflect different types of explicit and implicit threats to children's online personal data privacy.

Each focus group study contained four parts, including a warm-up and introduction session, a sharing of favourite apps, a walk-through of three hypothetical risk scenarios, and finally an open-ended session about issues not so far discussed.

Children were encouraged to share their personal experience related to the scenarios, by responding to questions like \textit{whether this has happened to them} and \textit{what they did}. This enabled us to listen to children's descriptions of their experiences not covered in the scenarios.

The hypothetical scenarios used in the study have been created into story cards that can be reused and shared in the classrooms or at home. They can be found in the Appendix towards the end of the report, accompanied by descriptions about the privacy risks related to each card. Parents and teachers can also download our story cards from \href{https://sites.google.com/view/koala-project-ox/toolkits}{our project web site} and use these cards to inspire conversations with children about online privacy issues.


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

R~Leo is an O-rich Asymptotic Giant Branch star (AGB) located at $\simeq70-85$~pc from Earth \citep{debeck_2010,ramstedt_2014}.
It has an effective temperature of $\simeq2500-3000$~K and an angular diameter of $\simeq0\farcs025-0\farcs030$ in the K-band \citep{perrin_1999,fedele_2005,wittkowski_2016}.
It pulsates with a period of $\simeq310$~days ejecting processed material with a low mass-loss rate of $\simeq1.0\times10^{-7}$~M$_\odot$~yr$^{-1}${} and a terminal velocity of $6-9$~km~s$^{-1}${} \citep{bujarrabal_1994,cernicharo_1997}.

The continuum at different wavelengths (optical, infrared, mm) together with the thermal and maser molecular emissions coming from this star have been analysed with different observing techniques \citep*[e.g.,][]{castelaz_1997,pardo_1998,ryde_1999,gonzalez-delgado_2003,soria-ruiz_2007}.
Some of these observations suggested the existence of complex structures in the innermost envelope.
\citet{cernicharo_1994} carried out lunar occultation observations of the SiO($v=1,J=2-1$) maser finding a low velocity structure within the 0\farcs5-sized region around the star that could be interpreted as a developing bipolar outflow or a rotating torus.
However, the limited spatial information perpendicular to the occultation direction achieved was insufficient to properly constrain this structure.
Additional asymmetries were found at larger scales by \citet{plez_1994} suggesting that the structure revealed by \citet{cernicharo_1994} could be part of a bigger one covering a significant fraction of the envelope.

In this Letter, we present new interferometer observations of CO and $^{29}$SiO toward the AGB star R~Leo carried out with the Atacama Large Millimeter Array (ALMA) with an angular resolution of up to 0\farcs026, comparable to the size of the star.

\section{Observations}
\label{sec:observations}

\begin{table*}
  \caption{Observations Summary\label{tab:table1}}
  \centering
  \begin{tabular}{c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c@{\hspace{0.75ex}}c}
    \hline
    Conf./Run & Date        & Phase & \multicolumn{2}{c}{PSF}   & W     & MRS & \multicolumn{2}{c}{Continuum} & \multicolumn{3}{c}{Gaussian Fit}\\
              &             &       & HPBW                      & P.A.  &     &           & Peak              & Integral      & Peak       & D. FWHM                  & D. P.A.\\
              &             &       & (\arcsec$\times$\arcsec)  & (deg) &     & (\arcsec) & (mJy~beam$^{-1}$)        & (mJy)         & (mJy~beam$^{-1}$) & (\arcsec$\times$\arcsec) & (deg)\\
    \hline
    C40-2     & 2017-Mar-22 & 0.80  & $1.560\times1.080^a$      & 80    & N   & 11.6      & $91.16\pm0.12$    & $105.8\pm0.6$ & ---        & ---                      & ---\\
    C40-5     & 2017-May-03 & 0.94  & $0.387\times0.359^b$      & 45    & N   & 4.5       & $91.10\pm0.16$    & $111.6\pm1.1$ & ---        & ---                      & ---\\
    C40-6     & 2016-Oct-01 & 0.25  & $0.186\times0.157^c$      & 32    & N   & 1.6       & $101.0\pm0.3$     & $140.2\pm2.1$ & ---        & ---                      & ---\\
    C43-9/A   & 2017-Sep-21 & 0.40  & $0.096\times0.035^d$      & 45    & N   & 0.4       & $68.19\pm0.20$    & $103.1\pm1.6$ & ---        & ---                      & ---\\
              &             &       & $0.071\times0.025^e$      & 32    & U   &           & $52.55\pm0.18$    & $103.3\pm1.8$ & ---        & ---                      & ---\\
    C43-9/B   & 2017-Oct-27 & 0.50  & $0.060\times0.046^f$      & 134   & N   & 0.6       & $45.10\pm0.21$    & $90\pm3$      & 37.8       & $0.042\times0.034$       & 173\\
              &             &       &                           &       &     &           &                   &               & 7.7        & $0.190\times0.021$       & 50 \\
              &             &       & $0.026\times0.026^g$      & ---   & U   &           & $21.38\pm0.16$    & $85\pm3$      & 16.5       & $0.045\times0.030$       & 155\\
              &             &       &                           &       &     &           &                   &               & 5.4        & $0.132\times0.022$       & 50\\
    \hline
  \end{tabular}
  \tablefoot{From left to right, the columns contain the ALMA configuration and the run for duplicated observations, the observation date, the optical pulsation phase (pulsation period~$=310$~days; JD$_\subscript{max}=2457896$ following AAVSO\footnote{\texttt{https://www.aavso.org/}}), the size and P.A. of the synthesised PSF, the visibilities weighting (N, ``natural'', means \texttt{robust}~$=1$ and U, ``uniform'', \texttt{robust}~$=0.1$ in GILDAS), the Maximum Recoverable Scale, the Continuum Peak for every configuration, the Continuum Integral above the $5\sigma$-level contours, the Peak Emission of the Gaussian components of the fit, and the FWHM and P.A. after deconvolution.
    If the Gaussian fit is not provided (---) the source can be considered as point-like.
    The errors of the Continuum Peak and the Continuum Integral do not consider the flux calibration uncertainties.
    The statistical errors of the Gaussian fit are always $\lesssim5\%$.
  (\textit{a-g}) The characteristic HPBWs, $\theta_\subscript{PSF}$, are 1\farcs30, 0\farcs37, 0\farcs17, 0\farcs058, 0\farcs042, 0\farcs052, and 0\farcs026, respectively.}
\end{table*}

We observed R~Leo with ALMA during Cycle 4 and 5 under the frame of project 2016.1.01202.S.
Array configurations C40-2, C40-5, C40-6, and C43-9 were used providing baselines from 15~m up to 13.9~km.
These observations give a complete view of the envelope of R~Leo at different scales along $\simeq1.25$ pulsation periods.
The observation details can be found in Table~\ref{tab:table1}.

Four spectral windows covered the frequency ranges $212.7-216.8$ and $227.5-231.5$~GHz with a channel width of 488~kHz.
The flux, bandpass, and pointing were calibrated in the usual way by observing J0854+2006 or J1058+0133.
We estimate a flux error $\simeq5\%$.
The phase calibrators and check sources J1002+1216, J1008+0621, and J0946+1017 were periodically observed.
R~Leo was observed twice with configuration C43-9.
In the first run (run A) the baseline of one spectral window was affected by strong spurious periodic spikes, the weather was worse, and the $uv$ plane coverage was limited.
The data quality highly increased during the second run (run B) and run A was deprecated.
The data calibration was done with the pipeline in CASA~4.7.2 \citep{mcmullin_2007}.

Mapping and data analysis were almost fully performed with GILDAS\footnote{\texttt{https://www.iram.fr/IRAMFR/GILDAS}}.
Imaging restoration was done adopting a robust parameter of 1.0 and 0.1 (natural and uniform weightings, hereafter).
We estimate a positional uncertainty for the highest angular resolution of $\simeq0\farcs005$ from the standard deviation of the emission peak position of the phase calibrator over time \citep{menten_2006,fonfria_2014}.

In this Letter, we have focused on the short scale observations to describe the stellar vicinity.
The analysis of the envelope through the larger scale maps is out of the scope of this work and will be presented elsewhere (Fonfr\'ia et al., \textit{in preparation}).
Due to the stellar pulsation phase incoherence and the long time lag ($\simeq1.1$~yr) the data taken with the C40-6 and C43-9 configurations have not been merged to prevent artefacts.

\section{Results and discussion}
\label{sec:results}

\subsection{The structure of the continuum emission}
\label{sec:continuum}

\begin{figure}[hbt!]
\centering
\includegraphics[width=0.475\textwidth]{f1.pdf}
\caption{Continuum emission observed with configuration C43-9, run B assuming natural weighting (\textit{lower}; $\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs052$) and uniform weighting (\textit{upper}; $\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs026$).
The lowest contours are at $\pm5\sigma$ (2.3\% and 3.75\% of the maximum for the lower and upper insets).
The rest are at 5, 10, 30, 50, 70, and 90\% of the maximum.
The grey crosses indicate where the emission peaks in the uniformly weighted map.
The star is plotted in the upper inset as a dashed ellipse. Its size is the deconvolved size of the compact emission component.}
\label{fig:f1}
\end{figure}

The observations show one continuum emission source in the region of the sky covered by the primary beam (Fig.~\ref{fig:f1}) peaking at $(\textnormal{RA},\textnormal{DEC})=(09^h47^m33.\!\!^s4915\pm0.\!\!^s0004,11^\circ25^\prime 42\farcs899\pm0\farcs005)$.

The highest resolution observations show a compact source surrounded by a faint extended brightness distribution elongated roughly along the Northeast-Southwest direction (Fig.~\ref{fig:f1}).
It is noteworthy that the Maximum Recoverable Scale (MRS; $\simeq 0\farcs6$) is similar to the size of the extended emission in Fig.~\ref{fig:f1}, which is $\simeq0\farcs35$ above the $5\sigma$ level.
Thus, the array could be filtering emission, preventing us to derive reliable brightness temperatures for the star.

The continuum emission can be described by two Gaussian-like components: one compact and one extended (Table~\ref{tab:table1}).
We can estimate the size of the continuum source at 7~mm at the same pulsation phase as during our observations ($\simeq 0.50$) from the results by \citet{matthews_2018} assuming first an uniform disk fit, using linear interpolations, and finally scaling the result to get a Gaussian size.
The result is $\simeq (0\farcs035\pm 0\farcs006)\times(0\farcs025\pm 0\farcs003)$, smaller than the size of our compact source ($0\farcs045\times0\farcs030$).
This only can occur if the star is surrounded by a continuum source that emits more at 1~mm than at 7~mm like a dusty shell \citep{norris_2012}.
All this is compatible with the extended continuum component (size~$\simeq 0\farcs05$) found by \citet{paladini_2017} in the mid-IR around the star.
For simplicity, we hereby consider this compact component as the star.

The extended Gaussian-like component is elongated along the Northeast-Southwest direction.
Its deconvolved size is smaller than the angular resolution provided by the C40-6 configuration ($\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs17$) and it is not evident in the corresponding map.
The ellipticity of this component is $e=1-\theta_\subscript{min}/\theta_\subscript{maj}\simeq0.83$, where $\theta_\subscript{maj}$ and $\theta_\subscript{min}$ are the deconvolved major and minor axes (Table~\ref{tab:table1}).
Since no extended emission has been discovered at longer wavelengths \citep{reid_2007,matthews_2018} the extended component would be produced by dust.

\subsection{Complex molecular line profiles in the vicinity of the star}
\label{sec:complexity}

The CO and $^{29}$SiO lines at the central pixel reveal profiles with emission and absorption components, suggesting the presence of a complex structure (Fig.~\ref{fig:f2}).
The CO($v=0,J=2-1$) and $^{29}$SiO($v=0,J=5-4$) lines observed with an angular resolution $\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs17$ (Table~\ref{tab:table1}) comprise a main emission contribution roughly centred at a systemic velocity of $-0.5$~km~s$^{-1}${} \citep{teyssier_2006} and a prominent red-shifted absorption of $\simeq25\%$ of the continuum peaking at $\simeq8$~km~s$^{-1}$.
The $^{29}$SiO line observed with $\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs052$ is dominated by a strong absorption peak at $\simeq-2.5$~km~s$^{-1}${} ($\simeq70\%$ of the continuum) while the CO line shows a clear inverse P-Cygni profile.
Line CO($v=1,J=2-1$) also shows a similar profile with its absorption component peaking at $\simeq10$~km~s$^{-1}$.
The $^{29}$SiO($v=1,J=5-4$) line lacks emission but shows a broad absorption peaking at $\simeq3-6$~km~s$^{-1}$.
All these absorptions extend up to velocities $\simeq10-15$~km~s$^{-1}$, higher than the terminal expansion velocity \citep*[$\simeq6-9$~km~s$^{-1}$;][]{bujarrabal_1994,cernicharo_1997,debeck_2010,ramstedt_2014}.
No absorption is observed for all these lines with $\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs37$ and 1\farcs30 except for $^{29}$SiO($v=1,J=5-4$).
This line displays blue- and red-shifted absorptions and a well differentiated, narrow emission component at the systemic velocity.
The large variations seen between the $v=1$ lines, expected to arise very close to the star, observed with $\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs17$ and 0\farcs052 indicate substantial excitation changes throughout a single pulsation period.

\begin{figure}[hbt!]
\centering
\includegraphics[width=0.475\textwidth]{f2.pdf}
\caption{Spectra in the central pixel of the CO($2-1$) and $^{29}$SiO($5-4$) lines ($v=0$ in blue and $v=1$ in red) acquired with different ALMA configurations.
  The number of pixels in the PSF is the same for every configuration.
  The grey vertical bands indicate the gas expansion velocity derived from single-dish observations.
  We have assumed a $v_\subscript{sys}=-0.5$~km~s$^{-1}${} \citep{teyssier_2006}.}
\label{fig:f2}
\end{figure}

All this suggest that
(1) there is gas moving away from us in front of the star producing the red-shifted absorptions,
(2) most of the emission comes from directions that do not enclose the star, and
(3) the upper level involved in the $^{29}$SiO($v=1,J=5-4$) line might be strongly drained, considering that this line forms in shells with temperatures $\gtrsim 1500$~K where the $^{29}$SiO($v=0,J=5-4$) do show a noticeable emission.
This effect could invert the populations of the levels involved in the $^{29}$SiO($v=1,J=6-5$) line, which might show maser emission \citep{cernicharo_1991,gonzalez-alfonso_1997}.

\subsection{Gas infall and hints of rotation}
\label{sec:torus}

\begin{figure*}[hbt!]
\centering
\includegraphics[width=\textwidth]{f3.pdf}
\caption{$^{29}$SiO and CO emission in the sub-arcsec region.
(\textit{Upper left}) $^{29}$SiO($v=0,J=5-4$) line ($\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs17$).
The purple contours are at levels 30, 40, and 50\% of the integrated emission peak.
The green dashed-dotted arrows describe the main directions of the total emission.
The solid blue, red, and cyan contours are the blue-shifted, red-shifted, and systemic velocity emissions ($[-15,-0.4]$, $[0.3,7.6]$, and $[-0.4,0.3]$~km~s$^{-1}$).
The dashed red ones describe the red-shifted absorption ($[7.6,15]$~km~s$^{-1}$), which matches spatially up with that of the $v=1$ line.
(\textit{Upper right}) CO($v=0,J=2-1$) line ($\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs026$).
The solid blue and red contours are the blue- and red-shifted emission ($[-15,-3.4]$ and $[1.2,8.4]$~km~s$^{-1}$).
The cyan contours are the emission around the systemic velocity ($[-3.4,1.2]$~km~s$^{-1}$; $v_\subscript{sys}=-0.5$~km~s$^{-1}$).
The dashed red contours are the red-shifted absorption.
The solid and dashed purple contours are the emission and absorption of the CO($v=1,J=2-1$) line ($[-3.3,6.9]$ and $[6.9,15]$~km~s$^{-1}$).
The thin yellow ellipse represents the star and the dashed-dotted pink arrows define the main directions of the continuum emission.
In both insets, the contours are at levels 70, 90, and 99\% of the peak emission/absorption.
The grey crosses position the continuum emission peak.
(\textit{Lower left} and \textit{right}) Position-velocity diagrams of the $^{29}$SiO($v=0,J=5-4$) and CO($v=0,J=2-1$) lines (\textit{left} and \textit{right}, respectively) along the corresponding main directions.
The horizontal dotted lines are the limits between the previous velocity ranges.
The yellow arrows (\textit{left}) point out features in the diagrams typical of the emission of rotating gas.
The blue dashed lines (\textit{right}) connect the peak emissions along the Northeast-Southwest and Southeast-Northwest revealing the dependence of the Doppler shift on the direction.
The contours are at 3, 7, 15, 30, 50, 70, and 90\% ($^{29}$SiO; $9.5\sigma\simeq3\%$) and 30, 50, 70, 87, and 97\% (CO; $5\sigma\simeq3.6\%$) of the peak emissions.}
\label{fig:f3}
\end{figure*}

The CO and $^{29}$SiO emission distribution seen in the vicinity of the star are shown in Fig.~\ref{fig:f3}.
The integrated emission of line $^{29}$SiO($v=0,J=5-4$) mapped with $\ensuremath{\theta_\subscript{PSF}}\simeq0\farcs17$ (Fig.~\ref{fig:f3}, upper left panel) shows a bipolar structure laying roughly along the East-West direction (P.A.~$\simeq100^\circ$).
Lower level contours ($\simeq30-50\%$ of the peak emission) are elongated along the perpendicular direction.
The East-West structure can be separated into different contributions in the velocity range $[-15,7.6]$~km~s$^{-1}$:
(1) a blue-shifted one centred at $\simeq0\farcs1$ to the East of the star,
(2) a red-shifted contribution centred at $\simeq0\farcs2$ to the West, and
(3) another contribution at the systemic velocity covering the same region of the envelope as the previous ones.
In the middle of the red- and blue-shifted contributions there is an absorption in the velocity range $[7.6,15]$~km~s$^{-1}$, also seen in the $^{29}$SiO($v=1,J=5-4$) line, which matches up with the star within the positional uncertainty.

The division of the CO($v=0,J=2-1$) line into velocity intervals reveals a more complex structure (Fig.~\ref{fig:f3}, upper right panel).
The emission at the systemic and positive velocities of the $v=0$ and 1 lines come from the western hemisphere of the envelope.
The red-shifted absorption of the $v=0$ line matches up with the star but it is slightly shifted toward the Southwest for the $v=1$ line.
The blue-shifted emission of the $v=1$ line is mostly absent.
The blue- and red-shifted emissions of the $v=0$ line and the emission of the $v=1$ line are overlapped to the S.
There is a CO emission deficit at velocities $\gtrsim 3$~km~s$^{-1}${} due Northeast.
The whole picture seems to be roughly symmetric around the Northeast-Southwest direction.
This main direction is compatible with that defined by the CO($v=0,J=2-1$) emission mapped with lower spatial resolutions (Fonfr\'ia et al, \textit{in preparation}) and by the continuum extended emission (Fig.~\ref{fig:f1}).
The position-velocity diagram along the Southeast-Northwest direction (Fig.~\ref{fig:f3}, lower right panel) is similar to those of the rotating disks in expansion found in several post-AGB stars except for the infalling gas \citep{bujarrabal_2016,bujarrabal_2017}.

The disagreement between the symmetry axes in the $^{29}$SiO and CO maps could be consequence of an inhomogeneous chemistry favoured by anisotropic physical and chemical conditions.
Excitation effects cannot be ruled out considering the high number of blended lines shown by SiO and its isotopologues displaying maser emission in this and other O-rich stars \citep*[e.g.,][]{gonzalez-alfonso_1997,pardo_1998}.

Two models can reasonably explain the observations:
\begin{enumerate}
\item Gas rotating around the axis described by the extended continuum emission (Northeast-Southwest direction; P.A.~$\simeq50^\circ$).
  Part of the gas responsible of the blue-shifted emission to the Southeast would leave the main stream during its spinning movement toward the Northwest due to instabilities or interactions with the recently ejected matter, falling onto the star with a Doppler velocity relative to the systemic one up to $\simeq10-15$~km~s$^{-1}${} (Fig.~\ref{fig:f3}, lower insets).
  The rest of the gas would keep revolving around the star at an average velocity of $\simeq6-8$~km~s$^{-1}$.
  A similar Doppler shift is expected for gas orbiting a star with a mass of $\simeq 0.5-1.0~M_\odot$ if an orbital inclination of $\simeq 25^\circ-50^\circ$ is assumed.
\item The previous rotating structure would exist as well but there would not be any leakage resulting in gas falling onto the star.
  The red-shifted absorptions would be produced by radial movements of photospheric layers (or convective cells) due to the stellar pulsation \citep{hinkle_1978,hinkle_1979a,hinkle_1979b}.
\end{enumerate}

Contrary to what \citet{vlemmings_2017} claimed to happen in the photosphere of W~Hya, where there is a bright spot suggested to be caused by the collision of infalling matter onto the star, we do not find any footprint of this phenomenon in our continuum maps.
However, the angular resolution of our observations is insufficient to explore the stellar photosphere and a similar bright spot could be highly diluted.
In fact, \citet{tatebe_2008} found evidences of a hot spot in the southern part of R~Leo that might have a counterpart in the mm range.

It can be argued that the Doppler velocities seen in our observations are related to random gas movements in the stellar photosphere or to effects of periodic shocks.
Several works based on VLBI observations of SiO masers have been published for a number of AGB stars such as R~Leo or TX Cam \citep*[e.g.,][]{cotton_2004,soria-ruiz_2007,gonidakis_2013}.
The maser spots are supposed to live enough time to trace the gas movement very close the photosphere, if a good sampling is adopted.
Those observations are constrained to the shell where the maser emission is higher, which covers the region ranging from 0\farcs030 to 0\farcs040.
They thus sample a region of the envelope considerably closer to the star than the one traced by our data.
Common conclusions are therefore difficult to draw.
Contrary to what happens with VLBI observations, the angular resolution of our observations is insufficient to describe random gas movements, which are cancelled out during the averaging performed by the PSF and contribute to broaden line widths.

A set of several clumps ejected with significant axial velocities resulting, for instance, from ballistic trajectories, could also be suggested to comprise the detected structure.
However, our observations show a highly continuous velocity field around the star difficult to explain assuming randomly ejected, discrete clumps with unconnected kinematical properties.
Additionally, there exist evidences of rotating molecular torus in the envelopes of other AGB stars such as L$_2$ Pup or R Dor \citep{kervella_2016,homan_2018}, suggesting that these structures might be more common than initially expected.

The existence of rotating gas in the surroundings of the star needs an angular momentum input, since the rotation of an AGB star is expected to be very low.
This extra angular momentum could be provided by a dipolar magnetic field \citep{matt_2000} or an additional body \citep{vlemmings_2018}.
The first scenario would require a plasma or a gas of charged particles to be coupled with.
No continuum emission typical of a plasma has been detected but dust grains are probably charged in the innermost envelope.
A dipolar magnetic field could force the radially ejected grains to drag the gas inducing circumstellar rotation.
Regarding the second scenario, no companion has been found in the optical and near-IR spectral ranges \citep{gatewood_1992} but this disagreement could be solved if it is a small, cold body such as a low mass star or a Jovian planet instead.
This idea was already invoked by \citet{wiesemeyer_2009b} to explain the polarisation fluctuations detected in the SiO maser emission \citep{wiesemeyer_2009a}, that are also formed very close to the stellar photosphere.
The existence of a planet orbiting an AGB star has been also suggested to be possible by \citet{kervella_2016}, who claim to have discovered what could be one of them in L$_2$ Pup.

\section{Summary and Conclusions}
\label{sec:summary}

We present new interferometer observations of the AGB star R~Leo carried out with ALMA with a maximum angular resolution of $\simeq0\farcs026$ and a MRS~$\simeq0\farcs6$.

The continuum emission can be described as a compact source surrounded by a highly collimated extended emission in the Northeast-Southwest direction.
The CO maps show a 0\farcs1-sized rotating torus-like structure with its axis lying along the Northeast-Southwest direction.
This structure could be explained by (1) either the action of a dipolar magnetic field on the radially expanding wind or (2) the presence of a small, cold companion (maybe a Jovian planet) that is orbiting the central star and injecting angular momentum to the gas.
This seemingly rotating structure is accompanied by evident red-shifted absorptions that could be produced by infalling gas coming from it or by the regular movement of the pulsating photosphere.
At larger scales $^{29}$SiO shows also a torus-like structure with its axis roughly oriented along the North-South direction unseen in the corresponding CO map.
Its origin is unclear but it may be produced by anisotropic physical and chemical conditions, or unexpected excitation effects.

These observations reveal a high complexity structure in the envelope of R~Leo tightly related to the changes undergone by this AGB star.
They will be addressed accurately in the near future through detailed analyses of these and new larger scale data.

\begin{acknowledgements}
We thank the European Research Council (ERC Grant 610256: NANOCOSMOS) and the Spanish MINECO/MICINN for funding support through grant AYA2016-75066-C-1-P and the ASTROMOL Consolider project CSD2009-00038.
MA thanks funding support from the Ram\'on y Cajal program of Spanish MINECO (RyC-2014-16277).
We also thank V. Bujarrabal for his valuable comments about the observations interpretation and the anonymous referee for his/her interesting comments.
\end{acknowledgements}