A numerical survey of relativistic rotating neutron star structures using the Hartle-Thorne formalism

A broad numerical survey of relativistic rotating neutron star structures was compiled using an exhaustive list of presently available equation of state models for neutron star matter. The structure parameters (spherical deformations in mass and radii, the moment of inertia and quadrupole moment, oblateness, and free precession) are calculated using the formalism proposed by Hartle and Thorne (1968). The results are discussed in relation to the relevant observational information. Binary pulsar data and X-ray burst sources provide information on the bulk properties of neutron stars, enabling the derivation of constraints that can be put on the structure of neutron stars and equation of state models.

A Triggering Mechanism for Enhanced Star-Formation in Colliding Galaxies

We propose a physical mechanism to explain the origin of the intense burst of massive-star formation seen in colliding/merging, gas-rich, field spiral galaxies. We explicitly take account of the different parameters for the two main mass components, H-2 and H I, of the interstellar medium within a galaxy and follow their consequent different evolution during a collision between two galaxies. We also note that, in a typical spiral galaxy-like our galaxy, the Giant Molecular Clouds (GMCs) are in a near-virial equilibrium and form the current sites of massive-star formation, but have a low star formation rate. We show that this star formation rate is increased following a collision between galaxies. During a typical collision between two field spiral galaxies, the H I clouds from the two galaxies undergo collisions at a relative velocity of approximately 300 km s-1. However, the GMCs, with their smaller volume filling factor, do not collide. The collisions among the H I clouds from the two galaxies lead to the formation of a hot, ionized, high-pressure remnant gas. The over-pressure due to this hot gas causes a radiative shock compression of the outer layers of a preexisting GMC in the overlapping wedge region. This makes these layers gravitationally unstable, thus triggering a burst of massive-star formation in the initially barely stable GMCs.The resulting value of the typical IR luminosity from the young, massive stars from a pair of colliding galaxies is estimated to be approximately 2 x 10(11) L., in agreement with the observed values. In our model, the massive-star formation occurs in situ in the overlapping regions of a pair of colliding galaxies. We can thus explain the origin of enhanced star formation over an extended, central area approximately several kiloparsecs in size, as seen in typical colliding galaxies, and also the origin of starbursts in extranuclear regions of disk overlap as seen in Arp 299 (NGC 3690/IC 694) and in Arp 244 (NGC 4038/39). Whether the IR emission from the central region or that from the surrounding extranuclear galactic disk dominates depends on the geometry and the epoch of the collision and on the initial radial gas distribution in the two galaxies. In general, the central starburst would be stronger than that in the disks, due to the higher preexisting gas densities in the central region. The burst of star formation is expected to last over a galactic gas disk crossing time approximately 4 x 10(7) yr. We can also explain the simultaneous existence of nearly normal CO galaxy luminosities and shocked H-2 gas, as seen in colliding field galaxies.This is a minimal model, in that the only necessary condition for it to work is that there should be a sufficient overlap between the spatial gas distributions of the colliding galaxy pair.

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.

Photometric calibration of the CCD camera of 1-m telescope at VBO

Calibration of the CCD camera of the 1-m telescope at the Vainu Bappu Observatory, Kavalur, to the BVR system is reported here based on the observations of stars in the 'dipper asterism' in the open cluster M 67 (NGC 2682). Transformations involving B and V have negligible colour terms, while those involving R are slightly colour dependent. The possibility of using scale-down R band fluxes to estimate the continuum flux at H-alpha is investigated by comparing the counts in R band with those through an interference filter centred at H-alpha. The scaling factor is found to remain constant over a wide range of colours. The sensitivity of the telescope-filter-CCD combination is estimated to be 2.0 per cent, 8.3 per cent and 9.7 per cent in B, V and R bands, respectively. The star F117 appears to be a small-amplitude (approximately 0.05 mag) variable.

Focal plane analysis of the Knox-Thompson algorithm in speckle interferometry

For the specific case of binary stars, this paper presents signal-to-noise ratio (SNR) calculations for the detection of the parity (the side of the brighter component) of the binary using the double correlation method. This double correlation method is a focal plane version of the well-known Knox-Thompson method used in speckle interferometry. It is shown that SNR for parity detection using double correlation depends linearly on binary separation. This new result was entirely missed by previous analytical calculations dealing with a point source. It is concluded that, for magnitudes relevant to the present day speckle interferometry and for binary separations close to the diffraction limit, speckle masking has better SNR for parity detection.

Reynolds stress and heat flux in spherical shell convection

Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Aims. We compute turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime in order to study whether restricted geometry introduces artefacts into the results. Methods. We employ direct numerical simulations of turbulent convection in spherical and Cartesian geometries. In order to alleviate the computational cost in the spherical runs and to reach as high spatial resolution as possible, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number or inverse Rossby number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in accordance with earlier numerical results, but in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong 'banana cells'. Their effect in the spherical case does not appear to be as large. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. Longitudinal heat flux is always in the retrograde direction. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor-Proudman balance.

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

Kinematics of the Cometary Globules in the Gum Nebula

We report a study of the kinematics of the cometary globules in the Gum Nebula using the J = 1 yields 0 transition line of (CO-12)O. A morphological center for the system with which 60 percent of the globules are associated is identified. It is shown that the observed radial velocities of the heads of the globules are consistent with an expansion of the system. Systematic velocity gradients are present along some of the tails. The estimated expansion age and the tail stretching age are both about a few million years, suggesting a common origin for the expansion and the formation of the tails. The presence of young stars of similar ages in some of the globules points to star formation triggered by the same cause. Possible scenarios are briefly discussed.

Magnetic and spin evolution of pulsars

We explore the consequences of the model of spin-down-induced flux expulsion for the magnetic field evolution in solitary as well as in binary neutron stars. The spin evolution of pulsars, allowing for their field evolution according to this model, is shown to be consistent with the existing observational constraints in both low- and high-mass X-ray binary systems. The contribution from pulsars recycled in massive binaries to the observed excess in the number of low-field (10(11)-10(12) G) solitary pulsars is argued to be negligible in comparison with that of normal pulsars undergoing a 'restricted' field decay predicted by the adopted field decay model. Magnetic fields of neutron stars born in close binaries with intermediate- or high-mass main-sequence companions are predicted to decay down to values as low as similar to 10(6) G, which would leave them unobservable as pulsars during most of their lifetimes. The post-recycling evolution of some of these systems can, however, account for the observed binary pulsars having neutron star or massive white dwarf companions. Pulsars recycled in the disc population low-mass binaries are expected to have residual fields greater than or similar to 10(8) G, while for those processed in globular clusters larger residual fields are predicted because of the lower field strength of the neutron star at the epoch of binary formation. A value of tau similar to 1-2 x 10(7) yr for the mean value of the Ohmic decay time-scale in the crusts of neutron stars is suggested, based on the consistency of the model predictions with the observed distribution of periods and magnetic fields in the single and binary pulsars.

Lopsidedness in WHISP galaxies I. Rotation curves and kinematic lopsidedness

The frequently observed lopsidedness of the distribution of stars and gas in disc galaxies is still considered as a major problem in galaxy dynamics. It is even discussed as an imprint of the formation history of discs and the evolution of baryons in dark matter haloes. Here, we analyse a selected sample of 70 galaxies from the Westerbork Hi Survey of Spiral and Irregular Galaxies. The Hi data allow us to follow the morphology and the kinematics out to very large radii. In the present paper, we present the rotation curves and study the kinematic asymmetry. We extract the rotation curves of the receding and approaching sides separately and show that the kinematic behaviour of disc galaxies can be classified into five different types: symmetric velocity fields where the rotation curves of the receding and approaching sides are almost identical; global distortions where the rotation velocities of the receding and approaching sides have an offset that is constant with radius; local distortions leading to large deviations in the inner and negligible deviations in the outer parts (and vice versa); and distortions that divide the galaxies into two kinematic systems that are visible in terms of the different behaviour of the rotation curves of the receding and approaching sides, which leads to a crossing and a change in side. The kinematic lopsidedness is measured from the maximum rotation velocities, averaged over the plateau of the rotation curves. This gives a good estimate of the global lopsidedness in the outer parts of the sample galaxies. We find that the mean value of the perturbation parameter denoting the lopsided potential as obtained from the kinematic data is 0.056. Altogether, 36% of the sample galaxies are globally lopsided, which can be interpreted as the disc responding to a halo that was distorted by a tidal encounter. In Paper II, we study the morphological lopsidedness of the same sample of galaxies.

Lopsidedness in WHISP galaxies II. Morphological lopsidedness

The distribution of stars and gas in many galaxies is asymmetric. This so-called lopsidedness is expected to significantly affect the dynamics and evolution of the disc, including the star formation activity. Here, we measure the degree of lopsidedness for the gas distribution in a selected sample of 70 galaxies from the Westerbork Hi Survey of Spiral and Irregular Galaxies. This complements our earlier work (Paper I) where the kinematic lopsidedness was derived for the same galaxies. The morphological lopsidedness is measured by performing a harmonic decomposition of the surface density maps. The amplitude of lopsidedness A(1), the fractional value of the first Fourier component, is typically quite high (about 0.1) within the optical disc and has a constant phase. Thus, lopsidedness is a common feature in galaxies and indicates a global mode. We measure A(1) out to typically one to four optical radii, sometimes even further. This is, on average, four times larger than the distance to which lopsidedness was measured in the past using near-IR as a tracer of the old stellar component, and therefore provides a new, more stringent constraint on the mechanism for the origin of lopsidedness. Interestingly, the value of A(1) saturates beyond the optical radius. Furthermore, the plot of A(1) versus radius shows fluctuations that we argue are due to local spiral features. We also try to explain the physical origin of this observed disc lopsidedness. No clear trend is found when the degree of lopsidedness is compared to a measure of the isolation or interaction probability of the sample galaxies. However, this does not rule out a tidal origin if the lopsidedness is long-lived. In addition, we find that the early-type galaxies tend to be more morphologically lopsided than the late-type galaxies. Both results together indicate that lopsidedness has a tidal origin.

Starbursts triggered by cloud compression in interacting galaxies

We propose a physical mechanism for the triggering of starbursts in interacting spiral galaxies by shock compression of the pre-existing disk giant molecular clouds (GMCs). We show that as a disk GMC tumbles into the central region of a galaxy following a galactic tidal encounter, it undergoes a radiative shock compression by the pre-existing high pressure of the central molecular intercloud medium. The shocked outer shell of a GMC becomes gravitationally unstable, which results in a burst of star formation in the initially stable GMC. In the case of colliding galaxies with physical overlap such as Arp 244, the cloud compression is shown to occur due to the hot, high-pressure remnant gas resulting from the collisions of atomic hydrogen gas clouds from the two galaxies. The resulting values of infrared luminosity agree with observations. The main mode of triggered star formation is via clusters of stars, thus we can naturally explain the formation of young, luminous star clusters observed in starburst galaxies.

Incremental-angle and angular velocity estimation using a star sensor

A method is described for estimating the incremental angle and angular velocity of a spacecraft using integrated rate parameters with the help of a star sensor alone. The chief advantage of this method is that the measured stars need not be identified, whereas the identification of the stars is necessary in earlier methods. This proposed estimation can be carried out with all of the available measurements by a simple linear Kalman filter, albeit with a time-varying sensitivity matrix. The residuals of estimated angular velocity by the proposed spacecraft incremental-angle and angular velocity estimation method are as accurate as the earlier methods. This method also enables the spacecraft attitude to be reconstructed for mapping the stars into an imaginary unit sphere in the body reference frame, which will preserve the true angular separation of the stars. This will pave the way for identification of the stars using any angular separation or triangle matching techniques applied to even a narrow field of view sensor that is made to sweep the sky. A numerical simulation for inertial as well as Earth pointing spacecraft is carried out to establish the results.