Modelling the secular evolution of migrating planet pairs

The subject of this paper is the secular behaviour of a pair of planets evolving under dissipative forces. In particular, we investigate the case when dissipative forces affect the planetary semimajor axes and the planets move inwards/outwards the central star, in a process known as planet migration. To perform this investigation, we introduce fundamental concepts of conservative and dissipative dynamics of the three-body problem. Based on these concepts, we develop a qualitative model of the secular evolution of the migrating planetary pair. Our approach is based on the analysis of the energy and the orbital angular momentum exchange between the two-planet system and an external medium; thus no specific kind of dissipative forces is invoked. We show that, under the assumption that dissipation is weak and slow, the evolutionary routes of the migrating planets are traced by the Mode I and Mode II stationary solutions of the conservative secular problem. The ultimate convergence and the evolution of the system along one of these secular modes of motion are determined uniquely by the condition that the dissipation rate is sufficiently smaller than the proper secular frequency of the system. We show that it is possible to reassemble the starting configurations and the migration history of the systems on the basis of their final states and consequently to constrain the parameters of the physical processes involved.

Hydrodynamical simulations of Galactic fountains - II. Evolution of multiple fountains

The ejection of gas out of the disc in late-type galaxies is related to star formation and is mainly due to the explosion of Type II supernovae (SN II). In a previous paper, we considered the evolution of a single Galactic fountain, that is, a fountain powered by a single SN cluster. Using three-dimensional hydrodynamical simulations, we studied in detail the fountain flow and its dependence with several factors, such as the Galactic rotation, the distance to the Galactic centre and the presence of a hot gaseous halo. As a natural followup, this paper investigates the dynamical evolution of multiple generations of fountains generated by similar to 100 OB associations. We have considered the observed size-frequency distribution of young stellar clusters within the Galaxy in order to appropriately fuel the multiple fountains in our simulations. Most of the results of the previous paper have been confirmed, like for example the formation of intermediate velocity clouds above the disc by the multiple fountains. Also, this work confirms the localized nature of the fountain flows: the freshly ejected metals tend to fall back close to the same Galactocentric region where they are delivered. Therefore, the fountains do not change significantly the radial profile of the disc chemical abundance. The multiple fountain simulations also allowed us to consistently calculate the feedback of the star formation on the halo gas. We found that the hot gas gains about 10 per cent of all the SN II energy produced in the disc. Thus, the SN feedback more than compensate for the halo radiative losses and allow for a quasi steady-state disc-halo circulation to exist. Finally, we have also considered the possibility of mass infall from the intergalactic medium and its interaction with the clouds that are formed by the fountains. Though our simulations are not suitable to reproduce the slow rotational pattern that is typically observed in the haloes around the disc galaxies, they indicate that the presence of an external gas infall may help to slow down the rotation of the gas in the clouds and thus the amount of angular momentum that they transfer to the coronal gas, as previously suggested in the literature.

Dynamical erosion of asteroid groups in the region of the Phocaea family

In a previous paper, the current state of knowledge of the region containing the Phocaea dynamical family was revised. Here, the dynamical evolution and possible origin of the Phocaea dynamical family and asteroid groups in the region are investigated. First, I study the case of asteroids at high eccentricity (e > 0.31). I find that these objects are unstable because of encounters with Mars on time-scales of up to 270 Myr. The minimum time needed by members of the Phocaea classical family to reach the orbital locations of these objects, 370 Myr, can be used to set a lower limit on the age of the Phocaea family.Next, attention is focused on the chaotic layer previously identified near the nu(6) secular resonance border. Using analytical and numerical tools, I find that the presence of the nu(6) secular resonance forces asteroids with vertical bar g-g(6)vertical bar < 2.55 arcsec yr(-1) to reach eccentricities high enough to allow them to experience deep, close encounters with Mars. Results of the analytical model of Yoshikawa and of my numerical simulations fully explain the low-inclination chaotic region found by Carruba.Finally, I investigate the long-term stability of the minor families and clumps identified in the previous paper, with particular emphasis on a clump only identifiable in the domain of proper frequencies (n, g, g - s) around (6246) Komurotoru. I find that while the clumps identified in the space of proper elements quickly disperse when the Yarkovsky effect is considered, the family around (19536) is still observable for time-scales of more than 50 Myr. The (6246) clump, characterized by its interaction with the nu(5) + nu(16) and 2 nu(6) - nu(16) secular resonances, is robust on time-scales of 50 Myr. I confirm that this group may be the first clump ever detected in the frequency domain that can be associated with a real collisional event.

Dynamical erosion of asteroid groups in the region of the Pallas family

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Theoretical models of planetary system formation II. Post-formation evolution

Aims. We extend the results of planetary formation synthesis by computing the long-term evolution of synthetic systems from the clearing of the gas disk into the dynamical evolution phase. Methods. We use the symplectic integrator SyMBA to numerically integrate the orbits of planets for 100 Myr, using populations from previous studies as initial conditions. Results. We show that within the populations studied, mass and semimajor axis distributions experience only minor changes from post-formation evolution. We also show that, depending upon their initial distribution, planetary eccentricities can statistically increase or decrease as a result of gravitational interactions. We find that planetary masses and orbital spacings provided by planet formation models do not result in eccentricity distributions comparable to observed exoplanet eccentricities, requiring other phenomena, such as stellar fly-bys, to account for observed eccentricities.

On gravitational fields of stars in f(R) gravity theories

Einstein's general relativity is a classical theory of gravitation: it is a postulate on the coupling between the four-dimensional, continuos spacetime and the matter fields in the universe, and it yields their dynamical evolution. It is believed that general relativity must be replaced by a quantum theory of gravity at least at extremely high energies of the early universe and at regions of strong curvature of spacetime, cf. black holes. Various attempts to quantize gravity, including conceptually new models such as string theory, have suggested that modification to general relativity might show up even at lower energy scales. On the other hand, also the late time acceleration of the expansion of the universe, known as the dark energy problem, might originate from new gravitational physics. Thus, although there has been no direct experimental evidence contradicting general relativity so far - on the contrary, it has passed a variety of observational tests - it is a question worth asking, why should the effective theory of gravity be of the exact form of general relativity? If general relativity is modified, how do the predictions of the theory change? Furthermore, how far can we go with the changes before we are face with contradictions with the experiments? Along with the changes, could there be new phenomena, which we could measure to find hints of the form of the quantum theory of gravity? This thesis is on a class of modified gravity theories called f(R) models, and in particular on the effects of changing the theory of gravity on stellar solutions. It is discussed how experimental constraints from the measurements in the Solar System restrict the form of f(R) theories. Moreover, it is shown that models, which do not differ from general relativity at the weak field scale of the Solar System, can produce very different predictions for dense stars like neutron stars. Due to the nature of f(R) models, the role of independent connection of the spacetime is emphasized throughout the thesis.

Influence of viscoelastic nature on the intermittent peel-front dynamics of adhesive tape

We investigate the influence of viscoelastic nature of the adhesive on the intermittent peel front dynamics by extending a recently introduced model for peeling of an adhesive tape. As time and rate-dependent deformation of the adhesives are measured in stationary conditions, a crucial step in incorporating the viscoelastic effects applicable to unstable intermittent peel dynamics is the introduction of a dynamization scheme that eliminates the explicit time dependence in terms of dynamical variables. We find contrasting influences of viscoelastic contribution in different regions of tape mass, roller inertia, and pull velocity. As the model acoustic energy dissipated depends on the nature of the peel front and its dynamical evolution, the combined effect of the roller inertia and pull velocity makes the acoustic energy noisier for small tape mass and low-pull velocity while it is burstlike for low-tape mass, intermediate values of the roller inertia and high-pull velocity. The changes are quantified by calculating the largest Lyapunov exponent and analyzing the statistical distributions of the amplitudes and durations of the model acoustic energy signals. Both single and two stage power-law distributions are observed. Scaling relations between the exponents are derived which show that the exponents corresponding to large values of event sizes and durations are completely determined by those for small values. Th scaling relations are found to be satisfied in all cases studied. Interestingly, we find only five types of model acoustic emission signals among multitude of possibilities of the peel front configurations.

Lopsidedness and sloshing in centres of advanced mergers of galaxies

We measure the non-axisymmetry in the luminosity distribution in the central few kpc of a sample of advanced mergers of galaxies, by analyzing their 2MASS images. All mergers show a high central asymmetry: the centres of isophotes show a striking sloshing pattern with a spatial variation of upto 30% within the central 1 kpc; and the Fourier amplitude for lopsidedness (m = 1) shows high values upto 0.2 within the central 5 kpc. The central asymmetry is estimated to be long-lived, lasting for ~ a few Gyr or ~ 100 local dynamical timescales. This will significantly affect the dynamical evolution of this region, by helping fuel the central active galactic nucleus, and also by causing the secular growth of the bulge driven by lopsidedness.

Globally coupled stochastic two-state oscillators: Fluctuations due to finite numbers

Infinite arrays of coupled two-state stochastic oscillators exhibit well-defined steady states. We study the fluctuations that occur when the number N of oscillators in the array is finite. We choose a particular form of global coupling that in the infinite array leads to a pitchfork bifurcation from a monostable to a bistable steady state, the latter with two equally probable stationary states. The control parameter for this bifurcation is the coupling strength. In finite arrays these states become metastable: The fluctuations lead to distributions around the most probable states, with one maximum in the monostable regime and two maxima in the bistable regime. In the latter regime, the fluctuations lead to transitions between the two peak regions of the distribution. Also, we find that the fluctuations break the symmetry in the bimodal regime, that is, one metastable state becomes more probable than the other, increasingly so with increasing array size. To arrive at these results, we start from microscopic dynamical evolution equations from which we derive a Langevin equation that exhibits an interesting multiplicative noise structure. We also present a master equation description of the dynamics. Both of these equations lead to the same Fokker-Planck equation, the master equation via a 1/N expansion and the Langevin equation via standard methods of Ito calculus for multiplicative noise. From the Fokker-Planck equation we obtain an effective potential that reflects the transition from the monomodal to the bimodal distribution as a function of a control parameter. We present a variety of numerical and analytic results that illustrate the strong effects of the fluctuations. We also show that the limits N -> infinity and t -> infinity(t is the time) do not commute. In fact, the two orders of implementation lead to drastically different results.

An evolutionary formalism for weak quantum measurements

Unitary evolution and projective measurement are fundamental axioms of quantum mechanics. Even though projective measurement yields one of the eigenstates of the measured operator as the outcome, there is no theory that predicts which eigenstate will be observed in which experimental run. There exists only an ensemble description, which predicts probabilities of various outcomes over many experimental runs. We propose a dynamical evolution equation for the projective collapse of the quantum state in individual experimental runs, which is consistent with the well-established framework of quantum mechanics. In case of gradual weak measurements, its predictions for ensemble evolution are different from those of the Born rule. It is an open question whether or not suitably designed experiments can observe this alternate evolution.

Sample-specific behavior in failure models of disordered medial

The concept ''sample-specific'' is suggested to describe the behavior of disordered media close to macroscopic failure. it is pointed out that the transition from universal scaling to sample-specific behavior may be a common phenomenon in failure models of disordered media. The dynamical evolution plays an important role in the transition.

The mechanism of energy storage by shearing the solar magnetic field

In this paper, a complete set of MHD equations have been solved by numerical calculations in an attempt to study the dynamical evolutionary processes of the initial equilibrium configuration and to discuss the energy storage mechanism of the solar atmosphere by shearing the magnetic field. The initial equilibrium configuration with an arch bipolar potential field obtained from the numerical solution is similar to the configuration in the vicinity of typical solar flare before its eruption. From the magnetic induction equation in the set of MHD equations and dealing with the non-linear coupling effects between the flow field and magnetic field, the quantitative relationship has been derived for their dynamical evolution. Results show that plasma shear motion at the bottom of the solar atmosphere causes the magnetic field to shear; meanwhile the magnetic field energy is stored in local regions. With the increase of time the local magnetic energy increases and it may reach an order of 4×10^25 J during a day. Thus the local storage of magnetic energy is large enough to trigger a big solar flare and can be considered as the energy source of solar flares. The energy storage mechanism by shearing the magnetic field can well explain the slow changes in solar active regions.

Time-Domain TeraHertz Spectroscopy and Observational Probes of Prebiotic Interstellar Gas and Ice Chemistry

Understanding the origin of life on Earth has long fascinated the minds of the global community, and has been a driving factor in interdisciplinary research for centuries. Beyond the pioneering work of Darwin, perhaps the most widely known study in the last century is that of Miller and Urey, who examined the possibility of the formation of prebiotic chemical precursors on the primordial Earth [1]. More recent studies have shown that amino acids, the chemical building blocks of the biopolymers that comprise life as we know it on Earth, are present in meteoritic samples, and that the molecules extracted from the meteorites display isotopic signatures indicative of an extraterrestrial origin [2]. The most recent major discovery in this area has been the detection of glycine (NH₂CH₂COOH), the simplest amino acid, in pristine cometary samples returned by the NASA STARDUST mission [3]. Indeed, the open questions left by these discoveries, both in the public and scientific communities, hold such fascination that NASA has designated the understanding of our "Cosmic Origins" as a key mission priority.

Despite these exciting discoveries, our understanding of the chemical and physical pathways to the formation of prebiotic molecules is woefully incomplete. This is largely because we do not yet fully understand how the interplay between grain-surface and sub-surface ice reactions and the gas-phase affects astrophysical chemical evolution, and our knowledge of chemical inventories in these regions is incomplete. The research presented here aims to directly address both these issues, so that future work to understand the formation of prebiotic molecules has a solid foundation from which to work.

From an observational standpoint, a dedicated campaign to identify hydroxylamine (NH₂OH), potentially a direct precursor to glycine, in the gas-phase was undertaken. No trace of NH₂OH was found. These observations motivated a refinement of the chemical models of glycine formation, and have largely ruled out a gas-phase route to the synthesis of the simplest amino acid in the ISM. A molecular mystery in the case of the carrier of a series of transitions was resolved using observational data toward a large number of sources, confirming the identity of this important carbon-chemistry intermediate B11244 as l-C₃H⁺ and identifying it in at least two new environments. Finally, the doubly-nitrogenated molecule carbodiimide HNCNH was identified in the ISM for the first time through maser emission features in the centimeter-wavelength regime.

In the laboratory, a TeraHertz Time-Domain Spectrometer was constructed to obtain the experimental spectra necessary to search for solid-phase species in the ISM in the THz region of the spectrum. These investigations have shown a striking dependence on large-scale, long-range (i.e. lattice) structure of the ices on the spectra they present in the THz. A database of molecular spectra has been started, and both the simplest and most abundant ice species, which have already been identified, as well as a number of more complex species, have been studied. The exquisite sensitivity of the THz spectra to both the structure and thermal history of these ices may lead to better probes of complex chemical and dynamical evolution in interstellar environments.

Optical colors of 56 near-Earth objects: trends with size and orbit

We present the results of BVRIZ photometry of 56 near-Earth objects (NEOs) obtained with the 1-m Jacobus Kapteyn telescope on La Palma during 2000 and 2001. Our sample includes many NEOs with particularly deep 1 - mum pyroxene/olivine absorption bands, similar to Q-type asteroids. We also classify three NEOs with particularly blue colors. No D-type asteroids were found, placing an upper limit of similar to2% on the fraction of the NEO population originating in the outer main belt or the Trojan clouds. The ratio of dark to bright objects in our sample was found to be 0.40, significantly higher than current theoretical predictions. As well as classifying the NEOs, we have investigated color trends with size and orbit. We see a general trend for larger silicate objects to have shallower absorption bands but find no significant difference in the distribution of taxonomic classes at small and large sizes. Our data clearly show that different taxonomic classes tend to occupy different regions of (a, e) space. By comparing our data with current model predictions for NEO dynamical evolution we see that Q- R-, and V-type NEOs tend to have orbits associated with "fast track" delivery from the main belt, whereas S-type NEOs tend to have orbits associated with "slow track" delivery. This outcome would be expected if space weathering occurs on time scales of >10(6) years. (C) 2003 Elsevier Science (USA). All rights reserved.

CH abundance gradient in TMC-1

Aims. The aim of this study is to examine if the well-known chemical gradient in TMC-1 is reflected in the amount of rudimentary forms of carbon available in the gas-phase. As a tracer we use the CH radical which is supposed to be well correlated with carbon atoms and simple hydrocarbon ions. Methods. We observed the 9-cm ?-doubling lines of CH along the dense filament of TMC-1. The CH column densities were compared with the total H2 column densities derived using the 2MASS NIR data and previously published SCUBA maps and with OH column densities derived using previous observations with Effelsberg. We also modelled the chemical evolution of TMC-1 adopting physical conditions typical of dark clouds using the UMIST Database for Astrochemistry gas-phase reaction network to aid the interpretation of the observed OH/CH abundance ratios. Results. The CH column density has a clear peak in the vicinity of the cyanopolyyne maximum of TMC-1. The fractional CH abundance relative to H2 increases steadily from the northwestern end of the filament where it lies around 1.0 × 10-8 , to the southeast where it reaches a value of 2.0 × 10-8. The OH and CH column densities are well correlated, and we obtained OH/CH abundance ratios of ~16–20. These values are clearly larger than what has been measured recently in diffuse interstellar gas and is likely to be related to C to CO conversion at higher densities. The good correlation between CH and OH can be explained by similar production and destruction pathways. We suggest that the observed CH and OH abundance gradients are mainly due to enhanced abundances in a low-density envelope which becomes more prominent in the southeastern part and seems to continue beyond the dense filament. Conclusions. An extensive envelope probably signifies an early stage of dynamical evolution, and conforms with the detection of a large CH abundance in the southeastern part of the cloud. The implied presence of other simple forms of carbon in the gas phase provides a natural explanation for the observation of “early-type” molecules in this region.