With the advent of the Large Synoptic Survey Telescope (LSST), time-domain astronomy will be faced with an unprecedented volume and rate of data. Real-time processing of variables and transients detected by such large-scale surveys is critical to identifying the more unusual events and allocating scarce follow-up resources efficiently. We develop an algorithm to identify these novel events within a given population of variable sources. We determine the distributions of magnitude changes (dm) over time intervals (dt) for a given passband f, pf(dm|dt), and use these distributions to compute the likelihood of a test source being consistent with the population, or an outlier. We demonstrate our algorithm by applying it to the DECam multi-band time-series data of more than 2000 variable stars identified by Saha et al. (2019) in the Galactic Bulge that are largely dominated by long-period variables and pulsating stars. Our algorithm discovers 18 outlier sources in the sample, including a microlensing event, a dwarf nova, and two chromospherically active RS CVn stars, as well as sources in the Blue Horizontal Branch region of the color-magnitude diagram without any known counterparts. We compare the performance of our algorithm for novelty detection with multivariate KDE and Isolation Forest on the simulated PLAsTiCC dataset. We find that our algorithm yields comparable results despite its simplicity. Our method provides an efficient way for flagging the most unusual events in a real-time alert-broker system.

The JUNO neutrino detector system is simulated using Monte-Carlo and analytical methods. A large number of proton decay events are also simulated. Preliminary results from this endeavor are presented in the present article.

A detailed numerical procedure has been developed to simulate the mechanical configurations and optical properties of Imaging Atmospheric Cerenkov Telescope systems. To test these procedures a few existing ACT arrays are simulated. First results from these simulations are presented.

IUCAA Digital Sampling Array Controller (IDSAC) is a flexible and generic yet powerful CCD controller which can handle a wide range of scientific detectors. Based on an easily scalable modular backplane architecture consisting of Single Board Controllers (SBC), IDSAC can control large detector arrays and mosaics. Each of the SBCs offers the full functionality required to control a CCD independently. The SBCs can be cold swapped without the need to reconfigure them. Each SBC can handle data from up to four video channels with or without dummy outputs at speeds up to 500 kilo Pixels Per Second (kPPS) Per Channel with a resolution of 16 bits. Communication with Linux based host computer is through a USB3.0 interface, with the option of using copper or optical fibers. A Field Programmable Gate Array (FPGA) is used as the master controller in each SBC which allows great flexibility in optimizing performance by adjusting gain, timing signals, bias levels, etc. using user-editable configuration files without altering the circuit topology. Elimination of thermal kTC noise is achieved via Digital Correlated Double Sampling (DCDS). We present the results of noise performance characterization of IDSAC through simulation, theoretical modeling, and actual measurements. The contribution of different types of noise sources is modeled using a tool to predict noise of a generic DCDS signal chain analytically. The analytical model predicts the net input referenced noise of the signal chain to be 5 electrons for 200k pixels per second per channel readout rate with 3 samples per pixel. Using a cryogenic test set up in the lab, the noise is measured to be 5.4 e (24.3 \muV), for the same readout configuration.

The Hyper-Kamiokande Water-Cherenkov Neutrino Detectors are simulated using Monte-Carlo and analytical methods. A few simple events are also simulated and these preliminary results are presented.

The present work deals with the late time evolution of the linear density contrast in the dark energy models reconstructed from the jerk parameter. It is found that the non-interacting models are favoured compared to the models where an interaction is allowed in the dark sector.

Core-collapse supernovae (CCSNe) are the extremely energetic deaths of massive stars. They play a vital role in the synthesis and dissemination of many heavy elements in the universe. In the past, CCSN nucleosynthesis calculations have relied on artificial explosion methods that do not adequately capture the physics of the innermost layers of the star. The PUSH method, calibrated against SN1987A, utilizes the energy of heavy-flavor neutrinos emitted by the proto-neutron star (PNS) to trigger parametrized explosions. This makes it possible to follow the consistent evolution of the PNS and to ensure a more accurate treatment of the electron fraction of the ejecta. Here, we present the Iron group nucleosynthesis results for core-collapse supernovae, exploded with PUSH, for two different progenitor series. Comparisons of the calculated yields to observational metal-poor star data are also presented. Nucleosynthesis yields will be calculated for all elements and over a wide range of progenitor masses. These yields can be immensely useful for models of galactic chemical evolution.

We report on a method, PUSH, for triggering core-collapse supernova (CCSN) explosions of massive stars in spherical symmetry. This method provides a framework to study many important aspects of core collapse supernovae: the effects of the shock passage through the star, explosive supernova nucleosynthesis and the progenitor-remnant connection. Here we give an overview of the method, compare the results to multi-dimensional simulations and investigate the effects of the progenitor and the equation of state on black hole formation.

The Cadmium Zinc Telluride Imager (CZTI) is a high energy, wide-field imaging instrument on AstroSat. CZT's namesake Cadmium Zinc Telluride detectors cover an energy range from 20 keV to > 200 keV, with 11% energy resolution at 60 keV. The coded aperture mask attains an angular resolution of 17' over a 4.6 deg x 4.6 deg (FWHM) field of view. CZTI functions as an open detector above 100 keV, continuously sensitive to GRBs and other transients in about 30% of the sky. The pixellated detectors are sensitive to polarisation above ~100 keV, with exciting possibilities for polarisation studies of transients and bright persistent sources. In this paper, we provide details of the complete CZTI instrument, detectors, coded aperture mask, mechanical and electronic configuration, as well as data and products.

We investigate the role of the spin angular momentum of astrophysical black holes in controlling the special relativistic chaotic dynamics of test particles moving under the influence of a post-Newtonian pseudo-Kerr black hole potential, along with a perturbative potential created by a asymmetrically placed (dipolar) halo. Proposing a Lyapunov-like exponent to be the effective measure of the degree of chaos observed in the system under consideration, it has been found that black hole spin anti-correlates with the degree of chaos for the aforementioned dynamics. Our findings have been explained applying the general principles of dynamical systems analysis.

We present the results of a search for gravitational-wave bursts associated with 137 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments during the fifth LIGO science run and first Virgo science run. The data used in this analysis were collected from 2005 November 4 to 2007 October 1, and most of the GRB triggers were from the Swift satellite. The search uses a coherent network analysis method that takes into account the different locations and orientations of the interferometers at the three LIGO-Virgo sites. We find no evidence for gravitational-wave burst signals associated with this sample of GRBs. Using simulated short-duration (<1 s) waveforms, we set upper limits on the amplitude of gravitational waves associated with each GRB. We also place lower bounds on the distance to each GRB under the assumption of a fixed energy emission in gravitational waves, with typical limits of D ~ 15 Mpc (E_GW^iso / 0.01 M_o c^2)^1/2 for emission at frequencies around 150 Hz, where the LIGO-Virgo detector network has best sensitivity. We present astrophysical interpretations and implications of these results, and prospects for corresponding searches during future LIGO-Virgo runs.

We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in the search. We report no signal detection from any of the targets and therefore interpret our results as upper limits on the gravitational wave signal strength. The most interesting limits are those for young pulsars. We present updated limits on gravitational radiation from the Crab pulsar, where the measured limit is now a factor of seven below the spin-down limit. This limits the power radiated via gravitational waves to be less than ~2% of the available spin-down power. For the X-ray pulsar J0537-6910 we reach the spin-down limit under the assumption that any gravitational wave signal from it stays phase locked to the X-ray pulses over timing glitches, and for pulsars J1913+1011 and J1952+3252 we are only a factor of a few above the spin-down limit. Of the recycled millisecond pulsars several of the measured upper limits are only about an order of magnitude above their spin-down limits. For these our best (lowest) upper limit on gravitational wave amplitude is 2.3x10^-26 for J1603-7202 and our best (lowest) limit on the inferred pulsar ellipticity is 7.0x10^-8 for J2124-3358.

The modulation of high-energy transients' (or steadily emitting sources') light curves due to the imperfect alignment of the detector's view axis with the spin axis in a spin-stabilized satellite is derived. It is shown how the orientation of the detector's view axis with respect to the satellite's spin axis may be estimated using observed light curves. The effects of statistical fluctuations are considered. Conversely, it is shown how the attitude of a spin-axis stabilized satellite as well as the unknown position of a celestial source of high-energy photons may be determined using a detector whose view-axis is intentionally kept inclined and is known accurately beforehand. The case of three-axes stabilized satellites is also discussed.

The value of the cosmological constant is one of the major puzzles of modern cosmology: it is tiny but nonzero. Sorkin predicted, from the Causet approach to quantum gravity, that the cosmological constant has quantum fluctuations. The predicted order of magnitude of the fluctuations agrees with the subsequently observed value of the cosmological constant. We had earlier developed an analogy between the cosmological constant of the Universe and the surface tension of fluid membranes. Here we demonstrate by computer simulations that the surface tension of a fluid membrane has statistical fluctuations stemming from its discrete molecular structure. Our analogy enables us to view these numerical experiments as probing a small and fluctuating cosmological constant. Deriving insights from our analogy, we show that a fluctuating cosmological constant is a GENERIC feature of quantum gravity models and is far more general than the specific context in which it was originally proposed. We pursue and refine the idea of a fluctuating cosmological constant and work towards making further testable predictions.

Recently it was shown that the inclusion of higher signal harmonics in the inspiral signals of binary supermassive black holes (SMBH) leads to dramatic improvements in parameter estimation with the Laser Interferometer Space Antenna (LISA). In particular, the angular resolution becomes good enough to identify the host galaxy or galaxy cluster, in which case the redshift can be determined by electromagnetic means. The gravitational wave signal also provides the luminosity distance with high accuracy, and the relationship between this and the redshift depends sensitively on the cosmological parameters, such as the equation-of-state parameter $w=p_{\rm DE}/\rho_{\rm DE}$ of dark energy. With a single binary SMBH event at $z < 1$ having appropriate masses and orientation, one would be able to constrain $w$ to within a few percent. We show that, if the measured sky location is folded into the error analysis, the uncertainty on $w$ goes down by an additional factor of 2-3, leaving weak lensing as the only limiting factor in using LISA as a dark energy probe.

The LISA Parameter Estimation (LISAPE) Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models, and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large, and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.

A detailed Monte-Carlo code has been developed from basic principles that simulates almost all of the basic photon and charged particle interactions. The code is used to derive the response functions of a high energy photon detector to incident beams of photons of various energies. The detector response matrices (DRMs) are calculated using this code. Deconvolution of an artificially generated spectrum is presented.

We present direct upper limits on gravitational wave emission from the Crab pulsar using data from the first nine months of the fifth science run of the Laser Interferometer Gravitational-wave Observatory (LIGO). These limits are based on two searches. In the first we assume that the gravitational wave emission follows the observed radio timing, giving an upper limit on gravitational wave emission that beats indirect limits inferred from the spin-down and braking index of the pulsar and the energetics of the nebula. In the second we allow for a small mismatch between the gravitational and radio signal frequencies and interpret our results in the context of two possible gravitational wave emission mechanisms.

It is generally believed that the angular resolution of the Laser Interferometer Space Antenna (LISA) for binary supermassive black holes (SMBH) will not be good enough to identify the host galaxy or galaxy cluster. This conclusion, based on using only the dominant harmonic of the binary SMBH signal, changes substantially when higher signal harmonics are included in assessing the parameter estimation problem. We show that in a subset of the source parameter space the angular resolution increases by more than a factor of 10, thereby making it possible for LISA to identify the host galaxy/galaxy cluster. Thus, LISA's observation of certain binary SMBH coalescence events could constrain the dark energy equation of state to within a few percent, comparable to the level expected from other dark energy missions.

Current expectations on the signal to noise ratios and masses of supermassive black holes which the Laser Interferometer Space Antenna (LISA) can observe are based on using in matched filtering only the dominant harmonic of the inspiral waveform at twice the orbital frequency. Other harmonics will affect the signal-to-noise ratio of systems currently believed to be observable by LISA. More significantly, inclusion of other harmonics in our matched filters would mean that more massive systems that were previously thought to be {\it not} visible in LISA should be detectable with reasonable SNRs. Our estimates show that we should be able to significantly increase the mass reach of LISA and observe the more commonly occurring supermassive black holes of masses $\sim 10^8M_\odot.$ More specifically, with the inclusion of all known harmonics LISA will be able to observe even supermassive black hole coalescences with total mass $\sim 10^8 M_\odot (10^9M_\odot)$ (and mass-ratio 0.1) for a low frequency cut-off of $10^{-4}{\rm Hz}$ $(10^{-5}{\rm Hz})$ with an SNR up to $\sim 60$ $(\sim 30)$ at a distance of 3 Gpc. This is important from the astrophysical viewpoint since observational evidence for the existence of black holes in this mass range is quite strong and binaries containing such supermassive black holes will be inaccessible to LISA if one uses as detection templates only the dominant harmonic.

The astronomically observed value of the cosmological constant is small but non-zero. This raises two questions together known as the cosmological constant problem a) why is lambda so nearly zero? b) why is lambda not EXACTLY zero? Sorkin has proposed that b) can be naturally explained as a one by square root N fluctuation by invoking discreteness of spacetime at the Planck scale due to quantum gravity. In this paper we shed light on these questions by developing an analogy between the cosmological constant and the surface tension of membranes. The ``cosmological constant problem'' has a natural analogue in the membrane context: the vanishingly small surface tension of fluid membranes provides an example where question a) above arises and is answered. We go on to find a direct analogue of Sorkin's proposal for answering question b) in the membrane context, where the discreteness of spacetime translates into the molecular structure of matter. We propose analogue experiments to probe a small and fluctuating surface tension in fluid membranes. A counterpart of dimensional reduction a la Kaluza-Klein and large extra dimensions also appears in the physics of fluid membranes.

We present Interplanetary Network (IPN) detection and localization information for 211 gamma-ray bursts (GRBs) observed as untriggered events by the Burst and Transient Source Experiment (BATSE), and published in catalogs by Kommers et al. (2001) and Stern et al. (2001). IPN confirmations have been obtained by analyzing the data from 11 experiments. For any given burst observed by BATSE and one other distant spacecraft, arrival time analysis (or ``triangulation'') results in an annulus of possible arrival directions whose half-width varies between 14 arcseconds and 5.6 degrees, depending on the intensity, time history, and arrival direction of the burst, as well as the distance between the spacecraft. This annulus generally intersects the BATSE error circle, resulting in a reduction of the area of up to a factor of ~650. When three widely separated spacecraft observed a burst, the result is an error box whose area is as much as 30000 times smaller than that of the BATSE error circle. Because the IPN instruments are considerably less sensitive than BATSE, they generally did not detect the weakest untriggered bursts, but did detect the more intense ones which failed to trigger BATSE when the trigger was disabled. In a few cases, we have been able to identify the probable origin of bursts as soft gamma repeaters. The vast majority of the IPN-detected events, however, are GRBs, and the confirmation of them validates many of the procedures utilized to detect BATSE untriggered bursts.

The existence of compact gravitational lenses, with masses around 0.5 (M_{\odot}), has been reported in the halo of the Milky Way. The nature of these dark lenses is as yet obscure, particularly because these objects have masses well above the threshold for nuclear fusion. In this work, we show that they find a natural explanation as being the evolutionary product of the metastable false vacuum domains (the so-called strange quark nuggets) formed in a first order cosmic quark-hadron transition.

Using the concept of open systems where the classical geometry is treated as the system and the quantum matter field as the environment, we derive a fluctuation-dissipation theorem for semiclassical cosmology. This theorem which exists under very general conditions for dissipations in the dynamics of the system, and the noise and fluctuations in the environment, can be traced to the formal mathematical relation between the dissipation and noise kernels of the influence functional depicting the open system, and is ultimately a consequence of the unitarity of the closed system. In particular, for semiclassical gravity, it embodies the backreaction effect of matter fields on the dynamics of spacetime. The backreaction equation derivable from the influence action is in the form of a Einstein-Langevin equation. It contains a dissipative term in the equation of motion for the dynamics of spacetime and a noise term related to the fluctuations of particle creation in the matter field. Using the well-studied model of a quantum scalar field in a Bianchi Type-I universe we illustrate how this Langevin equation and the noise term are derived and show how the creation of particles and the dissipation of anisotropy during the expansion of the universe can be understood as a manifestation of this fluctuation-dissipation relation.

A perturbed Reissner-Nordstr\"om-de Sitter solution is used to emphasize the nature of the singularity along the Cauchy horizon of a charged spherically symmetric black hole. For these solutions, conditions may prevail under which the mass function is bounded and yet the curvature scalar $R_{\alpha\beta\gamma\delta} R^{\alpha\beta\gamma\delta}$ diverges.