Astrometric detection involves a precise measurement of stellar positions, and is widely regarded as the leading concept presently ready to find earth-mass planets in temperate orbits around nearby sun-like stars. The TOLIMAN space telescope[39] is a low-cost, agile mission concept dedicated to narrow-angle astrometric monitoring of bright binary stars. In particular the mission will be optimised to search for habitable-zone planets around Alpha Centauri AB. If the separation between these two stars can be monitored with sufficient precision, tiny perturbations due to the gravitational tug from an unseen planet can be witnessed and, given the configuration of the optical system, the scale of the shifts in the image plane are about one millionth of a pixel. Image registration at this level of precision has never been demonstrated (to our knowledge) in any setting within science. In this paper we demonstrate that a Deep Convolutional Auto-Encoder is able to retrieve such a signal from simplified simulations of the TOLIMAN data and we present the full experimental pipeline to recreate out experiments from the simulations to the signal analysis. In future works, all the more realistic sources of noise and systematic effects present in the real-world system will be injected into the simulations.

Modern astronomy has been rapidly increasing our ability to see deeper into the universe, acquiring enormous samples of cosmic populations. Gaining astrophysical insights from these datasets requires a wide range of sophisticated statistical and machine learning methods. Long-standing problems in cosmology include characterization of galaxy clustering and estimation of galaxy distances from photometric colors. Bayesian inference, central to linking astronomical data to nonlinear astrophysical models, addresses problems in solar physics, properties of star clusters, and exoplanet systems. Likelihood-free methods are growing in importance. Detection of faint signals in complicated noise is needed to find periodic behaviors in stars and detect explosive gravitational wave events. Open issues concern treatment of heteroscedastic measurement errors and understanding probability distributions characterizing astrophysical systems. The field of astrostatistics needs increased collaboration with statisticians in the design and analysis stages of research projects, and to jointly develop new statistical methodologies. Together, they will draw more astrophysical insights into astronomical populations and the cosmos itself.

Cosmic voids and their corresponding redshift-aggregated projections of mass densities, known as troughs, play an important role in our attempt to model the large-scale structure of the Universe. Understanding these structures leads to tests comparing the standard model with alternative cosmologies, constraints on the dark energy equation of state, and provides evidence to differentiate among gravitational theories. In this paper, we extend the subspace-constrained mean shift algorithm, a recently introduced method to estimate density ridges, and apply it to 2D weak-lensing mass density maps from the Dark Energy Survey Y1 data release to identify curvilinear filamentary structures. We compare the obtained ridges with previous approaches to extract trough structure in the same data, and apply curvelets as an alternative wavelet-based method to constrain densities. We then invoke the Wasserstein distance between noisy and noiseless simulations to validate the denoising capabilities of our method. Our results demonstrate the viability of ridge estimation as a precursor for denoising weak lensing quantities to recover the large-scale structure, paving the way for a more versatile and effective search for troughs.

Big bang nucleosynthesis provides the earliest probe of standard model physics, at a time when the universe was less than a thousand seconds old. It determines the abundances of the lightest nuclides, which give rise to the subsequent history of the visible matter in the Universe. This work derives new $^7$Be(n,p)$^7$Li thermonuclear reaction rates based on all available experimental information. This reaction sensitively impacts the primordial abundances of $^{7}$Be and $^7$Li during big bang nucleosynthesis. We critically evaluate all available data and disregard experimental results that are questionable. For the nuclear model, we adopt an incoherent sum of single-level, two-channel R-matrix approximation expressions, which are implemented into a hierarchical Bayesian model, to analyze the remaining six data sets we deem most reliable. In the fitting of the data, we consistently model all known sources of uncertainty, including discrepant absolute normalizations of different data sets, and also take the variation of the neutron and proton channel radii into account, hence providing less biased estimates of the $^7$Be(n,p)$^7$Li thermonuclear rates. From the resulting posteriors, we extract R-matrix parameters ($E_r$, $\gamma^2_n$, $\gamma^2_p$) and derive excitation energies, partial and total widths. Our fit is sensitive to the contributions of the first three levels above the neutron threshold. Reaction rates were computed by integrating 10,000 samples of the reduced cross section. Our $^7$Be(n,p)$^7$Li thermonuclear rates have uncertainties between 1.5% and 2.0% at temperatures of $\leq$1 GK. We compare our rates to previous results and find that the $^7$Be(n,p)$^7$Li rates most commonly used in big bang simulations have too optimistic uncertainties.

The so-called ultraviolet (UV) upturn of elliptical galaxies is a phenomenon characterised by the up-rise of their fluxes in bluer wavelengths, typically in the 1,200-2,500A range. This work aims at estimating the rate of occurrence of the UV upturn over the entire red-sequence population of galaxies that show significant UV emission. This assessment is made considering it as function of three parameters: redshift, stellar mass, and -- what may seem counter-intuitive at first -- emission-line classification. We built a multiwavelength spectro-photometric catalogue from the Galaxy Mass Assembly survey, together with aperture-matched data from Galaxy Evolution Explorer Medium-Depth Imaging Survey (MIS) and Sloan Digital Sky Survey, covering the redshift range between 0.06 and 0.40. From this sample, we analyse the UV emission among UV bright galaxies, by selecting those that occupy the red-sequence locus in the (NUV-r) x (FUV-NUV) chart; then, we stratify the sample by their emission-line classes. To that end, we make use of emission-line diagnostic diagrams, focusing the analysis in retired/passive lineless galaxies. Then, a Bayesian logistic model was built to simultaneously deal with the effects of all galaxy properties (including emission-line classification or lack thereof). The main results show that retired/passive systems host an up-rise in the fraction of UV upturn or redshifts between 0.06 and 0.25, followed by an in-fall up to 0.35. Additionally, we show that the fraction of UV upturn hosts rises with increasing stellar mass.

In this paper we report the current status of a new research program. The primary goal of the "Vanishing & Appearing Sources during a Century of Observations" (VASCO) project is to search for vanishing and appearing sources using existing survey data to find examples of exceptional astrophysical transients. The implications of finding such objects extend from traditional astrophysics fields to the more exotic searches for evidence of technologically advanced civilizations. In this first paper we present new, deeper observations of the tentative candidate discovered by Villarroel et al. (2016). We then perform the first searches for vanishing objects throughout the sky by comparing 600 million objects from the US Naval Observatory Catalogue (USNO) B1.0 down to a limiting magnitude of $\sim 20 - 21$ with the recent Pan-STARRS Data Release-1 (DR1) with a limiting magnitude of $\sim$ 23.4. We find about 150,000 preliminary candidates that do not have any Pan-STARRS counterpart within a 30 arcsec radius. We show that these objects are redder and have larger proper motions than typical USNO objects. We visually examine the images for a subset of about 24,000 candidates, superseding the 2016 study with a sample ten times larger. We find about $\sim$ 100 point sources visible in only one epoch in the red band of the USNO which may be of interest in searches for strong M dwarf flares, high-redshift supernovae or other catagories of unidentified red transients.

We present a star/galaxy classification for the Southern Photometric Local Universe Survey (S-PLUS), based on a Machine Learning approach: the Random Forest algorithm. We train the algorithm using the S-PLUS optical photometry up to $r$=21, matched to SDSS/DR13, and morphological parameters. The metric of importance is defined as the relative decrease of the initial accuracy when all correlations related to a certain feature is vanished. In general, the broad photometric bands presented higher importance when compared to narrow ones. The influence of the morphological parameters has been evaluated training the RF with and without the inclusion of morphological parameters, presenting accuracy values of 95.0\% and 88.1\%, respectively. Particularly, the morphological parameter {\rm FWHM/PSF} performed the highest importance over all features to distinguish between stars and galaxies, indicating that it is crucial to classify objects into stars and galaxies. We investigate the misclassification of stars and galaxies in the broad-band colour-colour diagram $(g-r)$ versus $(r-i)$. The morphology can notably improve the classification of objects at regions in the diagram where the misclassification was relatively high. Consequently, it provides cleaner samples for statistical studies. The expected contamination rate of red galaxies as a function of the redshift is estimated, providing corrections for red galaxy samples. The classification of QSOs as extragalactic objects is slightly better using photometric-only case. An extragalactic point-source catalogue is provided using the classification without any morphology feature (only the SED information) with additional constraints on photometric redshifts and {\rm FWHM/PSF} values.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

In this paper we present a thorough discussion about the photometric redshift (photo-z) performance of the Southern Photometric Local Universe Survey (S-PLUS). This survey combines a 7 narrow + 5 broad passband filter system, with a typical photometric-depth of r$\sim$21 AB. For this exercise, we utilize the Data Release 1 (DR1), corresponding to 336 deg$^{2}$ from the Stripe-82 region. We rely on the \texttt{BPZ2} code to compute our estimates, using a new library of SED models, which includes additional templates for quiescent galaxies. When compared to a spectroscopic redshift control sample of $\sim$100k galaxies, we find a precision of $\sigma_{z}<$0.8\%, $<$2.0\% or $<$3.0\% for galaxies with magnitudes r$<$17, $<$19 and $<$21, respectively. A precision of 0.6\% is attained for galaxies with the highest \texttt{Odds} values. These estimates have a negligible bias and a fraction of catastrophic outliers inferior to 1\%. We identify a redshift window (i.e., 0.26$<z<$0.32) where our estimates double their precision, due to the simultaneous detection of two emission-lines in two distinct narrow-bands; representing a window opportunity to conduct statistical studies such as luminosity functions. We forecast a total of $\sim$2M, $\sim$16M and $\sim$32M galaxies in the S-PLUS survey with a photo-z precision of $\sigma_{z}<$1.0\%, $<$2.0\% and $<$2.5\% after observing 8000 $deg^{2}$. We also derive redshift Probability Density Functions, proving their reliability encoding redshift uncertainties and their potential recovering the $n(z)$ of galaxies at $z<0.4$, with an unprecedented precision for a photometric survey in the southern hemisphere.

This work determines the degree to which a standard Lambda-CDM analysis based on type Ia supernovae can identify deviations from a cosmological constant in the form of a redshift-dependent dark energy equation of state w(z). We introduce and apply a novel random curve generator to simulate instances of w(z) from constraint families with increasing distinction from a cosmological constant. After producing a series of mock catalogs of binned type Ia supernovae corresponding to each w(z) curve, we perform a standard Lambda-CDM analysis to estimate the corresponding posterior densities of the absolute magnitude of type Ia supernovae, the present-day matter density, and the equation of state parameter. Using the Kullback-Leibler divergence between posterior densities as a difference measure, we demonstrate that a standard type Ia supernova cosmology analysis has limited sensitivity to extensive redshift dependencies of the dark energy equation of state. In addition, we report that larger redshift-dependent departures from a cosmological constant do not necessarily manifest easier-detectable incompatibilities with the Lambda-CDM model. Our results suggest that physics beyond the standard model may simply be hidden in plain sight.

The new generation of deep photometric surveys requires unprecedentedly precise shape and photometry measurements of billions of galaxies to achieve their main science goals. At such depths, one major limiting factor is the blending of galaxies due to line-of-sight projection, with an expected fraction of blended galaxies of up to 50%. Current deblending approaches are in most cases either too slow or not accurate enough to reach the level of requirements. This work explores the use of deep neural networks to estimate the photometry of blended pairs of galaxies in monochrome space images, similar to the ones that will be delivered by the Euclid space telescope. Using a clean sample of isolated galaxies from the CANDELS survey, we artificially blend them and train two different network models to recover the photometry of the two galaxies. We show that our approach can recover the original photometry of the galaxies before being blended with $\sim$7% accuracy without any human intervention and without any assumption on the galaxy shape. This represents an improvement of at least a factor of 4 compared to the classical SExtractor approach. We also show that forcing the network to simultaneously estimate a binary segmentation map results in a slightly improved photometry. All data products and codes will be made public to ease the comparison with other approaches on a common data set.

Open clusters (OCs) are popular tracers of the structure and evolutionary history of the Galactic disk. The OC population is often considered to be complete within 1.8 kpc of the Sun. The recent Gaia Data Release 2 (DR2) allows the latter claim to be challenged. We perform a systematic search for new OCs in the direction of Perseus using precise and accurate astrometry from Gaia DR2. We implement a coarse-to-fine search method. First, we exploit spatial proximity using a fast density-aware partitioning of the sky via a k-d tree in the spatial domain of Galactic coordinates, (l, b). Secondly, we employ a Gaussian mixture model in the proper motion space to quickly tag fields around OC candidates. Thirdly, we apply an unsupervised membership assignment method, UPMASK, to scrutinise the candidates. We visually inspect colour-magnitude diagrams to validate the detected objects. Finally, we perform a diagnostic to quantify the significance of each identified overdensity in proper motion and in parallax space We report the discovery of 41 new stellar clusters. This represents an increment of at least 20% of the previously known OC population in this volume of the Milky Way. We also report on the clear identification of NGC 886, an object previously considered an asterism. This letter challenges the previous claim of a near-complete sample of open clusters up to 1.8 kpc. Our results reveal that this claim requires revision, and a complete census of nearby open clusters is yet to be found.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

We developed a hierarchical Bayesian framework to estimate S-factors and thermonuclear rates for the $^3$He(d,p)$^{4}$He reaction, which impacts the primordial abundances of $^3$He and $^7$Li. The available data are evaluated and all direct measurements are taken into account in our analysis for which we can estimate separate uncertainties for systematic and statistical effects. For the nuclear reaction model, we adopt a single-level, two-channel approximation of R-matrix theory, suitably modified to take the effects of electron screening at lower energies into account. Apart from the usual resonance parameters (resonance location and reduced widths for the incoming and outgoing reaction channel), we include for the first time the channel radii and boundary condition in the fitting process. Our new analysis of the $^3$He(d,p)$^{4}$He S-factor data results in improved estimates for the thermonuclear rates. This work represents the first nuclear rate evaluation using the R-matrix theory embedded into a hierarchical Bayesian framework, properly accounting for all known sources of uncertainty. Therefore, it provides a test bed for future studies of more complex reactions.

The $^3$H(d,n)$^4$He reaction has a large low-energy cross section and will likely be utilized in future commercial fusion reactors. This reaction also takes place during big bang nucleosynthesis. Studies of both scenarios require accurate and precise fusion rates. To this end, we implement a one-level, two-channel R-matrix approximation into a Bayesian model. Our main goals are to predict reliable astrophysical S-factors and to estimate R-matrix parameters using the Bayesian approach. All relevant parameters are sampled in our study, including the channel radii, boundary condition parameters, and data set normalization factors. In addition, we take uncertainties in both measured bombarding energies and S-factors rigorously into account. Thermonuclear rates and reactivities of the $^3$H(d,n)$^4$He reaction are derived by numerically integrating the Bayesian S-factor samples. The present reaction rate uncertainties at temperatures between $1.0$ MK and $1.0$ GK are in the range of 0.2% to 0.6%. Our reaction rates differ from previous results by 2.9% near 1.0 GK. Our reactivities are smaller than previous results, with a maximum deviation of 2.9% near a thermal energy of $4$ keV. The present rate or reactivity uncertainties are more reliable compared to previous studies that did not include the channel radii, boundary condition parameters, and data set normalization factors in the fitting. Finally, we investigate previous claims of electron screening effects in the published $^3$H(d,n)$^4$He data. No such effects are evident and only an upper limit for the electron screening potential can be obtained.

The dark ages of the Universe end with the formation of the first generation of stars residing in primeval galaxies. These objects were the first to produce ultraviolet ionizing photons in a period when the cosmic gas changed from a neutral state to an ionized one, known as Epoch of Reionization (EoR). A pivotal aspect to comprehend the EoR is to probe the intertwined relationship between the fraction of ionizing photons capable to escape dark haloes, also known as the escape fraction ($f_{esc}$), and the physical properties of the galaxy. This work develops a sound statistical model suitable to account for such non-linear relationships and the non-Gaussian nature of $f_{esc}$. This model simultaneously estimates the probability that a given primordial galaxy starts the ionizing photon production and estimates the mean level of the $f_{esc}$ once it is triggered. The model was employed in the First Billion Years simulation suite, from which we show that the baryonic fraction and the rate of ionizing photons appear to have a larger impact on $f_{esc}$ than previously thought. A naive univariate analysis of the same problem would suggest smaller effects for these properties and a much larger impact for the specific star formation rate, which is lessened after accounting for other galaxy properties and non-linearities in the statistical model.

We report a framework for spectroscopic follow-up design for optimizing supernova photometric classification. The strategy accounts for the unavoidable mismatch between spectroscopic and photometric samples, and can be used even in the beginning of a new survey -- without any initial training set. The framework falls under the umbrella of active learning (AL), a class of algorithms that aims to minimize labelling costs by identifying a few, carefully chosen, objects which have high potential in improving the classifier predictions. As a proof of concept, we use the simulated data released after the Supernova Photometric Classification Challenge (SNPCC) and a random forest classifier. Our results show that, using only 12\% the number of training objects in the SNPCC spectroscopic sample, this approach is able to double purity results. Moreover, in order to take into account multiple spectroscopic observations in the same night, we propose a semi-supervised batch-mode AL algorithm which selects a set of $N=5$ most informative objects at each night. In comparison with the initial state using the traditional approach, our method achieves 2.3 times higher purity and comparable figure of merit results after only 180 days of observation, or 800 queries (73% of the SNPCC spectroscopic sample size). Such results were obtained using the same amount of spectroscopic time necessary to observe the original SNPCC spectroscopic sample, showing that this type of strategy is feasible with current available spectroscopic resources. The code used in this work is available in the COINtoolbox: https://github.com/COINtoolbox/ActSNClass .

Machine learning techniques offer a precious tool box for use within astronomy to solve problems involving so-called big data. They provide a means to make accurate predictions about a particular system without prior knowledge of the underlying physical processes of the data. In this article, and the companion papers of this series, we present the set of Generalized Linear Models (GLMs) as a fast alternative method for tackling general astronomical problems, including the ones related to the machine learning paradigm. To demonstrate the applicability of GLMs to inherently positive and continuous physical observables, we explore their use in estimating the photometric redshifts of galaxies from their multi-wavelength photometry. Using the gamma family with a log link function we predict redshifts from the PHoto-z Accuracy Testing simulated catalogue and a subset of the Sloan Digital Sky Survey from Data Release 10. We obtain fits that result in catastrophic outlier rates as low as ~1% for simulated and ~2% for real data. Moreover, we can easily obtain such levels of precision within a matter of seconds on a normal desktop computer and with training sets that contain merely thousands of galaxies. Our software is made publicly available as an user-friendly package developed in Python, R and via an interactive web application (https://cosmostatisticsinitiative.shinyapps.io/CosmoPhotoz). This software allows users to apply a set of GLMs to their own photometric catalogues and generates publication quality plots with minimum effort from the user. By facilitating their ease of use to the astronomical community, this paper series aims to make GLMs widely known and to encourage their implementation in future large-scale projects, such as the Large Synoptic Survey Telescope.

Astronomical observations of extended sources, such as cubes of integral field spectroscopy (IFS), encode auto-correlated spatial structures that cannot be optimally exploited by standard methodologies. This work introduces a novel technique to model IFS datasets, which treats the observed galaxy properties as realizations of an unobserved Gaussian Markov random field. The method is computationally efficient, resilient to the presence of low-signal-to-noise regions, and uses an alternative to Markov Chain Monte Carlo for fast Bayesian inference, the Integrated Nested Laplace Approximation (INLA). As a case study, we analyse 721 IFS data cubes of nearby galaxies from the CALIFA and PISCO surveys, for which we retrieve the maps of the following physical properties: age, metallicity, mass and extinction. The proposed Bayesian approach, built on a generative representation of the galaxy properties, enables the creation of synthetic images, recovery of areas with bad pixels, and an increased power to detect structures in datasets subject to substantial noise and/or sparsity of sampling. A snippet code to reproduce the analysis of this paper is available in the COIN toolbox, together with the field reconstructions of the CALIFA and PISCO samples.

We introduce a value-added full-sky catalogue of galaxies, named as Galaxy List for the Advanced Detector Era, or GLADE. The purpose of this catalogue is to (i) help identifications of host candidates for gravitational-wave events, (ii) support target selections for electromagnetic follow-up observations of gravitational-wave candidates, (iii) provide input data on the matter distribution of the local universe for astrophysical or cosmological simulations, and (iv) help identifications of host candidates for poorly localised electromagnetic transients, such as gamma-ray bursts observed with the InterPlanetary Network. Both being potential hosts of astrophysical sources of gravitational waves, GLADE includes inactive and active galaxies as well. GLADE was constructed by cross-matching and combining data from five separate (but not independent) astronomical catalogues: GWGC, 2MPZ, 2MASS XSC, HyperLEDA and SDSS-DR12Q. GLADE is complete up to $d_L = 37^{+3}_{-4}$ Mpc in terms of the cumulative B-band luminosity of galaxies within luminosity distance $d_L$, and contains all of the brightest galaxies giving half of the total B-band luminosity up to $d_L = 91$ Mpc. As B-band luminosity is expected to be a tracer of binary neutron star mergers (currently the prime targets of joint GW+EM detections), our completeness measures can be used as estimations of completeness for containing all binary neutron star merger hosts in the local universe.

We present a review of data types and statistical methods often encountered in astronomy. The aim is to provide an introduction to statistical applications in astronomy for statisticians and computer scientists. We highlight the complex, often hierarchical, nature of many astronomy inference problems and advocate for cross-disciplinary collaborations to address these challenges.

Approximate Bayesian Computation (ABC) enables parameter inference for complex physical systems in cases where the true likelihood function is unknown, unavailable, or computationally too expensive. It relies on the forward simulation of mock data and comparison between observed and synthetic catalogues. Here we present cosmoabc, a Python ABC sampler featuring a Population Monte Carlo (PMC) variation of the original ABC algorithm, which uses an adaptive importance sampling scheme. The code is very flexible and can be easily coupled to an external simulator, while allowing to incorporate arbitrary distance and prior functions. As an example of practical application, we coupled cosmoabc with the numcosmo library and demonstrate how it can be used to estimate posterior probability distributions over cosmological parameters based on measurements of galaxy clusters number counts without computing the likelihood function. cosmoabc is published under the GPLv3 license on PyPI and GitHub and documentation is available at http://goo.gl/SmB8EX

We developed a Bayesian framework to determine in a robust way the relation between velocity dispersions and chemical abundances in a sample of stars. Our modelling takes into account the uncertainties in the chemical and kinematic properties. We make use of RAVE DR5 radial velocities and abundances together with Gaia DR1 proper motions and parallaxes (when possible, otherwise UCAC4 data is used). We found that, in general, the velocity dispersions increase with decreasing [Fe/H] and increasing [Mg/Fe]. A possible decrease in velocity dispersion for stars with high [Mg/Fe] is a property of a negligible fraction of stars and hardly a robust result. At low [Fe/H] and high [Mg/Fe] the sample is incomplete, affected by biases, and likely not representative of the underlying stellar population.

We invoke a Gaussian mixture model (GMM) to jointly analyse two traditional emission-line classification schemes of galaxy ionization sources: the Baldwin-Phillips-Terlevich (BPT) and $\rm W_{H\alpha}$ vs. [NII]/H$\alpha$ (WHAN) diagrams, using spectroscopic data from the Sloan Digital Sky Survey Data Release 7 and SEAGal/STARLIGHT datasets. We apply a GMM to empirically define classes of galaxies in a three-dimensional space spanned by the $\log$ [OIII]/H$\beta$, $\log$ [NII]/H$\alpha$, and $\log$ EW(H${\alpha}$), optical parameters. The best-fit GMM based on several statistical criteria suggests a solution around four Gaussian components (GCs), which are capable to explain up to 97 per cent of the data variance. Using elements of information theory, we compare each GC to their respective astronomical counterpart. GC1 and GC4 are associated with star-forming galaxies, suggesting the need to define a new starburst subgroup. GC2 is associated with BPT's Active Galaxy Nuclei (AGN) class and WHAN's weak AGN class. GC3 is associated with BPT's composite class and WHAN's strong AGN class. Conversely, there is no statistical evidence -- based on four GCs -- for the existence of a Seyfert/LINER dichotomy in our sample. Notwithstanding, the inclusion of an additional GC5 unravels it. The GC5 appears associated to the LINER and Passive galaxies on the BPT and WHAN diagrams respectively. Subtleties aside, we demonstrate the potential of our methodology to recover/unravel different objects inside the wilderness of astronomical datasets, without lacking the ability to convey physically interpretable results. The probabilistic classifications from the GMM analysis are publicly available within the COINtoolbox (https://cointoolbox.github.io/GMM\_Catalogue/).

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Two of the main problems encountered in the development and accurate validation of photometric redshift (photo-z) techniques are the lack of spectroscopic coverage in feature space (e.g. colours and magnitudes) and the mismatch between photometric error distributions associated with the spectroscopic and photometric samples. Although these issues are well known, there is currently no standard benchmark allowing a quantitative analysis of their impact on the final photo-z estimation. In this work, we present two galaxy catalogues, Teddy and Happy, built to enable a more demanding and realistic test of photo-z methods. Using photometry from the Sloan Digital Sky Survey and spectroscopy from a collection of sources, we constructed datasets which mimic the biases between the underlying probability distribution of the real spectroscopic and photometric sample. We demonstrate the potential of these catalogues by submitting them to the scrutiny of different photo-z methods, including machine learning (ML) and template fitting approaches. Beyond the expected bad results from most ML algorithms for cases with missing coverage in feature space, we were able to recognize the superiority of global models in the same situation and the general failure across all types of methods when incomplete coverage is convoluted with the presence of photometric errors - a data situation which photo-z methods were not trained to deal with up to now and which must be addressed by future large scale surveys. Our catalogues represent the first controlled environment allowing a straightforward implementation of such tests. The data are publicly available within the COINtoolbox (https://github.com/COINtoolbox/photoz_catalogues).

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

We revisit the present status of the stiffness of the supranuclear equations of state, particularly the solutions that increase the stiffness in the presence of hyperons, the putative transition to a quark matter phase and the robustness of massive compact star observations.

We developed a hierarchical Bayesian model (HBM) to investigate how the presence of Seyfert activity relates to their environment, herein represented by the galaxy cluster mass, $M_{200}$, and the normalized cluster-centric distance, $r/r_{200}$. We achieved this by constructing an unbiased sample of galaxies from the Sloan Digital Sky Survey, with morphological classifications provided by the Galaxy Zoo Project. A propensity score matching approach is introduced to control for the effects of confounding variables: stellar mass, galaxy colour, and star formation rate. The connection between Seyfert-activity and environmental properties in the de-biased sample is modelled within an HBM framework using the so-called logistic regression technique, suitable for the analysis of binary data (e.g., whether or not a galaxy hosts an AGN). Unlike standard ordinary least square fitting methods, our methodology naturally allows modelling the probability of Seyfert-AGN activity in galaxies on their natural scale, i.e. as a binary variable. Furthermore, we demonstrate how an HBM can incorporate information of each particular galaxy morphological type in a unified framework. In elliptical galaxies, our analysis indicates a strong correlation of Seyfert-AGN activity with $r/r_{200}$, and a weaker correlation with the mass of the host. In spiral galaxies these trends do not appear, suggesting that the link between Seyfert activity and the properties of spiral galaxies are independent of the environment.

The existence of multiple subclasses of type Ia supernovae (SNeIa) has been the subject of great debate in the last decade. One major challenge inevitably met when trying to infer the existence of one or more subclasses is the time consuming, and subjective, process of subclass definition. In this work, we show how machine learning tools facilitate identification of subtypes of SNeIa through the establishment of a hierarchical group structure in the continuous space of spectral diversity formed by these objects. Using Deep Learning, we were capable of performing such identification in a 4 dimensional feature space (+1 for time evolution), while the standard Principal Component Analysis barely achieves similar results using 15 principal components. This is evidence that the progenitor system and the explosion mechanism can be described by a small number of initial physical parameters. As a proof of concept, we show that our results are in close agreement with a previously suggested classification scheme and that our proposed method can grasp the main spectral features behind the definition of such subtypes. This allows the confirmation of the velocity of lines as a first order effect in the determination of SNIa subtypes, followed by 91bg-like events. Given the expected data deluge in the forthcoming years, our proposed approach is essential to allow a quick and statistically coherent identification of SNeIa subtypes (and outliers). All tools used in this work were made publicly available in the Python package Dimensionality Reduction And Clustering for Unsupervised Learning in Astronomy (DRACULA) and can be found within COINtoolbox (https://github.com/COINtoolbox/DRACULA).

The motion of a point like object of mass $M$ passing through the background potential of massive collisionless particles ($m << M$) suffers a steady deceleration named dynamical friction. In his classical work, Chandrasekhar assumed a Maxwellian velocity distribution in the halo and neglected the self gravity of the wake induced by the gravitational focusing of the mass $M$. In this paper, by relaxing the validity of the Maxwellian distribution due to the presence of long range forces, we derive an analytical formula for the dynamical friction in the context of the $q$-nonextensive kinetic theory. In the extensive limiting case ($q = 1$), the classical Gaussian Chandrasekhar result is recovered. As an application, the dynamical friction timescale for Globular Clusters spiraling to the galactic center is explicitly obtained. Our results suggest that the problem concerning the large timescale as derived by numerical $N$-body simulations or semi-analytical models can be understood as a departure from the standard extensive Maxwellian regime as measured by the Tsallis nonextensive $q$-parameter.

The first supernovae will soon be visible at the edge of the observable universe, revealing the birthplaces of Population III stars. With upcoming near-infrared missions, a broad analysis of the detectability of high-$z$ supernovae is paramount. We combine cosmological and radiation transport simulations, instrument specifications, and survey strategies to create synthetic observations of primeval core-collapse, Type IIn and pair-instability supernovae with the James Webb Space Telescope ($JWST$). We show that a dedicated observational campaign with the $JWST$ can detect up to $\sim 15$ pair-instability explosions, $\sim 300$ core-collapse supernovae, but less than one Type IIn explosion per year, depending on the Population III star formation history. Our synthetic survey also shows that $\approx 1-2 \times10^2$ supernovae detections, depending on the accuracy of the classification, are sufficient to discriminate between a Salpeter and flat mass distribution for high redshift stars with a confidence level greater than 99.5 per cent. We discuss how the purity of the sample affects our results and how supervised learning methods may help to discriminate between CC and PI SNe.

This paper is devoted to discuss the difference in the thermodynamic entropy budget {\it per baryon} in each type of stellar object found in Universe. We track and discuss the actual {\it decrease} of the stored baryonic thermodynamic entropy from the most primitive molecular cloud up to the final fate of matter in the black holes, passing through evolved states of matter as found in white dwarfs and neutron stars. We then discuss the case of actual stars of different masses throughout their {\it evolution}, clarifying the role of virial equilibrium condition for the decrease of the entropy and related issues. Finally, we discuss how gravity ultimately drives composition, hence structural changes along the stellar evolution all the way until the ultimate collapse to black holes, which may increase dramatically their entropy because of the gravitational contribution itself.

In this paper, the third in a series illustrating the power of generalized linear models (GLMs) for the astronomical community, we elucidate the potential of the class of GLMs which handles count data. The size of a galaxy's globular cluster population $N_{\rm GC}$ is a prolonged puzzle in the astronomical literature. It falls in the category of count data analysis, yet it is usually modelled as if it were a continuous response variable. We have developed a Bayesian negative binomial regression model to study the connection between $N_{\rm GC}$ and the following galaxy properties: central black hole mass, dynamical bulge mass, bulge velocity dispersion, and absolute visual magnitude. The methodology introduced herein naturally accounts for heteroscedasticity, intrinsic scatter, errors in measurements in both axes (either discrete or continuous), and allows modelling the population of globular clusters on their natural scale as a non-negative integer variable. Prediction intervals of 99% around the trend for expected $N_{\rm GC}$comfortably envelope the data, notably including the Milky Way, which has hitherto been considered a problematic outlier. Finally, we demonstrate how random intercept models can incorporate information of each particular galaxy morphological type. Bayesian variable selection methodology allows for automatically identifying galaxy types with different productions of GCs, suggesting that on average S0 galaxies have a GC population 35% smaller than other types with similar brightness.

Machine learning techniques offer a plethora of opportunities in tackling big data within the astronomical community. We present the set of Generalized Linear Models as a fast alternative for determining photometric redshifts of galaxies, a set of tools not commonly applied within astronomy, despite being widely used in other professions. With this technique, we achieve catastrophic outlier rates of the order of ~1%, that can be achieved in a matter of seconds on large datasets of size ~1,000,000. To make these techniques easily accessible to the astronomical community, we developed a set of libraries and tools that are publicly available.

We present AMADA, an interactive web application to analyse multidimensional datasets. The user uploads a simple ASCII file and AMADA performs a number of exploratory analysis together with contemporary visualizations diagnostics. The package performs a hierarchical clustering in the parameter space, and the user can choose among linear, monotonic or non-linear correlation analysis. AMADA provides a number of clustering visualization diagnostics such as heatmaps, dendrograms, chord diagrams, and graphs. In addition, AMADA has the option to run a standard or robust principal components analysis, displaying the results as polar bar plots. The code is written in R and the web interface was created using the Shiny framework. AMADA source-code is freely available at https://goo.gl/KeSPue, and the shiny-app at http://goo.gl/UTnU7I.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

Impulsive 30 THz continuum bursts have been recently observed in solar flares, utilizing small telescopes with a unique and relatively simple optical setup concept. The most intense burst was observed together with a GOES X2 class event on October 27, 2014, also detected at two sub-THz frequencies, RHESSI X-rays and SDO/HMI and EUV. It exhibits strikingly good correlation in time and in space with white light flare emission. It is likely that this association may prove to be very common. All three 30 THz events recently observed exhibited intense fluxes in the range of 104 solar flux units, considerably larger than those measured for the same events at microwave and sub-mm wavelengths. The 30 THz burst emission might be part of the same spectral burst component found at sub-THz frequencies. The 30 THz solar bursts open a promising new window for the study of flares at their origin

Revealing hidden patterns in astronomical data is often the path to fundamental scientific breakthroughs; meanwhile the complexity of scientific inquiry increases as more subtle relationships are sought. Contemporary data analysis problems often elude the capabilities of classical statistical techniques, suggesting the use of cutting edge statistical methods. In this light, astronomers have overlooked a whole family of statistical techniques for exploratory data analysis and robust regression, the so-called Generalized Linear Models (GLMs). In this paper -- the first in a series aimed at illustrating the power of these methods in astronomical applications -- we elucidate the potential of a particular class of GLMs for handling binary/binomial data, the so-called logit and probit regression techniques, from both a maximum likelihood and a Bayesian perspective. As a case in point, we present the use of these GLMs to explore the conditions of star formation activity and metal enrichment in primordial minihaloes from cosmological hydro-simulations including detailed chemistry, gas physics, and stellar feedback. We predict that for a dark mini-halo with metallicity $\approx 1.3 \times 10^{-4} Z_{\bigodot}$, an increase of $1.2 \times 10^{-2}$ in the gas molecular fraction, increases the probability of star formation occurrence by a factor of 75%. Finally, we highlight the use of receiver operating characteristic curves as a diagnostic for binary classifiers, and ultimately we use these to demonstrate the competitive predictive performance of GLMs against the popular technique of artificial neural networks.

We report the first analytical expression purely constructed by a machine to determine photometric redshifts ($z_{\rm phot}$) of galaxies. A simple and reliable functional form is derived using $41,214$ galaxies from the Sloan Digital Sky Survey Data Release 10 (SDSS-DR10) spectroscopic sample. The method automatically dropped the $u$ and $z$ bands, relying only on $g$, $r$ and $i$ for the final solution. Applying this expression to other $1,417,181$ SDSS-DR10 galaxies, with measured spectroscopic redshifts ($z_{\rm spec}$), we achieved a mean $\langle (z_{\rm phot} - z_{\rm spec})/(1+z_{\rm spec})\rangle\lesssim 0.0086$ and a scatter $\sigma_{(z_{\rm phot} - z_{\rm spec})/(1+z_{\rm spec})}\lesssim 0.045$ when averaged up to $z \lesssim 1.0$. The method was also applied to the PHAT0 dataset, confirming the competitiveness of our results when faced with other methods from the literature. This is the first use of symbolic regression in cosmology, representing a leap forward in astronomy-data-mining connection.

Besides its major objective tuned to the detection of the stellar galactic population the Gaia mission experiment will also observe a large number of galaxies. In this work we intend to evaluate the number and the characteristics of the galaxies that will effectively pass the onboard selection algorithm of Gaia. The detection of objects in Gaia will be performed in a section of the focal plane known as the Sky Mapper. Taking into account the Video Processing Algorithm criterion of detection and considering the known light profiles of discs and bulges galaxies we assess the number and the type of extra-galactic objects that will be observed by Gaia. We show that the stellar disk population of galaxies will be very difficult to be observed. On the contrary the spheroidal component of elliptical galaxies and bulges having higher central surface brightness and steeper brightness profile will be more easy to be detected. We estimate that most of the 20 000 elliptical population of nearby galaxies inside the local region up to 170 Mpc are in condition to be observed by Gaia. A similar number of bulges could also be observed although the low luminosity bulges should escape detection. About two thirds of the more distant objects up to 600 Mpc could also be detected increasing the total sample to half a million objects including ellipticals and bulges. The angular size of the detected objects will never exceed 4.72 arcsec which is the size of the largest transmitted windows. An heterogeneous population of elliptical galaxies and bulges will be observable by Gaia. This nearby Universe sample of galaxies should constitute a very rich and interesting sample to study their structural properties and their distribution.

We present a novel approach, based on robust principal components analysis (RPCA) and maximal information coefficient (MIC), to study the redshift dependence of halo baryonic properties. Our data are composed of a set of different physical quantities for primordial minihaloes: dark-matter mass ($M_{\mathrm{dm}}$), gas mass ($M_{\mathrm{gas}}$), stellar mass ($M_{\mathrm{star}}$), molecular fraction (${\mathrm{x_{mol}}}$), metallicity ($Z$), star formation rate (SFR) and temperature. We find that $M_{\mathrm{dm}}$ and $M_{\mathrm{gas}}$ are dominant factors for variance, particularly at high redshift. Nonetheless, with the emergence of the first stars and subsequent feedback mechanisms, ${\mathrm{x_{mol}}}$, SFR and $Z$ start to have a more dominant role. Standard PCA gives three principal components (PCs) capable to explain more than 97 per cent of the data variance at any redshift (two PCs usually accounting for no less than 92 per cent), whilst the first PC from the RPCA analysis explains no less than 84 per cent of the total variance in the entire redshift range (with two PCs explaining $ \gtrsim 95$ per cent anytime). Our analysis also suggests that all the gaseous properties have a stronger correlation with $M_{\mathrm{gas}}$ than with $M_{\mathrm{dm}}$, while $M_{\mathrm{gas}}$ has a deeper correlation with ${\mathrm{x_{mol}}}$ than with $Z$ or SFR. This indicates the crucial role of gas molecular content to initiate star formation and consequent metal pollution from Population III and Population II/I regimes in primordial galaxies. Finally, a comparison between MIC and Spearman correlation coefficient shows that the former is a more reliable indicator when halo properties are weakly correlated.

We present a fully self-consistent simulation of a synthetic survey of the furthermost cosmic explosions. The appearance of the first generation of stars (Population III) in the Universe represents a critical point during cosmic evolution, signaling the end of the dark ages, a period of absence of light sources. Despite their importance, there is no confirmed detection of Population III stars so far. A fraction of these primordial stars are expected to die as pair-instability supernovae (PISNe), and should be bright enough to be observed up to a few hundred million years after the big bang. While the quest for Population III stars continues, detailed theoretical models and computer simulations serve as a testbed for their observability. With the upcoming near-infrared missions, estimates of the feasibility of detecting PISNe are not only timely but imperative. To address this problem, we combine state-of-the-art cosmological and radiative simulations into a complete and self-consistent framework, which includes detailed features of the observational process. We show that a dedicated observational strategy using $\lesssim 8$ per cent of total allocation time of the James Webb Space Telescope mission can provide us up to $\sim 9-15$ detectable PISNe per year.

This paper consists in analyzing an action that describes boson and fermion fields minimally coupled to the gravity and a common matter field. The self-interaction potentials of the fields are not chosen a priori but from the Noether symmetry approach. The Noether forms of the potentials allow the boson field to play the role of dark energy and matter and the fermion field to behave as standard matter. The constant of motion and the cyclic variable associated with the Noether symmetry allow the complete integration of the field equations, whose solution produces a Universe with alternated periods of accelerated and decelerated expansion.

In this work it is investigated general models with interactions between two canonical scalar fields and between one non-canonical (tachyon-type) and one canonical scalar field. The potentials and couplings to the gravity are selected through the Noether symmetry approach. These general models are employed to describe interactions between dark energy and dark matter, with the fields being constrained by the astronomical data. The cosmological solutions of some cases are compared with the observed evolution of the late Universe.

Recently, it has been suggested that a critical electrical field arises during the gravitational collapse of massive stars leading to a vacuum polarization. This, in turn, leads to the necessity of a reexamination of the gravito-electrodynamical properties of compact stars of the class of neutron stars. Rotondo, Rueda, Ruffini and Xue claim to have proved the impossibility of local charge neutrality and then solved the coupled system of the general relativistic Thomas-Fermi-Einstein-Maxwell equations for the structure of neutron stars. Within the same approach of Avellar and Horvath (2012) we have calculated how the global neutrality hypothesis affects the order/disorder of these systems for a simple equation of state. We show the relative preference of local vs. global conservation in terms of the obtained information content of the systems under consideration.

The near-infrared spectral region is becoming a very useful wavelength range to detect and quantify the stellar population of galaxies. Models are developing to predict the contribution of TP-AGB stars, that should dominate the NIR spectra of populations 0.3 to 2 Gyr old. When present in a given stellar population, these stars leave unique signatures that can be used to detect them unambiguously. However, these models have to be tested in a homogeneous database of star-forming galaxies, to check if the results are consistent with what is found from different wavelength ranges. In this work we performed stellar population synthesis on the nuclear and extended regions of 23 star-forming galaxies to understand how the star-formation tracers in the near-infrared can be used in practice. The stellar population synthesis shows that for the galaxies with strong emission in the NIR, there is an important fraction of young/intermediate population contributing to the spectra, which is probably the ionisation source in these galaxies. Galaxies that had no emission lines measured in the NIR were found to have older average ages and less contribution of young populations. Although the stellar population synthesis method proved to be very effective to find the young ionising population in these galaxies, no clear correlation between these results and the NIR spectral indexes were found. Thus, we believe that, in practice, the use of these indexes is still very limited due to observational limitations.

The ESA Gaia mission, to be launched during 2013, will observe billions of objects, among which many galaxies, during its scanning of the sky. This will provide a large space-based dataset with unprecedented spatial resolution. Because of its natural Galactic and Astrometric priority, Gaia's observational strategy was optimized for point sources. Nonetheless, it is expected that 10^6 sources will be extragalactic, and a large portion of them will be angularly small galaxies. Although the mission was designed for point sources, a dedicated analysis of the raw data will allow the recovery of morphology of those objects at a 0.2" level. This may constitute a unique all-sky survey of such galaxies. We describe the conceptual design of the method we created for performing the morphological analysis of these objects as well as first results obtained from data simulations of low-resolution, highly binned, satellite data. Based on the obtained results we conclude that it is possible to push the limits of the Gaia space mission by analyzing galaxy morphology. (Abridged)

Structures in Warm Dark Matter (WDM) models are exponentially suppressed below a certain scale, characterized by the dark matter particle mass, $m_{\rm x}$. Since structures form hierarchically, the presence of collapsed objects at high-redshifts can set strong lower limits on $m_{\rm x}$. We place robust constraints on $m_{\rm x}$ using recent results from the {\it Swift} database of high-redshift gamma-ray bursts (GRBs). We parameterize the redshift evolution of the ratio between the cosmic GRB rate and star formation rate (SFR) as $\propto (1+z)^\alpha$, thereby allowing astrophysical uncertainties to partially mimic the cosmological suppression of structures in WDM models. Using a maximum likelihood estimator on two different $z>4$ GRB subsamples (including two bursts at $z>8$), we constrain $m_{\rm x} \gtrsim 1.6$-1.8 keV at 95% C.L., when marginalized over a flat prior in $\alpha$. We further estimate that 5 years of a SVOM-like mission would tighten these constraints to $m_{\rm x} \gtrsim 2.3 $ keV. Our results show that GRBs are a powerful probe of high-redshift structures, providing robust and competitive constraints on $m_{\rm x}$.

Ground- and space-based observations of solar flares from radio wavelengths to gamma-rays have produced considerable insights but raised several unsolved controversies. The last unexplored wavelength frontier for solar flares is in the range of submillimeter and infrared wavelengths. Here we report the detection of an intense impulsive burst at 30 THz using a new imaging system. The 30 THz emission exhibited remarkable time coincidence with peaks observed at microwave, mm/submm, visible, EUV and hard X-ray wavelengths. The emission location coincides with a very weak white-light feature, and is consistent with heating below the temperature minimum in the atmosphere. However, there are problems in attributing the heating to accelerated electrons. The peak 30 THz flux is several times larger than the usual microwave peak near 9 GHz, attributed to non-thermal electrons in the corona. The 30 THz emission could be consistent with an optically thick spectrum increasing from low to high frequencies. It might be part of the same spectral component found at sub-THz frequencies whose nature remains mysterious. Further observations at these wavelengths will provide a new window for flare studies.

Recent models show that TP-AGB stars should dominate the NIR spectra of populations 0.3 to 2 Gyr old, leaving unique signatures that can be used to detect young/intermediate stellar population in galaxies. However, no homogeneous database of star-forming galaxies is available in the NIR to fully explore these results. With this in mind, we study the NIR spectra of a sample of 23 HII and starburst galaxies, aimed at characterizing the most prominent spectral features and continuum shape in the 0.8-2.4 micron region of these objects. Spectral indices are derived for the relevant absorption lines/bands and a comparison with optical indices of the same sample available in the literature is made. We found no correlation between the optical and the NIR indexes. This is probably due to the differences in aperture between these two sets of data. That result is further supported by the absence or weakness of emission lines in the NIR for a subsample galaxies, while in the optical the emission lines are strong and clear, which means that the ionisation source in many of these galaxies is not nuclear, but circumnuclear or located in hot spots. We detected important signatures predicted for a stellar population dominated by the TP-AGBs, like CN 1.1 micron and CO 2.3 micron. In at least one galaxy (NGC 4102) the CN band at 1.4 micron was detected for the first time. We also detect TiO and ZrO bands that have never been reported before in extragalactic sources. The shape of the continuum emission is found to be strongly correlated to the presence/lack of emission lines. An observational template for the star-forming galaxies is derived to be used as a benchmark of stellar population(s) in starburst galaxies against which to compare near-IR spectroscopy of different types of galaxies, especially those with AGN activity and/or those at high-redshift.

In this work, we propose the use of Kernel Principal Component Analysis (KPCA) combined with k = 1 nearest neighbour algorithm (1NN) as a framework for supernovae (SNe) photometric classification. The classification is entirely based on information within the spectroscopic confirmed sample and each new light curve is classified one at a time. This allows us to update the principal component (PC) parameter space if a new spectroscopic light curve is available while also avoids the need of re-determining it for each individual new classification. We applied the method to different instances of the \textit{Supernova Photometric Classification Challenge} (SNPCC) data set. Our method provide good purity results in all data sample analysed, when SNR$\geq$5. As a consequence, we can state that if a sample as the post-SNPCC was available today, we would be able to classify $\approx 15%$ of the initial data set with purity $\gtrsim$ 90% (D$_{7}$+SNR3). Results from the original SNPCC sample, reported as a function of redshift, show that our method provides high purity (up to $\approx 97%$), specially in the range of $0.2\leq z < 0.4$, when compared to results from the SNPCC, while maintaining a moderate figure of merit ($\approx 0.25$). We also present results for SNe photometric classification using only pre-maximum epochs, obtaining 63% purity and 77% successful classification rates (SNR$\geq$5). Results are sensitive to the information contained in each light curve, as a consequence, higher quality data points lead to higher successful classification rates. The method is flexible enough to be applied to other astrophysical transients, as long as a training and a test sample are provided.

We study the statistical properties (such as shape and spin) of high-z halos likely hosting the first (PopIII) stars with cosmological simulations including detailed gas physics. In the redshift range considered ($11 < z < 16$) the average sphericity is $<s> = 0.3 \pm 0.1$, and for more than 90% of halos the triaxiality parameter is $T \lesssim 0.4$, showing a clear preference for oblateness over prolateness. Larger halos in the simulation tend to be both more spherical and prolate: we find $s \propto M_h^{\alpha_s}$ and $T \propto M_h^{\alpha_T}$, with $\alpha_s \approx 0.128$ and $\alpha_T= 0.276$ at z = 11. The spin distributions of dark matter and gas are considerably different at $z=16$, with the baryons rotating slower than the dark matter. At lower redshift, instead, the spin distributions of dark matter and gas track each other almost perfectly, as a consequence of a longer time interval available for momentum redistribution between the two components. The spin of both the gas and dark matter follows a lognormal distribution, with a mean value at z=16 of $<\lambda> =0.0184$, virtually independent of halo mass. This is in good agreement with previous studies. Using the results of two feedback models (MT1 and MT2) by McKee & Tan (2008) and mapping our halo spin distribution into a PopIII IMF, we find that at high-$z$ the IMF closely tracks the spin lognormal distribution. Depending on the feedback model, though, the distribution can be centered at $\approx 65 M_\odot$ (MT1) or $\approx 140 M_\odot$ (MT2). At later times, model MT1 evolves into a bimodal distribution with a second prominent peak located at $35-40 M_\odot$ as a result of the non-linear relation between rotation and halo mass. We conclude that the dark matter halo properties might be a key factor shaping the IMF of the first stars.

We construct a theoretical model to predict the number of orphan afterglows (OA) from gamma-ray bursts (GRBs) triggered by primordial metal-free (Pop III) stars expected to be observed by the Gaia mission. In particular, we consider primordial metal-free stars that were affected by radiation from other stars (Pop III.2) as a possible target. We use a semi-analytical approach that includes all relevant feedback effects to construct cosmic star formation history and its connection with the cumulative number of GRBs. The OA events are generated using the Monte Carlo method, and realistic simulations of Gaia's scanning law are performed to derive the observation probability expectation. We show that Gaia can observe up to 2.28 $\pm$ 0.88 off-axis afterglows and 2.78 $\pm$ 1.41 on-axis during the five-year nominal mission. This implies that a nonnegligible percentage of afterglows that may be observed by Gaia ($\sim 10%$) could have Pop III stars as progenitors.

We discuss the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive primordial stars. We construct a theoretical model to calculate the rate and detectability of these GRBs taking into account all important feedback and recent results from numerical simulations of pristine gas. We expect to observe a maximum of N $\lesssim$ 0.2 GRBs per year integrated over at z > 6 with \textit{Swift} and N $\lesssim$ 10 GRBs per year integrated over at z > 6 with EXIST.

Based on astrophysical constraints derived from Chandrasekhar's mass limit for white-dwarfs, we study the effects of the model on the parameters of unparticle-inspired gravity, on scales $\Lambda_U > 1 \; TeV$ and $d_U \approx 1$.

The motion of a point like object of mass M passing through the background potential of massive collisionless particles (m << M) suffers a steady deceleration named dynamic friction. In his classical work, Chandrasekhar assumed a Maxwellian velocity distribution in the halo and neglected the self gravity of the wake induced by the gravitational focusing of the mass M. In this paper, by relaxing the validity of the Maxwellian distribution due to the presence of long range forces, we derive an analytical formula for the dynamic friction in the context of the q-nonextensive kinetic theory. In the extensive limiting case (q = 1), the classical Gaussian Chandrasekhar result is recovered. As an application, the dynamic friction timescale for Globular Clusters spiraling to the galactic center is explicitly obtained. Our results suggest that the problem concerning the large timescale as derived by numerical N-body simulations or semi-analytical models can be understood as a departure from the standard extensive Maxwellian regime as measured by the Tsallis nonextensive q-parameter.

The origin of magnetic fields in astrophysical objects is a challenging problem in astrophysics. Throughout the years, many scientists have suggested that non-minimal gravitational-electromagnetic coupling (NMGEC) could be the origin of the ubiquitous astrophysical magnetic fields. We investigate the possible origin of intense magnetic fields $\sim 10^{15}-10^{16}$ by NMGEC near rotating neutron stars and black holes, connected with magnetars, quasars, and gamma-ray bursts. Whereas these intense magnetic fields are difficult to explain astrophysically, we find that they are easily explained by NMGEC.

Radio and optical observations of the evolution of flare-associated phenomena have shown an initial and rapid burst at 0.4 THz only followed subsequently by a localized chromospheric heating producing an H{\alpha} brightening with later heating of the whole active region. A major instability occurred several minutes later producing one impulsive burst at microwaves only, associated with an M2.0 GOES X-ray flare that exhibited the main H{\alpha} brightening at the same site as the first flash. The possible association between long-enduring time profiles at soft X-rays, microwaves, H{\alpha} and sub-THz wavelengths is discussed. In the decay phase the H{\alpha} movie shows a disrupting magnetic arch structure ejecting dark, presumably chromospheric, material upwards. The time sequence of events suggests genuine interdependent and possibly non-thermal instabilities triggering phenomena, with concurrent active region plasma heating and material ejection.

We calculate the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive first-generation (Population III; Pop III) stars. The Pop III GRBs could be super-energetic with the isotropic energy up to $E_{\rm iso} \gtrsim 10^{55-57}$ ergs, providing a unique probe of the high-redshift Universe. We consider both the so-called Pop III.1 stars (primordial) and Pop III.2 stars (primordial but affected by radiation from other stars). We employ a semi-analytical approach that considers inhomogeneous hydrogen reionization and chemical evolution of the intergalactic medium. We show that Pop III.2 GRBs occur more than 100 times more frequently than Pop III.1 GRBs, and thus should be suitable targets for future GRB missions. Interestingly, our optimistic model predicts an event rate that is already constrained by the current radio transient searches. We expect $\sim10-10^4$ radio afterglows above $\sim$ 0.3 mJy on the sky with $\sim 1$ year variability and mostly without GRBs (orphans), which are detectable by ALMA, EVLA, LOFAR, and SKA, while we expect to observe maximum of $N < 20$ GRBs per year integrated over at $z>6$ for Pop III.2 and $N < 0.08$ per year integrated over at $z>10$ for Pop III.1 with EXIST, and $N < 0.2$ for Pop III.2 GRBs per year integrated over at $z > 6$ with \textit{Swift}.

We propose a novel approach, based on Principal Components Analysis, to the use of Gamma-Ray Bursts (GRBs) as probes of cosmic star formation history (SFH) up to very high redshifts. The main advantage of such approach is to avoid the necessity of assuming an \textit{ad hoc} parameterization of the SFH. We first validate the method by reconstructing a known SFH from Monte Carlo-generated mock data. We then apply the method to the most recent \textit{Swift} data of GRBs with known redshift and compare it against the SFH obtained by independent methods. The main conclusion is that the level of star formation activity at $z \approx 9.4$ could have been already as high as the present-day one ($\approx 0.01 M_\odot$ yr$^{-1}$ Mpc$^{-3}$). This is a factor 3-5 times higher than deduced from high-$z$ galaxy searches through drop-out techniques. If true, this might alleviate the long-standing problem of a photon-starving reionization; it might also indicate that galaxies accounting for most of the star formation activity at high redshift go undetected by even the most deep searches.

Some observations of galaxies, and in particular dwarf galaxies, indicate a presence of cored density profiles in apparent contradiction with cusp profiles predicted by dark matter N-body simulations. We constructed an analytical model, using particle distribution functions (DFs), to show how a supernova (SN) explosion can transform a cusp density profile in a small-mass dark matter halo into a cored one. Considering the fact that a SN efficiently removes matter from the centre of the first haloes, we study the effect of mass removal through a SN perturbation in the DFs. We found that the transformation from a cusp into a cored profile is present even for changes as small as 0.5% of the total energy of the halo, that can be produced by the expulsion of matter caused by a single SN explosion.

A model-independent reconstruction of the cosmic expansion rate is essential to a robust analysis of cosmological observations. Our goal is to demonstrate that current data are able to provide reasonable constraints on the behavior of the Hubble parameter with redshift, independently of any cosmological model or underlying gravity theory. Using type Ia supernova data, we show that it is possible to analytically calculate the Fisher matrix components in a Hubble parameter analysis without assumptions about the energy content of the Universe. We used a principal component analysis to reconstruct the Hubble parameter as a linear combination of the Fisher matrix eigenvectors (principal components). To suppress the bias introduced by the high redshift behavior of the components, we considered the value of the Hubble parameter at high redshift as a free parameter. We first tested our procedure using a mock sample of type Ia supernova observations, we then applied it to the real data compiled by the Sloan Digital Sky Survey (SDSS) group. In the mock sample analysis, we demonstrate that it is possible to drastically suppress the bias introduced by the high redshift behavior of the principal components. Applying our procedure to the real data, we show that it allows us to determine the behavior of the Hubble parameter with reasonable uncertainty, without introducing any ad-hoc parameterizations. Beyond that, our reconstruction agrees with completely independent measurements of the Hubble parameter obtained from red-envelope galaxies.

An analytical solution for the discrepancy between observed core-like profiles and predicted cusp profiles in dark matter halos is studied. We calculate the distribution function for Navarro-Frenk-White halos and extract energy from the distribution, taking into account the effects of baryonic physics processes. We show with a simple argument that we can reproduce the evolution of a cusp to a flat density profile by a decrease of the initial potential energy.

We recently predicted the existence of random primordial magnetic fields (RPMF) in the form of randomly oriented cells with dipole-like structure with a cell size $L_0$ and an average magnetic field $B_0$. Here we investigate models for primordial magnetic field with a similar web-like structure, and other geometries, differing perhaps in $L_0$ and $B_0$. The effect of RPMF on the formation of the first galaxies is investigated. The filtering mass, $M_F$, is the halo mass below which baryon accretion is severely depressed. We show that these RPMF could influence the formation of galaxies by altering the filtering mass and the baryon gas fraction of a halo, $f_g$. The effect is particularly strong in small galaxies. We find, for example, for a comoving $B_0=0.1\muG$, and a reionization epoch that starts at $z_s=11$ and ends at $z_e=8$, for $L_0=100\,\text{pc}$ at $z=12$, the $f_g$ becomes severely depressed for $M<10^7\msun$, whereas for $B_0=0$ the $f_g$ becomes severely depressed only for much smaller masses, $M<10^5\msun$. We suggest that the observation of $M_F$ and $f_g$ at high redshifts can give information on the intensity and structure of primordial magnetic fields.

We performed stellar population synthesis on the nuclear and extended regions of NGC 1068 by means of near-infrared spectroscopy to disentangle their spectral energy distribution components. This is the first time that such a technique is applied to the whole 0.8 - 2.4 micron wavelength interval in this galaxy. NGC 1068 is one of the nearest and probably the most studied Seyfert 2 galaxy, becoming an excellent laboratory to study the interaction between black holes, the jets that they can produce and the medium in which they propagate. Our main result is that traces of young stellar population are found at ~ 100 south of the nucleus. The contribution of a power-law continuum in the centre is about 25%, which is expected if the light is scattered from a Seyfert 1 nucleus. We find peaks in the contribution of the featureless continuum about 100 - 150 pc from the nucleus on both sides. They might be associated with regions where the jet encounters dense clouds. Further support to this scenario is given by the peaks of hot dust distribution found around these same regions and the H2 emission line profile, leading us to propose that the peaks might be associate to regions where stars are being formed. Hot dust also has an important contribution to the nuclear region, reinforcing the idea of the presence of a dense, circumnuclear torus in this galaxy. Cold dust appears mostly in the south direction, which supports the view that the southwest emission is behind the plane of the galaxy and is extinguished very likely by dust in the plane. Intermediate age stellar population contributes significantly to the continuum, specially in the inner 200 pc.

We report the first simultaneous zJHK spectroscopy on the archetypical Seyfert 2 Galaxy NGC 1068 covering the wavelength region 0.9 to 2.4 micron. The slit, aligned in the NS direction and centred in the optical nucleus, maps a region 300 pc in radius at sub-arcsec resolution, with a spectral resolving power of 360 km s^-1. This configuration allow us to study the physical properties of the nuclear gas including that of the north side of the ionization cone, map the strong excess of continuum emission in the K-band and attributed to dust and study the variations, both in flux and profile, in the emission lines. Our results show that (1) Mid- to low-ionization emission lines are splitted into two components, whose relative strengths vary with the position along the slit and seem to be correlated with the jet. (2) The coronal lines are single-peaked and are detected only in the central few hundred of parsecs from the nucleus. (3) The absorption lines indicate the presence of intermediate age stellar population, which might be a significant contributor to the continuum in the NIR spectra. (4) Through some simple photoionization models we find photoionization as the main mechanism powering the emitting gas. (5) Calculations using stellar features point to a mass concentration inside the 100 - 200 pc of about 10^10 solar masses.

It is generally assumed that the magnetic fields of millisecond pulsars (MSPs) are $\sim 10^{8}$G. We argue that this may not be true and the fields may be appreciably greater. We present six evidences for this: (1) The $\sim 10^{8}$ G field estimate is based on magnetic dipole emission losses which is shown to be questionable; (2) The MSPs in low mass X-ray binaries (LMXBs) are claimed to have $< 10^{11}$ G on the basis of a Rayleygh-Taylor instability accretion argument. We show that the accretion argument is questionable and the upper limit $10^{11}$ G may be much higher; (3) Low magnetic field neutron stars have difficulty being produced in LMXBs; (4) MSPs may still be accreting indicating a much higher magnetic field; (5) The data that predict $\sim 10^{8}$ G for MSPs also predict ages on the order of, and greater than, ten billion years, which is much greater than normal pulsars. If the predicted ages are wrong, most likely the predicted $\sim 10^{8}$ G fields of MSPs are wrong; (6) When magnetic fields are measured directly with cyclotron lines in X-ray binaries, fields $\gg 10^{8}$ G are indicated. Other scenarios should be investigated. One such scenario is the following. Over 85% of MSPs are confirmed members of a binary. It is possible that all MSPs are in large separation binaries having magnetic fields $> 10^{8}$ G with their magnetic dipole emission being balanced by low level accretion from their companions.

We investigate the impact of the existence of a primordial magnetic field on the filter mass, characterizing the minimum baryonic mass that can form in dark matter (DM) haloes. For masses below the filter mass, the baryon content of DM haloes are severely depressed. The filter mass is the mass when the baryon to DM mass ratio in a halo is equal to half the baryon to DM ratio of the Universe. The filter mass has previously been used in semianalytic calculations of galaxy formation, without taking into account the possible existence of a primordial magnetic field. We examine here its effect on the filter mass. For homogeneous comoving primordial magnetic fields of $B_0 \sim 1$ or 2 nG and a reionization epoch that starts at a redshift $z_s=11$ and is completed at $z_r=8$, the filter mass is increased at redshift 8, for example, by factors 4.1 and 19.8, respectively. The dependence of the filter mass on the parameters describing the reionization epoch is investigated. Our results are particularly important for the formation of low mass galaxies in the presence of a homogeneous primordial magnetic field. For example, for $B_0\sim 1\nG$ and a reionization epoch of $z_s\sim 11$ and $z_r\sim7$, our results indicate that galaxies of total mass $M\sim5 \times 10^8\msun$ need to form at redshifts $z_F\gtrsim 2.0$, and galaxies of total mass $M\sim10^8\msun$ at redshifts $z_F\gtrsim 7.7$.

Microgauss magnetic fields are observed in all galaxies at low and high redshifts. The origin of these intense magnetic fields is a challenging question in astrophysics. We show here that the natural plasma fluctuations in the primordial universe (assumed to be random), predicted by the Fluctuation-Dissipation-Theorem, predicts $\sim 0.034 \mu G$ fields over $\sim 0.3$ kpc regions in galaxies. If the dipole magnetic fields predicted by the Fluctuation-Dissipation-Theorem are not completely random, microgauss fields over regions $\gtrsim 0.34$ kpc are easily obtained. The model is thus a strong candidate for resolving the problem of the origin of magnetic fields in $\lesssim 10^{9}$ years in high redshift galaxies.

Various authors have suggested that the gamma-ray burst (GRB) central engine is a rapidly rotating, strongly magnetized, $(\sim 10^{15}-10^{16}$ G) compact object. The strong magnetic field can accelerate and collimate the relativistic flow and the rotation of the compact object can be the energy source of the GRB. The major problem in this scenario is the difficulty of finding an astrophysical mechanism for obtaining such intense fields. Whereas, in principle, a neutron star could maintain such strong fields, it is difficult to justify a scenario for their creation. If the compact object is a black hole, the problem is more difficult since, according to general relativity it has "no hair" (i.e., no magnetic field). Schuster, Blackett, Pauli, and others have suggested that a rotating neutral body can create a magnetic field by non-minimal gravitational-electromagnetic coupling (NMGEC). The Schuster-Blackett form of NMGEC was obtained from the Mikhail and Wanas's tetrad theory of gravitation (MW). We call the general theory NMGEC-MW. We investigate here the possible origin of the intense magnetic fields $\sim 10^{15}-10^{16}$ G in GRBs by NMGEC-MW. Whereas these fields are difficult to explain astrophysically, we find that they are easily explained by NMGEC-MW. It not only explains the origin of the $\sim 10^{15}-10^{16}$G fields when the compact object is a neutron star, but also when it is a black hole.

(Abridged) The main purpose of this paper is to consider the contribution of all three non-thermal components to total mass measurements of galaxy clusters: cosmic rays, turbulence and magnetic pressures. To estimate the thermal pressure we used public XMM-\textit{Newton} archival data of 5 Abell clusters. To describe the magnetic pressure, we assume a radial distribution for the magnetic field, $B(r) \propto \rho_{g}^{\alpha}$, to seek generality we assume $\alpha$ within the range of 0.5 to 0.9, as indicated by observations and numerical simulations. For the turbulent component, we assumed an isotropic pressure, $P_{\rm turb} = {1/3}\rho_{\rm g}(\sigma_{r}^{2}+\sigma_{t}^{2})$. We also consider the contribution of cosmic ray pressure, $P_{cr}\propto r^{-0.5}$. It follows that a consistent description for the non-thermal component could yield variation in mass estimates that vary from 10% up to $\sim$30%. We verified that in the inner parts of cool-core clusters the cosmic ray component is comparable to the magnetic pressure, while in non cool-core cluster the cosmic ray component is dominant. For cool-core clusters the magnetic pressure is the dominant component, contributing with more than 50% of total mass variation due to non-thermal pressure components. However, for non cool-core clusters, the major influence comes from the cosmic ray pressure that accounts with more than 80% of total mass variation due to non-thermal pressure effects. For our sample, the maximum influence of the turbulent component to total mass variation can be almost 20%. We show that this analysis can be regarded as a starting point for a more detailed and refined exploration of the influence of non-thermal pressure in the intra-cluster medium (ICM).

We compare the results of the photometrical analysis of barred galaxies with those of a similar analysis from N-body simulations. The photometry is for a sample of nine barred galaxies observed in the J and Ks bands with the CANICA near infrared (NIR) camera at the 2.1-m telescope of the Observatorio Astrofisico Guillermo Haro (OAGH) in Cananea, Sonora, Mexico. The comparison includes radial ellipticity profiles and surface brightness (density for the N-body galaxies) profiles along the bar major and minor axes. We find very good agreement, arguing that the exchange of angular momentum within the galaxy plays a determinant role in the evolution of barred galaxies.

The structure of globular clusters and elliptical galaxies are described in an unified way through a new class of lowered models inspired on the nonextensive kinetic theory. These power law models are specified by a single parameter q which quantifies to what extent they depart from the class of lowered stellar distributions discussed by Michie and King. For q equal to unity, the Michie-King profiles are recovered. However, for q smaller than unity there is a gradual modification in the shape of the density profiles which depends on the degree of tidal damage imposed on the model, thereby also providing a good fit for globular clusters. It is also shown that a subclass of these models, those with a deeper potential and $q$ slightly less than unity, present a distribution resembling the de Vaucoulers $r^{1/4}$ profile which yields a good description of the structure of elliptical galaxies. This subset of models follows this trend, with a slight departure over nearly 10 orders of magnitudes.

(abridged) We suggest here that the large scale fields $\sim \mu$G, observed in galaxies at both high and low redshifts by Faraday Rotation Measurements (FRMs), have their origin in the electromagnetic fluctuations that naturally occurred in the dense hot plasma that existed just after the QHPT. We evolve the predicted fields to the present time. The size of the region containing a coherent magnetic field increased due to the fusion (polymerization) of smaller regions. Magnetic fields (MFs) $\sim 10 \mu G$ over a comoving $\sim 1$ pc region are predicted at redshift z $\sim 10$. These fields are orders of magnitude greater than those predicted in previous scenarios for creating primordial magnetic fields. Line-of-sight average magnetic fields (MFs) $\sim$ $10^{-2}$ $\mu$G, valid for FRMs, are obtained over a 1 Mpc comoving region at the redshift z $\sim$ 10. In the collapse to a galaxy (comoving size $\sim$ 30 kpc) at z $\sim$ 10, the fields are amplified to $\sim 10 \mu$G. This indicates that the MFs created immediately after the QHPT, predicted by the Fluctuation-Dissipation Theorem, could be the origin of the $\sim \mu G$ fields observed by FRMs in galaxies at both high and low redshifts.

We have obtained deep J and Ks images of a sample of nine barred galaxies in order to collect a reliable and homogeneous set of images to which N-body simulations of barred galaxies will be compared. The observations were performed using the new near-infrared camera available at the 2.1-m telescope of the Observatorio Astrofisico Guillermo Haro (OAGH) in Cananea, Sonora, Mexico. We present the results of surface photometry techniques applied to the observed images, as well as to the deprojected images. These results include radial profiles of surface brightness (elliptically averaged), colour, position angle, ellipticity and the b4 Fourier component. In addition, we present isophotal maps, colour maps, surface brightness profiles along the bar major and minor axes, characteristic radial scale-lengths and bar length estimates. We discuss how projection effects can influence these measurements and the uncertainties introduced by deprojecting galaxy images. We show that analytical expressions can be used to obtain reliable estimates of deprojected bar lengths, ellipticities and position angles directly from the observed images. These expressions are based on the assumption that the outer parts of the bar are vertically thin, as shown by theoretical work. The usefulness of our data in addressing issues on bar formation and evolution is also discussed. In particular, we present results showing a steep drop in the ellipticity profile, as expected for bar formation processes in which the dark matter halo plays a fundamental role. Furthermore, we show that the location of this drop is a good indicator of the end of the bar in strongly barred galaxies, as predicted by numerical models.

To better understand our Galaxy, we investigate the pertinency of describing the sys tem of nearby disk stars in terms of a two-components Schwarzschild velocity distributio n.Using the proper motion and parallax information of Hipparcos database, we determine t he parameters characterizing the local stellar velocity field of a sample of 22000 disk stars. The sample we use is essentially the same as the one described by the criteria ad opted to study the LSR and the stream motion of the nearby stellar population

In an effort to obtain further observational evidences for secular evolution processes in galaxies, as well as observational constraints to current theoretical models of secular evolution, we have used BVRI and Ks images of a sample of 18 barred galaxies to measure the lengths and colours of bars, create colour maps and estimate global colour gradients. In addition, applying a method we developed in a previous article, we could distinguish for 7 galaxies in our sample those whose bars have been recently formed from the ones with already evolved bars. We estimated an average difference in the optical colours between young and evolved bars that may be translated to an age difference of the order of 10 Gyr, meaning that bars may be long standing structures. Moreover, our results show that, on average, evolved bars are longer than young bars. This seems to indicate that, during its evolution, a bar grows longer by capturing stars from the disk, in agreement with recent numerical and analytical results.

Along with a brief analysis we present data obtained from BVRI and Ks images of a sample of 19 galaxies (18 barred and 1 unbarred) which will be further explored in a future paper. We measured the lengths and colors of the bars, created color maps and estimated global color gradients. Applying a method developed in a companion paper, we could distinguish for 7 galaxies in our sample those whose bars have been recently formed from the ones with already evolved bars. We estimated an average difference in the optical colors between young and evolved bars that may be translated to an age difference of the order of 10 Gyr, meaning that bars may be, at least in some cases, long standing structures. Moreover, our results show that, on average, evolved bars are longer than young bars. This seems to indicate that, during its evolution, a bar grows longer by capturing stars from the disk, in agreement with recent numerical and analytical results. Although the statistical significance of these results is low, and further studies are needed to confirm them, we discuss the implications from our results on the possibility of bars being a recurrent phenomenon. We also present isophotal contours for all our images as well as radial profiles of relevant photometric and geometric parameters.

The last few years witnessed a change of perspective in the context of galaxy formation and evolution. The major role played by bars in the history of these fundamental building blocks of our universe was finally realized. However, and despite the fact that one of the major concerns of any physical science is to measure timescales for natural phenomena, we have as yet no established way to estimate the ages of these important galactic components. In this contribution we outline the method we developed to estimate bar ages and briefly describe its first results.

In order to perform a detailed study of the stellar kinematics in the vertical axis of bars, we obtained high signal-to-noise spectra along the major and minor axes of the bars in a sample of 14 face-on galaxies, and used them to determine the line of sight stellar velocity distribution, parameterized as Gauss-Hermite series. With these data, we developed a diagnostic tool that allows one to distinguish between recently formed and evolved bars, as well as estimate their ages, assuming that bars form in vertically thin disks, recognizable by low values for the vertical velocity dispersion sigma_z. Through N-body realizations of bar unstable disk galaxies we could also check the time scales involved in the processes which give bars an important vertical structure. We show that sigma_z in evolved bars is roughly around 100 Km/s, which translates to a height scale of about 1.4 Kpc, giving support to scenarios in which bulges form through disk material. Furthermore, the bars in our numerical simulations have values for sigma_z generally smaller than 50 Km/s even after evolving for 2 Gyr, suggesting that a slow process is responsible for making bars as vertically thick as we observe. We verify theoretically that the Spitzer-Schwarzschild mechanism is quantitatively able to explain these observations if we assume that giant molecular clouds are twice as much concentrated along the bar as in the remaining of the disk.

The density profiles and other quantities of physical interest for spherically symmetric systems are computed by assuming that a collisionless stellar gas may relax to the non-Gaussian power law distribution suggested by the nonextensive kinetic theory. There are two different classes of solutions. The first class behaves like a subset of the polytropic Lane-Emden spheres, whereas the second one corresponds to a transition between two different polytropic indices. Unlike the isothermal Maxwellian sphere, the total mass and sizes of both classes are finite for a large range of the nonextensive q-parameter.

In an effort to elevate to higher grounds our understanding on the impact of the formation and evolution of bars in the formation and evolution of galaxies, we have developed a diagnostic tool to distinguish between recently formed and evolved bars. Our method was applied in the study of a sample of 14 galaxies and revealed that, apparently, AGN activity tends to appear in galaxies which have young bars rather than evolved bars. This suggests that the time scale for the fueling of AGN by bars is short, and may help to explain, for instance, why there is not a clear correlation between the presence of bars and AGN in galaxies.

We present BUDDA (BUlge/Disk Decomposition Analysis), a new code devoted to perform a two-dimensional bulge/disk decomposition directly from the images of galaxies. The bulge component is fitted with a generalized Sersic profile whereas disks have an exponential profile. No other components are included. Bars and other sub-structures, like lenses, rings, inner bars and inner disks, are studied with the residual images obtained through the subtraction of bulges and disks from the original images. This means that a detailed structural analysis of galaxies may be performed with a small number of parameters, and sub-structures may be directly studied with no a priori assumptions. To illustrate its usefulness, we present the results obtained with a sample of 51 mostly early-type galaxies (but covering the whole Hubble sequence). These results show some of the applications in which the code may be used: the determination of parameters for Fundamental Plane and structural studies, quantitative morphological classification of galaxies and identification and study of hidden sub-structures. We have determined the structural parameters of the galaxies in our sample and found many examples of hidden inner disks in ellipticals, secondary bars, nuclear rings and dust lanes in lenticulars and spirals, and also wrong morphological classification cases.

We have performed for the first time in the literature a 2D structural analysis on the barred lenticular and face--on galaxies NGC 4608 and NGC 5701. The results indicate that these galaxies either have never had large disks or their disks were almost completely destroyed by their strong bars, due to secular evolution processes. We discuss these surprising conclusions checking for signs of secular evolution, considering bar forming instabilities, and suggesting, based on N-body simulations, a new mechanism to form bars in spheroids, which includes non-spherical halos. Quantitative predictions from our new mechanism are compared with those from other recent models for bar formation and evolution.

To further enhance our understanding on the formation and evolution of bars in lenticular (S0) galaxies, we are undertaking a detailed photometric and spectroscopic study on a sample of 22 objects. Here we report the results of a 2D structural analysis on two barred face-on S0's, which indicate that presently these galaxies do not possess disks. We discuss two possibilities to explain these surprising results, namely strong secular evolution and bar formation without disks.

We present optical and IR spectroscopy of the first two months of evolution of the Type II SN 1999em. We combine these data with high-quality optical/IR photometry beginning only three days after shock breakout, in order to study the performance of the ``Expanding Photosphere Method'' (EPM) in the determination of distances. With this purpose we develop a technique to measure accurate photospheric velocities by cross-correlating observed and model spectra. The application of this technique to SN 1999em shows that we can reach an average uncertainty of 11% in velocity from an individual spectrum. Our analysis shows that EPM is quite robust to the effects of dust. In particular, the distances derived from the VI filters change by only 7% when the adopted visual extinction in the host galaxy is varied by 0.45 mag. The superb time sampling of the BVIZJHK light-curves of SN 1999em permits us to study the internal consistency of EPM and test the dilution factors computed from atmosphere models for Type II plateau supernovae. We find that, in the first week since explosion, the EPM distances are up to 50% lower than the average, possibly due the presence of circumstellar material. Over the following 65 days, on the other hand, our tests lend strong credence to the atmosphere models, and confirm previous claims that EPM can produce consistent distances without having to craft specific models to each supernova. This is particularly true for the VI filters which yield distances with an internal consistency of 4%. From the whole set of BVIZJHK photometry, we obtain an average distance of 7.5+/-0.5 Mpc, where the quoted uncertainty (7%) is a conservative estimate of the internal precision of the method obtained from the analysis of the first 70 days of the supernova evolution.

We present an analysis of the orbital properties in 9 intermediate-redshifts cluster of the CNOC1 survey and we compare them to a control sample of 12 nearby clusters. Similar to the nearby elliptical galaxies, the bulge-dominated galaxies in clusters at redshifts ~0.1-0.4 present orbits that are more eccentric than those for disk-dominated galaxies. However, the orbital segregation is less significant than that found for elliptical and spiral galaxies in nearby cluster. When galaxies are separated by colors - red galaxies with colors in the rest frame (U-V)_o > 1.4, and blue galaxies with (U-V)_o =< 1.4 - the strongest orbital segregation is found. Therefore, the segregation we found seems to modify more efficiently the star formation activity than the internal shape of the galaxies. When we compare the orbits of early-type galaxies at intermediate-redshift with those for z=0, they seem to develop significant changes getting much more eccentric. A different behavior is observed in the late-type galaxies, which present no-significant evolution in their orbit shapes.

A kinematical analysis applied to a sample of galaxy clusters indicates that the differences between the velocity distribution of elliptical and spiral galaxies are associated with the shape of their orbit families. The orbital anisotropies present on each morphological population could be measured with the use of a parameter which is the ratio of the radial and tangential velocity dispersions, and can be recovered through the observed line-of-sight velocity distribution. When a Gaussian velocity distribution is assumed, having different dispersions along the radial and tangential directions, we conclude that the orbits of elliptical galaxies in clusters are close to radial, while spirals have more circular shaped or isotropic orbits. Lenticulars galaxies shares an intermediate orbital parameter, between spirals and ellipticals.

A catalogue of binary galaxies with 621 pairs has been determined by applying a surface density enhancement procedure to {\it The Surface Photometry Catalogue of the ESO-Uppsala Galaxies}. The method does not require any redshift information. An additional restriction, based on objective criteria that take into account the completeness of the source catalogue, led to a sample of 189 isolated pairs which are listed. We have obtained the optical luminosity function of binary galaxies in the catalogue, from which we estimate that the luminosity density of galaxies in binaries is $\sim$4% of that found for field galaxies. The general properties of our sample are similar to those from CPG and CMG.