We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

The detection of several radio-loud narrow-line Seyfert 1 (NLS1) galaxies by the Fermi Gamma-Ray Space Telescope hints at the existence of a rare, new class of gamma-ray emitting active galactic nuclei with low black hole masses. Like flat spectrum radio quasars (FSRQs), their gamma-ray emission is thought to be produced via the external Compton mechanism whereby relativistic jet electrons upscatter a photon field external to the jet, e.g. from the accretion disc, broad line region (BLR) and dusty torus, to higher energies. Here we study the origin of the gamma-ray emission in the lowest-redshift candidate among the currently-known gamma-ray emitting NLS1s, 1H 0323+342, and take a new approach. We observationally constrain the external photon field using quasi-simultaneous near-IR, optical and X-ray spectroscopy. Applying a one-zone leptonic jet model, we simulate the range of jet parameters for which this photon field, when Compton scattered to higher energies, can explain the gamma-ray emission. We find that the site of the gamma-ray emission lies well within the BLR and that the seed photons mainly originate from the accretion disc. The jet power that we determine, $1.0 \times 10^{45}$ erg s$^{-1}$, is approximately half the accretion disc luminosity. We show that this object is not simply a low-mass FSRQ, its jet is intrinsically less powerful than predicted by scaling a typical FSRQ jet by black hole mass and accretion rate. That gamma-ray emitting NLS1s appear to host underpowered jets may go some way to explaining why so few have been detected to date.

The new multi-wavelength monitoring campaign on NGC 5548 shows clearly that the variability of the UV/optical lightcurves lags by progressively longer times at longer wavelengths, as expected from reprocessing of an optically thick disk, but that the timescales are longer than expected for a standard Shakura-Sunyaev accretion disc. We build a full spectral-timing reprocessing model to simulate the UV/optical lightcurves of NGC 5548. We show that disc reprocessing of the observed hard X-ray lightcurve produces optical lightcurves with too much fast variability as well as too short a lag time. Supressing the fast variability requires an intervening structure preventing the hard X-rays from illuminating the disc. We propose this is the disc itself, perhaps due to atomic processes in the UV lifting the photosphere, increasing the scale-height, making it less dense and less able to thermalise, so that it radiates low temperature Comptonised emission as required to produce the soft X-ray excess. The outer edge of the puffed-up Comptonised disc region emits FUV flux, and can illuminate the outer thin blackbody disc but while this gives reprocessed variable emission which is much closer to the observed UV and optical lightcurves, the light travel lags are still too short to match the data. We reverse engineer a solution to match the observations and find that the luminosity and temperature of the lagged emission is not consistent with material at the light travel lag distance responding to the irradiating flux (either FUV or X-ray) from the AGN. We conclude that the UV/optical lags of NGC 5548 are not the light travel time from X-ray reprocessing, nor the light travel time from FUV reprocessing, but instead could be the timescale for the outer blackbody disc vertical structure to respond to the changing FUV illumination.

Flat Spectrum Radio Quasars (FSRQs) are the most powerful relativistic jets seen from supermassive black holes (BHs) accreting via a radiatively efficient thin disc. Their high energy emission is well modelled by highly relativistic electrons in the jet Compton upscattering an external source of seed photons, primarily from the broad line region. Strong Doppler boosting by the jet bulk motion makes these FSRQs readily detectable by the Fermi Large Area Telescope. We combine jet spectral models with scaling relations for the jet physical parameters as a function of mass and accretion rate. This does not match well to the Gamma-ray loud Narrow Line Seyfert 1s, assuming their low BH masses are reliable, but is able to predict much of the spectral evolution observed along the Blazar sequence. We use these models in conjunction with cosmological simulations of efficiently accreting BH number densities, and find that they overpredict the observed number of FSRQs by 2 orders of magnitude if all of these objects produce a FSRQ jet. We can better reproduce the observed numbers if jets are only produced by high spin BHs and BH spin is built from chaotically aligned accretion episodes so that high spin BHs are rare. However, this does not reproduce the observed redshift and mass accretion rate distributions of the FSRQs. This may indicate a redshift dependence in accretion mode, with sustained alignment accretion episodes being more prevalent at higher redshift, or that there is some other trigger for FSRQ jets.

We present a detailed multi-wavelength study of an unobscured, highly super-Eddington Type-1 QSO RX J0439.6-5311. We combine the latest XMM-Newton observation with all archival data from infrared to hard X-rays. The optical spectrum is very similar to that of 1H 0707-495 in having extremely weak [O {\sc iii}] and strong Fe {\sc ii} emission lines, although the black hole mass is probably slightly higher at $5-10 \times10^{6}~\rm M_{\odot}$. The broadband SED is uniquely well-defined due to the extremely low Galactic and intrinsic absorption, so the bolometric luminosity is tightly constrained. The optical/UV accretion disc continuum is seen down to 900 \AA, showing that there is a standard thin disc structure down to $R \ge$ 190-380 $R_{\rm g}$ and determining the mass accretion rate through the outer disc. This predicts a much higher bolometric luminosity than observed, indicating that there must be strong wind and/or advective energy losses from the inner disc, as expected for a highly super-Eddington accretion flow. Significant outflows are detected in both the NLR and BLR emission lines, confirming the presence of a wind. We propose a global picture for the structure of a super-Eddington accretion flow where the inner disc puffs up, shielding much of the potential NLR material, and show how inclination angle with respect to this and the wind can explain very different X-ray properties of RX J0439.6-5311 and 1H 0707-495. Therefore, this source provides strong supporting evidence that `simple' and `complex' super-Eddington NLS1s can be unified within the same accretion flow scenario but with different inclination angles. We also propose that these extreme NLS1s could be the low-redshift analogs of weak emission-line quasars (WLQs).

Swift monitoring of NGC 4151 with ~6 hr sampling over a total of 69 days in early 2016 is used to construct light curves covering five bands in the X-rays (0.3-50 keV) and six in the ultraviolet (UV)/optical (1900-5500 A). The three hardest X-ray bands (>2.5 keV) are all strongly correlated with no measurable interband lag while the two softer bands show lower variability and weaker correlations. The UV/optical bands are significantly correlated with the X-rays, lagging ~3-4 days behind the hard X-rays. The variability within the UV/optical bands is also strongly correlated, with the UV appearing to lead the optical by ~0.5-1 day. This combination of >~3 day lags between the X-rays and UV and <~1 day lags within the UV/optical appears to rule out the "lamp-post" reprocessing model in which a hot, X-ray emitting corona directly illuminates the accretion disk, which then reprocesses the energy in the UV/optical. Instead, these results appear consistent with the Gardner & Done picture in which two separate reprocessings occur: first, emission from the corona illuminates an extreme-UV-emitting toroidal component that shields the disk from the corona; this then heats the extreme-UV component which illuminates the disk and drives its variability.

Complex narrow line Seyfert 1s (NLS1s), such as 1H0707-495, differ from simple NLS1s like PG1244+026 by showing stronger broad spectral features at Fe K and larger amplitude flux variability. These are correlated: the strongest Fe K features are seen during deep dips in the light curves of complex NLS1s. There are two competing explanations for these features, one where a compact X-ray source on the spin axis of a highly spinning black hole approaches the horizon and the consequent strong relativistic effects focus the intrinsic flux onto the inner edge of a thin disc, giving a dim, reflection dominated spectrum. The other is that the deep dips are caused by complex absorption by clumps close to the hard X-ray source. The reflection dominated model is able to reproduce the very short 30s soft lag from reverberation seen in the complex NLS1 1H0707-495. However, it does not explain the characteristic switch to hard lags on longer timescales. Instead, a full model of propagating fluctuations coupled to reverberation can explain the switch in the simple NLS1 PG1244+026 using a low spin black hole. However PG1244+026 has a longer reverberation lag of $\sim 200$s. Here we extend the successful propagation-reverberation model for the simple NLS1 PG1244+026 to include the effect of absorption from clumps in a turbulent region above the disk. The resulting occultations of the inner accretion flow can introduce additional hard lags when relativistic effects are taken into account. This dilutes the soft lag from reverberation and shifts it to higher frequencies, making a smooth transition between the 200s lags seen in simple NLS1s to the 30s lags in complex NLS1s. These two classes of NLS1 could then be determined by inclination angle with respect to a clumpy, probably turbulent, failed wind structure on the disc.

We study the origin of the soft X-ray excess seen in the 'simple' Narrow Line Seyfert 1 galaxy PG1244+026 using all available spectral-timing information. This object shows the standard switch between soft leading the hard band on long timescales, to the opposite behaviour on short timescales. This is interpreted as a combination of intrinsic fluctuations propagating down through the accretion flow giving the soft lead, together with reflection of the hard X-rays giving the soft lag. We build a full model of the spectral and time variability including both propagation and reflection, and compare our model with the observed power spectra, coherence, covariance, lag-frequency and lag-energy spectra. We compare models based on a separate soft excess component with those based on reflection dominated soft emission. Reflection dominated spectra have difficulty reproducing the soft lead at low frequency since reflection will always lag. They also suffer from high coherence and nearly identical hard and soft band power spectra in disagreement with the observations. This is a direct result of the power law and reflection components both contributing to the hard and soft energy bands, and the small radii over which the relativistically smeared reflection is produced transmitting too much high frequency power into the soft band. Conversely, we find the separate soft excess models (where the inner disc radius is $>6R_g$) have difficulty reproducing the soft lag at high frequency, as reflected flux does not contribute enough signal to overwhelm the soft lead. However, reflection should also be accompanied by reprocessing and this should add to the soft excess at low energies. This model can quantitatively reproduce the switch from soft lead to soft lag seen in the data and reproduces well the observed power spectra and other timing features which reflection dominated models cannot.

Using simulations of box/peanut- (B/P-) shaped bulges, we explore the nature of the X-shape of the Milky Way's bulge. An X-shape can be associated with a B/P-shaped bulge driven by a bar. By comparing in detail the simulations and the observations we show that the principal kinematic imprint of the X-shape is a minimum in the difference between the near and far side mean line-of-sight velocity along the minor axis. This minimum occurs at around |b| = 4{\deg}, which is close to the lower limit at which the X-shape can be detected. No coherent signature of an X-shape can be found in Galactocentric azimuthal velocities, vertical velocities, or any of the dispersions. After scaling our simulations, we find that a best fit to the BRAVA data leads to a bar angle of 15{\deg}. We also explore a purely geometric method for determining the distance to the Galactic Centre by tracing the arms of the X-shape. We find that we are able to determine this ill-known distance to an accuracy of about 5% per cent with sufficiently accurate distance measurements for the red clump stars in the arms.

The dramatic relativistic jets pointing directly at us in BL Lac objects can be well modelled by bulk motion beaming of synchrotron self-Compton emission powered by a low Eddington fraction accretion flow. Nearly 500 of these AGN are seen in the 2nd Fermi Large Area Telescope catalogue of AGN. We combine the jet models which describe individual spectra with the expected jet parameter scalings with mass and mass accretion rate to predict the expected number of Fermi detected sources given the number densities of AGN from cosmological simulations. We select only sources with Eddington scaled mass accretion rate $< 0.01$ (i.e. radiatively inefficient flows), and include cooling, orientation effects, and the effects of absorption from pair production on the extragalactic IR background. These models overpredict the number of Fermi detected BL Lacs by a factor of 1000! This clearly shows that one of the underlying assumptions is incorrect, almost certainly that jets do not scale simply with mass and accretion rate. The most plausible additional parameter which can affect the region producing the Fermi emission is black hole spin. We can reproduce the observed numbers of BL Lacs if such relativistic jets are only produced by the highest spin ($a_*>0.8$) black holes, in agreement with the longstanding spin-jet paradigm. This also requires that high spins are intrinsically rare, as predicted by the cosmological simulations for growing black hole mass via chaotic (randomly aligned) accretion episodes, where only the most massive black holes have high spin due to black hole-black hole mergers.

The X-ray spectra of black hole binaries in the low/hard state first harden as the flux decreases, then soften. This change in behaviour has been variously attributed to either the X-rays switching from being produced in the flow to being dominated by the jet, or to the flow switching seed photons from the disc to self generated seed photons from cyclo-synchrotron. Here we build a simple truncated disc, hot inner flow, plus standard conical synchrotron jet model to explore what this predicts for the X-ray emission mechanism as a function of mass accretion rate. We find that the change in X-ray spectral index can be quantitatively (not just qualitatively) explained by the seed photon switch in the hot flow i.e. this supports models where the X-rays are always produced by the hot flow. By contrast, standard conical jet models are as radiatively inefficient as the hot flow so there is no transition in X-ray production mechanism with $\dot{m}$. Including the effects of electron cooling allows the jet X-rays to drop more slowly with accretion rate and hence overtake the X-rays from the hot flow, however this produces a corresponding change in the radio-X-ray correlation, which is not observed. We argue that the unbroken radio-X-ray correlation down to quiescence rules out the jet transition model as an explanation for the trend in X-ray spectral index. Our favoured model is then a truncated disc with an inner, hot, radiatively inefficient flow which always dominates the hard X-rays, coupled to a conical synchrotron jet which produces the radio emission. However, even this has issues at low $\dot{m}$ as the low optical depth and high temperature of the flow means that the Compton spectrum is not well approximated by a power law. This shows the need for a more sophisticated model for the electron distribution in the hot flow.

It has recently been discovered that the Galactic bulge is X-shaped, with the two southern arms of the X both crossing the lines of sight at l=0 and |b|>4, hence producing a double red clump (RC) in the bulge CMD. Dynamical models predict the formation of X-shaped bulges, as extreme cases of boxy-peanut bulges. However, since X-shaped bulges were known to be present only in external galaxies, models have never been compared to 3D kinematical data for individual stars. We study the orbital motion of Galactic bulge stars, in the two arms of the X in the southern hemisphere. The goal is to provide observational constraints to bulge formation models that predict the formation of X-shapes through bar dynamical instabilities. Radial velocities have been obtained for a sample of 454 bulge giants, roughly equally distributed between the bright and faint RC, in a field at (l,b)=(0,-6). Proper motions were derived for all RC stars in the same field by combining images from two epochs obtained 11 years apart. The proper motions for the spectroscopic subsample are analyzed taking into account the radial velocities and metallicities measured from near-IR Calcium triplet lines, for a subsample of 352 stars. The radial velocity distribution of stars in the bright RC, tracing the closer overdensity of bulge stars, shows an excess of stars moving towards the Sun. Similarly, an excess of stars receding from the Sun is seen in the far overdensity, traced by faint RC stars. This is explained with the presence of stars on elongated orbits, most likely streaming along the arms of the X-shaped bulge. Proper motions for these stars are consistent with qualitative predictions of dynamical models of peanut-shaped bulges. Surprisingly, stars on elongated orbits have preferentially metal poor metallicities, while the metal rich ones, in both overdensities, are preferentially found in more axysimmetric orbits.

The long-term dynamics of Oort cloud comets are studied under the influence of both the radial and the vertical components of the Galactic tidal field. Sporadic dynamical perturbation processes are ignored, such as passing stars, since we aim to study the influence of just the axisymmetric Galactic tidal field on the cometary motion and how it changes in time. We use a model of the Galaxy with a disc, bulge and dark halo, and a local disc density, and disc scale length constrained to fit the best available observational constraints. By integrating a few million of cometary orbits over 1 Gyr, we calculate the time variable flux of Oort cloud comets that enter the inner Solar System, for the cases of a constant Galactic tidal field, and a realistically varying tidal field which is a function of the Sun's orbit. The applied method calculates the evolution of the comets by using first-order averaged mean elements. We find that the periodicity in the cometary flux is complicated and quasi-periodic. The amplitude of the variations in the flux are of order 30%. The radial motion of the Sun is the chief cause of this behaviour, and should be taken into account when the Galactic influence on the Oort cloud comets is studied.

It has been suggested that a resonance between a rotating bar and stars in the solar neighbourhood can produce the so called 'Hercules stream'. Recently, a second bar may have been identified in the Galactic centre, the so called 'long bar', which is longer and much flatter than the traditional Galactic bar, and has a similar mass. We looked at the dynamical effects of both bars, separately and together, on orbits of stars integrated backwards from local position and velocities, and a model of the Galactic potential which includes the bars directly. Both bars can produce Hercules like features, and allow us to measure the rotation rate of the bar(s). We measure a pattern speed, for both bars, of 1.87 +/- 0.02 times the local circular frequency. This is on par with previous measurements for the Galactic bar, although we do adopt a slightly different Solar motion. Finally, we identify a new kinematic feature in local velocity space, caused by the long bar, which is tempting to identify with the high velocity 'Arcturus' stream.

We examine the dynamical effects on disk stars of a "long bar" in the Milky Way by inserting a triaxial rotating bar into an axisymmetric disk+bulge+dark halo potential and integrating 3-D orbits of 104 tracer stars over a period of 2 Gyr. The long bar has been detected via "clump giants" in the IR by Lopez-Corredoira et al. (2007), and is estimated to have semi-major axes of (3.9 : 0.6 : 0.1) kpc and a mass of 6 10^9 Msun. We find such a structure has a slight impact on the inner disk-system, moving tracers near to the bar into the bar-region, as well as into the bulge. These effects are under continuing study.

If the origin of life and the evolution of observers on a planet is favoured by atypical properties of a planet's host star, we would expect our Sun to be atypical with respect to such properties. The Sun has been described by previous studies as both typical and atypical. In an effort to reduce this ambiguity and quantify how typical the Sun is, we identify eleven maximally-independent properties that have plausible correlations with habitability, and that have been observed by, or can be derived from, sufficiently large, currently available and representative stellar surveys. By comparing solar values for the eleven properties, to the resultant stellar distributions, we make the most comprehensive comparison of the Sun to other stars. The two most atypical properties of the Sun are its mass and orbit. The Sun is more massive than 95 -/+ 2% of nearby stars and its orbit around the Galaxy is less eccentric than 93 +/- 1% of FGK stars within 40 parsecs. Despite these apparently atypical properties, a chi^2 -analysis of the Sun's values for eleven properties, taken together, yields a solar chi^2 = 8.39 +/- 0.96. If a star is chosen at random, the probability that it will have a lower value (be more typical) than the Sun, with respect to the eleven properties analysed here, is only 29 +/- 11%. These values quantify, and are consistent with, the idea that the Sun is a typical star. If we have sampled all reasonable properties associated with habitability, our result suggests that there are no special requirements for a star to host a planet with life.