We report on high-resolution simulations that explore the orbital decay of massive black hole (MBH) pairs with masses between $10^5$ and $10^7 M_{\odot}$ embedded in a circumnuclear gas disk (CND). An adiabatic equation of state is adopted, with a range of adiabatic indices, which maintains a smooth flow. Mergers between MBHs in this mass range would be detectable by the upcoming Laser Inteferometer Space Antenna (LISA). The orbital evolution is followed from the CND scale ($100$~pc) down to separations of $0.1$--$0.01$~pc at which a circumbinary disk (CBD) could form. The decay is erratic and strongly dependent on the gas flow within the disk, that ultimately determines the net torques experienced by the sinking MBH. Overall, we can identify three different evolutionary stages: (i) an initially slow decay that leads to no significant change in the orbital angular momentum, resulting in some circularization; (ii) a fast migration phase in which the orbital angular momentum decreases rapidly; and (iii) a final, very slow decay phase, in which orbital angular momentum can even increase, and a CBD can form. The fast migration phase owes to disk-driven torques originating primarily from the co-orbital region of the secondary MBH, at a distance of 1--3 Hill radii. We find strong analogies with fast Type III migration for massive planets in protoplanetary disks. The CBD forms only when the decay rate becomes small enough to allow it enough time to carve a cavity around the primary MBH, at scales $\lesssim 1$~pc; when this happens, the MBH separation nearly stalls in our higher-resolution run. We suggest an empirically modified gap-opening criterion that takes into account such timescale effects as well as other deviations from standard assumptions made in the literature. [Abriged]

A Friedmann like cosmological model in Einstein-Cartan framework is studied when the torsion function is assumed to be proportional to a single $\phi(t)$ function coming just from the spin vector contribution of ordinary matter. By analysing four different types of torsion function written in terms of one, two and three free parameters, we found that a model with $\phi(t)=- \alpha H(t) \big({\rho_{m}(t)}/{\rho_{0c}}\big)^n$ is totally compatible with recent cosmological data, where $\alpha$ and $n$ are free parameters to be constrained from observations, $\rho_m$ is the matter energy density and $\rho_{0c}$ the critical density. The recent accelerated phase of expansion of the universe is correctly reproduced by the contribution coming from torsion function, with a deceleration parameter indicating a transition redshift of about $0.65$.

Light pollution poses a growing threat to optical astronomy, in addition to its detrimental impacts on the natural environment, the intangible heritage of humankind related to the contemplation of the starry sky and, potentially, on human health. The computation of maps showing the spatial distribution of several light pollution related functions (e.g. the anthropogenic zenithal night sky brightness, or the average brightness of the celestial hemisphere) is a key tool for light pollution monitoring and control, providing the scientific rationale for the adoption of informed decisions on public lighting and astronomical site preservation. The calculation of such maps from satellite radiance data for wide regions of the planet with sub-kilometric spatial resolution often implies a huge amount of basic pixel operations, requiring in many cases extremely large computation times. In this paper we show that, using adequate geographical projections, a wide set of light pollution map calculations can be reframed in terms of two-dimensional convolutions that can be easily evaluated using conventional fast Fourier-transform (FFT) algorithms, with typical computation times smaller than 10^-6 s per output pixel.

The light pollution levels experienced at any given site generally depend on a wide number of artificial light sources distributed throughout the surrounding territory. Since photons can travel long distances before being scattered by the atmosphere, any effective proposal for reducing local light pollution levels needs an accurate assessment of the relative weight of all intervening light sources, including those located tens or even hundreds of km away. In this paper we describe several ways of quantifying and visualizing these relative weights. Particular emphasis is made on the aggregate contribution of the municipalities, which are -- in many regions of the world -- the administrative bodies primarily responsible for the planning and maintenance of public outdoor lighting systems.

Several light pollution indicators are commonly used to monitor the effects of the transition from outdoor lighting systems based on traditional gas-discharge lamps to solid-state light sources. In this work we analyze a subset of these indicators, including the artificial zenithal night sky brightness in the visual photopic and scotopic bands, the brightness in the specific photometric band of the widely used Sky Quality Meter (SQM), and the top-of-atmosphere radiance detected by the VIIRS-DNB radiometer onboard the satellite Suomi-NPP. Using a single-scattering approximation in a layered atmosphere we quantitatively show that, depending on the transition scenarios, these indicators may show different, even opposite behaviors. This is mainly due to the combined effects of the changes in the sources' spectra and angular radiation patterns, the wavelength-dependent atmospheric propagation processes and the differences in the detector spectral sensitivity bands. It is suggested that the possible presence of this differential behavior should be taken into account when evaluating light pollution indicator datasets for assessing the outcomes of public policy decisions regarding the upgrading of outdoor lighting systems.

Under stable atmospheric conditions, the zenithal brightness of the urban sky varies throughout the night following the time course of the anthropogenic emissions of light. Different types of artificial light sources (e.g. streetlights, residential, and vehicle lights) present specific time signatures, and this feature makes it possible to estimate the amount of sky brightness contributed by each one of them. Our approach is based on transforming the time representation of the zenithal sky brightness into a modal coefficients one, in terms of the time course signatures of the sources. The modal coefficients, and hence the absolute and relative contributions of each type of source, can be estimated from the measured brightness by means of linear least squares fits. A method for determining the time signatures is described, based on wide-field time-lapse photometry of the urban nightscape. Our preliminary results suggest that artificial light leaking out of the windows of residential buildings may account for a significant share of the time-varying part of the zenithal sky brightness, whilst the contribution of the vehicle lights seems to be significantly smaller.

We study the orbital decay of a pair of massive black holes (BHs) with masses 5 * 10^5 and 10^7 M_sun, using hydrodynamical simulations of circumnuclear disks (CNDs) with the alternating presence of sub-grid physics such as radiative cooling, star formation, supernova feedback, BH accretion and feedback. In the absence of such processes, the orbit of the secondary BH decays over timescales of ~10 Myr to the center of the CND, where the primary BH resides. When strong dissipation operates in CNDs, fragmentation into massive objects the size of giant molecular clouds and with densities in the range 10^4 - 10^7 amu / cm^3 occurs, causing stochastic torques and hits that can eject the secondary BH from the midplane. Outside the plane, the low-density medium provides only weak drag, and the BH return is governed by inefficient dynamical friction. In rare cases, clump-BH interactions can lead to a faster decay. Feedback processes lead to outflows, but do not change significantly the overall density of the CND midplane. However, with a spherically distributed BH feedback a hot bubble is generated behind the secondary, which almost shuts off dynamical friction, a phenomenon we dub "wake evacuation", leading to delays in the decay of possibly ~0.3 Gyr. We discuss the non-trivial implications on the discovery space of the eLISA telescope. Our results suggest the largest uncertainty in predicting BH merger rates lies in the potentially wide variety of galaxy host systems, with different degrees of gas dissipation and heating, yielding decay timescales from ~10 to ~300 Myr.