Nuclei of Nearby Active Galaxies: NIR and Sub-mm Views and Finalizing the LINC-NIRVANA Fringe and Flexure Tracking System: Flexure and Temperature Behavior

The work presented in my thesis addresses the two cornerstones of modern astronomy: Observation and Instrumentation. Part I deals with the observation of two nearby active galaxies, the Seyfert 2 galaxy NGC 1433 and the Seyfert 1 galaxy NGC 1566, both at a distance of $\sim10$ Mpc, which are part of the Nuclei of Galaxies (NUGA) sample. It is well established that every galaxy harbors a super massive black hole (SMBH) at its center. Furthermore, there seems to be a fundamental correlation between the stellar bulge and SMBH masses. Simulations show that massive feedback, e.g., powerful outflows, in Quasi Stellar Objects (QSOs) has an impact on the mutual growth of bulge and SMBH. Nearby galaxies follow this relation but accrete mass at much lower rates. This gives rise to the following questions: Which mechanisms allow feeding of nearby Active Galactic Nuclei (AGN)? Is this feeding triggered by events, e.g., star formation, nuclear spirals, outflows, on $\sim500$ pc scales around the AGN? Does feedback on these scales play a role in quenching the feeding process? Does it have an effect on the star formation close to the nucleus? To answer these questions I have carried out observations with the Spectrograph for INtegral Field Observation in the Near Infrared (SINFONI) at the Very Large Telescope (VLT) situated on Cerro Paranal in Chile. I have reduced and analyzed the recorded data, which contain spatial and spectral information in the H-band ($1.45 \mic-1.85 \mic$) and K-band ($1.95 \mic-2.45 \mic$) on the central $10\arcsec\times10\arcsec$ of the observed galaxies. Additionally, Atacama Large Millimeter/Sub-millimeter Array (ALMA) data at $350$ GHz ($\sim0.87$ mm) as well as optical high resolution Hubble Space Telescope (HST) images are used for the analysis. For NGC 1433 I deduce from comparison of the distributions of gas, dust, and intensity of highly ionized emission lines that the galaxy center lies $\sim70$ pc north-northwest of the prior estimate. A velocity gradient is observed at the new center, which I interpret as a bipolar outflow, a circum nuclear disk, or a combination of both. At least one dust and gas arm leads from a $r\sim200$ pc ring towards the nucleus and might feed the SMBH. Two bright warm H$_2$ gas spots are detected that indicate hidden star formation or a spiral arm-arm interaction. From the stellar velocity dispersion (SVD) I estimate a SMBH mass of $\sim1.74\times10^7$ \msol. For NGC 1566 I observe a nuclear gas disk of $\sim150$ pc in radius with a spiral structure. I estimate the total mass of this disk to be $\sim5.4\times10^7$ \msol. What mechanisms excite the gas in the disk is not clear. Neither can the existence of outflows be proven nor is star formation detected over the whole disk. On one side of the spiral structure I detect a star forming region with an estimated star formation rate of $\sim2.6\times10^{-3}$ \msol\ yr$^{-1}$. From broad Br$\gamma$ emission and SVD I estimate a mean SMBH mass of $\sim5.3\times10^6$ \msol\ with an Eddington ratio of $\sim2\times10^{-3}$. Part II deals with the final tests of the Fringe and Flexure Tracker (FFTS) for LBT INterferometric Camera and the NIR/Visible Adaptive iNterferometer for Astronomy (LINC-NIRVANA) at the Large Binocular Telescope (LBT) in Arizona, USA, which I conducted. The FFTS is the subsystem that combines the two separate beams of the LBT and enables near-infrared interferometry with a significantly large field of view. The FFTS has a cryogenic system and an ambient temperature system which are separated by the baffle system. I redesigned this baffle to guarantee the functionality of the system after the final tests in the Cologne cryostat. The redesign did not affect any scientific performance of LINC-NIRVANA. I show in the final cooldown tests that the baffle fulfills the temperature requirement and stays $<110$ K whereas the moving stages in the ambient system stay $>273$ K, which was not given for the old baffle design. Additionally, I test the tilting flexure of the whole FFTS and show that accurate positioning of the detector and the tracking during observation can be guaranteed.

Very large array detection of radio recombination lines from the radio nucleus of NGC 253: Ionization by a weak active galactic nucleus, an obscured super star cluster, or a compact supernova remnant?

We have imaged the H92alpha and H75alpha radio recombination line (RRL) emissions from the starburst galaxy NGC 253 with a resolution of similar to4 pc. The peak of the RRL emission at both frequencies coincides with the unresolved radio nucleus. Both lines observed toward the nucleus are extremely wide, with FWHMs of similar to200 km s(-1). Modeling the RRL and radio continuum data for the radio nucleus shows that the lines arise in gas whose density is similar to10(4) cm(-3) and mass is a few thousand M., which requires an ionizing flux of (6-20) x 10(51) photons s(-1). We consider a supernova remnant (SNR) expanding in a dense medium, a star cluster, and also an active galactic nucleus (AGN) as potential ionizing sources. Based on dynamical arguments, we rule out an SNR as a viable ionizing source. A star cluster model is considered, and the dynamics of the ionized gas in a stellar-wind driven structure are investigated. Such a model is only consistent with the properties of the ionized gas for a cluster younger than similar to10(5) yr. The existence of such a young cluster at the nucleus seems improbable. The third model assumes the ionizing source to be an AGN at the nucleus. In this model, it is shown that the observed X-ray flux is too weak to account for the required ionizing photon flux. However, the ionization requirement can be explained if the accretion disk is assumed to have a big blue bump in its spectrum. Hence, we favor an AGN at the nucleus as the source responsible for ionizing the observed RRLs. A hybrid model consisting of an inner advection-dominated accretion flow disk and an outer thin disk is suggested, which could explain the radio, UV, and X-ray luminosities of the nucleus.

Isotropic heating of galaxy cluster cores via rapidly reorienting active galactic nucleus jets

Active galactic nucleus (AGN) jets carry more than sufficient energy to stave off catastrophic cooling of the intracluster medium (ICM) in the cores of cool-core clusters. However, in order to prevent catastrophic cooling, the ICM must be heated in a near-isotropic fashion and narrow bipolar jets with P-jet = 10(44-45) erg s(-1), typical of radio AGNs at cluster centers, are inefficient in heating the gas in the transverse direction to the jets. We argue that due to existent conditions in cluster cores, the supermassive black holes (SMBHs) will, in addition to accreting gas via radiatively inefficient flows, experience short stochastic episodes of enhanced accretion via thin disks. In general, the orientation of these accretion disks will be misaligned with the spin axis of the black holes (BHs) and the ensuing torques will cause the BH's spin axis (and therefore the jet axis) to slew and rapidly change direction. This model not only explains recent observations showing successive generations of jet-lobes-bubbles in individual cool-core clusters that are offset from each other in the angular direction with respect to the cluster center, but also shows that AGN jets can heat the cluster core nearly isotropically on the gas cooling timescale. Our model does require that the SMBHs at the centers of cool-core clusters be spinning relatively slowly. Torques from individual misaligned disks are ineffective at tilting rapidly spinning BHs by more than a few degrees. Additionally, since SMBHs that host thin accretion disks will manifest as quasars, we predict that roughly 1-2 rich clusters within z < 0.5 should have quasars at their centers.

Spitzer power-law active galactic nucleus candidates in the chandra deep field-north

We define a sample of 62 galaxies in the Chandra Deep Field-North whose Spitzer IRAC SEDs exhibit the characteristic power-law emission expected of luminous AGNs. We study the multiwavelength properties of this sample and compare the AGNs selected in this way to those selected via other Spitzer color-color criteria. Only 55% of the power-law galaxies are detected in the X-ray catalog at exposures of >0.5 Ms, although a search for faint emission results in the detection of 85% of the power-law galaxies at the ≥2.5 σ detection level. Most of the remaining galaxies are likely to host AGNs that are heavily obscured in the X-ray. Because the power-law selection requires the AGNs to be energetically dominant in the near- and mid-infrared, the power-law galaxies comprise a significant fraction of the Spitzer-detected AGN population at high luminosities and redshifts. The high 24 μm detection fraction also points to a luminous population. The power-law galaxies comprise a subset of color-selected AGN candidates. A comparison with various mid-infrared color selection criteria demonstrates that while the color-selected samples contain a larger fraction of the X-ray-luminous AGNs, there is evidence that these selection techniques also suffer from a higher degree of contamination by star-forming galaxies in the deepest exposures. Considering only those power-law galaxies detected in the X-ray catalog, we derive an obscured fraction of 68% (2 : 1). Including all of the power-law galaxies suggests an obscured fraction of <81% (4 : 1).

Gamma-ray luminosity function of gamma-ray bright AGNs

Detection of gamma-ray emissions from a class of active galactic nuclei (viz blazars), has been one of the important findings from the Compton Gamma-Ray Observatory (CGRO). However, their gamma-ray luminosity function has not-been well determined. Few attempts have been made in earlier works, where BL Lacs and Flat Spectrum Radio Quasars (FSRQs) have been considered as a single source class. In this paper, we investigated the evolution and gamma-ray luminosity function of FSRQs and BL Lacs separately. Our investigation indicates no evolution for BL Lacs, however FSRQs show significant evolution. Pure luminosity evolution is assumed for FSRQs and exponential and power law evolution models are examined. Due to the small number of sources, the low luminosity end index of the luminosity function for FSRQs is constrained with an upper limit. BL Lac luminosity function shows no signature of break. As a consistency check, the model source distributions derived from these luminosity functions show no significant departure from the observed source distributions.

Contribution from unresolved discrete sources to the extragalactic gamma-ray background (EGRB)

The origin of the extragalactic gamma-ray background (EGRB) is still an open question, even nearly forty years after its discovery. The emission could originate either from truly diffuse processes or from unresolved point sources. Although the majority of the 271 point sources detected by EGRET (Energetic Gamma Ray Experiment Telescope) are unidentified, of the identified sources, blazars are the dominant candidates. Therefore, unresolved blazars may be considered the main contributor to the EGRB, and many studies have been carried out to understand their distribution, evolution and contribution to the EGRB. Considering that gamma-ray emission comes mostly from jets of blazars and that the jet emission decreases rapidly with increasing jet to line-of-sight angle, it is not surprising that EGRET was not able to detect many large inclination angle active galactic nuclei (AGNs). Though Fermi could only detect a few large inclination angle AGNs during the first three months of its survey, it is expected to detect many such sources in the near future. Since non-blazar AGNs are expected to have higher density as compared to blazars, these could also contribute significantly to the EGRB. In this paper, we discuss contributions from unresolved discrete sources including normal galaxies, starburst galaxies, blazars and off-axis AGNs to the EGRB.

Two-temperature accretion around rotating black holes: a description of the general advective flow paradigm in the presence of various cooling processes to explain low to high luminous sources

We investigate viscous two-temperature accretion disc flows around rotating black holes. We describe the global solution of accretion flows with a sub-Keplerian angular momentum profile, by solving the underlying conservation equations including explicit cooling processes self-consistently. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. We focus on the set of solutions for sub-Eddington, Eddington and super-Eddington mass accretion rates around Schwarzschild and Kerr black holes with a Kerr parameter of 0.998. It is found that the flow, during its infall from the Keplerian to sub-Kepleria transition region to the black hole event horizon, passes through various phases of advection: the general advective paradigm to the radiatively inefficient phase, and vice versa. Hence, the flow governs a much lower electron temperature similar to 10(8)-10(9.5) K, in the range of accretion rate in Eddington units 0.01 less than or similar to (M) over dot less than or similar to 100, compared to the hot protons of temperature similar to 10(10.2)-10(11.8) K. Therefore, the solution may potentially explain the hard X-rays and gamma-rays emitted from active galactic nuclei (AGNs) and X-ray binaries. We then compare the solutions for two different regimes of viscosity. We conclude that a weakly viscous flow is expected to be cooling dominated, particularly at the inner region of the disc, compared to its highly viscous counterpart, which is radiatively inefficient. With all the solutions in hand, we finally reproduce the observed luminosities of the underfed AGNs and quasars (e. g. Sgr A*) to ultraluminous X-ray sources (e. g. SS433), at different combinations of input parameters, such as the mass accretion rate and the ratio of specific heats. The set of solutions also predicts appropriately the luminosity observed in highly luminous AGNs and ultraluminous quasars (e. g. PKS 0743-67).

A possible origin of viscosity in Keplerian accretion disks due to secondary perturbation: Turbulent transport without magnetic fields

The origin of hydrodynamic turbulence in rotating shear flow is a long standing puzzle. Resolving it is especially important in astrophysics when the flow's angular momentum profile is Keplerian which forms an accretion disk having negligible molecular viscosity. Hence, any viscosity in such systems must be due to turbulence, arguably governed by magnetorotational instability, especially when temperature T greater than or similar to 10(5). However, such disks around quiescent cataclysmic variables, protoplanetary and star-forming disks, and the outer regions of disks in active galactic nuclei are practically neutral in charge because of their low temperature, and thus are not expected to be coupled with magnetic fields enough to generate any transport due to the magnetorotational instability. This flow is similar to plane Couette flow including the Coriolis force, at least locally. What drives their turbulence and then transport, when such flows do not exhibit any unstable mode under linear hydrodynamic perturbation? We demonstrate that the three-dimensional secondary disturbance to the primarily perturbed flow that triggers elliptical instability may generate significant turbulent viscosity in the range 0.0001 less than or similar to nu(t) less than or similar to 0.1, which can explain transport in accretion flows.

Secondary Perturbation Effects in Keplerian Accretion Disks: Elliptical Instability

Origin of turbulence in cold accretion disks, particularly in 3D, which is expected to be hydrodynamic but not magnetohydrodynamic, is a big puzzle. While the flow must exhibit some turbulence in support of the transfer of mass inward and angular momentum outward, according to the linear perturbation theory it should always be stable. We demonstrate that the 3D secondary disturbance to the primarily perturbed disk which exhibits elliptical vortices into the system solves the problem. This result is essentially applicable to the outer region of accretion disks in active galactic nuclei where the gas is significantly cold and neutral in charge and the magnetic Reynolds number is smaller than 10^4.

Cause and Effect of Feedback: Multiphase Gas in Cluster Cores Heated by AGN Jets

Multiwavelength data indicate that the X-ray-emitting plasma in the cores of galaxy clusters is not cooling catastrophically. To a large extent, cooling is offset by heating due to active galactic nuclei (AGNs) via jets. The cool-core clusters, with cooler/denser plasmas, show multiphase gas and signs of some cooling in their cores. These observations suggest that the cool core is locally thermally unstable while maintaining global thermal equilibrium. Using high-resolution, three-dimensional simulations we study the formation of multiphase gas in cluster cores heated by collimated bipolar AGN jets. Our key conclusion is that spatially extended multiphase filaments form only when the instantaneous ratio of the thermal instability and free-fall timescales (t(TI)/t(ff)) falls below a critical threshold of approximate to 10. When this happens, dense cold gas decouples from the hot intracluster medium (ICM) phase and generates inhomogeneous and spatially extended Ha filaments. These cold gas clumps and filaments ``rain'' down onto the central regions of the core, forming a cold rotating torus and in part feeding the supermassive black hole. Consequently, the self-regulated feedback enhances AGN heating and the core returns to a higher entropy level with t(TI)/t(ff) > 10. Eventually, the core reaches quasi-stable global thermal equilibrium, and cold filaments condense out of the hot ICM whenever t(TI)/t(ff) less than or similar to 10. This occurs despite the fact that the energy from AGN jets is supplied to the core in a highly anisotropic fashion. The effective spatial redistribution of heat is enabled in part by the turbulent motions in the wake of freely falling cold filaments. Increased AGN activity can locally reverse the cold gas flow, launching cold filamentary gas away from the cluster center. Our criterion for the condensation of spatially extended cold gas is in agreement with observations and previous idealized simulations.

Turbulence and cooling in galaxy cluster cores

We study the interplay between turbulent heating, mixing, and radiative cooling in an idealized model of cool cluster cores. Active galactic nuclei (AGN) jets are expected to drive turbulence and heat cluster cores. Cooling of the intracluster medium (ICM) and stirring by AGN jets are tightly coupled in a feedback loop. We impose the feedback loop by balancing radiative cooling with turbulent heating. In addition to heating the plasma, turbulence also mixes it, suppressing the formation of cold gas at small scales. In this regard, the effect of turbulence is analogous to thermal conduction. For uniform plasma in thermal balance (turbulent heating balancing radiative cooling), cold gas condenses only if the cooling time is shorter than the mixing time. This condition requires the turbulent kinetic energy to be a parts per thousand(3) the plasma internal energy; such high velocities in cool cores are ruled out by observations. The results with realistic magnetic fields and thermal conduction are qualitatively similar to the hydrodynamic simulations. Simulations where the runaway cooling of the cool core is prevented due to mixing with the hot ICM show cold gas even with subsonic turbulence, consistent with observations. Thus, turbulent mixing is the likely mechanism via which AGN jets heat cluster cores. The thermal instability growth rates observed in simulations with turbulence are consistent with the local thermal instability interpretation of cold gas in cluster cores.

COOL CORE CYCLES: COLD GAS AND AGN JET FEEDBACK IN CLUSTER CORES

Using high-resolution 3D and 2D (axisymmetric) hydrodynamic simulations in spherical geometry, we study the evolution of cool cluster cores heated by feedback-driven bipolar active galactic nuclei (AGNs) jets. Condensation of cold gas, and the consequent enhanced accretion, is required for AGN feedback to balance radiative cooling with reasonable efficiencies, and to match the observed cool core properties. A feedback efficiency (mechanical luminosity approximate to epsilon(M) over dot(acc)c(2); where (M) over dot(acc). is the mass accretion rate at 1 kpc) as small as 6 x 10(-5) is sufficient to reduce the cooling/accretion rate by similar to 10 compared to a pure cooling flow in clusters (with M-200 less than or similar to 7 x 10(14) M-circle dot). This value is much smaller compared to the ones considered earlier, and is consistent with the jet efficiency and the fact that only a small fraction of gas at 1 kpc is accreted onto the supermassive black hole (SMBH). The feedback efficiency in earlier works was so high that the cluster core reached equilibrium in a hot state without much precipitation, unlike what is observed in cool-core clusters. We find hysteresis cycles in all our simulations with cold mode feedback: condensation of cold gas when the ratio of the cooling-time to the free-fall time (t(cool)/t(ff)) is less than or similar to 10 leads to a sudden enhancement in the accretion rate; a large accretion rate causes strong jets and overheating of the hot intracluster medium such that t(cool)/t(ff) > 10; further condensation of cold gas is suppressed and the accretion rate falls, leading to slow cooling of the core and condensation of cold gas, restarting the cycle. Therefore, there is a spread in core properties, such as the jet power, accretion rate, for the same value of core entropy t(cool)/t(ff). A smaller number of cycles is observed for higher efficiencies and for lower mass halos because the core is overheated to a longer cooling time. The 3D simulations show the formation of a few-kpc scale, rotationally supported, massive (similar to 10(11) M-circle dot) cold gas torus. Since the torus gas is not accreted onto the SMBH, it is largely decoupled from the feedback cycle. The radially dominant cold gas (T < 5 x 10(4) K; vertical bar v(r)vertical bar >vertical bar v(phi vertical bar)) consists of fast cold gas uplifted by AGN jets and freely infalling cold gas condensing out of the core. The radially dominant cold gas extends out to 25 kpc for the fiducial run (halo mass 7 x 10(14) M-circle dot and feedback efficiency 6 x 10(-5)), with the average mass inflow rate dominating the outflow rate by a factor of approximate to 2. We compare our simulation results with recent observations.

Computational studies of the transverse structure of AGN jets

Both the emission properties and the evolution of the radio jets of Active Galactic Nuclei are dependent on the magnetic (B) fields that thread them. A number of observations of AGN jets suggest that the B fields they carry have a significant helical component, at least on parsec scales. This thesis uses a model, first proposed by Laing and then developed by Papageorgiou, to explore how well the observed properties of AGN jets can be reproduced by assuming a helical B field with three parameters; pitch angle, viewing angle and degree of entanglement. This model has been applied to multifrequency Very Long Baseline Interferometry (VLBI) observations of the AGN jets of Markarian 501 and M87, making it possible to derive values for the helical pitch angle, the viewing angle and the degree of entanglement for these jets. Faraday rotation measurements are another important tool for investigating the B fields of AGN jets. A helical B field component should result in a systematic gradient in the observed Faraday rotation across the jet. Real observed radio images have finite resolution; typical beam sizes for cm-wavelength VLBI observations are often comparable to or larger than the intrinsic jet widths, raising questions about how well resolved a jet must be in the transverse direction in order to reliably detect transverse Faraday-rotation structure. This thesis presents results of Monte Carlo simulations of Faraday rotation images designed to directly investigate this question, together with a detailed investigation into the probabilities of observing spurious Faraday Rotation gradients as a result of random noise and finite resolution. These simulations clearly demonstrate the possibility of detecting transverse Faraday-rotation structures even when the intrinsic jet widths are appreciably smaller than the beam width.

The development of new methods for high resolution radio astronomy imaging

Very Long Baseline Interferometry (VLBI) polarisation observations of the relativistic jets from Active Galactic Nuclei (AGN) allow the magnetic field environment around the jet to be probed. In particular, multi-wavelength observations of AGN jets allow the creation of Faraday rotation measure maps which can be used to gain an insight into the magnetic field component of the jet along the line of sight. Recent polarisation and Faraday rotation measure maps of many AGN show possible evidence for the presence of helical magnetic fields. The detection of such evidence is highly dependent both on the resolution of the images and the quality of the error analysis and statistics used in the detection. This thesis focuses on the development of new methods for high resolution radio astronomy imaging in both of these areas. An implementation of the Maximum Entropy Method (MEM) suitable for multi-wavelength VLBI polarisation observations is presented and the advantage in resolution it possesses over the CLEAN algorithm is discussed and demonstrated using Monte Carlo simulations. This new polarisation MEM code has been applied to multi-wavelength imaging of the Active Galactic Nuclei 0716+714, Mrk 501 and 1633+382, in each case providing improved polarisation imaging compared to the case of deconvolution using the standard CLEAN algorithm. The first MEM-based fractional polarisation and Faraday-rotation VLBI images are presented, using these sources as examples. Recent detections of gradients in Faraday rotation measure are presented, including an observation of a reversal in the direction of a gradient further along a jet. Simulated observations confirming the observability of such a phenomenon are conducted, and possible explanations for a reversal in the direction of the Faraday rotation measure gradient are discussed. These results were originally published in Mahmud et al. (2013). Finally, a new error model for the CLEAN algorithm is developed which takes into account correlation between neighbouring pixels. Comparison of error maps calculated using this new model and Monte Carlo maps show striking similarities when the sources considered are well resolved, indicating that the method is correctly reproducing at least some component of the overall uncertainty in the images. The calculation of many useful quantities using this model is demonstrated and the advantages it poses over traditional single pixel calculations is illustrated. The limitations of the model as revealed by Monte Carlo simulations are also discussed; unfortunately, the error model does not work well when applied to compact regions of emission.

Synthetic X-ray spectra for simulations of the dynamics of an accretion flow irradiated by a quasar

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. It has been proposed that some of these outflows are driven off gas infalling towards the central supermassive black hole. We perform radiative transfer calculations to compute the gas ionization state and the emergent X-ray spectra for both two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of relatively dense (10 -10 g cm ) outflowing gas. This outflow causes strong absorption for observer orientation angles of ?35°. Particularly in 3D, the strength of this absorption varies significantly for different lines of sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features in the X-ray spectra since the ionization state is predicted to be very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a Compton-scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow are predicted to give rise to several emission features including a multicomponent Fe Ka emission complex for all observer orientations. For the hydrodynamical simulations considered, we predict both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. Higher resolution hydrodynamic simulations are needed to determine whether the outflows fragment on scales unresolved in our current study, which may yield the denser lower ionization material that could reconcile the models and the observations. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.