Radiative transfer is an important component of hydrodynamic simulations as it determines the thermal properties of a physical system. It is especially important in cases where heating and cooling regulate significant processes, such as in the collapse of molecular clouds, the development of gravitational instabilities in protostellar discs, disc-planet interactions, and planet migration. We compare two approximate radiative transfer methods which indirectly estimate optical depths within hydrodynamic simulations using two different metrics: (i) the gravitational potential and density of the gas (Stamatellos et al.), and (ii) the pressure scale-height (Lombardi et al.). We find that both methods are accurate for spherical configurations e.g. in collapsing molecular clouds and within clumps that form in protostellar discs. However, the pressure scale-height approach is more accurate in protostellar discs (low and high-mass discs, discs with spiral features, discs with embedded planets). We also investigate the $\beta$-cooling approximation which is commonly used when simulating protostellar discs, and in which the cooling time is proportional to the orbital period of the gas. We demonstrate that the use of a constant $\beta$ cannot capture the wide range of spatial and temporal variations of cooling in protostellar discs, which may affect the development of gravitational instabilities, planet migration, planet mass growth, and the orbital properties of planets.

The identification of on-going planet formation requires the finest angular resolutions and deepest sensitivities in observations inspired by state-of-the-art numerical simulations. Hydrodynamic simulations of planet-disk interactions predict the formation of circumplanetary disks (CPDs) around accreting planetary cores. These CPDs have eluded unequivocal detection -their identification requires predictions in CPD tracers. In this work, we aim to assess the observability of embedded CPDs with ALMA as features imprinted in the gas kinematics. We use 3D Smooth Particle Hydrodynamic (SPH) simulations of CPDs around 1 and 5 M_Jup planets at large stellocentric radii, in locally isothermal and adiabatic disks. The simulations are then connected with 3D radiative transfer for predictions in CO isotopologues. Observability is assessed by corrupting with realistic long baseline phase noise extracted from the recent HL Tau ALMA data. We find that the presence of a CPD produces distinct signposts: 1) compact emission separated in velocity from the overall circumstellar disk's Keplerian pattern, 2) a strong impact on the velocity pattern when the Doppler shifted line emission sweeps across the CPD location, and 3) a local increase in the velocity dispersion. We test our predictions with a simulation tailored for HD 100546 -which has a reported protoplanet candidate. We find that the CPDs are detectable in all 3 signposts with ALMA Cycle 3 capabilities for both 1 and 5 M_Jup protoplanets, when embedded in an isothermal disk.

The finding of residual gas in the large central cavity of the HD142527 disk motivates questions on the origin of its non-Keplerian kinematics, and possible connections with planet formation. We aim to understand the physical structure that underlies the intra-cavity gaseous flows, guided by new molecular-line data in CO(6-5) with unprecedented angular resolutions. Given the warped structure inferred from the identification of scattered-light shadows cast on the outer disk, the kinematics are consistent, to first order, with axisymmetric accretion onto the inner disk occurring at all azimuth. A steady-state accretion profile, fixed at the stellar accretion rate, explains the depth of the cavity as traced in CO isotopologues. The abrupt warp and evidence for near free-fall radial flows in HD 142527 resemble theoretical models for disk tearing, which could be driven by the reported low mass companion, whose orbit may be contained in the plane of the inner disk. The companion's high inclination with respect to the massive outer disk could drive Kozai oscillations over long time-scales; high-eccentricity periods may perhaps account for the large cavity. While shadowing by the tilted disk could imprint an azimuthal modulation in the molecular-line maps, further observations are required to ascertain the significance of azimuthal structure in the density field inside the cavity of HD142527.

We present the results of smoothed particle hydrodynamics (SPH) simulations of the disc around the young, eccentric stellar binary HD 104237. We find that the binary clears out a large cavity in the disc, driving a significant eccentricity at the cavity edge. This then precesses around the binary at a rate of $\dot{\varpi} = 0.48^{\circ}T_{\mathrm{b}}^{-1}$, which for HD 104237 corresponds to 40 years. We find that the accretion pattern into the cavity and onto the binary changes with this precession, resulting in a periodic accretion variability driven purely by the physical parameters of the binary and its orbit. For each star we find that this results in order of magnitude changes in the accretion rate. We also find that the accretion variability allows the primary to accrete gas at a higher rate than the secondary for approximately half of each precession period. Using a large number of 3-body integrations of test particles orbiting different binaries, we find good agreement between the precession rate of a test particle and our SPH disc precession. These rates also agree very well with the precession rates predicted by the analytic theory of Leung & Lee (2013), showing that their prescription can be accurately used to predict long-term accretion variability timescales for eccentric binaries accreting from a disc. We discuss the implications of our result, and suggest that this process provides a viable way of preserving unequal mass ratios in accreting eccentric binaries in both the stellar and supermassive black hole regimes

I present the results of 1D models of circumplanetary discs around planets on eccentric orbits. I use a classical viscous heating model to calculate emission fluxes at the wavelengths targeted by the NIRCam instrument on JWST, and compare the variability of this signal with the published NIRCam sensitivity specifications. This variability is theoretically detectable by JWST for a sufficiently viscous disc ($\alpha \sim 10^{-2}$) around a sufficiently eccentric planet ($e \sim 0.1-0.2$) and if the circumplanetary disc accretes material from its parent disc at a rate $\dot{M} \gtrsim 10^{-7}\, \mathrm{M}_{\odot}\,\mathrm{yr}^{-1}$. I discuss the limitations of the models used, and the implications of the result for probing the effectiveness of disc interactions for growing a planet's orbital eccentricity.

We present the results of high-resolution numerical simulations of gas clouds falling onto binary supermassive black holes to form circumbinary accretion discs, with both prograde and retrograde cloud orbits. We explore a range of clouds masses and cooling rates. We find that for low mass discs that cool fast enough to fragment, prograde discs are significantly shorter-lived than similar discs orbiting retrograde with respect to the binary. For fragmenting discs of all masses, we also find that prograde discs fragment across a narrower radial region. If the cooling is slow enough that the disc enters a self-regulating gravitoturbulent regime, we find that alignment between the disc and binary planes occurs on a timescale primarily dictated by the disc thickness. We estimate realistic cooling times for such discs, and find that in the majority of cases we expect fragmentation to occur. The longer lifetime of low-mass fragmenting retrograde discs allows them to drive significant binary evolution, and may provide a mechanism for solving the 'last parsec problem'.

The recent discovery of a number of circumbinary planets lends a new tool to astrophysicists seeking to understand how and where planet formation takes place. Of the increasingly numerous circumbinary systems, Kepler-16 is arguably the most dynamically interesting: it consists of a planet on an almost perfectly circular orbit (e = 0.0069) around a moderately eccentric binary (e = 0.16). We present high-resolution 3D smoothed-particle hydrodynamics simulations of a Kepler-16 analogue embedded in a circumbinary disc, and show that the planet's eccentricity is damped by its interaction with the protoplanetary disc. We use this to place a lower limit on the gas surface density in the real disc through which Kepler-16b migrated of \Sigma_min ~ 10 g cm^-2. This suggests that Kepler-16b, and other circumbinary planets, formed and migrated in relatively massive discs. We argue that secular evolution of circumbinary discs requires that these planets likely formed early on in the lifetime of the disc and migrated inwards before the disc lost a significant amount of its original mass.

We present high resolution 3-D simulations of the planet-disc interaction using smoothed particle hydrodynamics, to investigate the possibility of driving eccentricity growth by this mechanism. For models with a given disc viscosity (\alpha = 0.01), we find that for small planet masses (a few Jupiter masses) and canonical surface densities, no significant eccentricity growth is seen over the duration of our simulations. This contrasts with the limiting case of large planet mass (over twenty Jupiter masses) and extremely high surface densities, where we find eccentricity growth in agreement with previously published results. We identify the cause of this as being a threshold surface density for a given planet mass below which eccentricity growth cannot be excited by this method. Further, the radial profile of the disc surface density is found to have a stronger effect on eccentricity growth than previously acknowledged, implying that care must be taken when contrasting results from different disc models. We discuss the implication of this result for real planets embedded in gaseous discs, and suggest that the disc-planet interaction does not contribute significantly to observed exoplanet eccentricities.