%~Mouliné par MaN_auto v.0.27.3 2023-01-06 16:51:34
\documentclass[CRPHYS, Unicode, XML, Thematic, screen, HideArticleType]{cedram}
\makeatletter\def\cdr@extrapackages{upgreek}\makeatother

%\TopicEN{Exoplanets}
%\TopicFR{Exoplanètes}

\ItHasTeXPublished
%\usepackage{amsmath}
\usepackage[noadjust]{cite}
\usepackage{upgreek}
%Le caracture\upbeta ne fonctionne pas sur tralics et a dû être retiré dans la biblio pour beta...%%
%\hypersetup{ colorlinks=true, breaklinks=true, urlcolor= blue, linkcolor= black, citecolor=blue, }
%\newcommand{\upbeta}{\beta}

\def\thanksname{Note}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\graphicspath{{./figures/}}

\newcommand*{\mk}{\mkern -1mu}
\newcommand*{\Mk}{\mkern -2mu}
\newcommand*{\mK}{\mkern 1mu}
\newcommand*{\MK}{\mkern 2mu}

%\hypersetup{urlcolor=purple, linkcolor=blue, citecolor=red}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\RubricEF{Exoplanets}{Exoplanètes}


\title{Observations of circumstellar disks in scattered light with SPHERE at the VLT}
\alttitle{Observations des disques circumstellaires en lumière diffusée avec SPHERE au VLT}

\author{\firstname{Anthony} \lastname{Boccaletti}}
\address{LESIA, Observatoire de Paris, Universit{\'e} PSL, CNRS, Sorbonne Universit{\'e}, Sorbonne Paris Cit{\'e}, 5 place Jules Janssen, 92195 Meudon, France}
\email[A. Boccalettit]{anthony.boccaletti@obspm.fr}

\thanks{Follows up on a conference-debate of the French Academy of Sciences entitled ``Exoplanets: the new challenges'' held on 18 May 2021, visible via\\
\url{https://www.academie-sciences.fr/fr/Colloques-conferences-et-debats/exoplanetes.html}}

\thanks{\selectlanguage{french}\textbf{Note.} Fait suite à une conférence-débat de l’Académie des sciences intitulée \og~Exoplanètes~: les nouveaux défis~\fg{} tenue le 18 mai 2021, visible via\\
\url{https://www.academie-sciences.fr/fr/Colloques-conferences-et-debats/exoplanetes.html}}

\keywords{Circumstellar disks, planet formation, exoplanets, high contrast imaging, polarimetric imaging}
\altkeywords{Disques circumstellaires, formation des planètes, exoplanètes, imagerie haut contraste, imagerie polarimétrique}



%\subjclass{00X99}

\begin{abstract}
The technological developments initiated in the early 21$^{\rm st}$ century have led to the implementation of ``planet finders'' instruments on 8-m class telescopes which are in operation since 2014-2015. Such facilities are at the inception of significant progresses in the study of exoplanetary systems by enabling high contrast imaging of the environment of young nearby stars. One of the most productive science in this field is certainly the observations at high angular resolution of circumstellar disks, from the very young gas-rich protoplanetary disks all the way to more elusive dust-rich debris disks. In this paper, we summarize the main results obtained for these two categories of objects focusing on their morphological characteristics, as well as on the measurement of optical properties of the small dust particles to which the near IR spectral range is sensitive to. We mostly provide exemples based on observations performed with the SPHERE instrument at the VLT.
\end{abstract}


\begin{altabstract}
Les développements technologiques initiés au début du 21\textsuperscript{ème} siècle ont conduit à la mise en place d'instruments ``découvreurs de planètes'' sur des télescopes de la classe 8m en fonctionnement depuis 2014-2015. De telles installations sont à l'origine de progrès significatifs dans l'étude des systèmes exoplanétaires en permettant l'imagerie à haut contraste de l'environnement des jeunes étoiles proches. L'une des sciences les plus productives dans ce domaine est certainement l'observation à haute résolution angulaire des disques circumstellaires, depuis les très jeunes disques protoplanétaires riches en gaz jusqu'aux disques de débris plus insaisissables, riches en poussières. Dans cet article, nous résumons les principaux résultats obtenus pour ces deux catégories d'objets en nous concentrant sur leurs caractéristiques morphologiques, ainsi que sur la mesure des propriétés optiques des petites particules de poussière auxquelles la gamme spectrale du proche IR est sensible. Nous présentons surtout des exemples basés sur des observations réalisées avec l'instrument SPHERE au VLT.
\end{altabstract}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\dateposted{2023-03-21}
\begin{document}
\maketitle



%\selectlanguage{french}

\section{Introduction}
Understanding how and when planets form is a central question in the study of exoplanets. In the current paradigm, mostly based on Solar System observations, a stellar system emerges from the gravitational collapse of a gaseous cloud, which leaves a protostar in the center and a disk of gas and dust orbiting the star in a nearly flat distribution, owing to the conservation of angular momentum. Dust will grow from sub micron-size to form pebbles~\cite{Ida2016}, which after a series of accretion and collision interplays can produce km-size bodies, i.e. the planetesimals~\cite{Kokubo2000, Armitage2016}. Oligarchic growth~\cite{Kokubo1998} can explain how these planetesimals will end up in massive planetary cores (likely a few $M_{\oplus}$), and finally, these cores will accrete gas~\cite{Pollack1996} resulting in the formation of giant planets. This whole process of “core accretion” is estimated to last about 1\,Myr or so, after which the gas is dissipated by photo-evaporation~\cite{Hollenbach1994}, or integrated into the atmospheres of giant planets. This scenario (see details in~\cite{Raymond2020}) provides the big picture, but one lesson from exoplanets searches is clearly the diversity of their properties. For instance, some alternative models to the core accretion scenario have been proposed, like the gravitational instability~\cite{Boss1997}, which is able to generate giant planets more rapidly and directly from the collapse of the gas in the protoplanetary disk, in a sort of small-scale stellar formation mechanism. It is not unlikely that several mechanisms take place among stellar systems or even in a same system.

Two means can be considered to address the context of planet formation from the observation of the exoplanet population. First, we can study statistically the orbital properties (semi major axis, inclination, eccentricity), and the bulk properties (mass, radius, hence density) of known planets, and compare it to synthetic planet population models~\cite{Mordasini2009}. Another way consists in measuring their composition to derive chemical constraints on the abundance of elements in comparison with planet formation models~\cite{Oberg2011}. These two aspects require to characterize beforehand the environment in which planets form and evolve to provide inputs to these planet formation models.

In the early stage, a system contains large amount of gas and dust (with a typical 100:1 ratio) in the form of a protoplanetary disk, also known as a ``planet forming disk'', which can determine locally the composition of a planet, or even its fate given that planets interact with the gas, and migrate inward or outward~\cite{Kley2012}. Mature systems, which represent the vast majority of known exoplanets, have experienced significant dynamical evolution, so that the determination of the conditions of their formation is not straightforward. For this reason, observing the systems when they are just a few million year old is more relevant to put constraints on planetary formation. Protoplanetary disks have complex morphologies as a result of gas dynamics and size-dependent coupling with dust particles, eventually leading to planet formation. In a later stage, when the system evolves and the gas dissipates, planetesimals become the most abundant objects. Collisional cascades among the planetesimals produce second generation micron-sized grains, hence the term of ``debris disk''~\cite{Wyatt2008}. Like in the Solar System the planetesimals are distributed in the form of belts, similar to the Kuiper and the asteroid belts, and located at a few tens of astronomical units (AU) from the star.

Historically, the first disks were identified from the excess emission that the dust produces at infrared wavelengths, like the so-called Vega phenomenon~\cite{Aumann1984, Backman1993, Artymowicz1997}, the emission peak being proportional to the temperature of the grains, and thus to the stellocentric distance of the dust, together with its optical properties. The IRAS satellite has been instrumental in detecting such signatures (in excess compared to the blackbody-like emission of a single star, see Fig.~\ref{fig:sed}) at IR wavelengths of 10 to 100\,$\mu$m typically, and interpreted as either protoplanetary disks (near/mid IR emission, high fractional IR luminosity) or debris disks (far IR emission, low fractional IR luminosity). However, the spectral energy distribution (SED) of a system derived from photometric measurements carries important degeneracies related to its geometry (single belt, multiple belts, continuous disk with or without cavity) and to the dust properties (emission depends on temperature, size, shape, and composition of grains) \cite{Wyatt2008}.


\begin{figure*}[htbp]
\begin{center}
\includegraphics[trim={0cm 8cm 0cm 6cm},clip, height=6cm]{figures/SED_bPic_Malfait98.pdf}
\includegraphics[trim={2cm 8cm 2cm 7cm},clip, height=5cm]{figures/bPic_SmithTerrile87.pdf}
\end{center}
\caption{Spectral Energy Distribution of $\beta$ Pictoris measured with IRAS (left), and first image of the debris disk observed with a 2.5-m telescope at 850\,nm and with a 7$''$ diameter coronagraph~\cite{Smith1984} (right). The FoV is about 92$''$x82$''$.}\label{fig:sed}
\end{figure*}


Given the aforementioned characteristics, and because exoplanetary systems have AU-scale dimensions, imaging disks with high angular resolution and high sensitivity has been decisive to characterize them further, both in scattered light at visible and near IR wavelengths, as well as in thermal emission at millimetric wavelengths. The first image of a debris disk was obtained in the early 80s from the ground, around the famous star $\beta$ Pic, featuring a favorable combination of proximity, size, luminosity and inclination~\cite{Smith1984}. It took more than a decade to achieve similar level of performance on other systems, thanks to the HST~\cite{Weinberger1999, Augereau1999, Augereau2001, Schneider1999}, while on the ground the sensitivity has been a strong limitation~\cite{Boccaletti2003} until the implementation of adaptive optics on 8-m class telescopes. Taking advantage of an angular resolution of 40-50 mas corresponding to a few AU for the nearest stars like $\beta$ Pic, to a few tens of AU in the main young star associations at 100-150\,pc, the exploration of the brightest disks unveiled a rich morphology with spirals, gaps, and rings~\cite{Clampin2003, Quanz2011, Quanz2013, Boccaletti2009, Kalas2005, Kalas2006, Kalas2007, Krist2005, Schneider2014}. Unsurprisingly, this breakthrough in disk imaging coincides with the first obtained images of extrasolar giant planets~\cite{Chauvin2004, Lagrange2009, Lagrange2010, Marois2008, Marois2010, Kalas2008} (see~\cite{Chauvin2022} this proceeding for a review of exoplanet imaging results and projects).

A revolution came with extreme adaptive optics facilities (so called ``planet finders'') combined to optimized coronagraphs, and also because significant efforts were invested in the observational and data reduction strategies~\cite{Cantalloube2021}. The instruments SPHERE at the VLT~\cite{Beuzit2019}, GPI at Gemini~\cite{Macintosh2014}, and SCExAO at Subaru~\cite{Lozi2018} were prolific to resolve structures in protoplanetary and debris disks in scattered light, shedding light on how the matter is distributed in the first stages of planetary formation. The morphology and the fine structures distribution of millimeter dust and gas in protoplanetary disks is also revealed owing to radio-interferometry with ALMA which provides complementary informations about those exosystems with a similar angular resolution as SPHERE~\cite{Andrews2018, ALMApartnership2015}.

This paper focuses exclusively on the measured properties of circumstellar disks recently observed in scattered light, based in particular on examples from SPHERE observing programs, although not exhaustively.


\section{Detection methods}

Circumstellar disks are intrinsically angularly small, and faint, featuring a large contrast (>$10^{-3}$, defined as the ratio between the intensity of an object in the field and the host star) with respect to the central star, hence require large telescope diameters, high sensitivity and high dynamic to be observed. Therefore, the techniques for disk imaging and exoplanet imaging are identical, both being off-axis, resolved or unresolved, objects with respect to the optical axis defined by the star (assumed unresolved). In that respect, an extreme adaptive optics (AO) system like SPHERE is able to routinely achieve Strehl ratio larger than 80\% in the near IR (J, H and Ks bands) for bright stars (R<10) and can even provide a Strehl of 50\% at visible wavelengths for R<9~\cite{Sauvage2016}, allowing to deliver exquisite image quality with angular resolution of 15 to 50\,mas for wavelengths in the range 0.6-2.0\,$\mu$m. Such a low wavefront error level is necessary for a coronagraph to provide a large rejection of the starlight (see~\cite{GalicherMazoyer2022} this proceeding for a review of coronagraphic techniques). The SPHERE coronagraph~\cite{Carbillet2011} has a diameter of about 0.2$''$ which means that no off-axis object can be recovered below the Inner Working Angle of $\sim0.1''$. Raw contrasts measured with SPHERE are typically reaching $10^{-3}-10^{-4}$ in separation range $0.1-0.3''$, and better than $10^{-4}$ between $0.3''$ and the correction radius at $\sim0.8''$ (see~\cite[Fig.~10]{Beuzit2019}).

To reach a contrast that is compatible with exoplanet or disk imaging, the residual starlight left behind the coronagraph must be further removed or at least attenuated. This starlight pattern is a complex mixture of diffraction effects which can be static, highly variable, or quasi-static, the latter being the more difficult to get rid of. In conjunction with the improvement of AO systems, several techniques, combining observing strategies and numerical methods, were developed with the purpose to calibrate the stellar light, independently of the off-axis object to detect. For disk imaging, this can be accomplished in three different ways, taking advantage of data diversity introduced by a specific observing mode in order to disentangle the starlight from astrophysical sources around the star, and taking into account the properties of these sources.

The most obvious technique relies on the observation of one or several reference stars observed in the same configuration as the science target, usually in a short time frame to ensure a somehow similar starlight pattern. The so-called reference star differential imaging (RDI) has been very successful with the HST (Fig.~\ref{fig:Diff_Imaging}), in an atmosphere-free environment where the global wavefront is highly stable. In fact, following the discovery of the $\beta$ Pictoris debris disks, the first scattered light disk images at high angular resolution, were obtained with RDI~\cite{Weinberger1999, Augereau1999}. The method was generalized later to using a set of reference stars, with for instance linear combinations of frames and Principal Component Analysis, which allowed to discover a few more disks left undetected in archival HST data~\cite{Choquet2014,Choquet2016}. On the ground, because of the atmospheric turbulence, but also due to the inevitable gravity related variations of the optical wavefront, RDI has been much less effective until recently~\cite{Xie2022}. However, advanced numerical techniques, like artificial intelligence, are now being deployed to tackle this challenge since the amount of data, hence potentially of reference frames, is just gigantic.

A second technique, devised for exoplanet detection, is taking advantage of the field rotation at the focus arising from the alt-azimuthal mount configuration in large telescopes or rolls with space telescopes (Fig.~\ref{fig:Diff_Imaging}). If the wavefront errors are mostly coming from a pupil plane the net effect is that an off-axis object rotates in the image as the parallactic angle, while most of the starlight is fixed, or evolves differently (with the exception of the aberrations caused by high altitude winds). As such, Angular Differential Imaging (\cite[ADI]{Marois2006}) has been used for detecting all the directly-imaged exoplanets known today and comes with various flavors of data processing techniques~\cite{Lafreniere2007, Soummer2012}. An important aspect of ADI is that it implies a self-subtraction of the object to detect, as long as this object contributes, even in a small fraction, to the construction of the reference frame(s) to be subtracted to the data set. For a planet, which can be modeled as a point source, the self-subtraction can be taken into account to measure its photometry, but the problem is far more complex for an extended structure like a disk. In that case, the self-subtraction depends on the morphology of the object which is basically what we are trying to measure~\cite{Milli2012}. As a result, we have to assume a particular morphology and perform forward modeling to recover the disk geometric parameters and then its optical properties~\cite{Bhowmik2019}. While this is doable for debris disks which can be approximated as a belt (or multiple belts) and providing the inclination is relatively large, it is totally unpractical for protoplanetary disks. In practice, we are only measuring departures to an azimuthally-symmetrical geometry. To get around this problem some algorithms have been recently proposed to attempt preserving at most the flux of an extended object~\cite{Flasseur2021, Pairet2021}. Still, performing direct measurements of complex structures around stars with ADI remains challenging and could lead to wrong interpretations of the morphology: streamers can be seen at the edge of a highly inclined debris disk belt as an amplification in one direction of the disk halo~\cite{Thalmann2011}, or rings and spirals can be pinched in some directions as a result of the disk minor axis rotating slower than the outer part combined to the effect of the scattering phase function~\cite{Garufi2016}. Both ADI and RDI can be used with either IRDIS the camera of SPHERE and IFS the low resolution spectro-imager, providing spectral information from 0.95 to $2.3\,\mu$m.

Finally, the third technique benefits from the polarization of the light induced by the scattering of the starlight by dust grains, while the starlight itself is generally unpolarized. Therefore, filtering the linearly polarized signal in the image of an exosystem is highly efficient to reduce the contribution from the star and to reveal the dust contribution. Dual Polarization Imaging (DPI, also referred as Polarimetric Differential Imaging) is implemented in SPHERE IRDIS but is only available in broad band or narrow band filters. Unlike ADI, the morphology of the disk is fully preserved, but as a drawback the beam splitter and the polarizers implemented in SPHERE transmits only 50\% of the light, while the polarized signal of a disk is also intrinsically smaller (usually a few tens of \% of the unpolarized flux). DPI comes also with a series of observational constraints but combining DPI and pupil tracking can circumvent these issues, and provides simultaneously for the same data set the polarized and unpolarized image~\cite{Holstein2020, Boer2020}. Yet, DPI imaging has been extremely efficient for resolving structures in the vast majority of protoplanetary disks observed with SPHERE~\cite{Garufi2017}.

\begin{figure*}[htbp]
\begin{center}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=7cm]{figures/HD141_RDI_ACS.jpg}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=7cm]{figures/bpic_ADI_ACS.jpg}
\end{center}
\caption{Illustrations of the RDI (left) and ADI (right) methods from the same instrument: ACS on HST. The ACS PSF being strongly aberrated, the coronagraphic images that are displayed here are very extended. The left panel shows an observation of HD\,141569 and a reference star together with the subtraction revealing the gas-rich debris disk~\cite{Clampin2003}. The right panel is an observation of the disk of $\beta$ Pic observed at two roll angles, demonstrating the stability of the coronagraphic diffraction pattern with respect to the disk itself which rotates with the field~\cite{Golimowski2006}.}\label{fig:Diff_Imaging}
\end{figure*}



\section{Characterization of protoplanetary disks}

Protoplanetary disk is a generic term which corresponds to the phase between the formation of the star and the transition to the debris disk phase when the gas is dissipated. This period probably lasts a few Myr during which the system can undergo significant evolution which impacts the disk structuring. Unresolved photometry in the near to far IR has provided in the past some hints about the structures in these disks, in particular the presence of cavity or gaps, which shows up as features in the SED of these systems owing to the temperature-distance dependence. Imaging is able to retrieve much more details in the morphology, and allows to leave some degeneracies when modeling the SED. All disks observed in scattered light have structures. If they don't, it is essentially because they are too small and/or too faint given the sensitivity~\cite{Garufi2016}. Structures are easier to detect in polarimetry if the disk is nearly face-on, while at high inclinations the effect due to the phase function or the spatial distribution of the dust can be degenerated. In the case of protoplanetary disks, combining observations in the near IR and millimeter wavelengths is crucial to test formation theories because the two regimes are sensitive to different sizes of dust particles (sub-micron, and millimeter respectively), the distribution of which can differ during the evolution of a system. For instance, small particles are efficiently coupled to the gas while the largest are not. Therefore, we can observe different distributions and dynamics with SPHERE and ALMA. Below, we discuss the main structures that we manage to observe in scattered light.


\subsection{Spirals}
Spirals can be the result of density waves being launched by a planet orbiting the star and propagating in the gas~\cite{Dong2015a}. Near IR observations in scattered light are sensitive to micron-size grains which are strongly coupled to the gas, so these images are effectively tracing the perturbation in the gas distribution. One of the first protoplanetary disk observed with SPHERE, MWC\,758~\cite{Benisty2015}, has two prominent and almost symmetrical spirals which totally dominate the scattered light image. These patterns can be interpreted using the linear theory of wave density~\cite{Rafikov2002} as being induced by a gravitational perturber in Keplerian rotation around the star, so presumably a planet. The linear relation to describe the shape of a spiral in polar coordinates $\theta(r)$ ($r$ being the radial distance to the star), as a function of the planet location ($r_c$, $\theta_0$) and disk properties reads as follows~\cite{Muto2012}:
\begin{equation}\label{eq:muto}
\theta(r) = \theta_0-\frac{\text{sgn}(r-r_c)}{h_c} \times \left[\left(\frac{r}{r_c}\right)^{1+\beta} \left\{ \frac{1}{1+\beta}-\frac{1}{1-\alpha+\beta} \left(\frac{r}{r_c}\right)^{-\alpha} \right\} \right. \left. -\left(\frac{1}{1+\beta} - \frac{1}{1-\alpha+\beta} \right)\right]
\end{equation}


The linear theory is convenient in the context of direct imaging as it provides an analytical expression of the spiral shape being a function on the planet location, the angular velocity and the disk aspect ratio (scale height divided by the radius), hence related to the gas temperature. The properties of the disk are contained in the parameter $\alpha$ and $\beta$, which are the exponents of the power law describing the variation of the angular velocity $\Omega \propto r^{-\alpha}$, and the sound speed $c_s \propto r^{-\beta}$. The fit of this relationship in MWC\,758 (Fig.~\ref{fig:diskproto}), but also HD\,135344B~\cite{Stolker2016} images, leads to several solutions with either inner or outer planets. The case with inner planets yields large aspect ratios corresponding to unphysical temperatures, so that the scenario with outer planets is usually favored. The same conclusion is reached based on hydrodynamical models in which a pair of a so-called ``grand-design'' spirals can be generated by a single perturber~\cite{Dong2015a}. In complement, direct imaging in unpolarized light with SPHERE allows to put strong constraints on the presence of giant planets (undetected in polarized light) in these systems with more sensitivity at large separations, precisely where the presence of planets is favored by the analysis of spiral arms. However, no planet was discovered yet in these systems. For instance, the SPHERE data achieve detection limits as low as 1$\sim$2 M$_J$ for HD\,135344B~\cite{Maire2017} and MWC\,758~\cite{Boccaletti2021}.

As of today we never identified a bona-fide planet related to a spiral arm, with the exception of HD\,100453~\cite{Wagner2015, Dong2016}, in which the perturber is not a planet but a low-mass star. In the context of findings connection between spirals and planets, the even more complex system AB\,Aur is unique to understand this connection, because it is one out of the two cases where a pair of spirals is detected both in scattered light and in CO gas with ALMA~\cite{Tang2017}. Two independent observations reported hints of the presence of two protoplanets possibly associated to the spiral arms. First, with SPHERE in polarimetry, the eastern spiral features a twist which matches relatively well the classical spiral model, this twist being the place where the inward and the outward spiral components are launched~\cite{Boccaletti2020} (Fig.~\ref{fig:structuresproto}). Second, based on SCExAO data, a diffuse region at the tip of the western spiral has been modeled as a combination of a protoplanet emission and scattering from the disk~\cite{Currie2022}, the location of which is also in line with ALMA observations~\cite{Tang2017}. The disk around HD\,169142 was also scrutinized with SPHERE with very aggressive processing to remove the starlight which allowed to detect or re-detect several blobs inside the complex spiral pattern of this disk. Since these disk features also move at Keplerian velocity~\cite{Gratton2019} the connection to planets is still to be demonstrated.

Dynamically, being waves, spirals are expected to rotate at the same angular velocity as the source, and as such can appear to move at a super-Keplerian speed outwards of the launching site, or sub-Keplerian inwards. However, it is challenging to obtain a precise measurement, especially for outer planets simply because the orbital period can be as large as a few hundred years, while an inner planet may produce a visible motion in a few years. We note that in the case of MWC\,758 the extreme reproducibility of SPHERE images and the knowledge on the field orientation allowed to derive a rotation of $0.22^{\circ}$ consistent with an outer perturber at 172\,AU~\cite{Ren2020} still to be discovered. Follow up of other systems, in particular AB\,Aur may allow to better constrain the origin of spirals.

Most of the disk with spirals were observed face-on because it is the most efficient case in DPI with SPHERE, but in a few cases, based on a comparison to radiative transfer models, we were able to identify spirals in more inclined disks, like for RY~Lup~\cite{Langlois2018}. Interestingly, for RY~Lup and also T~Cha~\cite{Pohl2017}, the direction of polarization is not fully azimuthal, indicating that multiple scattering occurs within the dusty disk.

Overall, the connection between spirals and planets is not fully confirmed. Spirals can also arise from gravitational instabilities which can develop when the disk is self-gravitating~\cite{Dong2015b}, under the conditions that the disk is more massive than a quarter of the star's mass. The aforementioned systems are not expected to be that massive. However, some disks like AB\,Aur (Fig.~\ref{fig:structuresproto}) and HD\,100546 show a variety of spirals morphologies which may have a different origin, possibly self-gravity~\cite{Boccaletti2020, Garufi2016}.

\begin{figure*}[htbp]
\begin{center}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, width=\linewidth]{figures/montage_disk_proto.jpg}
\end{center}
\caption{Images obtained with SPHERE illustrating main features observed in protoplanetary disks in scattered light (from left to right): the spiral arm of MWC\,758~\cite{Benisty2015}, the rings/gaps in TW\,Hya~\cite{vanBoekel2017}, and the shadows in HD\,135344\,B~\cite{Stolker2016}.}\label{fig:diskproto}
\end{figure*}


\subsection{Rings and cavities}
Cavities and rings in systems featuring IR excess were very soon suspected from the dip in the excess at near IR wavelengths. The dust depletion at short stellocentric distances can be the result of photo-evaporation, dust evolution, or gravitational sweeping by massive bodies like planets, the latter being frequently involved as an argument to target such systems in high-contrast imaging surveys~\cite{Meshkat2015}. Cavities are found in most protoplanetary and transition disks observed in scattered light with SPHERE. Like spirals, the rings and cavities are easier to detect at moderate disk inclinations with either DPI and/or ADI, as showcase with the images in both polarized and unpolarized intensity of HD\,97048~\cite{Ginski2016} and RXJ\,1615~\cite{deBoer2018}. Their size can be different than in the millimeter regime indicative of dust segregation which may be explained by dust trapping by planets at the gap edges. With sizes in the range $\sim200-400$\,AU, and inclinations of $40-50^{\circ}$, the fitting of ellipses to the multiple rings detected in these systems allows to directly measure the scale height which can be as large as 15 to 20\,AU (measured at 100\,AU). Such large values are related to the fact that scattered light is sensitive to the surface of the disk (corresponding to optical depth $\tau$ = 1). In nearly face-on geometries, alternation of rings and gaps was also evidenced owing to differential polarimetry in the disk around TW\,Hya~\cite{vanBoekel2017} or HD\,169142.~\cite{Ligi2017}. In that cases, surface density can be inferred from surface brightness radial variations. The gaps in TW\,Hya (Fig.~\ref{fig:diskproto}) are found to be relatively shallow which would imply a few to tens Earth-masses perturbers if they have a gravitational origin. Since such planets are not yet detected, other scenarios involving the accumulation of dust particles at snowlines are still plausible~\cite{Zhang2015}.

There are a few examples for which a cavity can be directly connected to one or more perturbers, most are low mass stars like for HD\,142527~\cite{Claudi2018} and GG Tau~\cite{Keppler2020} (Fig.~\ref{fig:structuresproto}). In the latter, the remarkable precision delivered by SPHERE allowed to interpret the many fine structures at the edges of the cavity with hydrodynamical models as being triggered by periodic perturbations of the central binary star. The case of PDS\,70~\cite{Keppler2018, Muller2018, Haffert2019, Mesa2019} is certainly the most relevant in this category since at least two giant planets have been identified in the cavity orbiting respectively at 22\,AU and 35\,AU while the cavity is 55\,AU wide. These very young objects still in a formation phase~\cite{Benisty2021} are presumably sculpting the disk.

\subsection{Shadows} 
Direct imaging and in particular DPI has also revealed shadow patterns in some of the brightest protoplanetary disks. Shadows are systematically associated to large near IR excess which is in line with the fact that they take origins in the inner part of these systems due to inner belts casting shadows in the outer parts. One of the first example was found in the disk of HD\,142527 which features a very large cavity and two localized centro-symmetrical shadows in the outer belt at the northern and the southern tips. They are observed both in the visible~\cite{Hunziker2021} and near IR~\cite{Avenhaus2014} and explained using radiative transfer with an optically thick disk close to the star and warped by 70$^{\circ}$ with respect to the outer one~\cite{Marino2015}. Similar conclusions were reached for the shadows of HD\,100453~\cite{Benisty2017} with the particularity that these shadows are apparently connected to two spirals. Although these spirals are likely gravitationally triggered by an outer M star companion, the shadows may also play a particular role here. Finally, shadows can have much more complex structures than just two anti-symmetrical dark regions casted on an outer belt, as in SAO\,206906 where several shadows feature azimuthal variations at different epochs (Fig.~\ref{fig:diskproto}), reminiscent of material in rotation at much closer orbital radius~\cite{Stolker2016}.


\begin{figure*}[htbp]
\begin{center}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=4.5cm]{figures/HD169142_blobs.jpg}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=4.5cm]{figures/GGTau.jpg}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=4.5cm]{figures/ABAur_twist.jpg}
\end{center}
\caption{Typical small-scale structures: blobs in HD\,169142~\cite{Gratton2019}, cavity and multiple spirals in GG\,Tau~\cite{Keppler2020}, spiral twist in AB\,Aur~\cite{Boccaletti2020}.}\label{fig:structuresproto}
\end{figure*}

\section{Characterization of debris disks}
The debris disk phase is referring to the next stage in the evolution (the typical age is 10$\sim$15 Myr) of planetary systems when most of the gas is dissipated and large bodies are formed. The system becomes dominated by gravity so that collisions among planetesimals produce ``second-generation'' dust particles down to micron sizes. As a result, the dust is mostly generated in one or several belts where the planetesimals collide. Because debris disks are intrinsically faint, polarimetric imaging is less efficient than for protoplanetary disks but remains the only option for disks seen at low inclination. On the contrary, non polarimetric imaging is also used more frequently as long as the disk morphology is easier to describe with parametric models and the use of forward modelling. This observing mode strongly bias the detection towards highly inclined disks. High contrast imaging can essentially constrain their morphology (possibly with temporal resolution) and the optical properties of dust particles. Hereafter, we summarize the main physical characteristics of debris disks that are derived from observations.

\subsection{Morphology and interaction with planets}
For most known cases, debris disks are dust belts typically located at a few tens of AU from the star. The disk around HR\,4796 as observed with SPHERE~\cite{Milli2017} gives a clear example of such an evolved system. It is relatively straightforward to estimate its size, its inclination (the angle with respect to a face-on orientation) and its position angle (rotation in the plane of the sky with respect to the North). As such, direct imaging of debris disks provides unambiguous measurements of the locations of planetesimals, which once compared to previous estimations made from SEDs, lead to large discrepancies, as a clear proof that the latter is subject to degeneracies between temperature and dust location. Belts in debris disks are usually rather sharp, especially at their inner boundary which calls for the presence of a perturber shepherding the edge of the disk due to its chaotic zone~\cite{Wisdom1980}, although several planets would probably be needed to explain the large cavity observed in debris disks. No such connection between a planet and the sharp edge of a disk has been evidenced so far. On the contrary, the outer intensity slope of debris disk belts usually goes like $r^{-4}$ ($r$ being the physical separation to the star), which is the characteristics of dust particles blown away by the radiation pressure. A sharper slope may indicate truncation by an object outwards of the belt like for HR\,4796~\cite{Milli2017}.

Generally speaking, it is tantalizing to seek for any departure to this simple belt geometry since any non centro-symmetrical features could be the result of gravitational perturbations by planets, and could thus be used as an indirect signature of their presence. Indeed, the spatial structure and then the scattered light image of a debris disk can be strongly affected by the presence of planets as illustrated in some theoretical works~\cite{Lee2016}, depending on the disk inclination and planet's eccentricity. A case study is the system $\beta$ Pic, an edge-on debris disk, in which a warp has been identified in the first AO+coronagraph images~\cite{Mouillet1997}, and interpreted as the presence of a $\sim10\,M_J$ giant planet at $\sim10$\,AU on an inclined orbit (4-5$^{\circ}$). The planet was confirmed a decade later~\cite{Lagrange2010}, the properties of which are matching very well the predictions. While mostly all directly imaged planets are found in young systems, which also contain a dust component, only one other example, (with the exception of the much younger PDS\,70 system), features a planet and a debris disk in the same image: HD\,106906~\cite{Kalas2015}. However, this planet is found at a larger angular separation from the star (corresponding to 650\,AU in projection), but the fact that the disk belt is resolved allows to infer a much larger physical distance of 2000-3000\,AU~\cite{Lagrange2016} (if the planet and the disk were coplanar) suggesting an ejection mechanism.

An interesting outcome of SPHERE observations is the case of debris disks around HIP\,73145~\cite{Feldt2017}, NZ\,Lup~\cite{Boccaletti2019}, and HIP\,67497~\cite{Bonnefoy2017} which features at least two or multiple belts, together with the more complex HD\,141569 in which many rings, ringlets and broken structures have been detected~\cite{Perrot2016}. One plausible explanation is that such multiple-ring geometries are generated by instabilities rather than shaped by planets, for instance based on the so-called ``photo electric'' effect~\cite{Lyra2013}, which implies that the gas compound is not negligible in these system with gas to dust ratio close to unity. CO gas was detected at least in HIP\,73145 and HD\,141569, with the latter being categorized as a gas-rich debris disk and showing a very strong match between the gas and dust distribution~\cite{Difolco2020, Singh2021}. To explain this it has been shown that collisions within planetesimals in a belt can release both second generation dust and gas explaining their colocation~\cite{Kral2019}.

On the other hand we expect massive collisions between very large asteroids or small planets to occur in a young planetary system, leaving a reservoir of dust at the place of the release, which would be detected in scattered light images as brightness asymmetries~\cite{Jackson2014}. Such a mechanism has been posited to explain asymmetric features in several debris disks although without firm conclusions~\cite{Olofsson2020, Singh2021}.

\begin{figure*}[htbp]
\begin{center}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=5cm]{figures/HR4796_ZIMPOL.jpg}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=5cm]{figures/HD106906_Lagrange.jpg}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=4.7cm]{figures/HD141569_rings.jpg}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=4.7cm]{figures/HD32297_Itot_Ipol.jpg}
\end{center}
\caption{Exemple of debris images illustrating some of their morphological characteristics (from left to right, top to bottom): the planetesimals belt of HR\,4796 observed in visible polarimetry~\cite{Milli2019}, the ring and the planet of HD\,106906~\cite{Lagrange2016}, the multiple rings and ringlets in the gas-rich debris disk of HD\,141569~\cite{Perrot2016}, and the highly inclined disk of HD\,32297~\cite{Bhowmik2019} in total intensity and in polarimetry showing halo beyond the birth ring (dashed line).}\label{fig:diskdebris_phase}
\end{figure*}


\subsection{Scattering phase function}
In the generic case of a debris disk with a single belt seen at a given inclination, the intensity (relative to the star) along the belt varies with the scattering angle. This scattering phase function can be measured to derive optical properties of the dust grains as long as we assume the disk to be infinitely thin, with an azimuthally constant surface density, although we can expect departure from this simple hypothesis~\cite{Olofsson2020}. Both non polarimetric and polarimetric data are useful in this case in providing complementary informations, with the problematic that non polarimetric data are affected by self-subtraction modifying the disk intensity, so special care is needed in that case, for instance using forward modeling as long as the disk geometry can be described by a few parameters, or ideally, RDI processing. Because it is bright and seen at a favorable inclination, the disk around HR\,4796 has been extensively studied in this context~\cite{Milli2017,Milli2019}. The scattering phase function can be measured at small scattering angles (down to 13.6$^{\circ}$) which is decisive to test various dust properties (Fig.~\ref{fig:diskdebris_phase}). On the contrary, the polarized scattering phase function is rather flat at small scattering angles. The phase functions measured for HR\,4796 and other disks like GSC\,07396-00759~\cite{Adam2021}, HD\,32297~\cite{Bhowmik2019} and HD\,15115~\cite{Engler2019} show that Mie scattering of compact grains are not compatible with the observations, favoring particle aggregates as opposed to compact spheres. In the case of HR\,4796, to account for the amount of back scattering, the modeling of the scattering phase function suggests that grains are bigger than the observing wavelength. In a general case, it is found difficult to reproduce the scattering phase function with the standard Henyey--Greenstein function but it often requires multiple components. Equations~\ref{eq:hgItot} and~\ref{eq:hgpol} provide the analytical expression of these functions in total intensity and polarized intensity respectively, where $g$ is the aninostropic scattering factor, and $\theta$ is the scattering angle.
\begin{equation}\label{eq:hgItot}
f_{Itot}(\theta) = \frac{1 - g^2}{4\pi(1+g^2+2g\times \cos(\theta))^{3/2} }
\end{equation}

\begin{equation}\label{eq:hgpol}
f_{Ipol}(\theta) = f_{Itot}(\theta) \times \frac{1-\cos(\theta)}{1+\cos(\theta)}
\end{equation}

For illustration, a model of a disk image in scattered light, the SPF of which is defined by Eq. \ref{eq:hgItot}, is shown in Fig.~\ref{fig:modelspf}, for several inclinations and anisotropic scattering factors.

\begin{figure*}[htbp]
\begin{center}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=5cm]{figures/PhaseFunction_exemples.pdf}
\includegraphics[trim={0cm 1.5cm 0cm 2cm},clip, height=5cm]{figures/SPF_SPFP.pdf}
\end{center}
\caption{Left: disk models in total intensity for several $g$ values and inclinations (i=0$^{\circ}$, 30$^{\circ}$, 70$^{\circ}$ from top to bottom). Right: scattering phase functions and polarized scattering phase functions for several $g$ values.}\label{fig:modelspf}
\end{figure*}


Dust production is also a feature of comets in the Solar System which can be studied from a different perspective, sometimes in-situ, so there has been some attempts to compare the properties of dust in these two kind of objects to identify common optical properties hence to put constraints on the dust particles in debris disks in terms of size, shape and composition~\cite{Levasseur2020}.


\subsection{Spectroscopy}
SPHERE and other planet finders provide spectro-imaging capabilities at low resolution in the near IR extensively used for exoplanet characterization. In the case of disks, the spectral information has been mostly useful to process the data and reveal substructures~\cite{Gratton2019} since retrieving the spectral information in non polarimetric data involves a careful modeling to circumvent the problematic of self-subtraction, and necessarily suffer from a lower signal to noise ratio than in imaging. This has been doable in a very few cases like the bright disk of HD\,32297~\cite{Bhowmik2019} for which the negative slope (blue color) of the spectrum is compatible with sub-micron dust particles (Fig.~\ref{fig:diskdebris_phase}). This in line with the global picture of debris disks around A type stars in which small dust grains are put on very eccentric orbits or even blown out by the radiation pressure. Given the very high inclination of this disk, the spectrum is likely dominated by these grains that are seen along the line of sight everywhere in the disk image. On the contrary, HR\,4796 has a red spectrum indicative of bigger grains (micron size) in agreement with the size inferred from the scattered phase function~\cite{Milli2017} but at odd with respect to the SED. The flux measurements in scattered light either from photometry or from spectroscopy are usually difficult to reconcile with thermal emission derived from SED modeling for instance or from ALMA observations. More work towards combining multiwavelength observations is highly desirable.


\begin{figure*}[htbp]
\begin{center}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=5cm]{figures/HR4796_phasefunction.jpg}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=5cm]{figures/HD32297_spectrum.jpg}
\end{center}
\caption{Exemple of a scattering phase function (left) measured in total intensity in the disk of HR\,4796~\cite{Milli2017}, and of a near IR spectrum of HD\,32297~\cite{Bhowmik2019}.}\label{fig:diskdebris_phase}
\end{figure*}


\subsection{Variability and late type stars}
The exquisite resolution of SPHERE allows to obtain consistent measurements at several epochs to perform dynamical studies in debris disks. SPHERE has revealed an unexpected phenomenon in the one of a kind disk around the M-type star AU\,Mic which is young and ideally located at about 10\,pc. Large structures forming some arches a few AU wide were observed during SPHERE commissioning in 2014 on one side of this edge-on debris disk~\cite{Boccaletti2015}. Not only the structures were all recovered in older data from HST (2010, 2011) but they featured an apparent motion, the interpretation of which is complicated by the edge-on orientation. Follow up observations in the next years confirmed the evidence of the motion, and discovered other smaller structures in different places of the disk, all being moving over time~\cite{Boccaletti2018} (Fig.~\ref{fig:aumic}). Two characteristics were important to understand the behavior of the structures: first they are moving on straight trajectories (linear with time), and the most distant (in apparent projection) are even faster than the escape velocity. All together, some dust particles are blown away from the system, a particularity that was only expected for A-type star owing to the radiation pressure. For AU\,Mic, being a very active M-type star, the motion of the dust is very likely correlated with its strong stellar winds~\cite{Sezestre2017}. In addition, since we observe these structures in specific regions of the disk, it is not a uniform mechanism, so another object is necessarily involved to break the symmetry, possibly being a planet or a local reservoir of dust. Interestingly, the presence of at least two planets was confirmed with transit and radial velocity~\cite{Plavchan2020, Martioli2020}, although they are orbiting too close to the star to be involved in the production of the fast-moving dust features, because of timescale reasons. Whether the same phenomenon can be observed in other M-star systems is still to be determined.

\begin{figure*}[htbp]
\begin{center}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=5cm]{figures/AUMic_southeast_epochs.jpg}
\includegraphics[trim={0cm 0cm 0cm 0cm},clip, height=5cm]{figures/AUMic_model.jpg}
\end{center}
\caption{The edge-on disk of AU\,Mic observed at four epochs (left) showing the fast motion of five extended dust structures (SE1 to SE5). The star is at the right edge of each panel. The right panel displays a possible model with a planet at 34\,au~\cite{Boccaletti2018}.}\label{fig:aumic}
\vspace*{15pt}
\end{figure*}


\section{Conclusion}
Improving the angular resolution of instruments generally comes with higher contrasts. In that respect, SPHERE has been able to achieve significantly larger and more stable contrasts than the previous instruments generation. At the same time, some dedicated observing modes were specifically designed for disk science, allowing imaging in total intensity and polarimetry with a rather large field of view both in the visible and near IR, as well as spectro-imaging at low spectral resolution.

These unique capabilities have literally boosted the science on disks at very high angular resolution, and then also provided a mean to characterize in the very details the environment in which planet form. Emblematic structures like spirals, gaps, rings or shadows, were observed in planet forming disks with a precision never achieved before, leading to the identification of clear categories of objects likely reminiscent of different evolutionary states. For the first time, we have been able to witness directly the link between disks and exoplanets but objects like PDS\,70 are rare since only a few are observable with SPHERE.

For more mature disks, providing some assumptions on the morphology we are now able to infer optical properties of the dust particles based on spectral reflectance measurements and phase function measurements. Here also both polarimetric and non polarimetric imaging were find crucial to bring the interpretation of debris disks data to a higher level than before. Still, the comparison to models is not fully satisfactory and some works are definitely needed toward reconciling the theory of dust scattering with observations. The extreme stability of SPHERE has enabled the investigation of morphological variability, which at least was possible in two cases: First, when the source of variability is located at short physical distances to be able to detect Keplerian motion in just a few months, as we observed with some shadows in young disks. Second, for the particular case of a very active flaring star, like AU\,Mic.

As of now, we understand the limitations of SPHERE and other equivalent instruments, so we can foresee the next steps in the form of upgrades or future instruments on the Extremely Large Telescope. An upgrade of SPHERE is proposed to boost exoplanet and disk science, in particular by pushing the contrast further in to detect and to characterize planets and dust structures at close angular separations corresponding to physical separations where we expect planets to form. A second objective is to allow the observations of fainter and redder stars to observe more PDS\,70 and AU\,Mic analogs and also to reinforce the synergy with ALMA observations sensitive to sub millimeter dust grains.

A completely new field in disk science at high angular resolution and high contrast is now being opened with the James Webb Space Telescope, in particular in the $5-30\,\mu$m spectral range. By providing different types of observables, JWST is intermediate, hence complementary, to scattered light observations in the visible and near IR, and thermal emissions from the dust and gas observed at sub-millimeter. Altogether, these three different regimes will certainly revolutionize our understanding of the birthplace of planets.

\section*{Conflicts of interest}
The author has no conflict of interest to declare.

%\nocite{*}


\bibliographystyle{crunsrt}
\bibliography{crphys20221188
}
\end{document}
