CIGALE is a powerful multiwavelength spectral energy distribution (SED) fitting code for extragalactic studies. However, the current version of CIGALE is not able to fit X-ray data, which often provide unique insights into AGN intrinsic power. We develop a new X-ray module for CIGALE, allowing it to fit SEDs from the X-ray to infrared (IR). We also improve the AGN fitting of CIGALE from UV-to-IR wavelengths. We implement a modern clumpy two-phase torus model, SKIRTOR. To account for moderately extincted type 1 AGNs, we implement polar-dust extinction. We publicly release the source code (named X-CIGALE). We test X-CIGALE with X-ray detected AGNs in SDSS, COSMOS, and AKARI-NEP. The fitting quality (as indicated by reduced $\chi^2$) is good in general, indicating that X-CIGALE is capable of modelling the observed SED from X-ray to IR. We discuss constrainability and degeneracy of model parameters in the fitting of AKARI-NEP, for which excellent mid-IR photometric coverage is available. We also test fitting a sample of AKARI-NEP galaxies for which only X-ray upper limits are available from Chandra observations, and find that the upper limit can effectively constrain the AGN SED contribution for some systems. Finally, using X-CIGALE, we assess the ability of Athena to constrain the AGN activity in future extragalactic studies.

There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z = 1-4 $ galaxies using the SKIRT simulated spectral energy distributions of galaxies sampled from the EAGLE simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR$-{\rm M}_{\star}$ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm \mu m}$, MIPS - 24, 70 and 160 ${\rm \mu m}$ or SPIRE - 250, 350, 500 ${\rm \mu m}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by -0.5 dex at $z \sim 4$ but overpredict the measurements by 0.3 dex at $z \sim 1$. Relations relying on dust-corrected rest-frame UV luminosities, are flatter since they overpredict/underpredict SFRs for low/high star forming objects and yield deviations from the intrinsic relation from 0.10 dex to -0.13 dex at $z \sim 4$. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs.

We present a detailed study of ionized outflows in a large sample of ~650 hard X-ray detected AGN. Using optical spectroscopy from the BAT AGN Spectroscopic Survey (BASS) we are able to reveal the faint wings of the [OIII] emission lines associated with outflows covering, for the first time, an unexplored range of low AGN bolometric luminosity at low redshift (z~0.05). We test if and how the incidence and velocity of ionized outflow is related to AGN physical parameters: black hole mass, gas column density, Eddington Ratio, [OIII], X-ray, and bolometric luminosities. We find a higher occurrence of ionized outflows in type 1.9 (55%) and type 1 AGN (46%) with respect to type 2 AGN (24%). While outflows in type 2 AGN are evenly balanced between blue and red velocity offsets with respect to the [OIII] narrow component, they are almost exclusively blueshifted in type 1 and type 1.9 AGN. We observe a significant dependence between the outflow occurrence and accretion rate, which becomes relevant at high Eddington ratios (> -1.7). We interpret such behaviour in the framework of covering factor-Eddington ratio dependence. We don't find strong trends of the outflow maximum velocity with AGN physical parameters, as an increase with bolometric luminosity can be only identified when including samples of AGN at high luminosity and high redshift taken from literature.

The scatter (${\rm\sigma_{\text{sSFR}}}$) of the specific star formation rates (sSFRs) of galaxies is a measure of the diversity in their star formation histories (SFHs) at a given mass. In this paper we employ the EAGLE simulations to study the dependence of the ${\rm \sigma_{\text{sSFR}}}$ of galaxies on stellar mass (${\rm M_{\star}}$) through the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation in $ {\rm z \sim 0-4}$. We find that the relation evolves with time, with the dispersion depending on both stellar mass and redshift. The models point to an evolving U-shape form for the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation with the scatter being minimal at a characteristic mass $M^{\star}$ of ${\rm 10^{9.5}}$ ${\rm M_{\odot}}$ and increasing both at lower and higher masses. This implication is that the diversity of SFHs increases towards both at the low- and high-mass ends. We find that active galactic nuclei feedback is important for increasing the ${\rm \sigma_{\text{sSFR}}}$ for high mass objects. On the other hand, we suggest that SNe feedback increases the ${\rm \sigma_{\text{sSFR}}}$ of galaxies at the low-mass end. We also find that excluding galaxies that have experienced recent mergers does not significantly affect the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation. Furthermore, we employ the combination of the EAGLE simulations with the radiative transfer code SKIRT to evaluate the effect of SFR/stellar mass diagnostics in the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation and find that the ${\rm SFR/M_{\star}}$ methodologies (e.g. SED fitting, UV+IR, UV+IRX-$\beta$) widely used in the literature to obtain intrinsic properties of galaxies have a large effect on the derived shape and normalization of the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Recent observations which resolved the mid-infrared (MIR) emission of nearby active galactic nuclei (AGN), surprisingly revealed that their dust emission appears prominently extended in the polar direction, at odds with the expectations from the canonical dusty torus. This polar dust, tentatively associated with dusty winds driven by radiation pressure, is found to have a major contribution to the MIR flux from scales of a few to hundreds of parsecs. When facing a potential change of paradigm, case studies of objects with the best intrinsic resolution are essential. One such source with a clear detection of polar dust is a nearby, well-known AGN in the Circinus galaxy. In the first paper, we successfully explained the peculiar MIR morphology of Circinus observed on large, tens of parsec scales with a model consisting of a compact dusty disc and an extended hollow dusty cone. In this work, we further refine the model on smaller, parsecs scales to test whether it can also explain the MIR interferometric data. We find that a model composed of a thin dusty disc seen almost edge-on and a polar outflow in the form of a hyperboloid shell can reproduce well the VLTI/MIDI observations at all wavelengths, baselines and position angles. In contrast, while providing a good fit to the integrated MIR spectrum, the dusty torus model fails to reproduce the spatially resolved interferometric data. We put forth the disc$+$hyperboloid wind model of Circinus AGN as a prototype for the dust structure in the AGN population with polar dust.

We quantify the luminosity contribution of active galactic nuclei (AGN) to the 12 $\mu$m, mid-infrared (MIR; 5-38 $\mu$m), and the total IR (5-1000 $\mu$m) emission in the local AGN detected in the all-sky 70-month Swift/Burst Alert Telescope (BAT) ultra hard X-ray survey. We decompose the IR spectral energy distributions (SEDs) of 587 objects into AGN and starburst components using AGN torus and star-forming galaxy templates. This enables us to recover the AGN torus emission also for low-luminosity end, down to $\log (L_{14-150}/{\rm erg}~{\rm s}^{-1}) \simeq 41$, which typically have significant host galaxy contamination. We find that the luminosity contribution of the AGN to the 12 $\mu$m, the MIR, and the total IR band is an increasing function of the 14-150 keV luminosity. We also find that for the most extreme cases, the IR pure-AGN emission from the torus can extend up to 90 $\mu$m. The obtained total IR AGN luminosity through the IR SED decomposition enables us to estimate the fraction of the sky obscured by dust, i.e., the dust covering factor. We demonstrate that the median of the dust covering factor is always smaller than that of the X-ray obscuration fraction above the AGN bolometric luminosity of $\log (L_{\rm bol}/{\rm erg}~{\rm s}^{-1}) \simeq 42.5$. Considering that X-ray obscuration fraction is equivalent to the covering factor coming from both the dust and gas, it indicates that an additional neutral gas component, along with the dusty torus, is responsible for the absorption of X-ray emission.

Active galactic nuclei (AGN) are one of the most luminous objects in the Universe, emitting powerful continuum and line emission across all wavelength bands. They represent an important link in the investigations of the galaxy evolution and cosmology. The resolving of the AGN inner structure is still a difficult task with current instruments, therefore the spectroscopy and spectropolarimetry are crucial tools to investigate these objects and their components, such as the properties of the supermassive black hole, the broad line region, and the dusty torus. In this review, we present the results of the project "Astrophysical spectroscopy of extragalactic objects", from the observations, data processing and analysis, to the modelling of different regions in AGN.

Recent high angular resolution observations resolved for the first time the mid-infrared (MIR) structure of nearby active galactic nuclei (AGN). Surprisingly, they revealed that a major fraction of their MIR emission comes from the polar regions. This is at odds with the expectation based on AGN unification, which postulates a dusty torus in the equatorial region. The nearby, archetypical AGN in the Circinus galaxy offers one of the best opportunities to study the MIR emission in greater detail. New, high quality MIR images obtained with the upgraded VISIR instrument at the Very Large Telescope show that the previously detected bar-like structure extends up to at least 40 pc on both sides of the nucleus along the edges of the ionization cone. Motivated by observations across a wide wavelength range and on different spatial scales, we propose a phenomenological dust emission model for the AGN in the Circinus galaxy consisting of a compact dusty disk and a large-scale dusty cone shell, illuminated by a tilted accretion disk with an anisotropic emission pattern. Undertaking detailed radiative transfer simulations, we demonstrate that such a model is able to explain the peculiar MIR morphology and account for the entire IR spectral energy distribution. Our results call for caution when attributing dust emission of unresolved sources entirely to the torus and warrant further investigation of the MIR emission in the polar regions of AGN.

Polarization is an important tool to further the understanding of interstellar dust and the sources behind it. In this paper we describe our implementation of polarization that is due to scattering of light by spherical grains and electrons in the dust Monte Carlo radiative transfer code SKIRT. In contrast to the implementations of other Monte Carlo radiative transfer codes, ours uses co-moving reference frames that rely solely on the scattering processes. It fully supports the peel-off mechanism that is crucial for the efficient calculation of images in 3D Monte Carlo codes. We develop reproducible test cases that push the limits of our code. The results of our program are validated by comparison with analytically calculated solutions. Additionally, we compare results of our code to previously published results. We apply our method to models of dusty spiral galaxies at near-infrared and optical wavelengths. We calculate polarization degree maps and show them to contain signatures that trace characteristics of the dust arms independent of the inclination or rotation of the galaxy.

The primary source of emission of active galactic nuclei (AGN), the accretion disk, is surrounded by an optically and geometrically thick dusty structure ("the so-called dusty torus"). The infrared radiation emitted by the dust is nothing but a reprocessed fraction of the accretion disk emission, so the ratio of the torus to the AGN luminosity ($L_{\text{torus}}/L_{\text{AGN}}$) should correspond to the fraction of the sky obscured by dust, i.e. the covering factor. We undertook a critical investigation of the $L_{\text{torus}}/L_{\text{AGN}}$ as the dust covering factor proxy. Using state-of-the-art 3D Monte Carlo radiative transfer code, we calculated a grid of spectral energy distributions (SEDs) emitted by the clumpy two-phase dusty structure. With this grid of SEDs, we studied the relation between $L_{\text{torus}}/L_{\text{AGN}}$ and the dust covering factor for different parameters of the torus. We found that in the case of type 1 AGNs the torus anisotropy makes $L_{\text{torus}}/L_{\text{AGN}}$ underestimate low covering factors and overestimate high covering factors. In type 2 AGNs $L_{\text{torus}}/L_{\text{AGN}}$ always underestimates covering factors. Our results provide a novel easy-to-use method to account for anisotropy and obtain correct covering factors. Using two samples from the literature, we demonstrated the importance of our result for inferring the obscured AGN fraction. We found that after the anisotropy is properly accounted for, the dust covering factors show very weak dependence on $L_{\text{AGN}}$, with values in the range of $\approx0.6-0.7$. Our results also suggest a higher fraction of obscured AGNs at high luminosities than those found by X-ray surveys, in part owing to the presence of a Compton-thick AGN population predicted by population synthesis models.

Hot, Dust-Obscured Galaxies (Hot DOGs), selected from the WISE all sky infrared survey, host some of the most powerful Active Galactic Nuclei (AGN) known, and might represent an important stage in the evolution of galaxies. Most known Hot DOGs are at $z> 1.5$, due in part to a strong bias against identifying them at lower redshift related to the selection criteria. We present a new selection method that identifies 153 Hot DOG candidates at $z\sim 1$, where they are significantly brighter and easier to study. We validate this approach by measuring a redshift $z=1.009$, and an SED similar to higher redshift Hot DOGs for one of these objects, WISE J1036+0449 ($L_{\rm\,Bol}\simeq 8\times 10^{46}\rm\,erg\,s^{-1}$), using data from Keck/LRIS and NIRSPEC, SDSS, and CSO. We find evidence of a broadened component in MgII, which, if due to the gravitational potential of the supermassive black hole, would imply a black hole mass of $M_{\rm\,BH}\simeq 2 \times 10^8 M_{\odot}$, and an Eddington ratio of $\lambda_{\rm\,Edd}\simeq 2.7$. WISE J1036+0449 is the first Hot DOG detected by NuSTAR, and the observations show that the source is heavily obscured, with a column density of $N_{\rm\,H}\simeq(2-15)\times10^{23}\rm\,cm^{-2}$. The source has an intrinsic 2-10 keV luminosity of $\sim 6\times 10^{44}\rm\,erg\,s^{-1}$, a value significantly lower than that expected from the mid-infrared/X-ray correlation. We also find that the other Hot DOGs observed by X-ray facilities show a similar deficiency of X-ray flux. We discuss the origin of the X-ray weakness and the absorption properties of Hot DOGs. Hot DOGs at $z\lesssim1$ could be excellent laboratories to probe the characteristics of the accretion flow and of the X-ray emitting plasma at extreme values of the Eddington ratio.

We present an analysis of 43 years (1972 to 2015) of spectroscopic observations of the Seyfert 1 galaxy NGC 5548. This includes 12 years of new unpublished observations (2003 to 2015). We compiled about 1600 H$\beta$ spectra and analyzed the long-term spectral variations of the 5100 \AA\ continuum and the H$\beta$ line. Our analysis is based on standard procedures including the Lomb-Scargle method, which is known to be rather limited to such heterogeneous data sets, and new method developed specifically for this project that is more robust and reveals a $\sim$5700 day periodicity in the continuum light curve, the H$\beta$ light curve and the radial velocity curve of the red wing of the H$\beta$ line. The data are consistent with orbital motion inside the broad emission line region of the source. We discuss several possible mechanisms that can explain this periodicity, including orbiting dusty and dust-free clouds, a binary black hole system, tidal disruption events, and the effect of an orbiting star periodically passing through an accretion disc.

If the existence of an obscuring circumnuclear region around the innermost regions of active galactic nuclei (AGN) has been observationally proven, its geometry remains highly uncertain. The morphology usually adopted for this region is a toroidal structure, but other alternatives, such as flared disks, can be a good representative of equatorial outflows. Those two geometries usually provide very similar spectroscopic signatures, even when they are modeled under the assumption of fragmentation. In this lecture note, we show that the resulting polarization signatures of the two models, either a torus or a flared disk, are quite different from each other. We use a radiative transfer code that computes the 2000 - 8000 angstrom polarization of the two morphologies in a clumpy environment, and show that varying the sizes of a toroidal region has deep impacts onto the resulting polarization, while the polarization of flared disks is independent of the outer radius. Clumpy flared disks also produce higher polarization degrees (about 10 % at best) together with highly variable polarization position angles.

The narrow component of the iron K$\alpha$ is an almost ubiquitous feature in the X-ray spectra of active galactic nuclei (AGN) and is believed to originate in neutral material, possibly located in the molecular torus. This would imply a tight connection between the Fe K$\alpha$ equivalent width (EW) and the physical properties of the torus. In a recent work we have shown that the decrease of the covering factor of the torus with the luminosity, as expected by luminosity-dependent unification models, would be able to explain the decrease of Fe K$\alpha$ EW with the luminosity (i.e., the X-ray Baldwin effect). Recent developments in the study of the mid-IR (MIR) spectrum of AGN allow important parameters of the torus to be deduced, such as its covering factor ($f_{\rm\,obs}$) and equatorial column density ($N_{\rm\,H}^{\rm\,T}$), by applying clumpy torus models. Using XMM-Newton/EPIC observations of a sample of 24 type-I AGN, we investigate the relation between the physical parameters of the torus obtained by recent MIR works and the properties of the Fe K$\alpha$ line. We correct the values of the Fe K$\alpha$ EW by taking the inclination angle, the photon index, the equatorial column density, and half-opening angle of the torus into account using a physical torus model of X-ray reprocessed radiation. We find that the relation between Fe K$\alpha$ EW and $f_{\rm\,obs}$ shows a slope that is consistent with the expected value, albeit with a low statistical significance. A trend that is consistent with the theoretical prediction is also found when comparing the Fe K$\alpha$ EW to $N_{\rm\,H}^{\rm\,T}$. Our work seems to confirm that the bulk of the narrow Fe K$\alpha$ line is produced by the same material responsible for the MIR emission.

A crucial aspect of 3D Monte Carlo radiative transfer is the choice of the spatial grid used to partition the dusty medium. We critically investigate the use of octree grids in Monte Carlo dust radiative transfer, with two different octree construction algorithms (regular and barycentric subdivision) and three different octree traversal algorithms (top-down, neighbour list, and the bookkeeping method). In general, regular octree grids need higher levels of subdivision compared to the barycentric grids for a fixed maximum cell mass threshold criterion. The total number of grid cells, however, depends on the geometry of the model. Surprisingly, regular octree grid simulations turn out to be 10 to 20% more efficient in run time than the barycentric grid simulations, even for those cases where the latter contain fewer grid cells than the former. Furthermore, we find that storing neighbour lists for each cell in an octree, ordered according to decreasing overlap area, is worth the additional memory and implementation overhead: using neighbour lists can cut down the grid traversal by 20% compared to the traditional top-down method. In conclusion, the combination of a regular node subdivision and the neighbour list method results in the most efficient octree structure for Monte Carlo radiative transfer simulations.

We modeled the AGN dusty torus as a clumpy two-phase medium, with high-density clumps embedded in a low-density interclump dust. To obtain spectral energy distributions and images of the torus at different wavelengths, we employed the 3D Monte Carlo radiative transfer code SKIRT. Apart from the grid of two-phase models, we calculated the corresponding sets of clumps-only models and models with a smooth dust distribution for comparison. We found that the most striking feature of the two-phase model is that it might offer a natural solution to the common issue reported in a number of papers -- the observed excess of the near-infrared emission.

We investigated the gravitational microlensing of active galactic nucleus dusty tori in the case of lensed quasars in the infrared domain. The dusty torus is modeled as a clumpy two-phase medium. To obtain spectral energy distributions and images of tori at different wavelengths, we used the 3D Monte Carlo radiative transfer code SKIRT. A ray-shooting technique has been used to calculate microlensing magnification maps. We simulated microlensing by the stars in the lens galaxy for different configurations of the lensed system and different values of the torus parameters, in order to estimate (a) amplitudes and timescales of high magnification events, and (b) the influence of geometrical and physical properties of dusty tori on light curves in the infrared domain. We found that, despite their large size, dusty tori could be significantly affected by microlensing in some cases, especially in the near-infrared domain (rest-frame). The very long time-scales of such events, in the range from several decades to hundreds of years, are limiting the practical use of this method to study the properties of dusty tori. However, our results indicate that, when studying flux ratios between the images in different wavebands of lensed quasars, one should not disregard the possibility that the near and mid-infrared flux ratios could be under the influence of microlensing.

We investigated the emission of active galactic nuclei dusty tori in the infrared domain, with a focus on the 10 micron silicate feature. We modeled the dusty torus as a clumpy two-phase medium with high-density clumps and a low-density medium filling the space between the clumps. We employed a three-dimensional radiative transfer code to obtain spectral energy distributions and images of tori at different wavelengths. We calculated a grid of models for different parameters and analyzed the influence of these parameters on the shape of the mid-infrared emission. A corresponding set of clumps-only models and models with a smooth dust distribution is calculated for comparison. We found that the dust distribution, the optical depth and a random arrangement of clumps in the innermost region, all have an impact on the shape and strength of the silicate feature. The 10 micron silicate feature can be suppressed for some parameters, but models with smooth dust distribution are also able to produce a wide range of the silicate feature strength.

We study the photocenter position variability due to variations in the quasar inner structure. We consider variability in the accretion disk emissivity and torus structure variability due to different illumination by the central source. We discuss possible detection of these effects by Gaia. Observations of the photocenter variability in two AGNs, SDSS J121855+020002 and SDSS J162011+1724327 have been reported and discussed. With investigation of the variations in the quasar inner structure we explore how much this effect can affect the position determination and whether it can be (or not) detected with Gaia mission. We used (a) a model of a relativistic disk, including the perturbation that can increase brightness of a part of the disk, and consequently offset the photocenter position, and (b) a model of a dusty torus which absorbs and re-emits the incoming radiation from accretion disk. We estimated the value of the photocenter offset due to these two effects. We found that perturbations in the inner structure can significantly offset the photocenter. It depends on the characteristics of perturbation and accretion disk and structure of the torus. In the case of two considered QSOs the observed photocenter offsets cannot be explained by variations in the accretion disk and other effects should be considered. We discussed possibility of exploding stars very close to the AGN source, and also possibility that there are two variable sources in the center of these two AGNs that may indicate a binary super-massive black hole system on a kpc (pc) scale. The Gaia mission seems to be very perspective, not only for astrometry, but also for exploring the inner structure of AGNs. We conclude that variations in the quasar inner structure can affect the observed photocenter (up to several mas). There is a chance to observe such effect in the case of bright and low-redshifted QSOs.

We investigate the emission of active galactic nuclei (AGN) dusty tori in the infrared domain. Following theoretical predictions coming from hydrodynamical simulations, we model the dusty torus as a 3D two-phase medium with high-density clumps and low-density medium filling the space between the clumps. Spectral energy distributions (SED) and images of the torus at different wavelengths are obtained using 3D Monte Carlo radiative transfer code SKIRT. Our approach of generating clumpy structure allows us to model tori with single clumps, complex structures of merged clumps or interconnected sponge-like structure. A corresponding set of clumps-only models and models with smooth dust distribution is calculated for comparison. We found that dust distribution, optical depth, clump size and their actual arrangement in the innermost region, all have an impact on the shape of near- and mid-infrared SED. The 10 micron silicate feature can be suppressed for some parameters, but models with smooth dust distribution are also able to produce a wide range of the silicate feature strength. Finally, we find that having the dust distributed in a two-phase medium, might offer a natural solution to the lack of emission in the near-infrared, compared to observed data, which affects clumpy models currently available in the literature.

We present an updated version of SKIRT, a 3D Monte Carlo radiative transfer code developed to simulate dusty galaxies. The main novel characteristics of the SKIRT code are the use of a stellar foam to generate random positions, an efficient combination of eternal forced scattering and continuous absorption, and a new library approach that links the radiative transfer code to the DustEM dust emission library. This approach enables a fast, accurate and self-consistent calculation of the dust emission of arbitrary mixtures of transiently heated dust grains and polycyclic aromatic hydrocarbons, even for full 3D models containing millions of dust cells. We have demonstrated the accuracy of the SKIRT code through a set of simulations based on the edge-on spiral galaxy UGC 4754. The models we ran were gradually refined from a smooth, 2D, LTE model to a fully 3D model that includes NLTE dust emission and a clumpy structure of the dusty ISM. We find that clumpy models absorb UV and optical radiation less efficiently than smooth models with the same amount of dust, and that the dust in clumpy models is on average both cooler and less luminous. Our simulations demonstrate that, given the appropriate use of optimization techniques, it is possible to efficiently and accurately run Monte Carlo radiative transfer simulations of arbitrary 3D structures of several million dust cells, including a full calculation of the NLTE emission by arbitrary dust mixtures.

Here we show that in the case when double peaked emission lines originate from outer parts of accretion disk, their variability could be caused by perturbations in the disk emissivity. In order to test this hypothesis, we introduced a model of disk perturbing region in the form of a single bright spot (or flare) by a modification of the power law disk emissivity in appropriate way. The disk emission was then analyzed using numerical simulations based on ray-tracing method in Kerr metric and the corresponding simulated line profiles were obtained. We applied this model to the observed H-beta line profiles of 3C 390.3 (observed in the period 1995-1999), and estimated the parameters of both, accretion disk and perturbing region. Our results show that two large amplitude outbursts of the H-beta line observed in 3C 390.3 could be explained by successive occurrences of two bright spots on approaching side of the disk. These bright spots are either moving, originating in the inner regions of the disk and spiralling outwards by crossing small distances during the period of several years, or stationary. In both cases, their widths increase with time, indicating that they most likely decay.

Interstellar dust grains efficiently absorb and scatter UV and optical radiation in galaxies, and therefore can significantly affect the apparent structure of spiral galaxies. We discuss the effect of dust attenuation on the observed structural properties of bulges and discs. We also present some first results on modelling the dust content of edge-on spiral galaxies using both optical and Herschel far-infrared data. Both of these results demonstrate the complex interplay of dust and star light in spiral galaxies.