Global lopsided instability in a purely stellar galactic disc

It is shown that pure exponential discs in spiral galaxies are capable of supporting slowly varying discrete global lopsided modes, which can explain the observed features of lopsidedness in the stellar discs. Using linearized fluid dynamical equations with the softened self-gravity and pressure of the perturbation as the collective effect, we derive self-consistently a quadratic eigenvalue equation for the lopsided perturbation in the galactic disc. On solving this, we find that the ground-state mode shows the observed characteristics of the lopsidedness in a galactic disc, namely the fractional Fourier amplitude A(1), increases smoothly with the radius. These lopsided patterns precess in the disc with a very slow pattern speed with no preferred sense of precession. We show that the lopsided modes in the stellar disc are long-lived because of a substantial reduction (approximately a factor of 10 compared to the local free precession rate) in the differential precession. The numerical solution of the equations shows that the groundstate lopsided modes are either very slowly precessing stationary normal mode oscillations of the disc or growing modes with a slow growth rate depending on the relative importance of the collective effect of the self-gravity. N-body simulations are performed to test the spontaneous growth of lopsidedness in a pure stellar disc. Both approaches are then compared and interpreted in terms of long-lived global m = 1 instabilities, with almost zero pattern speed.

Acceleration of the Cosmological Expansion as an Effect of Inhomogeneities

The cosmological observations of light from type Ia supernovae, the cosmic microwave background and the galaxy distribution seem to indicate that the expansion of the universe has accelerated during the latter half of its age. Within standard cosmology, this is ascribed to dark energy, a uniform fluid with large negative pressure that gives rise to repulsive gravity but also entails serious theoretical problems. Understanding the physical origin of the perceived accelerated expansion has been described as one of the greatest challenges in theoretical physics today. In this thesis, we discuss the possibility that, instead of dark energy, the acceleration would be caused by an effect of the nonlinear structure formation on light, ignored in the standard cosmology. A physical interpretation of the effect goes as follows: due to the clustering of the initially smooth matter with time as filaments of opaque galaxies, the regions where the detectable light travels get emptier and emptier relative to the average. As the developing voids begin to expand the faster the lower their matter density becomes, the expansion can then accelerate along our line of sight without local acceleration, potentially obviating the need for the mysterious dark energy. In addition to offering a natural physical interpretation to the acceleration, we have further shown that an inhomogeneous model is able to match the main cosmological observations without dark energy, resulting in a concordant picture of the universe with 90% dark matter, 10% baryonic matter and 15 billion years as the age of the universe. The model also provides a smart solution to the coincidence problem: if induced by the voids, the onset of the perceived acceleration naturally coincides with the formation of the voids. Additional future tests include quantitative predictions for angular deviations and a theoretical derivation of the model to reduce the required phenomenology. A spin-off of the research is a physical classification of the cosmic inhomogeneities according to how they could induce accelerated expansion along our line of sight. We have identified three physically distinct mechanisms: global acceleration due to spatial variations in the expansion rate, faster local expansion rate due to a large local void and biased light propagation through voids that expand faster than the average. A general conclusion is that the physical properties crucial to account for the perceived acceleration are the growth of the inhomogeneities and the inhomogeneities in the expansion rate. The existence of these properties in the real universe is supported by both observational data and theoretical calculations. However, better data and more sophisticated theoretical models are required to vindicate or disprove the conjecture that the inhomogeneities are responsible for the acceleration.

Non-Gaussianities and Preheating from N-flation and Magnetogenesis via rotating cosmic string loops

In this thesis we examine multi-field inflationary models of the early Universe. Since non-Gaussianities may allow for the possibility to discriminate between models of inflation, we compute deviations from a Gaussian spectrum of primordial perturbations by extending the delta-N formalism. We use N-flation as a concrete model; our findings show that these models are generically indistinguishable as long as the slow roll approximation is still valid. Besides computing non-Guassinities, we also investigate Preheating after multi-field inflation. Within the framework of N-flation, we find that preheating via parametric resonance is suppressed, an indication that it is the old theory of preheating that is applicable. In addition to studying non-Gaussianities and preheatng in multi-field inflationary models, we study magnetogenesis in the early universe. To this aim, we propose a mechanism to generate primordial magnetic fields via rotating cosmic string loops. Magnetic fields in the micro-Gauss range have been observed in galaxies and clusters, but their origin has remained elusive. We consider a network of strings and find that rotating cosmic string loops, which are continuously produced in such networks, are viable candidates for magnetogenesis with relevant strength and length scales, provided we use a high string tension and an efficient dynamo.

Studies of the star-forming structures in the dense interstellar medium: a view by dust extinction

New stars in galaxies form in dense, molecular clouds of the interstellar medium. Measuring how the mass is distributed in these clouds is of crucial importance for the current theories of star formation. This is because several open issues in them, such as the strength of different mechanism regulating star formation and the origin of stellar masses, can be addressed using detailed information on the cloud structure. Unfortunately, quantifying the mass distribution in molecular clouds accurately over a wide spatial and dynamical range is a fundamental problem in the modern astrophysics. This thesis presents studies examining the structure of dense molecular clouds and the distribution of mass in them, with the emphasis on nearby clouds that are sites of low-mass star formation. In particular, this thesis concentrates on investigating the mass distributions using the near infrared dust extinction mapping technique. In this technique, the gas column densities towards molecular clouds are determined by examining radiation from the stars that shine through the clouds. In addition, the thesis examines the feasibility of using a similar technique to derive the masses of molecular clouds in nearby external galaxies. The papers presented in this thesis demonstrate how the near infrared dust extinction mapping technique can be used to extract detailed information on the mass distribution in nearby molecular clouds. Furthermore, such information is used to examine characteristics crucial for the star formation in the clouds. Regarding the use of extinction mapping technique in nearby galaxies, the papers of this thesis show that deriving the masses of molecular clouds using the technique suffers from strong biases. However, it is shown that some structural properties can still be examined with the technique.

Formation of structure in dark energy cosmologies

Acceleration of the universe has been established but not explained. During the past few years precise cosmological experiments have confirmed the standard big bang scenario of a flat universe undergoing an inflationary expansion in its earliest stages, where the perturbations are generated that eventually form into galaxies and other structure in matter, most of which is non-baryonic dark matter. Curiously, the universe has presently entered into another period of acceleration. Such a result is inferred from observations of extra-galactic supernovae and is independently supported by the cosmic microwave background radiation and large scale structure data. It seems there is a positive cosmological constant speeding up the universal expansion of space. Then the vacuum energy density the constant describes should be about a dozen times the present energy density in visible matter, but particle physics scales are enormously larger than that. This is the cosmological constant problem, perhaps the greatest mystery of contemporary cosmology. In this thesis we will explore alternative agents of the acceleration. Generically, such are called dark energy. If some symmetry turns off vacuum energy, its value is not a problem but one needs some dark energy. Such could be a scalar field dynamically evolving in its potential, or some other exotic constituent exhibiting negative pressure. Another option is to assume that gravity at cosmological scales is not well described by general relativity. In a modified theory of gravity one might find the expansion rate increasing in a universe filled by just dark matter and baryons. Such possibilities are taken here under investigation. The main goal is to uncover observational consequences of different models of dark energy, the emphasis being on their implications for the formation of large-scale structure of the universe. Possible properties of dark energy are investigated using phenomenological paramaterizations, but several specific models are also considered in detail. Difficulties in unifying dark matter and dark energy into a single concept are pointed out. Considerable attention is on modifications of gravity resulting in second order field equations. It is shown that in a general class of such models the viable ones represent effectively the cosmological constant, while from another class one might find interesting modifications of the standard cosmological scenario yet allowed by observations. The thesis consists of seven research papers preceded by an introductory discussion.

Possible use of self-calibration to reduce systematic uncertainties in determining distance-redshift relation via gravitational radiation from merging binaries

By observing mergers of compact objects, future gravity wave experiments would measure the luminosity distance to a large number of sources to a high precision but not their redshifts. Given the directional sensitivity of an experiment, a fraction of such sources (gold plated) can be identified optically as single objects in the direction of the source. We show that if an approximate distance-redshift relation is known then it is possible to statistically resolve those sources that have multiple galaxies in the beam. We study the feasibility of using gold plated sources to iteratively resolve the unresolved sources, obtain the self-calibrated best possible distance-redshift relation and provide an analytical expression for the accuracy achievable. We derive the lower limit on the total number of sources that is needed to achieve this accuracy through self-calibration. We show that this limit depends exponentially on the beam width and give estimates for various experimental parameters representative of future gravitational wave experiments DECIGO and BBO.

The slowing down of galaxy disks in dissipationless minor mergers

We have investigated the impact of dissipationless minor galaxy mergers on the angular momentum of the remnant. Our simulations cover a range of initial orbital characteristics, and the system consists of a massive galaxy with a bulge and disk merging with a much less massive (one-tenth or one-twentieth) gasless companion that has a variety of morphologies (disk-or elliptical-like) and central baryonic mass concentrations. During the process of merging, the orbital angular momentum is redistributed into the internal angular momentum of the final system; the internal angular momentum of the primary galaxy can increase or decrease depending on the relative orientation of the orbital spin vectors (direct or retrograde), while the initially nonrotating dark matter halo always gains angular momentum. The specific angular momentum of the stellar component always decreases independently of the orbital parameters or morphology of the satellite, the decrease in the rotation velocity of the primary galaxy is accompanied by a change in the anisotropy of the orbits, and the ratio of rotation speed to velocity dispersion of the merger remnant is lower than the initial value, not only because of an increase in the dispersion but also of the slowing-down of the disk rotation. We briefly discuss several astrophysical implications of these results, suggesting that minor mergers do not cause a "random walk" process of the angular momentum of the stellar disk component of galaxies, but rather a steady decrease. Minor mergers may play a role in producing the large scatter observed in the Tully-Fisher relation for S0 galaxies, as well as in the increase of the velocity dispersion and the decrease in upsilon/sigma at large radii as observed in S0 galaxies.

Undamped Oscillations of Collisionless Stellar Systems: Spheres,Spheroids and Discs.

The collisionless Boltzmann equation governing self-gravitating systems such as galaxies has recently been shown to admit exact oscillating solutions with planar and spherical symmetry. The relation of the spherically symmetric solutions to the Virial theorem, as well as generalizations to non-uniform spheres, uniform spheroids and discs form the subject of this paper. These models generalize known families of static solutions. The case of the spheroid is worked out in some detail. Quasiperiodic as well as chaotic time variation of the two axes is demonstrated by studying the surface of section for the associated Hamiltonian system with two degrees of freedom. The relation to earlier work and possible implications for the general problem of collisionless relaxation in self gravitating systems are also discussed.

Effect of a QSO on the interstellar medium of its host galaxy

We consider the hydrodynamic evolution of gas in the interstellar medium of the host galaxy of a quasar due to Compton heating by the QSO radiation. We show that a Lagrangean formulation of the problem is necessary. It is found that the "hydrodynamic time scale" becomes important compared to the Compton heating time scale. We also relax the "single fluid" approximation by considering the existence of clouds and taking into account the mass loss from stars. The results predict star burst activity, and thus we explain the blue colors of the active galaxies.

Embracing Integrability in Stellar and Planetary Dynamics

Hamiltonian systems in stellar and planetary dynamics are typically near integrable. For example, Solar System planets are almost in two-body orbits, and in simulations of the Galaxy, the orbits of stars seem regular. For such systems, sophisticated numerical methods can be developed through integrable approximations. Following this theme, we discuss three distinct problems. We start by considering numerical integration techniques for planetary systems. Perturbation methods (that utilize the integrability of the two-body motion) are preferred over conventional "blind" integration schemes. We introduce perturbation methods formulated with Cartesian variables. In our numerical comparisons, these are superior to their conventional counterparts, but, by definition, lack the energy-preserving properties of symplectic integrators. However, they are exceptionally well suited for relatively short-term integrations in which moderately high positional accuracy is required. The next exercise falls into the category of stability questions in solar systems. Traditionally, the interest has been on the orbital stability of planets, which have been quantified, e.g., by Liapunov exponents. We offer a complementary aspect by considering the protective effect that massive gas giants, like Jupiter, can offer to Earth-like planets inside the habitable zone of a planetary system. Our method produces a single quantity, called the escape rate, which characterizes the system of giant planets. We obtain some interesting results by computing escape rates for the Solar System. Galaxy modelling is our third and final topic. Because of the sheer number of stars (about 10^11 in Milky Way) galaxies are often modelled as smooth potentials hosting distributions of stars. Unfortunately, only a handful of suitable potentials are integrable (harmonic oscillator, isochrone and Stäckel potential). This severely limits the possibilities of finding an integrable approximation for an observed galaxy. A solution to this problem is torus construction; a method for numerically creating a foliation of invariant phase-space tori corresponding to a given target Hamiltonian. Canonically, the invariant tori are constructed by deforming the tori of some existing integrable toy Hamiltonian. Our contribution is to demonstrate how this can be accomplished by using a Stäckel toy Hamiltonian in ellipsoidal coordinates.

Compressive tidal heating of a disk galaxy in a rich cluster

The restricted three-body method is used to model the effect of the mean tidal field of a cluster of galaxies on the internal dynamics of a disk galaxy falling into the cluster for the first time. In the model adopted the galaxy experiences a tidal field that is compressive within the core of the cluster. The planar random velocities of all components in the disk increase after the galaxy passes through the core of the cluster. The low-velocity dispersion gas clouds experience a relatively larger increase in random velocity than the hotter stellar components. The increase in planar velocities results in a strong anisotropy between the planar and vertical velocity dispersions. It is argued that this will make the disk unstable to the 'fire-hose instability' which leads to bending modes in the disk and which will thicken the disk slightly. The mean tidal fields in rich clusters were probably stronger during the epoch of cluster formation and relaxation than they are in present-day relaxed clusters.

Results from PAMELA, ATIC and FERMI: Pulsars or dark matter?

It is well known that dark matter dominates the dynamics of galaxies and clusters of galaxies. Its constituents remain a mystery despite an assiduous search for them over the past three decades. Recent results from the satellite-based PAMELA experiment show an excess in the positron fraction at energies between 10 and 100 GeV in the secondary cosmic ray spectrum. Other experiments, namely ATIC, HESS and FERMI, show an excess in the total electron (e(+) + e(-)) spectrum for energies greater than 100 GeV. These excesses in the positron fraction as well as the electron spectrum can arise in local astrophysical processes like pulsars, or can be attributed to the annihilation of the dark matter particles. The latter possibility gives clues to the possible candidates for the dark matter in galaxies and other astrophysical systems. In this article, we give a report of these exciting developments.

Flat rotation curves: A result of the helical inverse cascade in turbulent media

We have modeled the rotation curves of 21 galaxies observed by Amram et al. (1992), by combining the effects of rigid rotation, gravity, and turbulence. The main motivation behind such modeling is to study the formation of coherent structures in turbulent media and explore its role in the large-scale structures of the universe. The values of the parameters such as mass, turbulent velocity, and angular velocity derived from the rotation curve fits are in good agreement with those derived from the prevalent models.

Star formation in giant extragalactic H II regions

Star formation properties in Giant Extragalactic H II Regions (GEHRs) are investigated using optical photometry and evolutionary population synthesis models. Photometric data in $BVR$ bands and in the emission line of H-alpha are obtained by CCD imaging at Vainu Bappu Observatory, Kavalur. Aperture photometry is performed for 180 GEHRs in galaxies NGC 1365, 1566, 2366, 2903, 2997, 3351, 4303, 4449, 4656 and 5253. Thirty six of these GEHRs having published spectroscopic data are studied for star formation properties. The population synthesis model is constructed based on Maeder's stellar evolutionary and Kurucz stellar atmosphere models, to synthesize observational quantities of embedded clusters in GEHRs. The observed H-alpha luminosity is a measure of the number of massive stars while the contribution to BVR bands is from intermediate mass (5-15 solar mass) stars when the cluster is young and from evolving supergiants when the cluster is old (age >/= 6~Myr). Differential reddening between gas and embedded stars is essential to constrain the dereddened cluster colors within the range of youngest clusters. Obscuring dust closely associated with gas, which is distributed in filaments and clumps, as in the case of 30 Doradus, is the most likely configuration giving rise to net reduction of extinction towards stars. The fraction of the stellar photons escaping the nebula unattenuated is estimated to be 50%. GEHRs are rarely found to be simple systems containing stars from single generation. In the present sample such regions in addition to being older than 3~Myr, have their Lyman continuum luminosity reduced by as much as 60%, compared to the observed $B$ band luminosity for a normal IMF. The missing ionizing photons may be escaping the nebula, leading to the ionization of extra-H II region ionized medium. Co-existence of young (age