Thermodynamical and dynamical evolution of a self-gravitating object and the influence of ionization

We use a simple equation of state, in which the adiabatic index depends on opacity and ionization and we integrate the dynamical and thermodynamical equations for the gravitational collapse of a typical solar composition protocloud, up to the virialization of the energies. Following the evolution of the thermal energy and ionization fraction, violent bounces are obtained at the sudden hardening of the equation of state, when the material becomes ionized.We also suggest a mechanism to explain the onset of protostellar winds.We introduce radiation losses in the model, and integrate again the modified equations, studying the evolution of a 1.1M protocloud. The object's effective temperature stays in a confined small zone of the IR region throughout its fast (40 yr) evolution and its luminosity oscillates and decreases from 5000L to 500L . The radius starts from 35 AU and shrinks down to 140R , before a physical instability gives birth to a strong shock wave with consequent mass loss.     

The dynamical evolution of substructure

The evolution of substructure embedded in non-dissipative dark haloes is studied through N -body simulations of isolated systems, both in and out of initial equilibrium, complementing cosmological simulations of the growth of structure. We determine by both analytic calculations and direct analysis of the N -body simulations the relative importance of various dynamical processes acting on the clumps, such as the removal of material by global tides, clump–clump heating, clump–clump merging and dynamical friction. The ratio of the internal clump velocity dispersion to that of the dark halo is an important parameter; as this ratio approaches a value of unity, heating by close encounters between clumps becomes less important, while the other dynamical processes continue to increase in importance. Our comparison between merging and disruption processes implies that spiral galaxies cannot be formed in a protosystem that contains a few large clumps, but can be formed through the accretion of many small clumps; elliptical galaxies form in a more clumpy environment than do spiral galaxies. Our results support the idea that the central cusp in the density profiles of dark haloes is the consequence of self-limiting merging of small, dense haloes. This implies that the collapse of a system of clumps/substructure is not sufficient to form a cD galaxy, with an extended envelope; plausibly, subsequent accretion of large galaxies is required. The post-collapse system is in general triaxial, with rounder systems resulting from fewer, but more massive, clumps. Persistent streams of material from disrupted clumps can be found in the outer regions of the final system, and at an overdensity of around 0.75, can cover 10 to 30 per cent of the sky.     

The dynamical evolution of Saturn's E-ring

The dynamical evolution of fine dust particles ejected from Enceladus and subsequently electrically charged within the Saturnian magnetosphere is studied. It is shown that the gyro-phase drift, which is radially outwards due to the strong radial temperature and density gradients in the magnetospheric plasma, is, by far, the fastest transport mechanism of these grains. Maintenance of the E-ring in a steady state throughout the age of the solar system would need a mass loss from Enceladus of about 2 parts in 1000.     

The dynamical evolution of stellar superclusters

Recent images taken with the Hubble Space Telescope ( HST ) of the interacting disc galaxies NGC 4038/4039 (the Antennae) reveal clusters of many dozens and possibly hundreds of young compact massive star clusters within projected regions spanning about 100 to 500 pc. It is shown here that a large fraction of the individual star clusters merge within a few tens to a hundred Myr. Bound stellar systems with radii of a few hundred parsecs, masses ≲ 10⁹ M⊙ and relaxation times of 10¹¹ − 10¹² yr may form from these. These spheroidal dwarf galaxies contain old stars from the pre-merger galaxy and much younger stars formed in the massive star clusters, and possibly from later gas accretion events. The possibility that star formation in the outer regions of gas-rich tidal tails may also lead to superclusters is raised. The mass-to-light ratio of these objects is small, because they contain an insignificant amount of dark matter. After many hundred Myr such systems may resemble dwarf spheroidal satellite galaxies with large apparent mass-to-light ratios, if tidal shaping is important.     

The evolution of the lunar orbit revisited. I

After recalling the contribution of Halley, J. Kepler, and G. Darwin to our understanding of the secular acceleration of the Moon, we establish a set of differential equations for the variation of the semi-major axis, and the inclination of the Moon on the maximum area plane. These equations are obtained without expanding the disturbing function, due to the tidal bulge, in term of the elliptic elements. The equations thus obtained are simple enough to allow us a qualitative discussion of the solution, followed by a numerical integration.The results obtained show the Moon was in the distant past in a retrograde orbit, approaching the Earth, its inclination increasing towards 90°; once after a closer approach to the Earth, the Moon receeded and it will finally reach an equilibrium point, the orbital and the equatorial planes being blended.The solution of the equations appears as a fascicle of curves, becoming extremely dense as we come nearer to the present. Owing to the high sensitivity of the solution to the initial conditions, a weak disturbance added to our modeled forces may lead to a past situation very different from the conclusion drawn by Goldreich (1966) and MacDonald (1964); the minimal approach distance could be greater than 10 Earth's radii.     

The evolution of the lunar orbit revisited,II

We present here a model for the tidal evolution of an isolated two-body system. Equations are derived, including the dissipation in the planet as in the satellite, in a frequency dependent lag model. The set of differential equations obtained is still valid for large eccentricity, as well as for all inclinations. The reference plane chosen enables us to study the evolution for both the orbital plane and the equatorial plane.The results obtained show the Moon, after having approached the Earth with small variations for the inclination and the eccentricity, exhibits strong increase for the two parameters in the vicinity of the closest approach. In every case the eccentricity tends towards the value 1, whereas the variations of the in clinations are dependent on the magnitude of the dissipation in the satellite.Some qualitative results are also investigated for the final behaviour of satellites such as Triton and the Galilean satellites.     

On the evolution of the lunar orbit

It is generally accepted that the Earth-Moon separation is at present increasing due to tidal dissipation. Values for the corresponding lunar deceleration and the related slowing of the Earth's rotation are obtained from astronomical observations and by studies of ancient eclipses. Extrapolation of these values leads to a close approach of the Earth and Moon 1–3 b.y. BP. However, justification for such an extrapolation is required. It has been hypothesized that periodicities in the Precambrian stromatolites can be used to determine the number of solar days in a lunar month prior to 500 m.y. BP. These data combined with dynamic constraints on the number of solar days in a lunar month indicate a close approach of the Earth and Moon at 2.85 ± 0.25 b.y. BP. It is suggested that the mare volcanism on the Moon and high-temperature Archean volcanism on the Earth prior to this date were caused by tidal heating. It is also suggested that the strong tidal heating during a close approach could have contributed to the formation of the first living organisms.     

Modeling the dynamical profile of a BL Lacertae object

A simple dynamical model for a BL Lacertae active galaxy is presented. The model consists of a logarithmic potential with an additional term representing internal perturbations. The time independent and the evolving model are investigated. In both cases we search for regular and chaotic motion and study the velocity distribution near the centre of the system. Numerical calculations suggest that responsible for the chaotic phenomena is the internal perturbation, the flattening parameter and the dense nucleus. The radius of the nucleus also affects the maximum velocity in the central regions of the galaxy. Our numerical outcomes are supported by theoretical arguments and analytical calculations. A linking of our numerical outcomes to observational data is also presented. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The collisional and dynamical evolution of the main-belt and NEA size distributions

The size distribution of main belt of asteroids is determined primarily by collisional processes. Large asteroids break up and form smaller asteroids in a collisional cascade, with the outcome controlled by the strength-size relationship for asteroids. In addition to collisional processes, the non-collisional removal of asteroids from the main belt (and their insertion into the near-Earth asteroid (NEA) population) is critical, and involves several effects: strong resonances increase the orbital eccentricity of asteroids and cause them to enter the inner planet region; chaotic diffusion by numerous weak resonances causes a slow leak of asteroids into the Mars- and Earth-crossing populations; and the Yarkovsky effect, a radiation force on asteroids, is the primary process that drives asteroids into these resonant escape routes. Yarkovsky drift is size-dependent and can modify the main-belt size distribution. The NEA size distribution is primarily determined by its source, the main-belt population, and by the size-dependent processes that deliver bodies from the main belt. All of these effects are simulated in a numerical collisional evolution model that incorporates removal by non-collisional processes. We test our model against a wide range of observational constraints, such as the observed main-belt and NEA size distributions, the number of asteroid families, the preserved basaltic crust of Vesta and its large south-pole impact basin, the cosmic ray exposure ages of meteorites, and the cratering records on asteroids. We find a strength-size relationship for main-belt asteroids and non-collisional removal rates from the main belt such that our model fits these constraints as best as possible within the parameter space we explore. Our results are consistent with other independent estimates of strength and removal rates.     

On the evolution of the orbit of Nereid

The evolution of the orbit of the second satellite of Neptune (Nereid) on time intervals up to 500 thousand years is investigated. The methods used by authors earlier are supplemented with the possibility of considering the influence of the attraction of the internal satellite Triton on the evolution of the orbit of the external satellite Nereid.     

A minimal dynamical model for tidal synchronization and orbit circularization

We study tidal synchronization and orbit circularization in a minimal model that takes into account only the essential ingredients of tidal deformation and dissipation in the secondary body. In previous work we introduced the model (Escribano et al. in Phys. Rev. E, 78:036216, 2008); here we investigate in depth the complex dynamics that can arise from this simplest model of tidal synchronization and orbit circularization. We model an extended secondary body of mass m by two point masses of mass m/2 connected with a damped spring. This composite body moves in the gravitational field of a primary of mass M ≫ m located at the origin. In this simplest case oscillation and rotation of the secondary are assumed to take place in the plane of the Keplerian orbit. The gravitational interactions of both point masses with the primary are taken into account, but that between the point masses is neglected. We perform a Taylor expansion on the exact equations of motion to isolate and identify the different effects of tidal interactions. We compare both sets of equations and study the applicability of the approximations, in the presence of chaos. We introduce the resonance function as a resource to identify resonant states. The approximate equations of motion can account for both synchronization into the 1:1 spin-orbit resonance and the circularization of the orbit as the only true asymptotic attractors, together with the existence of relatively long-lived metastable orbits with the secondary in p:q (p and q being co-prime integers) synchronous rotation.     

Chemical and dynamical evolution of spiral galaxies

We investigate the consequences of the hypothesis of the secular evolution (growth of the bulge from disc material via a bar and temporal evolution of the Hubble sequence) on the chemical evolution of a galaxy. We present the first dynamical and chemical results of our 3D tree-SPH simulations. This revised version was published online in August 2006 with corrections to the Cover Date.     

Chemical and dynamical evolution of galactic discs

The relative roles of star formation and viscosly-induced radial flows in galactic discs are discussed. It is shown that the present-day distributions of stars, gas, and metals in galactic spirals need not reflect initial conditions but may instead indicate a cooperation between star forming and viscous processes over the disc lifetime.     

Chemical and dynamical evolution of spiral galaxies

We present the very first results of a new 3D numerical model for the formation and evolution of spiral galaxies along the Hubble sequence. We take into account the hydrodynamical properties of the gas with an SPH method while we use a tree code for the gravitational forces of the dark matter and stars. The chemical evolution is also fully included, with both SNe Ia and SNe II explosions being followed, and this will allows us to predict abundances of various chemical species, abundance ratios and their radial distributions. This revised version was published online in September 2006 with corrections to the Cover Date.     

On the dynamical stability of the very low-mass object Gliese 22 Bb

Gliese 22 is a hierarchical triple red dwarf system formed by two close components, Aa and Ab, and a distant component, B, which is moving around the center of mass of the first two.The possible existence of a fourth very low-mass object (15 Jupiter mass) orbiting around component B was reported by Docobo et al. (Docobo, J.A. et al. [2007]. IAU Commun. 26, 3-4). In this probable scenario with four bodies, component B would be in reality two: star Ba and the new object, Bb.Two full three-dimensional accurate (circular and elliptical) solutions for the orbit of Bb have been obtained, along with an improved arrangement of the system masses. In addition, such a multiple system is analyzed by means of a (2 + 2)-body model considering its evolution during 10 Myr. In particular, we have studied its apsidal motion in order to eventually find any evidence of chaotic behavior. The nature of the new object as a giant planet or a brown star is also discussed.     

The dynamical evolution of uniformly rotating asteroids subject to YORP

A detailed derivation is given of the effect of solar radiation on the rotational dynamics of asteroids, commonly called the YORP effect. The current derivation goes beyond previous discussions published in the literature and provides a comprehensive secular dynamical analysis of the effect of solar radiation torques acting on a uniformly rotating body, and the evolution of its rotation state over time. Our predicted model has the global radiation properties of the asteroid as explicit parameters, and hence can be specified independent of these parameters. The resulting secular equations for the rotation rate and rotation pole are characterized by three parameters of the body's shape and explicitly includes the effect of thermal inertia on the evolution of these rotation state parameters. With this detailed model, in conjunction with estimated asteroid shapes and poles, we compute the expected YORP torques and dynamic response of several asteroids and the change in rotation rate for specific shapes as a function of obliquity. Finally, we define a convenient dimensionless parameter that is only a function of the body geometry and that can be used to characterize the effects of YORP.     

Long term evolution of comet Halleys orbit

The aim of the paper is to study the long term evolution of comet Halleys orbit taking into account small errors in the initial conditions. Recent papers deal with mapping methods to model cometary dynamics; (e.g. Petrosky and Broucke, 1987 and Chirikov and Vecheslavov, 1986). They will be discussed critically and compared with our own results. We then tested the model using numerical integration methods. For the moment we limited our calculation to 2.10⁵ years, but a 10⁶ year integration is still in progress. We show the expected dynamical evolution of Hallyes orbit taking into account also smaller and larger errors of the initial conditions (nongravitational effects are only roughly estimated). Finally we discuss alsothe controversal opinions concerning the role of the planets (especially the earth).