Source author record

Stefano Facchini

Stefano Facchini appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

30 published item(s)

preprint2026arXiv

Circumplanetary Disk Candidate in the Disk of HD 163296 Traced by Localized Emission from Simple Organics

Atacama Large Millimeter/submillimeter Array observations suggest that the disc of HD 163296 is being actively shaped by embedded, yet unseen protoplanets, as indicated by numerous gas and dust substructures consistent with planet-disc interaction models. We report the first detection of simple organic molecules, HCN and C2H, tracing a candidate circumplanetary disc (CPD) in the HD 163296 system, located at an orbital radius of $R=88\pm7$ au and azimuth $ϕ=46\pm3^\circ$ (or $R=0.75''$, $\rm{PA}=350^\circ$ in projected sky coordinates), and originating near the midplane of the circumstellar disc. The signature is localised but spectrally resolved, and it overlaps with a previously reported planet candidate, P94, identified through kinematic perturbations traced by CO lines. We propose a scenario in which the observed chemical anomalies arise from increased heating driven by the forming planet and ongoing accretion through its CPD, facilitating the thermal desorption of species that would otherwise remain frozen out in the disc midplane, and potentially triggering the activation barriers of chemical reactions that lead to enhanced molecular production. Based on a first-order dynamical analysis of the HCN spectrum from the CPD--isolated with a 7$σ$ significance--we infer an upper limit on the planet mass of 1.8 $M_{\rm Jup}$, consistent with predictions from CO kinematics and constraints from direct imaging studies. By comparing the CPD sizes derived from our models with theoretical expectations where the CPD radius corresponds to roughly one-third of the planet's Hill radius, we favor CPD gas temperatures $T > 150$ K, planet masses $M_{\rm p} < 1.0$ $M_{\rm Jup}$, and CPD radii $R_{\rm CPD} < 2$ au.

preprint2026arXiv

Water vapor emission at the warm cavity wall of the HD 100546 disk as revealed by ALMA

We present spatially resolved ALMA observations of the water line at 183 GHz in the disk around the Herbig star HD 100546. The water vapor emission peaks at the inner edge of the warm dust cavity, located ~15 au from the central star. We attribute this to thermal desorption at the water snowline, shifted outward at the dust cavity wall directly heated by the intense radiation. This represents the first spatially resolved image of the water snowline using ALMA observations of the main water isotopologue in a protoplanetary disk. The water emission morphology peaking inside the first dust ring is consistent with previous ALMA detections of oxygen-bearing complex organic molecules in the disk, including thermally desorbed methanol. These findings signal that warm cavities of transition disks provide ideal targets to directly reconstruct the spatial distribution of water vapor and the snowline location with ALMA, and directly connect water vapor emission to ice desorption of complex organic species.

preprint2025arXiv

Brightness variability in polar circumbinary disks

In binary systems with a strongly misaligned disk, the central binary stars can travel a significant vertical distance above and below the disk's orbital plane. This can cause large changes in illumination of the disk over the course of the binary orbital period. We use both analytic and radiative transfer models to examine the effect of changes in stellar illumination on the appearance of the disk, particularly in the case of the polar disk HD 98800B. We find that the observed flux from the disk can vary significantly over the binary orbital period, producing a periodically varying lightcurve which peaks twice each binary orbit. The amount of flux variation is strongly influenced by the disk geometry. We suggest that these flux variations produce several observable signatures, and that these observables may provide constraints on different properties of the disk such as its vertical structure, geometry, and cooling rate.

preprint2024arXiv

Mapping the Vertical Gas Structure of the Planet-hosting PDS 70 Disk

PDS 70 hosts two massive, still-accreting planets and the inclined orientation of its protoplanetary disk presents a unique opportunity to directly probe the vertical gas structure of a planet-hosting disk. Here, we use high-spatial-resolution (${\approx}$0."1;10 au) observations in a set of CO isotopologue lines and HCO$^+$ J=4-3 to map the full 2D $(r,z)$ disk structure from the disk atmosphere, as traced by $^{12}$CO, to closer to the midplane, as probed by less abundant isotopologues and HCO$^+$. In the PDS 70 disk, $^{12}$CO traces a height of $z/r\approx0.3$, $^{13}$CO is found at $z/r\approx0.1$, and C$^{18}$O originates at, or near, the midplane. The HCO$^+$ surface arises from $z/r\approx0.2$ and is one of the few non-CO emission surfaces constrained with high fidelity in disks to date. In the $^{12}$CO J=3-2 line, we resolve a vertical dip and steep rise in height at the cavity wall, making PDS 70 the first transition disk where this effect is directly seen in line emitting heights. In the outer disk, the CO emission heights of PDS 70 appear typical for its stellar mass and disk size and are not substantially altered by the two inner embedded planets. By combining CO isotopologue and HCO$^+$ lines, we derive the 2D gas temperature structure and estimate a midplane CO snowline of ${\approx}$56-85 au. This implies that both PDS 70b and 70c are located interior to the CO snowline and are likely accreting gas with a high C/O ratio of ${\approx}$1.0, which provides context for future planetary atmospheric measurements from, e.g., JWST, and for properly modeling their formation histories.

preprint2022arXiv

A new planet candidate detected in a dust gap of the disc around HD 163296 through localised kinematic signatures. An observational validation of the Discminer

We report the robust detection of coherent, localised deviations from Keplerian rotation possibly associated with the presence of two giant planets embedded in the disc around HD 163296. The analysis is performed using the DISCMINER channel-map modelling framework on $^{12}$CO $J=2-1$ DSHARP data. Not only orbital radius, but also azimuth of the planets are retrieved by our technique. One of the candidate planets, detected at $R=94\pm6$ au, $ϕ=50\pm3^\circ$ (P94), is near the centre of one of the gaps in dust continuum emission, and is consistent with a planet mass of 1 $\rm M_{\rm Jup}$. The other planet, located at $R=261\pm4$ au, $ϕ=57\pm1^\circ$ (P261), is in the region where a velocity kink was previously observed in $^{12}$CO channel maps. Also, we provide a simultaneous description of the height and temperature of the upper and lower emitting surfaces of the disc, and propose the line width as a solid observable to track gas substructure. Using azimuthally averaged line width profiles we detect gas gaps at $R=38$ au, $R=88$ au, and $R=136$ au, closely matching the location of their dust and kinematical counterparts. Furthermore, we observe strong azimuthal asymmetries in line widths around the gas gap at $R=88$ au, possibly linked to turbulent motions driven by the P94 planet. Our results confirm that the DISCMINER is capable of finding localised, otherwise unseen velocity perturbations thanks to its robust statistical framework, but also that it is well suited for studies of the gas properties and vertical structure of protoplanetary discs.

preprint2022arXiv

An APEX search for carbon emission from NGC 1977 proplyds

We used the Atacama Pathfinder Experiment (APEX) telescope to search for CI 1-0 (492.16GHz) emission towards 8 proplyds in NGC 1977, which is an FUV radiation environment two orders of magnitude weaker than that irradiating the Orion Nebular Cluster (ONC) proplyds. CI is expected to enable us to probe the wind launching region of externally photoevaporating discs. Of the 8 targets observed, no 3$σ$ detections of the CI line were made despite reaching sensitivities deeper than the anticipated requirement for detection from prior APEX CI observations of nearby discs and models of external photoevaporation of quite massive discs. By comparing both the proplyd mass loss rates and CI flux constraints with a large grid of external photoevaporation simulations, we determine that the non-detections are in fact fully consistent with the models if the proplyd discs are very low mass. Deeper observations in CI and probes of the disc mass with other tracers (e.g. in the continuum and CO) can test this. If such a test finds higher masses, this would imply carbon depletion in the outer disc, as has been proposed for other discs with surprisingly low CI fluxes, though more massive discs would also be incompatible with models that can explain the observed mass loss rates and CI non-detections. The expected remaining lifetimes of the proplyds are estimated to be similar to those of proplyds in the ONC at 0.1Myr. Rapid destruction of discs is therefore also a feature of common, intermediate UV environments.

preprint2022arXiv

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): A Panchromatic View of DO Tau's Complex Kilo-au Environment

While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales.

preprint2022arXiv

Gas and Dust Shadows in the TW Hydrae Disk

We present new observations of CO J=2-1 emission from the protoplanetary disk around TW Hya. Emission is detected out to 240 au (4") and found to exhibit azimuthal variations up to 20% beyond 180 au (3"), with the west side of the disk brighter than the east. This asymmetry is interpreted as tracing the shadow previously seen in scattered light. A reanalysis of the multi-epoch observations of the dust shadow in scattered light from Debes et al. (2017) suggests that an oscillatory motion would provide a better model of the temporal evolution of the dust shadow rather than orbital motion. Both models predict an angular offset between the dust shadow and the gas shadow of up to ~100 deg. We attribute this offset to the finite rate at which dust grains and gas molecules can exchange heat, dominated by the collisional rate between gas molecules and dust grains, $t_{\rm coll}$. The angular offsets derived are equivalent to collisional timescales that range from the near instantaneous up to $t_{\rm coll}$ ~ 10 years, depending on whether a straight or a curved dust shadow, as suggested by HST observations reported by Debes et al. (2017), is adopted. The inferred range of $t_{\rm coll}$ are consistent with those predictions based on representative gas densities, temperatures, gas-to-dust ratios and grain sizes. These results represent the first time empirical constraints can be placed on $t_{\rm coll}$.

preprint2022arXiv

High-resolution ALMA observations of transition disk candidates in Lupus

Transition disks with small inner dust cavities are interesting targets for the study of disk clearing mechanisms. Such disks have been identified through a deficit in the infrared part of their SED, but spatially resolved millimeter imaging is required to confirm the presence of an inner dust cavity. We use high-resolution ALMA observations of 30 mas resolution in Band 6 continuum and $^{12}$CO 2--1 emission of 10 transition disk candidates in the Lupus star forming region, in order to confirm the presence of inner dust cavities and infer the responsible mechanism. The continuum data are analyzed using visibility modeling and the SEDs are compared with radiative transfer models. Out of the six transition disk candidates selected from their SED, only one disk revealed an inner dust cavity of 4 au in radius. Three of the other disks are highly inclined, which limits the detectability of an inner dust cavity but it is also demonstrated to be the possible cause for the infrared deficit in their SED. The two remaining SED-selected disks are very compact, with dust radii of only $\sim$3 au. From the four candidates selected from low-resolution images, three new transition disks with large inner cavities $>$20 au are identified, bringing the total number of transition disks with large cavities in Lupus to 13. SED-selected transition disks with small cavities are biased towards highly inclined and compact disks, which casts doubt on the use of their occurrence rates in estimating dispersal timescales of photoevaporation. Using newly derived disk dust masses and radii, we re-evaluate the size-luminosity and $M_{\rm dust}-M_{\rm star}$ relations. These relations can be understood if the bright disks are dominated by disks with substructure whereas faint disks are dominated by drift-dominated disks. (Abridged)

preprint2022arXiv

Molecules with ALMA at Planet-forming Scales (MAPS). A Circumplanetary Disk Candidate in Molecular Line Emission in the AS 209 Disk

We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a $117\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using $^{12}$CO and near-infrared scattered light observations, and is associated with localized velocity perturbations in $^{12}$CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the $^{13}$CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is $\gtrsim35$ K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3σ$ flux density $\lesssim26.4~μ$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-4}$. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion.

preprint2022arXiv

The ODYSSEUS Survey. Motivation and First Results: Accretion, Ejection, and Disk Irradiation of CVSO 109

The Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program of low-mass pre-main-sequence stars, coupled with forthcoming data from ALMA and JWST, will provide the foundation to revolutionize our understanding of the relationship between young stars and their protoplanetary disks. A comprehensive evaluation of the physics of disk evolution and planet formation requires understanding the intricate relationships between mass accretion, mass outflow, and disk structure. Here we describe the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey and present initial results of the classical T Tauri Star CVSO 109 in Orion OB1b as a demonstration of the science that will result from the survey. ODYSSEUS will analyze the ULLYSES spectral database, ensuring a uniform and systematic approach in order to (1) measure how the accretion flow depends on the accretion rate and magnetic structures, (2) determine where winds and jets are launched and how mass-loss rates compare with accretion, and (3) establish the influence of FUV radiation on the chemistry of the warm inner regions of planet-forming disks. ODYSSEUS will also acquire and provide contemporaneous observations at X-ray, optical, NIR, and millimeter wavelengths to enhance the impact of the ULLYSES data. Our goal is to provide a consistent framework to accurately measure the level and evolution of mass accretion in protoplanetary disks, the properties and magnitudes of inner-disk mass loss, and the influence of UV radiation fields that determine ionization levels and drive disk chemistry.

preprint2020arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT): I. CO, CS, CN, and H2CO around DG Tau B

The chemical composition of planets is inherited by the distribution of the various molecular species in the protoplanetary disk at thetime of their formation. As of today, only a handful of disks has been imaged in multiple spectral lines with high spatial resolution. As part of a small campaign devoted to the chemical characterization of disk-outflow sources in Taurus, we report on new ALMA Band 6 observations with 20 au resolution toward the embedded young star DG Tau B. Images of the continuum emission reveals a dust disk with rings and, putatively, a leading spiral arm. The disk, as well as the prominent outflow cavities, are detected in CO, H2CO, CS, and CN while they remain undetected in SO2, HDO, and CH3OH. From the absorption of the back-side outflow, we inferred that the disk emission is optically thick in the inner 50 au. This morphology explains why no line emission is detected from this inner region and poses some limitations toward the calculation of the dust mass and the characterization of the inner gaseous disk. The H2CO and CS emission from the inner 200 au is mostly from the disk and their morphology is very similar. The CN emission significantly differs from the other two molecules as it is observed only beyond 150 au. This ring-like morphology is consistent with previous observations and the predictions of thermochemical disk models. Finally, we constrained the disk-integrated column density of all molecules. In particular, we found that the CH3OH/H2CO ratio must be smaller than 2, making the methanol non-detection still consistent with the only such a ratio available from the literature (1.27 in TW Hya).

preprint2020arXiv

Visualizing the Kinematics of Planet Formation

A stunning range of substructures in the dust of protoplanetary disks is routinely observed across a range of wavelengths. These gaps, rings and spirals are highly indicative of a population of unseen planets, hinting at the possibility of current observational facilities being able to capture planet-formation in action. Over the last decade, our understanding of the influence of a young planet on the dynamical structure of its parental disk has progressed significantly, revealing a host of potentially observable features which would betray the presence of a deeply embedded planet. In concert, recent observations have shown that subtle perturbations in the kinematic structure of protoplanetary disks are found in multiple sources, potentially the characteristic disturbances associated with embedded planets. In this work, we review the theoretical background of planet-disk interactions, focusing on the kinematical features, and the current methodologies used to observe these interactions in spatially and spectrally resolved observations. We discuss the potential pit falls of such kinematical detections of planets, providing best-practices for imaging and analysing interferometric data, along with a set of criteria to use as a benchmark for any claimed detection of embedded planets. We finish with a discussion on the current state of simulations in regard to planet-disk interactions, highlighting areas of particular interest and future directions which will provide the most significant impact in our search for embedded planets. This work is the culmination of the 'Visualizing the Kinematics of Planet Formation' workshop, held in October 2019 at the Center for Computational Astrophysics at the Flatiron Institute in New York City.

preprint2016arXiv

External photoevaporation of protoplanetary discs in sparse stellar groups: the impact of dust growth

We estimate the mass loss rates of photoevaporative winds launched from the outer edge of protoplanetary discs impinged by an ambient radiation field. We focus on mild/moderate environments (the number of stars in the group/cluster is N ~ 50), and explore disc sizes ranging between 20 and 250 AU. We evaluate the steady-state structures of the photoevaporative winds by coupling temperature estimates obtained with a PDR code with 1D radial hydrodynamical equations. We also consider the impact of dust dragging and grain growth on the final mass loss rates. We find that these winds are much more significant than have been appreciated hitherto when grain growth is included in the modelling: in particular, mass loss rates > 1e-8 M_sun/yr are predicted even for modest background field strengths ( ~ 30 G_0) in the case of discs that extend to R > 150 AU. Grain growth significantly affects the final mass loss rates by reducing the average cross section at FUV wavelengths, and thus allowing a much more vigorous flow. The radial profiles of observable quantities (in particular surface density, temperature and velocity patterns) indicate that these winds have characteristic features that are now potentially observable with ALMA. In particular, such discs should have extended gaseous emission that is dust depleted in the outer regions, characterised by a non-Keplerian rotation curve, and with a radially increasing temperature gradient.

preprint2016arXiv

Grand challenges in protoplanetary disc modelling

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

preprint2016arXiv

Photochemical-dynamical models of externally FUV irradiated protoplanetary discs

There is growing theoretical and observational evidence that protoplanetary disc evolution may be significantly affected by the canonical levels of far ultraviolet (FUV) radiation found in a star forming environment, leading to substantial stripping of material from the disc outer edge even in the absence of nearby massive stars. In this paper we perform the first full radiation hydrodynamic simulations of the flow from the outer rim of protoplanetary discs externally irradiated by such intermediate strength FUV fields, including direct modelling of the photon dominated region (PDR) which is required to accurately compute the thermal properties. We find excellent agreement between our models and the semi-analytic models of Facchini et al. (2016) for the profile of the flow itself, as well as the mass loss rate and location of their "critical radius". This both validates their results (which differed significantly from prior semi-analytic estimates) and our new numerical method, the latter of which can now be applied to elements of the problem that the semi--analytic approaches are incapable of modelling. We also obtain the composition of the flow, but given the simple geometry of our models we can only hint at some diagnostics for future observations of externally irradiated discs at this stage. We also discuss the potential for these models as benchmarks for future photochemical-dynamical codes.

preprint2016arXiv

Proplyds around a B1 star - 42 Orionis in NGC 1977

We present the discovery of seven new proplyds (i.e. sources surrounded by cometary Halpha emission characteristic of offset ionization fronts) in NGC 1977, located about 30' north of the Orion Nebula Cluster at a distance of ~400 pc. Each of these proplyds are situated at projected distances 0.04-0.27 pc from the B1V star 42 Orionis (c Ori), which is the main source of UV photons in the region. In all cases the ionization fronts of the proplyds are clearly pointing toward the common ionizing source, 42 Ori, and 6 of the 7 proplyds clearly show tails pointing away from it. These are the first proplyds to be found around a B star, with previously known examples instead being located around O stars, including those in the Orion Nebula Cluster around $θ^1$ Ori C. The radii of the offset ionization fronts in our proplyds are between ~200 and 550 AU; two objects also contain clearly resolved central sources that we associate with disks of radii 50-70 AU. The estimated strength of the FUV radiation field impinging on the proplyds is around 10-30 times less than that incident on the classic proplyds in the Orion Nebula Cluster. We show that the observed proplyd sizes are however consistent with recent models for FUV photoevaporation in relatively weak FUV radiation fields.

preprint2016arXiv

Quantum walking in curved spacetime

A discrete-time Quantum Walk (QW) is essentially a unitary operator driving the evolution of a single particle on the lattice. Some QWs admit a continuum limit, leading to familiar PDEs (e.g. the Dirac equation). In this paper, we study the continuum limit of a wide class of QWs, and show that it leads to an entire class of PDEs, encompassing the Hamiltonian form of the massive Dirac equation in $(1+1)$ curved spacetime. Therefore a certain QW, which we make explicit, provides us with a unitary discrete toy model of a test particle in curved spacetime, in spite of the fixed background lattice. Mathematically we have introduced two novel ingredients for taking the continuum limit of a QW, but which apply to any quantum cellular automata: encoding and grouping.

preprint2016arXiv

Quantum walking in curved spacetime: $(3+1)$ dimensions, and beyond

A discrete-time Quantum Walk (QW) is essentially an operator driving the evolution of a single particle on the lattice, through local unitaries. Some QWs admit a continuum limit, leading to familiar PDEs (e.g. the Dirac equation). Recently it was discovered that prior grouping and encoding allows for more general continuum limit equations (e.g. the Dirac equation in $(1+1)$ curved spacetime). In this paper, we extend these results to arbitrary space dimension and internal degree of freedom. We recover an entire class of PDEs encompassing the massive Dirac equation in $(3+1)$ curved spacetime. This means that the metric field can be represented by a field of local unitaries over a lattice.

preprint2016arXiv

The origin of the eccentricity of the hot Jupiter in CI Tau

Following the recent discovery of the first radial velocity planet in a star still possessing a protoplanetary disc (CI Tau), we examine the origin of the planet's eccentricity (e $\sim 0.3$). We show through long timescale ($10^5$ orbits) simulations that the planetary eccentricity can be pumped by the disc, even when its local surface density is well below the threshold previously derived from short timescale integrations. We show that the disc may be able to excite the planet's orbital eccentricity in $<$ a Myr for the system parameters of CI Tau. We also perform two planet scattering experiments and show that alternatively the observed planet may plausibly have acquired its eccentricity through dynamical scattering of a migrating lower mass planet, which has either been ejected from the system or swallowed by the central star. In the latter case the present location and eccentricity of the observed planet can be recovered if it was previously stalled within the disc's magnetospheric cavity.

preprint2015arXiv

A tidal encounter caught in the act: modelling a star-disc fly-by in the young RW Aurigae system

RW Aurigae (RW Aur) is a binary star system with a long molecular arm trailing the primary star. Cabrit et al. (2006) noted the resemblance between this extended structure and the tidal arm stripped from the primary star in the simulations of star-disc encounters by Clarke & Pringle (1993). In this paper we use new hydrodynamical models and synthetic observations to fit many of the parameters of RW Aur. Using hydrodynamic models we find that the morphological appearance of RW Aur can be indeed explained by a tidal encounter with the secondary star. We reproduce all the major morphological and kinematic features of the system. Using radiative transfer calculations, we find that synthetic CO and dust continuum observations of our hydrodynamic models agree well with observations. We reproduce all the main features of the line profiles, from emission fluxes to the optical depth of the different components of the system. The agreement between observations and simulations thus lends strong support to the hypothesis of a tidal encounter scenario. Finally, we propose a possible solution for the origin of the dimming of the primary star observed in 2010/2011 by Rodriguez et al. (2013).

preprint2014arXiv

Discrete Lorentz covariance for Quantum Walks and Quantum Cellular Automata

We formalize a notion of discrete Lorentz transforms for Quantum Walks (QW) and Quantum Cellular Automata (QCA), in (1 + 1)-dimensional discrete spacetime. The theory admits a diagrammatic representation in terms of a few local, circuit equivalence rules. Within this framework, we show the first-order-only covariance of the Dirac QW. We then introduce the Clock QW and the Clock QCA, and prove that they are exactly discrete Lorentz covariant. The theory also allows for non-homogeneous Lorentz transforms, between non-inertial frames.

preprint2014arXiv

Probing the presence of planets in transition discs' cavities via warps: the case of TW Hya

We are entering the era in which observations of protoplanetary discs properties can indirectly probe the presence of massive planets or low mass stellar companions interacting with the disc. In particular, the detection of warped discs can provide important clues to the properties of the star-disc system. In this paper we show how observations of warped discs can be used to infer the dynamical properties of the systems. We concentrate on circumbinary discs, where the mass of the secondary can be planetary. First, we provide some simple relations that link the amplitude of the warp in the linear regime to the parameters of the system. Secondly, we apply our method to the case of TW Hya, a transition disc for which a warp has been proposed based on spectroscopic observations. Assuming values for the disc and stellar parameters from observations, we conclude that, in order for a warp induced by a planetary companion to be detectable, the planet mass should be large ($M_{\rm p} \approx 10 - 14M_{\rm J}$) and the disc should be viscous ($α\approx 0.15 - 0.25$). We also apply our model to LkCa 15 and T Cha, where a substellar companion has been detected within the central cavity of the transition discs.

preprint2014arXiv

Quantum Circuits for the Unitary Permutation Problem

We consider the Unitary Permutation problem which consists, given $n$ unitary gates $U_1, \ldots, U_n$ and a permutation $σ$ of $\{1,\ldots, n\}$, in applying the unitary gates in the order specified by $σ$, i.e. in performing $U_{σ(n)}\ldots U_{σ(1)}$. This problem has been introduced and investigated by Colnaghi et al. where two models of computations are considered. This first is the (standard) model of query complexity: the complexity measure is the number of calls to any of the unitary gates $U_i$ in a quantum circuit which solves the problem. The second model provides quantum switches and treats unitary transformations as inputs of second order. In that case the complexity measure is the number of quantum switches. In their paper, Colnaghi et al. have shown that the problem can be solved within $n^2$ calls in the query model and $\frac{n(n-1)}2$ quantum switches in the new model. We refine these results by proving that $n\log_2(n) +Θ(n)$ quantum switches are necessary and sufficient to solve this problem, whereas $n^2-2n+4$ calls are sufficient to solve this problem in the standard quantum circuit model. We prove, with an additional assumption on the family of gates used in the circuits, that $n^2-o(n^{7/4+ε})$ queries are required, for any $ε>0$. The upper and lower bounds for the standard quantum circuit model are established by pointing out connections with the permutation as substring problem introduced by Karp.

preprint2014arXiv

The theory of globulettes: candidate precursors of brown dwarfs and free floating planets in H II regions

Large numbers of small opaque dust clouds - termed 'globulettes' by Gahm et al - have been observed in the H II regions surrounding young stellar clusters. With masses typically in the planetary (or low mass brown dwarf) regime, these objects are so numerous in some regions (e.g. the Rosette) that, if only a small fraction of them could ultimately collapse, then they would be a very significant source of free floating planets. Here we review the properties of globulettes and present a theoretical framework for their structure and evolution. We demonstrate that their interior structure is well described by a pressure confined isothermal Bonnor-Ebert sphere and that the observed mass-radius relation (mass approximately proportional to the radius squared) is a systematic consequence of a column density threshold below which components of the globulette are not identified. We also find that globulettes with this interior structure are very stable against collapse within H II regions. We follow Gahm et al in assuming that globulettes are detached from the tips of pillars protruding in from the swept up shell that borders the expanding H II region and produce a model for their dynamics, finding that globulettes will eventually impact the shell. We derive an expression for the time it takes to do so and show that dissipation of energy via dust cooling allows all globulettes to survive this encounter and escape into the wider ISM. Once there the ambient pressure drops and they disperse on timescales around 30-300 kyr and should be observable using ALMA out to distances of order a parsec. Since we find that globulettes are stable, the only route via which they might still form brown dwarfs or planets is during their collision with the shell or some other violent perturbative event.

preprint2013arXiv

Decoupled Quantum Walks, models of the Klein-Gordon and wave equations

Decoupling a vectorial PDE consists in solving the system for each component, thereby obtaining scalar PDEs that prescribe the evolution of each component independently. We present a general approach to decoupling of Quantum Walks, again defined as a procedure to obtain an evolution law for each scalar component of the QW, in such a way that it does not depend on the other components. In particular, the method is applied to show the relation between the Dirac (or Weyl) Quantum Walk in three space dimensions with (or without) mass term, and the Klein-Gordon (or wave) equation.

preprint2013arXiv

Wave-like warp propagation in circumbinary discs I. Analytic theory and numerical simulations

In this paper we analyse the propagation of warps in protostellar circumbinary discs. We use these systems as a test environment in which to study warp propagation in the bending-wave regime, with the addition of an external torque due to the binary gravitational potential. In particular, we want to test the linear regime, for which an analytic theory has been developed. In order to do so, we first compute analytically the steady state shape of an inviscid disc subject to the binary torques. The steady state tilt is a monotonically increasing function of radius. In the absence of viscosity, the disc does not present any twist. Then, we compare the time-dependent evolution of the warped disc calculated via the known linearised equations both with the analytic solutions and with full 3D numerical simulations, which have been performed with the PHANTOM SPH code using 2 million particles. We find a good agreement both in the tilt and in the phase evolution for small inclinations, even at very low viscosities. Moreover, we have verified that the linearised equations are able to reproduce the diffusive behaviour when α > H/R, where α is the disc viscosity parameter. Finally, we have used the 3D simulations to explore the non-linear regime. We observe a strongly non-linear behaviour, which leads to the breaking of the disc. Then, the inner disc starts precessing with its own precessional frequency. This behaviour has already been observed with numerical simulations in accretion discs around spinning black holes. The evolution of circumstellar accretion discs strongly depends on the warp evolution. Therefore the issue explored in this paper could be of fundamental importance in order to understand the evolution of accretion discs in crowded environments, when the gravitational interaction with other stars is highly likely, and in multiple systems.

preprint2013arXiv

Wave-like warp propagation in circumbinary discs II. Application to KH 15D

KH 15D is a protostellar binary system that shows a peculiar light curve. In order to model it, a narrow circumbinary precessing disc has been invoked, but a proper dynamical model has never been developed. In this paper, we analytically address the issue of whether such a disc can rigidly precess around KH 15D, and we relate the precessional period to the main parameters of the system. Then, we simulate the disc's dynamics by using a 1D model developed in a companion paper, such that the warp propagates into the disc as a bending wave, which is expected to be the case for protostellar discs. The validity of such an approach has been confirmed by comparing its results with full 3D SPH simulations on extended discs. In the present case, we use this 1D code to model the propagation of the warp in a narrow disc. If the inner truncation radius of the disc is set by the binary tidal torques at {\sim} 1 AU, we find that the disc should extend out to 6-10 AU (depending on the models), and is therefore wider than previously suggested. Our simulations show that such a disc does reach an almost steady state, and then precesses as a rigid body. The disc displays a very small warp, with a tilt inclination that increases with radius in order to keep the disc in equilibrium against the binary torque. However, for such wider discs, the presence of viscosity leads to a secular decay of the tilt on a timescale of {\approx} 3000 (α/0.05)^(-1) years, where α is the disc viscosity parameter. The presence of a third body (such as a planet), orbiting at roughly 10 AU might simultaneously explain the outer truncation of the disc and the maintenance of the tilt for a prolonged time.

preprint2012arXiv

Quantum computation with programmable connections between gates

A new model of quantum computation is considered, in which the connections between gates are programmed by the state of a quantum register. This new model of computation is shown to be more powerful than the usual quantum computation, e. g. in achieving the programmability of permutations of N different unitary channels with 1 use instead of N uses per channel. For this task, a new elemental resource is needed, the "quantum switch", which can be programmed to switch the order of two channels with a single use of each one.

preprint2010arXiv

No-signaling, entanglement-breaking, and localizability in bipartite channels

A bipartite quantum channel represents the interaction between systems, generally allowing for exchange of information. A special class of bipartite channels are the no-signaling ones, which do not allow communication. In Ref. [1] it has been conjectured that all no-signaling channels are mixtures of entanglement-breaking and localizable channels, which require only local operations and entanglement. Here we provide the general realization scheme, giving a counterexample to the conjecture.