Observational data suggest that a belt of planetesimals is expected close to the snow line in protoplanetary disks. Assuming there is such a belt in TRAPPIST-1 system, we examine possibilities of water delivery to the planets via planetesimals from the belt. The study is accomplished by numerical simulations of dynamical evolution of a hypothetical planetesimal belt. Our results show that the inner part of the belt is dynamically unstable and planetesimals located in this region are quickly scattered away, with many of them entering the region around the planets. The main dynamical mechanism responsible for the instability are close encounters with the outermost planet Trappist-1h. A low-order mean-motion resonance 2:3 with Trappist-1h, located in the same region, also contributes to the objects transport. In our nominal model, the planets have received non-negligible amount of water, with the smallest amount of $15$\% of the current Earth's water amount (EWA) being delivered to the planet 1b, while the planets Trappist-1e and Trappist-1g have received more than $60$\% of the EWA. We have found that while the estimated efficiency of water transport to the planets is robust, the amount of water delivered to each planet may vary significantly depending on the initial masses and orbits of the planets. The estimated dynamical "half-lives" have shown that the impactors' source region should be emptied in less then 1~Myr. Therefore, the obtained results suggest that transport of planetesimals through the system preferably occurs during an early phase of the planetary system evolution.

In this paper we present the updated empirical radio surface-brightness-to-diameter ($\Sigma$--$D$) relation for Galactic supernova remnants (SNRs) calibrated using $110$ SNRs with reliable distances. We apply orthogonal fitting procedure and kernel density smoothing in $\Sigma-D$ plane and compare the results with the latest theoretical $\Sigma-D$ relations derived from simulations of radio evolution of SNRs. We argue that the best agreement between the empirical and simulated $\Sigma-D$ relations is achieved if the mixed-morphology SNRs and SNRs of both, low brightness and small diameter, are filtered out from the calibration sample. The distances to $5$ newly discovered remnants and $27$ new candidates for shell SNRs are estimated from our full and filtered calibration samples.

The problem of the extent of habitable zones in different kinds of galaxies is one of the outstanding challenges for contemporary astrobiology. In the present study, we investigate habitability in a large sample of simulated galaxies from the Illustris Project in order to at least roughly quantify the hospitality to life of different galactic types. The pioneering study of Dayal et al. (2015) is critically examined and some of its results are amended. In particular, we find a tentative evidence for a second mode of galactic habitability comprising metal-rich dwarfs similar to IC 225, LMC or M32. The role of the galactic environment and the observation selection effects is briefly discussed and prospects for further research on the topic outlined.

Galaxies represent the main form of organization of matter in our universe. Therefore, they are of obvious interest for the new multidisciplinary field of astrobiology. In particular, to study habitability of galaxies represents one of the main emerging challenges of theoretical and numerical astrobiology. Its theoretical underpinnings are, however, often confused and vague. Here we present a systematic attempt to list and categorize major causal factors playing a role in emergent habitability of galaxies. Furthermore, we argue that the methodology of cosmological merger trees is particularly useful in delineating what are systematic and lawful astrobiological properties of galaxies at present epoch vs. those which are product of historical contingency and, in particular, interaction with wider extragalactic environment. Employing merger trees extracted from cosmological N-body simulations as a new and promising research method for astrobiology has been pioneered by Stanway et al. (2018). We analyse the general issue of applicability of merger trees and present preliminary results on a set of trees extracted from the Illustris Project. In a sense, this approach is directly complementary to using large-scale cosmological simulations to study habitable zones of individual galaxies with high mass/spatial resolution; taken together, they usher a new era of synergy and synthesis between cosmology and astrobiology.

In a recent paper in this journal, Lingam and Loeb (2018) develop an excellent heuristic for searches for biosignatures vs. technosignatures. We consider two ways in which their approach could be extended and sharpened, with focus on durability of technosignatures. We also note an important consequence of the adopted heuristic which offers strong support to the ideas of the Dysonian SETI.

Context. Astrobiological evolution of the Milky Way (or the shape of its "astrobiological landscape") has emerged as one of the key research topics in recent years. In order to build precise, quantitative models of the Galactic habitability, we need to account for two opposing tendencies of life and intelligence in the most general context: the tendency to spread to all available ecological niches (conventionally dubbed "colonization") and the tendency to succumb to various types of existential catastrophes ("catastrophism"). These evolutionary tendencies have become objects of study in fields such as ecology, macroevolution, risk analysis, and futures studies, while a serious astrobiological treatment has been lacking so far. Aims. Our aim is to numerically investigate the dynamics of opposed processes of expansion (panspermia, colonization) and extinction (catastrophic mechanisms) of life in the Galaxy. Methods. We employ a new type of numerical simulation based on 1D probabilistic cellular automaton with very high temporal resolution, in order to study astrobiological dynamics. Results. While the largest part of the examined parameter space shows very low habitability values, as expected, the remaining part has some observationally appealing features that imply, among other things, a reduction in the amount of fine-tuning necessary for resolving the Fermi paradox. Conclusions. Advanced aspects of Galactic habitability are amenable to precision studies using massive parallel computer simulations. There are regions of parameter space corresponding to a quasi-stationary state satisfying observable constraints and possessing viable SETI targets.

We construct the most complete sample of supernova remnants (SNRs) in any galaxy - the Large Magellanic Cloud (LMC) SNR sample. We study their various properties such as spectral index ($\alpha$), size and surface-brightness. We suggest an association between the spatial distribution, environment density of LMC SNRs and their tendency to be located around supergiant shells. We find evidence that the 16 known type Ia LMC SNRs are expanding in a lower density environment compared to the Core-Collapse (CC) type. The mean diameter of our entire population (74) is 41 pc, which is comparable to nearby galaxies. We didn't find any correlation between the type of SN explosion, ovality or age. The $N(<D)$ relationship of $a={0.96}$ implies that the randomised diameters are readily mimicking such an exponent. The rate of SNe occurring in the LMC is estimated to be $\sim$1 per 200 yr. The mean $\alpha$ of the entire LMC SNR population is $\alpha=-0.52$, which is typical of most SNRs. However, our estimates show a clear flattening of the synchrotron $\alpha$ as the remnants age. As predicted, our CC SNRs sample are significantly brighter radio emitters than the type Ia remnants. We also estimate the $\Sigma - D$ relation for the LMC to have a slope $\sim3.8$ which is comparable with other nearby galaxies. We also find the residency time of electrons in the galaxy ($4.0-14.3$ Myr), implying that SNRs should be the dominant mechanism for the production and acceleration of CRs.

We analyze the influence of fractal structure of the interstellar matter (ISM) density on the parameter values for the radio surface brightness to diameter ($\Sigma-D$) relation for supernovae remnants (SNRs). We model a dense ISM as a molecular cloud with fractal density structure. SNRs are modelled as spheres of different radius scattered in the modelled ISM. The surface brightness of the SNRs is calculated from the simple relation $\Sigma \propto \rho^{0.5}D^{-3.5}$ and also from the parametrized more general form $\Sigma \propto \rho^{\eta}D^{-\beta_0}$. Our results demonstrate that empirical $\Sigma-D$ slopes that are steeper than the ones derived from theory, might be partly explained with the fractal structure of the ambient medium into which SNRs expand. The slope of the $\Sigma-D$ relation steepens if the density of the regions where SNRs are formed is higher. The simple geometrical effects combined with the fractal structure of the ISM can contribute to a steeper empirical $\Sigma-D$ slopes, especially for older remnants, and this is more pronounced if $\Sigma$ has a stronger dependence on ambient density.

We present an isolated Milky Way-like simulation in GADGET2 N-body SPH code. The Galactic disk star formation rate (SFR) surface densities and stellar mass indicative of Solar neighbourhood are used as thresholds to model the distribution of stellar mass in life friendly environments. SFR and stellar component density are calculated averaging the GADGET2 particle properties on a 2D grid mapped on the Galactic plane. The peak values for possibly habitable stellar mass surface density move from $10$ to $15$ kpc cylindrical galactocentric distance in $10$ Gyr simulated time span. At $10$ Gyr the simulation results imply the following. Stellar particles which have spent almost all of their life time in habitable friendly conditions reside typically at $\sim16$ kpc from Galactic centre and are $\sim 3$ Gyr old. Stellar particles that have spent $\ge 90 \%$ of their $4-5$ Gyr long life time in habitable friendly conditions, are also predominantly found in the outskirts of the Galactic disk. Less then $1 \%$ of these particles can be found at a typical Solar system galactocentric distance of $8-10$ kpc. Our results imply that the evolution of an isolated spiral galaxy is likely to result in galactic civilizations emerging at the outskirts of the galactic disk around stellar hosts younger than the Sun.

We present updated empirical radio surface-brightness-to-diameter ($\Sigma-D$) relation for supernova remnants (SNRs) in our Galaxy. Our original calibration sample of Galactic SNRs with independently determined distances (Pavlovic et al. 2013, hereafter Paper I) is reconsidered and updated with data which became available in the past two years. The orthogonal fitting procedure and probability-density-function-based (PDF) method are applied to the calibration sample in the $\log \Sigma - \log D$ plane. Orthogonal regression keeps $\Sigma-D$ and $D-\Sigma$ relations invariant. Our previous Monte Carlo simulations verified that the slopes of the empirical $\Sigma-D$ relation should be determined by using orthogonal regression. Updated calibration sample contains 65 shell SNRs. 6 new Galactic SNRs are added to the sample from Paper I, one is omitted and distances are changed for 10 SNRs. The slope derived here is slightly steeper ($\beta \approx 5.2$) than $\Sigma-D$ slope in Paper I ($\beta \approx 4.8$). The PDF method relies on data points density maps which can provide more reliable calibrations that preserve more information contained in the calibration sample. We estimate distances to five new faint Galactic SNRs discovered for the first time by Canadian Galactic Plane Survey and obtained distances of 2.3, 4.0, 1.3, 2.9 and 4.7 kiloparsecs for G108.5+11.0, G128.5+2.6, G149.5+3.2, G150.8+3.8 and G160.1$-$1.1, respectively. The updated empirical relation is used to estimate distances of 160 shell Galactic SNRs and new results change their distance scales up to 15 per cent, compared to results from Paper I. The PDF calculation can provide even few times higher or lower values in comparison with orthogonal fit, as it uses totally different approach. However, in average, this difference is 32, 24 and 18 per cent for mode, median and mean distances.

Search for extraterrestrial life and intelligence constitutes one of the major endeavors in science, but has yet been quantitatively modeled only rarely and in a cursory and superficial fashion. We argue that probabilistic cellular automata (PCA) represent the best quantitative framework for modeling astrobiological history of the Milky Way and its Galactic Habitable Zone. The relevant astrobiological parameters are to be modeled as the elements of the input probability matrix for the PCA kernel. With the underlying simplicity of the cellular automata constructs, this approach enables a quick analysis of large and ambiguous input parameters' space. We perform a simple clustering analysis of typical astrobiological histories and discuss the relevant boundary conditions of practical importance for planning and guiding actual empirical astrobiological and SETI projects. In addition to showing how the present framework is adaptable to more complex situations and updated observational databases from current and near-future space missions, we demonstrate how numerical results could offer a cautious rationale for continuation of practical SETI searches.

We present an extended analysis of the relation between radio surface brightness and diameter -- the so-called $\Sigma-D$ relation for planetary nebulae (PNe). We revise our previous derivation of the theoretical $\Sigma-D$ relation for the evolution of bremsstrahlung surface brightness in order to include the influence of the fast wind from the central star. Different theoretical forms are derived: $\Sigma \propto D^{-1}$ for the first and second phases of evolution and $\Sigma\propto D^{-3}$ for the final stage of evolution. Also, we analyzed several different Galactic PN samples. All samples are influenced by severe selection effects, but Malmquist bias seems to be less influential here than in the supernova remnant (SNR) samples. We derived empirical $\Sigma-D$ relations for 27 sample sets using 6 updated PN papers from which an additional 21 new sets were extracted. Twenty four of these have a trivial form of $\beta \approx 2$. However, we obtain one empirical $\Sigma-D$ relation that may be useful for determining distances to PNe. This relation is obtained by extracting a recent nearby (< 1 kpc) Galactic PN sample.