Context. In most current debris disc models, the dynamical and the collisional evolutions are studied separately with N-body and statistical codes, respectively, because of stringent computational ...constraints. In particular, incorporating collisional effects (especially destructive collisions) into an N-body scheme has proven a very arduous task because of the exponential increase of particles it would imply. Aims. We present here LIDT-DD, the first code able to mix both approaches in a fully self-consistent way. Our aim is for it to be generic enough to be applied to any astrophysical case where we expect dynamics and collisions to be deeply interlocked with one another: planets in discs, violent massive breakups, destabilized planetesimal belts, bright exozodiacal discs, etc. Methods. The code takes its basic architecture from the LIDT3D algorithm for protoplanetary discs, but has been strongly modified and updated to handle the very constraining specificities of debris disc physics: high-velocity fragmenting collisions, radiation-pressure affected orbits, absence of gas that never relaxes initial conditions, etc. It has a 3D Lagrangian-Eulerian structure, where grains of a given size at a given location in a disc are grouped into super-particles or tracers whose orbits are evolved with an N-body code and whose mutual collisions are individually tracked and treated using a particle-in-a-box prescription designed to handle fragmenting impacts. To cope with the wide range of possible dynamics for same-sized particles at any given location in the disc, and in order not to lose important dynamical information, tracers are sorted and regrouped into dynamical families depending on their orbits. A complex reassignment routine that searches for redundant tracers in each family and reassignes them where they are needed, prevents the number of tracers from diverging. Results. The LIDT-DD code has been successfully tested on simplified cases for which robust results have been obtained in past studies: we retrieve the classical features of particle size distributions in unperturbed discs and the outer radial density profiles in ~r-1.5 outside narrow collisionally active rings as well as the depletion of small grains in dynamically cold discs. The potential of the new code is illustrated with the test case of the violent breakup of a massive planetesimal within a debris disc. Preliminary results show that we are able for the first time to quantify the timescale over which the signature of such massive break-ups can be detected. In addition to studying such violent transient events, the main potential future applications of the code are planet and disc interactions, and more generally, any configurations where dynamics and collisions are expected to be intricately connected.
Context. Debris discs are traditionally studied using two distinct types of numerical models: statistical particle-in-a-box codes to study their collisional and size distribution evolution, and ...dynamical N-body models to study their spatial structure. The absence of collisions in N-body codes is in particular a major shortcoming, as collisional processes are expected to significantly alter the results obtained from pure N-body runs. Aims. We present a new numerical model, to study the spatial structure of perturbed debris discs in both a dynamical and collisional steady-state. We focus on the competing effects of gravitational perturbations by a massive body (planet or star), the collisional production of small grains, and the radiation pressure placing these grains in possibly dynamically unstable regions. Methods. We consider a disc of parent bodies in a dynamical steady-state, from which small radiation-pressure-affected grains are released in a series of runs, each corresponding to a different orbital position of the perturber, where particles are assigned a collisional destruction probability. These collisional runs produce successive position maps that are then recombined, following a complex procedure, to generate surface density profiles for each orbital position of the perturbing body. Results. We apply our code to the case of a circumprimary disc in a binary. We find pronounced structures inside and outside the dynamical stability regions. For low eB, the disc’s structure is time varying, with spiral arms in the dynamically “forbidden” region precessing with the companion star. For high eB, the disc is strongly asymmetric but time invariant, with a pronounced density drop in the binary’s periastron direction.
Aims.
We present new optical and near-infrared images of the debris disk around the F-type star HD 114082 in the Scorpius-Centaurus OB association. We obtained direct imaging observations and ...analyzed the TESS photometric time series data of this target with the goal of searching for planetary companions to HD 114082 and characterizing the morphology of the debris disk and the scattering properties of dust particles.
Methods.
HD 114082 was observed with the VLT/SPHERE instrument in different modes – the IRDIS camera in the
K
band (2.0–2.3 μm) together with the IFS in the
Y
,
J
, and
H
bands (0.95–1.66 μm) using the angular differential imaging technique as well as IRDIS in the
H
band (1.5–1.8 μm) and ZIMPOL in the I_PRIME band (0.71–0.87 μm) using the polarimetric differential imaging technique. To constrain the basic geometrical parameters of the disk and the scattering properties of dust grains, scattered light images were fitted with a 3D model for single scattering in an optically thin dust disk using a Markov chain Monte Carlo approach. We performed aperture photometry to derive the scattering and polarized phase functions, the polarization fraction, and the spectral scattering albedo for the dust particles in the disk. This method was also used to obtain the reflectance spectrum of the disk and, in turn, to retrieve the disk color and study the dust reflectivity in comparison to the debris disk HD 117214. We also performed the modeling of the HD 114082 light curve measured by TESS using models for planet transit and stellar activity to put constraints on the radius of the detected planet and its orbit. Last, we searched for additional planets in the system by combining archival radial velocity data, astrometry, and direct imaging.
Results.
The debris disk HD 114082 appears as an axisymmetric debris belt with a radius of ~0.37″ (35 au), an inclination of ~83°, and a wide inner cavity. Dust particles in HD 114082 have a maximum polarization fraction of ~17% and a higher reflectivity when compared to the debris disk HD 117214. This high reflectivity results in a spectral scattering albedo of ~0.65 for the HD 114082 disk at near-infrared wavelengths. The disk reflectance spectrum exhibits a red color at the position of the planetesimal belt and shows no obvious features, whereas that of HD 117214 might indicate the presence of CO
2
ice. The analysis of TESS photometric data reveals a transiting planetary companion to HD 114082 with a radius of ~1
R
Jup
on an orbit with a semimajor axis of 0.7 ± 0.4 au. No additional planet was detected in the system when we combined the SPHERE images with constraints from astrometry and radial velocity. We reach deep sensitivity limits down to ~5
M
Jup
at 50 au and ~10
M
Jup
at 30 au from the central star.
Context. Planetary systems hold the imprint of the formation and of the evolution of planets especially at young ages, and in particular at the stage when the gas has dissipated leaving mostly ...secondary dust grains. The dynamical perturbation of planets in the dust distribution can be revealed with high-contrast imaging in a variety of structures. Aims. SPHERE, the high-contrast imaging device installed at the VLT, was designed to search for young giant planets in long period, but is also able to resolve fine details of planetary systems at the scale of astronomical units in the scattered-light regime. As a young and nearby star, NZ Lup was observed in the course of the SPHERE survey. A debris disk had been formerly identified with HST/NICMOS. Methods. We observed this system in the near-infrared with the camera in narrow and broad band filters and with the integral field spectrograph. High contrasts are achieved by the mean of pupil tracking combined with angular differential imaging algorithms. Results. The high angular resolution provided by SPHERE allows us to reveal a new feature in the disk which is interpreted as a superimposition of two belts of planetesimals located at stellocentric distances of ~85 and ~115 au, and with a mutual inclination of about 5°. Despite the very high inclination of the disk with respect to the line of sight, we conclude that the presence of a gap, that is, a void in the dust distribution between the belts, is likely. Conclusions. We discuss the implication of the existence of two belts and their relative inclination with respect to the presence of planets.
Context. Debris disks are the intrinsic by-products of the star and planet formation processes. Most likely due to instrumental limitations and their natural faintness, little is known about debris ...disks around low mass stars, especially when it comes to spatially resolved observations. Aims. We present new VLT/SPHERE IRDIS dual-polarization imaging (DPI) observations in which we detect the dust ring around the M2 spectral type star TWA 7. Combined with additional angular differential imaging observations we aim at a fine characterization of the debris disk and setting constraints on the presence of low-mass planets. Methods. We modeled the SPHERE DPI observations and constrain the location of the small dust grains, as well as the spectral energy distribution of the debris disk, using the results inferred from the observations, and performed simple N-body simulations. Results. We find that the dust density distribution peaks at ~0.72′′ (25 au), with a very shallow outer power-law slope, and that the disk has an inclination of ~13° with a position angle of ~91° east of north. We also report low signal-to-noise ratio detections of an outer belt at a distance of ~1.5′′ (~52 au) from the star, of a spiral arm in the southern side of the star, and of a possible dusty clump at 0.11′′. These findings seem to persist over timescales of at least a year. Using the intensity images, we do not detect any planets in the close vicinity of the star, but the sensitivity reaches Jovian planet mass upper limits. We find that the SED is best reproduced with an inner disk at ~0.2′′ (~7 au) and another belt at 0.72′′ (25 au). Conclusions. We report the detections of several unexpected features in the disk around TWA 7. A yet undetected 100 Solar Mass planet with a semi-major axis at 20−30 au could possibly explain the outer belt as well as the spiral arm. We conclude that stellar winds are unlikely to be responsible for the spiral arm.
Context. The existence of dusty debris disks around a large fraction of solar type main-sequence stars, inferred from excess far-IR and submillimetre emission compared to that expected from stellar ...photospheres, suggests that leftover planetesimal belts analogous to the asteroid- and comet reservoirs of the solar system are common. Aims. Sensitive submillimetre observations are essential to detect and characterise cold extended dust originating from collisions of small bodies in disks, belts, or rings at Kuiper-belt distances (30–50 AU or beyond). Measurements of the flux densities at these wavelengths will extend existing IR photometry and permit more detailed modelling of the Rayleigh-Jeans tail of the disks spectral energy distribution (SED), effectively constraining dust properties and disk extensions. By observing stars spanning from a few up to several hundred Myr, the evolution of debris disks during crucial phases of planet formation can be studied. Methods. We observed 22 exo-Kuiper-belt candidates at 870 μm, as part of a large programme with the LABOCA bolometer at the APEX telescope. Dust masses (or upper limits) were calculated from integrated 870 μm fluxes, and fits to the SED of detected sources revealed the fractional dust luminosities $f_{\mathrm{dust}}$, dust temperatures $T_{\mathrm{dust}}$, and power-law exponents β of the opacity law. Results. A total of 10 detections with at least 3σ significance were made, out of which five (HD 95086, HD 131835, HD 161868, HD 170773, and HD 207129) have previously never been detected at submillimetre wavelengths. Three additional sources are marginally detected with > 2.5σ significance. The best-fit β parameters all lie between 0.1 and 0.8, in agreement with previous results indicating the presence of significantly larger grains than those in the ISM. From our relatively small sample we estimate $f_{\mathrm{dust}}$ $\propto$ $t^{-\alpha}$, with α ~ 0.8–2.0, and identify an evolution of the characteristic radial dust distance $R_{\mathrm{dust}}$ that is consistent with the $t^{1/3}$ increase predicted from models of self-stirred collisions in debris disks.
Context. The vertical thickness of debris discs is often used as a measure of these systems' dynamical excitation, and as clues to the presence of hidden massive perturbers such as planetary embryos. ...However, this argument might be flawed because the observed dust should be naturally placed on inclined orbits by the combined effect of radiation pressure and mutual collisions. Aims. We critically reinvestigate this issue and numerically estimate the “natural” vertical thickness of a collisionally evolving disc, in the absence of any additional perturbing body. Methods. We use a deterministic collisional code, to follow the dynamical evolution of a population of indestructible test grains suffering mutual inelastic impacts. Grain differential sizes as well as the effect of radiation pressure are taken into account. Results. We find that, under the coupled effect of radiation pressure and collisions, grains naturally acquire inclinations of a few degrees. The disc is stratified with respect to grain sizes, the smallest grains having the largest vertical dispersion and the largest being clustered closer to the midplane. Conclusions. Debris discs should have a minimum “natural” observed aspect ratio $h_{\mathrm{min}}$ ~ 0.04±0.02 from visible to mid-IR wavelengths, where the flux is dominated by the smallest bound grains. These values are comparable to the estimated thicknesses of several vertically resolved debris discs, as illustrated by the specific example of AU Mic. For all systems with h ~ $h_{\mathrm{min}}$, the presence (or absence) of embedded perturbing bodies cannot be inferred from the vertical dispersion of the disc.
Innovative therapeutic strategies are needed to diminish the impact of harmful immunosuppression in transplantation. Dendritic cell (DC)‐based therapy is a promising approach for induction of ...antigen‐specific tolerance. Using a heart allograft model in rats, we analyzed the immunoregulatory mechanisms by which injection of autologous tolerogenic DCs (ATDCs) plus suboptimal immunosuppression promotes indefinite graft survival. Surprisingly, we determined that Interferon‐gamma (IFNG), a cytokine expected to be propathogenic, was threefold increased in the spleen of tolerant rats. Importantly, its blockade led to allograft rejection Mean Survival Time (MST) = 25.6 ± 4 days, showing that IFNG plays a critical role in immunoregulatory mechanisms triggered by ATDCs. IFNG was expressed by TCRαβ+CD3+CD4−CD8−NKRP1− cells (double negative T cells, DNT), which accumulated in the spleen of tolerant rats. Interestingly, ATDCs specifically induced IFNG production by DNT cells. ATDCs expressed the cytokinic chain Epstein‐Barr virus‐induced gene 3 (EBI3), an IL‐12 family member. EBI3 blockade or knock‐down through siRNA completely abolished IFNG expression in DNT cells. Finally, EBI3 blockade in vivo led to allograft rejection (MST = 36.8 ± 19.7 days), demonstrating for the first time a role for EBI3 in transplantation tolerance. Taken together our results have important implications in the rationalization of DC‐based therapy in transplantation as well as in the patient immunomonitoring follow‐up.
The authors describe how the injection of autologous tolerogenic dendritic cells triggers Interferon‐gamma and Epstein‐Barr virus‐induced gene 3‐dependent immunoregulatory mechanisms leading to allograft tolerance in a rat model. See editorial by Zhang and Thomson on page 2003.
There has been a renewed interest for metabolism in the computational biology community, leading to an avalanche of papers coming from methodological network analysis as well as experimental and ...theoretical biology. This paper is meant to serve as an initial guide for both the biologists interested in formal approaches and the mathematicians or computer scientists wishing to inject more realism into their models. The paper is focused on the structural aspects of metabolism only. The literature is vast enough already, and the thread through it difficult to follow even for the more experienced worker in the field. We explain methods for acquiring data and reconstructing metabolic networks, and review the various models that have been used for their structural analysis. Several concepts such as modularity are introduced, as are the controversies that have beset the field these past few years, for instance, on whether metabolic networks are small-world or scale-free, and on which model better explains the evolution of metabolism. Clarifying the work that has been done also helps in identifying open questions and in proposing relevant future directions in the field, which we do along the paper and in the conclusion.
Context. Planet perturbations have been often invoked as a potential explanation for many spatial structures that have been imaged in debris discs. So far this issue has been mostly investigated with ...pure N-body numerical models, which neglect the crucial effect collisions within the disc can have on the disc’s response to dynamical perturbations. Aims. We numerically investigate how the coupled effect of collisions and radiation pressure can affect the formation and survival of radial and azimutal structures in a disc perturbed by a planet. We consider two different set-ups: a planet embedded within an extended disc and a planet exterior to an inner debris ring. One important issue we want to address is under which conditions a planet’s signature can be observable in a collisionally active disc. Methods. We use our DyCoSS code, which is designed to investigate the structure of perturbed debris discs at dynamical and collisional steady-state, and derive synthetic images of the system in scattered light. The planet’s mass and orbit, as well as the disc’s collisional activity (parameterized by its average vertical optical depth τ0) are explored as free parameters. Results. We find that collisions always significantly damp planet-induced spatial structures. For the case of an embedded planet, the planet’s signature, mostly a density gap around its radial position, should remain detectable in head-on images if Mplanet ≥ MSaturn. If the system is seen edge-on, however, inferring the presence of the planet is much more difficult, as only weak asymmetries remain in a collisionally active disc, although some planet-induced signatures might be observable under very favourable conditions. For the case of an inner ring and an external planet, planetary perturbations cannot prevent collision-produced small fragments from populating the regions beyond the ring. The radial luminosity profile exterior to the ring is in most cases close to the one it should have in the absence of the external planet. The most significant signature left by a Jovian planet on a circular orbit are precessing azimutal structures that can be used to indirectly infer its presence. For a planet on an eccentric orbit, we show that the ring becomes elliptic and that the well known pericentre glow effect is visible despite of collisions and radiation pressure, but that detecting such features in observed discs is not an unambiguous indicator of the presence of an outer planet.