We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates.

We construct a SCUBA-2 450-$\mu$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$\sigma$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $\mu$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$\mu$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$\beta_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$\mu$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$.

The streaming instability is a popular candidate for planetesimal formation by concentrating dust particles to trigger gravitational collapse. However, its robustness against physical conditions expected in protoplanetary disks is unclear. In particular, particle stirring by turbulence may impede the instability. To quantify this effect, we develop the linear theory of the streaming instability with external turbulence modelled by gas viscosity and particle diffusion. We find the streaming instability is sensitive to turbulence, with growth rates becoming negligible for alpha-viscosity parameters $\alpha \gtrsim \mathrm{St} ^{1.5}$, where $\mathrm{St}$ is the particle Stokes number. We explore the effect of non-linear drag laws, which may be applicable to porous dust particles, and find growth rates are modestly reduced. We also find that gas compressibility increase growth rates by reducing the effect of diffusion. We then apply linear theory to global models of viscous protoplanetary disks. For minimum-mass Solar nebula disk models, we find the streaming instability only grows within disk lifetimes beyond $\sim 10$s of AU, even for cm-sized particles and weak turbulence ($\alpha\sim 10^{-4}$). Our results suggest it is rather difficult to trigger the streaming instability in non-laminar protoplanetary disks, especially for small particles.

The $M_{BH}$ - $\sigma_{\star}$ relation is considered a result of co-evolution between the host galaxies and their super-massive black holes. For elliptical-bulge hosting inactive galaxies, this relation is well established, but there is still a debate whether active galaxies follow the same relation. In this paper, we estimate black hole masses for a sample of 19 local luminous AGNs (LLAMA) in order to test their location on the $M_{BH}$ - $\sigma_{\star}$ relation. Super-massive black hole masses ($M_{BH}$) were derived from the broad-line based relations for H$\alpha$, H$\beta$ and Pa$\beta$ emission line profiles for the Type 1 AGNs. We compare the bulge stellar velocity dispersion ($\sigma_{\star}$) as determined from the Ca II triplet (CaT) with the dispersion measured from the near-infrared CO (2-0) absorption features for each AGN and find them to be consistent with each other. We apply an extinction correction to the observed broad line fluxes and we correct the stellar velocity dispersion by an average rotation contribution as determined from spatially resolved stellar kinematic maps. The H$\alpha$-based black hole masses of our sample of AGNs were estimated in the range 6.34 $\leq$ $\log{M_{BH}}$ $\leq$ 7.75 M$_\odot$ and the $\sigma_{\star CaT}$ estimates range between 73 $\leq$ $\sigma_{\star CaT}$ $\leq$ 227 km s$^{-1}$. From the so-constructed $M_{BH}$ - $\sigma_{\star}$ relation for our Type 1 AGNs, we estimate the black hole masses for the Type 2 AGNs and the inactive galaxies in our sample. In conclusion, we find that our sample of local luminous AGNs is consistent with the $M_{BH}$ - $\sigma_{\star}$ relation of lower luminosity AGNs and inactive galaxies, after correcting for dust extinction and the rotational contribution to the stellar velocity dispersion.

We discuss the ability of a dark fluid becoming relevant around the time of matter radiation equality to significantly relieve the tension between local measurements of the Hubble constant and CMB inference, within the $\Lambda$CDM model. We show that the gravitational impact of acoustic oscillations in the dark fluid balance the effects on the CMB and result in an improved fit to CMB measurements themselves while simultaneously raising the Hubble constant. The required balance favors a model where the fluid is a scalar field that converts its potential to kinetic energy around matter radiation equality which then quickly redshifts away. We derive the requirements on the potential for this conversion mechanism and find that a simple canonical scalar with two free parameters for its local slope and amplitude robustly improves the fit to the combined data by $\Delta\chi^2 \approx 12.7$ over $\Lambda$CDM. We uncover the CMB polarization signatures that can definitively test this scenario with future data.

We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar bar which is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with HII regions. We model the kinematics of the ring using the velocity field of the CO (2--1) emission and stars and find it is consistent with a rotating disk. The outer regions of the disk, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$_{\sun}$ yr$^{-1}$. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disk, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$_{\sun}$ yr$^{-1}$ consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region.

We discuss the difference between various gauge-invariant quantities typically used in single-field inflation, namely synchronous $\zeta_s$, comoving $\zeta_c$, and unitary $\zeta_u$ curvatures. We show that conservation of $\zeta_c$ outside the horizon is quite restrictive on models as it leads to conservation of $\zeta_s$ and $\zeta_u$, whereas the reverse does not hold. We illustrate the consequence of these differences with two inflationary models: ultra-slow-roll (USR) and braiding-ultra-slow-roll (BUSR). In USR, we show that out of the three curvatures, only $\zeta_s$ is conserved outside the horizon, and we connect this result to the concepts of separate universe and the usage of the $\delta N$ formalism. We find that even though $\zeta_s$ is conserved, there is still a mild violation of the separate universe approximation in the continuity equation. Nevertheless, the $\delta N$ formalism can still be applied to calculate the primordial power spectrum of some gauge-invariant quantities such as $\zeta_u$, although it breaks down for others such as the uniform-density curvature. In BUSR, we show that both $\zeta_u$ and $\zeta_s$ are conserved outside the horizon, but take different values. Additionally, since $\zeta_u\not=\zeta_c$ we find that the prediction for observable curvature fluctuations after inflation does not reflect $\zeta_c$ at horizon crossing during inflation and moreover involves not just $\zeta_u$ at that epoch but also the manner in which the braiding phase ends.

In the innermost regions of protoplanerary discs, the solid-to-gas ratio can be increased considerably by a number of processes, including photoevaporative and particle drift. MHD disc models also suggest the existence of a dead-zone at $R\lesssim 10$ AU, where the regions close to the midplane remain laminar. In this context, we use two-fluid hydrodynamical simulations to study the interaction between a low-mass planet ($\sim 1.7 \;{\rm M_\oplus}$) on a fixed orbit and an inviscid pebble-rich disc with solid-to-gas ratio $\epsilon\ge 0.5$. For pebbles with Stokes numbers St=0.1, 0.5, multiple dusty vortices are formed through the Rossby Wave Instability at the planet separatrix. Effects due to gas drag then lead to a strong enhancement in the solid-to-gas ratio, which can increase by a factor of $\sim 10^3$ for marginally coupled particles with St=0.5. As in streaming instabilities, pebble clumps reorganize into filaments that may plausibly collapse to form planetesimals. When the planet is allowed to migrate in a MMSN disc, the vortex instability is delayed due to migration but sets in once inward migration stops due a strong positive pebble torque. Again, particle filaments evolving in a gap are formed in the disc while the planet undergoes an episode of outward migration. Our results suggest that vortex instabilities triggered by low-mass planets could play an important role in forming planetesimals in pebble-rich, inviscid discs, and may significantly modify the migration of low-mass planets. They also imply that planetary dust gaps may not necessarily contain planets if these migrated away.

Enhancing the local dust-to-gas ratio in protoplanetary discs is a necessary first step to planetesimal formation. In laminar discs, dust settling is an efficient mechanism to raise the dust-to-gas ratio at the disc midplane. However, turbulence, if present, can stir and lift dust particles, which ultimately hinders planetesimal formation. In this work, we study dust settling in protoplanetary discs with hydrodynamic turbulence sustained by the vertical shear instability. We perform axisymmetric numerical simulations to investigate the effect of turbulence, particle size, and solid abundance or metallicity on dust settling. We highlight the positive role of drag forces exerted onto the gas by the dust for settling to overcome the vertical shear instability. In typical disc models we find particles with a Stokes number $\sim 10^{-3}$ can sediment to $\lesssim 10\%$ of the gas scale-height, provided that $\Sigma_\mathrm{d}/\Sigma_\mathrm{g}\gtrsim 0.02$-$0.05$, where $\Sigma_\mathrm{d,g}$ are the surface densities in dust and gas, respectively. This coincides with the metallicity condition for small particles to undergo clumping via the streaming instability. Super-solar metallicities, at least locally, are thus required for a self-consistent picture of planetesimal formation. Our results also imply that dust rings observed in protoplanetary discs should have smaller scale-heights than dust gaps, provided that the metallicity contrast between rings and gaps exceed the corresponding contrast in gas density.

Gap-opening planets can generate dust-trapping vortices that may explain some of the latest discoveries of high-contrast crescent-shaped dust asymmetries in transition discs. While planet-induced vortices were previously thought to have concentrated shapes, recent computational work has shown that these features naturally become much more elongated in the gas when simulations account for the relatively long timescale over which planets accrete their mass. In this work, we conduct two-fluid hydrodynamical simulations of vortices induced by slowly-growing Jupiter-mass planets in discs with very low viscosity ($\alpha = 3 \times 10^{-5}$). We simulate the dust dynamics for four particle sizes spanning 0.3 mm to 1 cm in order to produce synthetic ALMA images. In our simulations, we find that an elongated vortex still traps dust, but not directly at its center. With a flatter pressure bump and disruptions from the planet's overlapping spiral density waves, the dust instead circulates around the vortex. This motion (1) typically carries the peak off-center, (2) spreads the dust out over a wider azimuthal extent $\geq 180^{\circ}$, (3) skews the azimuthal profile towards the front of the vortex, and (4) can also create double peaks in newly-formed vortices. In particular, we expect that the most defining observational signature, a peak offset of more than $30^{\circ}$, should be detectable $>30\%$ of the time in observations with a beam diameter of at most the planet's separation from its star.

We discuss the phenomenological imprints of modifications to gravity in the early universe with a specific focus on the time of recombination. We derive several interesting results regarding the effect that such modifications have on cosmological observables, especially on the driving and phasing of acoustic oscillations, observed in the CMB and BAO, as well as the weak gravitational lensing of the CMB and of galaxy shapes. This widens the pool of measurements that can be used to test gravity with present and future surveys, in particular realizing the full constraining power of the structure of the primary peaks of the CMB spectrum. We investigate whether such a phenomenology can relax tensions between cosmological measurements and find that a modification of the gravitational constant at recombination would help in reconciling measurements of the CMB with local measurements of the Hubble constant.

Large-scale, dust-trapping vortices may account for observations of asymmetric protoplanetary discs. Disc vortices are also potential sites for accelerated planetesimal formation by concentrating dust grains. However, in 3D discs vortices are subject to destructive `elliptic instabilities', which reduces their viability as dust traps. The survival of vortices in 3D accretion discs is thus an important issue to address. In this work, we perform shearing box simulations to show that disc self-gravity enhances the survival of 3D vortices, even when self-gravity is weak in the classic sense (e.g. with a Toomre $Q\simeq5$). We find a 3D, self-gravitating vortex can grow on secular timescales in spite of the elliptic instability. The vortex aspect-ratio decreases as it strengthens, which feeds the elliptic instability. The result is a 3D vortex with a turbulent core that persists for $\sim 10^{3}$ orbits. We find when gravitational and hydrodynamic stresses become comparable, the vortex may undergo episodic bursts, which we interpret as interaction between elliptic and gravitational instabilities. We estimate the distribution of dust particles in self-gravitating, turbulent vortices. Our results suggest large-scale vortices in protoplanetary discs are more easily observed at large radii.

It is expected that a pressure bump can be formed at the inner edge of a dead-zone, and where vortices can develop through the Rossby Wave Instability (RWI). It has been suggested that self-gravity can significantly affect the evolution of such vortices. We present the results of 2D hydrodynamical simulations of the evolution of vortices forming at a pressure bump in self-gravitating discs with Toomre parameter in the range $4-30$. We consider isothermal plus non-isothermal disc models that employ either the classical $\beta$ prescription or a more realistic treatment for cooling. The main aim is to investigate whether the condensating effect of self-gravity can stabilize vortices in sufficiently massive discs. We confirm that in isothermal disc models with ${\cal Q} \gtrsim 15$, vortex decay occurs due to the vortex self-gravitational torque. For discs with $3\lesssim {\cal Q} \lesssim 7$, the vortex develops gravitational instabilities within its core and undergoes gravitational collapse, whereas more massive discs give rise to the formation of global eccentric modes. In non-isothermal discs with $\beta$ cooling, the vortex maintains a turbulent core prior to undergoing gravitational collapse for $\beta \lesssim 0.1$, whereas it decays if $\beta \ge 1$. In models that incorpore both self-gravity and a better treatment for cooling, however, a stable vortex is formed with aspect ratio $\chi \sim 3-4$. Our results indicate that self-gravity significantly impacts the evolution of vortices forming in protoplanetary discs, although the thermodynamical structure of the vortex is equally important for determining its long-term dynamics.

Protoplanets may be born into dust-rich environments if planetesimals formed through streaming or gravitational instabilities, or if the protoplanetary disc is undergoing mass loss due to disc winds or photoevaporation. Motivated by this possibility, we explore the interaction between low mass planets and dusty protoplanetary discs with focus on disc-planet torques. We implement Lin & Youdin's newly developed, purely hydrodynamic model of dusty gas into the PLUTO code to simulate dusty protoplanetary discs with an embedded planet. We find that for imperfectly coupled dust and high metallicity, e.g. Stokes number $10^{-3}$ and dust-to-gas ratio $\Sigma_\mathrm{d}/\Sigma_\mathrm{g}=0.5$ , a `bubble' develops inside the planet's co-orbital region, which introduces unsteadiness in the flow. The resulting disc-planet torques sustain large amplitude oscillations that persists well beyond that in simulations with perfectly coupled dust or low dust-loading, where co-rotation torques are always damped. We show that the desaturation of the co-rotation torques by finite-sized particles is related to potential vorticity generation from the misalignment of dust and gas densities. We briefly discuss possible implications for the orbital evolution of protoplanets in dust-rich discs. We also demonstrate Lin & Youdin's dust-free framework reproduces previous results pertaining to dusty protoplanetary discs, including dust-trapping by pressure bumps, dust settling, and the streaming instability.

We explore the relationship between X-ray absorption and optical obscuration within the BAT AGN Spectroscopic Survey (BASS) which has been collecting and analyzing the optical and X-ray spectra for 641 hard X-ray selected ($E>14$ keV) active galactic nuclei (AGN). We use the deviation from a linear broad H$\alpha$-to-X-ray relationship as an estimate of the maximum optical obscuration towards the broad line region and compare the $A_{\rm V}$ to the hydrogen column densities ($N_{\rm H}$) found through systematic modeling of their X-ray spectra. We find that the inferred columns implied by $A_{\rm V}$ towards the broad line region (BLR) are often orders of magnitude less than the columns measured towards the X-ray emitting region indicating a small scale origin for the X-ray absorbing gas. After removing 30\% of Sy 1.9s that potentially have been misclassified due to outflows, we find that 86\% (164/190) of the Type 1 population (Sy 1--1.9) are X-ray unabsorbed as expected based on a single obscuring structure. However, 14\% (26/190), of which 70\% (18/26) are classified as Sy 1.9, are X-ray absorbed, suggesting the broad line region itself is providing extra obscuration towards the X-ray corona. The fraction of X-ray absorbed Type 1 AGN remains relatively constant with AGN luminosity and Eddington ratio, indicating a stable broad line region covering fraction.

Using new APEX and JCMT spectroscopy of the CO 2-1 line, we undertake a controlled study of cold molecular gas in moderately luminous Active Galactic Nuclei (AGN) and inactive galaxies from the Luminous Local AGN with Matched Analogs (LLAMA) survey. We use spatially resolved infrared photometry of the LLAMA galaxies from 2MASS, WISE, IRAS & Herschel, corrected for nuclear emission using multi-component spectral energy distribution (SED) fits, to examine the dust-reprocessed star-formation rates (SFRs), molecular gas fractions and star formation efficiencies (SFEs) over their central 1 - 3 kpc. We find that the gas fractions and central SFEs of both active and inactive galaxies are similar when controlling for host stellar mass and morphology (Hubble type). The equivalent central molecular gas depletion times are consistent with the discs of normal spiral galaxies in the local Universe. Despite energetic arguments that the AGN in LLAMA should be capable of disrupting the observable cold molecular gas in their central environments, our results indicate that nuclear radiation only couples weakly with this phase. We find a mild preference for obscured AGN to contain higher amounts of central molecular gas, which suggests a connection between AGN obscuration and the gaseous environment of the nucleus. Systems with depressed SFEs are not found among the LLAMA AGN. We speculate that the processes that sustain the collapse of molecular gas into dense pre-stellar cores may also be a prerequisite for the inflow of material on to AGN accretion disks.

In a complete sample of local 14-195 keV selected AGNs and inactive galaxies, matched by their host galaxy properties, we study the spatially resolved stellar kinematics and luminosity distributions at near-infrared wavelengths on scales of 10-150 pc, using SINFONI on the VLT. In this paper, we present the first half of the sample, which comprises 13 galaxies, 8 AGNs and 5 inactive galaxies. The stellar velocity fields show a disk-like rotating pattern, for which the kinematic position angle is in agreement with the photometric position angle obtained from large scale images. For this set of galaxies, the stellar surface brightness of the inactive galaxy sample is generally comparable to the matched sample of AGN but extends to lower surface brightness. After removal of the bulge contribution, we find a nuclear stellar light excess with an extended nuclear disk structure, and which exhibits a size-luminosity relation. While we expect the excess luminosity to be associated with a dynamically cooler young stellar population, we do not typically see a matching drop in dispersion. This may be because these galaxies have pseudo-bulges in which the intrinsic dispersion increases towards the centre. And although the young stars may have an impact in the observed kinematics, their fraction is too small to dominate over the bulge and compensate the increase in dispersion at small radii, so no dispersion drop is seen. Finally, we find no evidence for a difference in the stellar kinematics and nuclear stellar luminosity excess between these active and inactive galaxies.

Small solids embedded in gaseous protoplanetary disks are subject to strong dust-gas friction. Consequently, tightly-coupled dust particles almost follow the gas flow. This near conservation of dust-to-gas ratio along streamlines is analogous to the near conservation of entropy along flows of (dust-free) gas with weak heating and cooling. We develop this thermodynamic analogy into a framework to study dusty gas dynamics in protoplanetary disks. We show that an isothermal dusty gas behaves like an adiabatic pure gas; and that finite dust-gas coupling may be regarded as an effective heating/cooling. We exploit this correspondence to deduce that 1) perfectly coupled, thin dust layers cannot cause axisymmetric instabilities; 2) radial dust edges are unstable if the dust is vertically well-mixed; 3) the streaming instability necessarily involves a gas pressure response that lags behind dust density; 4) dust-loading introduces buoyancy forces that generally stabilizes the vertical shear instability associated with global radial temperature gradients. We also discuss dusty analogs of other hydrodynamic processes (e.g. Rossby wave instability, convective overstability, and zombie vortices), and how to simulate dusty protoplanetary disks with minor tweaks to existing codes for pure gas dynamics.

The presence of a giant planet in a low-viscosity disc can create a gap edge in the disc's radial density profile sharp enough to excite the Rossby Wave Instability. This instability may evolve into dust-trapping vortices that might explain the "banana-shaped" features in recently observed asymmetric transition discs with inner cavities. Previous hydrodynamical simulations of planet-induced vortices have neglected the timescale of hundreds to thousands of orbits to grow a massive planet to Jupiter-size. In this work, we study the effect of a giant planet's runaway growth timescale on the lifetime and characteristics of the resulting vortex. For two different planet masses (1 and 5 Jupiter masses) and two different disc viscosities ($\alpha$=3$\times 10^{-4}$ and 3$\times10^{-5}$), we compare the vortices induced by planets with several different growth timescales between 10 and 4000 planet orbits. In general, we find that slowly-growing planets create significantly weaker vortices with lifetimes and surface densities reduced by more than $50\%$. For the higher disc viscosity, the longest growth timescales in our study inhibit vortex formation altogether. Additionally, slowly-growing planets produce vortices that are up to twice as elongated, with azimuthal extents well above $180^{\circ}$ in some cases. These unique, elongated vortices likely create a distinct signature in the dust observations that differentiates them from the more concentrated vortices that correspond to planets with faster growth timescales. Lastly, we find that the low viscosities necessary for vortex formation likely prevent planets from growing quickly enough to trigger the instability in self-consistent models.

We discuss the environment of local hard X-ray selected active galaxies, with reference to two independent group catalogues. We find that the fraction of these AGN in S0 host galaxies decreases strongly as a function of galaxy group size (halo mass) - which contrasts with the increasing fraction of galaxies of S0 type in denser environments. However, there is no evidence for an environmental dependence of AGN in spiral galaxies. Because most AGN are found in spiral galaxies, this dilutes the signature of environmental dependence for the population as a whole. We argue that the differing results for AGN in disk-dominated and bulge-dominated galaxies is related to the source of the gas fuelling the AGN, and so may also impact the luminosity function, duty cycle, and obscuration. We find that there is a significant difference in the luminosity function for AGN in spiral and S0 galaxies, and tentative evidence for some difference in the fraction of obscured AGN.

In recent years, direct detection, indirect detection, and collider experiments have placed increasingly stringent constraints on particle dark matter, exploring much of the parameter space associated with the WIMP paradigm. In this paper, we focus on the subset of WIMP models in which the dark matter annihilates in the early universe through couplings to either the Standard Model $Z$ or the Standard Model Higgs boson. Considering fermionic, scalar, and vector dark matter candidates within a model-independent context, we find that the overwhelming majority of these dark matter candidates are already ruled out by existing experiments. In the case of $Z$ mediated dark matter, the only scenarios that are not currently excluded are those in which the dark matter is a fermion with an axial coupling and with a mass either within a few GeV of the $Z$ resonance ($m_{\rm DM} \simeq m_Z/2$) or greater than 200 GeV, or with a vector coupling and with $m_{\rm DM} > 6$ TeV. Several Higgs mediated scenarios are currently viable if the mass of the dark matter is near the Higgs pole ($m_{\rm DM} \simeq m_H/2$). Otherwise, the only scenarios that are not excluded are those in which the dark matter is a scalar (vector) heavier than 400 GeV (1160 GeV) with a Higgs portal coupling, or a fermion with a pseudoscalar (CP violating) coupling to the Standard Model Higgs boson. With the exception of dark matter with a purely pseudoscalar coupling to the Higgs, it is anticipated that planned direct detection experiments will probe nearly the entire range of models considered in this study.

We use high spectral resolution (R > 8000) data covering 3800-13000\r{A} to study the physical conditions of the broad line region (BLR) of nine nearby Seyfert 1 galaxies. Up to six broad HI lines are present in each spectrum. A comparison - for the first time using simultaneous optical to near-infrared observations - to photoionisation calculations with our devised simple scheme yields the extinction to the BLR at the same time as determining the density and photon flux, and hence distance from the nucleus, of the emitting gas. This points to a typical density for the HI emitting gas of 10$^{11}$cm$^{-3}$ and shows that a significant amount of this gas lies at regions near the dust sublimation radius, consistent with theoretical predictions. We also confirm that in many objects the line ratios are far from case B, the best-fit intrinsic broad-line H$\alpha$/H$\beta$ ratios being in the range 2.5-6.6 as derived with our photoionization modeling scheme. The extinction to the BLR, based on independent estimates from HI and HeII lines, is A$_V$ $\le$ 3 for Seyfert 1-1.5s, while Seyfert 1.8-1.9s have A$_V$ in the range 4-8. A comparison of the extinction towards the BLR and narrow line region (NLR) indicates that the structure obscuring the BLR exists on scales smaller than the NLR. This could be the dusty torus, but dusty nuclear spirals or filaments could also be responsible. The ratios between the X-ray absorbing column N$_H$ and the extinction to the BLR are consistent with the Galactic gas-to-dust ratio if N$_H$ variations are considered.

Low mass, self-gravitating accretion disks admit quasi-steady,`gravito-turbulent' states in which cooling balances turbulent viscous heating. However, numerical simulations show that gravito-turbulence cannot be sustained beyond dynamical timescales when the cooling rate or corresponding turbulent viscosity is too large. The result is disk fragmentation. We motivate and quantify an interpretation of disk fragmentation as the inability to maintain gravito-turbulence due to formal secondary instabilities driven by: 1) cooling, which reduces pressure support; and/or 2) viscosity, which reduces rotational support. We analyze the axisymmetric gravitational stability of viscous, non-adiabatic accretion disks with internal heating, external irradiation, and cooling in the shearing box approximation. We consider parameterized cooling functions in 2D and 3D disks, as well as radiative diffusion in 3D. We show that generally there is no critical cooling rate/viscosity below which the disk is formally stable, although interesting limits appear for unstable modes with lengthscales on the order of the disk thickness. We apply this new linear theory to protoplanetary disks subject to gravito-turbulence modeled as an effective viscosity, and cooling regulated by dust opacity. We find that viscosity renders the disk beyond $\sim 60$AU dynamically unstable on radial lengthscales a few times the local disk thickness. This is coincident with the empirical condition for disk fragmentation based on a maximum sustainable stress. We suggest turbulent stresses can play an active role in realistic disk fragmentation by removing rotational stabilization against self-gravity, and that the observed transition in behavior from gravito-turbulent to fragmenting may reflect instability of the gravito-turbulent state itself.

We quantify the thermodynamic requirement for the Vertical Shear Instability and evaluate its relevance to realistic protoplanetary disks as a potential route to hydrodynamic turbulence.

In order to understand the interactions between active galactic nuclei (AGN) and star formation during the evolution of galaxies, we investigate 142 galaxies detected in both X-ray and 70{\mu}m observations in the COSMOS (Cosmic Evolution Survey) field. All of our data are obtained from the archive, X-ray point source catalogs from Chandra and XMM-Newton observations; far-infrared 70{\mu}m point source catalog from Spitzer-MIPS observations. Although the IRAC [3.6{\mu}m]-[4.5{\mu}m] vs. [5.8{\mu}m]-[8.0{\mu}m] colours of our sample indicate that only ~63% of our sources would be classified as AGN, the ratio of the rest-frame 2-10 keV luminosity to the total infrared luminosity (8-1000{\mu}m) shows that all of the sample has comparatively higher X-ray luminosity than that expected from pure star-forming galaxies, suggesting the presence of an AGN in all of our sources. From the analysis of the X-ray hardness ratio, we find that sources with both 70{\mu}m and X-ray detection tend to have a higher hardness ratio relative to the whole X-ray selected source population, suggesting the presence of more X-ray absorption in the 70{\mu}m detected sources. In addition, we find that the observed far-infrared colours of 70{\mu}m detected sources with and without X-ray emission are similar, suggesting the far-infrared emission could be mainly powered by star formation.

We discuss the dense molecular gas in central regions of nearby Seyfert galaxies, and report new arcsec resolution observations of HCN(1-0) and HCO$^+$(1-0) for 3 objects. In NGC 3079 the lines show complex profiles as a result of self-absorption and saturated continuum absorption. H$^{13}$CN reveals the continuum absorption profile, with a peak close to the galaxy's systemic velocity that traces disk rotation, and a second feature with a blue wing extending to $-350$km s$^{-1}$ that most likely traces a nuclear outflow. The morphological and spectral properties of the emission lines allow us to constrain the dense gas dynamics. We combine our kinematic analysis for these 3 objects, as well as another with archival data, with a previous comparable analysis of 4 other objects, to create a sample of 8 Seyferts. In 7 of these, the emission line kinematics imply thick disk structures on radial scales of $\sim$100pc, suggesting such structures are a common occurrence. We find a relation between the circumnuclear LHCN and Mdyn that can be explained by a gas fraction of 10% and a conversion factor {\alpha}HCN $\sim$ 10 between gas mass and HCN luminosity. Finally, adopting a different perspective to probe the physical properties of the gas around AGN, we report on an analysis of molecular line ratios which indicates that the clouds in this region are not self-gravitating.

The optical classification of a Seyfert galaxy and whether it is considered X-ray absorbed are often used interchangeably. But there are many borderline cases and also numerous examples where the optical and X-ray classifications appear to be in conflict. In this article we re-visit the relation between optical obscuration and X-ray absorption in AGNs. We make use of our "dust color" method (Burtscher et al. 2015) to derive the optical obscuration A_V and consistently estimated X-ray absorbing columns using 0.3--150 keV spectral energy distributions. We also take into account the variable nature of the neutral gas column N_H and derive the Seyfert sub-classes of all our objects in a consistent way. We show in a sample of 25 local, hard-X-ray detected Seyfert galaxies (log L_X / (erg/s) ~ 41.5 - 43.5) that there can actually be a good agreement between optical and X-ray classification. If Seyfert types 1.8 and 1.9 are considered unobscured, the threshold between X-ray unabsorbed and absorbed should be chosen at a column N_H = 10^22.3 / cm^2 to be consistent with the optical classification. We find that N_H is related to A_V and that the N_H/A_V ratio is approximately Galactic or higher in all sources, as indicated previously. But in several objects we also see that deviations from the Galactic ratio are only due to a variable X-ray column, showing that (1) deviations from the Galactic N_H/A_V can simply be explained by dust-free neutral gas within the broad line region in some sources, that (2) the dust properties in AGNs can be similar to Galactic dust and that (3) the dust color method is a robust way to estimate the optical extinction towards the sublimation radius in all but the most obscured AGNs.

We suggest stellar oscillations are responsible for the strange radio behaviors of Anomalous X-ray pulsars and soft Gamma-ray repeaters (AXP/SGRs), within the framework of both solid quark star model and magnetar model. In solid quark star model, the extra voltage provided by oscillations activates the star from under death line to above death line. In magnetar model, oscillations enlarge the radio beam so that increase the possibility to detect it. Later radio emission decays and vanishes as oscillations damp.

The vertical shear instability (VSI) offers a potential hydrodynamic mechanism for angular momentum transport in protoplanetary disks (PPDs). The VSI is driven by a weak vertical gradient in the disk's orbital motion, but must overcome vertical buoyancy, a strongly stabilizing influence in cold disks, where heating is dominated by external irradiation. Rapid radiative cooling reduces the effective buoyancy and allows the VSI to operate. We quantify the cooling timescale $t_c$ needed for efficient VSI growth, through a linear analysis of the VSI with cooling in vertically global, radially local disk models. We find the VSI is most vigorous for rapid cooling with $t_c<\Omega_\mathrm{K}^{-1}h|q|/(\gamma -1)$ in terms of the Keplerian orbital frequency, $\Omega_\mathrm{K}$; the disk's aspect-ratio, $h\ll1$; the radial power-law temperature gradient, $q$; and the adiabatic index, $\gamma$. For longer $t_c$, the VSI is much less effective because growth slows and shifts to smaller length scales, which are more prone to viscous or turbulent decay. We apply our results to PPD models where $t_c$ is determined by the opacity of dust grains. We find that the VSI is most effective at intermediate radii, from $\sim5$AU to $\sim50$AU with a characteristic growth time of $\sim30$ local orbital periods. Growth is suppressed by long cooling times both in the opaque inner disk and the optically thin outer disk. Reducing the dust opacity by a factor of 10 increases cooling times enough to quench the VSI at all disk radii. Thus the formation of solid protoplanets, a sink for dust grains, can impede the VSI.

We investigate properties of Galactic microlensing events in which a stellar object is lensed by a neutron star. For an all-sky photometric microlensing survey, we determine the number of lensing events caused by $\sim10^{5}$ potentially-observable radio pulsars to be $\sim0.2\ \rm{yr^{-1}}$ for $10^{10}$ background stellar sources. We expect a few detectable events per year for the same number of background sources from an astrometric microlensing survey. We show that such a study could lead to precise measurements of radio pulsar masses. For instance, if a pulsar distance could be constrained through radio observations, then its mass would be determined with a precision of $\sim10\%$. We also investigate the time-scale distributions for neutron star events, finding that they are much shorter than had been previously thought. For photometric events towards the Galactic centre that last $\sim15$ days, around $7\%$ will have a neutron star lens. This fraction drops rapidly for longer time-scales. Away from the bulge region we find that neutron stars will contribute $\sim40\%$ of the events that last less than $\sim10$ days. These results are in contrast to earlier work which found that the maximum fraction of neutron star events would occur on time-scales of hundreds of days.

We describe a complete volume limited sample of nearby active galaxies selected by their 14-195keV luminosity, and outline its rationale for studying the mechanisms regulating gas inflow and outflow. We describe also a complementary sample of inactive galaxies, selected to match the AGN host galaxy properties. The active sample appears to have no bias in terms of AGN type, the only difference being the neutral absorbing column which is two orders of magnitude greater for the Seyfert 2s. In the luminosity range spanned by the sample, log L_{14-195keV} [erg/s] = 42.4-43.7, the optically obscured and X-ray absorbed fractions are 50-65%. The similarity of these fractions to more distant spectroscopic AGN samples, although over a limited luminosity range, suggests that the torus does not strongly evolve with redshift. Our sample confirms that X-ray unabsorbed Seyfert 2s are rare, comprising not more than a few percent of the Seyfert 2 population. At higher luminosities, the optically obscured fraction decreases (as expected for the increasing dust sublimation radius), but the X-ray absorbed fraction changes little. We argue that the cold X-ray absorption in these Seyfert 1s can be accounted for by neutral gas in clouds that also contribute to the broad line region (BLR) emission; and suggest that a geometrically thick neutral gas torus co-exists with the BLR and bridges the gap to the dusty torus.

We combine two approaches to isolate the AGN luminosity at near-infrared wavelengths and relate the near-IR pure AGN luminosity to other tracers of the AGN. Using integral-field spectroscopic data of an archival sample of 51 local AGNs, we estimate the fraction of non-stellar light by comparing the nuclear equivalent width of the stellar 2.3 micron CO absorption feature with the intrinsic value for each galaxy. We compare this fraction to that derived from a spectral decomposition of the integrated light in the central arc second and find them to be consistent with each other. Using our estimates of the near-IR AGN light, we find a strong correlation with presumably isotropic AGN tracers. We show that a significant offset exists between type 1 and type 2 sources in the sense that type 1 sources are 7 (10) times brighter in the near-IR at log L_MIR = 42.5 (log L_X = 42.5). These offsets only becomes clear when treating infrared type 1 sources as type 1 AGNs. All AGNs have very red near-to-mid-IR dust colors. This, as well as the range of observed near-IR temperatures, can be explained with a simple model with only two free parameters: the obscuration to the hot dust and the ratio between the warm and hot dust areas. We find obscurations of A_V (hot) = 5 - 15 mag for infrared type 1 sources and A_V (hot) = 15 - 35 mag for type 2 sources. The ratio of hot dust to warm dust areas of about 1000 is nicely consistent with the ratio of radii of the respective regions as found by infrared interferometry.

Several observations of transition discs show lopsided dust-distributions. A potential explanation is the formation of a large-scale vortex acting as a dust-trap at the edge of a gap opened by a giant planet. Numerical models of gap-edge vortices have thus far employed locally isothermal discs, but the theory of this vortex-forming or `Rossby wave' instability was originally developed for adiabatic discs. We generalise the study of planetary gap stability to non-isothermal discs using customised numerical simulations of disc-planet systems where the planet opens an unstable gap. We include in the energy equation a simple cooling function with cooling timescale $t_c=\beta\Omega_k^{-1}$, where $\Omega_k$ is the Keplerian frequency, and examine the effect of $\beta$ on the stability of gap edges and vortex lifetimes. We find increasing $\beta$ lowers the growth rate of non-axisymmetric perturbations, and the dominant azimuthal wavenumber $m$ decreases. We find a quasi-steady state consisting of one large-scale, over-dense vortex circulating the outer gap edge, typically lasting $O(10^3)$ orbits. Vortex lifetimes were found to generally increase with cooling times up to an optimal value, beyond which vortex lifetimes decrease. This non-monotonic dependence is qualitatively consistent with recent studies using strictly isothermal discs that vary the disc aspect ratio.

We describe a new mechanism that leads to the destabilisation of non-axisymmetric waves in astrophysical discs with an imposed radial temperature gradient. This might apply, for example, to the outer parts of protoplanetary discs. We use linear density wave theory to show that non-axisymmetric perturbations generally do not conserve their angular momentum in the presence of a forced temperature gradient. This implies an exchange of angular momentum between linear perturbations and the background disc. In particular, when the disturbance is a low-frequency trailing wave and the disc temperature decreases outwards, this interaction is unstable and leads to the growth of the wave. We demonstrate this phenomenon through numerical hydrodynamic simulations of locally isothermal discs in 2D using the FARGO code and in 3D with the ZEUS-MP and PLUTO codes. We consider radially structured discs with a self-gravitating region which remains stable in the absence of a temperature gradient. However, when a temperature gradient is imposed we observe exponential growth of a one-armed spiral mode (azimuthal wavenumber $m=1$) with co-rotation radius outside the bulk of the spiral arm, resulting in a nearly-stationary one-armed spiral pattern. The development of this one-armed spiral does not require the movement of the central star, as found in previous studies. Because destabilisation by a forced temperature gradient does not explicitly require disc self-gravity, we suggest this mechanism may also affect low-frequency one-armed oscillations in non-self-gravitating discs.

We study the magnetospheric activity in the polar cap region of pulsars under stellar oscillations. The toroidal oscillation of the star propagates into the magnetosphere, which provides additional voltage due to unipolar induction, changes Goldreich-Julian charge density from the traditional value due to rotation, and hence, influences particle acceleration. We present a general solution of the effect of oscillations within the framework of the inner vacuum gap model, and consider three different inner gap modes controlled by curvature radiation, inverse Compton scattering, and two photon annihilation, respectively. With different pulsar parameters and oscillation amplitudes, one of three modes would play a dominant role in defining the gap properties. When the amplitude of oscillation exceeds a critical value, mode changing would occur. Oscillations also lead to change of the size of the polar cap. As applications, we show the inner gap properties under oscillations in both normal pulsars and anomalous X-ray pulsars / soft gamma-ray repeaters (AXPs/SGRs). We interpret the onset of radio emission after glitches/flares in AXPs/SGRs as due to oscilation-driven magnetic activities in these objects, within the framework of both the magnetar model and the solid quark star model. Within the magnetar model, radio activation may be caused by the enlargement of the effective polar cap angle and the radio emission beam due to oscillation; whereas within the solid quark star angle, it may be caused by activation of the pulsar inner gap from below the radio emission death line due to an oscillation-induced voltage enhancement. The model can also explain the glitch-induced radio profile change observed in PSR J1119-6127.

We analyse the 2-dimensional distribution and kinematics of the stars as well as molecular and ionised gas in the central few hundred parsecs of 5 active and 5 matched inactive galaxies. The equivalent widths of the Br-gamma line indicate there is no on-going star formation in their nuclei, although recent (terminated) starbursts are possible in the active galaxies. The stellar velocity fields show no signs of non-circular motions, while the 1-0S(1) H_2 kinematics exhibit significant deviations from simple circular rotation. In the active galaxies the H_2 kinematics reveal inflow and outflow superimposed on disk rotation. Steady-state circumnuclear inflow is seen in three AGN, and hydrodynamical models indicate it can be driven by a large scale bar. In three of the five AGN, molecular outflows are spatially resolved. The outflows are oriented such that they intersect, or have an edge close to, the disk - which may be the source of molecular gas in the outflow. The relatively low speeds imply the gas will fall back onto the disk; and with moderate outflow rates, they will have only a local impact on the host galaxy. H_2 was detected in two inactive galaxies. These exhibit chaotic circumnuclear dust morphologies and have molecular structures that are counter-rotating with respect to the main gas component, which could lead to gas inflow in the near future. In our sample, all four galaxies with chaotic dust morphology in the circumnuclear region exist in moderately dense groups with 10-15 members where accretion of stripped gas can easily occur.

Magneto-rotational instability (MRI) and gravitational instability (GI) are the two principle routes to turbulent angular momentum transport in accretion disks. Protoplanetary disks may develop both. This paper aims to reinvigorate interest in the study of magnetized massive protoplanetary disks, starting from the basic issue of stability. The local linear stability of a self-gravitating, uniformly magnetized, differentially rotating, three-dimensional stratified disk subject to axisymmetric perturbations is calculated numerically. The formulation includes resistivity. It is found that the reduction in the disk thickness by self-gravity can decrease MRI growth rates; the MRI becomes global in the vertical direction, and MRI modes with small radial length scales are stabilized. The maximum vertical field strength that permits the MRI in a strongly self-gravitating polytropic disk with polytropic index $\Gamma=1$ is estimated to be $B_{z,\mathrm{max}} \simeq c_{s0}\Omega\sqrt{\mu_0/16\pi G} $, where $c_{s0}$ is the midplane sound speed and $\Omega$ is the angular velocity. In massive disks with layered resistivity, the MRI is not well-localized to regions where the Elsasser number exceeds unity. For MRI modes with radial length scales on the order of the disk thickness, self-gravity can enhance density perturbations, an effect that becomes significant in the presence of a strong toroidal field, and which depends on the symmetry of the underlying MRI mode. In gravitationally unstable disks where GI and MRI growth rates are comparable, the character of unstable modes can transition smoothly between MRI and GI. Implications for non-linear simulations are discussed briefly.

We present a rare observation of direct magnetic interaction between an activating filament and a coronal hole (CH). The filament was a quiescent one located at the northwest of the CH. It underwent a nonradial activation, during which filament material constantly fell and intruded into the CH. As a result, the CH was clearly destroyed by the intrusion. Brightenings appeared at the boundaries and in the interior of the CH, meanwhile, its west boundaries began to retreat and the area gradually shrank. It is noted that the CH went on shrinking after the end of the intrusion and finally disappeared entirely. Following the filament activation, three coronal dimmings (D1-D3) were formed, among which D1 and D2 persisted throughout the complete disappearance of the CH. The derived coronal magnetic configuration shows that the filament was located below an extended loop system which obviously linked D1 to D2. By comparison with this result of extrapolation, our observations imply that the interaction between the filament and the CH involved direct intrusion of the filament material to the CH and the disappearance of the CH might be due to interchange reconnection between the expanding loop system and the CH's open field.

Gap formation by giant planets in self-gravitating disks may lead to a gravitational edge instability (GEI). We demonstrate this GEI with global 3D and 2D self-gravitating disk-planet simulations using the ZEUS, PLUTO and FARGO hydrodynamic codes. High resolution 2D simulations show that an unstable outer gap edge can lead to outwards orbital migration. Our results have important implications for theories of giant planet formation in massive disks.

Vortex formation through the Rossby wave instability (RWI) in protoplanetary discs has been invoked to play a role in planet formation theory, and suggested to explain the observation of large dust asymmetries in several transitional discs. However, whether or not vortex formation operates in layered accretion discs, i.e. models of protoplanetary discs including dead zones near the disc midplane --- regions that are magnetically inactive and the effective viscosity greatly reduced --- has not been verified. As a first step toward testing the robustness of vortex formation in layered discs, we present non-linear hydrodynamical simulations of global 3D protoplanetary discs with an imposed kinematic viscosity that increases away from the disc midplane. Two sets of numerical experiments are performed: (i) non-axisymmetric instability of artificial radial density bumps in viscous discs; (ii) vortex-formation at planetary gap edges in layered discs. Experiment (i) shows that the linear instability is largely unaffected by viscosity and remains dynamical. The disc-planet simulations also show the initial development of vortices at gap edges, but in layered discs the vortices are transient structures which disappear well into the non-linear regime. We suggest that the long term survival of columnar vortices, such as those formed via the RWI, requires low effective viscosity throughout the vertical extent of the disc, so such vortices do not persist in layered discs.

Composite dark matter candidates, which can arise from new strongly-coupled sectors, are well-motivated and phenomenologically interesting, particularly in the context of asymmetric generation of the relic density. In this work, we employ lattice calculations to study the electromagnetic form factors of electroweak-neutral dark-matter baryons for a three-color, QCD-like theory with Nf = 2 and 6 degenerate fermions in the fundamental representation. We calculate the (connected) charge radius and anomalous magnetic moment, both of which can play a significant role for direct detection of composite dark matter. We find minimal Nf dependence in these quantities. We generate mass-dependent cross-sections for dark matter-nucleon interactions and use them in conjunction with experimental results from XENON100, excluding dark matter candidates of this type with masses below 10 TeV.

The Atacama Large Millimeter Array (ALMA) has been returning images of transitional disks in which large asymmetries are seen in the distribution of mm-sized dust in the outer disk. The explanation in vogue borrows from the vortex literature by suggesting that these asymmetries are the result of dust trapping in giant vortices, excited via Rossby wave instability (RWI) at planetary gap edges. Due to the drag force, dust trapped in vortices will accumulate in the center, and diffusion is needed to maintain a steady state over the lifetime of the disk. While previous work derived semi-analytical models of the process, in this paper we provide analytical steady-state solutions. Exact solutions exist for certain vortex models. The solution is determined by the vortex rotation profile, the gas scale height, the vortex aspect ratio, and the ratio of dust diffusion to gas-dust friction. In principle, all these quantities can be derived from observations, which would give validation of the model, also giving constrains on the strength of the turbulence inside the vortex core. Based on our solution, we derive quantities such as the gas-dust contrast, the trapped dust mass, and the dust contrast at the same orbital location. We apply our model to the recently imaged Oph IRS 48 system, finding values within the range of the observational uncertainties.

It has been suggested that long-period giant planets, such as HD 95086b and HR 8799bcde, may have formed through gravitational instability of protoplanetary discs. However, self-gravitating disc-satellite interaction can lead to the formation of a gravitationally unstable gap. Such an instability significantly affects the orbital migration of gap-opening perturbers in massive discs. We use 2D hydrodynamical simulations to examine the role of planet mass on the gravitational stability of gaps and its impact on orbital migration. We consider giant planets with planet-to-star mass ratio q=0.0003 to q=0.003, in a self-gravitating disc with disc-to-star mass ratio M_d/M_*=0.08, aspect ratio h=0.05, and Keplerian Toomre parameter Q = 1.5 at 2.5 times the planet's initial orbital radius. Fixed-orbit simulations show that all planet masses we consider open gravitationally unstable gaps, but the instability is stronger and develops sooner with increasing planet mass. The disc-on-planet torques typically become more positive with increasing planet mass. In freely-migrating simulations, we observe faster outward migration with increasing planet mass, but only for planet masses capable of opening unstable gaps early on. For q=0.0003, the planet undergoes rapid inward type III migration before it can open a gap. For q=0.0013 we find it is possible to balance the tendency for inward migration by the positive torques due to an unstable gap, but only for a few 10's of orbital periods. We find the unstable outer gap edge can trigger outward type III migration, sending the planet to twice its initial orbital radius on dynamical timescales. This triggered type III migration may be a way of bringing giant planets to large orbital distances. On the other hand, our results also imply gap-opening is not a sufficient condition for the in situ formation of giant planets on wide orbits through disc fragmentation.

Numerical simulations are presented to study the stability of gaps opened by giant planets in 3D self-gravitating disks. In weakly self-gravitating disks, a few vortices develop at the gap edge and merge on orbital time-scales. The result is one large but weak vortex with Rossby number -0.01. In moderately self-gravitating disks, more vortices develop and their merging is resisted on dynamical time-scales. Self-gravity can sustain multi-vortex configurations, with Rossby number -0.2 to -0.1, over a time-scale of order 100 orbits. Self-gravity also enhances the vortex vertical density stratification, even in disks with initial Toomre parameter of order 10. However, vortex formation is suppressed in strongly self-gravitating disks and replaced by a global spiral instability associated with the gap edge which develops during gap formation.

Astrophysical disks with localized radial structure, such as protoplanetary disks containing dead zones or gaps due to disk-planet interaction, may be subject to the non-axisymmetric Rossby wave instability (RWI) that lead to vortex-formation. The linear instability has recently been demonstrated in three-dimensional (3D) barotropic disks. It is the purpose of this study to generalize the 3D linear problem to include an energy equation, thereby accounting for baroclinity in three-dimensions. Linear stability calculations are presented for radially structured, vertically stratified, geometrically-thin disks with non-uniform entropy distribution in both directions. Polytropic equilibria are considered but adiabatic perturbations assumed. The unperturbed disk has a localized radial density bump making it susceptible to the RWI. The linearized fluid equations are solved numerically as a partial differential equation eigenvalue problem. Emphasis on the ease of method implementation is given. It is found that when the polytropic index is fixed and adiabatic index increased, non-uniform entropy has negligible effect on the RWI growth rate, but pressure and density perturbation magnitudes near a pressure enhancement increases away from the midplane. The associated meridional flow is also qualitatively changed from homentropic calculations. Meridional vortical motion is identified in the nonhomentropic linear solution, as well as in a nonlinear global hydrodynamic simulation of the RWI in an initially isothermal disk evolved adiabatically. Numerical results suggest buoyancy forces play an important role in the internal flow of Rossby vortices.

The linear Rossby wave instability (RWI) in global, 3D polytropic discs is revisited with a much simpler numerical method than that previously employed by the author. The governing partial differential equation is solved with finite differences in the radial direction and spectral collocation in the vertical direction. RWI modes are calculated subject to different upper disc boundary conditions. These include free surface, solid boundaries and variable vertical domain size. Boundary conditions that oppose vertical motion increase the instability growth rate by a few per cent. The magnitude of vertical flow throughout the fluid column can be affected but the overall flow pattern is qualitatively unchanged. Numerical results support the notion that the RWI is intrinsically two dimensional. This implies that inconsistent upper disc boundary conditions, such as vanishing enthalpy perturbation, may inhibit the RWI in 3D.

Numerical simulations of global three-dimensional (3D), self-gravitating discs with a gap opened by an embedded planet are presented. The simulations are customised to examine planetary gap stability. Previous results, obtained by Lin & Papaloizou from two-dimensional (2D) disc models, are reproduced in 3D. These include (i) the development of vortices associated with local vortensity minima at gap edges and their merging on dynamical timescales in weakly self-gravitating discs, (ii) the increased number of vortices as the strength of self-gravity is increased and their resisted merging, and (iii) suppression of the vortex instability and development of global spiral arms associated with local vortensity maxima in massive discs. The vertical structure of these disturbances are examined. In terms of the relative density perturbation, the vortex disturbance has weak vertical dependence when self-gravity is neglected. Vortices become more vertically stratified with increasing self-gravity. This effect is seen even when the unperturbed region around the planet's orbital radius has a Toomre stability parameter ~10. The spiral modes display significant vertical structure at the gap edge, with the midplane density enhancement being several times larger than that near the upper disc boundary. However, for both instabilities the vertical Mach number is typically a few per cent,and on average vertical motions near the gap edge do not dominate horizontal motions.

Numerical calculations of the linear Rossby wave instability (RWI) in global three-dimensional (3D) disks are presented. The linearized fluid equations are solved for vertically stratified, radially structured disks with either a locally isothermal or polytropic equation of state, by decomposing the vertical dependence of the perturbed hydrodynamic quantities into Hermite and Gegenbauer polynomials, respectively. It is confirmed that the RWI operates in 3D. For perturbations with vertical dependence assumed above, there is little difference in growth rates between 3D and two-dimensional (2D) calculations. Comparison between 2D and 3D solutions of this type suggest the RWI is predominantly a 2D instability and that three-dimensional effects, such as vertical motion, to be interpreted as a perturbative consequence of the dominant 2D flow. The vertical flow around co-rotation, where vortex-formation is expected, is examined. In locally isothermal disks the expected vortex center remains in approximate vertical hydrostatic equilibrium. For polytropic disks the vortex center has positive vertical velocity, whose magnitude increases with decreasing polytropic index $n$.

Protoplanetary discs may become dynamically unstable due to structure induced by an embedded giant planet. In this thesis, I discuss the stability of such systems and explore the consequence of instability on planetary migration. I begin with non-self-gravitating, low viscosity discs and show that giant planets induce shocks inside its co-orbital region, leading to a profile unstable to vortex formation around a potential vorticity minimum. This instability is commonly known as the vortex or Rossby wave instability. Vortex-planet interaction lead to episodic phases of migration, which can be understood in the framework of type III migration. I then examine the effect of disc self-gravity on gap stability. The linear theory of the Rossby wave instability is extended to include disc gravity, which shows that self-gravity is effective at stabilising the vortex instability at small azimuthal wavenumber. This is consistent with the observation that more vortices develop with increasing disc mass in hydrodynamic simulations. Vortices in self-gravitating discs also resist merging, and is most simply understood as pair-vortices undergoing mutual horsehoe turns upon encounter. I show that in sufficiently massive discs vortex modes are suppressed. Instead, global spiral instabilities develop which are associated with a potential vorticity maximum at the gap edge. These edge modes can be physically understood as a result of unstable interaction between the gap edge and the exterior disc through gravity. I show the spiral arms can provide a positive torque on the planet, leading to fast migration outwards. I confirm the above results, obtained from razor-thin disc models, persist in three-dimensions.

We present numerical simulations of disc-planet interactions where the planet opens a gravitationally unstable gap in an otherwise gravitationally stable disc. In our disc models, where the outer gap edge can be unstable to global spiral modes, we find that as we increase the surface density scale the gap becomes more unstable and the planet migrates outwards more rapidly. We show that the positive torque is provided by material brought into the planet's coorbital region by the spiral arms. This material is expected to execute horseshoe turns upon approaching the planet and hence torque it. Our results suggest that standard type II migration, applicable to giant planets in non-self-gravitating viscous discs, is likely to be significantly modified in massive discs when gravitational instabilities associated with the gap occur.

Young protostars embedded in circumstellar discs accrete from an angular momentum-rich mass reservoir. Without some braking mechanism, all stars should be spinning at or near break-up velocity. In this paper, we perform simulations of the self-gravitational collapse of an isothermal cloud using the ORION adaptive mesh refinement code and investigate the role that gravitational torques might play in the spin-down of the dense central object. While magnetic effects likely dominate for low mass stars, high mass and Population III stars might be less well magnetised. We find that gravitational torques alone prevent the central object from spinning up to more than half of its breakup velocity, because higher rotation rates lead to bar-like deformations that enable efficient angular momentum transfer to the surrounding medium. We also find that the long-term spin evolution of the central object is dictated by the properties of the surrounding disc. In particular, spiral modes with azimuthal wavenumber $m=2$ couple more effectively to its spin than the lopsided $m=1$ mode, which was found to inhibit spin evolution. We suggest that even in the absence of magnetic fields, gravitational torques may provide an upper limit on stellar spin, and that moderately massive circumstellar discs can cause long-term spin down.

We study the effect of disc self-gravity on vortex-forming instabilities associated with gaps opened by a Saturn mass planet in a protoplanetary disc. It is shown analytically and confirmed through linear calculations that vortex modes with low azimuthal mode number $,m,$ are stabilised by increasing self-gravity if the basic state is fixed. However, linear calculations show that the combined effect of self-gravity through the background and through the linear response shifts the most unstable vortex mode to higher $m.$ Nonlinear hydrodynamic simulations of planetary gaps show more vortices develop with increasing strength of self-gravity. For sufficiently large disc mass the vortex instability is suppressed and replaced by a new global instability, consistent with analytical expectations. In the nonlinear regime, vortex merging is increasingly delayed as the disc mass increases and multiple vortices may persist until the end of simulations. With self-gravity, the post-merger vortex is localised in azimuth and has similar structure to a Kida-like vortex. This is unlike the case without self-gravity where vortices merge to form a single vortex extended in azimuth. We also performed a series of supplementary simulations of co-orbital Kida-like vortices and found that self-gravity enables such vortices to execute horseshoe turns upon encountering each other. As a result vortex merging is avoided on time-scales where it would occur without self-gravity. Thus we suggest that mutual repulsion of self-gravitating vortices in a rotating flow is responsible for the delayed vortex merging above. The effect of self-gravity on vortex-induced migration is briefly discussed. We found that when self-gravity is included, the vortex-induced type III migration of Lin & Papaloizou (2010) is delayed but the extent of migration is unchanged.

We study the stability of gaps opened by a giant planet in a self-gravitating protoplanetary disc. We find a linear instability associated with both the self-gravity of the disc and local vortensity maxima which coincide with gap edges. For our models, these edge modes develop and extend to twice the orbital radius of a Saturn mass planet in discs with disc-to-star mass ratio >0.06, corresponding to a Toomre Q < 1.5 at the outer disc boundary. Unlike the local vortex-forming instabilities associated with gap edges in weakly or non-self-gravitating low viscosity discs, the edge modes are global and exist only in sufficiently massive discs, but for the typical viscosity values adopted for protoplanetary discs. Analytic modelling and linear calculations show edge modes may be interpreted as a localised disturbance associated with a gap edge inducing activity in the extended disc, through the launching of density waves excited at Lindblad resonances. Nonlinear hydrodynamic simulations are performed to investigate the evolution of edge modes in disc-planet systems. The form and growth rates of unstable modes are consistent with linear theory. Their dependence on viscosity and gravitational softening is also explored. We also performed a first study of the effect of edge modes on planetary migration. We found that if edge modes develop, then the average disc-on-planet torque becomes more positive with increasing disc mass. In simulations where the planet was allowed to migrate, although a fast type III migration could be seen that was similar to that seen in non-self-gravitating discs, we found that it was possible for the planet to interact gravitationally with the spiral arms associated with an edge mode and that this could result in the planet being scattered outwards. Thus orbital migration is likely to be complex and non monotonic in massive discs of the type we consider.

We study the type III migration of a Saturn mass planet in low viscosity discs. The planet is found to experience cyclic episodes of rapid decay in orbital radius, each amounting to a few Hill radii. We find this to be due to the scattering of large- scale vortices present in the disc. The origin and role of vortices in the context of type III migration is explored. It is shown through numerical simulations and semi- analytical modelling that spiral shocks induced by a sufficiently massive planet will extend close to the planet orbital radius. The production of vortensity across shock tips results in thin high vortensity rings with a characteristic width of the local scale height. For planets with masses equal to and above that of Saturn, the rings are co-orbital features extending the entire azimuth. Linear stability analysis show there exists unstable modes that are localised about local vortensity minima which coincide with gap edges. Simulations show that vortices are non-linear a outcome. We used hydrodynamic simulations to examine vortex-planet interactions. Their effect is present in discs with kinematic viscosity less than about an order of magnitude smaller than the typically adopted value of \nu = 10^{-5}\Omega_pr_p(0)^2, where r_p(0) and \Omega_p are the initial orbital radius and angular velocity of the planet respectively. We find that the magnitude of viscosity affects the nature of type III migration but not the extent of the orbital decay. The role of vortices as a function of initial disc mass is also explored and it is found that the amount of orbital decay during one episode of vortex-planet interaction is independent of initial disc mass. We incorporate the concept of the co-orbital mass deficit in the analysis of our results and link it to the presence of vortices at gap edges.

It is believed that some stars have two or more convection zones in close proximity near to the stellar photosphere. These zones are separated by convectively stable regions that are relatively narrow. Due to the close proximity of these regions it is important to construct mathematical models to understand the transport and mixing of passive and dynamic quantities. One key quantity of interest is a magnetic field, a dynamic vector quantity, that can drastically alter the convectively driven flows, and have an important role in coupling the different layers. In this paper we present the first investigation into the effect of an imposed magnetic field in such a geometry. We focus our attention on the effect of field strength and show that, while there are some similarities with results for magnetic field evolution in a single layer, new and interesting phenomena are also present in a three layer system.