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.

Hypersonic stagnation‐point boundary layers with massive blowing in the presence of a magnetic field

The effect of massive blowing rates on the steady laminar hypersonic boundary‐layer flow of an electrically conducting fluid in the stagnation region of an axisymmetric body with an applied magnetic field has been studied. The governing equations have been solved numerically by combining the implicit finite‐difference scheme with the quasi‐linearization technique. It is observed that the effect of massive blowing rates is to remove the viscous layer away from the boundary, whereas the effect of the magnetic field is just the opposite. It is also found that the velocity overshoot increases with blowing rates and also with magnetic field. The effect of the variation of the density‐viscosity product across the boundary layer is strong only when the blowing rate is small, but for the massive blowing rate the effect is negligible.

Infrared divergences and confinement in massive (2+1)-dimensional QED

In (2+1)-dimensional quantum electrodynamics with massless photons and massive matter fields, it is shown that the mass renormalization of the latter is infrared divergent at one loop. This result remains unchanged at two loops. A simple argument based on a similar divergence of the Coulomb potential leads us to conjecture that charged states are not observable in this model. This argument holds in 1+1 dimensions also.

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

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.

Swing amplification of nonaxisymmetric perturbations in stars and gas in a sheared galactic disk

We present a study of the growth of local, nonaxisymmetric perturbations in gravitationally coupled stars and gas in a differentially rotating galactic disk. The stars and gas are treated as two isothermal fluids of different velocity dispersions, with the stellar velocity dispersion being greater than that for the gas. We examine the physical effects of inclusion of a low-velocity dispersion component (gas) on the growth of non-axisymmetric perturbations in both stars and gas, as done for the axisymmetric case by Jog & Solomon. The amplified perturbations in stars and gas constitute trailing, material, spiral features which may be identified with the local spiral features seen in all spiral galaxies. The formulation of the two-fluid equations closely follows the one-fluid treatment by Goldreich & Lynden-Bell. The local, linearized perturbation equations in the sheared frame are solved to obtain the results for a temporary growth via swing amplification. The problem is formulated in terms of five dimensionless parameters-namely, the Q-factors for stars and gas, respectively; the gas mass fraction; the shearing rate in the galactic disk; and the length scale of perturbation. By using the observed values of these parameters, we obtain the amplifications and the pitch angles for features in stars and gas for dynamically distinct cases, as applicable for different regions of spiral galaxies. A real galaxy consisting of stars and gas may display growth of nonaxisymmetric perturbations even when it is stable against axisymmetric perturbations and/or when either fluid by itself is stable against non-axisymmetric perturbations. Due to its lower velocity dispersion, the gas exhibits a higher amplification than do the stars, and the amplified gas features are slightly more tightly wound than the stellar features. When the gas contribution is high, the stellar amplification and the range of pitch angles over which it can occur are both increased, due to the gravitational coupling between the two fluids. Thus, the two-fluid scheme can explain the origin of the broad spiral arms in the underlying old stellar populations of galaxies, as observed by Schweizer and Elmegreen & Elmegreen. The arms are predicted to be broader in gas-rich galaxies, as is indeed seen for example in M33. In the linear regime studied here, the arm contrast is shown to increase with radius in the inner Galaxy, in agreement with observations of external galaxies by Schweizer. These results follow directly due to the inclusion of gas in the problem.

New Massive Gravity and AdS(4) Counterterms

We show that the recently proposed Dirac-Born-Infeld extension of new massive gravity emerges naturally as a counterterm in four-dimensional anti-de Sitter space (AdS(4)). The resulting on-shell Euclidean action is independent of the cutoff at zero temperature. We also find that the same choice of counterterm gives the usual area law for the AdS(4) Schwarzschild black hole entropy in a cutoff-independent manner. The parameter values of the resulting counterterm action correspond to a c = 0 theory in the context of the duality between AdS(3) gravity and two-dimensional conformal field theory. We rewrite this theory in terms of the gauge field that is used to recast 3D gravity as a Chern-Simons theory.

Unsteady laminar compressible swirling flow with massive blowing

THE study of swirling boundary layers is of considerable importance in many rotodynamic machines such as rockets, jet engines, swirl generators, swirl atomizers, arc heaters, etc. For example, the introduction of swirl in a flow acceleration device such as a nozzle in a rocket engine promises efficient mass flow control. In nuclear rockets, swirl is used to retain the uranium atoms in the rocket chamber. With these applications in mind, Back1 and Muthanna and Nath2 have obtained the similarity solutions for a low-speed three-dimensional steady laminar compressible boundary layer with swirl inside an axisymmetric surface of variable cross section. The aim of the present analysis is to study the effect of massive blowing rates on the unsteady laminar swirling compressible boundary-layer flow of an axisymmetric body of arbitrary cross section when the freestream velocity and blowing rate vary with time. The type of swirl considered here is that of a free vortex superimposed on the longitudinal flow of a compressible fluid with variable properties. The analysis is applicable to external flow over a body as well as internal flow along a surface. For the case of external flow, strong blowing can have significant use in cooling the surface of hypervelocity vehicles, particularly when ablation occurs under large aerodynamic or radiative heating, but there may not be such an important application of strong blowing in the case of internal flow. The governing partial differential equations have been solved numerically using an implicit finite difference scheme with a quasilinearization technique.3 High temperature gas effects, such as radiation, dissociation, and ionization, etc., are not investigated. The nomenclature is usually that of Ref. 4 and is listed in the full paper.

Temperature profiles of accretion discs around rapidly rotating strange stars in general relativity: A comparison with neutron stars

We compute the temperature profiles of accretion discs around rapidly rotating strange stars, using constant gravitational mass equilibrium sequences of these objects, considering the full effect of general relativity. Beyond a certain critical value of stellar angular momentum (J), we observe the radius ( $r_{\rm orb}$) of the innermost stable circular orbit (ISCO) to increase with J (a property seen neither in rotating black holes nor in rotating neutron stars). The reason for this is traced to the crucial dependence of ${\rm d}r_{\rm orb}/{\rm d}J$ on the rate of change of the radial gradient of the Keplerian angular velocity at $r_{\rm orb}$ with respect to J. The structure parameters and temperature profiles obtained are compared with those of neutron stars, as an attempt to provide signatures for distinguishing between the two. We show that when the full gamut of strange star equation of state models, with varying degrees of stiffness are considered, there exists a substantial overlap in properties of both neutron stars and strange stars. However, applying accretion disc model constraints to rule out stiff strange star equation of state models, we notice that neutron stars and strange stars exclusively occupy certain parameter spaces. This result implies the possibility of distinguishing these objects from each other by sensitive observations through future X-ray detectors.

Establishing a Relation between the Mass and the Spin of Stellar-Mass Black Holes

Stellar mass black holes (SMBHs), forming by the core collapse of very massive, rapidly rotating stars, are expected to exhibit a high density accretion disk around them developed from the spinning mantle of the collapsing star. A wide class of such disks, due to their high density and temperature, are effective emitters of neutrinos and hence called neutrino cooled disks. Tracking the physics relating the observed (neutrino) luminosity to the mass, spin of black holes (BHs) and the accretion rate ((M) over dot) of such disks, here we establish a correlation between the spin and mass of SMBHs at their formation stage. Our work shows that spinning BHs are more massive than nonspinning BHs for a given (M) over dot. However, slowly spinning BHs can turn out to be more massive than spinning BHs if (M) over dot at their formation stage was higher compared to faster spinning BHs.

Nucleosynthesis in the accretion disks of type II collapsars

We investigate nucleosynthesis inside the gamma-ray burst (GRB) accretion disks formed by the Type II collapsars. In these collapsars, the core collapse of massive stars first leads to the formation of a proto-neutron star. After that, an outward moving shock triggers a successful supernova. However, the supernova ejecta lacks momentum and within a few seconds the newly formed neutron star gets transformed to a stellar mass black hole via massive fallback. The hydrodynamics of such an accretion disk formed from the fallback material of the supernova ejecta has been studied extensively in the past. We use these well-established hydrodynamic models for our accretion disk in order to understand nucleosynthesis, which is mainly advection dominated in the outer regions. Neutrino cooling becomes important in the inner disk where the temperature and density are higher. The higher the accretion rate (M) over dot is, the higher the density and temperature are in the disks. We deal with accretion disks with relatively low accretion rates: 0.001 M-circle dot s(-1) less than or similar to (M) over dot less than or similar to 0.01 M-circle dot s(-1) and hence these disks are predominantly advection dominated. We use He-rich and Si-rich abundances as the initial condition of nucleosynthesis at the outer disk, and being equipped with the disk hydrodynamics and the nuclear network code, we study the abundance evolution as matter inflows and falls into the central object. We investigate the variation in the nucleosynthesis products in the disk with the change in the initial abundance at the outer disk and also with the change in the mass accretion rate. We report the synthesis of several unusual nuclei like P-31, K-39, Sc-43, Cl-35 and various isotopes of titanium, vanadium, chromium, manganese and copper. We also confirm that isotopes of iron, cobalt, nickel, argon, calcium, sulphur and silicon get synthesized in the disk, as shown by previous authors. Much of these heavy elements thus synthesized are ejected from the disk via outflows and hence they should leave their signature in observed data.

Nucleosynthesis in the outflows associated with accretion disks of type II collapsars

We investigate nucleosynthesis inside the outflows from gamma-ray burst (GRB) accretion disks formed by the Type II collapsars. In these collapsars, massive stars undergo core collapse to form a proto-neutron star initially, and a mild supernova (SN) explosion is driven. The SN ejecta lack momentum, and subsequently this newly formed neutron star gets transformed to a stellar mass black hole via massive fallback. The hydrodynamics and the nucleosynthesis in these accretion disks have been studied extensively in the past. Several heavy elements are synthesized in the disk, and much of these heavy elements are ejected from the disk via winds and outflows. We study nucleosynthesis in the outflows launched from these disks by using an adiabatic, spherically expanding outflow model, to understand which of these elements thus synthesized in the disk survive in the outflow. While studying this, we find that many new elements like isotopes of titanium, copper, zinc, etc., are present in the outflows. Ni-56 is abundantly synthesized in most of the cases in the outflow, which implies that the outflows from these disks in a majority of cases will lead to an observable SN explosion. It is mainly present when outflow is considered from the He-rich, Ni-56/Fe-54-rich zones of the disks. However, outflow from the Si-rich zone of the disk remains rich in silicon. Although emission lines of many of these heavy elements have been observed in the X-ray afterglows of several GRBs by Chandra, BeppoSAX, XMM-Newton, etc., Swift seems to have not yet detected these lines.

A DYNAMO MODEL OF MAGNETIC ACTIVITY IN SOLAR-LIKE STARS WITH DIFFERENT ROTATIONAL VELOCITIES

We attempt to provide a quantitative theoretical explanation for the observations that Ca II H/K emission and X-ray emission from solar-like stars increase with decreasing Rossby number (i.e., with faster rotation). Assuming that these emissions are caused by magnetic cycles similar to the sunspot cycle, we construct flux transport dynamo models of 1M(circle dot) stars rotating with different rotation periods. We first compute the differential rotation and the meridional circulation inside these stars from a mean-field hydrodynamics model. Then these are substituted in our dynamo code to produce periodic solutions. We find that the dimensionless amplitude f(m) of the toroidal flux through the star increases with decreasing rotation period. The observational data can be matched if we assume the emissions to go as the power 3-4 of f(m). Assuming that the Babcock-Leighton mechanism saturates with increasing rotation, we can provide an explanation for the observed saturation of emission at low Rossby numbers. The main failure of our model is that it predicts an increase of the magnetic cycle period with increasing rotation rate, which is the opposite of what is found observationally. Much of our calculations are based on the assumption that the magnetic buoyancy makes the magnetic flux tubes rise radially from the bottom of the convection zone. Taking into account the fact that the Coriolis force diverts the magnetic flux tubes to rise parallel to the rotation axis in rapidly rotating stars, the results do not change qualitatively.

Large-Scale Multiuser SM-MIMO Versus Massive MIMO

Spatial modulation (SM) is attractive for multiantenna wireless communications. SM uses multiple transmit antenna elements but only one transmit radio frequency (RF) chain. In SM, in addition to the information bits conveyed through conventional modulation symbols (e.g., QAM), the index of the active transmit antenna also conveys information bits. In this paper, we establish that SM has significant signal-to-noise (SNR) advantage over conventional modulation in large-scale multiuser (multiple-input multiple-output) MIMO systems. Our new contribution in this paper addresses the key issue of large-dimension signal processing at the base station (BS) receiver (e.g., signal detection) in large-scale multiuser SM-MIMO systems, where each user is equipped with multiple transmit antennas (e.g., 2 or 4 antennas) but only one transmit RF chain, and the BS is equipped with tens to hundreds of (e.g., 128) receive antennas. Specifically, we propose two novel algorithms for detection of large-scale SM-MIMO signals at the BS; one is based on message passing and the other is based on local search. The proposed algorithms achieve very good performance and scale well. For the same spectral efficiency, multiuser SM-MIMO outperforms conventional multiuser MIMO (recently being referred to as massive MIMO) by several dBs. The SNR advantage of SM-MIMO over massive MIMO can be attributed to: (i) because of the spatial index bits, SM-MIMO can use a lower-order QAM alphabet compared to that in massive MIMO to achieve the same spectral efficiency, and (ii) for the same spectral efficiency and QAM size, massive MIMO will need more spatial streams per user which leads to increased spatial interference.

In search of slow-moving ionizing massive particles

Many particles proposed by theories, such as GUT monopoles, nuclearites and 1/5 charge superstring particles, can be categorized as Slow-moving, Ionizing, Massive Particles (SIMPs).

Detailed calculations of the signal-to-noise ratios in vanous acoustic and mechanical methods for detecting such SIMPs are presented. It is shown that the previous belief that such methods are intrinsically prohibited by the thermal noise is incorrect, and that ways to solve the thermal noise problem are already within the reach of today's technology. In fact, many running and finished gravitational wave detection ( GWD) experiments are already sensitive to certain SIMPs. As an example, a published GWD result is used to obtain a flux limit for nuclearites.

The result of a search using a scintillator array on Earth's surface is reported. A flux limit of 4.7 x 10^(-12) cm^(-2)sr^(-1)s^(-1) (90% c.l.) is set for any SIMP with 2.7 x 10^(-4) less than β less than 5 x 10^(-3) and ionization greater than 1/3 of minimum ionizing muons. Although this limit is above the limits from underground experiments for typical supermassive particles (10^(16)GeV), it is a new limit in certain β and ionization regions for less massive ones (~10^9 GeV) not able to penetrate deep underground, and implies a stringent limit on the fraction of the dark matter that can be composed of massive electrically and/ or magnetically charged particles.

The prospect of the future SIMP search in the MACRO detector is discussed. The special problem of SIMP trigger is examined and a circuit proposed, which may solve most of the problems of the previous ones proposed or used by others and may even enable MACRO to detect certain SIMP species with β as low as the orbital velocity around the earth.