Using an idealized set-up, we investigate the dynamical role of cosmic rays (CRs) in the early stages of galactic outflows for galaxies of halo masses $10^{8}$, $10^{11}$ and $10^{12}$ $M_\odot$. The outflow is launched from a central region in the galactic disk where we consider three different constant star formation rates ($0.1$, $1$, and $10$ $M_{\odot}$ yr$^{-1}$) over a dynamical timescale of $50$ Myr. We determine the temperature distribution of the gas and find that CRs can reduce the temperature of the shocked gas, which is consistent with previous results. However, we show that CRs do not have any noticeable effect on the mass loading by the outflow. We find that CRs can reduce the size of the outflow, which contradicts previous claims of efficient dynamical impact of CRs; however, it is consistent with earlier theoretical models of cosmic ray driven blastwave as well as stellar wind. We discuss the dependence of our results on CR injection prescriptions and compare them with earlier studies. We conclude that in the early stages of galactic outflows the dynamical role of CRs is not important.

Different measurements of the Hubble constant ($H_{0}$) are not consistent and a tension between the CMB based methods and cosmic distance ladder based methods has been observed. Measurements from various distance based methods are also inconsistent. To aggravate the problem, same cosmological probe (Type Ia SNe for instance) calibrated through different methods also provide different value of $H_{0}$. We compare various distance ladder based methods through the already available unique data obtained from Hubble Space Telescope (HST). Our analysis is based on parametric (T-test) as well as non-parametric statistical methods such as the Mann-Whitney U test and Kolmogorov-Smirnov test. Our results show that different methods provide different values of $H_0$ and the differences are statistically significant. The biases in the calibration would not account for these differences as the data has been taken from a single telescope with common calibration scheme. The unknown physical effects or issues with the empirical relations of distance measurement from different probes could give rise to these differences.

Low energy ground-based cosmic ray air shower experiments generally have energy threshold in the range of a few tens to a few hundreds of TeV. The shower observables are measured indirectly with an array of detectors. The atmospheric absorption of low energy secondaries limits their detection frequencies at the Earth's surface. However, due to selection effects, a tiny fraction of low energy showers, which are produced in the lower atmosphere can reach the observational level. But, due to less information of shower observables, the reconstruction of these showers are arduous. Hence, it is believed that direct measurements by experiments aboard on satellites and balloon flights are more reliable at low energies. Despite having very small efficiency ($\sim$0.1%) at low energies, the large acceptance ($\sim$5 m$^2$sr) of GRAPES-3 experiment allows observing primary cosmic rays down below to $\sim$1 TeV and opens up the possibility to measure primary energy spectrum spanning from a few TeV to beyond cosmic ray knee (up to 10$^{16}$ eV), covering five orders of magnitude. The GRAPES-3 energy threshold for primary protons through Monte Carlo simulations are calculated, which gives reasonably good agreement with data. Furthermore, the total efficiencies and acceptance are also calculated for protons primaries. The ability of GRAPES-3 experiment to cover such a broader energy range may provide a unique handle to bridge the energy spectrum between direct measurements at low energies and indirect measurements at ultra-high energies.

Motivated by previous studies in literature about the potential importance of relativistic corrections to galaxy cluster hydrostatic masses, we calculate the masses of 12 relaxed clusters (with Chandra X-ray data) using the Tolman-Oppenheimer-Volkov (TOV) equation of hydrostatic equilibrium and the ideal gas equation of state. Analytical formulae for gas density and temperature profiles for these clusters, previously derived by Vikhlinin et al (astro-ph/0507092) were used to obtain these masses. We compare the TOV-based masses with those obtained using the corresponding Newtonian equation of hydrostatic equilibrium. We find that the fractional relative difference between the two masses are negligible, corresponding to $\sim \mathcal{O}(10^{-5})$.

Although galaxy clusters have proved to be wonderful laboratories for testing a plethora of modified gravity theories and other exotic alternatives to $\Lambda$CDM, until a few years ago, there was only one paper (from 1974), which obtained a limit on graviton mass of $\mathcal{O}(10^{-29})$ eV with clusters. To rectify this, in the last few years multiple works have obtained different bounds on graviton mass using single cluster data as well as stacking galaxy catalogs. We review these recent limits on graviton mass using galaxy clusters obtained using disparate methods.

The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $\gamma$-ray sources in the multi-TeV energy range. The GRAPES-3 an EAS array located at Ooty in India (11.4$^{\circ}$N, 76.7$^{\circ}$E, 2200 m altitude) is designed to study $\gamma$-rays in the TeV-PeV energy range. It comprises of a dense array of 400 plastic scintillators deployed over an area of 25000 m$^2$ and a large area (560 m$^2$) muon telescope. A new statistical method allowed real time determination of the propagation delay of each detector in the GRAPES-3 array. The shape of shower front is known to be curved and here the details of a new method developed for accurate measurement of the shower front curvature is presented. These two developments have led to a sizable improvement in the angular resolution of GRAPES-3 array. It is shown that the curvature depends on the size and age of an EAS. By employing two different techniques, namely, the odd-even and the left-right methods, independent estimates of the angular resolution are obtained. The odd-even method estimates the best achievable resolution of the array. For obtaining the angular resolution, the left-right method is used after implementing the size and age dependent curvature corrections. A comparison of the angular resolution as a function of EAS energy by these two methods shows them be virtually indistinguishable. The angular resolution of GRAPES-3 array is 47$^{\prime}$ for energies E$>$5 TeV and improves to 17$^{\prime}$ at E$>$100 TeV and finally approaching 10$^{\prime}$ at E$>$500 TeV.

Cosmic ray (CR) sources leave signatures in the isotopic abundances of CRs. Current models of Galactic CRs that consider supernovae (SNe) shocks as the main sites of particle acceleration cannot satisfactorily explain the higher ${\rm ^{22}Ne/^{20}Ne}$ ratio in CRs compared to the interstellar medium. Although stellar winds from massive stars have been invoked, their contribution relative to SNe ejecta has been taken as a free parameter. Here we present a theoretical calculation of the relative contributions of wind termination shocks (WTSs) and SNe shocks in superbubbles, based on the hydrodynamics of winds in clusters, the standard stellar mass function, and stellar evolution theory. We find that the contribution of WTSs towards the total CR production is at least $25\%$, which rises to $\gtrsim 50\%$ for young ($\lesssim 10$ Myr) clusters, and explains the observed $^{22}{\rm Ne}/^{20} {\rm Ne}$ ratio. We argue that since the progenitors of apparently isolated supernovae remnants (SNRs) are born in massive star clusters, both WTS and SNe shocks can be integrated into a combined scenario of CRs being accelerated in massive clusters. This scenario is consistent with the observed ratio of SNRs to $\gamma$-ray bright ($L_\gamma \gtrsim 10^{35}$ erg s$^{-1}$) star clusters, as predicted by star cluster mass function. Moreover, WTSs can accelerate CRs to PeV energies, and solve other longstanding problems of the standard supernova paradigm of CR acceleration.

We report the results obtained from a detailed timing and spectral studies of Be/X-ray binary pulsar 2S 1417-624 using data from Swift and NuSTAR observatories. The observations were carried out at the peak of a giant outburst of the pulsar in 2018. X-ray pulsations at 17.475 s were detected in the source light curves up to ~79 keV. The evolution of the pulse profiles with energy was found to be complex. A four-peaked profile at lower energies gradually evolved into a double-peak structure at higher energies. The pulsed fraction of the pulsar, calculated from the NuSTAR observation was found to follow an anti-correlation trend with luminosity as observed during previous giant X-ray outburst studies in 2009. The broadband spectrum of the pulsar is well described by a composite model consisting of a cutoff power law model modified with the interstellar absorption, a thermal blackbody component with a temperature of $\approx$1 keV and a Gaussian function for the 6.4 keV iron emission line. Though the pulsar was observed at the peak of the giant outburst, there was no signature of presence of any cyclotron line feature in the spectrum. The radius of the blackbody emitting region was estimated to be $\approx$2 km, suggesting that the most probable site of its origin is the stellar surface of the neutron star. Physical models were also explored to understand the emission geometry of the pulsar and are discussed in the paper.

The GRAPES-3 muon telescope located in Ooty, India records 4x10^9 muons daily. These muons are produced by interaction of primary cosmic rays (PCRs) in the atmosphere. The high statistics of muons enables GRAPES-3 to make precise measurement of various sun-induced phenomenon including coronal mass ejections (CME), Forbush decreases, geomagnetic storms (GMS) and atmosphere acceleration during the overhead passage of thunderclouds. However, the understanding and interpretation of observed data requires Monte Carlo (MC) simulation of PCRs and subsequent development of showers in the atmosphere. CORSIKA is a standard MC simulation code widely used for this purpose. However, these simulations are time consuming as large number of interactions and decays need to be taken into account at various stages of shower development from top of the atmosphere down to ground level. Therefore, computing resources become an important consideration particularly when billion of PCRs need to be simulated to match the high statistical accuracy of the data. During the GRAPES-3 simulations, it was observed that over 60% of simulated events don't really reach the Earth's atmosphere. The geomagnetic field (GMF) creates a threshold to PCRs called cutoff rigidity Rc, a direction dependent parameter below which PCRs can't reach the Earth's atmosphere. However, in CORSIKA there is no provision to set a direction dependent threshold. We have devised an efficient method that has taken into account of this Rc dependence. A reduction by a factor ~3 in simulation time and ~2 in output data size was achieved for GRAPES-3 simulations. This has been incorporated in CORSIKA version v75600 onwards. Detailed implementation of this along the potential benefits are discussed in this work.

Cosmic rays (CRs) are frequently modeled as an additional fluid in hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations of astrophysical flows. The standard CR two-fluid model is described in terms of three conservation laws (expressing conservation of mass, momentum and total energy) and one additional equation (for the CR pressure) that cannot be cast in a satisfactory conservative form. The presence of non-conservative terms with spatial derivatives in the model equations prevents a unique weak solution behind a shock. We investigate a number of methods for the numerical solution of the two-fluid equations and find that, in the presence of shock waves, the results depend on the choice of the numerical methods (spatial reconstruction, time stepping, and the CFL number) and the adopted discretization. Nevertheless, all methods converge to a unique result only if the energy partition between the thermal and non-thermal fluids at the shock is prescribed a priori. This highlights the closure problem of the two-fluid equations at shocks. We suggest a robust method where the solutions are insensitive to the numerical method. Comparison with the currently used methods, critical test problems, and future directions are discussed.

We present new limits on the graviton mass using a sample of 12 relaxed galaxy clusters, for which temperature and gas density profiles were derived by Vikhlinin et al (astro-ph/0507092) using Chandra X-ray observations. These limits can be converted to a bound on the graviton mass, assuming a non-zero graviton mass would lead to a Yukawa potential at these scales. For this purpose, we first calculate the total dynamical mass from the hydrostatic equilibrium equation in Yukawa gravity and then compare it with the corresponding mass in Newtonian gravity. We calculate a 90 % c.l. lower/upper limit on the graviton Compton wavelength/ mass for each of the 12 clusters in the sample. The best limit is obtained for Abell 2390, corresponding to $\lambda_g > 3.58\times 10^{19}$ km or $m_g<3.46 \times 10^{-29}$ eV. This is the first proof of principles demonstration of setting a limit on the graviton mass using a sample of related galaxy clusters with X-ray measurements and can be easily applied to upcoming X-ray surveys such as eRosita.

In the last few years, there has been a resurgence of interest in obtaining observational bounds on the graviton mass, following the detection of gravitational waves, because of the versatility of massive graviton theories in resolving multiple problems in cosmology and fundamental physics. In this work, we apply the method proposed in Rana et al.(arXiv:1801.03309), which consists of looking for Yukawa-like fall off in the gravitational potential, to stacked galaxy cluster catalogs from three disparate surveys. These include catalogs from 2500 sq. degree SPT-SZ survey, the Planck all-sky SZ catalog, and a redMaPPer selected catalog from 10,000 sq. degree of SDSS-DR8 data. The 90\% c.l. limits which we obtained on the graviton mass using SPT, Planck and SDSS are: $m_g < 4.73 \times 10^{-30}$ eV, $3.0 \times 10^{-30}$ eV, and $1.27 \times 10^{-30}$ eV respectively; or in terms of Compton wavelength are $\lambda_g >2.62 \times 10^{20}$ km, $4.12 \times 10^{20}$ km, $9.76 \times 10^{20}$ km. These limits are about five times more stringent than the previous best bound from galaxy clusters.

We present the results obtained from timing and spectral studies of the accretion powered X-ray pulsar 2S~1417-624 during a giant outburst in 2009 by using {\it Rossi X-ray Timing Explorer (RXTE)} observations. X-ray pulsations were detected in the light curves obtained from all epochs of observations. The pulsar was found to be spinning-up during the outburst. The pulse profiles were observed to be strongly dependent on photon energy and luminosity. A double peaked profile at lower luminosity evolved into a triple peaked profile at the peak of the outburst which is further reverted back to a double peaked structure during the decay of the outburst. An anti-correlation was also observed between the pulse fraction and the source flux. The 3-70 keV energy spectrum of pulsar was well described with a power law modified with high energy cutoff model along with an iron fluorescence line at 6.4 keV. Based on the evolution of pulse profile, pulse fraction and spectral parameters across observed luminosity, we interpret our results in terms of changes in the pulsar beam configuration from sub-critical to super-critical regimes.

This work looks at some definitive signatures of magnetars, in particular of period closures accompanied by a decline of X-ray radiation in two models. We review some of the previous works which are based on the well known dynamo model in which the star is born with a period of a few milliseconds at high temperatures. In such a convection regime the dynamo mechanism can amplify the the magnetic fields to the magnetar value. This is in contrast to a screened core model which posits that a high density phase transition occurs in the inner core of magnetars that dynamically aligns all the neutron magnetic moments producing a large magnetic field in the core. The accompanying change of flux gives rise to shielding or screening currents in the surrounding high conductivity plasma that do not permit the field to exit to the surface. Ambipolar diffusion then transports the field to the crust dissipating energy in neutrinos and X-rays. The up-welling field cleaves the crust resulting in flares and X-ray radiation from ohmic dissipation in the crust till the screening currents are spent and the surface polar field attains its final value. In the dynamo model the polar magnetic field decreases with time whereas in our screened model it increases to its final value. One consequence of this is that in the latter model, as a function of time and period, the ratio of the dipole radiation loss, $\dot E$ to the X-ray luminosity, $ L_X$, is a 'U' curve, indicating that it is the exponential decline in $L_X$, that brings closure to the periods that are observed for magnetars.

We present detailed spectral and timing studies using a NuSTAR observation of GX 1+4 in October 2015 during an intermediate intensity state. The measured spin period of 176.778 s is found to be one of the highest values since its discovery. In contrast to a broad sinusoidal-like pulse profile, a peculiar sharp peak is observed in profiles below ~25 keV. The profiles at higher energies are found to be significantly phase-shifted compared to the soft X-ray profiles. Broadband energy spectra of GX 1+4, obtained from NuSTAR and Swift observations, are described with various continuum models. Among these, a two component model consisting of a bremsstrahlung and a blackbody component is found to best-fit the phase-averaged and phase-resolved spectra. Physical models are also used to investigate the emission mechanism in the pulsar, which allows us to estimate the magnetic field strength to be in $\sim$(5-10)$\times$10$^{12}$ G range. Phase-resolved spectroscopy of NuSTAR observation shows a strong blackbody emission component in a narrow pulse phase range. This component is interpreted as the origin of the peculiar peak in the pulse profiles below $\le$25 keV. The size of emitting region is calculated to be $\sim$400 m. The bremsstrahlung component is found to dominate in hard X-rays and explains the nature of simple profiles at high energies.

X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models.

We use $1$D and $3$D two-fluid cosmic ray (CR) hydrodynamic simulations to investigate the role of CRs in the vicinity of a compact young star cluster. We model a self-gravitating cloud (density profile $\rho \propto r^{-1}$), include important thermal and non-thermal processes, and explore two different CR injection scenarios. We show that if internal shocks in the wind-driving region are the main site for CR acceleration, then the resulting $\gamma$-ray luminosity ($L_{\rm \gamma}$) can reach $\approx 5\%$ of the mechanical luminosity ($L_{\rm w}$), independent of the fraction of wind energy ($\sim 1-20\%$) injected into CRs. In contrast, if the forward/reverse shock of a bubble is the injection site then $L_{\rm \gamma}$ increases linearly with the CR injection fraction, as expected analytically. We find that the X-ray luminosity ($L_{\rm x}$) in the forward/reverse shock injection scenario is $\gtrsim 10^{-3} L_{\rm w}$, which is $\sim 10$ times larger than in the central wind-driving injection case. We predict the corresponding range of the synchrotron radio luminosity. We show how multi-wavelength observations can constrain the CR parameters. Comparing the predicted multi-wavelength luminosities with those of 30 Doradus we identify the reverse shock as the most probable CR injection site, and that thermal conduction is important. We do not find significant dynamical impact of CRs in our models.

The GRAPES-3 muon telescope in Ooty, India had claimed detection of a 2 hour (h) high-energy ($\sim$20 GeV) burst of galactic cosmic-rays (GCRs) through a $>$50$\sigma$ surge in GeV muons, was caused by reconnection of the interplanetary magnetic field (IMF) in the magnetosphere that led to transient weakening of Earth's magnetic shield. This burst had occurred during a G4-class geomagnetic storm (storm) with a delay of $\frac{1}{2}$h relative to the coronal mass ejection (CME) of 22 June 2015 (Mohanty et al., 2016). However, recently a group interpreted the occurrence of the same burst in a subset of 31 neutron monitors (NMs) to have been the result of an anisotropy in interplanetary space (Evenson et al., 2017) in contrast to the claim in (Mohanty et al., 2016). A new analysis of the GRAPES-3 data with a fine 10.6$^{\circ}$ angular segmentation shows the speculation of interplanetary anisotropy to be incorrect, and offers a possible explanation of the NM observations. The observed 28 minutes (min) delay of the burst relative to the CME can be explained by the movement of the reconnection front from the bow shock to the surface of Earth at an average speed of 35 km/s, much lower than the CME speed of 700 km/s. This measurement may provide a more accurate estimate of the start of the storm.

Observations strongly suggest that filaments in galactic molecular clouds are in a non-thermal state. As a simple model of a filament we study a two-dimensional system of self-gravitating point particles by means of numerical simulations of the dynamics, with various methods: direct $N$-body integration of the equations of motion, particle-in-cell simulations and a recently developed numerical scheme that includes multiparticle collisions in a particle-in-cell approach. Studying the collapse of Gaussian overdensities we find that after the damping of virial oscillations the system settles in a non-thermal steady state whose radial density profile is similar to the observed ones, thus suggesting a dynamical origin of the non-thermal states observed in real filaments. Moreover, for sufficiently cold collapses the density profiles are anticorrelated with the kinetic temperature, i.e., exhibit temperature inversion, again a feature that has been found in some observations of filaments. The same happens in the state reached after a strong perturbation of an initially isothermal cylinder. Finally, we discuss our results in the light of recent findings in other contexts (including non-astrophysical ones) and argue that the same kind of non-thermal states may be observed in any physical system with long-range interactions.

We study the origin of the non-thermal profiles observed in filamentary structures in galactic molecular clouds by means of numerical dynamical simulations. We find that such profiles are intrinsic features of the end products of dissipationless collapse in cylindrical symmetry. Moreover, for sufficiently cold initial conditions, we obtain end states characterized by markedly anticorrelated radial density and temperature profiles. Gravitational, dissipationless dynamics alone is thus sufficient to reproduce, at least qualitatively, many of the properties of the observed non-thermal structures.

In this paper we present the design and measured performance of a novel cryogenic motor based on a superconducting magnetic bearing (SMB). The motor is tailored for use in millimeter-wave half-wave plate (HWP) polarimeters, where a HWP is rapidly rotated in front of a polarization analyzer or polarization-sensitive detector. This polarimetry technique is commonly used in cosmic microwave background (CMB) polarization studies. The SMB we use is composed of fourteen yttrium barium copper oxide (YBCO) disks and a contiguous neodymium iron boron (NdFeB) ring magnet. The motor is a hollow-shaft motor because the HWP is ultimately installed in the rotor. The motor presented here has a 100 mm diameter rotor aperture. However, the design can be scaled up to rotor aperture diameters of approximately 500 mm. Our motor system is composed of four primary subsystems: (i) the rotor assembly, which includes the NdFeB ring magnet, (ii) the stator assembly, which includes the YBCO disks, (iii) an incremental encoder, and (iv) the drive electronics. While the YBCO is cooling through its superconducting transition, the rotor is held above the stator by a novel hold and release mechanism (HRM). The encoder subsystem consists of a custom-built encoder disk read out by two fiber optic readout sensors. For the demonstration described in this paper, we ran the motor at 50 K and tested rotation frequencies up to approximately 10 Hz. The feedback system was able to stabilize the the rotation speed to approximately 0.4%, and the measured rotor orientation angle uncertainty is less than 0.15 deg. Lower temperature operation will require additional development activities, which we will discuss.

We present a comprehensive timing and spectral studies of Be/X-ray binary pulsar EXO 2030+375 using extensive Rossi X-ray Timing Explorer observations from 1995 till 2011, covering numerous Type I and 2006 Type II outbursts. Pulse profiles of the pulsar were found to be strongly luminosity dependent. At low luminosity, the pulse profile consisted of a main peak and a minor peak that evolved into a broad structure at high luminosity with a significant phase shift. A narrow and sharp absorption dip, also dependent on energy and luminosity, was detected in the pulse profile. Comparison of pulse profiles showed that the features at a particular luminosity are independent of type of X-ray outbursts. This indicates that the emission geometry is solely a function of mass accretion rate. The broadband energy spectrum was described with a partial covering high energy cutoff model as well as a physical model based on thermal and bulk Comptonization in accretion column. We did not find any signature of cyclotron resonance scattering feature in the spectra obtained from all the observations. A detailed analysis of spectral parameters showed that, depending on source luminosity, the power-law photon index was distributed in three distinct regions. It suggests the phases of spectral transition from sub-critical to super-critical regimes in the pulsar as proposed theoretically. A region with constant photon index was also observed in ~(2-4) x 10^37 erg/s range, indicating critical luminosity regime in EXO 2030+375.

Using analytic methods and $1$-D two-fluid simulations, we study the effect of cosmic rays (CRs) on the dynamics of interstellar superbubbles (ISBs) driven by multiple supernovae (SNe)/stellar winds in OB associations. In addition to CR advection and diffusion, our models include thermal conduction and radiative cooling. We find that CR injection at the reverse shock or within a central wind-driving region can affect the thermal profiles of ISBs and hence their X-ray properties. Even if a small fraction ($10-20\%$) of the total mechanical power is injected into CRs, a significant fraction of the ram pressure at the reverse shock can be transferred to CRs. The energy transfer becomes efficient if (1) the reverse shock gas Mach number exceeds a critical value ($M_{\rm th}\gtrsim 12$) and (2) the CR acceleration time scale $\tau_{\rm acc}\sim \kappa_{\rm cr}/v^2$ is shorter than the dynamical time, where $\kappa_{\rm cr}$ is CR diffusion constant and $v$ is the upstream velocity. We show that CR affected bubbles can exhibit a volume averaged hot gas temperature $1-5\times10^{6}$ K, lower by a factor of $2-10$ than without CRs. Thus CRs can potentially solve the long-standing problem of the observed low ISB temperatures.

Gamow-Teller (GT) strength distributions (B(GT)) in electron-capture (EC) daughters stemming from the parent ground state are computed with the shell-model in the full pf-shell space, with quasi-particle random-phase approximation (QRPA) in the formalism of Krumlinde and M\"oller and with an Approximate Method (AM) for assigning an effective B(GT). These are compared to data available from decay and charge-exchange (CE) experiments across titanium isotopes in the pf-shell from A=43 to A=62, the largest set available for any chain of isotopes in the pf-shell. The present study is the first to examine B(GT) and the associated EC rates across a particular chain of isotopes with the purpose of examining rate sensitivities as neutron number increases. EC rates are also computed for a wide variety of stellar electron densities and temperatures providing concise estimates of the relative size of rate sensitivities for particular astrophysical scenarios. This work underscores the astrophysical motivation for CE experiments in inverse kinematics for nuclei away from stability at the luminosities of future Radioactive Ion Beam Facilities.

Assuming the Central Limit Theorem, experimental uncertainties in any data set are expected to follow the Gaussian distribution with zero mean. We propose an elegant method based on Kolmogorov-Smirnov statistic to test the above; and apply it on the measurement of Hubble constant which determines the expansion rate of the Universe. The measurements were made using Hubble Space Telescope. Our analysis shows that the uncertainties in the above measurement are non-Gaussian.

Energetic winds and radiation from massive star clusters push the surrounding gas and blow superbubbles in the interstellar medium (ISM). Using 1-D hydrodynamic simulations, we study the role of radiation in the dynamics of superbubbles driven by a young star cluster of mass $10^{6}$ M$_{\odot}$. We have considered a realistic time evolution of the mechanical power as well as radiation power of the star cluster, and detailed heating and cooling processes. We find that the ratio of the radiation pressure on the shell (shocked ISM) to the thermal pressure ($\sim10^{7}$ K) of the shocked wind region is almost independent of the ambient density, and it is greater than unity before $\lesssim 1$ Myr. We explore the parameter space of density and dust opacity of the ambient medium, and find that the size of the hot gas ($\sim$ 10$^{7}$ K) cavity is insensitive to the dust opacity ($\sigma_{d}\approx(0.1-1.5)\times 10^{-21}$ cm$^{2}$), but the structure of the photoionized ($\sim10^4$ K) gas depends on it. Most of the radiative losses occur at $\sim10^{4}$ K, with sub-dominant losses at $\lesssim 10^3$ K and $\sim10^{6}-10^{8}$ K. The superbubbles can retain as high as $\sim 10\%$ of its input energy, for an ambient density of $10^{3}\,m{\rm_{H}\,cm^{-3}}$. We discuss the role of ionization parameter and recombination-averaged density in understanding the dominant feedback mechanism. Finally, we compare our results with the observations of 30 Doradus.

The nature of random errors in any data set that is Gaussian is a well established fact according to the Central Limit Theorem. Supernovae type Ia data have played a crucial role in major discoveries in cosmology. Unlike in laboratory experiments, astronomical measurements can not be performed in controlled situations. Thus, errors in astronomical data can be more severe in terms of systematics and non-Gaussianity compared to those of laboratory experiments. In this paper, we use the Kolmogorov-Smirnov statistic to test non-Gaussianity in high-z supernovae data. We apply this statistic to four data sets, i.e., Gold data(2004), Gold data(2007), Union2 catalogue and the Union2.1 data set for our analysis. Our results shows that in all four data sets the errors are consistent with the Gaussian distribution.

We use the $\Delta_{\chi^2}$ statistic introduced in \cite{gup08,gup10} to study directional dependence, in the high-z supernovae data. This dependence could arise due to departures from the cosmological principle or from direction dependent statistical systematics in the data. We apply our statistic to the gold data set from \cite{rie04} and \cite{rie07}, and Union2 catalogue from \cite{aman10}. Our results show that all the three data sets show a weak but consistent direction dependence. In 2007 data errors are Gaussian, however other two data sets show non-Gaussian features.

Temperature inversions occur in nature, e.g., in the solar corona and in interstellar molecular clouds: somewhat counterintuitively, denser parts of the system are colder than dilute ones. We propose a simple and appealing way to spontaneously generate temperature inversions in systems with long-range interactions, by preparing them in inhomogeneous thermal equilibrium states and then applying an impulsive perturbation. In similar situations, short-range systems would typically relax to another thermal equilibrium, with uniform temperature profile. By contrast, in long-range systems, the interplay between wave-particle interaction and spatial inhomogeneity drives the system to nonequilibrium stationary states that generically exhibit temperature inversion. We demonstrate this mechanism in a simple mean-field model and in a two-dimensional self-gravitating system. Our work underlines the crucial role the range of interparticle interaction plays in determining the nature of steady states out of thermal equilibrium.

Aims. Forbush decrease (FD) is a transient decrease followed by a gradual recovery in the observed galactic cosmic ray intensity. We seek to understand the relationship between the FDs and near-Earth interplanetary magnetic field (IMF) enhancements associated with solar coronal mass ejections (CMEs). Methods. We use muon data at cutoff rigidities ranging from 14 to 24 GV from the GRAPES-3 tracking muon telescope to identify FD events. We select those FD events that have a reasonably clean profile, and magnitude > 0.25%. We use IMF data from ACE/WIND spacecrafts. We look for correlations between the FD profile and that of the one hour averaged IMF. We ask if the diffusion of high energy protons into the large scale magnetic field is the cause of the lag observed between the FD and the IMF. Results. The enhancement of the IMF associated with FDs occurs mainly in the shock-sheath region, and the turbulence level in the magnetic field is also enhanced in this region. The observed FD profiles look remarkably similar to the IMF enhancement profiles. The FDs typically lag the IMF enhancement by a few hours. The lag corresponds to the time taken by high energy protons to diffuse into the magnetic field enhancement via cross-field diffusion. Conclusions. Our findings show that high rigidity FDs associated with CMEs are caused primarily by the cumulative diffusion of protons across the magnetic field enhancement in the turbulent sheath region between the shock and the CME.

Cosmic ray Forbush decreases (FDs) are usually thought to be due to Earth-directed coronal mass ejections (CMEs) from the Sun and their associated shocks. When CMEs and their shocks reach the Earth, they cause magnetic field compressions. We seek to understand the relation between these magnetic field compressions and FDs at rigidities between 12 and 42 GV using data from the GRAPES-3 instrument at Ooty. We find that the shapes of the Forbush decrease profiles show a startling similarity to that of the magnetic field compression in the near-Earth IP medium. We seek to understand the implications of this interesting result.

Temperature inversion due to velocity filtration, a mechanism originally proposed to explain the heating of the solar corona, is demonstrated to occur also in a simple paradigmatic model with long-range interactions, the Hamiltonian mean-field model. Using molecular dynamics simulations, we show that when the system settles into an inhomogeneous quasi-stationary state in which the velocity distribution has suprathermal tails, the temperature and density profiles are anticorrelated: denser parts of the system are colder than dilute ones. We argue that this may be a generic property of long-range interacting systems.

The temperature in the crust of an accreting neutron star, which comprises its outermost kilometer, is set by heating from nuclear reactions at large densities, neutrino cooling, and heat transport from the interior. The heated crust has been thought to affect observable phenomena at shallower depths, such as thermonuclear bursts in the accreted envelope. Here we report that cycles of electron capture and its inverse, $\beta^-$ decay, involving neutron-rich nuclei at a typical depth of about 150 m, cool the outer neutron star crust by emitting neutrinos while also thermally decoupling the surface layers from the deeper crust. This Urca mechanism has been studied in the context of white dwarfs and Type Ia supernovae, but hitherto was not considered in neutron stars, because previous models computed the crust reactions using a zero-temperature approximation and assumed that only a single nuclear species was present at any given depth. This thermal decoupling means that X-ray bursts and other surface phenomena are largely independent of the strength of deep crustal heating. The unexpectedly short recurrence times, of the order of years, observed for very energetic thermonuclear superbursts are therefore not an indicator of a hot crust, but may point instead to an unknown local heating mechanism near the neutron star surface.

We seek to identify the primary agents causing Forbush decreases (FDs) observed at the Earth in high rigidity cosmic rays. In particular, we ask if such FDs are caused mainly by coronal mass ejections (CMEs) from the Sun that are directed towards the Earth, or by their associated shocks. We use the muon data at cutoff rigidities ranging from 14 to 24 GV from the GRAPES-3 tracking muon telescope to identify FD events. We select those FD events that have a reasonably clean profile, and can be reasonably well associated with an Earth-directed CME and its associated shock. We employ two models: one that considers the CME as the sole cause of the FD (the CME-only model) and one that considers the shock as the only agent causing the FD (the shock-only model). We use an extensive set of observationally determined parameters for both these models. The only free parameter in these models is the level of MHD turbulence in the sheath region, which mediates cosmic ray diffusion (into the CME, for the CME-only model and across the shock sheath, for the shock-only model). We find that good fits to the GRAPES-3 multi-rigidity data using the CME-only model require turbulence levels in the CME sheath region that are only slightly higher than those estimated for the quiet solar wind. On the other hand, reasonable model fits with the shock-only model require turbulence levels in the sheath region that are an order of magnitude higher than those in the quiet solar wind. This observation naturally leads to the conclusion that the Earth-directed CMEs are the primary contributors to FDs observed in high rigidity cosmic rays.

This paper presents a systematic evaluation of the ability of theoretical models to reproduce experimental Gamow-Teller transition strength distributions measured via (n,p)-type charge-exchange reactions at intermediate beam energies. The focus is on transitions from stable nuclei in the pf shell (45<A<64). The impact of deviations between experimental and theoretical Gamow-Teller strength distributions on derived stellar electron-capture rates at densities and temperatures of relevance for Type Ia and Type II supernovae is investigated. The theoretical models included in the study are based on the shell-model, using the KB3G and GXPF1a interactions, and quasiparticle random-phase approximation (QRPA) using ground-state deformation parameters and masses from the finite-range droplet model.

The location of electron capture heat sources in the crust of accreting neutron stars depends on the masses of extremely neutron-rich nuclei. We present first results from a new implementation of the time-of-flight technique to measure nuclear masses of rare isotopes at the National Superconducting Cyclotron Laboratory. The masses of 16 neutron-rich nuclei in the scandium -- nickel range were determined simultaneously, improving the accuracy compared to previous data in 12 cases. The masses of $^{61}${V}, $^{63}${Cr}, $^{66}${Mn}, and $^{74}${Ni} were measured for the first time with mass excesses of $-30.510(890)$ MeV, $-35.280(650)$ MeV, $-36.900(790)$ MeV, and $-49.210(990)$ MeV, respectively. With the measurement of the $^{66}$Mn mass, the locations of the two dominant electron capture heat sources in the outer crust of accreting neutron stars that exhibit superbursts are now experimentally constrained. We find that the location of the $^{66}$Fe$\rightarrow^{66}$Mn electron capture transition occurs significantly closer to the surface than previously assumed because our new experimental Q-value is 2.1 MeV (2.6$\sigma$) smaller than predicted by the FRDM mass model.

Two new statistics, namely $\Delta_\chi^2$ and $\Delta_\chi$, based on extreme value theory, were derived in \cite{gupta08,gupta10}. We use these statistics to study direction dependence in the HST key project data which provides the most precise measurement of the Hubble constant. We also study the non-Gaussianity in this data set using these statistics. Our results for $\Delta_\chi^2$ show that the significance of direction dependent systematics is restricted to well below one $\sigma$ confidence limit, however, presence of non-Gaussian features is subtle. On the other hand $\Delta_\chi$ statistic, which is more sensitive to direction dependence, shows direction dependence systematics to be at slightly higher confidence level, and the presence of non-Gaussian features at a level similar to the $\Delta_\chi^2$ statistic.

We present a catalog of compact sources derived from the QUaD Galactic Plane Survey. The survey covers ~800 square degrees of the inner galaxy (|b|<4 degrees) in Stokes I, Q and U parameters at 100 and 150 GHz, with angular resolution 5 and 3.5 arcminutes. 505 unique sources are identified in I, of which 239 are spatially matched between frequency bands, with 50 (216) detected at 100 (150) GHz alone; 182 sources are identified as ultracompact HII (UCHII) regions. Approximating the distribution of total intensity source fluxes as a power-law, we find a slope of $\gamma_{S,100}=-1.8\pm0.4$ at 100 GHz, and $\gamma_{S,150}=-2.2\pm0.4$ at 150 GHz. Similarly, the power-law index of the source two-point angular correlation function is $\gamma_{\theta,100}=-1.21\pm0.04$ and $\gamma_{\theta,150}=-1.25\pm0.04$. The total intensity spectral index distribution peaks at $\alpha_{I}\sim0.25$, indicating that dust emission is not the only source of radiation produced by these objects between 100 and 150 GHz; free-free radiation is likely significant in the 100 GHz band. Four sources are detected in polarized intensity P, of which three have matching counterparts in I. Three of the polarized sources lie close to the galactic center, Sagittarius A*, Sagittarius B2 and the Galactic Radio Arc, while the fourth is RCW 49, a bright HII region. An extended polarized source, undetected by the source extraction algorithm on account of its $\sim0.5^{\circ}$ size, is identified visually, and is an isolated example of large-scale polarized emission oriented distinctly from the bulk galactic dust polarization.

We present a survey of ~800 square degrees of the galactic plane observed with the QUaD telescope. The primary product of the survey are maps of Stokes I, Q and U parameters at 100 and 150 GHz, with spatial resolution 5 and 3.5 arcminutes respectively. Two regions are covered, spanning approximately 245-295 and 315-5 degrees in galactic longitude l, and -4<b<+4 degrees in galactic latitude b. At 0.02 degree square pixel size, the median sensitivity is 74 and 107 kJy/sr at 100 GHz and 150 GHz respectively in I, and 98 and 120 kJy/sr for Q and U. In total intensity, we find an average spectral index of 2.35+/-0.01 (stat) +/- 0.02 (sys) for |b|<1 degree, indicative of emission components other than thermal dust. A comparison to published dust, synchrotron and free-free models implies an excess of emission in the 100 GHz QUaD band, while better agreement is found at 150 GHz. A smaller excess is observed when comparing QUaD 100 GHz data to WMAP 5-year W band; in this case the excess is likely due to the wider bandwidth of QUaD. Combining the QUaD and WMAP data, a two-component spectral fit to the inner galactic plane (|b|<1 degree) yields mean spectral indices of -0.32 +/- 0.03 and 2.84 +/- 0.03; the former is interpreted as a combination of the spectral indices of synchrotron, free-free and dust, while the second is attributed largely to the thermal dust continuum. In the same galactic latitude range, the polarization data show a high degree of alignment perpendicular to the expected galactic magnetic field direction, and exhibit mean polarization fraction 1.38+/-0.08 (stat) +/-0.1 (sys) % at 100 GHz and 1.70 +/- 0.06 (stat) +/- 0.1 (sys) % at 150GHz. We find agreement in polarization fraction between QUaD 100 GHz and WMAP W band, the latter giving 1.1 +/- 0.4 %.

We present results of five years of optical (UBVRI) observations of the very-high-energy gamma-ray blazar 1ES 1011+496 at the MDM Observatory. We calibrated UBVRI magnitudes of five comparison stars in the field of the object. Most of our observations were done during moderately faint states of 1ES 1011+496 with R > 15.0. The light curves exhibit moderate, closely correlated variability in all optical wavebands on time scales of a few days. A cross-correlation analysis between optical bands does not show significant evidence for time lags. We find a positive correlation (Pearson's r = 0.57; probability of non-correlation P(>r) ~ 4e-8) between the R-band magnitude and the B - R color index, indicating a bluer-when-brighter trend. Snap-shot optical spectral energy distributions (SEDs) exhibit a peak within the optical regime, typically between the V and B bands. We find a strong (r = 0.78; probability of non-correlation P (>r) ~ 1e-15) positive correlation between the peak flux and the peak frequency, best fit by a relation $\nu F_{\nu}^{\rm pk} \propto \nu_{\rm pk}^k$ with k = 2.05 +/- 0.17. Such a correlation is consistent with the optical (synchrotron) variability of 1ES 1011+496 being primarily driven by changes in the magnetic field.

Centaurus (Cen) A represents one of the best candidates for an isolated, compact, highly polarized source that is bright at typical cosmic microwave background (CMB) experiment frequencies. We present measurements of the 4 degree by 2 degree region centered on Cen A with QUaD, a CMB polarimeter whose absolute polarization angle is known to 0.5 degrees. Simulations are performed to assess the effect of misestimation of the instrumental parameters on the final measurement and systematic errors due to the field's background structure and temporal variability from Cen A's nuclear region are determined. The total (Q, U) of the inner lobe region is (1.00 +/- 0.07 (stat.) +/- 0.04 (sys.), -1.72 +/- 0.06 +/- 0.05) Jy at 100 GHz and (0.80 +/- 0.06 +/- 0.06, -1.40 +/- 0.07 +/- 0.08) Jy at 150 GHz, leading to polarization angles and total errors of -30.0 +/- 1.1 degrees and -29.1 +/- 1.7 degrees. These measurements will allow the use of Cen A as a polarized calibration source for future millimeter experiments.

We have gathered optical photometry data from the literature on a large sample of Swift-era gamma-ray burst (GRB) afterglows including GRBs up to September 2009, for a total of 76 GRBs, and present an additional three pre-Swift GRBs not included in an earlier sample. Furthermore, we publish 840 additional new photometry data points on a total of 42 GRB afterglows, including large data sets for GRBs 050319, 050408, 050802, 050820A, 050922C, 060418, 080413A and 080810. We analyzed the light curves of all GRBs in the sample and derived spectral energy distributions for the sample with the best data quality, allowing us to estimate the host galaxy extinction. We transformed the afterglow light curves into an extinction-corrected z=1 system and compared their luminosities with a sample of pre-Swift afterglows. The results of a former study, which showed that GRB afterglows clustered and exhibited a bimodal distribution in luminosity space, is weakened by the larger sample. We found that the luminosity distribution of the two afterglow samples (Swift-era and pre-Swift) are very similar, and that a subsample for which we were not able to estimate the extinction, which is fainter than the main sample, can be explained by assuming a moderate amount of line-of-sight host extinction. We derived bolometric isotropic energies for all GRBs in our sample, and found only a tentative correlation between the prompt energy release and the optical afterglow luminosity at one day after the GRB in the z=1 system. A comparative study of the optical luminosities of GRB afterglows with echelle spectra (which show a high number of foreground absorbing systems) and those without reveals no indication that the former are statistically significantly more luminous. (abridged)

We evaluate the contribution of cosmic microwave background (CMB) polarization spectra to cosmological parameter constraints. We produce cosmological parameters using high-quality CMB polarization data from the ground-based QUaD experiment and demonstrate for the majority of parameters that there is significant improvement on the constraints obtained from satellite CMB polarization data. We split a multi-experiment CMB dataset into temperature and polarization subsets and show that the best-fit confidence regions for the LCDM 6-parameter cosmological model are consistent with each other, and that polarization data reduces the confidence regions on all parameters. We provide the best limits on parameters from QUaD EE/BB polarization data and we find best-fit parameters from the multi-experiment CMB dataset using the optimal pivot scale of k_p=0.013 Mpc-1 to be {omch2, ombh2, H_0, A_s, n_s, tau}= {0.113, 0.0224, 70.6, 2.29 times 10^-9, 0.960, 0.086}.

We revise and extend the extreme value statistic, introduced in \cite{gup08}, to study directional dependence in the high redshift supernova data; arising either from departures from the cosmological principle or due to direction dependent statistical systematics in e data. We introduce a likelihood function that analytically marginalises over the Hubble constant, and use it to extend our previous statistic. We also introduce a new statistic that is sensitive to direction dependence arising from living off-centre inside a large void, as well as previously mentioned reasons for anisotropy. We show that for large data sets this statistic has a limiting form that can be computed analytically. We apply our statistics to the gold data sets from \cite{rie04} and \cite{rie07}, as in our previous work. Our revision and extension of previous statistic shows that 1) the effect of marginalsing over Hubble constant instead of using its best fit value has only a marginal effect on our results. However, correction of errors in our previous work reduce the level of non-Gaussianity in the 2004 gold data that was found in our earlier work. The revised results for the 2007 gold data show that the data is consistent with isotropy and Gaussianity. Our second statistic confirms these results.

We present an improved analysis of the final dataset from the QUaD experiment. Using an improved technique to remove ground contamination, we double the effective sky area and hence increase the precision of our CMB power spectrum measurements by ~30% versus that previously reported. In addition, we have improved our modeling of the instrument beams and have reduced our absolute calibration uncertainty from 5% to 3.5% in temperature. The robustness of our results is confirmed through extensive jackknife tests and by way of the agreement we find between our two fully independent analysis pipelines. For the standard 6-parameter LCDM model, the addition of QUaD data marginally improves the constraints on a number of cosmological parameters over those obtained from the WMAP experiment alone. The impact of QUaD data is significantly greater for a model extended to include either a running in the scalar spectral index, or a possible tensor component, or both. Adding both the QUaD data and the results from the ACBAR experiment, the uncertainty in the spectral index running is reduced by ~25% compared to WMAP alone, while the upper limit on the tensor-to-scalar ratio is reduced from r < 0.48 to r < 0.33 (95% c.l). This is the strongest limit on tensors to date from the CMB alone. We also use our polarization measurements to place constraints on parity violating interactions to the surface of last scattering, constraining the energy scale of Lorentz violating interactions to < 1.5 x 10^{-43} GeV (68% c.l.). Finally, we place a robust upper limit on the strength of the lensing B-mode signal. Assuming a single flat band power between l = 200 and l = 2000, we constrain the amplitude of B-modes to be < 0.57 micro-K^2 (95% c.l.).

The thermal conductivity of a dense {\it multi-component} plasma is critical to the modeling of accreting neutron stars. To this end, we perform large-scale molecular dynamics simulations to calculate the static structure factor of the dense multi-component plasma in the neutron star crust from near the photosphere-ocean boundary to the vicinity of the neutron drip point. The structure factors are used to validate a microscopic linear mixing rule that is valid for arbitrarily complex plasmas over a wide range of Coulomb couplings. The microscopic mixing rule in turn implies and validates the linear mixing rule for the equation of state properties and also the linear mixing rule for the electrical and thermal conductivities of dense multi-component plasmas. To make our result as useful as possible, for the specific cases of electrical and thermal conductivities, we provide a simple analytic fit that is valid for arbitrarily complex multi-component plasmas over a wide range of Coulomb couplings. We compute the thermal conductivity for a representative compositional profile of the outer crust of an accreting neutron star in which hundreds of nuclear species can be present. We utilize our results to re-examine the so-called impurity parameter formalism as used to characterize impure plasmas.

We present measurements of the cosmic microwave background (CMB) radiation temperature anisotropy in the multipole range 2000<ell<3000 from the QUaD telescope's second and third observing seasons. After masking the brightest point sources our results are consistent with the primary LCDM expectation alone. We estimate the contribution of residual (un-masked) radio point sources using a model calibrated to our own bright source observations, and a full simulation of the source finding and masking procedure. Including this contribution slightly improves the chi^2. We also fit a standard SZ template to the bandpowers and see no strong evidence of an SZ contribution, which is as expected for sigma_8 approx 0.8.

A new instrument with a dual-beam H-alpha Doppler system is being developed at the Udaipur Solar Observatory (USO) in order to improve the quality and quantity of data on quiet, activated and erupting filaments and prominences on the Sun, especially those associated with geo-effective coronal mass ejections. These data can be potentially used to construct three-dimensional topology of erupting filaments as they leave the surface of the Sun and can be compared with multi-wavelength data obtained from space missions such as STEREO, SOHO, and Hinode. The characterization of various optical components for the instrument is being carried out, and some preliminary results are described in the paper.

We describe the QUaD experiment, a millimeter-wavelength polarimeter designed to observe the Cosmic Microwave Background (CMB) from a site at the South Pole. The experiment comprises a 2.64 m Cassegrain telescope equipped with a cryogenically cooled receiver containing an array of 62 polarization-sensitive bolometers. The focal plane contains pixels at two different frequency bands, 100 GHz and 150 GHz, with angular resolutions of 5 arcmin and 3.5 arcmin, respectively. The high angular resolution allows observation of CMB temperature and polarization anisotropies over a wide range of scales. The instrument commenced operation in early 2005 and collected science data during three successive Austral winter seasons of observation.

We constrain parity-violating interactions to the surface of last scattering using spectra from the QUaD experiment's second and third seasons of observations by searching for a possible systematic rotation of the polarization directions of CMB photons. We measure the rotation angle due to such a possible "cosmological birefringence" to be 0.55 deg. +/- 0.82 deg. (random) +/- 0.5 deg. (systematic) using QUaD's 100 and 150 GHz TB and EB spectra over the multipole range 200 < l < 2000, consistent with null, and constrain Lorentz violating interactions to < 2^-43 GeV (68% confidence limit). This is the best constraint to date on electrodynamic parity violation on cosmological scales.

We report results from the second and third seasons of observation with the QUaD experiment. Angular power spectra of the Cosmic Microwave Background are derived for both temperature and polarization at both 100 GHz and 150 GHz, and as cross frequency spectra. All spectra are subjected to an extensive set of jackknife tests to probe for possible systematic contamination. For the implemented data cuts and processing technique such contamination is undetectable. We analyze the difference map formed between the 100 and 150 GHz bands and find no evidence of foreground contamination in polarization. The spectra are then combined to form a single set of results which are shown to be consistent with the prevailing LCDM model. The sensitivity of the polarization results is considerably better than that of any previous experiment -- for the first time multiple acoustic peaks are detected in the E-mode power spectrum at high significance.

In this paper we present a parameter estimation analysis of the polarization and temperature power spectra from the second and third season of observations with the QUaD experiment. QUaD has for the first time detected multiple acoustic peaks in the E-mode polarization spectrum with high significance. Although QUaD-only parameter constraints are not competitive with previous results for the standard 6-parameter LCDM cosmology, they do allow meaningful polarization-only parameter analyses for the first time. In a standard 6-parameter LCDM analysis we find the QUaD TT power spectrum to be in good agreement with previous results. However, the QUaD polarization data shows some tension with LCDM. The origin of this 1 to 2 sigma tension remains unclear, and may point to new physics, residual systematics or simple random chance. We also combine QUaD with the five-year WMAP data set and the SDSS Luminous Red Galaxies 4th data release power spectrum, and extend our analysis to constrain individual isocurvature mode fractions, constraining cold dark matter density, alpha(cdmi)<0.11 (95 % CL), neutrino density, alpha(ndi)<0.26 (95 % CL), and neutrino velocity, alpha(nvi)<0.23 (95 % CL), modes. Our analysis sets a benchmark for future polarization experiments.

We seek to estimate the average level of MHD turbulence near coronal mass ejection (CME) fronts as they propagate from the Sun to the Earth. We examine the cosmic ray data from the GRAPES-3 tracking muon telescope at Ooty, together with the data from other sources for three well observed Forbush decrease events. Each of these events are associated with frontside halo Coronal Mass Ejections (CMEs) and near-Earth magnetic clouds. In each case, we estimate the magnitude of the Forbush decrease using a simple model for the diffusion of high energy protons through the largely closed field lines enclosing the CME as it expands and propagates from the Sun to the Earth. We use estimates of the cross-field diffusion coefficient $D_{\perp}$ derived from published results of extensive Monte Carlo simulations of cosmic rays propagating through turbulent magnetic fields. Our method helps constrain the ratio of energy density in the turbulent magnetic fields to that in the mean magnetic fields near the CME fronts. This ratio is found to be $\sim$ 2% for the 11 April 2001 Forbush decrease event, $\sim$ 6% for the 20 November 2003 Forbush decrease event and $\sim$ 249% for the much more energetic event of 29 October 2003.

In our calculation of neutron star crust heating we include several key new model features. In earlier work electron capture (EC) only allowed neutron emission from the daughter ground-state; here we calculate, in a deformed QRPA model, EC decay rates to all states in the daughter that are allowed by Gamow-Teller selection rules and energetics. The subsequent branching ratios between the 1n,...,xn channels and the competing $\gamma$-decay are calculated in a Hauser-Feshbach model. Since EC accesses excited states, many more neutrons are emitted in our calculation than in previous work, leading to accelerated reaction flows. In our multi-component plasma model a single (EC,xn) reaction step can produce several neutron-deficient nuclei, each of which can further decay by (EC,xn). Hence, the neutron emission occurs more continuously with increasing depth as compared to that in a one-component plasma model.

We derive upper limits on the ratio $f_{GRB/CCSN}(z) \equiv R_{GRB}(z)/R_{CCSN}(z) \equiv f_{GRB/CCSN}(0)(1+z)^\alpha$, the ratio of the rate, $R_{GRB}$, of long-duration Gamma Ray Bursts (GRBs) to the rate, $R_{CCSN}$, of core-collapse supernovae (CCSNe) in the Universe ($z$ being the cosmological redshift and $\alpha\geq 0$), by using the upper limit on the diffuse TeV--PeV neutrino background given by the AMANDA-II experiment in the South Pole, under the assumption that GRBs are sources of TeV--PeV neutrinos produced from decay of charged pions produced in $p\gamma$ interaction of protons accelerated to ultrahigh energies at internal shocks within GRB jets. For the assumed ``concordance model'' of cosmic star formation rate, $R_{SF}$, with $R_{CCSN}(z) \propto R_{SF}(z)$, our conservative upper limits are $f_{GRB/CCSN}(0)\leq 5.0\times10^{-3}$ for $\alpha=0$, and $f_{GRB/CCSN}(0)\leq 1.1\times10^{-3}$ for $\alpha=2$, for example. These limits are already comparable to (and, for $\alpha\geq 1$ already more restrictive than) the current upper limit on this ratio inferred from other astronomical considerations, thus providing a useful independent probe of and constraint on the CCSN-GRB connection. Non-detection of a diffuse TeV--PeV neutrino background by the up-coming IceCube detector in the South pole after three years of operation, for example, will bring down the upper limit on $f_{GRB/CCSN}(0)$ to below few $\times10^{-5}$ level, while a detection will confirm the hypothesis of proton acceleration to ultrahigh energies in GRBs and will potentially also yield the true rate of occurrence of these events in the Universe.

We perform a self-consistent calculation of the thermal structure in the crust of a superbursting neutron star. In particular, we follow the nucleosynthetic evolution of an accreted fluid element from its deposition into the atmosphere down to a depth where the electron Fermi energy is 20 MeV. We include temperature-dependent continuum electron capture rates and realistic sources of heat loss by thermal neutrino emission from the crust and core. We show that, in contrast to previous calculations, electron captures to excited states and subsequent gamma-emission significantly reduce the local heat loss due to weak-interaction neutrinos. Depending on the initial composition these reactions release up to a factor of 10 times more heat at densities < 10^{11} g/cc than obtained previously. This heating reduces the ignition depth of superbursts. In particular, it reduces the discrepancy noted by Cumming et al. between the temperatures needed for unstable 12C ignition on timescales consistent with observations and the reduction in crust temperature from Cooper pair neutrino emission.

The most detailed constraints on the accelerating expansion of the universe and details of nature of dark energy are derived from the high redshift supernova data, assuming that the errors in the measurements are Gaussian in nature. There is a possibility that there are direction dependent systematics in the data, either due to uncorrected, known physical processes or because there are tiny departures from the cosmological principle, making the universe slightly anisotropic. To investigate this possibility we introduce a statistic based on the extreme value theory and apply it to the gold data set from Riess et al. (2004). Our analysis indicates a systematic, direction dependent non-gaussianity at about one sigma level.

The Galactic high-mass X-ray binary and jet source (microquasar) LSI +61 303 has recently been detected at TeV gamma-ray energies by the MAGIC telescope. We have applied a time-dependent leptonic jet model to the broadband spectral energy distribution and suggested (though not unambiguously detected) orbital modulation of the very high energy gamma-ray emission of this source. Our model takes into account time dependent electron injection and acceleration, and the adiabatic and radiative cooling of non-thermal electrons. It includes synchrotron, synchrotron self-Compton and external inverse Compton (with seed photons from the companion star), as well as gamma-gamma absorption of gamma-rays by starlight photons. The model can successfully reproduce the available multiwavelength observational data. Our best fit to the SED indicates that a magnetic field of B_0 ~ 5 X 10^3 G at ~ 10^3 R_g is required, and electrons need to be accelerated out to TeV energies (gamma_2 = 10^6) with a nonthermal injection spectrum with a spectral index of q = 1.7, indicating the operation of acceleration mechanisms beyond the standard first-order Fermi mechanism at relativistic or non-relativistic shocks. The orbital modulation of the VHE gamma-ray emission can be explained solely by the geometrical effect of changes in the relative orientation of the stellar companion with respect to the compact object and jet as it impacts the position and depth of the gamma-gamma absorption trough. Such a scenario predicts a trend of spectral hardening during VHE gamma-ray low orbital phases.

We calculate the general solutions for a warm inflationary scenario with weak dissipation, reviewing the dissipative dynamics of the two-fluid system, and calculate the bispectrum of the gravitational field fluctuations generated in the case where dissipation of the vacuum potential during inflation is the mechanism for structure formation, but is the sub-dominant effect in the dynamics of the scalar field during slow-roll. The bispectrum is non-zero because of the self-interaction of the scalar field. We compare the predictions with both those of standard, or `supercooled', inflationary models, and warm inflation models with strong dissipation and consider the detectability of these levels of non-Gaussianity in the bispectrum of the cosmic microwave background. We find that the levels of non--Gaussianity for warm and supercooled inflation are an order of magnitude different.

Jet models for the high-energy emission of Galactic X-ray binary sources have regained significant interest with detailed spectral and timing studies of the X-ray emission from microquasars, the recent detection by the HESS collaboration of very-high-energy gamma-rays from the microquasar LS~5039, and the earlier suggestion of jet models for ultraluminous X-ray sources observed in many nearby galaxies. Here we study the synchrotron and Compton signatures of time-dependent electron injection and acceleration, adiabatic and radiative cooling, and different jet geometries in the jets of Galactic microquasars. Synchrotron, synchrotron-self-Compton, and external-Compton radiation processes with soft photons provided by the companion star and the accretion disk are treated. An analytical solution is presented to the electron kinetic equation for general power-law geometries of the jets for Compton scattering in the Thomson regime. We pay particular attention to predictions concerning the rapid flux and spectral variability signatures expected in a variety of scenarios, making specific predictions concerning possible spectral hysteresis, similar to what has been observed in several TeV blazars. Such predictions should be testable with dedicated monitoring observations of Galactic microquasars and ultraluminous X-ray sources using Chandra and/or XMM-Newton.

We show how the primordial bispectrum of density perturbations from inflation may be characterised in terms of manifestly gauge-invariant cosmological perturbations at second order. The primordial metric perturbation, zeta, describing the perturbed expansion of uniform-density hypersurfaces on large scales is related to scalar field perturbations on unperturbed (spatially-flat) hypersurfaces at first- and second-order. The bispectrum of the metric perturbation is thus composed of (i) a local contribution due to the second-order gauge-transformation, and (ii) the instrinsic bispectrum of the field perturbations on spatially flat hypersurfaces. We generalise previous results to allow for scale-dependence of the scalar field power spectra and correlations that can develop between fields on super-Hubble scales.

We present a search for non-Gaussianity in the WMAP first-year data using the two-point correlation function of maxima and minima in the temperature map. We find evidence for non-Gaussianity on large scales, whose origin appears to be associated with unsubstracted foregrounds, but which is not entirely clear. The signal appears to be associated most strongly with cold spots, and is more pronounced in the Southern galactic hemisphere. Removal of the region of sky near the galactic plane, or filtering out large-scale modes removes the signal. Analysis of individual frequency maps shows strongest signal in the 41GHz Q band. A study of difference maps tests the hypothesis that the non-Gaussianity is due to residual foregrounds and noise, but shows no significant detection. We suggest that the detection is due to large-scale residual foregrounds affecting more than one frequency band, but a primordial contribution from the Cosmic Microwave Background cannot be excluded.

We investigate the case of a homogeneous tachyon field coupled to gravity in a spatially flat Friedman-Robertson-Walker spacetime. Assuming the field evolution to be exponentially decaying with time we solve the field equations and show that, under certain conditions, the scale factor represents an accelerating universe, following a phase of decelerated expansion. We make use of a model of dark energy (with p=-\rho) and dark matter (p=0) where a single scalar field (tachyon) governs the dynamics of both the dark components. We show that this model fits the current supernova data as well as the canonical \LambdaCDM model. We give the bounds on the parameters allowed by the current data.

We study the evolution of the cosmological perturbations after inflation in curvaton models where the non-relativistic curvaton decays into both radiation and a cold dark matter component. We calculate the primordial curvature and correlated isocurvature perturbations inherited by the radiation and cold dark matter after the curvaton has decayed. We give the transfer coefficient in terms of the initial curvaton density relative to the curvaton decay rate.

The generic prediction for the perturbations generated during slow--roll, single--field inflation, as they appear in the cosmic microwave background (CMB), is a flat, close--to--Gaussian spectrum. We calculate the general solution for a warm inflationary scenario with weak dissipation, reviewing the dissipative dynamics of the two--fluid system, and calculate the bispectrum of the gravitational field fluctuations generated in the case where dissipation of the vacuum potential during inflation is the mechanism for structure formation, but is the sub--dominant effect in the dynamics of the scalar field.

We calculate the bispectrum of the gravitational field fluctuations generated during warm inflation, where dissipation of the vacuum potential during inflation is the mechanism for structure formation. The bispectrum is non--zero because of the self--interaction of the scalar field. We compare the predictions with those of standard, or `supercooled', inflationary models, and consider the detectability of these levels of non--Gaussianity in the bispectrum of the cosmic microwave background. We find that the levels of non--Gaussianity for warm and supercooled inflation are comparable, and over--ridden by the contribution to the bispectrum due to other physical effects. We also conclude that the resulting bispectrum values will be undetectable in the cosmic microwave background for both the MAP and Planck Surveyor satellites.

The statistical properties of a map of the primary fluctuations in the cosmic microwave background (CMB) may be specified to high accuracy by a few thousand power spectra measurements, provided the fluctuations are gaussian, yet the number of parameters relevant for the CMB is probably no more than about 10-20. There is consequently a large degree of redundancy in the power spectrum data. In this paper, we show that the MOPED data compression technique can reduce the CMB power spectrum measurements to about 10-20 numbers (one for each parameter), from which the cosmological parameters can be estimated virtually as accurately as from the complete power spectrum. This offers opportunities for very fast parameter estimation from real and simulated CMB skies, with accurate likelihood calculations at Planck resolution being speeded up by a factor of around five hundred million.

We present the analysis of 65 hours of high speed photometric observations of HD 12098 taken from State Observatory, Naini Tal and Gurushikhar Observatory, Mt.Abu on sixteen nights spanning from November 1999 to November 2000. HD 12098 is the first rapidly oscillating Ap star discovered from the `Naini Tal-Cape survey for northern hemisphere roAp stars'. It is the 32nd in the complete list. HD 12098 exhibits one predominant mode of oscillation at nu_1 = 2.1738 mHz. The second-most significant frequency in our data is at nu_2 = 2.1641 mHz with a 1 cycle/day alias ambiguity. We argue that nu_2 is a rotational sidelobe of nu_1, rather than an independent pulsation mode. Evidence for the presence of two other frequencies at 2.1807 and 2.3056 mHz is also presented.

We compute precise predictions for the two-point correlation function of local maxima (or minima) in the temperature of the microwave background, under the assumption that it is a random gaussian field. For a given power spectrum and peak threshold there are no adjustable parameters, and since this analysis does not make the small-angle approximation of Heavens & Sheth (1999), it is essentially complete. We find oscillatory features which are absent in the temperature autocorrelation function, and we also find that the small-angle approximation to the peak-peak correlation function is accurate to better than 0.01 on all scales. These high-precision predictions can form the basis of a sensitive test of the gaussian hypothesis with upcoming all-sky microwave background experiments MAP and Planck, affording a thorough test of the inflationary theory of the early Universe. To illustrate the effectiveness of the technique, we apply it to simulated maps of the microwave sky arising from the cosmic string model of structure formation, and compare with the bispectrum as a non-gaussian discriminant. We also show how peak statistics can be a valuable tool in assessing and statistically removing contamination of the map by foreground point sources.