Smoothed particle hydrodynamics simulations of gas and dust mixtures

We present a 'two-fluid' implementation of dust in smoothed particle hydrodynamics (SPH) in the test particle limit. The scheme is able to handle both short and long stopping times and reproduces the short friction time limit, which is not properly handled in other implementations. We apply novel tests to verify its accuracy and limitations, including multi-dimensional tests that have not been previously applied to the drag-coupled dust problem and which are particularly relevant to self-gravitating protoplanetary discs. Our tests demonstrate several key requirements for accurate simulations of gas-dust mixtures. Firstly, in standard SPH particle jitter can degrade the dust solution, even when the gas density is well reproduced. The use of integral gradients, a Wendland kernel and a large number of neighbours can control this, albeit at a greater computational cost. Secondly, when it is necessary to limit the artificial viscosity we recommend using the Cullen & Dehnen (2010) switch, since the alternative, using ${\alpha} \sim 0.1$, can generate a large velocity noise up to ${{\sigma}_v} \lesssim 0.3 c_s$ in the dust particles. Thirdly, we find that an accurate dust density estimate requires $>400$ neighbours, since, unlike the gas, the dust particles do not feel regularization forces. This density noise applies to all particle-based two-fluid implementations of dust, irrespective of the hydro solver and could lead to numerically induced fragmentation. Although our tests show accurate dusty gas simulations are possible, care must be taken to minimize the contribution from numerical noise

A dusty origin for the correlation between protoplanetary disc accretion rates and dust masses

ABSTRACT Recent observations have uncovered a correlation between the accretion rates (measured from the UV continuum excess) of protoplanetary discs and their masses inferred from observations of the submm continuum. While viscous evolution models predict such a correlation, the predicted values are in tension with data obtained from the Lupus and Upper Scorpius star-forming regions; for example, they underpredict the scatter in accretion rates, particularly in older regions. Here, we argue that since the submm observations trace the discs’ dust, by explicitly modelling the dust grain growth, evolution, and emission, we can better understand the correlation. We show that for turbulent viscosities with α ≲ 10−3, the depletion of dust from the disc due to radial drift means we can reproduce the range of masses and accretion rates seen in the Lupus and Upper Sco data sets. One consequence of this model is that the upper locus of accretion rates at a given dust mass does not evolve with the age of the region. Moreover, we find that internal photoevaporation is necessary to produce the lowest accretion rates observed. In order to replicate the correct dust masses at the time of disc dispersal, we favour relatively low photoevaporation rates ≲ 10−9 M⊙ yr−1 for most sources but cannot discriminate between EUV or X-ray-driven winds. A limited number of sources, particularly in Lupus, are shown to have higher masses than predicted by our models which may be evidence for variations in the properties of the dust or dust trapping induced in substructures.</jats:p

Evidence of a past disc-disc encounter: HV and DO Tau

© 2018 The Author(s). Published by Oxford University Press on behalf of The Royal Astronomical Society. Theory and observations suggest that star formation occurs hierarchically due to the fragmentation of giant molecular clouds. In this case we would expect substructure and enhanced stellar multiplicity in the primordial cluster. This substructure is expected to decay quickly in most environments, however historic stellar encounters might leave imprints in a protoplanetary disc (PPD) population. In a low-density environment such as Taurus, tidal tails from violent star-disc or disc-disc encounters might be preserved over timescales sufficient to be observed. In this work, we investigate the possibility that just such an event occurred between HV Tau C (itself a component of a triple system) and DO Tau ~0.1 Myr ago, as evidenced by an apparent 'bridge' structure evident in the 160 μ m emission. By modelling the encounter using smoothed particle hydrodynamics (SPH) we reproduce the main features of the observed extended structure ('V'-shaped emission pointing west of HV Tau and a tail-like structure extending east of DO Tau). We suggest that HV Tau and DO Tau formed together in a quadruple system on a scale of ~5000 au (0.025 pc)

Barriers and facilitators to health screening in men: A systematic review.

RATIONALE: Men have poorer health status and are less likely to attend health screening compared to women. OBJECTIVE: This systematic review presents current evidence on the barriers and facilitators to engaging men in health screening. METHODS: We included qualitative, quantitative and mixed-method studies identified through five electronic databases, contact with experts and reference mining. Two researchers selected and appraised the studies independently. Data extraction and synthesis were conducted using the 'best fit' framework synthesis method. RESULTS: 53 qualitative, 44 quantitative and 6 mixed-method studies were included. Factors influencing health screening uptake in men can be categorized into five domains: individual, social, health system, healthcare professional and screening procedure. The most commonly reported barriers are fear of getting the disease and low risk perception; for facilitators, they are perceived risk and benefits of screening. Male-dominant barriers include heterosexual -self-presentation, avoidance of femininity and lack of time. The partner's role is the most common male-dominant facilitator to screening. CONCLUSIONS: This systematic review provides a comprehensive overview of barriers and facilitators to health screening in men including the male-dominant factors. The findings are particularly useful for clinicians, researchers and policy makers who are developing interventions and policies to increase screening uptake in men

A super-resolution analysis of the DSHARP survey: Substructure is common in the inner 30 au

The DSHARP survey evidenced the ubiquity of substructure in the mm dust distribution of large, bright protoplanetary discs. Intriguingly, these datasets have yet higher resolution information that is not recovered in a CLEAN image. We first show that the intrinsic performance of the CLEAN algorithm is resolution-limited. Then analyzing all 20 DSHARP sources using the 1D, super-resolution code Frankenstein (frank), we accurately fit the 1D visibilities to a mean factor of 4.3 longer baseline than the Fourier transform of the CLEAN images and a factor of 3.0 longer baseline than the transform of the CLEAN component models. This yields a higher resolution brightness profile for each source, identifying new substructure interior to 30 au in multiple discs; resolving known gaps to be deeper, wider, and more structured; and known rings to be narrower and brighter. Across the survey, high contrast gaps are an average 14% wider and 44% deeper in the frank profiles relative to CLEAN, and high contrast rings are an average 26% narrower. Categorizing the frank brightness profiles into trends, we find that the relative scarcity of features interior to 30 au in the survey's CLEAN images is an artifact of resolving power, rather than an intrinsic rarity of inner disc (or compact disc) substructure. Finally the rings in the frank profiles are narrower than the previously inferred deconvolved widths, indicating smaller alpha / St ratios in the local gas disc

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.This work has been supported by the DISCSIM project, grant agreement 341137 funded by the European Research Council under ERC-2013-ADG, and from STFC through grant ST/L000636/1. This work used the DIRAC Shared Memory Processing and Data Analytic systems, both at the University of Cambridge and operated, respectively, by the COSMOS Project at the Department of Applied Mathematics and Theoretical Physics and the Cambridge High Performance Computing Service, on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk). This equipment was funded by BIS National E-infrastructure capital grants ST/J005673/1 and ST/K001590/1, STFC capital grants ST/H008586/1, ST/H008861/1 and ST/H00887X/1, STFC DiRAC Operations grant ST/K00333X/1, and STFC DiRAC Operations grant ST/K00333X/1. DiRAC is part of the National E-Infrastructure.This is the final version of the article. It first appeared from Oxford University Press via https://doi.org/10.1093/mnrasl/slw18

On the Origin of the Spiral Morphology in the Elias 2-27 Circumstellar Disk

The young star Elias 2-27 has recently been observed to posses a massive circumstellar disk with two prominent large-scale spiral arms. In this Letter, we perform three-dimensional Smoothed Particle Hydrodynamics simulations, radiative transfer modeling, synthetic ALMA imaging, and an unsharped masking technique to explore three possibilities for the origin of the observed structures - an undetected companion either internal or external to the spirals, and a self-gravitating disk. We find that a gravitationally unstable disk and a disk with an external companion can produce morphology that is consistent with the observations. In addition, for the latter, we find that the companion could be a relatively massive planetary-mass companion (≲10-13 M Jup ) and located at large radial distances (between ≈300-700 au). We therefore suggest that Elias 2-27 may be one of the first detections of a disk undergoing gravitational instabilities, or a disk that has recently undergone fragmentation to produce a massive companion.We acknowledge support from the DISCSIM project, grant agreement 341137 under ERC-2013-ADG. F.M. acknowledges support from The Leverhulme Trust. This Letter uses the following ALMA data: ADS/JAO.ALMA# 2013.1.00498.S. This work used the Darwin DiRAC HPC cluster at the University of Cambridge and was undertaken on the Cambridge COSMOS SMP system, part of the STFC DiRAC HPC Facility supported by BIS NeI capital grant ST/J005673/1 and STFC grants ST/H008586/1, ST/K00333X/1

Eccentricity evolution during planet-disc interaction

During the process of planet formation, the planet-disc interactions might excite (or damp) the orbital eccentricity of the planet. In this paper, we present two long (t ~ 3 × 10 5 orbits) numerical simulations: (a) one (with a relatively light disc, M d /M p = 0.2), where the eccentricity initially stalls before growing at later times and (b) one (with a more massive disc, M d /M p = 0.65) with fast growth and a late decrease of the eccentricity. We recover the well-known result that a more massive disc promotes a faster initial growth of the planet eccentricity. However, at late times the planet eccentricity decreases in the massive disc case, but increases in the light disc case. Both simulations show periodic eccentricity oscillations superimposed on a growing/decreasing trend and a rapid transition between fast and slow pericentre precession. The peculiar and contrasting evolution of the eccentricity of both planet and disc in the two simulations can be understood by invoking a simple toy model where the disc is treated as a second point-like gravitating body, subject to secular planet-planet interaction and eccentricity pumping/damping provided by the disc. We show how the counterintuitive result that the more massive simulation produces a lower planet eccentricity at late times can be understood in terms of the different ratios of the disc-to-planet angular momentum in the two simulations. In our interpretation, at late times the planet eccentricity can increase more in low-mass discs rather than in high-mass discs, contrary to previous claims in the literature.This work has been supported by the DISCSIM project, grant agreement 341137 funded by the European Research Council under ERC-2013-ADG. JT acknowledges support from STFC through grant ST/L000636/1. This work used the Wilkes GPU cluster at the University of Cambridge High Performance Computing Service (http://www.hpc.cam.ac.uk/), provided by Dell Inc., NVIDIA and Mellanox, and part funded by STFC with industrial sponsorship from Rolls Royce and Mitsubishi Heavy Industries. We also thank the MIAPP for hosting us for the ‘Protoplanetary Disks and Planet Formation and Evolution’ topical workshop held in Munich during June 2017. All the figures were generated with the python-based package matplotlib (Hunter 2007)

The extremely truncated circumstellar disc of V410 X-ray 1: A precursor to TRAPPIST-1?

Protoplanetary discs around brown dwarfs and very low mass stars offer some of the best prospects for forming Earth-sized planets in their habitable zones. To this end, we study the nature of the disc around the very low mass star V410 X-ray 1, whose SED is indicative of an optically thick and very truncated dust disc, with our modelling suggesting an outer radius of only 0.6 au. We investigate two scenarios that could lead to such a truncation, and find that the observed SED is compatible with both. The first scenario involves the truncation of both the dust and gas in the disc, perhaps due to a previous dynamical interaction or the presence of an undetected companion. The second scenario involves the fact that a radial location of 0.6 au is close to the expected location of the H$_2$O snowline in the disc. As such, a combination of efficient dust growth, radial migration, and subsequent fragmentation within the snowline leads to an optically thick inner dust disc and larger, optically thin outer dust disc. We find that a firm measurement of the CO $J=2$--1 line flux would enable us to distinguish between these two scenarios, by enabling a measurement of the radial extent of gas in the disc. Many models we consider contain at least several Earth-masses of dust interior to 0.6 au, suggesting that V410 X-ray 1 could be a precursor to a system with tightly-packed inner planets, such as TRAPPIST-1

Severe Epstein-Barr virus infection in primary immunodeficiency and the normal host

Epstein–Barr virus (EBV) infection is ubiquitous in humans, but the majority of infections have an asymptomatic or self-limiting clinical course. Rarely, individuals may develop a pathological EBV infection with a variety of life threatening complications (including haemophagocytosis and malignancy) and others develop asymptomatic chronic EBV viraemia. Although an impaired ability to control EBV infection has long been recognised as a hallmark of severe T-cell immunodeficiency, the advent of next generation sequencing has identified a series of Primary Immunodeficiencies in which EBV-related pathology is the dominant feature. Chronic active EBV infection is defined as chronic EBV viraemia associated with systemic lymphoproliferative disease, in the absence of immunodeficiency. Descriptions of larger cohorts of patients with chronic active EBV in recent years have significantly advanced our understanding of this clinical syndrome. In this review we summarise the current understanding of the pathophysiology and natural history of these diseases and clinical syndromes, and discuss approaches to the investigation and treatment of severe or atypical EBV infection