In this work, we address the effects of a phenomenon known as Verlinde gravity. Here we show that its effect over the planets and the Moon in our solar system is quite negligible. We find that the Verlinde gravity effects on the orbits of planets are at least 10 times smaller than the precision with which we can determine the Sun's mass, and the one on the orbit of the Moon is about 100 times smaller than the precision with which we can determine the Earth's mass. These results let us infer that statements in the literature that Verlinde gravity is ruled out by the observed motion of planets in our solar system aren't correct.

The goal of this research is to study how the fragmentation of planetary embryos can affect the physical and dynamical properties of terrestrial planets around solar-type stars. Our work focuses on the formation and evolution of planets and water delivery in the Habitable Zone. We distinguish class A and class B HZ planets, which have an accretion seed initially located inside and beyond the snowline, respectively. We develop an Nbody integrator that incorporates fragmentation and hit-and-run collisions, which is called D3 N-body code. From this, we perform 46 numerical simulations of planetary accretion in systems that host two gaseous giants like Jupiter and Saturn. We compare two sets of 23 N-body simulations, one of which includes a realistic collisional treatment and the other one models all impacts as perfect mergers. The final masses of the HZ planets formed in runs with fragmentation are about 15\%-20\% smaller than those obtained without fragmentation. As for the class A HZ planets, those formed in simulations without fragmentation experience very significant increases in mass respect to their initial values, while the growth of those produced in runs with fragmentation is less relevant. We remark that the fragments play a secondary role in the masses of the class A HZ planets, providing less than 30\% of their final values. In runs without fragmentation, the final fraction of water of the class A HZ planets keeps the initial value since they do not accrete water-rich embryos. In runs with fragmentation, the final fraction of water of such planets strongly depends on the model used to distribute the water after each collision. The class B HZ planets do not show significant differences concerning their final water contents in runs with and without fragmentation. From this, we find that the collisional fragmentation is not a barrier to the survival of water worlds in the HZ.

In the present research, we study the effects of a single giant planet in the dynamical evolution of water-rich embryos and planetesimals, located beyond the snow line of systems around Sun-like stars, in order to determine what kind of terrestrial-like planets could be formed in the habitable zone (hereafter HZ) of these systems. To do this, we carry out N-body simulations of planetary accretion, considering that the gas has been already dissipated from the disk and a single giant planet has been formed beyond the snow line of the system, at 3 au. We find that a giant planet with a value of mass between Saturn-mass and Jupiter-mass, represents a limit from which the amount of water-rich embryos that moves inward from beyond the snow line starts to decrease. From this, our research suggests that giant planets more massive than one Jupiter-mass become efficient dynamical barriers to inward-migrating water-rich embryos. Moreover, we infer that the number of these embryos that survive in the HZ significantly decreases for systems that host a giant planet more massive than one Jupiter-mass. This result has important consequences concerning the formation of terrestrial-like planets in the HZ with very high water contents and could provide a selection criteria in the search of potentially habitable exoplanets in systems that host a gaseous giant around solar-type stars.

Several observational works have shown the existence of Jupiter-mass planets covering a wide range of semi-major axes around Sun-like stars. We aim to analyse the planetary formation processes around Sun-like stars that host a Jupiter-mass planet at intermediate distances ranging from $\sim$1 au to 2 au. Our study focusses on the formation and evolution of terrestrial-like planets and water delivery in the habitable zone (HZ) of the system. Our goal is also to analyse the long-term dynamical stability of the resulting systems. A semi-analytic model was used to define the properties of a protoplanetary disk that produces a Jupiter-mass planet around the snow line, which is located at $\sim$2.7 au for a solar-mass star. Then, it was used to describe the evolution of embryos and planetesimals during the gaseous phase up to the formation of the Jupiter-mass planet, and we used the results as the initial conditions to carry out N-body simulations of planetary accretion. Our simulations produce three different classes of planets in the HZ: 'water worlds', with masses between 2.75 $M_{\oplus}$ and 3.57 $M_{\oplus}$ and water contents of 58% and 75% by mass, terrestrial-like planets, with masses ranging from 0.58 $M_{\oplus}$ to 3.8 $M_{\oplus}$ and water contents less than 1.2% by mass, and 'dry worlds', simulations of which show no water. A relevant result suggests the efficient coexistence in the HZ of a Jupiter-mass planet and a terrestrial-like planet with a percentage of water by mass comparable to the Earth. Moreover, our study indicates that these planetary systems are dynamically stable for at least 1 Gyr. Systems with a Jupiter-mass planet located at 1.5 au - 2 au around solar-type stars are of astrobiological interest. These systems are likely to harbour terrestrial-like planets in the HZ with a wide diversity of water contents.

We have developed a computationally competitive N-body model of a previrialized aggregation of galaxies in a flat LambdaCDM universe to assess the role of the multiple mergers that take place during the formation stage of such systems in the configuration of the remnants assembled at their centres. An analysis of a suite of 48 simulations of low-mass forming groups (of about 1E13 solar masses) demonstrates that the gravitational dynamics involved in their hierarchical collapse is capable of creating realistic first-ranked galaxies without the aid of dissipative processes. Our simulations indicate that the brightest group galaxies (BGGs) constitute a distinct population from other group members, sketching a scenario in which the assembly path of these objects is dictated largely by the formation of their host system. We detect significant differences in the distribution of Sersic indices and total magnitudes, as well as a luminosity gap between BGGs and the next brightest galaxy that is positively correlated with the total luminosity of the parent group. Such gaps arise from both the grow of BGGs at the expense of lesser companions and the decrease in the relevance of second-ranked objects in equal measure. This results in a dearth of intermediate-mass galaxies which explains the characteristic central dip detected in their luminosity functions in dynamically young galaxy aggregations. The fact that the basic global properties of our BGGs define a thin mass fundamental plane strikingly similar to that followed giant early-type galaxies in the local universe reinforces confidence in the results obtained.

This paper reports one recent result from a set of pre-virialized galaxy group simulations that are being used in an investigation of measurement techniques for the quantity of intragroup light (IGL). We present evidence that the binding energy of the stellar material stripped from the galaxies is essentially uncorrelated with the local mass density. This suggests that IGL detection methods based on the distribution of luminosity perform poorly in detecting the unbound stars.

The reader can find in the literature a lot of different techniques to study the dynamics of a given system and also, many suitable numerical integrators to compute them. Notwithstanding the recent work of Maffione et al. (2011a) for mappings, a detailed comparison among the widespread indicators of chaos in a general system is still lacking. Such a comparison could lead to select the most efficient algorithms given a certain dynamical problem. Furthermore, in order to choose the appropriate numerical integrators to compute them, more comparative studies among numerical integrators are also needed. This work deals with both problems. We first extend the work of Maffione et al. (2011) for mappings to the 2D H\'enon & Heiles (1964) potential, and compare several variational indicators of chaos: the Lyapunov Indicator (LI); the Mean Exponential Growth Factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI) and its generalized version, the Generalized Alignment Index (GALI); the Fast Lyapunov Indicator (FLI) and its variant, the Orthogonal Fast Lyapunov Indicator (OFLI); the Spectral Distance (D) and the Dynamical Spectras of Stretching Numbers (SSNs). We also include in the record the Relative Lyapunov Indicator (RLI), which is not a variational indicator as the others. Then, we test a numerical technique to integrate Ordinary Differential Equations (ODEs) based on the Taylor method implemented by Jorba & Zou (2005) (called taylor), and we compare its performance with other two well-known efficient integrators: the Prince & Dormand (1981) implementation of a Runge-Kutta of order 7-8 (DOPRI8) and a Bulirsch-St\"oer implementation. These tests are run under two very different systems from the complexity of their equations point of view: a triaxial galactic potential model and a perturbed 3D quartic oscillator.

The aim of this research work is to compare the reliability of several variational indicators of chaos on mappings. The Lyapunov Indicator (LI); the Mean Exponential Growth factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI); the Fast Lyapunov Indicator (FLI); the Dynamical Spectra of stretching numbers (SSN) and the corresponding Spectral Distance (D); and the Relative Lyapunov Indicator (RLI), which is based on the evolution of the difference between two close orbits, have been included. The experiments presented herein allow us to reliably suggest a group of chaos indicators to analyze a general mapping. We show that a package composed of the FLI and the RLI (to analyze the phase portrait globally) and the MEGNO and the SALI (to analyze orbits individually) is good enough to make a description of the systems' dynamics.

We present an adaptation of the standard scenario of disk-galaxy formation to the concordant LCDM cosmology aimed to derive analytical expressions for the scale length and rotation speed of present-day disks that form within four different, cosmologically motivated protogalactic dark matter halo-density profiles. We invoke a standard galaxy-formation model that includes virial equilibrium of spherical dark halos, specific angular momentum conservation during gas cooling, and adiabatic halo response to the gas inflow. The mean mass-fraction and mass-to-light ratio of the central stellar disk are treated as free parameters whose values are tuned to match the zero points of the observed size-luminosity and circular speed-luminosity relations of galaxies. We supply analytical formulas for the characteristic size and rotation speed of disks built inside Einasto r^{1/6}, Hernquist, Burkert, and Navarro-Frenk-White dark matter halos. These expressions match simultaneously the observed zero points and slopes of the different correlations that can be built in the RVL space of disk galaxies from plausible values of the galaxy- and star-formation efficiencies.