Nuclear starburst discs (NSDs) are very compact star-forming regions in the centers of galaxies that have been studied as a possible origin for the absorbing gas around a central active galactic nucleus. NSDs may be most relevant at $z\sim 1$ when obscured accretion onto supermassive black holes (SMBHs) is common. This paper describes the characteristics of the stellar remnants of NSDs at $z=0.01$, taking into account the evolution from $z=1$. Using a stellar synthesis model, the colours, masses, and luminosities of the stellar remnants are computed for a suite of 192 two-dimensional NSD models. These properties are compared to observations of local nuclear star clusters (NSCs), and a good match is found between the predicted and observed properties. Dynamical effects will likely cause the final remnant to be a rotating, nearly spherical distribution. In addition, $\approx 20$% of the NSD remnants have half-light radii <~ 10 pc, consistent with NSCs hosted in both late-type and early-type galaxies, and all the remnants follow similar size-luminosity relationships as observed in nearby NSCs. NSDs require the presence of a central SMBH and the most massive and compact stellar remnants are associated with the most massive SMBHs, although stellar clusters with a variety of sizes can be produced by all considered SMBH masses. Overall, NSDs at $z\sim 1$ appear to be a promising origin for the $\gg 1$ Gyr NSC population in early- and late-type galaxies with large SMBHs.

A significant number of active galactic nuclei (AGNs) are observed to be hidden behind dust and gas. The distribution of material around AGNs plays an important role in modeling the cosmic X-ray background (CXB), especially the fraction of Type-2 AGNs ($f_2$). One of the possible explanations for the obscuration in Seyfert galaxies at intermediate redshift is dusty starburst discs. We compute the 2D hydrostatic structure of 768 nuclear starburst discs (NSDs) under various physical conditions and also the distribution of column density along the line of sight ($N_{\text{H}}$) associated with these discs. Then, the $N_{\text{H}}$ distribution is evolved with redshift by using the redshift dependent distribution function of input parameters. The $f_2$ shows a strong positive evolution up to $z=2$, but only a weak level of enhancement at higher $z$. The Compton-thin and Compton-thick AGN fractions associated with these starburst regions increase as $\propto (1+z)^{\delta}$ where the $\delta$ is estimated to be 1.12 and 1.45, respectively. The reflection parameter $R_f$ associated with column density $N_{\text{H}} \geq 10^{23.5}$ cm$^{-2}$ extends from 0.13 at $z=0$ to 0.58 at $z=4$. A CXB model employing this evolving $N_{\text{H}}$ distribution indicates more compact ($R_{\text{out}}<120$ pc) NSDs provide a better fit to the CXB. In addition to "Seyfert-like" AGNs obscured by nuclear starbursts, we predict that 40 to 60 per cent of quasars must be Compton-thick to produce the peak of the CXB spectrum within observational uncertainty. The predicted total number counts of AGNs in 8-24 keV band are in fair agreement with observations from NuSTAR.

Nuclear starburst discs (NSDs) are star-forming discs that may be residing in the nuclear regions of active galaxies at intermediate redshifts. One dimensional (1D) analytical models developed by Thompson et al. (2005) show that these discs can possess an inflationary atmosphere when dust is sublimated on parsec scales. This make NSDs a viable source for AGN obscuration. We model the two dimensional (2D) structure of NSDs using an iterative method in order to compute the explicit vertical solutions for a given annulus. These solutions satisfy energy and hydrostatic balance, as well as the radiative transfer equation. In comparison to the 1D model, the 2D calculation predicts a less extensive expansion of the atmosphere by orders of magnitude at the parsec/sub-parsec scale, but the new scale-height $h$ may still exceed the radial distance $R$ for various physical conditions. A total of 192 NSD models are computed across the input parameter space in order to predict distributions of a line of sight column density $N_H$. Assuming a random distribution of input parameters, the statistics yield 56% of Type 1, 23% of Compton-thin Type 2s (CN), and 21% of Compton-thick (CK) AGNs. Depending on a viewing angle ($\theta$) of a particular NSD (fixed physical conditions), any central AGN can appear to be Type 1, CN, or CK which is consistent with the basic unification theory of AGNs. Our results show that $\log[N_H(\text{cm}^{-2})]\in$ [23,25.5] can be oriented at any $\theta$ from 0$^\circ$ to $\approx$80$^\circ$ due to the degeneracy in the input parameters.

Obscured AGNs provide an opportunity to study the material surrounding the central engine. Geometric and physical constraints on the absorber can be deduced from the reprocessed AGN emission. In particular, the obscuring gas may reprocess the nuclear X-ray emission producing a narrow Fe K$\alpha$ line and a Compton reflection hump. In recent years, models of the X-ray reflection from an obscuring torus have been computed; however, although the reflecting gas may be dusty, the models do not yet take into account the effects of dust on the predicted spectrum. We study this problem by analyzing two sets of models, with and without the presence of dust, using the one dimensional photo-ionization code Cloudy. The calculations are performed for a range of column densities ($22 <{\rm log}[N_H(\rm cm^{-2})]< 24.5$ ) and hydrogen densities ( $6 <{\rm log}[n_H(\rm cm^{-3})]< 8$). The calculations show the presence of dust can enhance the Fe K$\alpha$ equivalent width (EW) in the reflected spectrum by factors up to $\approx$ 8 for Compton thick (CT) gas and a typical ISM grain size distribution. The enhancement in EW with respect to the reflection continuum is due to the reduction in the reflected continuum intensity caused by the anisotropic scattering behaviour of dust grains. This effect will be most relevant for reflection from distant, predominately neutral gas, and is a possible explanation for AGNs which show a strong Fe K$\alpha$ EW and a relatively weak reflection continuum. Our results show it is an important to take into account dust while modeling the X-ray reflection spectrum, and that inferring a CT column density from an observed Fe K$\alpha$ EW may not always be valid. Multi-dimensional models are needed to fully explore the magnitude of the effect.