Recent detection of the neutrino events IceCube-170922A, 13 muon-neutrino events observed in 2014-2015 and IceCube-141209A by IceCube observatory from the Blazars, namely TXS 0506+056, PKS 0502+049/TXS 0506+056 and GB6 J1040+0617 respectively in the state of enhanced gamma-ray emission, indicates the acceleration of cosmic rays in the blazar jets. The photo-meson ($p\gamma$) interaction cannot explain the IceCube observations of 13 neutrino events. The non-detection of broadline emission in the optical spectra of the IceCube blazars, however, question the hadronuclear (pp) interaction interpretation through relativistic jet meets with high density cloud. In this work, we investigate the proton blazar model in which the non-relativistic protons that come into existence under the charge neutrality condition of the blazar jet can offer sufficient target matter for $pp$ interaction with shock-accelerated protons, to describe the observed high-energy gamma-rays and neutrino signal from the said blazars. Our findings suggest that the model can explain consistently the observed electromagnetic spectrum in combination with appropriate number of neutrino events from the corresponding blazars.

Recent detection of the neutrino event, IceCube-170922A by IceCube observatory from the Blazar TXS 0506+056 in the state of enhanced gamma ray emission indicates for acceleration of cosmic rays in the blazar jet. The non-detection of the broadline emission in the optical spectrum of TXS 0506+056 and other BL Lac objects suggests that external photons emissions are weak and hence photo-meson (p-gamma) interaction may not be a favored mechanism for high energy neutrino production. The lack of broadline signatures also creates doubt about the presence of a high density cloud in the vicinity of the super-massive black hole (SMBH) of TXS 0506+056 and consequently raised question on hadronuclear (pp) interaction interpretation like relativistic jet meets with high density cloud. Here we demonstrate that non-relativistic protons in the proton blazar model, those come into existence under charge neutrality condition of the blazar jet, offer sufficient target matter for pp-interaction with shock accelerated protons and consequently the model can describe consistently the observed high energy gamma rays and neutrino signal from the blazar TXS 0506+056.

The galactic cosmic rays are generally believed to be originated in supernova remnants (SNRs), produced in diffusive shock acceleration (DSA) process in supernova blast waves driven by expanding SNRs. One of the key unsettled issue in SNR origin of cosmic ray model is the maximum attainable energy by a cosmic ray particle in the supernova shock. Recently it has been suggested that an amplification of effective magnetic field strength at the shock may take place in young SNRs due to growth of magnetic waves induced by accelerated cosmic rays and as a result the maximum energy achieved by cosmic rays in SNR may reach the knee energy instead of $\sim 200$ TeV as predicted earlier under normal magnetic field situation. In the present work we investigate the implication of such maximum energy scenarios on TeV gamma rays and neutrino fluxes from the molecular clouds interacting with the SNR W28. The authors compute the gamma-ray and neutrino flux assuming two different values for the maximum energy reached by cosmic rays in the SNR, from CR interaction in nearby molecular clouds. Both protons and nuclei are considered as accelerated particles and as target material. Our findings suggest that the issue of the maximum energy of cosmic rays in SNRs can be observationally settled by the upcoming gamma-ray experiment the Large High Altitude Air Shower Observatory (LHAASO). The estimated neutrino fluxes from the molecular clouds are , however, out of reach of the present/near future generation of neutrino telescopes.

Space-time geometry of halo region in spiral galaxies is derived using the observed flat galactic rotation curve feature and considering the presence of cold dark matter in the galaxy. Gravitational deflection angle due to the derived space time metric has been evaluated and implications of the findings are discussed. Exploiting cosmological observations about the dark matter content of the universe relative to the luminous matter, the radial extent of spiral galaxies, such as that of Milky Way, has been predicted which may be observationally tested in future.

The standard perception is that the detection of high energy (TeV energies and above) neutrinos from an astrophysical object is a conclusive evidence for the presence of hadronic cosmic rays at the source. In the present work we demonstrate that TeV neutrinos can also be originated from energetic electrons via electromagnetic interactions in different potential cosmic ray sources with flux levels comparable to that of the hadronic originated neutrinos at high energies. Our findings thus imply that at least a part of the neutrinos observed by Icecube observatory may be originated from energetic electrons. The present analysis further suggests that only a combine study of TeV gamma rays and neutrinos over a wide energy range from an astrophysical object can unambiguously identify the nature of their parents, hadrons or leptons.

It is widely believe that galactic cosmic rays are originated in supernova remnants (SNRs) where they are accelerated by diffusive shock acceleration process at supernova blast waves driven by expanding SNRs. In recent theoretical developments of the diffusive shock acceleration theory in SNRs, protons are expected to accelerate in SNRs at least up to the knee energy. If SNRs are true generator of cosmic rays, they should accelerate not only protons but also heavier nuclei with right proportion and the maximum energy of heavier nuclei should be atomic mass (Z) times that of protons. In this work we investigate the implications of acceleration of heavier nuclei in SNRs on energetic gamma rays those are produced in hadronic interaction of cosmic rays with ambient matter. Our findings suggest that the energy conversion efficiency has to be nearly double for the mixed cosmic ray composition instead of pure protons to explain the observation and secondly the gamma ray flux above few tens of TeV would be significantly higher if cosmic rays particles can attain energies Z times of the knee energy in lieu of 200 TeV, as suggested earlier for non-amplified magnetic fields. The two stated maximum energy paradigm will be discriminated in future by the upcoming gamma ray experiments like Cherenkov Telescope array (CTA).

The possibility of estimating the mass composition of primary cosmic rays above the knee of its energy spectrum through the study of high energy gamma rays, muons and neutrinos produced in the interactions of cosmic rays with the solar ambient matter and radiation has been explored. It is found that the theoretical fluxes of TeV gamma rays, muons and neutrinos from a region around $15^{o}$ of the Sun are sensitive to mass composition of cosmic rays in the PeV energy range. The experimental prospects for detection of such TeV gamma rays/neutrinos by future experiments are discussed.

We construct a self-consistent relativistic Newtonian analogue corresponding to gravitational static spherical symmetric spacetime geometries, staring directly from a generalized scalar relativistic gravitational action in Newtonian framework, which gives geodesic equations of motion identical to those of the parent metric. Consequently, the derived velocity-dependent relativistic scalar potential, which is a relativistic generalization of Newtonian gravitational potential, exactly reproduces the relativistic gravitational features corresponding to any static spherical symmetric spacetime geometry in its entirety, including all the experimentally tested gravitational effects in the weak field up to the present. This relativistic analogous potential is expected to be quite useful in studying wide range of astrophysical phenomena, especially in strong field gravity.

This work reports the first discovery of (negative) spectral lags in the X-ray emission below 10 keV from the gamma ray binary $LSI +61^{o} $ 303 during large flaring episodes using the Rossi X-ray Timing Explorer (RXTE) observations. It is found from the RXTE data that during the flares, low energy (3-5 KeV) variations lead the higher energy (8-10 keV) variations by few tens of seconds whereas no significant time lag is observed during the non-flaring states. The observed spectral lag features for flaring events suggest that inverse Compton scattering may be operative, at least in some part of the system. Another possibility is that the sites of particle acceleration may be different for flaring and non-flaring events such as in the microquasar model the flaring radiation may come from hot spots sitting above the black hole while steady state emissions are due to the jets.

The effects of several dark energy models on gravitational time delay of particles with non-zero mass are investigated and analytical expressions for the same are obtained at the first order accuracy. Also the expression for gravitational time delay under the influence of conformal gravity potential that well describes the flat rotation curve of spiral galaxies is derived. The findings suggest that i) the conformal gravity description of dark matter reduces the net time delay in contrast to the effect of normal dark matter and therefore in principle the models can be discriminated using gravitational time delay observations and ii)the effect of dark energy/flat rotation curve may be revealed from high precision measurements of gravitational time delay of particles involving megaparsec and beyond distance-scale.

In this work we examine with the help of Monte Carlo simulation whether a consistent primary energy spectrum of cosmic rays emerges from both the experimentally observed total charged particles and muon size spectra of cosmic ray extensive air showers considering primary composition may or may not change beyond the knee of the energy spectrum. It is found that EAS-TOP observations consistently infer a knee in the primary energy spectrum provided the primary is pure unchanging iron whereas no consistent primary spectrum emerges from simultaneous use of the KASCADE observed total charged particle and muon spectra. However, it is also found that when primary composition changes across the knee the estimation of spectral index of total charged particle spectrum is quite tricky, depends on the choice of selection of points near the knee in the size spectrum.

We formulate a generic Newtonian like analogous potential for static spherically symmetric general relativistic (GR) spacetime, and subsequently derived proper Newtonian like analogous potential corresponding to Janis-Newman-Winicour (JNW) and Reissner-Nordstr\"{o}m (RN) spacetimes, both exhibiting naked singularities. The derived potentials found to reproduce the entire GR features including the orbital dynamics of the test particle motion and the orbital trajectories, with precise accuracy. The nature of the particle orbital dynamics including their trajectory profiles in JNW and RN geometries show altogether different behavior with distinctive traits as compared to the nature of particle dynamics in Schwarzschild geometry. Exploiting the Newtonian like analogous potentials, we found that the radiative efficiency of a geometrically thin and optically thick Keplerian accretion disk around naked singularities corresponding to both JNW and RN geometries, in general, is always higher than that for Schwarzschild geometry. The derived potentials would thus be useful to study astrophysical processes, especially to investigate more complex accretion phenomena in AGNs or in XRBs in the presence of naked singularities and thereby exploring any noticeable differences in their observational features from those in the presence of BHs to ascertain outstanding debatable issues relating to gravity - whether the end state of gravitational collapse in our physical Universe renders black hole (BH) or naked singularity.

The effect of dark matter/energy on gravitational time advancement (negative effective time delay) has been investigated considering few dark energy/matter models including cosmological constant. It is found that dark energy gives only (positive) gravitational time delay irrespective of the position of the observer whereas pure Schwarzschild geometry leads to gravitational time advancement when the observer is situated at relatively stronger gravitational field point in the light trajectory. Consequently, there will be no time advancement effect at all at radial distances where gravitational field due to dark energy is stronger than the gravitational field of Schwarzschild geometry.

Primary energy spectrum of cosmic rays exhibits a knee at about $3$ PeV where a change in the spectral index occurs. Despite many efforts the origin of such a feature of the spectrum is not satisfactorily solved yet. Here it is proposed that the steepening of the spectrum beyond the knee may be a consequence of mass distribution of progenitor of cosmic ray source. The proposed speculative model can account all the major observed features about cosmic rays without invoking any fine tuning to match flux or spectra at any energy point. The prediction of the proposed model regarding primary composition scenario beyond the knee is quite different from most of the prevailing models of the knee and thereby can be discriminated from precise experimental measurement of the primary composition.

Distinguishing $\gamma$-ray and hadron initiated extensive air showers (EAS) based on lateral distribution of electrons has been studied by detailed Monte Carlo (MC) simulations. The possibility of using the local age parameter (LAP) of EAS for the gamma-hadron separation has been explored. It is found that separating $\gamma$-ray and hadron induced EAS on the basis of LAP can be useful for surface detector experiments those have no reliable muon measurement facilities.

Based on the conserved Hamiltonian for a test particle, we have formulated a Newtonian analogue of Kerr spacetime in the `low energy limit of the test particle motion' that, in principle, can be comprehensively used to describe general relativistic (GR) features of Kerr spacetime, however, with less accuracy for high spin. The derived potential, which has an explicit velocity dependence, contains the entire relativistic features of corresponding spacetime including the frame dragging effect, unlike other prevailing pseudo-Newtonian potentials (PNPs) for the Kerr metric where such an effect is either totally missing or introduced in a ad hoc manner. The particle dynamics with this potential precisely reproduce the GR results within a maximum ~ 10 % deviation in energy for a particle orbiting circularly in the vicinity of a rapidly corotating black hole. GR epicyclic frequencies are also well reproduced with the potential, though with a relatively higher percentage of deviation. For counterrotating cases, the obtained potential replicate the GR results with precise accuracy. The Kerr-Newtonian potential also approximates the radius of marginally stable and marginally bound circular orbits with reasonable accuracy for a < 0.7. Importantly, the derived potential can imitate the experimentally tested GR effects like perihelion advancement and bending of light with reasonable accuracy. The formulated Kerr-Newtonian potential thus can be useful to study complex accreting plasma dynamics and its implications around rotating BHs in the Newtonian framework, avoiding GR gas dynamical equations.

The late time accelerated expansion of the Universe demands that even in local galactic-scales it is desirable to study astrophysical phenomena, particularly relativistic accretion related phenomena in massive galaxies or in galaxy mergers and the dynamics of the kiloparsecs-scale structure and beyond, in the local-galaxies in Schwarzschild-de Sitter (SDS) background, rather than in Schwarzschild or Newtonian paradigm. Owing to the complex and nonlinear character of the underlying magnetohydrodynamical equations in general relativistic (GR) regime, it is quite useful to have an Newtonian analogous potential containing all the important GR features that allows to treat the problem in Newtonian framework for study of accretion and its related processes. From the principle of conserved Hamiltonian of the test particle motion, here, a three dimensional Newtonian analogous potential has been obtained in spherical geometry corresponding to SDS/Schwarzschild anti-de Sitter (SADS) spacetime, that reproduces almost all of the GR features in its entirety with remarkable accuracy. The derived potential contains an explicit velocity dependent term of the test particle that renders an approximate relativistic modification of Newtonian like potential. The complete orbital dynamics around SDS geometry and the epicyclic frequency corresponding to SDS metric have been extensively studied in the Newtonian framework using the derived potential. Applying the derived analogous potential it is found that the current accepted value of $\Lambda \sim 10^{-56} \rm cm^{-2}$ moderately influences both sonic radius as well as Bondi accretion rate, especially for spherical accretion with smaller values of adiabatic constant and temperature, which might have interesting consequences on the stability of accretion disk in AGNs/radio galaxies.

Conformal gravity theory can explain observed flat rotation curves of galaxies without invoking hypothetical dark matter. Within this theory, we obtain a generic formula for the sizes of galaxies exploiting the stability criterion of circular orbits. It is found that different galaxies have different finite sizes uniquely caused by the assumed quadratic potential of cosmological origin. Observations on where circular orbits might actually terminate could thus be very instructive in relation to the galactic sizes predicted here.

Analyzing simulated EAS events generated with the CORSIKA code, the characteristics of lateral distribution of electrons in EAS around the knee energy region of the primary energy spectrum have been studied and compared with experimental observations. The differences between the EGS4 and the NKG output of CORSIKA in respect to electron radial density distribution have been investigated. The relation between lateral and longitudinal age parameters has been studied after introducing the notion of the local age parameter that reflects the profile of the lateral distribution of electrons in EAS. The present analysis motivates the inclusion of the lateral shower age in a multiparameter study of EAS to extract information on hadronic interactions and primary composition.

After discussing some subtle issues concerning the computation of deflection angle, a general but simple expression of bending angle of light rays in weak deflection limit has been presented for a general static and spherically symmetric space-time. In this context the importance of proper choice of polar axis has been highlighted. Applying the prescribed method the explicit expression of bending angle up to an accuracy of second order in mass is obtained for the Schwarzschild geometry without restricting the source and the observer to be at infinite distance away from the lens.

We study low energy hadronic interaction models based on BESS observed cosmic ray proton and antiproton spectra at medium high altitude. Among the three popular low energy interaction models, we find that FLUKA reproduces results of BESS observations on secondary proton spectrum reasonably well over the whole observed energy range, the model UrQMD works well at relatively higher energies whereas spectrum obtained with GHEISHA differs significantly from the measured spectrum. Simulated antiproton spectrum with FLUKA, however, exhibits significant deviations from the BESS observation wheras UrQMD and GHEISHA reproduce the BESS observations within the experimental error.

Recent studies suggest that pulsars could be strong sources of TeV muon neutrinos provided positive ions are accelerated by pulsar polar caps to PeV energies. In such a situation muon neutrinos are produced through the delta resonance in interactions of pulsar accelerated ions with its thermal radiation field. High energy gamma rays also should be produced simultaneously in pulsar environment as both charged and neutral pions are generated in the interactions of energetic hadrons with the ambient photon fields. Here we estimate TeV gamma ray flux at Earth from few nearby young pulsars. When compared with the observations we find that proper consideration of the effect of polar cap geometry in flux calculation is important. Incorporating such an effect we obtain the (revised) event rates at Earth due to few potential nearby pulsars. The results suggest that pulsars are unlikely to be detected by the upcoming neutrino telescopes. We also estimate TeV gamma ray and neutrino fluxes from pulsar nebulae for the adopted model of particle acceleration.

Contribution of nearby pulsars to the cosmic rays observed at Earth has been studied. It is found that the experimental bound on amplitude of cosmic ray anisotropy may produce significant constraint on the efficiency of converting pulsar rotational energy to emitted particles kinetic energy. Cosmic ray fluxes from two well known nearby gamma ray pulsars, namely the Vela and Geminga pulsars, are estimated. The analysis suggests that observed bound on cosmic ray anisotropy restricts the contributions of the Vela and the Geminga pulsars to at most 1 % of the observed cosmic rays below the knee.

We study gravitational lensing by the Gibbons-Maeda-Garfinkle-Horowitz-Strominger (GMGHS) charged black hole of heterotic string theory and obtain the angular position and magnification of the relativistic images. Modeling the supermassive central object of the galaxy as a GMGHS black hole, we estimate the numerical values of different strong-lensing parameters. We find that there is no significant string effect present in the lensing observables in the strong-gravity scenario.