Researcher profile

Julien Milli

Julien Milli contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) III: The vertical structure of debris disks

Debris disks -- collisionally sustained belts of dust and sometimes gas around main sequence stars -- are remnants of planet formation processes and are found in systems ${\gtrsim}10$ Myr old. Millimeter-wavelength observations are particularly important, as the grains probed by these observations are not strongly affected by radiation pressure and stellar winds, allowing them to probe the dynamics of large bodies producing dust. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is analyzing high-resolution observations of 24 debris disks to enable the characterization of debris disk substructures across a large sample for the first time. For the most highly inclined disks, it is possible to recover the vertical structure of the disk. We aim to model and analyze the most highly inclined systems in the ARKS sample in order to uniformly extract the vertical dust distributions for a sample of well-resolved debris disks. We employed both parametric and nonparametric methods to constrain the vertical dust distributions for the most highly inclined ARKS targets. We find a broad range of aspect ratios, revealing a wide diversity in vertical structure, with a range of best-fit parametric values of $0.0026 \leq h_{\rm HWHM} \leq 0.193$ and a median best-fit value of $h_{\rm HWHM}=0.021$. The results obtained by nonparametric modeling are generally consistent with the parametric modeling results. We find that five of the 13 disks are consistent with having total disk masses less than that of Neptune (17 $M_{\oplus}$), assuming stirring by internal processes (self-stirring and collisional and frictional damping). Furthermore, most systems show a significant preference for a Lorentzian vertical profile rather than a Gaussian.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) IX: Gas-driven origin for the continuum arc in the debris disc of HD 121617

Debris discs were long considered to be largely gas-free environments governed by collisional fragmentation, gravitational stirring, and radiative forces. Recent CO detections show that gas is present, but its abundance and origin remain uncertain. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) revealed a narrow gas and dust ring in the disc HD 121617 with an asymmetric arc 40% brighter than the rest of the ring. We aim to constrain the total gas mass in HD 121617 assuming the dust arc is produced by hydrodynamical gas-dust interactions. We used the Dusty FARGO-ADSG code, modelling dust as Lagrangian particles, including radiation pressure and dust feedback, and varying the total gas mass. Simulations were compared to observations using radiative transfer. An unstable gas ring creates a size-dependent radial and azimuthal dust trap whose efficiency depends on gas mass. Two models, with 50 and 5 Earth masses of gas, reproduce both the ALMA band 7 arc and the outward offset of the VLT/SPHERE scattered-light ring via gas drag and radiation pressure. We infer a conservative gas-mass range of 2.5 to 250 Earth masses. If the ALMA asymmetry is caused by gas drag, the required gas mass compared with the observed CO implies substantial H2, consistent with primordial gas. HD 121617 would then be a hybrid disc between protoplanetary and debris stages. Since a planet could also create an arc, future observations are needed to distinguish these scenarios.

preprint2022arXiv

Polarization from Aligned Dust Grains in the $β$ Pic Debris Disk

We present 870 $μ$m ALMA polarization observations of thermal dust emission from the iconic, edge-on debris disk $β$ Pic. While the spatially resolved map does not exhibit detectable polarized dust emission, we detect polarization at the $\sim$3$σ$ level when averaging the emission across the entire disk. The corresponding polarization fraction is $P_\textrm{frac}$ = $0.51 \pm 0.19$%. The polarization position angle $χ$ is aligned with the minor axis of the disk, as expected from models of dust grains aligned via radiative alignment torques (RAT) with respect to a toroidal magnetic field ($B$-RAT) or with respect to the anisotropy in the radiation field ($k$-RAT). When averaging the polarized emission across the outer versus inner thirds of the disk, we find that the polarization arises primarily from the SW third. We perform synthetic observations assuming grain alignment via both $k$-RAT and $B$-RAT. Both models produce polarization fractions close to our observed value when the emission is averaged across the entire disk. When we average the models in the inner versus outer thirds of the disk, we find that $k$-RAT is the likely mechanism producing the polarized emission in $β$ Pic. A comparison of timescales relevant to grain alignment also yields the same conclusion. For dust grains with realistic aspect ratios (i.e., $s > 1.1$), our models imply low grain-alignment efficiencies.

preprint2022arXiv

Redundant Apodized Pupils (RAP) for high-contrast imagers robust to segmentation-due aberrations and island effects

The imaging and characterization of a larger range of exoplanets, down to young Jupiters and exo-Earths will require accessing very high contrasts at small angular separations with an increased robustness to aberrations, three constraints that drive current instrumentation development. This goal relies on efficient coronagraphs set up on extremely large diameter telescopes such as the Thirty Meter Telescope (TMT), the Giant Magellan Telescope (GMT), or the Extremely Large Telescope (ELT). However, they tend to be subject to specific aberrations that drastically deteriorate the coronagraph performance: their primary mirror segmentation implies phasing errors or even missing segments, and the size of the telescope imposes large spiders, generating low-wind effect as already observed on the Very Large Telescope (VLT)/SPHERE instrument or at the Subaru telescope, or adaptive-optics-due petaling, studied in simulations in the ELT case. The ongoing development of coronagraphs has then to take into account their sensitivity to such errors. We propose an innovative method to generate coronagraphs robust to primary mirror phasing errors and low-wind and adaptive-optics-due petaling effect. This method is based on the apodization of the segment or petal instead of the entire pupil, this apodization being then repeated to mimic the pupil redundancy. We validate this so-called Redundant Apodized Pupil (RAP) method on a James Webb Space Telescope-like pupil composed of 18 hexagonal segments segments to align, and on the VLT architecture in the case of residual low-wind effect.

preprint2022arXiv

The vertical structure of debris disks and the impact of gas

The vertical structure of debris disks provides clues about their dynamical evolution and the collision rate of the unseen planetesimals. Thanks to the ever-increasing angular resolution of contemporary instruments and facilities, we are beginning to constrain the scale height of a handful of debris disks, either at near-infrared or millimeter wavelengths. Nonetheless, this is often done for individual targets only. We present here the geometric modeling of eight disks close to edge-on, all observed with the same instrument (SPHERE) and using the same mode (dual-beam polarimetric imaging). Motivated by the presence of CO gas in two out of the eight disks, we then investigate the impact that gas can have on the scale height by performing N-body simulations including gas drag and collisions. We show that gas can quickly alter the dynamics of particles (both in the radial and vertical directions), otherwise governed by gravity and radiation pressure. We find that, in the presence of gas, particles smaller than a few tens of microns can efficiently settle toward the midplane at the same time as they migrate outward beyond the birth ring. For second generation gas ($M_\mathrm{gas} \leq 0.1$ $M_\oplus$), the vertical settling should be best observed in scattered light images compared to observations at millimeter wavelengths. But if the gas has a primordial origin ($M_\mathrm{gas} \geq 1$ $M_\oplus$), the disk will appear very flat both at near-infrared and sub-mm wavelengths. Finally, far beyond the birth ring, our results suggest that the surface brightness profile can be as shallow as $\sim -2.25$.

preprint2020arXiv

A low-mass stellar companion to the young variable star RZ Psc

RZ Psc is a young Sun-like star with a bright and warm infrared excess that is occasionally dimmed significantly by circumstellar dust structures. Optical depth arguments suggest that the dimming events do not probe a typical sight line through the circumstellar dust, and are instead caused by structures that appear above an optically thick mid-plane. This system may therefore be similar to systems where an outer disk is shadowed by material closer to the star. Here we report the discovery that RZ Psc hosts a 0.12$M_\odot$ companion at a projected separation of 23au. We conclude that the disk must orbit the primary star. While we do not detect orbital motion, comparison of the angle of linear polarization of the primary with the companion's on-sky position angle provides circumstantial evidence that the companion and disc may not share the same orbital plane. Whether the companion severely disrupts the disc, truncates it, or has little effect at all, will require further observations of both the companion and disc.

preprint2020arXiv

Increasing the raw contrast of VLT/SPHERE with the dark hole technique. I. Simulations and validation on the internal source

Context. Since 1995 and the first discovery of an exoplanet orbiting a main-sequence star, 4000 exoplanets have been discovered using several techniques. However, only a few of these exoplanets were detected through direct imaging. Indeed, the imaging of circumstellar environments requires high-contrast imaging facilities and accurate control of wavefront aberrations. Ground-based planet imagers such as VLT/SPHERE or Gemini/GPI have already demonstrated great performance. However, their limit of detection is hampered by suboptimal correction of aberrations unseen by adaptive optics (AO). Aims. Instead of focusing on the phase minimization of the pupil plane as in standard AO, we aim to directly minimize the stellar residual light in the SPHERE science camera behind the coronagraph to improve the contrast as close as possible to the inner working angle. Methods. We propose a dark hole (DH) strategy optimized for SPHERE. We used a numerical simulation to predict the global improvement of such a strategy on the overall performance of the instrument for different AO capabilities and particularly in the context of a SPHERE upgrade. Then, we tested our algorithm on the internal source with the AO in closed loop. Results. We demonstrate that our DH strategy can correct for aberrations of phase and amplitude. Moreover, this approach has the ability to strongly reduce the diffraction pattern induced by the telescope pupil and the coronagraph, unlike methods operating at the pupil plane. Our strategy enables us to reach a contrast of 5e-7 at 150 mas from the optical axis in a few minutes using the SPHERE internal source. This experiment establishes the grounds for implementing the algorithm on sky in the near future.

preprint2020arXiv

Linking studies of tiny meteoroids, zodiacal dust, cometary dust and circumstellar disks

Tiny meteoroids entering the Earth's atmosphere and inducing meteor showers have long been thought to originate partly from cometary dust. Together with other dust particles, they form a huge cloud around the Sun, the zodiacal cloud. From our previous studies of the zodiacal light, as well as other independent methods (dynamical studies, infrared observations, data related to Earth's environment), it is now established that a significant fraction of dust particles entering the Earth's atmosphere comes from Jupiter-family comets (JFCs). This paper relies on our understanding of key properties of the zodiacal cloud and of comet 67P/Churyumov-Gerasimenko, extensively studied by the Rosetta mission to a JFC. The interpretation, through numerical and experimental simulations of zodiacal light local polarimetric phase curves, has recently allowed us to establish that interplanetary dust is rich in absorbing organics and consists of fluffy particles. The ground-truth provided by Rosetta presently establishes that the cometary dust particles are rich in organic compounds and consist of quite fluffy and irregular aggregates. Our aims are as follows: (1) to make links, back in time, between peculiar micrometeorites, tiny meteoroids, interplanetary dust particles, cometary dust particles, and the early evolution of the Solar System, and (2) to show how detailed studies of such meteoroids and of cometary dust particles can improve the interpretation of observations of dust in protoplanetary and debris disks. Future modeling of dust in such disks should favor irregular porous particles instead of more conventional compact spherical particles.

preprint2020arXiv

Multiband GPI Imaging of the HR 4796A Debris Disk

We have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron.

preprint2020arXiv

Pushing Point spread function reconstruction to the next level. Application to SPHERE/ZIMPOL

Point spread function (PSF) reconstruction (PSF-R) is a well established technique to determine reliably and accurately the PSF from Adaptive Optics (AO) control loop data. We have successfully applied this technique to improve the precision on photometry and astrometry to observation of NGC6121 obtained with SPHERE/ZIMPOL as it will be presented in a forthcoming letter. Firstly, we present the methodology we followed to reconstruct the PSF combining pupil-plane and focal-plane measurements using using our PSF-R method PRIME (Beltramo-Martin et al. 2019), with upgrade of both the model and best-fitting steps compared to previous papers. Secondly, we highlight that PRIME allows to maintain the PSF fitting residual below 0.2% over 2 hours of observation and using only 30 s of AO telemetry, which may have important consequences for telemetry storage for PSF-R purpose on future 30-40 m class telescopes. Finally, we deploy PRIME in a more realistic regime using faint stars so as to identify the precision needed on the initial guess parameters to ensure the convergence towards the optimal solution.

preprint2020arXiv

The wind-driven halo in high-contrast images I: analysis from the focal plane images of SPHERE

Context. The wind driven halo is a feature observed within the images delivered by the latest generation of ground-based instruments equipped with an extreme adaptive optics system and a coronagraphic device, such as SPHERE at the VLT. This signature appears when the atmospheric turbulence conditions are varying faster than the adaptive optics loop can correct. The wind driven halo shows as a radial extension of the point spread function along a distinct direction (sometimes referred to as the butterfly pattern). When present, it significantly limits the contrast capabilities of the instrument and prevents the extraction of signals at close separation or extended signals such as circumstellar disks. This limitation is consequential because it contaminates the data a substantial fraction of the time: about 30% of the data produced by the VLT/SPHERE instrument are affected by the wind driven halo.Aims. This paper reviews the causes of the wind driven halo and presents a method to analyze its contribution directly from the scientific images. Its effect on the raw contrast and on the final contrast after post-processing is demonstrated.Methods. We used simulations and on-sky SPHERE data to verify that the parameters extracted with our method are capable of describing the wind driven halo present in the images. We studied the temporal, spatial and spectral variation of these parameters to point out its deleterious effect on the final contrast.Results. The data driven analysis we propose does provide information to accurately describe the wind driven halo contribution in the images. This analysis justifies why this is a fundamental limitation to the final contrast performance reached.Conclusions. With the established procedure, we will analyze a large sample of data delivered by SPHERE in order to propose, in the future, post-processing techniques tailored to remove the wind driven halo.