Boxy, peanut- or X-shaped "bulges" are observed in a large fraction of barred galaxies viewed in, or close to, edge-on projection, as well as in the Milky Way. They are the product of dynamical instabilities occurring in stellar bars, which cause the latter to buckle and thicken vertically. Recent studies have found nearby galaxies that harbour two such features arising at different radial scales, in a nested configuration. In this paper we explore the formation of such double peanuts, using a collisionless N-body simulation of a pure disc evolving in isolation within a live dark matter halo, which we analyse in a completely analogous way to observations of real galaxies. In the simulation we find a stable double configuration consisting of two X/peanut structures associated to the same galactic bar - rotating with the same pattern speed - but with different morphology, formation time, and evolution. The inner, conventional peanut-shaped structure forms early via the buckling of the bar, and experiences little evolution once it stabilises. This feature is consistent in terms of size, strength and morphology, with peanut structures observed in nearby galaxies. The outer structure, however, displays a strong X, or "bow-tie", morphology. It forms just after the inner peanut, and gradually extends in time (within 1 to 1.5 Gyr) to almost the end of the bar, a radial scale where ansae occur. We conclude that, although both structures form, and are dynamically coupled to, the same bar, they are supported by inherently different mechanisms.

In the MW bulge, metal-rich stars form a strong bar and are more peanut-shaped than metal-poor stars. It has been recently claimed that this behavior is driven by the initial in-plane radial velocity dispersion of these populations, rather than by their initial vertical random motions. This has led to the suggestion that a thick disc is not necessary to explain the characteristics of the MW bulge. We rediscuss this issue by analyzing two dissipationless N-body simulations of boxy/peanut (b/p)-shaped bulges formed from composite stellar discs, made of kinematically cold and hot stellar populations, and we conclude that initial vertical random motions are as important as in-plane random motions in determining the relative contribution of cold and hot disc populations with height above the plane, the metallicity and age trends. Previous statements emphasizing the dominant role of in-plane motions in determining these trends are not confirmed. However, differences exist in the morphology and strength of the resulting b/p-shaped bulges: a model where disc populations have initially only different in-plane random motions, but similar thickness, results into a b/p bulge where all populations have a similar peanut shape, independently on their initial kinematics, or metallicity. We discuss the reasons behind these differences, and also predict the signatures that these two extreme initial conditions would leave on the vertical age and metallicity gradients of disc stars, outside the bulge region. We conclude that a metal-poor, kinematically (radial and vertical) hot component, that is a thick disc, is necessary in the MW before bar formation, supporting the scenario traced in previous works. [abridged]

Type-III S0 galaxies present tight scaling relations between their surface brightness photometric and structural parameters. Several evolutionary models have been proposed for the formation of Type-III S0 galaxies but the observations of are usually limited to the local Universe. We study the evolution of the photometric and structural scaling relations found between the parameters of the surface brightness profiles in the rest-frame R-band of Type-III S0 galaxies with z and the possible differences between the rest-frame (B-R) colours of the inner and outer disc profiles. We make use of a sample of 14 Type-III E/S0--S0 galaxies at 0.2<z<0.6 to study if the correlations found in local Type-III S0 galaxies were present ~6 Gyr ago. We analyse the distribution of the surface brightness characteristic parameters as a function of the stellar mass and if there is a significant change with z. We also derive their rest-frame (B-R) colour profiles and we compare these results with the predictions from a grid of SSP models. We find that the inner and outer scale-lengths of Type-III S0 galaxies at 0.4<z<0.6 follow compatible trends and scaling relations with those observed in local S0 galaxies. We do not detect any significant differences between the location of Rbreak between z~0.6 and z=0 for a fixed stellar mass of the object, whereas the surface brightness at the break radius is ~1.5 mag arcsec-2 dimmer in the local Universe than at z~0.6 for a fixed stellar mass. In contrast to Type-II profiles, the Type-III surface brightness profiles of S0 galaxies present compatible Rbreak values and scaling relations during the last 6 Gyr. This result and the similarity of the colours of the inner and outer discs point to a highly scalable and stable formation process, probably more related to gravitational and dynamical processes than to the evolution of stellar populations (abridged).

The nature, size and orientation of the dominant structural components in the Milky Way's inner ~4 kpc - specifically the bulge and bar - have been the subject of conflicting interpretations in the literature. We present a different approach to inferring the properties of the long bar which extends beyond the inner bulge, via the information encoded in the Galaxy's X/peanut (X/P)-shaped structure. We perform a quantitative analysis of the X/P feature seen in WISE wide-field imaging at 3.4 $\mu$m and 4.6 $\mu$m. We measure the deviations of the isophotes from pure ellipses, and quantify the X/P structure via the radial profile of the Fourier n=6 harmonic (cosine term $B_6$). In addition to the vertical height and integrated `strength' of the X/P instability, we report an intrinsic radius of $R_{\Pi;{\rm int}} = 1.67\pm0.27$ kpc, and an orientation angle of $\alpha = {37^\circ}^{+7^\circ}_{-10^\circ}$ with respect to our line-of-sight to the Galactic Centre. Based on X/P structures observed in other galaxies, we make three assumptions: (i) the peanut is intrinsically symmetric, (ii) the peanut is aligned with the long Galactic bar, and (iii) their sizes are correlated. Thus the implication for the Galactic bar is that it is oriented at the same $37^\circ$ angle and has an expected radius of ~4.2 kpc, but possibly as low as ~3.2 kpc. We further investigate how the Milky Way's X/P structure compares with other analogues, and find that the Galaxy is broadly consistent with our recently established scaling relations, though with a moderately stronger peanut instability than expected. We additionally perform a photometric decomposition of the Milky Way's major axis surface brightness profile, accounting for spiral structure, and determine an average disc scale length of $h = 2.54\pm0.16$ kpc in the WISE bands, in good agreement with the literature.

The controversy about the origin of the structure of S0--E/S0 galaxies may be due to the difficulty of comparing surface brightness profiles with different depths, photometric corrections and PSF effects (almost always ignored). We aim to quantify the properties of Type-III (anti-truncated) discs in a sample of S0 galaxies at 0.2<z<0.6. In this paper, we present the sample selection and describe in detail the methods to robustly trace the structure in their outskirts and correct for PSF effects. We have selected and classified a sample of 150 quiescent galaxies at 0.2<z<0.6 in the GOODS-N field. We perform a quantitative structural analysis of 44 S0-E/S0 galaxies. We corrected their surface brightness profiles for PSF distortions and analysed the biases in the structural and photometric parameters when the PSF correction is not applied. Additionally, we have developed Elbow, an automatic statistical method to determine whether a possible break is significant - or not - and its type and made it publicly available. We found 14 anti-truncated S0-E/S0 galaxies in the range 0.2<z<0.6 (~30% of the final sample). This fraction is similar to the those reported in the local Universe. In our sample, ~25% of the Type-III breaks observed in PSF-uncorrected profiles are artifacts, and their profiles turn into a Type I after PSF correction. PSF effects also soften Type-II profiles. We found that the profiles of Type-I S0 and E/S0 galaxies of our sample are compatible with the inner profiles of the Type-III, in contrast with the outer profiles. We have obtained the first robust and reliable sample of 14 anti-truncated S0--E/S0 galaxies beyond the local Universe, in the range 0.2<z<0.6. PSF effects significantly affect the shape of the surface brightness profiles in galaxy discs even in the case of the narrow PSF of HST/ACS images, so future studies on the subject should make an effort to correct them.

Selected from a sample of nine, isolated, dwarf early-type galaxies (ETGs) having the same range of kinematic properties as dwarf ETGs in clusters, we use CG 611 (LEDA 2108986) to address the Nature versus Nurture debate regarding the formation of dwarf ETGs. The presence of faint disk structures and rotation within some cluster dwarf ETGs has often been heralded as evidence that they were once late-type spiral or dwarf irregular galaxies prior to experiencing a cluster-induced transformation into an ETG. However, CG 611 also contains significant stellar rotation (~20 km/s) over its inner half light radius, R_(e,maj)=0.71 kpc, and its stellar structure and kinematics resemble those of cluster ETGs. In addition to hosting a faint young nuclear spiral within a possible intermediate-scale stellar disk, CG 611 has accreted an intermediate-scale, counter-rotating gas disk. It is therefore apparent that dwarf ETGs can be built by accretion events, as opposed to disk-stripping scenarios. We go on to discuss how both dwarf and ordinary ETGs with intermediate-scale disks, whether under (de)construction or not, are not fully represented by the kinematic scaling S_0.5=sqrt{ 0.5(V_rot)^2 + sigma^2 }, and we also introduce a modified spin-ellipticity diagram, lambda(R)-epsilon(R), with the potential to track galaxies with such disks.

Aims. ESO 243-49 HLX-1, otherwise known as HLX-1, is an intermediate mass black hole (IMBH) candidate located 8" (3.7 Kpc) from the centre of the edge-on S0 galaxy ESO 243-49. How the black hole came to be associated with this galaxy, and the nature of the environment in which it resides, are still unclear. Using multi-wavelength observations we investigate the nature of the medium surrounding HLX-1, search for evidence of past mergers with ESO 243-49 and constrain parameters of the galaxy. Methods. We reduce and analyse integral field unit observations of ESO 243-49 that were taken with the MUSE instrument on the VLT. Using complementary multi-wavelength data, including X-Shooter, HST, Swift, Chandra and ATCA data, we further examine the vicinity of HLX-1. We additionally examine the nature of the host galaxy and estimate the mass of the central supermassive black hole (SMBH) in ESO 243-49. Results. No evidence for a recent minor-merger that could result in the presence of the IMBH is discerned, but the data are compatible with a scenario in which minor mergers may have occurred in the history of ESO 243-49. The MUSE data reveal a rapidly rotating disc in the centre of the galaxy, around the SMBH. The mass of the SMBH at the centre of ESO 243-49 is estimated to be 0.5-23 $\times$ 10$^7$ M$_\odot$. Studying the spectra of HLX-1, that were taken in the low/hard state, we determine H$_\alpha$ flux variability to be at least a factor 6, compared to observations taken during the high/soft state. This H$_\alpha$ flux variability over one year indicates that the line originates close to the IMBH, excluding the possibility that the line emanates from a surrounding nebula or a star cluster. The large variability associated with the X-ray states of HLX-1 confirms that the H$_\alpha$ line is associated with the object and therefore validates the distance to HLX-1.

Aims: This work is the first stage of a campaign to search for IMBHs, in low luminosity AGN (LLAGN) and dwarf galaxies. An additional, and equally important, aim of this pilot study is to investigate the consistency between the predictions of several BH scaling relations and the fundamental plane of black hole activity (FP-BH). Methods: We use X-ray and radio luminosity relations in accreting BHs, along with the latest scaling relations between the mass of the central black hole (MBH) and the properties of its host spheroid, to predict MBH in seven LLAGN, that were previously reported to be in the IMBH regime. Namely, we use the recently re-evaluated MBH - Msph (Msph : spheroid absolute magnitude at 3.6 microm) scaling relation for spiral galaxies, the MBH - nsph (nsph : major axis Sersic index of the bulge) relation, the MBH - PA (PA: pitch angle) relation and the FP-BH for weakly accreting BHs, to independently estimate MBH in all seven galaxies. Results: We find that all LLAGN in our list have low-mass central black holes - with log MBH/Msol~6.5 on average - but they are, not IMBHs. All 4 methods used predicted consistent BH masses, in the 1sigma range. Furthermore, we report that, in contrast to previous classification, galaxy NGC 4470 is bulge-less, and cast doubts on the AGN classification of NGC 3507. Conclusions: We find that our latest, state-of-the-art techniques for bulge magnitude & Sersic index computations and the most recent updates of the MBH - Lsph, MBH - nsph , MBH - PA relations and the FP-BH produce consistent results in the low mass regime. We establish a multiple method approach for predicting BH masses in the regime where their spheres-of- grav. influence cannot be spatially resolved. Our approach mitigates against outliers from any one relation and provides a more robust average prediction. We will use our new method to revisit more IMBH candidates in LLAGN

While spiral and lenticular galaxies have large-scale disks extending beyond their bulges, and most local early-type galaxies with 10^{10} < M_*/M_Sun < 2x10^{11} contain a disk (e.g., ATLAS^3D), the early-type galaxies do possess a range of disk sizes. The edge-on, `intermediate-scale' disk in the `disky elliptical' galaxy NGC 1271 has led to some uncertainty as to what its spheroidal component is. Walsh et al. reported a directly measured black hole mass of 3x10^9 M_Sun for this galaxy; which they remarked was an order of magnitude greater than what they expected based on their derivation of the host spheroid's luminosity. Our near-infrared image analysis supports a small embedded disk within a massive spheroidal component with M_{sph,*} = (0.9+/-0.2)x10^{11} M_Sun (using M_*/L_H = 1.4 from Walsh et al.). This places NGC 1271 just 1.6-sigma above the near-linear M_bh-M_{sph,*} relation for early-type galaxies. Therefore, past speculation that there may be a systematic difference in the black hole scaling relations between compact massive early-type galaxies with intermediate-scale disks, i.e. ES galaxies such as NGC 1271, and early-type galaxies with either no substantial disk (E) or a large-scale disk (S0) is not strongly supported by NGC 1271. We additionally (i) show how ES galaxies fit naturally in the (`bulge'-to-total)-(morphological type) diagram, while noting a complication with recent revisions to the Hubble-Jeans tuning-fork diagram, (ii) caution about claims of over-massive black holes in other ES galaxies if incorrectly modelled as S0 galaxies, and (iii) reveal that the compact massive spheroid in NGC 1271 has properties similar to bright bulges in other galaxies which have grown larger-scale disks.

I introduce $Profiler$, a new, user-friendly program written in $Python$ and designed to analyse the radial surface brightness profiles of galaxies. With an intuitive graphical user interface, $Profiler$ can accurately model a wide range of galaxies and galaxy components, such as elliptical galaxies, the bulges of spiral and lenticular galaxies, nuclear sources, discs, bars, rings, spiral arms, etc., with a variety of parametric functions routinely employed in the field (S\'ersic, core-S\'ersic, exponential, Gaussian, Moffat and Ferrers). In addition to these, $Profiler$ can employ the broken exponential model (relevant for disc truncations or antitruncations) and two special cases of the edge-on disc model: namely along the major axis (in the disc plane) and along the minor axis (perpendicular to the disc plane). $Profiler$ is optimised to work with galaxy light profiles obtained from isophotal measurements which capture radial gradients in the ellipticity, position angle and Fourier harmonic profiles of the isophotes, and are thus often better at capturing the total light than two-dimensional image-fitting programs. Additionally, the one-dimensional approach is generally less computationally expensive and more stable. In $Profiler$, the convolution of either circular or elliptical models with the point spread function is performed in two-dimensions, and offers a choice between Gaussian, Moffat or a user-provided data vector (a table of intensity values as a function of radius) for the point spread function. I demonstrate $Profiler$'s features and operation by decomposing three case-study galaxies: the cored elliptical galaxy NGC 3348, the nucleated dwarf Seyfert I galaxy Pox 52, and NGC 2549, a structurally complex, double-barred galaxy which also displays a Type II truncated disc viewed edge-on. $Profiler$ is freely available at: https://github.com/BogdanCiambur/PROFILER

X-shaped or peanut-shaped (X/P) bulges are observed in more than 40% of (nearly) edge-on disc galaxies, though to date a robust method to quantify them is lacking. Using Fourier harmonics to describe the deviation of galaxy isophotes from ellipses, we demonstrate with a sample of 11 such galaxies (including NGC 128) that the sixth Fourier component ($B_6$) carries physical meaning by tracing this X/P structure. We introduce five quantitative diagnostics based on the radial $B_6$ profile, namely: its `peak' amplitude ($\Pi_{\rm max}$); the (projected major-axis) `length' where this peak occurs ($R_{\Pi,max}$); its vertical `height' above the disc plane ($z_{\Pi,max}$); the $B_6$ profile's integrated `strength' ($S_{\Pi}$); and the $B_6$ peak `width' ($W_{\Pi}$). We also introduce different `classes' of $B_6$ profile shape. Furthermore, we convincingly detect and measure the properties of multiple (nested) X/P structures in individual galaxies which additionally display the signatures of multiple bars in their surface brightness profiles, thus consolidating further the scenario in which peanuts are associated with bars. We reveal that the peanut parameter space (`length', `strength' and `height') for real galaxies is not randomly populated, but the 3 metrics are inter-correlated (both in kpc and disc scale-length $h$). Additionally, the X/P `length' and `strength' appear to correlate with $v_{rot}/\sigma_*$, lending further support to the notion that peanuts `know' about the galactic disc in which they reside. Such constraints are important for simulations, as they provide a direct link between peanuts and their host disc. Our diagnostics reveal a spectrum of X/P properties and could provide a means of distinguishing between different peanut formation scenarios. Moreover, nested peanuts, as remnants of bar buckling events, can provide insights into the disc and bar instability history.

The mass scaling relation between supermassive black holes and their host spheroids has previously been described by a quadratic or steeper relation at low masses (10^5 < M_bh/M_sun < 10^7). How this extends into the realm of intermediate mass black holes (10^2 < M_bh/M_sun < 10^5) is not yet clear, although for the barred Sm galaxy LEDA 87300, Baldassare et al. have recently reported a nominal virial mass M_bh=5x10^4 M_sun residing in a `spheroid' of stellar mass equal to 6.3x10^8 M_sun. We point out, for the first time, that LEDA 87300 therefore appears to reside on the near-quadratic M_bh-M_{sph,*} relation. However, Baldassare et al. modelled the bulge _and_ bar as the single spheroidal component of this galaxy. Here we perform a 3-component bulge+bar+disk decomposition and find a bulge luminosity which is 7.7 times fainter than the published `bulge' luminosity. After correcting for dust we find that M_bulge=0.9x10^8 M_sun, and M_bulge/M_disk=0.04 - which is now in accord with ratios typically found in Scd-Sm galaxies. We go on to discuss slight revisions to the stellar velocity dispersion (40+/-11 km/s) and black hole mass (M_bh=2.9x10^4 M_sun) and show that LEDA 87300 remains consistent with both the M_bh-sigma relation and the near-quadratic M_bh-M_{sph,*} relation when using the reduced bulge mass. LEDA 87300 therefore offers the first support for the rapid but regulated (near-quadratic) growth of black holes, relative to their host bulge/spheroid, extending into the domain of intermediate mass black holes.

This work introduces a new fitting formalism for isophotes which enables more accurate modelling of galaxies with non-elliptical shapes, such as disk galaxies viewed edge-on or galaxies with X-shaped/peanut bulges. Within this scheme, the angular parameter which defines quasi-elliptical isophotes is transformed from the commonly used, but inappropriate, polar co-ordinate to the `eccentric anomaly'. This provides a superior description of deviations from ellipticity, better capturing the true isophotal shape. Furthermore, this makes it possible to accurately recover both the surface brightness profile, using the correct azimuthally-averaged isophote, and the two-dimensional model of any galaxy: the hitherto ubiquitous, but artificial, cross-like features in residual images are completely removed. The formalism has been implemented into the IRAF tasks $Ellipse$ and $Bmodel$ to create the new tasks `$Isofit$', and `$Cmodel$'. The new tools are demonstrated here with application to five galaxies, chosen to be representative case-studies for several areas where this technique makes it possible to gain new scientific insight. Specifically: properly quantifying boxy/disky isophotes via the fourth harmonic order in edge-on galaxies, quantifying X-shaped/peanut bulges, higher-order Fourier moments for modelling bars in disks, and complex isophote shapes. Higher order (n > 4) harmonics now become meaningful and may correlate with structural properties, as boxyness/diskyness is known to do. This work also illustrates how the accurate construction, and subtraction, of a model from a galaxy image facilitates the identification and recovery of over-lapping sources such as globular clusters and the optical counterparts of X-ray sources.

We have compared the observed distribution of galaxies in the plane of specific star formation rate versus stellar mass with the predictions of the Garching semi-analytic model at redshifts 0, 1 and 2. The goal is to test whether the implementation of radio mode AGN feedback, which is responsible for terminating the formation of stars in high mass galaxies, provides an adequate match to current high-redshift observations. The fraction of quenched galaxies as a function of stellar mass in the models is in good agreement with data at z=0 and z=1. By z=2, there are too few quenched galaxies with low stellar masses in the models. At z=2, the population of galaxies with no ongoing star formation is clearly separated from the `main sequence' of star-forming galaxies in the data. This is not found in the models, because z=2 galaxies with stellar masses less than 10^11 solar masses are predicted to host black holes with relatively low masses (less than 10^8 solar masses). The current implementation of radio mode feedback from such black holes reduces the cooling rates from the surrounding halo, but does not generate sufficient energy to stop star formation entirely. We suggest that the models may be brought into better agreement with the data if black hole growth is triggered by disc instabilities in addition to major mergers, and if feedback mechanisms associated with the formation of galactic bulges act to quench star formation in galaxies.