Reprocessed emission line profiles from dense clouds in geometrically thick accretion engines

The central engines of active galactic nuclei (AGN) contain cold, dense material as well as hot X-ray-emitting gas. The standard paradigm for the engine geometry is a cold thin disc sandwiched between hot X-ray coronae. Strong support for this geometry in Seyferts comes from the study of fluorescent iron line profiles, although the evidence is not ubiquitously airtight. The thin disc model of line profiles in AGN and in X-ray binaries should still be benchmarked against other plausible possibilities. One proposed alternative is an engine consisting of dense clouds embedded in an optically thin, geometrically thick X-ray-emitting engine. This model is also motivated by studies of geometrically thick engines such as advection-dominated accretion flows (ADAFs). Here we compute the reprocessed iron line profiles from dense clouds embedded in geometrically thick, optically thin X-ray-emitting discs near a Schwarzschild black hole. We consider a range of cloud distributions and disc solutions, including ADAFs, pure radial infall and bipolar outflows. We find that such models can reproduce line profiles similar to those from geometrically thin, optically thick discs and might help alleviate some of the problems encountered from the latter. Thus, independent of thin discs, thick disc engines can also exhibit iron line profiles if embedded dense clouds can survive long enough to reprocess radiation.     

Simulations of minor mergers – II. The phase-space structure of thick discs

We analyse the phase-space structure of simulated thick discs that are the result of a 5:1 mass-ratio merger between a disc galaxy and a satellite. Our main goal is to establish what would be the imprints of a merger origin for the Galactic thick disc. We find that the spatial distribution predicted for thick-disc stars is asymmetric, seemingly in agreement with recent observations of the Milky Way thick disc. Near the Sun, the accreted stars are expected to rotate more slowly, to have broad velocity distributions and to occupy preferentially the wings of the line-of-sight velocity distributions. The majority of the stars in our model thick discs have low eccentricity orbits (in clear reference to the pre-existing heated disc) which give rise to a characteristic (sinusoidal) pattern for their line-of-sight velocities as a function of galactic longitude. The z -component of the angular momentum of thick-disc stars provides a clear discriminant between stars from the pre-existing disc and those from the satellite, particularly at large radii. These results are robust against the particular choices of initial conditions made in our simulations.     

0422−476: a shell galaxy with azimuthally distributed shells

A morphological and two-colour charge-coupled device photometry study of the shell galaxy 0422−476, one of the richest known azimuthally distributed shell galaxies, is presented. Taking this galaxy as a prototype, a general method for reducing observational data of these objects is defined and quantitative parameters for use in further theoretical studies are derived.
According to some recent models (e.g. that of Thomson and Wright in 1990), the shells in such a galaxy could be density waves induced in a thick disc population of dynamically cold stars by a weak interaction with another galaxy. In 0422−476 there is no evidence of either a conventional exponential disc or a thick disc. Although it is not possible to rule out the weak interaction model, the observations continue to favour the merger model (e.g. that of Quinn from 1984).     

Stability of general relativistic Miyamoto–Nagai galaxies

The stability of a recently proposed general relativistic model of galaxies is studied in some detail. This model is a general relativistic version of the well-known Miyamoto–Nagai model that represents well a thick galactic disc. The stability of the disc is investigated under a general first-order perturbation keeping the space–time metric frozen (no gravitational radiation is taken into account). We find that the stability is associated with the thickness of the disc. We find that flat galaxies have more non-stable modes than the thick ones, i.e. flat galaxies have a tendency to form more complex structures like rings, bars and spiral arms.     

Long-term optical and X-ray variability of the Seyfert galaxy Markarian 79

We present the results of concurrent X-ray and optical monitoring of the Seyfert 1 galaxy Mrk 79 over a period of more than 5 yr. We find that on short to medium time-scales (days to a few tens of days) the 2–10 keV X-ray and optical u - and V -band fluxes are significantly correlated, with a delay between the bands consistent with 0 d. We show that most of these variations may be well reproduced by a model where the short-term optical variations originate from reprocessing of X-rays by an optically thick accretion disc. The optical light curves, however, also display long time-scale variations over thousands of days, which are not present in the X-ray light curve. These optical variations must originate from an independent variability mechanism and we show that they can be produced by variations in the (geometrically) thin disc accretion rate as well as by varying reprocessed fractions through changes in the location of the X-ray corona.     

An irradiated accretion disc in the narrow-line Seyfert 1 RE J1034+396?

We model the optical to X-ray continuum spectrum of the narrow-line Seyfert 1 galaxy RE J1034+396. We show that the flat optical spectrum is consistent with emission from an irradiated accretion disc. The X-ray emission can be modelled with a disc blackbody and a Comptonized component. The temperature at the inner edge of the disc     Using this constraint, we show that the flat optical spectrum is consistent with emission from the irradiatively heated outer part of the accretion disc. We constrain the outer radius of the optically thick disc     and the inner radius of the irradiation-dominated region     . Our optical and X-ray spectral fits indicate a mass     , and do not rule out a low (i.e. face-on) inclination angle for the system.     

Radial scalelengths of the galactic thin and thick disc with 2MASS data

This paper presents a global analysis of the 2MASS (Two Micron All Sky Survey) data as observed in seven fields at different galactic latitudes in our Galaxy. The data allow the preliminary determination of the scale parameters, which lead to strong constraints on the radial and vertical structure of the galactic thin and thick disc. The interpretation of star counts and colour distributions of stars in the near-infrared with the synthetic stellar population model gives strong evidence that the galactic thin disc density scalelength ( h _R ) is rather short (2.8±0.3 kpc). The galactic thick disc population is revisited in the light of new data. We find the thick disc to have a local density of 3.5±2.0 per cent of the thin disc, exponential scaleheight ( h _z ) of 860±200 pc and exponential scalelength ( h _R ) of 3.7±_0.5^0.8 kpc.     

Disc–halo interaction – I. Three-dimensional evolution of the Galactic disc

The results of a three-dimensional model for disc–halo interaction are presented here. The model considers explicitly the input of energy and mass by isolated and correlated supernovae in the disc. Once disrupted by the explosions, the disc never returns to its initial state. Instead it approaches a state where a thin H  i disc is formed in the Galactic plane, overlaid by thick H  i and H  ii gas discs with scaleheights of 500 pc and 1–1.5 kpc, respectively. The upper parts of the thick H  ii disc (the diffuse ionized medium) act as a disc–halo interface, and its formation and stability are directly correlated to the supernova rate per unit area in the simulated disc.     

Star formation and chemical evolution of damped Lyman α systems

In this paper, we investigate the star formation and chemical evolution of damped Lyman α systems (DLAs) based on the disc galaxy formation model developed by Mo, Mao & White. We propose that the DLAs are the central galaxies of less-massive dark haloes present at redshifts z ∼3, and they should inhabit haloes of moderately low circular velocity. The empirical Schmidt law of star formation rates, and closed box model of chemical evolution that an approximation known as instantaneous recycling is assumed, are adopted. In our models, when the predicted distribution of metallicity for DLAs is calculated, two cases are considered. One is that, using the closed-box model, empirical Schmidt law and star formation time, the distribution of metallicity can be directly calculated. The other is that, when the simple gravitational instability of a thin isothermal gas disc as first discussed by Toomre is considered, the star formation occurs only in the region where the surface density of gas satisfies the critical value, not everywhere of a gas disc. In this case, we first obtain the region where the star formation can occur by assuming that the disc has a flat rotation curve and rotational velocity is equal to the circular velocity of the surrounding dark matter halo, and then calculate the metallicity distribution as in case one. We assume that star formation in each DLA lasts for a period of 1 Gyr from redshifts z =3. There is only one output parameter in our models, i.e. the stellar yield, which relates to the time of star formation history and is obtained by normalizing the predicted distribution of metallicity to the mean value of 1/13 Z_⊙ as presented by Pettini et al.. The predicted metallicity distribution is consistent with the current (rather limited) observational data. A random distribution of galactic discs is taken into account.     

Chemical evolution with radial mixing

Models of the chemical evolution of our Galaxy are extended to include radial migration of stars and flow of gas through the disc. The models track the production of both iron and α-elements. A model is chosen that provides an excellent fit to the metallicity distribution of stars in the Geneva–Copenhagen survey (GCS) of the solar neighbourhood and a good fit to the local Hess diagram. The model provides a good fit to the distribution of GCS stars in the age–metallicity plane, although this plane was not used in the fitting process. Although this model's star formation rate is monotonically declining, its disc naturally splits into an α-enhanced thick disc and a normal thin disc. In particular, the model's distribution of stars in the ([O/Fe], [Fe/H]) plane resembles that of Galactic stars in displaying a ridge line for each disc. The thin-disc's ridge line is entirely due to stellar migration, and there is the characteristic variation of stellar angular momentum along it that has been noted by Haywood in survey data. Radial mixing of stellar populations with high  σ _z   from inner regions of the disc to the solar neighbourhood provides a natural explanation of why measurements yield a steeper increase of  σ _z   with age than predicted by theory. The metallicity gradient in the interstellar medium is predicted to be steeper than in earlier models, but appears to be in good agreement with data for both our Galaxy and external galaxies. The models are inconsistent with a cut-off in the star formation rate at low gas surface densities. The absolute magnitude of the disc is given as a function of time in several photometric bands, and radial colour profiles are plotted for representative times.     

A grid of chemical evolution models as a tool to interpret spiral and irregular galaxies data

We present a generalization of the multiphase chemical evolution model (CEM) applied to a wide set of theoretical galaxies with different masses and evolutionary rates. This generalized set of models has been computed using the so-called universal rotation curve from Persic, Salucci & Steel to calculate the radial mass distribution of 44 theoretical protogalaxies. This distribution is a fundamental input which, besides its own effect on the galaxy evolution, defines the characteristic collapse time-scale or gas infall rate on to the disc. We have adopted 10 sets of values, between 0 and 1, for the molecular cloud and star formation efficiencies, as corresponding to their probability nature, for each one of the radial distributions of total mass. Thus, we have constructed a biparametric grid of models, depending on those efficiency sets and on the rotation velocity, whose results are valid in principle for any spiral or irregular galaxy. The model results provide the time-evolution of different regions of the disc and the halo along galactocentric distance, measured by the gas (atomic and molecular) and stellar masses, the star formation rate (SFR) and chemical abundances of 14 elements, for a total of 440 models. This grid may be used to estimate the evolution of a given galaxy for which only present time information, such as radial distributions of elemental abundances, gas densities and/or star formation, which are the usual observational constraints of chemical evolution models (CEMs), is available.     

Contrasting the chemical evolution of the Milky Way and Andromeda

The chemical evolution history of a galaxy hides clues about how it formed and has been changing through time. We have studied the chemical evolution history of the Milky Way (MW) and Andromeda (M31) to find which are common features in the chemical evolution of disc galaxies as well as which are galaxy-dependent. We use a semi-analytic multizone chemical evolution model. Such models have succeeded in explaining the mean trends of the observed chemical properties in these two Local Group spiral galaxies with similar mass and morphology. Our results suggest that while the evolution of the MW and M31 shares general similarities, differences in the formation history are required to explain the observations in detail. In particular, we found that the observed higher metallicity in the M31 halo can be explained by either (i) a higher halo star formation efficiency (SFE), or (ii) a larger reservoir of infalling halo gas with a longer halo formation phase. These two different pictures would lead to (i) a higher [O/Fe] at low metallicities, or (ii) younger stellar populations in the M31 halo, respectively. Both pictures result in a more massive stellar halo in M31, which suggests a possible correlation between the halo metallicity and its stellar mass.     

A dimensional study of disc galaxies

We present a highly simplified model of the dynamical structure of a disc galaxy where only two parameters fully determine the solution, mass and angular momentum. We show through simple physical scalings that once the mass has been fixed, the angular momentum parameter λ is expected to regulate such critical galactic disc properties as colour, thickness of the disc and bulge-to-disc ratio. It is, hence, expected to be the determinant physical ingredient resulting in a given Hubble type. A simple analytic estimate of λ for an observed system is provided. An explicit comparison of the distribution of several galactic parameters against both Hubble type and λ is performed using observed galaxies. Both such distributions exhibit highly similar characteristics for all galactic properties studied, suggesting λ as a physically motivated classification parameter for disc galaxies.     

Narrow‐Line Seyfert 1 Galaxies with high‐energy sharp spectral drops: reflection‐versus absorption models

With the launch of XMM‐Newton in 1999, two Narrow‐Line Seyfert 1 Galaxies (NLS1s) have been detected (IRAS 13224–3809 and 1H 0707–495) showing sharp spectral drops at energies equal or above the neutral Fe K edge at 7.1 keV without any narrow Fe K reemission. In this paper I summarize our present knowledge on the observed properties of sharp high‐energy spectral drops. I list the problems presently arising from the reflection dominated and the optically thick disc models. Finally, I present an alternative solution which consists of a combination of the accretion disc model and the reflection dominated model. This might solve the problems of the standard accretion disc model. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Three-dimensional evolution of worms and chimneys in the Galactic disc

Three-dimensional modelling of the flow of gas and plasma in a section of the Galaxy has been carried out to study the evolution and formation of Galactic chimneys and worms. It is found that clustered supernovae located on either side of the Galactic plane are sources for the formation of well-collimated chimneys, having widths of ∼     . The thick gas disc may have a role in the collimation of chimneys. Channel maps of disc gas, obtained from the simulations, show the presence of sheet-like structures running perpendicular to the Galactic plane and resembling worms. Worms are believed to result from the break-up of the shells and supershells. However, the simulations show that although some worms correlate well with the debris of broken shells/supershells, others do not. They are cold gas that has been accelerated in the disc and rise on to the thick gas disc.     

Chemo-spectrophotometric evolution of spiral galaxies – III. Abundance and colour gradients in discs

We study the relations between luminosity and chemical-abundance profiles of spiral galaxies, using detailed models for the chemical and spectrophotometric evolution of galactic discs. The models are 'calibrated' on the Milky Way disc and are successfully extended to other discs with the help of simple 'scaling' relations, obtained in the framework of semi-analytic models of galaxy formation. We find that our models exhibit oxygen abundance gradients that increase in absolute value with decreasing disc luminosity (when expressed in dex kpc⁻¹) and are independent of disc luminosity (when expressed in dex scalelength⁻¹), both in agreement with observations. We notice an important strong correlation between abundance gradient and disc scalelength. These results support the idea of 'homologous evolution' of galactic discs.     

On the geometry of broad emission region in quasars

We study the geometry of the Hβ broad emission region by comparing the M _BH values derived from Hβ through the virial relation with those obtained from the host galaxy luminosity in a sample of 36 low-redshift  ( z ∼ 0.3)  quasars. This comparison lets us infer the geometrical factor f needed to deproject the line-of-sight velocity component of the emitting gas. The wide range of f values we found, together with the strong dependence of f on the observed linewidth, suggests that a disc-like model for the broad-line region is preferable to an isotropic model, both for radio-loud and radio-quiet quasars. We examined similar observations of the C  iv line and found no correlation in the width of the two lines. Our results indicate that an inflated disc broad-line region, in which the Carbon line is emitted in a flat disc while Hβ is produced in a geometrically thick region, can account for the observed differences in the width and shape of the two emission lines.     

The effect of stellar feedback on the formation and evolution of gas and dust tori in AGN

Recently, the existence of geometrically thick dust structures in active galactic nuclei (AGN) has been directly proven with the help of interferometric methods in the mid-infrared. The observations are consistent with a two-component model made up of a geometrically thin and warm central disc, surrounded by a colder, fluffy torus component. Within the framework of an exploratory study, we investigate one possible physical mechanism, which could produce such a structure, namely the effect of stellar feedback from a young nuclear star cluster on the interstellar medium in centres of AGN. The model is realized by numerical simulations with the help of the hydrodynamics code tramp . We follow the evolution of the interstellar medium by taking discrete mass-loss and energy ejection due to stellar processes, as well as optically thin radiative cooling into account. In a post-processing step, we calculate observable quantities like spectral energy distributions (SEDs) and surface brightness distributions with the help of the radiative transfer code mc3d . The interplay between injection of mass, supernova explosions and radiative cooling leads to a two-component structure made up of a cold geometrically thin, but optically thick and very turbulent disc residing in the vicinity of the angular momentum barrier, surrounded by a filamentary structure. The latter consists of cold long radial filaments flowing towards the disc and a hot tenuous medium in between, which shows both inwards and outwards directed motions. With the help of this modelling, we are able to reproduce the range of observed neutral hydrogen column densities of a sample of Seyfert galaxies as well as the relation between them and the strength of the silicate 10 μm spectral feature. Despite being quite crude, our mean Seyfert galaxy model is even able to describe the SEDs of two intermediate type Seyfert galaxies observed with the Spitzer Space Telescope .     

Disc-like objects in hierarchical hydrodynamical simulations: comparison with observations

We present results from a careful and detailed analysis of the structural and dynamical properties of a sample of 29 disc-like objects identified at z =0 in three AP3M–SPH fully consistent cosmological simulations. These simulations are realizations of a CDM hierarchical model, in which an inefficient Schmidt-law-like algorithm to model the stellar formation process has been implemented. We focus on properties that can be constrained with available data from observations of spiral galaxies, namely the bulge and disc structural parameters and the rotation curves. Comparison with data from Broeils, de Jong and Courteau gives satisfactory agreement, in contrast with previous findings using other codes. This suggests that the stellar formation implementation we have used has succeeded in forming compact bulges that stabilize disc-like structures in the violent phases of their assembly, while in the quiescent phases the gas has cooled and collapsed in accord with the Fall & Efstathiou standard model of disc formation.     

Black hole spin and radio loudness in a Λ cold dark matter universe

We use a combination of a cosmological N -body simulation of the concordance Λ cold dark matter paradigm and a semi-analytic model of galaxy formation to investigate the spin development of central supermassive black holes (BHs) and its relation to the BH host galaxy properties. In order to compute BH spins, we use the α model of Shakura & Sunyaev and consider the King et al. warped disc alignment criterion. The orientation of the accretion disc is inferred from the angular momentum of the source of accreted material, which bears a close relationship to the large-scale structure in the simulation. We find that the final BH spin depends almost exclusively on the accretion history and only weakly on the warped disc alignment. The main mechanisms of BH spin-up are found to be gas cooling processes and disc instabilities, a result that is only partially compatible with Monte Carlo models where the main spin-up mechanisms are major mergers and disc instabilities; the latter results are reproduced when implementing randomly oriented accretion discs in our model. Regarding the BH population, we find that more massive BHs, which are hosted by massive ellipticals, have higher spin values than less massive BHs, hosted by spiral galaxies. We analyse whether gas accretion rates and BH spins can be used as tracers of the radio loudness of active galactic nuclei (AGN). We find that the current observational indications of an increasing trend of radio-loud AGN fractions with stellar and BH mass can be easily obtained when placing lower limits on the BH spin, with a minimum influence from limits on the accretion rates; a model with random accretion disc orientations is unable to reproduce this trend. Our results favour a scenario where the BH spin is a key parameter to separate the radio-loud and radio-quiet galaxy populations.