Habitability has been generally defined as the capability of an environment to support life. Ecologists have been using Habitat Suitability Models (HSMs) for more than four decades to study the habitability of Earth from local to global scales. Astrobiologists have been proposing different habitability models for some time, with little integration and consistency between them and different in function to those used by ecologists. In this white paper, we suggest a mass-energy habitability model as an example of how to adapt and expand the models used by ecologists to the astrobiology field. We propose to implement these models into a NASA Habitability Standard (NHS) to standardize the habitability objectives of planetary missions. These standards will help to compare and characterize potentially habitable environments, prioritize target selections, and study correlations between habitability and biosignatures. Habitability models are the foundation of planetary habitability science. The synergy between the methods used by ecologists and astrobiologists will help to integrate and expand our understanding of the habitability of Earth, the Solar System, and exoplanets.

Rings around giant exoplanets (hereafter 'exorings') are still a missing planetary phenomenon among the vast number of discovered planets. Despite the fact there exist a large number of methods for identifying and characterizing these exorings, none of them has been successful to date. Most of those efforts focus on the photometric signatures produced by rings around transiting exoplanets; thus, little interest has been intended for the detectable signatures that non-transiting ringed planets might cause owing to the excess of scattered starlight from both their atmosphere and the considerably large surface of their (hypothetical) rings-system. This extra scattering produced by exorings would occur at an orbital location defined here as the `the summer solstice' of a stellar light curve. In this letter, we develop a first-order model to estimate the photometric signatures of non-transiting exorings and predict their detectability by using present and future facilities. We also show how, besides the discovery itself, our model can be used to constrain orbital and physical parameters of planet-rings systems.

We present the results from four stellar occultations by (486958) Arrokoth, the flyby target of the New Horizons extended mission. Three of the four efforts led to positive detections of the body, and all constrained the presence of rings and other debris, finding none. Twenty-five mobile stations were deployed for 2017 June 3 and augmented by fixed telescopes. There were no positive detections from this effort. The event on 2017 July 10 was observed by SOFIA with one very short chord. Twenty-four deployed stations on 2017 July 17 resulted in five chords that clearly showed a complicated shape consistent with a contact binary with rough dimensions of 20 by 30 km for the overall outline. A visible albedo of 10% was derived from these data. Twenty-two systems were deployed for the fourth event on 2018 Aug 4 and resulted in two chords. The combination of the occultation data and the flyby results provides a significant refinement of the rotation period, now estimated to be 15.9380 $\pm$ 0.0005 hours. The occultation data also provided high-precision astrometric constraints on the position of the object that were crucial for supporting the navigation for the New Horizons flyby. This work demonstrates an effective method for obtaining detailed size and shape information and probing for rings and dust on distant Kuiper Belt objects as well as being an important source of positional data that can aid in spacecraft navigation that is particularly useful for small and distant bodies.

Close-in giant planets represent the most significant evidence of planetary migration. If large exomoons form around migrating giant planets which are more stable (e.g. those in the Solar System), what happens to these moons after migration is still under intense research. This paper explores the scenario where large regular exomoons escape after tidal-interchange of angular momentum with its parent planet, becoming small planets by themselves. We name this hypothetical type of object a \textit{ploonet}. By performing semi-analytical simulations of tidal interactions between a large moon with a close-in giant, and integrating numerically their orbits for several Myr, we found that in $\sim$50 per cent of the cases a young ploonet may survive ejection from the planetary system, or collision with its parent planet and host star, being in principle detectable. Volatile-rich ploonets are dramatically affected by stellar radiation during both planetocentric and siderocentric orbital evolution, and their radius and mass change significantly due to the sublimation of most of their material during time-scales of hundred of Myr. We estimate the photometric signatures that ploonets may produce if they transit the star during the phase of evaporation, and compare them with noisy lightcurves of known objects ("Kronian" stars and non-periodical dips in dusty lightcurves). Additionally, the typical transit timing variations (TTV) induced by the interaction of a ploonet with its planet are computed. We find that present and future photometric surveys' capabilities can detect these effects and distinguish them from those produced by other nearby planetary encounters.

Exoplanet discoveries have motivated numerous efforts to find unseen populations of exomoons, yet they have been unsuccessful. A plausible explanation is that most discovered planets are located on close-in orbits, which would make their moons prone to tidal evolution and orbital detachment. In recent models of tidally-driven migration of exomoons, evolving planets might prevent what was considered their most plausible fate (i.e. colliding against their host planet), favouring scenarios where moons are pushed away and reach what we define as the "satellite tidal orbital parking" distance ($a_\mathrm{stop}$), which is often within the critical limit for unstable orbits and depends mainly on the system's initial conditions: mass-ratio, semi-major axes, and rotational rates. By using semi-analytical calculations and numerical simulations, we calculate $a_\mathrm{stop}$ for different initial system parameters and constrain the transit detectability of exomoons around close-in planets. We found that systems with $M_\mathrm{m}/M_\mathrm{p} \geq 10^{-4}$, which are less likely to form, are also stable and detectable with present facilities (e.g. Kepler and TESS) through their direct and secondary effects in planet+moon transit, as they are massive, oversized, and migrate slowly. In contrast, systems with lower moon-to-planet mass ratios are "ephemeral" and hardly detectable. Moreover, any detection, confirmation, and full characterisation would require both the short cadence capabilities of TESS and high photometric sensitivity of ground-based observatories. Finally, despite the shortage of discovered long-period planets in currently available databases, the tidal migration model adopted in this work supports the idea that they are more likely to host the first detectable exomoon.

During lunar eclipse of January 21, 2019 a meteoroid impacted the Moon producing a visible light flash. The impact was witnessed by casual observers offering an opportunity to study the phenomenon from multiple geographical locations. We use images and videos collected by observers in 7 countries to estimate the location, impact parameters (speed and incoming direction) and energy of the meteoroid. Using parallax, we achieve determining the impact location at lat. $-29.43^{+0.30}_{-0.21}$, lon. $-67.89^{+0.07}_{-0.09}$ and geocentric distance as 356553 km. After devising and applying a photo-metric procedure for measuring flash standard magnitudes in multiple RGB images having different exposure times, we found that the flash, had an average G-magnitude $\langle G\rangle = 6.7\pm0.3$. We use gravitational ray tracing (GRT) to estimate the orbital properties and likely radiant of the impactor. We find that the meteoroid impacted the moon with a speed of $14^{+7}_{-6}$ km/s (70% C.L.) and at a shallow angle, $\theta < 38.2$ degrees. Assuming a normal error for our estimated flash brightness, educated priors for the luminous efficiency and object density, and using the GRT-computed probability distributions of impact speed and incoming directions, we calculate posterior probability distributions for the kinetic energy (median $K_{\rm med}$ = 0.8 kton), body mass ($M_{\rm med}$ = 27 kg) and diameter ($d_{\rm med}$ = 29 cm), and crater size ($D_{\rm med}$ = 9 m). If our assumptions are correct, the crater left by the impact could be detectable by prospecting lunar probes. These results arose from a timely collaboration between professional and amateur astronomers which highlight the potential importance of citizen science in astronomy.

Despite the success of discovering transiting exoplanets, several recently observed objects (e.g. KIC-8462852, J1407 and PDS-110) exhibit unconventional transit signals, whose appropriate interpretation in terms of a spherical single body has been challenging, if not impossible. In the aforementioned examples the presence of a ring-like structure has been proposed for explaining the unusual data. Thus, in this paper we delve into the dynamics of a tilted exoring disturbed by a third close companion, and the role that the Lidov-Kozai mechanism may have to explain irregular and anomalous transit signals of ringed planets, as well as the ring's early stages. To that end, we performed numerical simulations and semi-analytical calculations to assess the ring's dynamical and morphological properties, and their related transit observables. We found that tilted ringed structures undergo short-term changes in shape and orientation that are manifested as strong variations of transit depth and contact times, even between consecutive eclipses. Any detected anomaly in transit characteristics may lead to a miscalculation of the system's properties (planetary radius, semi-major axis, stellar density and others). Moreover, oscillating ring-like structures may account for the strangeness of some light-curve features in already known and future discovered exoplanets.

Hypothetical exomoons around close-in giant planets may migrate inwards and/or outwards in virtue of the interplay of the star, planet and moon tidal interactions. These processes could be responsible for the disruption of lunar systems, the collision of moons with planets or could provide a mechanism for the formation of exorings. Several models have been developed to determine the fate of exomoons when subject to the tidal effects of their host planet. None of them have taken into account the key role that planetary evolution could play in this process. In this paper we put together numerical models of exomoon tidal-induced orbital evolution, results of planetary evolution and interior structure models, to study the final fate of exomoons around evolving close-in gas giants. We have found that planetary evolution significantly affects not only the time-scale of exomoon migration but also its final fate. Thus, if any change in planetary radius, internal mass distribution and rotation occurs in time-scales lower or comparable to orbital evolution, exomoon may only migrate outwards and prevent tidal disruption or a collision with the planet. If exomoons are discovered in the future around close-in giant planets, our results may contribute to constraint planetary evolution and internal structure models.

Compressible, Riemann S-type ellipsoids can emit gravitational waves (GWs) with a chirp-like behavior (hereafter chirping ellipsoids, CELs). We show that the GW frequency-amplitude evolution of CELs (mass $\sim 1$~M$_\odot$, radius $\sim10^3$~km, polytropic equation of state with index $n\approx 3$) is indistinguishable from that emitted by double white dwarfs (DWDs) and by extreme mass-ratio inspirals (EMRIs) composed of an intermediate-mass (e.g.~$10^3~M_\odot$) black hole and a planet-like (e.g.~$10^{-4}~M_\odot$) companion, in a specific frequency interval within the detector sensitivity band in which the GWs of all these systems are quasi-monochromatic. We estimate that for reasonable astrophysical assumptions, the rates in the local Universe of CELs, DWDs and EMRIs in the mass range considered here, are very similar, posing a detection-degeneracy challenge for space-based GW detectors. The astrophysical implications of this CEL-binary detection degeneracy by space-based GW-detection facilities, are outlined.

Every year, a few metre-sized meteoroids impact the atmosphere of the Earth. Most (if not all) of them are undetectable before the impact. Therefore, predicting where and how they will fall seems to be an impossible task. In this letter we show compelling evidence that we can constrain in advance, the dynamical and geometrical conditions of an impact. For this purpose, we analyse the well-documented case of the Chelyabinsk (Russia) impact and the more recent and smaller Vi\~nales (Cuba) event, whose conditions we estimate and provide here. After using the {\em Gravitational Ray Tracing} algorithm (GRT) to "predict" the impact conditions of the aforementioned events, we find that the speed, incoming direction and (marginally) the orbital elements of the corresponding meteoroids could be constrained in advance, starting only on one hand, with the geographical location and time of the impact, and on the other hand, with the distribution in configuration space of Near Earth Objects (NEOs). Any improvement in our capability to predict or at least to constrain impact properties of medium-sized and large meteoroids, will help us to be better prepared for its potentially damaging effects.

The discovery of planets in binaries is one the most interesting outcomes of planetary research. With the growing number of discoveries has also grown the interest on describing their formation, long-term evolution and potential habitability. In this work we revisit the dynamics of planets in S-type binary systems. For that purpose we develop explicit formulas for the secularized octupolar Hamiltonian, coupled with general relativistic corrections and non-conservative interactions. We implemented those formulas in an open-source package \texttt{SecDev3B}, that can be used to reproduce our results or test improved versions of the models. In order to test it, we study the long-term dynamical evolution of S-type binary planets during the pre-main-sequence phase of stellar evolution. During that phase, stellar radius significantly changes in timescales similar to secular timescales. We hypothesize that when close-encounters between the planet and its host star happens (e.g. via Lidov-Kozai effect), particularities in the secular formalism plus changes in stellar radius may alter significantly the dynamical evolution. We study the well-known binary planet HD 80606b and found that an octupolar expansion of the conservative Hamiltonian is required to properly predict its dynamical evolution. We also apply the dynamical model, enriched with results coming from stellar evolutionary models, to demonstrate that in S-type systems around low-mass stars, with relative high inclinations ($\itot\ge 60^\circ$), moderate eccentricities ($0.2\le e\le 0.4$) and planets located around 1 AU, the evolution of stellar radius during the first few hundreds of Myr, alters significantly the timescales of dynamical evolution.

We present the first systematic exploration of earth tides-seismicity correlation in northwestern South America, with a special emphasis in Colombia. For this purpose, we use a dataset of ~167,000 earthquakes, gathered by the Colombian Seismological Network between 1993 and 2017. Most of the events are intermediate-depth earthquakes from the Bucaramanga seismic nest and the Cauca seismic cluster. For this purpose, we implemented a novel approach for the calculation of tidal phases that considers the relative positions of the Earth-Moon-Sun system at the time of the events. After applying the standard Schuster test to the whole dataset and to several earthquake samples (classified by time, location, magnitude and depth), we found strong correlation anomalies with the diurnal and monthly components of the tide (global log(p) values around -7.0 for the diurnal constituent and -12.1 for the monthly constituent), especially for the intermediate depth events. These anomalies suggest that around 16% of the deep earthquakes in Colombia may be triggered by tides, especially when the monthly phase is between 350$^\circ$-10$^\circ$. We attribute our positive results, which favor the tidal-triggering hypothesis, in contrast to previous negative ones to: 1) the size of our dataset, and 2) the method we used to calculate tidal phases. Anyone willing to reproduce our results or to apply our methodology to custom datasets can use the public information system tQuakes that we developed for this work.

With the advent of more and deeper sky surveys, the discovery of interstellar small objects entering into the Solar System has been finally possible. In October 19, 2017, using observations of the Pan-STARRS survey, a fast moving object, now officially named 1I/2017 U1 ('Oumuamua), was discovered in a heliocentric unbound trajectory suggesting an interstellar origin. Assessing the provenance of interstellar small objects is key for understanding their distribution, spatial density and the processes responsible for their ejection from planetary systems. However, their peculiar trajectories place a limit on the number of observations available to determine a precise orbit. As a result, when its position is propagated $\sim 10^5-10^6$ years backward in time, small errors in orbital elements become large uncertainties in position in the interstellar space. In this paper we present a general method for assigning probabilities to nearby stars of being the parent system of an observed interstellar object. We describe the method in detail and apply it for assessing the origin of 'Oumuamua. A preliminary list of potential progenitors and their corresponding probabilities is provided. In the future, when further information about the object and/or the nearby stars be refined, the probabilities computed with our method can be updated. We provide all the data and codes we developed for this purpose in the form of an open source {\tt C/C++/Python package}, {\bf\tt iWander} which is publicly available at http://github.com/seap-udea/iWander.

Tunguska and Chelyabinsk impact events occurred inside a geographical area of only 3.4\% of the Earth's surface. Although two events hardly constitute a statistically significant demonstration of a geographical pattern of impacts, their spatial coincidence is at least tantalizing. To understand if this concurrence reflects an underlying geographical and/or temporal pattern, we must aim at predicting the spatio-temporal distribution of meteoroid impacts on Earth. For this purpose we designed, implemented and tested a novel numerical technique, the "Gravitational Ray Tracing" (GRT) designed to compute the relative impact probability (RIP) on the surface of any planet. GRT is inspired by the so-called ray-casting techniques used to render realistic images of complex 3D scenes. In this paper we describe the method and the results of testing it at the time of large impact events. Our findings suggest a non-trivial pattern of impact probabilities at any given time on Earth. Locations at $60-90\deg$ from the apex are more prone to impacts, especially at midnight. Counterintuitively, sites close to apex direction have the lowest RIP, while in the antapex RIP are slightly larger than average. We present here preliminary maps of RIP at the time of Tunguska and Chelyabinsk events and found no evidence of a spatial or temporal pattern, suggesting that their coincidence was fortuitous. We apply the GRT method to compute theoretical RIP at the location and time of 394 large fireballs. Although the predicted spatio-temporal impact distribution matches marginally the observed events, we successfully predict their impact speed distribution.

The discovery of a planet around the closest star to our Sun, Proxima Centauri, represents a quantum leap in the testability of exoplanetary models. Unlike any other discovered exoplanet, models of Proxima b could be contrasted against near future telescopic observations and far future in-situ measurements. In this paper we aim at predicting the planetary radius and the magnetic properties (dynamo lifetime and magnetic dipole moment) of Proxima b analogues (solid planets with masses of $\sim 1-3\,M_\oplus$, rotation periods of several days and habitable conditions). For this purpose we build a grid of planetary models with a wide range of compositions and masses. For each point in the grid we run the planetary evolution model developed in \citet{Zuluaga2013}. Our model assumes small orbital eccentricity, negligible tidal heating and earth-like radiogenic mantle elements abundances. We devise a statistical methodology to estimate the posterior distribution of the desired planetary properties assuming simple prior distributions for the orbital inclination and bulk composition. Our model predicts that Proxima b would have a mass $1.3\leq M_{\rm p}\leq 2.3\,M_{\oplus}$ and a radius $R_{\rm p}=1.4^{+0.3}_{-0.2}\,R_{\oplus}$. In our simulations, most Proxima b analogues develop intrinsic dynamos that last for $\geq$4 Gyr (the estimated age of the host star). If alive, the dynamo of Proxima b have a dipole moment ${\cal M}_{\rm dip}>0.32^{\times 2.3}_{\div 2.9}{\cal M}_{\rm dip,\oplus}$. These results are not restricted to Proxima b but they also apply to earth-like planets having similar observed properties.

In this work we find evidence that the object is of cometary origin.

We undertook observations with the Green Bank Telescope, simultaneously with the 300m telescope in Arecibo, as a follow-up of a possible flare of radio emission from Ross 128. We report here the non-detections from the GBT observations in C band (4-8 GHz), as well as non-detections in archival data at L band (1.1-1.9 GHz). We suggest that a likely scenario is that the emission comes from one or more satellites passing through the same region of the sky.

Potentially habitable planets within the habitable zone of M-dwarfs are affected by tidal interaction. We studied the tidal evolution in GJ 667C using a numerical code we call TIDEV. We reviewed the problem of the dynamical evolution focusing on the effects that a rheological treatment, different compositions and the inclusion of orbital perturbations, have on the spin-down time and the probability to be trapped in a low spin-orbit resonance. Composition have a strong effect on the spin-down time, changing, in some cases, by almost a factor of 2 with respect to the value estimated for a reference Earth-like model. We calculated the time to reach a low resonance value (3:2) for the configuration of 6 planets. Capture probabilities are affected when assuming different compositions and eccentricities variations. We chose planets b and c to evaluate the probabilities of capture in resonances below 5:2 for two compositions: Earth-like and Waterworld planets. We found that perturbations, although having a secular effect on eccentricities, have a low impact on capture probabilities and noth- ing on spin-down times. The implications of the eccentricity variations and actual habitability of the GJ 667C system are discussed.

The discovery of many planets using the Kepler telescope includes ten planets orbiting eight binary stars. Three binaries, Kepler-16, Kepler-47, and Kepler-453, have at least one planet in the circumbinary habitable-zone (BHZ). We constrain the level of high-energy radiation and the plasma environment in the BHZ of these systems. With this aim, BHZ limits in these Kepler binaries are calculated as a function of time, and the habitability lifetimes are estimated for hypothetical terrestrial planets and/or moons within the BHZ. With the time-dependent BHZ limits established, a self-consistent model is developed describing the evolution of stellar activity and radiation properties as proxies for stellar aggression toward planetary atmospheres. Modeling binary stellar rotation evolution, including the effect of tidal interaction between stars in binaries is key to establishing the environment around these systems. We find that Kepler-16 and its binary analogs provide a plasma environment favorable for the survival of atmospheres of putative Mars-sized planets and exomoons. Tides have modified the rotation of the stars in Kepler-47 making its radiation environment less harsh in comparison to the solar system. This is a good example of the mechanism first proposed by Mason et al. Kepler-453 has an environment similar to that of the solar system with slightly better than Earth radiation conditions at the inner edge of the BHZ. These results can be reproduced and even reparametrized as stellar evolution and binary tidal models progress, using our online tool http://bhmcalc.net.

The discovery of rings around extrasolar planets ("exorings") is one of the next breakthroughs in exoplanetary research. Previous studies have explored the feasibility of detecting exorings with present and future photometric sensitivities by seeking anomalous deviations in the residuals of a standard transit light curve fit, at the level of ~100 ppm for Kronian rings. In this work, we explore two much larger observational consequences of exorings: (1) the significant increase in transit depth that may lead to the misclassification of ringed planetary candidates as false-positives and/or the underestimation of planetary density; and (2) the so-called "photo-ring" effect, a new asterodensity profiling effect, revealed by a comparison of the light curve derived stellar density to that measured with independent methods (e.g., asteroseismology). While these methods do not provide an unambiguous detection of exorings, we show that the large amplitude of these effects combined with their relatively simple analytic description, makes them highly suited to large-scale surveys to identify candidate ringed planets worthy of more detailed investigation. Moreover, these methods lend themselves to ensemble analyses seeking to uncover evidence of a population of ringed planets. We describe the method in detail, develop the basic underlying formalism and test it in the parameter space of rings and transit config- uration. We discuss the prospects of using this method for the first systematic search of exoplanetary rings in the Kepler database and provide a basic computational code for implementing it.

The diversity and quantity of moons in the Solar System suggest a manifold population of natural satellites exist around extrasolar planets. Of peculiar interest from an astrobiological perspective, the number of sizable moons in the stellar habitable zones may outnumber planets in these circumstellar regions. With technological and theoretical methods now allowing for the detection of sub-Earth-sized extrasolar planets, the first detection of an extrasolar moon appears feasible. In this review, we summarize formation channels of massive exomoons that are potentially detectable with current or near-future instruments. We discuss the orbital effects that govern exomoon evolution, we present a framework to characterize an exomoon's stellar plus planetary illumination as well as its tidal heating, and we address the techniques that have been proposed to search for exomoons. Most notably, we show that natural satellites in the range of 0.1 - 0.5 Earth mass (i) are potentially habitable, (ii) can form within the circumplanetary debris and gas disk or via capture from a binary, and (iii) are detectable with current technology.

Binaries could provide the best niches for life in the galaxy. Though counterintuitive, this assertion follows directly from stellar tidal interaction theory and the evolution of lower mass stars. There is strong evidence that chromospheric activity of rapidly rotating young stars may be high enough to cause mass loss from atmospheres of potentially habitable planets. The removal of atmospheric water is most critical. Tidal breaking in binaries could help reduce magnetic dynamo action and thereby chromospheric activity in favor of life. We call this the Binary Habitability Mechanism (BHM), that we suggest allows for water retention at levels comparable to or better than Earth. We discuss novel advantages that life may exploit, in these cases, and suggest that life may even thrive on some circumbinary planets. We find that while many binaries do not benefit from BHM, high quality niches do exist for various combinations of stars between 0.55 and 1.0 solar masses. For a given pair of stellar masses, BHM operates only for certain combinations of period and eccentricity. Binaries having a solar-type primary seem to be quite well suited niches having wide and distant habitable zones with plentiful water and sufficient light for photosynthetic life. We speculate that, as a direct result of BHM, conditions may be suitable for life on several planets and possibly even moons of giant planets orbiting some binaries. Lower mass combinations, while more restrictive in parameter space, provide niches lasting many billions of years and are rich suppliers of photosynthetic photons. We provide a publicly available web-site (\href{http://bit.ly/BHM-calculator}{http://bit.ly/BHM-calculator}, \href{http://bit.ly/BHM-calculator-mirror}{http://bit.ly/BHM-calculator-mirror}), which calculates the BHM effects presented in this paper.

We present the results of measuring the geocentric lunar distance using what we propose is the simplest method to achieve a precise result. Although lunar distance has been systematically measured to a precision of few millimeters using powerful lasers and retroreflectors installed on the moon by the Apollo missions, the method devised and applied here can be readily used by nonscientist citizens (e.g. amateur astronomers or students) and it requires only a good digital camera. After launching a citizen science project called the Aristarchus Campaign, intended to involve astronomy enthusiasts in scientific measurement of the Lunar Eclipse of 15 April 2014, we compiled and measured a series of pictures obtained by one of us (J.C. Figueroa). These measurements allowed us to estimate the lunar distance to a precision of 3%. We describe here how to perform the measurements and the method to calculate from them the geocentric lunar distance using only the pictures time stamps and a precise measurement of the instantaneous lunar apparent diameter. Our aim here is not to provide any improved measurement of a well-known astronomical quantity, but rather to demonstrate how the public could be engaged in scientific endeavors and how using simple instrumentation and readily available technological devices such as smartphones and digital cameras, any person can measure the local Universe as ancient astronomers did.

In this paper we discuss and illustrate the hypothesis that life substantially alters the state of a planetary environment and therefore, modifies the limits of the HZ as estimated for an uninhabited planet. This hypothesis lead to the introduction of the Habitable Zone for Inhabited planets (hereafter InHZ), defined here as the region where the complex interaction between life and its abiotic environment is able to produce plausible equilibrium states with the necessary physical conditions for the existence and persistence of life itself. We support our hypothesis of an InHZ with three theoretical arguments, multiple evidences coming from observations of the Earth system, several conceptual experiments and illustrative numerical simulations. Conceptually the diference between the InHZ and the Abiotic HZ (AHZ) depends on unique and robust properties of life as an emergent physical phenomenon and not necesarily on the particular life forms bearing in the planet. Our aim here is to provide conceptual basis for the development of InHZ models incorporating consistently life-environment interactions. Although previous authors have explored the effects of life on habitability there is a gap in research developing the reasons why life should be systematically included at determining the HZ limits. We do not provide here definitive limits to the InHZ but we show through simple numerical models (as a parable of an inhabited planet) how the limits of the AHZ could be modified by including plausible interactions between biota and its environment. These examples aim also at posing the question that if limits of the HZ could be modified by the presence of life in those simple dynamical systems how will those limits change if life is included in established models of the AHZ.

With most planets and planetary candidates detected in the stellar habitable zone (HZ) being super-Earths and gas giants, rather than Earth-like planets, we naturally wonder if their moons could be habitable. The first detection of such an exomoon has now become feasible, and due to observational biases it will be at least twice as massive as Mars. But formation models predict moons can hardly be as massive as Earth. Hence, a giant planet's magnetosphere could be the only possibility for such a moon to be shielded from cosmic and stellar high-energy radiation. Yet, the planetary radiation belt could also have detrimental effects on exomoon habitability. We here synthesize models for the evolution of the magnetic environment of giant planets with thresholds from the runaway greenhouse (RG) effect to assess the habitability of exomoons. For modest eccentricities, we find that satellites around Neptune-sized planets in the center of the HZ around K dwarf stars will either be in an RG state and not be habitable, or they will be in wide orbits where they will not be affected by the planetary magnetosphere. Saturn-like planets have stronger fields, and Jupiter-like planets could coat close-in habitable moons soon after formation. Moons at distances between about 5 and 20 planetary radii from a giant planet can be habitable from an illumination and tidal heating point of view, but still the planetary magnetosphere would critically influence their habitability.

There is a group of newly recognized asteroids in the main belt that are exhibiting cometary characteristics. We will call them Asteroidal Belt Comets or ABCs. An analysis of their orbital and physical properties has resulted in the following conclusion: (1) We define the detached group, DG, as objects sublimating water that have aphelion distances Q < 4.5 AU. The detached group contains all the ABCs traditionally recognized, plus a few other members like 2P and 107P. With the above definition there are 11 members of the ABC group: 2P, 107P, 133P, 176P, 233P, 238P, C/2008 R1, C/2010 R2, 2011 CR42, 3200, and 300163 = 2006 VW139. (2) Active objects sublimate ices except collisioned asteroids, CA, P/2010 A2, 596 Scheila, and P/2012 F5 Gibbs, and 3200 Phaethon who was active due to sputtering. (3)We have calibrated the black body (color) temperature of comets vs perihelion distance, R. We find T = 325+-5 K/SQRT(R). (4) Using a mathematical model of the thermal wave we calculate the thickness of the dust layer on comet nuclei as 2.0+-0.5 m for comet 107P, 4.7+-1.2 m for comet 133P and 1.9+-0.5 m for a sample of 9 comets. (5) We have located three ABCs in an evolutionary diagram of Remaining Revolutions, RR, vs Water-Budget Age, WB-AGE. ABCs lie together in the upper right hand corner of the diagram. (6) RR vs WB-AGE defines the region of the graveyard of comets, objects with 1000 cy < WB-AGE (cy stands for comet years). Five members belong to the graveyard, 107P, 133P, 2006 VW139, D/1891 W1 Blanpain and 3200 Phaeton. Thus we propose that the asteroidal belt contains an enormous graveyard of ancient dormant and extinct rocky comets, that turn on (are rejuvenated), in response to a diminution of their perihelion distance.

We report a mechanism capable of reducing (or increasing) stellar activity in binary stars, thereby potentially enhancing (or destroying) circumbinary habitability. In single stars, stellar aggression towards planetary atmospheres causes mass-loss, which is especially detrimental for late-type stars, because habitable zones are very close and activity is long lasting. In binaries, tidal rotational breaking reduces magnetic activity, thus reducing harmful levels of XUV radiation and stellar mass-loss that are able to erode planetary atmospheres. We study this mechanism for all confirmed circumbinary (p-type) planets. We find that main sequence twins provide minimal flux variation and in some cases improved environments, if the stars rotationally synchronize within the first Gyr. Solar-like twins, like Kepler 34 and Kepler 35, provide low habitable zone XUV fluxes and stellar wind pressures. These wide, moist, habitable zones may potentially support multiple habitable planets. Solar-type stars with lower mass companions, like Kepler 47, allow for protected planets over a wide range of secondary masses and binary periods. Kepler 38 and related binaries are marginal cases. Kepler 64 and analogues have dramatically reduced stellar aggression due to synchronization of the primary, but are limited by the short lifetime. Kepler 16 appears to be inhospitable to planets due to extreme XUV flux. These results have important implications for estimates of the number of stellar systems containing habitable planets in the Galaxy and allow for the selection of binaries suitable for follow-up searches for habitable planets.

Magnetic protection of potentially habitable planets plays a central role in determining their actual habitability and/or the chances of detecting atmospheric biosignatures. We develop here a thermal evolution model of potentially habitable Earth-like planets and super-Earths. Using up-to-date dynamo scaling laws we predict the properties of core dynamo magnetic fields and study the influence of thermal evolution on their properties. The level of magnetic protection of tidally locked and unlocked planets is estimated by combining simplified models of the planetary magnetosphere and a phenomenological description of the stellar wind. Thermal evolution introduces a strong dependence of magnetic protection on planetary mass and rotation rate. Tidally locked terrestrial planets with an Earth-like composition would have early dayside magnetospause distances between 1.5 and 4.0 Rp, larger than previously estimated. Unlocked planets with periods of rotation ~1 day are protected by magnetospheres extending between 3 and 8 Rp. Our results are robust against variations in planetary bulk composition and uncertainties in other critical model parameters. For illustration purposes the thermal evolution and magnetic protection of the potentially habitable super-Earths GL 581d, GJ 667Cc and HD 40307g were also studied. Assuming an Earth-like composition we found that the dynamos of these planets are already extinct or close to being shut down. While GL 581d is the best protected, the protection of HD 40307g cannot be reliably estimated. GJ 667Cc, even under optimistic conditions, seems to be severely exposed to the stellar wind and, under the conditions of our model, has probably suffered massive atmospheric losses.

A ballistic reconstruction of a meteoroid orbit can be made if enough information is available about its trajectory inside the atmosphere. A few methods have been devised in the past and used in several cases to trace back the origin of small impactors. On February 15, 2013, a medium-sized meteoroid hit the atmosphere in the Chelyabinsk region of Russia, causing damage in several large cities. The incident, the largest registered since the Tunguska event, was witnessed by many thousands and recorded by hundreds of amateur and public video recording systems. The amount and quality of the information gathered by those systems is sufficient to attempt a reconstruction of the trajectory of the impactor body in the atmosphere, and from this the orbit of the body with respect to the Sun. Using amateur and public footage taken in four different places close to the event, we have determined precisely the properties of the entrance trajectory and the orbit of the Chelyabinsk event impactor. We found that the object entered the atmosphere at a velocity ranging from 16.0 to 17.4 km/s in a grazing trajectory, almost directly from the east, with an azimuth of velocity vector of 285$^o$, and with an elevation of 15.8$^o$ with respect to the local horizon. The orbit that best fits the observations has, at a 95% confidence level, a semi-major axis a = 1.26$\pm$0.05 AU, eccentricity e = 0.44$\pm$0.03, argument of perihelion $\omega$=95.5$^o\pm2^o$ and longitude of ascending node $\Omega$= 326.5$^o\pm0.3^o$. Using these properties the object can be classified as belonging to the Apollo family of asteroids. The absolute magnitude of the meteoroid was H= 25.8, well below the threshold for its detection and identification as a Potential Hazardous Asteroid (PHA). This result would imply that present efforts intended to detect and characterize PHAs are incomplete.

In February 15 2013 a medium-sized meteoroid impacted the atmosphere in the region of Chelyabinsk, Russia. After its entrance to the atmosphere and after travel by several hun- dred of kilometers the body exploded in a powerful event responsible for physical damages and injured people spread over a region enclosing several large cities. We present in this letter the results of a preliminary reconstruction of the orbit of the Chelyabinsk meteoroid. Using evidence gathered by one camera at the Revolution Square in the city of Chelyabinsk and other videos recorded by witnesses in the close city of Korkino, we calculate the trajectory of the body in the atmosphere and use it to reconstruct the orbit in space of the meteoroid previous to the violent encounter with our planet. In order to account for the uncertainties implicit in the determination of the trajectory of the body in the atmosphere, we use Monte Carlo methods to calculate the most probable orbital parameters. We use this result to classify the meteoroid among the near Earth asteroid families finding that the parent body belonged to the Apollo asteroids. Although semimajor axis and inclination of the preliminary orbit computed by us are uncertain, the rest of orbital elements are well constrained in this preliminary reconstruction.

We present a model for the evolution of the magnetic properties of habitable terrestrial planets and their effects on the protection of planetary atmosphere against the erosive action of stellar wind. Using up-to-date thermal evolution models and dynamo scaling laws we predict the evolution of the planetary dipole moment as a function of planetary mass and rotation rate. Combining these results with models for the evolution of the stellar wind, stellar XUV fluxes and planetary exosphere characteristics, we determine the properties of the magnetosphere and the exobase radius in order to estimate the level of atmospheric mass losses. We use this model to evaluate the magnetic protection of the potentially habitable super-Earths GJ 667Cc, Gl 581d and HD 85512b. We confirm that Earth-like planets, even under the highest attainable magnetic field strengths, will lose a significant fraction of their atmospheric volatiles if they are tidally locked in the habitable zone of dM stars, or even if having N/O-rich atmospheres they are in habitable zones closer than $\sim$ 0.8 AU. Similar mass-dependent inner limits have been found for super-Earths $M_p\gtrsim 3 M_\oplus$ that in any case seem to have better chances of preserving their atmospheres even if they are tidally locked. We predict that the atmosphere of GJ 667Cc has probably already been obliterated and it is presently uninhabitable. On the other hand, our model predicts that the atmospheres of Gl 581d and HD 85512b would be well protected by intrinsic magnetic fields, even under the worst expected conditions of stellar aggression. (abrigded abstract).

Planetary magnetic fields could impact the evolution of planetary atmospheres and have a role in the determination of the required conditions for the emergence and evolution of life (planetary habitability). We study here the role of rotation in the evolution of dynamo-generated magnetic fields in massive earth-like planets, Super Earths (1-10 $M_\oplus$). Using the most recent thermal evolution models of Super Earths (Gaidos et al. 2010; Tachinami et al. 2011) and updated scaling laws for convection-driven dynamos, we predict the evolution of the local Rossby number. This quantity is one of the proxies for core magnetic field regime, i.e. non-reversing dipolar, reversing dipolar and multipolar. We study the dependence of the local Rossby number and hence the core magnetic field regime on planetary mass and rotation rate. Previous works have focused only on the evolution of core magnetic fields assuming rapidly rotating planets i.e. planets in the dipolar regime. In this work we go further, including the effects of rotation in the evolution of planetary magnetic field regime and obtaining global constraints to the existence of intense protective magnetic fields in rapidly and slowly rotating Super Earths. We find that the emergence and continued existence of a protective planetary magnetic field is not only a function of planetary mass but also depend on rotation rate. Low-mass Super Earths ($M\lesssim 2 M_\oplus$) develop intense surface magnetic fields but their lifetimes will be limited to 2-4 Gyrs for rotational periods larger than 1-4 days. On the other hand and also in the case of slowly rotating planets, more massive Super Earths ($M\gtrsim 2 M_\oplus$) have weak magnetic fields but their dipoles will last longer. Finally we analyze tidally locked Super Earths inside and outside the habitable zone of GKM stars.

We study the constraints on neutrino masses that could be derived from the observation of a Galactic supernova neutrino signal with present and future neutrino detectors. Our analysis is based on a recently proposed method that uses the full statistics of neutrino events and does not depend on particular astrophysical assumptions. The statistical approach, originally justified mainly in terms of intuitive reasoning, is put on a more solid basis by means of Bayesian inference reasoning. Theoretical uncertainties in the neutrino signal time profiles are estimated by applying the method to two widely different supernova models. Present detectors can reach a sensitivity down to 1 eV. This is better than limits from tritium $\beta$-decay experiments, competitive with the most conservative results from neutrinoless double $\beta$-decay, less precise but less dependent from prior assumptions than cosmological bounds. Future megaton water Cerencov detectors will allow for about a factor of two improvement. However, they will not be competitive with the next generation of laboratory experiments.

We study the sensitivity to neutrino masses of a Galactic supernova neutrino signal as could be measured with the detectors presently in operation and with future large volume water \v{C}erencov and scintillator detectors. The analysis uses the full statistics of neutrino events. The method proposed uses the principles of Bayesian inference reasoning and has shown a remarkably independence of astrophysical assumptions. We show that, after accounting for the uncertainties in the detailed astrophysical description of the neutrino signal and taking into account the effects of neutrino oscillations in the supernova mantle, detectors presently in operation can have enough sensitivity to reveal a neutrino mass (or to set upper limits) at the level of 1 eV. This is sensibly better than present results from tritium $\beta$-decay experiments, competitive with the most conservative limits from neutrinoless double $\beta$-decay and less precise but remarkably less dependent from prior assumptions than cosmological measurements. Future megaton water \v{C}erencov detectors and large volume scintillator detectors will allow for about a factor of two improvement in the sensitivity; however, they will not be competitive with the next generation of tritium $\beta$-decay and neutrinoless double $\beta$-decay experiments. Using the codes developed to perform the generation of synthetic supernova signals and their analysis we created a computer package, SUNG (SUpernova Neutrino Generation tool, http://urania.udea.edu.co/sungweb), aimed to offer a general purpose solution to perform calculations in supernova neutrino studies.

A new method to study the effects of neutrino masses on a supernova neutrino signal is proposed. The method relies exclusively on the analysis of the full statistics of neutrino events, it is independent of astrophysical assumptions, and does not require the observation of any additional phenomenon to trace possible delays in the neutrino arrival times. The sensitivity of the method to the sub-eV neutrino mass range, defined as the capability of disentangling at 95% c.l. the case $m_\nu=1$eV from $m_\nu=0$, is tested by analyzing a set of synthetic neutrino samples modeled according to the signal that could be detected at SuperKamiokande. For a supernova at the Galactic center success is achieved in more than 50% of the cases. It is argued that a future Galactic supernova yielding several thousands of inverse $\beta$ decays might provide enough information to explore a neutrino mass range somewhat below 1 eV.

A convincing explanation for the observed pulsar large peculiar velocities is still missing. We argue that any viable particle physics solution would most likely involve the resonant production of a non-interacting neutrino $\nu_s$ of mass $m_{\nu_s}\sim 20$--50 keV. We propose a model where anisotropic magnetic field configurations strongly bias the resonant spin flavour precession of tau antineutrinos into $\nu_s$. For internal magnetic fields $B_{int} \gsim 10^{15}$ G a $\bar\nu_\tau$-$\nu_s$ transition magnetic moment of the order of $10^{-12}$ Bohr magnetons is required. The asymmetric emission of $\nu_s$ from the core can produce sizeable natal kicks and account for recoil velocities of several hundred kilometers per second.