We study the dynamical evolution of Jupiter and Saturn embedded in a gaseous, solar-nebula-type disc by means of hydrodynamics simulations with the FARGO2D1D code. We study the evolution for different initial separations of the planets' orbits, Delta a_SJ , to investigate whether they become captured in mean motion resonance (MMR) and the direction of the subsequent migration of the planet (inward or outward). We also provide an assessment of the planet's orbital dynamics at different epochs of Saturn's growth. We find that the evolution of initially compact orbital configurations is dependent on the value of Delta a_SJ . This implies that an evolution as that proposed in the Grand Tack model depends on the precise initial orbits of Jupiter and Saturn and on the timescales for their formation. Capture in the 1:2 MMR and inward or (nearly) stalled migration are highly favoured. Within its limits, our work suggests that the reversed migration, associated with the resonance capture of Jupiter and Saturn, may be a low probability evolutionary scenario, so that other planetary systems with giant planets are not expected to have experienced a Grand Tack-like evolutionary path.

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.

As part of our multi-observatory, multi-filter campaign, we present \rmi color observations of 82 Near-Earth Objects (NEOs) obtained with the RATIR instrument on the 1.5m robotic telescope at the San Pedro Martir's National Observatory in Mexico. Our project is particularly focused on rapid response observations of small ($\lesssim 850$ m) NEOs. The rapid response and the use of spectrophotometry allows us to constrain the taxonomic classification of NEOs with high efficiency. Here we present the methodology of our observations and our result, suggesting that the ratio of C-type to S-type asteroids in a size range of $\sim$30-850m is 1.1, which is in accordance with our previous results. We also find that 10$\%$ of all NEOs in our sample are neither C- nor S-type asteroids

We present a photometric follow-up of transiting exoplanets HAT-P-3b and TrES-3b, observed by using several optical and near-infrared filters, with four small-class telescopes (D = 36--152cm) in the Northern Hemisphere. Two of the facilities present their first scientific results. New 10 HAT-P-3b light curves and new 26 TrES-3b light curves are reduced and combined by filter in order to improve the quality of the photometry. Combined light curves fitting is carried out independently by using two different analysis packages, allowing the corroboration of the orbital and physical parameters in the literature. Results find no differences in the relative radius with the observing filter. In particular, we report for HAT-P-3b a first estimation of the planet-to-star radius Rp/R* = 0.1112+0.0025-0.0026 in the B band which is coherent with values found in the VRIz'JH filters. Concerning TrES-3b, we derive a value for the orbital period of P = 1.3061862+-0.0000001 days which shows no linear variations over nine years of photometric observations.

Using the Wide Field Camera for the United Kingdom Infrared Telescope, we measure the near-infrared colors of near-Earth objects (NEOs) in order to put constraints on their taxonomic classifications. The rapid-response character of our observations allows us to observe NEOs when they are close to the Earth and bright. Here we present near-infrared color measurements of 86 NEOs, most of which were observed within a few days of their discovery, allowing us to characterize NEOs with diameters of only a few meters. Using machine-learning methods, we compare our measurements to existing asteroid spectral data and provide probabilistic taxonomic classifications for our targets. Our observations allow us to distinguish between S-complex, C/X-complex, D-type, and V-type asteroids. Our results suggest that the fraction of S-complex asteroids in the whole NEO population is lower than the fraction of ordinary chondrites in the meteorite fall statistics. Future data obtained with UKIRT will be used to investigate the significance of this discrepancy.

With the use of long-term numerical simulations, we study the evolution and orbital behavior of cometary nuclei in cold Kuiper belt-like debris disks under the gravitational influence of dwarf planets (DPs); we carry out these simulations with and without the presence of a Neptune-like giant planet. This exploratory study shows that in the absence of a giant planet, 10 DPs are enough to induce strong radial and vertical heating on the orbits of belt particles. On the other hand, the presence of a giant planet close to the debris disk, acts as a stability agent reducing the radial and vertical heating. With enough DPs, even in the presence of a Neptune-like giant planet some radial heating remains; this heating grows steadily, re-filling resonances otherwise empty of cometary nuclei. Specifically for the solar system, this secular process seems to be able to provide material that, through resonant chaotic diffusion, increase the rate of new comets spiraling into the inner planetary system, but only if more than the $\sim10$ known DP sized objects exist in the trans-Neptunian region.

Three optical telescopes located at the San Pedro M\'artir National Observatory were used for the first time to obtain multi-filter defocused photometry of the transiting extrasolar planets WASP-39b and WASP-43b. We observed WASP-39b with the 2.12m telescope in the U filter for the first time, and additional observations were carried out in the R and I filters using the 0.84m telescope. WASP-43b was observed in VRI with the same instrument, and in the i filter with the robotic 1.50m telescope. We reduced the data using different pipelines and performed aperture photometry with the help of custom routines, in order to obtain the light curves. The fit of the light curves (1.5--2.5mmag rms), and of the period analysis, allowed a revision of the orbital and physical parameters, revealing for WASP-39b a period ($4.0552947 \pm 9.65 \times 10^{-7}$ days) which is $3.084 \pm 0.774$ seconds larger than previously reported. Moreover, we find for WASP-43b a planet/star radius ($0.1738 \pm 0.0033$) which is $0.01637 \pm 0.00371$ larger in the i filter with respect to previous works, and that should be confirmed with additional observations. Finally, we confirm no evidence of constant period variations in WASP-43b.

The orbital elements of comet Halley are known to a very high precision, suggesting that the calculation of its future dynamical evolution is straightforward. In this paper we seek to characterize the chaotic nature of the present day orbit of comet Halley and to quantify the timescale over which its motion can be predicted confidently. In addition, we attempt to determine the timescale over which its present day orbit will remain stable. Numerical simulations of the dynamics of test particles in orbits similar to that of comet Halley are carried out with the Mercury 6.2 code. On the basis of these we construct survival time maps to assess the absolute stability of Halley's orbit, frequency analysis maps, to study the variability of the orbit and we calculate the Lyapunov exponent for the orbit for variations in initial conditions at the level of the present day uncertainties in our knowledge of its orbital parameters. On the basis of our calculations of the Lyapunov exponent for comet Halley, the chaotic nature of its motion is demonstrated. The e-folding timescale for the divergence of initially very similar orbits is approximately 70 years. The sensitivity of the dynamics on initial conditions is also evident in the self-similarity character of the survival time and frequency analysis maps in the vicinity of Halley's orbit, which indicates that, on average, it is unstable on a timescale of hundreds of thousands of years. The chaotic nature of Halley's present day orbit implies that a precise determination of its motion, at the level of the present day observational uncertainty, is difficult to predict on a timescale of approximately 100 years. Furthermore, we also find that the ejection of Halley from the solar system or its collision with another body could occur on a timescale as short as 10,000 years.

To date, 17 circumbinary planets have been discovered. In this paper, we focus our attention on the stability of the Kepler circumbinary planetary systems with only one planet, i.e. Kepler-16, Kepler-34, Kepler-35, Kepler-38, Kepler-64 and Kepler-413. In addition to their intrinsic interest, the study of such systems is an opportunity to test our understanding of planetary system formation and evolution around binaries. The investigation is done by means of numerical simulations. We perform numerical integrations of the full equations of motion of each system with the aim of checking the stability of the planetary orbit. The investigation of the stability of the above systems consists of three numerical experiments. In the first one we perform a long term (1Gyr) numerical integration of the nominal solution of the six Kepler systems under investigation. In the second experiment, we look for the critical semimajor axis of the six planetary orbits, and finally, in the third experiment, we construct two dimensional stability maps on the eccentricity-pericentre distance plane. Additionally, using numerical integrations of the nominal solutions we checked if this solutions were close to the exact resonance.

We study the distribution of intragroup dark matter in small groups of dark matter halos of galaxy-like size in a LCDM cosmology. These groups are identified using a physical criterion and may be an appropriate representation of small galaxy groups. We quantify the amount of intragroup dark matter and characterize its distribution. We find that compact associations of halos, as well as those intermediate and loose groups, have rather flat intragroup dark matter profiles with logarithmic slopes of gamma ~0 and ~-0.2$, respectively. Hence, the intragroup dark matter of these halo systems does not follow the same cuspy tendency that halos of galaxies have. In intermediate and loose galaxy-size halo associations the intragroup matter tends to be <50% that of the total mass of the group, and in compact associations is <20% within their group radius.

We suggest a new web-based approach for browsing and visualizing data produced by a network of telescopes, such as those of the ongoing TAOS and the forthcoming TAOS II projects. We propose a modern client-side technology and we present two examples based on two software packages developed for different kinds of server- side database approaches. In spite our examples are specific for the browsing of TAOS light curves, the software is coded in a way to be suitable for the use in several types of astronomical projects.

A preliminary data analysis of the stellar light curves obtained by the robotic telescopes of the TAOS project is presented. We selected a data run relative to one of the stellar fields observed by three of the four TAOS telescopes, and we investigate the common trend and the correlation between the light curves. We propose two ways to remove these trends and show the preliminary results. A project aimed at flagging interesting behaviors, such as stellar variability, and to set up an automated follow-up with the San Pedro M\'artir Facilities is on the way.

We study the effect of a massive planetesimal disk on the dynamical stability of the outer planets assuming, as has been suggested recently, that these were initially locked in a compact and multiresonant configuration as a result of gas-driven migration in a protoplanetary disk. The gravitational interaction among all bodies in our simulations is included self-consistently using the Mercury6.5 code. Several initial multiresonant configurations and planetesimal disk models are considered. Under such conditions a strong dynamical instability, manifested as a rapid giant planet migration and planetesimal disk dispersal, develops on a timescale of less than 40 Myr in most cases. Dynamical disk heating due to the gravitational interactions among planetesimals leads to more frequent interactions between the planetesimals and the ice giants Uranus and Neptune, in comparison to models in which planetesitmal-planetesimal interactions are neglected. On account of the rapid evolution of the multiresonant configurations obtained with fully self-consistent simulations, our results are inconsistent with the dynamical instability origin of the Late Heavy Bombardment as currently considered by the Nice model for the Solar System.

Measurements conducted with the ASPERA-4 instrument and the magnetometer of the Venus Express spacecraft show that the dynamic pressure of planetary O+ ion fluxes measured in the Venus wake can be significantly larger than the local magnetic pressure and, as a result, those ions are not being driven by magnetic forces but by the kinetic energy of the solar wind. Beams of planetary O+ ions with those properties have been detected in several orbits of the Venus Express through the wake as the spacecraft traverses by the noon-midnight plane along its near polar trajectory. The momentum flux of the O+ ions leads to superalfvenic flow conditions. It is suggested that such O+ ion beams are produced in the vicinity of the magnetic polar regions of the Venus ionosphere where the solar wind erodes the local plasma leading to plasma channels that extend downstream from those regions.

It has been suggested that the ejection to interplanetary space of terrestrial crustal material, accelerated in a large impact, may result in the interchange of biological material between Earth and other Solar System bodies. In this paper, we analyze the fate of debris ejected from Earth by means of numerical simulations of the dynamics of a large collection of test particles. This allows us to determine the probability and conditions for the collision of ejecta with other planets of the Solar System. We also estimate the amount of particles falling-back to Earth as a function of time after being ejected. We find that, in general, the collision rates of Earth ejecta with Venus and the Moon, as well as the fall-back rates, are consistent with results reported in the literature. By considering a larger number of particles than in all previous calculations we have also determined directly the collision probability with Mars and, for the first time, computed collision probabilities with Jupiter. We find that the collision probability with Mars is greater than values determined from collision cross section estimations previously reported.

We study some of the effects that magnetic turbulent fluctuations have on the dynamics of pickup O+ ions in the magnetic polar regions of the Mars ionosheath. In particular we study their effect on the bulk velocity profiles of ions as a function of altitude over the magnetic poles, in order to compare them with recent Mars Express data; that indicate that their average velocity is very low and essentially in the anti-sunward direction. We find that, while magnetic field fluctuations do give rise to deviations from simple ExB-drift gyromotion, even fluctuation amplitudes much greater than those of in situ measurements are {\it not} able to reproduce the vertical velocity profile of O+ ions. We conclude that other physical mechanisms, different from a pure charged particle dynamics, are acting on pickup ions at the Martian terminator. A possibility being a viscous-like interaction between the Solar Wind and the Martian ionosphere at low altitudes.

We study the dynamics of the upsilon-Andromedae planetary system proposed by Curiel et al. (2011). We focus on the study of the 3:1 Mean Motion Resonance between upsilon-Andromedae-d and the recently discovered upsilon-Andromedae-e (hereafter upsilon-And-d and upsilon-And-e). Numerical simulations of the dynamics of the four planet system are conducted. The previously reported apsidal resonance between upsilon-And-c and upsilon-And-d is confirmed. In addition, we find that upsilon-And-d and upsilon-And-e are also in an apsidal resonance condition. Our results further indicate that the upsilon-Andromedae planetary system configuration is in the middle of a stability island in the semimajor axis-eccentricity domain. Additionally, we performed numerical integrations of the planetary configuration over 500 Myr and found it to be stable. We conclude that, within the uncertainties in the value of the orbital parameters, it is likely that upsilon-Andromedae planetary system will remain stable for a long timescale.

We compare the results of the numerical simulation of the viscous-like interaction of the solar wind with the plasma tail of a comet, with velocities of H2O+ ions in the tail of comet Swift-Tuttle determined by means of spectroscopic, ground based observations. Our aim is to constrain the value of the basic parameters in the viscous-like interaction model: the effective Reynolds number of the flow and the interspecies coupling timescale. We find that in our simulations the flow rapidly evolves from an arbitrary initial condition to a quasi-steady state for which there is a good agreement between the simulated tailward velocity of H2O+ ions and the kinematics derived from the observations. The fiducial case of our model, characterized by a low effective Reynolds number (Re \approx 20 and selected on the basis of a comparison to in situ measurements of the plasma flow at comet Halley, yields an excellent fit to the observed kinematics. Given the agreement between model and observations, with no ad hoc assumptions, we believe that this result suggests that viscous-like momentum transport may play an important role in the interaction of the solar wind and the cometary plasma environment.

We model the interaction of the solar wind with the plasma tail of a comet by means of numerical simulations, taking into account the effects of viscous-like forces.A 2D hydrodynamical, two species, finite difference code has been developed for the solution of the time dependent continuity, momentum and energy conservation equations, as applied to the problem at hand. Velocity, density and temperature profiles across the tail are obtained. Several cases with different flow parameters are considered in order to study the relative importance of viscous-like effects and the coupling between species on the flow dynamics. Assuming a Mach number equal to 2 for the incident solar wind as it flows past the comet nucleus, the flow exhibits three transitions with location and properties depending on the Reynolds number for each species and on the ratio of the timescale for inter-species coupling to the crossing time of the free flowing solar wind. By comparing our results with the measurements taken in situ by the Giotto spacecraft during its flyby of comet Halley we constrain the flow parameters for both plasmas. In the context of our approximations, we find that our model is qualitatively consistent with the in situ measurements as long as the Reynolds number of the solar wind protons and of cometary H2O+ ions is low, less than 100, suggesting that viscous-like momentum transport processes may play an important role in the interaction of the solar wind and the plasma environment of comets.

The density profiles of dark halos are usually inferred from the rotation curves of disk galaxies based on the assumption that the gas is a good tracer of the gravitational potential of the galaxies. Some authors have suggested that magnetic pinching forces could alter significantly the rotation curves of spiral galaxies. In contrast to other studies which have concentrated in the vertical structure of the disk, here we focus on the problem of magnetic confinement in the radial direction to bound the magnetic effects on the HI rotation curves. It is shown that azimuthal magnetic fields can hardly speed up the HI disk of galaxies as a whole. In fact, based on virial constraints we show that the contribution of galactic magnetic fields to the rotation curves cannot be larger than ~10 km/s at the outermost point of HI detection, if the galaxies did not contain dark matter at all, and up to 20 km/s in the conventional dark halo scenario. The procedure to estimate the maximum effect of magnetic fields is general and applicable to any particular galaxy disk. The inclusion of the surface terms, namely the intergalactic (thermal, magnetic or ram) pressure, does not change our conclusions. Other problems related with the magnetic alternative to dark halos are highlighted. The relevance of magnetic fields in the cuspy problem of dark halos is also discussed.

We examine the linear stability of a flow threaded by a weak, vertical magnetic field in a disk with a keplerian rotation profile and a vertical stratification of the ionization degree as that predicted for vast portions of protoplanetary disks. A quasi-global analysis is carried out, where the form of the perturbations in the vertical direction is determined. Considering the ohmic magnetic diffusivity of the gas, the conditions leading to the magnetorotational instability are analyzed as a function of the diffusivity at the disk surfaces, its vertical profile and the strength of the unperturbed magnetic field. For typical conditions believed to prevail in protoplanetary disks at radial distances between 0.1 and 10 AU, where the so-called dead zone is proposed to exist, we find that generally the instability is damped. This implies that, if the MRI is considered the only possible source of turbulence in protoplanetary disks, no viscous angular momentum transport occursat those radii.