**A cosmological surprise: the universe accelerates.**

Bruno Leibundgut & Jesper Sollerman

Cosmology is in turmoil. The standard model of a few years ago has recently been abandoned and replaced by new ideas. The reasons for this dramatic change are new measurements of the geometry and the matter contents of the universe. The new model implies a dynamical age of the universe that accommodates the oldest known stellar objects, but raises the need for a dark energy component, which is not readily explained within the current particle physics theories.

Most current cosmological models are based on the Big Bang theory, in which the universe started in a hot and dense state. Since then, the cosmic expansion has led to adiabatic cooling, and the subsequent condensations of the matter have formed stars, galaxies and clusters of galaxies. The Big Bang theory predicts not only the universal expansion, but also the baryonic matter content and a relic radiation from the original hot phase. All of these predictions have been observed. The detection of the cosmic microwave radiation and the tiny temperature fluctuations in it are often cited as the most spectacular success of this theory (e.g. Peebles 2001). An extension to the Big Bang theory is the proposed inflationary phase at the very earliest times of the expansion. The inflation would be driven by energy, which emerges from the decay of a particle field and produces a universe many times larger than a simple linear expansion. Inflation predicts the seeds for the growth of large-scale structure in the universe. Since the universe would be inflated several orders of magnitudes, it would essentially have a flat space structure (Guth 1997). Among the open questions in this picture are the matter/energy content of this universe and its future.

The average energy density determines the fate of a homogeneous and isotropic universe governed by gravitational forces. It was thought that the matter content governs the expansion of the universe today, while at early phases radiation dominated. In this simple picture the geometry of the universe is directly coupled to the matter density. A flat universe implies a density of about $8\times10^{-27}$ kg m$^{-3}$, the critical density. For an average density below this value, the universe would be open and expand forever, while at a higher density it would be closed and eventually re-collapse.

Just a few years ago the total energy content of the universe was unknown. All indications were that baryonic matter makes up only about 5$\%$ of the critical density (Bahcall et al. 1999). To explain the kinematics of galaxies and clusters of galaxies an additional mass component often referred to as "missing" or "dark" mass had to be introduced. The search for direct evidence for the dark matter has, despite substantial efforts, been futile so far. Nevertheless, from observations of the evolution and mass determinations of clusters of galaxies, it appears that dark and baryonic matter could explain about 30$\%$ of the critical density. Should the universe be flat, as required if inflation is correct, then some other energy component must contribute to the average density.

For a cosmologist, there are only few parameters needed to describe the universe. All models are based on Einstein's theory of general relativity. The world models are characterized by two parameters: the current rate and the deceleration of the expansion. The first parameter is called the Hubble constant after Edwin Hubble, who discovered the cosmic expansion in 1929. The other parameter describes the change of the expansion and depends on the energy density and the curvature of the universe. The contributions to the density are expressed as fractions of the critical density and are denoted by the Greek letter $\Omega$, e.g. $\Omega_{\rm M}$, for the matter density. The expansion itself is typically measured by the redshift. This is the ratio of the scale factor at two different times of the expansion and observationally a shift of spectral features to longer wavelengths. Hubble's law states that for small distances the redshift is proportional to the distance.

At large look-back times and distances the linearity of Hubble's law breaks down and the distances depend on the energy density of the universe. The various constituents, typically matter and radiation are considered, contribute in different ways to the energy density. Radiation ceased to be gravitationally important at a redshift of about 1000, a time from which we can only measure the cosmic microwave background radiation. Another component is the famous cosmological constant introduced by Albert Einstein to reconcile the solutions of his equations with a static universe. He later abandoned this term, when Edwin Hubble discovered the general expansion of the universe. For many decades the cosmological constant was not considered in the world models as there was not obvious reason to include it and as it was not possible to connect it to any particle theory. In modern terms, it represents the contribution of the vacuum energy (Carroll et al. 1992).

The last three years have seen truly exciting progress in observational cosmology. The flat geometry of the universe has been confirmed by balloon-borne experiments that measured the fluctuations in the cosmic microwave background. The physical scale of these fluctuations can be determined from first principles and the measurement of their angular extend on the sky gives a direct indication of the geometry of the universe. All experiments agree that the universe is most likely flat (de Bernardis et al. 2000, Balbi et al. 2000).

A complementary approach to determine the geometry is through the measurement of distances. This approach actually measures the deceleration due to the gravitational attraction, which slows down the cosmic expansion. The amount of deceleration hence directly determines the average density.

Cosmological distances are, however, notoriously difficult to attain. The long and arduous history of the determination of the Hubble constant painfully reflects this difficulty. Even today, it appears that the value of the Hubble constant is uncertain to about 10$\%$. The easiest way to measure distances is through standard candles, i.e. objects with identical absolute luminosity. Many candidates for standard candles have been proposed. Only few astronomical objects have turned out to be suitable.

A prime candidate for such a standard candle is a certain type of stellar explosions. Known to astronomers as Type Ia Supernovae these explosions occur at the end stages of stellar evolution, when low-mass stars exhaust their fuel and start to contract and cool down. The resulting white dwarfs are so compact and dense that they are supported by electron pressure. It can be shown that there is an upper limit to the mass supported by the pressure of the degenerate electrons. This Chandrasekhar limit, named after Subrahmanyan Chandrasekhar, is near 1.4 times the mass of the Sun. A single, isolated white dwarf will not change its mass, but a white dwarf in a double star system may acquire mass from its companion. If this process is efficient enough, the white dwarf will turn itself into a thermonuclear bomb that will burn carbon and oxygen explosively to nuclear statistical equilibrium. For the density and pressure in these explosions this is mostly a radioactive nickel isotope ($^{56}$Ni), which decays through gamma-decay to $^{56}$Co and further to stable $^{56}$Fe. The energy is deposited in the star and blows it apart. Such an explosion can outshine the light of a whole galaxy made up of some 10$^{10}$ stars (Figure 1). The Chandrasekhar mass is a natural limit and makes it conceivable that all Type Ia Supernovae are similar. Such a uniform constellation assures that there are only small differences between individual explosions, a prime condition for a standard candle.

# Observing supernovae

Supernovae are extremely rare. A galaxy like our Milky Way may produce a Type Ia Supernova only every 400 years. Hence, one has to observe a large number of galaxies for quite some time to detect a supernova. The sample of nearby supernovae (out to a redshift of about z=0.1) is still rather limited. In fact, well-sampled light curves (e.g. Figure 1) are an exception and it has taken over a decade of dedicated effort to collect a significant sample of objects. In the process a lot has been learnt about Type Ia Supernovae and their physics (for a review see Leibundgut 2000).

Not all Type Ia Supernovae are identical. However, there is a way to normalise their peak luminosity according to their light curves. More luminous supernovae display a slower temporal evolution. In this way, it is possible to normalise all objects to the same peak luminosity and make them exquisite standard candles.

One way to empirically test the quality of a standard candle is to look at the distance vs. redshift diagram (Figure 2). For a redshift below z=0.1 the scatter around the line of linear expansion is extremely small and proves that, at least in the nearby universe, the supernovae can be used as distance indicators. Direct determination of the absolute luminosity of Type Ia Supernovae has recently been obtained for the 10 nearest supernovae which exploded within the last 60 years. The peak luminosity of these objects is indeed very uniform.

Distant supernovae, like the example in Figure 1, are more difficult to observe. Since supernovae are such rare events a large volume has to be surveyed to discover a sufficient number for co-ordinated follow-up observations. It was soon recognised that such a project would exceed the possibilities of a single observatory or group of astronomers. In fact, this experiment makes use of all the largest telescopes available. A typical campaign involves the European Southern Observatory VLT (Figure 3), the Keck telescopes, the Canada-France-Hawaii telescope, the Cerro Tololo Inter-American Observatory and the Hubble Space Telescope. The supernovae are observed for about two months with all available telescopes. The resulting light curves, like the one shown in Figure 1, are analysed to derive the peak brightness and plot them in diagrams like Figure 2.

With the model lines calibrated by the nearby supernovae the distant explosions tell us how much the universe has expanded between then and now. From Figure 2 it is obvious that the distances are larger than what would be expected in a freely coasting, i.e. empty, universe. Even worse would be a fit for a flat and matter-filled Einstein-de Sitter model. The distant supernovae are simply too faint, i.e. too distant to be compatible with the old paradigm of a universe filled with matter and radiation only. There has to be a component to the energy density that has accelerated the expansion over the last $\sim6\times10^{9}$ years (Riess et al. 1998, Perlmutter et al. 1999, Hogan et al. 1999).

An obvious candidate is the cosmological constant. With the supernova result it appears that it should be re-introduced to explain the data. Another possibility is that the cosmological constant is indeed zero, but that there is a particle field that through its decay acts like a cosmological constant. It has now come to be known as 'quintessence' (Ostriker and Steinhardt 2001). Independent of which explanation is correct we can designate the energy density of this "dark energy" as $\Omega_{\Lambda}$. Figure 4 shows the probability distribution between the matter density $\Omega_{\rm M}$ and $\Omega_{\Lambda}$ as the supernova distances define them. It is obvious that models with $\Omega_{\Lambda}$=0 are excluded at the $>95\%$ level for all models with a positive matter content. This comes from the fact that an accelerating, negative pressure component is required, which can not be achieved with matter or radiation.

A nice feature of world models with these parameters would be that the dynamical age of the universe is no longer in conflict with the ages of the oldest stars in our Milky Way. We have marked the isochrones of the dynamical age for the combination of $\Omega_{\rm M}$ and $\Omega_{\Lambda}$.

The result from the distant supernovae that the expansion of the universe is not decelerating but in fact accelerates is certainly surprising. We should answer the obvious question how secure is this result is. There are a few other possible explanations for the faintness of the distant supernovae. They are an intrinsic evolution of the peak brightness of the supernovae, unrecognised dust, or gravitational lensing. Evolution has been the downfall of all previously proposed distance indicators. Take for example regular galaxies. Looking back for several 109 years means that a galaxy is composed of stars, which are on average younger than the ones currently in our Milky Way. Since a galaxy with young stars has a larger fraction of short-lived, massive stars, which are brighter and bluer, it will change its luminosity over time. Could Type Ia Supernovae also suffer from evolutionary effects? At first glance, such an evolution is less likely. According to the current models they are explosions of stars, which always end up in the same configuration, a white dwarf at the Chandrasekhar limit. However, even though the bomb always has the same mass, its composition may vary slightly. The explosive carbon burning in a white dwarf depends on the mixture of carbon and oxygen. One could imagine that the distant supernovae have had a different chemical composition at the explosion than their nearby counterparts. Also, stellar evolution predicts that younger white dwarfs are the descendants from predominantly more massive stars. Thus, the parent population of the distant supernovae could be different from the one observed nearby. Unlike for the nearby supernovae we can not use the linear expansion law to check for consistency. What has been tried in the last few years is to make sure that the distant objects have the same appearance as the nearby ones. This program is still under way, but first results are already becoming available. The observed spectra of the two populations are so far indistinguishable (Figure 5). There is further a very characteristic light curve shape with a second maximum of Type Ia Supernovae in the near-infrared, which by now has also been observed for the distant supernovae. On the other hand, there appears to be a trend to bluer colours for distant supernovae, which, if confirmed, may hint at evolution.

Interstellar dust also dims astronomical objects. Galactic dust distributed in the plane of the Milky Way is the reason we can not observe the Galactic Centre in optical light. But galactic dust not only dims the objects, it also preferably scatters blue light and makes objects appear redder. If the intrinsic colour of an object is know, the observed light can be corrected according to the galactic absorption law. The colour evolution of Type Ia Supernovae is in fact very uniform and by always obtaining colour information for the distant supernovae we can check for dust and correct the brightness, if necessary. But there is nothing that says that the galactic dust is representative throughout the universe. For example, larger dust particles scatter differently so that the supernovae still would be dimmed but not as reddened. Observing programs to check this possibility have been carried out last year and are currently being analysed.

Along the light path of any distant object lie other massive objects. Gravitational lensing of distant objects has been observed in several forms. Background objects are distorted by the potential wells through which light has to travel. Large arcs of distorted background galaxies are now regularly observed in massive clusters of galaxies. Another effect of lensing is amplification or de-amplification of the light. It turns out that distant objects are on average de-amplified, i.e. statistically appear fainter. This would be an obvious explanation of the dimness of the distant supernovae. However, the effect is not large enough, even if all mass would be concentrated in compact objects, and can not explain the distant supernovae. This result is based on model calculations, as this is the only effect that can not be observed with the supernova light itself.

There is one signature, which would almost unambiguously prove that the universe has been accelerating over about half its age. This is the deceleration during the early phases of the expansion. The cosmological constant does not change over time, but since the density decreases the gravitational action should have been much larger during the first few 10$^{9}$ years. In Figure 2 this would show up as a non-monotonic evolution. The critical redshift range is around z=1.2, a range that has so far not been explored by the observations. This is stretching the current capabilities of any existing telescope. The supernovae become exceedingly faint and due to the redshift have to be observed at near-infrared wavelengths where the night sky is much brighter than in the optical. Nevertheless, the next observing campaigns will target exactly this redshift range to test whether the acceleration of the universal expansion is indeed the correct interpretation for the dimness of the distant supernovae.

There is a big problem with the cosmological constant. In modern particle theories it is associated with the energy of the vacuum, but these theories also predict a value for the vacuum energy which deviates by more than 50 orders of magnitudes from the cosmological observations. Although this has been a problem all along, the supernova measurements have acerbated this problem, by requiring a non-zero but small ($<1$) cosmological constant. For these reasons many theorists currently favour quintessence models. But this requires introducing a new particle and corresponding potentials, which have to be fine-tuned to have an action like the observed one.

But there is more on the horizon. To distinguish between a cosmological constant and quintessence the time variability of the acceleration should be checked. Tracing the supernova distances in detail can do this. A large sample of supernovae out to a redshift of about z=1.5 is needed. Proposals to obtain such samples have already been made.

Type Ia Supernovae together with the recent results on the cosmic microwave background and the masses of clusters of galaxies are consistent with a flat universe, in which about 30$\%$ is gravitating matter and 70$\%$ is contributed by "dark energy" (cosmological constant or quintessence). Only about 5$\%$ of the total energy stem from baryonic matter. Another 5$\%$ may be contributed by massive neutrinos.

The opinions of cosmologists currently range from visions of "precision cosmology" to worries about the fact that we have to add new constituents to the universe for which we have currently no explanation at all. This is not necessarily a contradiction. Observers have been furnished with tools over the last decade, which allow them to probe many of the cosmological questions in much more detail and with much higher precision. On the other hand, these new results have shown that our picture of the universe was incomplete and will need further scrutiny.


author: T. Bhowmik, A. Boccaletti, P. Thébault, Q. Kral, J.Mazoyer[^1], J. Milli, A.L. Maire, R. G. van Holstein, J.-C. Augereau, P. Baudoz, M. Feldt, R. Galicher , T. Henning, A.-M. Lagrange, J. Olofsson, E. Pantin, C. Perrot
bibliography: all.bib
subtitle: Further evidence of small dust grains
title: Spatially resolved spectroscopy of the debris disk HD 32297[^2][^3]

Spectro-photometry of debris disks in total intensity and polarimetry can provide new insight into the properties of the dust grains therein (size distribution and optical properties). We aim to constrain the morphology of the highly inclined debris disk HD 32297. We also intend to obtain spectroscopic and polarimetric measurements to retrieve information on the particle size distribution within the disk for certain grain compositions.

We observed HD 32297 with SPHERE in Y, J, and H bands in total intensity and in J band in polarimetry. The observations are compared to synthetic models of debris disks and we developed methods to extract the photometry in total intensity overcoming the data-reduction artifacts, namely the self-subtraction. The spectro-photometric measurements averaged along the disk mid-plane are then compared to model spectra of various grain compositions.

These new images reveal the very inner part of the system as close as 0.15$''$. The disk image is mostly dominated by the forward scattering making one side (half-ellipse) of the disk more visible, but observations in total intensity are deep enough to also detect the back side for the very first time. The images as well as the surface brightness profiles of the disk rule out the presence of a gap as previously proposed. We do not detect any significant asymmetry between the northeast and southwest sides of the disk. The spectral reflectance features a "gray to blue" color which is interpreted as the presence of grains far below the blowout size. The presence of sub-micron grains in the disk is suspected to be the result of gas drag and/or "avalanche mechanisms". The blue color of the disk could be further investigated with additional total intensity and polarimetric observations in K and H bands respectively to confirm the spectral slope and the fraction of polarization.

# Introduction

Debris disks are a class of circumstellar disks in which the dust content is thought to be continuously replenished by collisions of planetesimals, usually distributed in a single belt or in multiple belts. The dust present in these disks can be observed by scattered light imaging and/or thermal emission. Morphological and photo-spectroscopic analysis provide constraints on the spatial distribution and physical properties of grains present in a planetary system. Even though debris disks are believed to be depleted in gas, a number of them are significantly gas rich, very likely because of second-generation production of gas released when planetesimals collide and produce the observed dust . In addition, debris disks frequently show identifiable structures like blobs , warps , multiple belts , asymmetries , and other irregularities in the disk. These structures might be sculpted by transient breakups of massive bodies , stellar flybys , stellar companions , or interactions with the inter-stellar medium (ISM) , or be owing to the presence of gas within the disk. One of the most frequently invoked explanations for the presence of most of these structures is, however, the perturbing effect of a planet . In this scenario, structures can be the telltale signature of undetected planets, such as in the case of $\beta$ Pictoris b, which was first inferred indirectly from a warp in the disk and later imaged directly confirming the prediction .

HD 32297 is an A-type star ($V=8.14\pm 0.01,  J=7.687\pm0.024, H=7.624\pm0.051, K=7.594\pm0.018$) located at 133 pc . The age of the star is estimated to be $\gtrsim$`<!-- -->`{=html}15 Myr and $<$`<!-- -->`{=html}30 Myr . A fractional luminosity of L$_{IR}$/L$_\star$ $\geq$ 2.7$\times10^{-3}$ found with the Infrared Astronomical Satellite provides evidence for the presence of cold dust arranged in a belt. This high fractional luminosity makes it one of the brightest debris disks known to date . The disk was first resolved in scattered light by the Hubble Space Telescope (HST) NICMOS as an edge-on system and detected up to radial distances of 400 au . A significant surface brightness asymmetry (southwest ansae brighter than northeast ansae) was reported and attributed tentatively to the presence of an unseen planet. detected the disk at larger separations (400 to 1680 au) from the ground with Keck in the R band and measured a blue color when comparing to the HST near-infrared (NIR) data.

In apparent contradiction to the scattered light observations of and , mid-infrared (MIR) observations published by and show that the northeast side appears brighter than the southwest side beyond $0.75''$, although the impact of the angular resolution and signal to noise ratio (S/N) in these images could be questionable. and observed a similar brightness asymmetry to and in millimeter emission and K band respectively, but still with low angular resolution.

At large distances, as observed with Keck and confirmed with HST/NICMOS, the disk appears bowed, which has been interpreted as the possible interaction of small dust grains residing at large distances with the ISM . confirmed the bow shape of the disk at the L$'$ band. Recently, showed that similar structures could also be the result of an interaction between a planet in an eccentric orbit and planetesimals perturbed by radiation pressure. The presence of small dust grains in a large halo has been confirmed with very deep HST/STIS observations .

With the advance of high-contrast imaging, , and were able to resolve the disk in the NIR as an inclined belt ($\sim$`<!-- -->`{=html}88$^{\circ}$) located at $\sim$`<!-- -->`{=html}130 au (after correcting for the new Gaia distance). According to , there is no significant brightness asymmetry between the two sides of the disk in H and K bands as observed with NACO, while suggest the southwest side to be brighter than the northeast side at $r= 35-80$ au.

High-contrast polarimetry combined with total intensity in the NIR was first achieved by with Subaru/HiCIAO observations, and with the aim to break degeneracies on the geometrical parameters and grain properties. They reported a gap in the total intensity H band data, which was not visible in their polarimetric data.

Finally, the most recent observations were carried out with ALMA by who concluded that at millimeter wavelengths the disk is composed of a planetesimal belt with an inner edge at 78 au and a outer edge at 122 au, and an extended halo up to 440 au. The presence of millimeter grains in the halo complicates the understanding of the large bow, which is expected to be populated with loosely bound grains placed on very eccentric orbits by stellar radiation pressure that should be sensitive to interactions with the ISM. We note however that there are some other alternatives and the ALMA halo may be explained as the presence of a scattered disk .

Furthermore, CO gas emission was detected with ALMA corresponding to a total mass of $\sim 7\times10^{-2}$ $M_{\oplus}$ derived from an optically thin CO isotopolog . This large quantity of gas could interact and drag the dust in this system . The smallest bound grains would be pushed backwards and the unbound grains would be slowed down and may accumulate in greater quantity than in gas-depleted systems.

Modeling of the spectral energy distribution (SED) suggests that the main belt is populated with sub-micron grains , which is in agreement with the color index measured between visible and NIR by . The presence of such sub-micron grains is unexpected, because they should be smaller than the blowout size due to radiation pressure, and should therefore be ejected on very short timescales . Moreover, taking into account the size of the belt inferred by and using photometry from Herschel, were able to draw constraints on the grain composition. These latter authors concluded that the system is made of an inner warm belt at $\geq$`<!-- -->`{=html}1.1 au, and an outer ring at 110 au populated with $\geq2.2\muup$m, high-porosity grains consistent with cometary-like composition. However, this result was questioned by who favored pure icy grains instead.

In this paper we present new high-contrast imaging observations of HD 32297 with Spectro-Polarimetric High-contrast Exoplanet Research at VLT in Chile, in the NIR from Y to H band complemented with polarimetric observations in the J band. With the new observations at NIR wavelengths, we resolve the disk at separations as close as 0.15$''$ and perform a spectral analysis. We briefly discuss the observational technique and the data reduction for total intensity and polarimetric data in Sect. . The disk features seen in the images are described in Sect. . The procedure of the modeling and parametric study is presented in Sect. . In Sect. , we describe how the photometry is retrieved from the images in the various spectral channels, for both total intensity and polarimetry. The resulting spectrum is compared to various grain models in Sect. . Finally, we discuss the implication of the colors in terms of grain sizes in Sect. .

# Observation and data reduction

```latex
\begin{table*}
\centering
\caption{SPHERE observation log}
\begin{tabular}{c c c c c c c c c c}     % 7 columns
\hline\hline                             % To combine 4 columns into a single one
Data UT & prog. ID & Filter & PC & Field rotation ($^{\circ}$) & DIT (s)  & T$_{\mathrm{exp}}$ (s) & Seeing ($''$) & $\tau_0$ (ms) & TN ($^{\circ}$)  \\
\hline
   2016-12-19 & 098.C-0686(A) & IRDIS-BB\_H & &25.16  & 64 & 7168 &0.72 & 8.8 & -1.75 \\
   2016-12-19 & 098.C-0686(A) & IFS-YJ & &25.42 & 64  & 7168 &0.72 & 8.7 & -1.75 \\
   2016-12-16 & 098.C-0686(B) & IRDIS-DPI-BB\_J & 10& stabilized & 64 & 5120 & 0.47& 8.0& -1.7\\
\hline
\end{tabular}
\tablefoot{From left to right: the observation data, program ID, filter combination, the number of polarimetric cycle (PC), the total field rotation in degrees, the individual integration time of each frame (DIT) in seconds, the true time in seconds ($T_{exp}$), the DIMM seeing in arcseconds, $\tau_0$ the correlation time in milleseconds and the true north correction angle in degrees (TN)}
\label{table:obs_log}
\end{table*}
```

## SPHERE

SPHERE is an instrument installed at the VLT for high-contrast direct imaging of giant planets around young nearby stars . Its combination of extreme adaptive optics and advanced coronagraphy has allowed observation of circumstellar disks from the ground. SPHERE consists of three instruments: the Infra-Red Dual-beam Imager and Spectrograph , the Integral Field Spectrograph , and the Zurich IMaging POLarimeter . IRDIS is a dual band imager using two narrowband or broadband filters in the Y, J, H, or K bands . The IFS is a spectro-imager which delivers 39 simultaneous images across the YJ (IRDIFS mode) or YJH (IRDIFS-ext mode) bands. These two instruments can be used for parallel observations with the IRDIFS mode.

## Observations

HD 32297 was observed with SPHERE on December 19$^{}$, 2016, in the IRDIFS mode using pupil stabilization in order to take advantage of angular differential imaging for calibration of stellar residuals during post processing. An apodized Lyot coronagraph N_ALC_YJH_S with a diameter of 185 mas allows attenuation of the starlight in the AO corrected radius (0.84$''$ at 1.65$\muup$m).

For IRDIS the observations were performed in the broadband H filter (1.625 $\muup$m central wavelength, 0.29 $\muup$m filter width). IRDIS has a field of view (FOV) of $11''\times11''$ and pixel size of 12.25$\pm$`<!-- -->`{=html}0.02 mas. For IFS we used the YJ mode (0.95-1.35 $\muup$m), which provides a spectral resolution R$\simeq$`<!-- -->`{=html}54. The IFS has a FOV of $1.73'' \times 1.73''$ and a pixel size of 7.46$\pm$`<!-- -->`{=html}0.02 mas . The atmospheric conditions were good with seeing below 0.8$''$ and correlation time of about $8-9$ ms. The observing sequence is as follows: point-spread function (PSF; with the star outside the coronagraphic mask with a neutral density ND2), star center , long science coronagraphic exposures, second PSF and sky background. Further details on the observation log are provided in Table .

Furthermore, HD 32297 was observed with IRDIS dual polarimetric imaging (DPI) mode a few days apart in field stabilized mode to measure the polarized flux of the disk. Several polarization cycles were taken in the J band (1.25 $\muup$m), each consisting of exposures at four orientations of the half wave plate ($0^{\circ}, 22.5^{\circ}, 45^{\circ}, 67.5^{\circ}$) to measure the full Stokes parameters. The observing sequence is similar to that of the IRDIFS mode for IRDIS; the two filters of IRDIS are replaced with polarizers to split the polarization into two orthogonal directions .

## IRDIFS data reduction

The preliminary data reduction, including flat-field corrections, sky and dark subtractions, star-centering using waffle pattern, bad-pixel removal, distortion correction , and wavelength calibration, is done at the SPHERE Data Centre[^4] using the data reduction and handling (DRH) pipeline for the total intensity data obtained in IRDIFS mode. This step provides a data cube which is further processed with ADI techniques based on several algorithms, such as Karhunen-Loève Image Projection , classical angular differential imaging , and template locally optimized combination of images , using two pipelines, namely the SpeCal pipeline and another one developed by .

Karhunen-Loève Image Projection is an algorithm based on the principle component analysis (PCA). The algorithm involves reformatting the science data into a covariance matrix to remove redundancy. Eigenvalues and eigenvectors are calculated for the covariance matrix which is used to form the KL (Karhunen-Loève) basis, which itself is truncated to a given number of modes . The reference image is built by projecting the science data onto the truncated KL basis and then subtracting them out frame by frame. The result is derotated according to the parallactic angle variation as in any other ADI techniques allowing reduction of the stellar halo. The final image is normalized to the maximum of the measured PSF.

In this paper, we work with the KLIP reduced data truncated to ten modes using the pipeline developed by as it provides the optimal S/N.

## Dual polarimetric imaging data reduction

The stokes vectors Q and U are calculated from the polarization cycle observed through each half-wave plate position. To mitigate the instrumental polarization we use the double subtraction method . The final Q and U vectors can be written as follows:

$$\label{eq1}
\begin{aligned}
Q &= \frac{Q^+ - Q^-}{2}\\
U &= \frac{U^+ - U^-}{2}
\end{aligned}
,$$

where $Q^+, U^+, Q^-, U^{-}$ are obtained from observations through the half wave plate positions at $0^{\circ}, 22.5^{\circ}, 45^{\circ}, 67.5^{\circ}$.

We retrieve the $Q_\phi$ and $U_\phi$ azimuthal Stokes vectors which are expressed in polar coordinates as explained in . The $Q_\phi$ and the $U_\phi$ vectors can be written as:

$$\label{eq2}
\begin{aligned}
Q_\phi &= -Q \cos 2 \phi{+ U} \sin 2 \phi\\
U_\phi &= -Q \sin 2 \phi- U \cos 2 \phi
\end{aligned}
,$$

where $\phi$ is the azimuthal angle with respect to the center of the star. We recover the disk signal in the $Q_\phi$ vector and the noise in the $U_\phi$ vector under the assumption that the disk is optically thin and undergoes only single scattering. Further correction to instrumental offset is not done as the $U_{\phi}$ image is dominated by noise and has very low true disk signal.

# Morphology of the disk

## General description

The full extension of the disk in total intensity (H band) and polarimetry (J band) is shown in the S/N map in Fig. , while a smaller FOV is displayed in Fig. for IRDIS (H band), IRDIS-DPI (J band), and IFS (YJ band, collapsed image). The S/N map is the ratio between the reduced image and the azimuthal standard deviation of the same image

The disk of HD 32297 is relatively bright compared to the stellar residuals. We detect the disk extending to the same distance as observed by HST-NICMOS . The disk is almost edge-on and the shape appears globally symmetrical in both the northeast (NE) and southwest (SW) sides as seen in both the total intensity and polarimetric images. Some local asymmetries between the NE and SW sides of the disk are discussed in detail in the photometric analysis in Sect. .

The IRDIS images feature a concavity at large separations towards the northwest, which becomes significant at a stellocentric distance of about $2''$. This pattern is observed in both the IRDIS BB_H image and the DPI BB_J image (as well observed in Fig.). This concavity or bow shape was also observed with HST/NICMOS as mentioned in , and attributed to the interaction with the ISM of small particles on very eccentric orbits.

On closer inspection at a smaller scale, inside $1''$ the disk takes a half-elliptical shape indicative of an inclined ring (likely the planetesimal belt), as suggested first by and . The high quality of SPHERE images now allows us to resolve this ellipse, with the innermost part of the disk observed at a stellocentric radial distance as close as 0.15$''$ . We measure the position of the ring ansa at $0.8-0.9''$. The asymmetry with respect to the major axis is reminiscent of forward scattering suggesting that the bright part is the front side of the disk. Interestingly, the S/N map (Fig. ) built from the IRDIS BB_H image suggests that the back side of the disk is also detected which is also visible in the IRDIS BB_H intensity image (Fig.). This is discussed in Sect. .

Even though the photometry of the disk is impacted by the so-called ADI self-subtraction , the intensity along the disk varies monotonically, and so does not feature any sign of a gap contrary to the observation reported in . Aditionally, the stellar residual halo in the SPHERE images is much lower than in HiCIAO observations, confidently ruling out such a gap. This is confirmed with our photometric analysis discussed further in Sect. .

The IFS image features a highly symmetrical ring on each side of the minor axis, inside the achievable FOV corresponding to stellocentric distances of $<0.9''$.

In polarimetry, the intensity along the disk varies differently than in total intensity because the scattering phase function is vastly different. As a result, we observe a peak of intensity located at $\sim 0.8''$ and associated to the location of the disk ansa. Due to the absence of self-subtraction in DPI, strong signal is observed at small angles \< 0.15$''$ (Fig. ). At this stage, it is not yet clear whether this signal is produced by stellar residual or by the disk itself, although it is moderately visible in the S/N map (Fig. ). There are two possible explanations for this signal to be a true disk signal, one being a strong forward scattering peak at small angular separation, and another being a detection of a potential inner belt.

## Position angle of the disk

The position angle (PA) is determined by a method developed for edge-on disks as presented in . The science data are first rotated to an initial guess of the PA, and a Gaussian profile is fitted perpendicular to the disk mid-plane in a range of angular separations from which the spine of the disk and the local slope, globally or separately for the NE and SW sides, are derived. This process is repeated until the slope reaches a minimum corresponding to a horizontal disk and providing a measurement of the actual disk PA. The average PA retrieved for both total intensity and polarimetric data is $47.60^{\circ}\pm0.2^{\circ}$ and $47.50^{\circ}\pm0.15^{\circ}$ , respectively. The errors consist of measurement error for the applied method and the additional TN uncertainty of 0.1$^{\circ}$.

The local minima in the spine of the disk can determine the position of the ansae . For higher accuracy, a simple ellipse fit to the spine of the disk can give the position of the ansae corresponding to the semi-major axes of the ellipse. Additionally, the inclination can be derived from the semi-major and semi-minor axes of the fitted ellipse. The ellipse fitting the total intensity data is centered at (-0.06$''$,-0.008$''$) and (-0.02$''$,-0.015$''$) for the polarimetric data. The plots are presented in Fig. . We do this ellipse fit to the spine measured on all the images reduced by the methods, TLOCI , KLIP 3,5,10 modes , classical ADI , No ADI , and DPI $Q_{\phi}$ , and the dispersion between these measurements is used as error. The position of the ansae measured from the spine is $126.4\pm{12.8}$ au and the inclination is $88.4^{\circ}\pm{0.6}^{\circ}$.

# Geometrical modeling

## Modeling total intensity with a single phase function

Angular differential imaging techniques induce biases to the disk photometry due to self-subtraction. To overcome this problem and recover the unbiased photometry of the disk we proceed with forward modeling, which involves the generation of synthetic images from a model, given some parameters, undergoing a similar post processing to that used on the original data. We used the GRaTer model to generate synthetic images of debris disks . The model assumes a ring of planetesimals releasing dust from a collisional cascade and located at a distance $r$ from the star. The position $R_0$ is where the dust density peaks corresponding to the position of the ansae and we then impose that the dust density distribution scales radially as ${\scaleto{r^{\alpha_{in}}}{9pt}}$ inwards and ${\scaleto{r^{\alpha_{out}}}{9pt}}$ outwards. The vertical distribution is fixed to be a Gaussian function, the height of which is controlled by the parameter $H_0$ obtained at a radius $R_0$.

From this geometrical prescription, GRaTer calculates, for a given model, the resulting scattered light image taking into account the scattering angle $\theta$, which depends on the inclination and PA. The phase function is given by the Henyey Greenstein (HG) function for total intensity

$$\label{e3}
\begin{aligned}
f_I (\theta) = \dfrac{1-g^2}{4\pi(1+g^2-2g\cos{\theta})^{3/2}}
\end{aligned}
,$$

where $g$ is an anisotropic scattering factor that we leave as a free parameter. For $0<g<1,$ the scattering by dust particles is predominantly forward (isotropic if $g=0$), and conversely backward for $-1<g<0$.

The free parameters of our model are listed below. The initial guesses are obtained from previous work on this system and first-order estimations from the SPHERE images, totalling 23040 models.

-   inclination i ($^{\circ}$): 87.5, 88.0, 88.5, 89.0;

-   position of the ansae $R_0$ (au): 115, 120, 125, 130,135, 140,145,150;

-   power-law index $\alpha_{in}$: 2, 5, 8, 10;

-   power-law index $\alpha_{out}$: -4, -5, -6, -7,-8 ;

-   HG parameter $g$: 0.4, 0.5, 0.6, 0.7, 0.8, 0.9;

-   disk aspect ratio $h = H_{0}/R_{0}$ : 0.01, 0.015, 0.02, 0.025, 0.03, 0.035.

The PA is kept constant at 47.6$^{\circ}$  as a result of the analysis presented in Sect. . Each GRaTer model is first injected into an empty data cube void of noise and convolved with the measured PSF. To account for the ADI photometric bias in the forward modeling, a given GRaTer model defined by a set of parameters is projected onto the KL basis which was formed and used for the data. Each model is then normalized to the maximum PSF intensity as formerly done for the data (Sect. ). Looking for the best-fit model implies that we compare the science image with a series of model images in an appropriate region encompassing the pixels containing disk signal. A numerical mask $\sim0.15''\times 3''$, aligned with the disk mid-plane, is applied, which is compared to observations in a $\chi^2$ fashion. The central part with a stellocentric distance lower than $\sim0.15''$ is removed from the mask. The northeast and the southwest parts are analyzed separately. A comparison between the science image and the noiseless model masked with an effective aperture is credible as the disk has high S/N assuming that there is no over-subtraction due to the combination of noise and disk. Also, in a comparison between forward modeling with noiseless models and the injected model into the science image at a different PA with a certain flux resulted in similar $\chi^2$ values. Therefore, we refrain from doing the latter due to reduced computational momentum.

The reduced $\chi^2$ is calculated between the science image ($S_{i,j}$) and the models ($M_{i,j}$) at $i,j^{th}$ pixel and summed over the number of pixels in the mask ($N_{data}$). The calculation is described as follows:

$$\label{eq4}
\begin{aligned}
\chi_\nu^2 =\dfrac{1}{\nu}  \sum_{i,j=1}^{N_{data}}\Bigg( \dfrac{S_{i,j}-a . M_{i,j}(p)}{\sigma_{i,j}}\Bigg)^2,
\end{aligned}$$

where $\nu$ is the degree of freedom equal to $N_{data} - N_{param}$, where $N_{param}$ is the number of free parameters and $p$ the parameter space explored in these models. The noise term ($\sigma_{i,j}$) is derived from the azimuthal standard deviation in the image masking the disk. In high-contrast imaging, the noise has a spatial structure related to the stellar residuals, varying with wavelengths, while $\chi^2$ minimization applies to Gaussian errors and linear models. These two conditions are not rigorously met in our forward modeling technique, which imposes some limitations to this approach. The parameter $a$ is the scaling factor between the data and the model.

To identify the best models, we select 1$\%$ of those with the lowest reduced $\chi^2$ values. This approach is considered instead of taking models corresponding to the 1$\sigma$ deviation of the reduced $\chi^2$ distribution using the theoretical threshold of $\sqrt{2\nu}$ . This is because there are very few models falling into this latter category, which puts overly strong restrictions on the measurements of error bars. Also, this threshold is theoretically applicable for Gaussian errors and linear models which, as discussed earlier, our models do not follow. As a final output we derive the histogram for each parameter, for the given set of best models. A Gaussian profile is fitted to the histograms and Table provides the peak of the Gaussian and the $1\sigma$ deviation from the measured peak as errors. Albiet, in two scenarios a Gaussian cannot be fitted. The first is when the distribution of the histogram is flat. The second is when the number of values explored for the parameter are less than four. The best-fit model taken for creating the residual image as shown in Fig. is the model which has the minimum value of the reduced $\chi^2$ (${\chi^{2}=5.26}$).

|  | NE side | SW side | Full disk |
|:--:|:--:|:--:|:--:|
| i ($^{\circ}$) | 88.2$\pm{0.4}$ | 88.3${\pm{0.3}}$ | 88.3${\pm{0.3}}$ (88.5) |
| $R_0$ (au) |  136.8${\pm{10.0}}$ | 129.4${\pm{8.4}}$ | 134.7${\pm{9.3}}$ (130) |
| $\alpha_{in}$ |  6.00${\pm{4.00}}$ |  6.00${\pm{4.00}}$ | 6.00${\pm{4.00}}$ (10) |
| $\alpha_{out}$ | -6.0${\pm{0.9}}$ | -5.7$\pm{0.8}$ | -6.17${\pm{1.17}}$ (-6) |
| $g$ | 0.55${\pm{0.14}}$ | 0.49${\pm{0.14}}$ | 0.55$\pm{0.13}$ (0.6) |
| $h$ |  0.026${\pm{0.005}}$ | 0.024$\pm{0.006}$ | 0.020${\pm{0.006}}$ (0.020) |

Parameters that provide the best GRaTer model fitting of the IRDIS BB_H science image with the one HG phase function

The inclination we find, $88.3^{\circ}\pm0.3^{\circ}$, is compatible with previous measurements but achieves a higher accuracy. The position of the ansae $R_0$ is ${134.7\pm9.3}$ au. The variation of the dust density inwards from the belt ($\alpha_{in}$) is relatively difficult to constrain for inclined disks; only limited values are explored for this parameter and evidently the retrieved value of ${6.0\pm4.0}$ is not constrained. On the other hand, the values for $\alpha_{out}$, $g,$ and $h$ are relatively well constrained with values of ${-6.17\pm1.17}$, $0.55\pm0.13,$ and ${0.020\pm0.006}$, respectively.

From Fig. , one can see that the residuals are still relatively large, which indicates that the adopted model does not perfectly explain the data. Using two HG phase functions instead of one has been proved effective to better fit inclined disks . Similarly, we also try to improve our model by using two HG phase functions to model the disk.

## Modeling total intensity with two HG phase functions

The previous section uses a single HG function to model the scattering phase function but allowed to explore a large range if parameters. In the total intensity IRDIS image, a faint signature of the back side is visible (Fig. ). In trying to better model the back scattered grains in the disk, we adopt a new phase function based on two HG functions as shown below.

$$\label{eq5}
\begin{aligned}
f_T (\theta) = w_1.f_I (g_1,\theta) + (1-w_1).f_I (g_2,\theta)
\end{aligned}
,$$

where $g_2$ is assumed to be negative as it models the backward scattering component of the disk.

We reduced the parameter space guided by our results from the previous section and thus the explored parameters for this case are:

-   inclination i ($^{\circ}$): 88.0, 88.5;

-   position of the ansae $R_0$ (au):120, 125, 130, 135, 140, 145;

-   power-law index $\alpha_{in}$: 2, 5, 8, 10;

-   power-law index $\alpha_{out}$: -5, -6, -7;

-   first HG parameter g$_1$: 0.4, 0.5, 0.6, 0.7;

-   second HG parameter g$_2$:-0.5, -0.4, -0.3, -0.2, -0.1;

-   disk aspect ratio $h$ : 0.02, 0.025, 0.03;

-   weight $w_1$: 0.80, 0.83, 0.87, 0.90, 0.93.

The range of weights $w_1$ is chosen to encompass the values derived for and another inclined debris disk .

With the combination of the above listed parameters, 43200 models are created, for which the reduced $\chi^2$ is measured and the best models are identified following the same approach as described in Sect. . The parameters giving the best-fit model for this case are provided in Table .

|  | NE side | SW side | Full disk |
|:--:|:--:|:--:|:--:|
| i ($^{\circ}$) | 88.0 | 88.2${\pm{0.2}}$ | 88.2${\pm{0.2}}$ (88.0) |
| $R_0$ (au) | 132.4${\pm{8.3}}$ | 128.5${\pm{8.6}}$ | 134.4${\pm{8.5}}$ (140) |
| $\alpha_{in}$ | 5.2${\pm{2.8}}$ | 6.0$\pm{4.0}$ | 6.0${\pm{4.0}}$ (2) |
| $\alpha_{out}$ | -6.0${\pm{1.0}}$ |  -6.0${\pm{1.0}}$ |  -6.0${\pm{1.0}}$ (-6) |
| $g_1$ |  0.68${\pm{0.05}}$ | 0.69$\pm{0.06}$ |  0.69${\pm{0.06}}$ (0.7) |
| $g_2$ | -0.3${\pm{0.2}}$ | -0.3${\pm{0.2}}$ | -0.3${\pm{0.2}}$ (-0.4) |
| $w_1$ | 0.82$\pm{0.04}$ | 0.80${\pm{0.06}}$ |  0.81${\pm{0.05}}$ (0.80) |
| $h$ | 0.025${\pm{0.005}}$ | 0.022$\pm{0.003}$ |  0.022${\pm{0.002}}$ (0.020) |

Parameters that provide the best GRaTer models fitting the IRDIS BB_H KLIP science image with two HG phase functions

The inclination is very well constrained to ${88.2\pm{0.2}}$$^{\circ}$, but we highlight the fact that we considered only two values to explore the parameter space as this is relatively well-constrained from previous studies as well as in the previous section (Sect. ). The position of the ansae and the inner power-law index are found to be ${R_0=134.4\pm{8.5}}$ au and ${\alpha_{in}=6.0\pm{4.0}}$, and ${\alpha_{out}}$ is ${-6.0\pm{1.0}}$ which are within the error bars retrieved in the previous section. The value of ${h = 0.022\pm{0.002}}$ is also consistent with results from the previous section.

We find ${g_1 = 0.69\pm{0.06}}$, ${g_2=-0.3\pm{0.2,}}$ and ${w_{1} = 0.81\pm{0.05}}$. These values do not agree with the best fits in K band in , where a similar HG function to Eq. was used to model the phase component with $g_1 = 0.96$, $g_2=-0.1,$ and $w_{1} = 0.9$. It should however be considered that our data have better S/N compared to the observations of and we observe the back side of the disk in our data.

In Fig. , we plot our best-fit model together with the data and the residuals found after subtracting them both. The reduced $\chi^2$ value for the best-fit model is 3.29. We find that our model provides a 31% better match to the data compared to models with single HG phase function. Irrespective of the improvement, some residuals remain, which may be due to limitations that are combination of several factors. The simplicity of our model with only eight parameters cannot reproduce the complexity of the disk; the variation of the PSF during observation, the possibility of over-subtraction and indirect self-subtraction in the forward modeling , along with nonlinear terms ignored in this study can all add to limitation of our best-fit model.

Independently fitting the disk observed in each spectral channel of IFS with the same set of parameters ($R_0$, $\alpha_{in}$, $\alpha_{out}$, h) but $g_1$, $g_2,$ and $w_1$ being free parameters, we found that the variation of the anisotropic scattering factors is no larger than 8%. This gives confidence in selecting the same value for all spectral channels. We note that the disk is observed only partly in the IFS FOV, and therefore any degeneracy between morphological and phase parameters is not considered in this test.

## Modeling polarimetric images

Polarimetric observations complement total intensity data and depend differently on the morphological parameters. For instance, as shown in , the phase function in polarimetry peaks close to $\sim90^{\circ}$ phase angle as opposed to the total intensity phase function, which reaches a maximum at small phase angles. As a result, the ansae has stronger signatures in polarimetry than in total intensity. Modeling the polarimetric image is therefore crucial to derive the geometry of the debris disk .

We use GRaTer to create a grid of 32256 geometrical models with a polarised phase function $f_P (\theta),$ given below:

$$\label{eq6}
\begin{aligned}
f_{P}(\theta) = f_{I}(\theta)\dfrac{1-\cos^2 \theta}{1+\cos^2 \theta}
\end{aligned}
,$$

where $\theta$ is the scattering angle. The phase function used is the combination of the HG function (Eq. ) and Rayleigh scattering, to account for the angular dependence of linear polarization due to single scattering of an optically thin disk.

The back side of the disk is not visible in the polarimetric image, and therefore we used a single HG function to construct the phase function in the GRaTer models. The free parameters are as follows:

-   inclination i ($^{\circ}$): 87.5, 88.0, 88.5, 89.0;

-   position of the ansae $R_0$ (au): 115, 120, 125, 130,135, 140,145,150;

-   power-law index $\alpha_{in}$: 2, 5, 8, 10;

-   power-law index $\alpha_{out}$: -3, -4, -5, -6, -7,-8;

-   HG parameter $g$: 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.99;

-   disk aspect ratio $h$: 0.01, 0.015, 0.02, 0.025, 0.03, 0.035.

We applied the same procedure as in Sect. to measure the reduced $\chi^2$ and derive the range of best models. The value of the reduced $\chi^2$ is 3.52 for the best-fit model.

|  | NE side | SW side | Full disk |
|:--:|:--:|:--:|:--:|
| i ($^{\circ}$) | 88.7$\pm{0.3}$ | 88.5$\pm{0.3}$ | 88.6$\pm{0.3}$ (88.5) |
| $R_0$ (au) | 135.8${\pm{9.8}}$ | 127.3${\pm{9.1}}$ | 127.9${\pm{8.0}}$ (125) |
| $\alpha_{in}$ | 7.4${\pm{3.2}}$ | 8.1${\pm{3.1}}$ | 8.1$\pm{3.2}$ (10) |
| $\alpha_{out}$ | -4.1$\pm{1.0}$ |  -3.7${\pm{0.7}}$ | -3.9${\pm{0.8}}$ (-4) |
| $g$ | 0.85${\pm{0.09}}$ | 0.88$\pm{0.05}$ |  0.84${\pm{0.08}}$ (0.8) |
| $h$ | 0.023$\pm{0.006}$ | 0.022$\pm{0.006}$ | 0.022 ${\pm{0.006}}$ (0.020) |

Parameters that provide the best GRaTer model fitting of the IRDIS BB_J science image

The parameters for the best-fit model of the polarimetric image are given in Table . We find that the inclination of the disk is 88.6$^{\circ}$$\pm{0.3}$$^{\circ}$and the position of the ansae (${127.9\pm{8.0}}$ au) is consistent with the one obtained in total intensity within error bars. Regarding the slopes of the surface density profile, polarimetry favors a steep inner edge ($\alpha_{in} =8.1\pm{3.2}$) while this parameter was essentially unconstrained in total intensity. On the contrary, the outer slope ${\alpha_{out}=-3.9\pm{0.8}}$ is flatter.

The anisotropic scattering factor is significantly larger in polarimetry ${(0.84\pm{0.08}})$ as compared to intensity irrespective of the number of $g$ parameters we considered. Therefore, the dust grains which are probed in these data are more prone to forward scattering in polarimetry. It should be noted that the in case of strong forward scattering, the polarised phase function peaks at a smaller scattering angle as well as at the ansae , which is visible in our best-fit model as seen in Fig. . Figure plots the scattering phase functions for total intensity and polarimetric best-fit models. The scattering angles probed for HD 32297 are between 6$^{\circ}$ and 175$^{\circ}$. The phase function in polarimetry peaks at scattering angles smaller than $\sim90^{\circ}$ as seen in Fig.  while in total intensity we observe an increase beyond 110$^{\circ}$ due to the back-scattering for a double HG function.

Finally, although the residuals displayed in Fig. are much lower than for the modeling of total intensity, there is still some intensity left near the ansae, indicating that the model does not perfectly reproduce the disk. The residuals could also be an indication of our preference of larger $g$ values over a possibility of an inner component at a separation of \< 40 AU. Modeling the polarimetric observation with consecutive inner and outer belts is beyond the scope of this paper

# Photometry and analysis

## Self-subtraction profile and photometry for the total intensity

There are several possible approaches to retrieving the spectrophotometry of the disk. We compare two methods which we use to retrieve the surface brightness of the disk and discuss the limitations associated with each of them.

-   *Method 1*: The first method is similar to the one used to measure the photometry of the HD 32297 disk in and HD 15115 in . It was developed for highly inclined debris disks and relies on estimating the self-subtraction caused by the ADI process. The edge-on geometry allows simplification of this calculation to a 1D problem. Given the best fit model fitting the data identified in Sect. , we first extracted the radial profile of the GRaTer model (convolved with the PSF) and its associated KLIP image. The profiles for both the model and the KLIP image are measured in the same numerical mask as for the $\chi^2$ minimization, and we average the flux in nonoverlapping concentric arcs of $0.15''$ vertical width, and four-pixel ($\sim0.05''$) length.

    The self-subtraction is derived from the ratio of these two profiles (Fig. , left) and can be significant at short separations when using KLIP (about 60 at $0.2''$ in the H band). The unbiased surface brightness (in counts) of the disk is then obtained by multiplying the self-subtraction profile with the science (KLIP) image profile to compensate for the ADI effect, in each spectral channel of the IFS and each filter of IRDIS.

    Self-subtraction could also be deduced directly from the ratio of the images of the GRaTer model (PSF convolved) and its KLIP version instead of using profiles. used this method in the case of a less inclined disk and using a specific ADI process (masked classical ADI), which minimizes the self-subtraction beforehand. Here, in the case of a highly inclined disk together with KLIP processing, this method would lead to the self-subtraction profile presented in the right panel of Fig. , which clearly cannot be used for photometric correction. The large variations are attributed to the strong positive to negative fluctuations resulting from the KLIP processing. Therefore, we conclude that this method of self-subtraction measurement should be avoided for KLIP-processed data.

-   *Method 2*: The second solution does not require evaluation of the self-subtraction, but instead uses one step of the modeling when the reduced $\chi^2$ is calculated. For the minimum value of reduced $\chi^2$, the scaling parameter $a$ in Eq. directly provides the scaling factor between the KLIP image of the best-fit model and the data. Therefore, contrary to method 1, the surface brightness profile is evaluated in the scaled GRaTer model (PSF convolved) image instead of the real disk image.

For further calculation we first obtain the stellar flux by integrating over a masked PSF which contains 99.99% of its total flux. The radius obtained for the mask is 0.4$''$ for the PSF obtained by IRDIS and 0.3$''$ for IFS. Surface brightness profiles are then converted to magnitude arcsec$^{-2}$ taking into account the pixel size and normalizing with respect to the stellar flux. The error bars linearly combine two terms, the dispersion of the stellar+background residual intensity as measured in the mask rotated by 90$^{\circ}$ relative to the real disk for each radial bin, and the accuracy of the photometric extraction. The latter is estimated with a fake disk (same parameters as the best-fit model) injected into the raw data at 90$^{\circ}$ from the real disk at a contrast level of $5\times10^{-4}$ and processed in the same way. The photometric extraction is found to be consistent with this initial contrast within $\sim3.3\%$.

The two methods provide very similar results as seen in Fig. . Hence, we derived surface brightness profiles for all spectral channels of IFS and IRDIS. Flowcharts representing the calculations using either method are shown in Fig . = \[rectangle, rounded corners, minimum width=3cm, minimum height=1cm,text centered, draw=black, fill=blue!10\] = \[rectangle, rounded corners, minimum width=3cm, minimum height=1cm,text centered, draw=black, fill=red!10\] = \[thick,-\>,\>=stealth\]

At this stage, and contrary to , we do not confirm the presence of a dip or break in the surface brightness profiles near 0.75$''$ on the NE side or 0.65$''$ in the SW, and therefore we rule out the presence of a gap (Fig. ).

## Photometry for polarimetric images

In polarimetry, the disk photometry is not affected by self-subtraction and therefore can be directly measured from the image. Contrast is calculated directly from the science DPI BB_J image, instead of a model, using the same process as explained in method 2 of Sect. .

As a result of the scattering angle dependence (Sect. ), the slope of the surface brightness profile is clearly different in polarimetry compared to total intensity, with a less steep decrease from 0.1$''$ to about $0.7''$ (Fig. ). At the location of the ansae ($\sim0.9''$), the surface brightness shows a peak instead of a break. As advocated in for the case of the inclined debris disk around HD 15115, this dependence of the phase function provides a better sensitivity for polarimetric data to pinpoint the inner edge or the ansae of a dust belt, or to reveal multiple structures. In the case of HD 32297, the polarimetric surface brightness does not show any particular signs of such multiple belts. The polarimetric surface brightness is about two magnitudes fainter than in total intensity for a stellocentric distance of 0.5$''$, which translates to a polarized fraction of about 15%.

## Average spectral reflectance

To derive the reflectance of the grains as a function of wavelength, we converted the surface brightness profiles into contrast with respect to the star and averaged the values between $0.2''$ and $0.8''$, for all spectral channels of IFS and IRDIS, as well as for polarimetric data. Even though, the intensity decreases by $\sim$ 1  mag/arcsec$^2$ between 0.2 and 0.8$''$ for all spectral channels, this separation range is chosen as it corresponds to the minimal and maximal distances where the disk signal starts to dominate over the stellar halo ($>0.2''$) while encompassing the IFS field of view. The reflectance spectra obtained with the two methods of Sect. are plotted in Fig. , in which the error bars are obtained from the averaged errors of the surface brightness between 0.2$''$ and 0.8$''$ for each wavelength. The main characteristic of the reflectance, irrespective of the method used, is a slow decreasing trend with wavelength which gives the disk a gray to blue color in the YJH spectral range (Fig. ). Independently of any assumptions on the grain properties, by fitting a straight line to the spectra, we measured a contrast arcsec$^{-2}$ variation of $-0.013\pm{0.002}$ per $\muup$m.

Method 1 (Fig. , left) clearly produces more dispersion in the YJ band as a result of biases introduced by the self-subtraction estimation. Since method 2 (Fig. , right) is based on photometry extracted from models, it provides a smoother spectrum, but a slightly lower reflectance ($\sim$`<!-- -->`{=html}12 to 20%) than method 1. For further assurance, we checked the consistency of the two methods by injecting a fake disk 90$^{\circ}$to the real disk PA and performing consecutive photometry for the real disk and the fake disk with both methods. We find that the uncertainty on the measurement between the real disk and the fake disk in Method 1 is twice that found when using Method 2. We therefore chose to use Method 2 for subsequent analyses.

Using previous measurements from the literature to confirm the trend of the reflectance on a larger spectral range would be valuable, but photometric measurements are usually derived from various methods, and at various locations on the disk image, which makes the comparison difficult. However, in the case of HD 32297, Ks band data were obtained by some of us using similar (but not quite identical) methods . In these data, the disk is only detected at stellocentric distances of 0.5-0.8$''$ for both the NE and SW sides and the surface brightness is $\sim15\pm{0.5}$ mag/arcsec$^2$ which translates to a contrast of $(1.09\pm{0.5})\times10^{-3}$ arcsec$^{-2}$.

Considering the variation of the disk intensity in SPHERE images between the two ranges of separations $0.2''-0.8''$ and $0.5''-0.8''$ at H band, we can extrapolate this Ks band contrast to $(1.99\pm{0.5})\times 10^{-3}$. Therefore, at first order, the photometry in the Ks band confirms the spectral slope of the reflectance (Fig.). Observing the disk at the Ks band with SPHERE would further validate the current value of the slope found.

The globally blue behavior of the spectrum is a key element that can help to constrain the size distribution of the grains composing the HD 32297 disk.

# Grain modeling

Now that we have constrained the geometry and morphology of the disk, we return to the GRaTer radiative transfer code to constrain the particle size distribution (PSD) of the dust grains. We take as a reference morphology ($R_0\pm{dr}$) the best fit obtained for the "full disk" fit displayed in Table and produce synthetic spectra of the disk as a function of grain-related parameters, which we then compare to the observed spectrum in the NIR.

We consider three different grain compositions: astro-silicates, porous astro-silicates with 80$\%$ porosity, and a mixture with 50$\%$ water ice and 50$\%$ astro-silicates. The two crucial free parameters are the minimum grain size $s_{min}$ and the index $\kappa$ of the power-law that the PSD is assumed to follow ($dn(s) \propto s^{\kappa}ds$). We explore values of $s_{min}$ between 0.1 and 10 $\muup$m with an increment of 0.1 $\muup$m and values of $\kappa$ ranging from -5.0 to -3.0 with an increment of 0.1.

Each model spectrum is interpolated on the forty wavelength channels (39 for IFS and 1 for IRDIS) and globally scaled to the data using the same $\chi^2$ minimization framework as in Sect. .

For all wavelength channels ($i$), the reduced $\chi^2$ is measured between the average spectral reflectance $D_{i}$ and each grain model spectrum $M_{i}$ as given below

$$\label{eq7}
\begin{aligned}
\chi_\nu^2 =\dfrac{1}{\nu}  \sum_{i=1}^{40}\Bigg( \dfrac{D_{i}- M_{i}(p)}{\sigma_{i}}\Bigg)^2
\end{aligned}
,$$

where $\sigma_{i}$ is given by the error bars for each data point in Fig. and $p$ is the parameter space. There are 40 wavelength channels ($n_{\rm data}$) and two parameters ($n_{\rm param}$; minimum grain size and distribution index), and therefore there are 38 degrees of freedom; $\nu=38$. The best models are selected as per 1$\sigma$ deviation of the reduced $\chi^2$ distribution, which is given by $\chi_{\nu,th}^2 = \chi_{\nu,min}^2 + \Delta{\chi_\nu^2}$, with $\Delta{\chi_\nu^2} = \sqrt{2\nu}$.

The reflectance spectrum corresponding to the best fits obtained for the three considered compositions is displayed in Fig. (right) and the best-fit parameters are provided in Table . An important result is that, for all considered compositions, $s_{min}$ is well below the blow-out limit size $s_{blow}$. We indeed always have $s_{min}/s_{blow} \leq 0.085$, where we derive $s_{blow}$ in gm/cm$^3$ using the prescription by :

$$\label{eq7a}
\begin{aligned}
    s_{blow}=0.8\dfrac{{\rm L}_*}{{\rm M}_*}\dfrac{2.7}{\rho}
    \end{aligned}
,$$

where L$_*$ and M$_*$ are the stellar luminosity and mass expressed in solar values, and $\rho$ is the bulk density of the material (given in Table ). We take L$_*= (8.4\pm{0.2})$ $L_{\odot}$ and M$_* =1.8$M$_\odot$ .

```latex
\begin{table*}[h]\centering
\caption{Parameters of the grains and their size distribution that generate the best fit to the spectrum}
\begin{tabular}{c | c|c |c | c | c | c |c  }     % 7 columns
\hline \hline

  Grain type & volume ratio &density g/cm$^3$& $\kappa$ & $s_{min}$ & $s_{blow}$ & $s_{min}/s_{blow}$ & $\chi^2$\\
\hline
  Astro-silicate          & -   & 2.7  & $-3.79\pm{0.34}$ & $0.30\pm{0.02}$ & 3.7 & 0.081 & $0.471\pm{0.084}$ \\
  Vacuum+astro-silicate  & 4:1 & 0.54 & $-4.24\pm{0.60}$ & $1.59\pm{0.11}$ & 18.7 & 0.085 & $0.628\pm{0.114}$ \\
  Astro-silicate+water ice & 1:1 & 1.85 & $-4.65\pm{0.35}$ & $0.44\pm{0.03}$ & 5.4 & 0.081 & $0.502\pm{0.091}$ \\
\hline \hline

\end{tabular}
\label{table:grain}

\end{table*}
```

As for the best fit of the slope (Fig., right) of the size distribution, we find $\kappa=-3.79\pm{0.34}$ for astro-silicates, which is relatively close to the slope expected for collisional steady states . For the other two compositions, we find slightly steeper PSDs, with $\kappa\sim-4.5$. For all three slopes, it is important to stress that the geometrical cross section, and thus the flux, is dominated by the smallest grains in the PSDs, that is, those close to $s_{min}$ . Also, $s_{min}=2.2\muup$m, the minimum size found by , also corresponds to $s<s_{blow}$ grains given the very high 90% porosity they assumed. To further stress the crucial role of unbound $s\leq s_{blow}$ grains, we display in Fig. (left) the best fits that would be obtained when considering a PSD stopping at $s_{min}=s_{blow}$ and a canonical slope with $\kappa = -3.5$. As can be clearly seen, in the absence of the unbound grains no satisfying match can be obtained, in particular regarding the blue slope of the reflectance spectrum.

These results, especially those regarding the minimum grain size, only weakly depend on the morphology assumed for the disk. Taking a size for the planetesimal belt other than that of Table indeed leads to relatively similar results. This is expected, as the slope of the reflectance spectrum is essentially imposed by the size-dependence of the scattering coefficient $Q_{sca}$, which does not depend on the location of the grains with respect to the star. The scattering anisotropy parameters $g$ do depend on this location, but under the assumption made in Sect. that there is no size dependence for $g$, this would not translate into a different slope of the reflectance spectra for a given PSD. Therefore, we note that at this stage it is safe to derive synthetic spectra without a dependence of anisotropic scattering parameter.

# Discussion

## Comparison of geometrical parameters to millimeter observations by ALMA

From the analysis presented in Sect. we find that the disk is best described as a relatively narrow belt peaking at $132.3\pm{6.2}$ au (according to polarimetric data), which agrees with previous scattered light observations . HD 32297 was also observed recently with ALMA in the dust continuum at 1.3 mm by MG18. However, the beam size was about ${0.76''\times0.51''}$ which is more than ten times larger than the SPHERE angular resolution, meaning that the ALMA resolution is significantly poorer.

While considering a geometrical model that is similar to GRaTer, MG18 concluded that the disk is rather broad, extending from $78\pm8$ au to $122\pm3$ au, with a surface density rising as $r^2$. The peak in density occurs slightly closer-in than in the SPHERE images, but this discrepancy could be the result of the angular resolution. In any case, the radial dependency of the surface density, even if relatively poorly constrained in our case, is considerably steeper with SPHERE.

Beyond the planetesimal ring at 122 au, MG18 observed a halo extending out to $440\pm{32}$ au, where the surface density decreases as $r^{-6}$. This is consistent with the SPHERE image in which we observe dust scattering as far as 3.3$''$ (equivalent to 440 au), with a comparable radial decrease in density. MG18 do not report any asymmetry in this outer part while the aforementioned concavity is obvious at shorter wavelengths. This could again be a resolution effect or because of the sensitivity of ALMA to bigger grains. Nevertheless, subtracting an axi-symmetrical disk model from the ALMA image leaves residuals co-located with the region where the concavity is detected by SPHERE and HST. The presence of millimeter grains at such stellocentric distances would, however, be at odds with the expected mechanism creating such concavity. For the same reason, our value for the disk inclination is significantly more precise ($88.2\pm{0.3}$$^{\circ}$) than the one derived by MG18 ($83.6^{\circ}$ $^{+4.6}_{-0.4}$).

In order to check the compatibility of the ALMA and SPHERE models, we used the MG18 parameters $\alpha_{in}$, $\alpha_{out}$, $R_{in}$, $R_{out}$, and took the other parameters ($g_1,g_2,w_1,h$) from our best-fit values (Table. ). As a result, we obtained a much larger $\chi^{2}$ value of 8.5 for total intensity compared to our best fit. For polarimetry, $\chi^2=3.6,$ which is not far from the value we find for our best-fit model. This is because $R_0=122$ au, as used in MG18, is close to that of our polarimetric best-fit model $R_0=125$ au. We also explored models with $g_1,g_2,w_1$ as free parameters corresponding to the parameter space of Sect for total intensity and Sect for polarimetry in order to restrict possible degeneracies between dust density distribution corresponding to the ALMA values and phase function. The attempt resulted in larger $\chi^{2}$ (8.3 for total intensity and 3.6 for polarimetry) values compared to our best fit indicating that the ALMA model does not match the SPHERE image. The disagreement between the models derived from these two instruments could stem from one of two factors, or both. The first is the angular resolution as mentioned above. The second is that ALMA and SPHERE probe different grain size of which the dynamics can be governed by different processes resulting in two distinct spatial distributions.

|         Parameters         | Total intensity | Polarimetry |
|:--------------------------:|:---------------:|:-----------:|
| Inclination i ($^{\circ}$) |      88.0       |    88.5     |
|           $R_0$            |       122       |     122     |
|       $\alpha_{in}$        |        2        |     10      |
|       $\alpha_{out}$       |       -6        |     -4      |
|          $R_{in}$          |       78        |     78      |
|         $R_{out}$          |       440       |     440     |
|           $g_1$            |       0.7       |     0.8     |
|           $g_2$            |      -0.4       |     0.0     |
|           $w_1$            |      0.80       |    1.00     |
|        $h=H_{0}/R$         |      0.020      |    0.020    |
|         $\chi^{2}$         |      8.53       |    3.64     |

Parameters used to create models comparable to millimeter observations by ALMA. Here, $R_{in}$ and $R_{out}$ are the inner and outer edges of the disk

## Disk color and sub-micron grains

### Quantity of sub-micron grains produced naturally in debris discs

Our new observations and analysis confirm two striking characteristics of the HD 32297 disk: the blue color of the spectrum in the NIR and the significant presence of tiny grains much smaller than $s_{blow}$. Moreover, we confirm that there is an intrinsic coupling between these two characteristics, as was also inferred for the HD 15115 or AU Mic systems .

This link between a blue NIR spectrum and sub-micron grains has been quantitatively investigated in the recent study by . This numerical exploration shows that for bright debris disks with high fractional luminosity $f_d$, a collisional cascade at steady-state can "naturally" produce a level of unbound sub-micron grains that is high enough to lead to a blue slope of the spectrum in the NIR. This is because for bright and dense disks, the drop in grain number density at the $s= s_{blow}$ frontier, which is to a first order $\propto 1/f_d$, is much less pronounced than for fainter systems. For a very bright disk with $f_d=5\times 10^{-3}$ comparable to that of HD 32297, found a profile of the relative NIR $L_d/L*$ spectrum that is qualitatively similar to the one obtained here (see Fig. 13 of that paper). However, the blue slope they obtained is not as steep as in the present case, with a flux ratio between the $\lambda=1\muup$m and $\lambda=1.6\muup$m fluxes that is $\sim 1.1$, as compared to $\sim 1.6$ here.

This could indicate an additional source of sub-micron grains, which cannot be explained by the steady-state collisional evolution of the system. One possible cause could be the so-called collisional "avalanche" mechanism , initiated by the break-up of a large planetesimal closer to the star, which releases large amounts of unbound dust grains that then trigger a collisional chain-reaction as they sandblast at very high velocity through a dense outer disk. The ideal case for an avalanche-producing system is a double-belt configuration, with an inner belt (where the large planetesimal breaks up) at $\sim$`<!-- -->`{=html}1-10 au with $f_d\gtrsim 10^{-4}$, and a bright outer belt with $f_d\gtrsim 10^{-3}$ .This could match the structure of HD 32297, for which an inner belt of brightness $f_d\sim6\times 10^{-4}$ has been inferred by , even though the reality of this inner belt is still debated . Observational confirmation of an inner belt would need to achieve contrasts significantly higher (a factor of 10) than those currently feasible with SPHERE. In this case, the level of $s\leq s_{blow}$ grains would vary stochastically, on a timescale $t_{av}$ that is roughly a third of the typical dynamical timescale in the disk . This would, however, correspond here to $t_{av}\sim 300-400\,$yrs, much too long to be observationally monitored. Moreover, it is not guaranteed that the rate at which large planetesimals break up in the inner regions is high enough for such an event to be likely to be witnessed .

### Effect due to gas on the presence of sub-micron grains

Another possibility is that the system is able to retain $s\leq s_{blow}$ grains significantly longer than the radiation pressure blow-out time. If there is enough gas in the system, gas drag could act to significantly increase the time for an unbound grain to leave the system. To check that, we first compute the stopping time for the case where grains are bound to their host star, equal to

$$\label{eqst}
\begin{aligned}
    T_{\rm sb} \sim 2 \left( \frac{\rho}{1.5\,{{\rm g/cm}^3}} \right) \left( \frac{s}{1\, \mu{\rm m}} \right) \left( \frac{M_{\rm gas}}{0.1\, {\rm M}_\oplus} \right)^{-1},
    \end{aligned}$$

where we assume M$_\star=1.8$ M$_\odot$, $R=130$ au, $\Delta R=50$ au to be consistent with results from Sect. . We also fix the gas temperature to 30K and its mean molecular weight to 28 based on , where they show that the gas mass is dominated by CO rather than carbon in HD 32297. Accounting for the observed neutral and ionized carbon , the total gas mass barely goes above 0.1 M$_\oplus$. If accounting for potential CO$_2$ or water being released from planetesimals at the same time as CO (and thus producing extra oxygen not coming from CO and some extra hydrogen), and assuming a solar-system comet-like composition , the total gas mass could go up to 0.5 M$_\oplus$[^5].

Therefore, small bound grains close to the blow-out limit (i.e., 1-10 $\mu$m, see Table ) will have a stopping time close to 1 and will be affected by gas drag over a few orbital periods before they have time to collisionally deplete. The smallest bound grains will likely (depending on the gas-pressure gradient) move outwards before being collisionally destroyed and will therefore be present for longer than usually assumed in a standard size distribution .

For unbound grains, Eq.  needs to be adjusted. The velocity of unbound grains can reach $\sqrt{2(\beta-1)}$ of the Keplerian velocity at which they are released initially, which would increase the velocity difference between gas and dust, hence entering the Stokes regime of drag . On top of that, the stopping time of Eq.  is calculated over an orbital period and unbound grains travel almost radially over $\Delta R$, the width of the disk. Therefore, we calculate a new dimensionless Stokes stopping time for unbound grains $T_{\rm su}$ , scaled by the crossing time over $\Delta R$ and therefore equal to $T_{\rm sb} [2 c_s/v_K] [R/\Delta R$\], where $v_K$ is the Keplerian velocity and $c_s$ the sound speed:

$$\label{eqst2}
\begin{aligned}
    T_{\rm su} \sim 0.03 \left( \frac{\rho}{1.5\,{{\rm g/cm}^3}} \right) \left( \frac{s}{0.1\, \mu{\rm m}} \right) \left( \frac{M_{\rm gas}}{0.1\, {\rm M}_\oplus} \right)^{-1},
    \end{aligned}$$

meaning that 0.1 $\mu$m grains will be slowed down significantly before they have time to leave the disk[^6]. The exact orbit of unbound grains interacting with gas is time dependent and is not derived here as it goes beyond the scope of this paper, but generally speaking, an unbound grain would start on a very hyperbolic orbit and eventually be circularized around the star. These grains will then accumulate before being destroyed collisionally. Therefore, we find tantalizing evidence that gas observed in this system may be able to explain the blue color of the disk by allowing small unbound grains to be present for longer.

## Point-source detection limit

We do not detect any point source in the entire IRDIS field of view. The contrast curves at 5$\sigma$ are measured with the SpeCal pipeline for KLIP-reduced data (Fig. ). IRDIS provides contrasts of about $10^{-5}$ at 0.5$''$ and $10^{-6}$ at 1$''$. The IFS contrasts are similar or slightly better at $<0.4''$ and then degrade for larger separations. For this reason, we considered the IRDIS contrast curve only to derive the limit of detection in terms of mass. Figure displays the limit of detection using the COND model , and assuming two possible ages of the system: 10 Myr and 30 Myr. We note that if the system were 10 Myr old we should be able to detect a planet of 3 Jupiter masses at a separation of 0.5$''$ (and respectively 1.3$''$ for 30 Myr). The values below 1 M$_\mathrm{J}$ are unreliable in the framework of the COND model. Since the noise is estimated azimuthally in SpeCal, the disk itself contributes 40% to the contrast curve in the range $0.2''-1.2''$.

# Conclusions

The findings of this paper can be summarized as follows:

-   We observed the debris disk of HD 32297 in the NIR in the Y, J, and H bands out to stellocentric distances of 3.3$''$, and for the first time as close as 0.15$''$. We obtained both total intensity and polarimetric images as well.

-   At large separations, the disk is characterized by a concavity as reported by . At shorter separations ($<1''$), a bow-like shape is reminiscent of a very inclined belt of which we see mostly one side (northwest) due to the forward scattering by the grains. Noticeably, we were able to detect the back side of the disk which we modeled using two HG phase functions. This feature is not observed in polarimetric data possibly due to low S/N.

-   Upon first inspection, the disk appears to be symmetrical in NE and SW sides and has no gapped structure in contrast to the claims of . This is confirmed unambiguously with our photometric study in which the surface brightness profiles do not show any significant brightness asymmetry between the two sides, or any gap.

-   We present two methods for extracting the photometry of inclined disks observed in total intensity. The first method includes the estimation of the ADI self-subtraction in model images and accounting for this bias into the data. We find that this method can induce some irregularities depending on the accuracy of measurement of the ADI self-subtraction. In the second method we calculate a scaling factor between the ADI processed data and its corresponding model. The scaled model is used to measure the photometry instead of the data. As a drawback any departure from the model is not represented in the measurements.

-   Comparing total intensity and polarimetry in the J band we derived a polarisation fraction of about $15\%$ which is in accordance with other debris disks.

-   From photometric measurements obtained in 40 spectral channels we obtained an average spectral reflectance and conclude that the disk is 'gray to blue' color in the YJH spectral range. Using a radiative transfer module in GRaTer we were able to compare this measured reflectance with those expected for a variety of grain sizes and compositions. We found that irrespective of the composition, grains should be significantly smaller than the corresponding blowout size (sub-micron size for astrosilicates).

-   Finally, we discussed that the presence of the small grains and the associated blue color of the disk can originate from a combination of several physical processes, including steady-state collisional evolution and the avalanche process. Given the amount of gas in this system , we also found that the gas drag can retain smaller unbound grains over a longer timescale.

HD 32297 is amongst the very few known bright and extended debris disks with gas. The SPHERE observations are of unprecedented quality allowing the detection of this disk at high S/N in all spectral channels, and strong constraints to be derived on the grain properties. Confirming the trend of the spectral reflectance would require additional SPHERE observations in the K band in total intensity as well as polarimetric data in the H and K bands. It would be interesting to perform the very same type of observations and data analysis for other gas-rich debris disks and investigate if they share similarities with HD 32297 as an attempt to understand whether the presence of gas can fully explain the dust size distribution.

[^1]: Sagan NHFP Fellow

[^2]: Based on data collected at the European Southern Observatory, Chile under the programs 098.C-0686(A) and 098.C-0686(B)

[^3]: Reduced data (FITS file) of the S/N maps are only available on CDS via anonymous ftp cdsarc.u-strasbg.fr(126.96.36.199) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/

[^4]: http://sphere.osug.fr

[^5]: We note that if the gas were primordial, the total gas mass (accounting for extra H$_2$) could go up to $10^2$ M$_\oplus$, but this is probably not the case as shown in

[^6]: If grains become slow enough because of gas drag and come back to the Epstein regime, we calculate that for $\beta$ values between 1 and 5 , the unbound stopping times (over the crossing time) are roughly one to two orders of magnitude larger than derived in Eq. , meaning that even in the Epstein regime, grains of 0.1 $\mu$m will also have stopping times close to 1-10, i.e., gas will have time to substantially brake unbound grains before they leave the disk.