Bursts of gamma ray showers have been observed in coincidence with downward propagating negative leaders in lightning flashes by the Telescope Array Surface Detector (TASD). The TASD is a 700~square kilometer cosmic ray observatory located in southwestern Utah, U.S.A. In data collected between 2014 and 2016, correlated observations showing the structure and temporal development of three shower-producing flashes were obtained with a 3D lightning mapping array, and electric field change measurements were obtained for an additional seven flashes, in both cases co-located with the TASD. National Lightning Detection Network (NLDN) information was also used throughout. The showers arrived in a sequence of 2--5 short-duration ($\le$10~$\mu$s) bursts over time intervals of several hundred microseconds, and originated at an altitude of $\simeq$3--5 kilometers above ground level during the first 1--2 ms of downward negative leader breakdown at the beginning of cloud-to-ground lightning flashes. The shower footprints, associated waveforms and the effect of atmospheric propagation indicate that the showers consist primarily of downward-beamed gamma radiation. This has been supported by GEANT simulation studies, which indicate primary source fluxes of $\simeq$$10^{12}$--$10^{14}$ photons for $16^{\circ}$ half-angle beams. We conclude that the showers are terrestrial gamma-ray flashes (TGFs), similar to those observed by satellites, but that the ground-based observations are more representative of the temporal source activity and are also more sensitive than satellite observations, which detect only the most powerful TGFs.

The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere. Cosmic ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors. We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data. Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.

IGR J14257-6117 is an unclassified source in the hard X-ray catalogues. Optical follow-ups suggest it could be a Cataclysmic Variable of the magnetic type. We present the first high S/N X-ray observation performed by \XMM\ at 0.3--10 keV, complemented with 10--80 keV coverage by \Swift/BAT, aimed at revealing the source nature. We detected for the first time a fast periodic variability at 509.5\,s and a longer periodic variability at 4.05\,h, ascribed to the white dwarf (WD) spin and binary orbital periods, respectively. These unambiguously identify IGR J14257-6117 as a magnetic CV of the Intermediate Polar (IP) type. The energy resolved light curves at both periods reveal amplitudes decreasing with increasing energy, with the orbital modulation reaching $\sim100\%$ in the softest band. The energy spectrum shows optically thin thermal emission with an excess at the iron complex, absorbed by two dense media (${\rm N_{H}\sim10^{22-23}\,cm^{-2}}$), partially covering the X-ray source. These are likely localised in the magnetically confined accretion flow above the WD surface and at the disc rim, producing the energy dependent spin and orbital variabilities, respectively. IGR J14257-6117, joins the group of strongest orbitally modulated IPs now counting four systems. Drawing similarities with low-mass X-ray binaries displaying orbital dips, these IPs should be seen at large orbital inclinations allowing azimuthally extended absorbing material fixed in the binary frame to intercept the line of sight. For IGR J14257-6117, we estimate ($50^o\,\lesssim\,i\,\lesssim\,70^o$). Whether also the mass accretion rate plays a role in the large orbital modulations in IPs cannot be established with the present data.

One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured number of muons with the MC prediction. In this paper, we present the results of studies of muon densities in UHE extensive air showers obtained by analyzing the signal of surface detector stations which should have high $\it{muon \, purity}$. The muon purity of a station will depend on both the inclination of the shower and the relative position of the station. In 7 years' data from the Telescope Array experiment, we find that the number of particles observed for signals with an expected muon purity of $\sim$65% at a lateral distance of 2000 m from the shower core is $1.72 \pm 0.10{\rm (stat.)} \pm 0.37 {\rm (syst.)}$ times larger than the MC prediction value using the QGSJET II-03 model for proton-induced showers. A similar effect is also seen in comparisons with other hadronic models such as QGSJET II-04, which shows a $1.67 \pm 0.10 \pm 0.36$ excess. We also studied the dependence of these excesses on lateral distances and found a slower decrease of the lateral distribution of muons in the data as compared to the MC, causing larger discrepancy at larger lateral distances.

We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources.

An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson likelihood ratio. The center of maximum significance is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies $10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$ eV. The post-trial probability of this energy anisotropy, appearing by chance anywhere on an isotropic sky, is found by Monte Carlo simulation to be $9$$\times$$10^{-5}$ ($3.74$$\sigma_{global}$).

The X-ray spectra of intermediate polars can be modelled to give a direct measurement of white dwarf mass. Here we fit accretion column models to NuSTAR spectra of three intermediate polars; V709 Cas, NY Lup and V1223 Sgr in order to determine their masses. From fits to 3-78 keV spectra, we find masses of $M_{\rm WD}=0.88^{+0.05}_{-0.04}M_{\odot}$, $1.16^{+0.04}_{-0.02}M_{\odot}$ and $0.75\pm0.02M_{\odot}$ for V709 Cas, NY Lup and V1223 Sgr, respectively. Our measurements are generally in agreement with those determined by previous surveys of intermediate polars, but with typically a factor $\sim2$ smaller uncertainties. This work paves the way for an approved NuSTAR Legacy Survey of white dwarf masses in intermediate polars.

The energy of the ultrahigh energy spectral cutoff was measured, integrating over the northern hemisphere sky, by the Telescope Array (TA) collaboration, to be $10^{19.78 \pm 0.06}$ eV, in agreement with the High Resolution Fly's Eye (HiRes) experiment, whereas the Pierre Auger experiment, integrating over the southern hemisphere sky, measured the cutoff to be at 10$^{19.62 \pm 0.02}$ eV. An 11% energy scale difference between the TA and Auger does not account for this difference. However, in comparing the spectra of the Telescope Array and Pierre Auger experiments in the band of declination common to both experiments ($-15.7^{\circ} < \delta < 24.8^{\circ}$) we have found agreement in the energy of the spectral cutoff. While the Auger result is essentially unchanged, the TA cutoff energy has changed to $10^{19.59 \pm 0.06}$ eV. In this paper we argue that this is an astrophysical effect.

Compared to mass transfer in cataclysmic variables, the nature of accretion in symbiotic binaries in which red giants transfer material to white dwarfs (WDs) has been difficult to uncover. The accretion flows in a symbiotic binary are most clearly observable, however, when there is no quasi-steady shell burning on the WD to hide them. RT Cru is the prototype of such non-burning symbiotics, with its hard (\delta-type) X-ray emission providing a view of its innermost accretion structures. In the past 20 years, RT Cru has experienced two, similar optical brightening events, separated by \sim 4000 days and with amplitudes of \Delta V \sim 1.5 mag. After Swift became operative, the BAT detector revealed a hard X-ray brightening event almost in coincidence with the second optical peak. Spectral and timing analyses of multi-wavelength observations that we describe here, from NuSTAR, Suzaku, Swift XRT+BAT+UVOT and optical photometry and spectroscopy, indicate that accretion proceeds through a disk that reaches down to the white dwarf surface. The scenario in which a massive, magnetic WD accretes from a magnetically truncated accretion disk is not supported. For example, none of our data show the minute-time-scale periodic modulations (with tight upper limits from X-ray data) expected from a spinning, magnetic WD. Moreover, the similarity of the UV and X-ray fluxes, as well as the approximate constancy of the hardness ratio within the BAT band, indicate that the boundary layer of the accretion disk remained optically thin to its own radiation throughout the brightening event, during which the rate of accretion onto the WD increased to 6.7$\times$ 10$^{-9}$ M$_{\odot}$ yr$^{-1}$ (d/2 kpc)$^2$. (Abridged abstract version, see full version in the pdf)

In magnetically accreting white dwarfs, the height above the white dwarf surface where the standing shock is formed is intimately related with the accretion rate and the white dwarf mass. However, it is difficult to measure. We obtained new data with NuSTAR and Swift that, together with archival Chandra data, allow us to constrain the height of the shock in the intermediate polar EX Hya. We conclude that the shock has to form at least at a distance of about one white dwarf radius from the surface in order to explain the weak Fe K{\alpha} 6.4 keV line, the absence of a reflection hump in the high-energy continuum, and the energy dependence of the white dwarf spin pulsed fraction. Additionally, the NuSTAR data allowed us to measure the true, uncontaminated hard X-ray (12-40 keV) flux, whose measurement was contaminated by the nearby galaxy cluster Abell 3528 in non-imaging X-ray instruments.

WD 1145+017 was observed from 2016 November through 2017 June for the purpose of further characterizing the transit behavior of the dusty debris clouds orbiting this white dwarf. The optical observations were carried out with a small ground-based telescope run by an amateur astronomer, and covered 53 different nights over the 8-month interval. We have found that the optical activity has increased to the highest level observed since its discovery with Kepler K2, with approximately 17% of the optical flux extinguished per orbit. The source exhibits some transits with depths of up to 55% and durations as long as two hours. The dominant period of the orbiting dust clouds during 2016-2017 is 4.49126 hours. We present 'waterfall' images for the entire 2016-2017 and 2015-2016 observing seasons. In addition, the white dwarf was observed with the Chandra X-ray Observatory for 10-ksec on each of four different occasions, separated by about a month each. The upper limit on the average X-ray flux from WD 1145+017 is ~5 x 10^{-15} ergs/cm^2/s (unabsorbed over the range 0.1-100 keV), which translates to an upper limit on the X-ray luminosity, Lx, of ~2 x 10^{28} ergs/s. If Lx = G Mwd Mdot/Rwd, where Mwd and Rwd are the mass and radius of the white dwarf, and Mdot is the accretion rate, then Mdot < 2 x 10^{11} g/s. This is just consistent with the value of Mdot that is inferred from the level of dust activity.

The Telescope Array (TA) experiment is located in the western desert of Utah, USA, and observes ultra high energy cosmic rays (UHECRs) in the Northern hemisphere. At the highest energies, $E>10$~EeV, the shape of cosmic ray energy spectrum may carry an imprint of the source density distribution along the line of sight different in different directions of the sky. In this study, we search for such directional variations in the shape of the energy spectrum using events observed with the Telescope Array's surface detector. We divide the TA field of view into two nearly equal-exposure regions: the "on-source" region which we define as $\pm 30^\circ$ of the supergalactic plane containing mostly nearby structures, and the complementary "off-source" region where the sources are further away on average. We compare the UHECR spectra in these regions by fitting them to the broken power law and comparing the resulting parameters. We find that the off-source spectrum has an earlier break at highest energies. The chance probability to obtain such or larger difference in statistically equivalent distributions is estimated as $6.2\pm1.1\times10^{-4}$ ($3.2\sigma$) by a Monte-Carlo simulation. The observed difference in spectra is in a reasonable quantitative agreement with a simplified model that assumes that the UHECR sources trace the galaxy distribution from the 2MRS catalogue, primary particles are protons and the magnetic deflections can be neglected.

We present timing and spectral analysis of a sample of seven hard X-ray selected Cataclysmic Variable candidates based on simultaneous X-ray and optical observations collected with XMM-Newton , complemented with Swift/BAT and INTEGRAL/IBIS hard X-ray data and ground-based optical photometry. For six sources, X-ray pulsations are detected for the first time in the range $\rm \sim296-6098\,s$, identifying them as members of the magnetic class. Swift J0927.7-6945, Swift J0958.0-4208, Swift J1701.3-4304, Swift J2113.5+5422, and possibly PBC J0801.2-4625, are Intermediate Polars (IPs), while Swift J0706.8+0325 is a short (1.7 h) orbital period Polar, the 11$^{\rm th}$ hard X-ray selected identified so far. X-ray orbital modulation is also observed in Swift J0927.7-6945 (5.2 h) and Swift J2113.5+5422 (4.1 h). Swift J1701.3-4304 is discovered as the longest orbital period (12.8 h) deep eclipsing IP. The spectra of the magnetic systems reveal optically thin multi-temperature emission between 0.2 and 60 keV. Energy dependent spin pulses and the orbital modulation in Swift J0927.7-6945 and Swift J2113.5+5422 are due to intervening local high density absorbing material ($\rm N_H\sim10^{22-23}\,cm^{-2}$). In Swift J0958.0-4208 and Swift J1701.3-4304, a soft X-ray blackbody (kT$\sim$50 and $\sim$80 eV) is detected, adding them to the growing group of "soft" IPs. White dwarf masses are determined in the range $\rm \sim0.58-1.18\,M_{\odot}$, indicating massive accreting primaries in five of them. Most sources accrete at rates lower than the expected secular value for their orbital period. Formerly proposed as a long-period (9.4 h) novalike CV, Swift J0746.3-1608 shows peculiar spectrum and light curves suggesting either an atypical low-luminosity CV or a low mass X-ray binary.

GK Per, a classical nova of 1901, is thought to undergo variable mass accretion on to a magnetized white dwarf (WD) in an intermediate polar system (IP). We organized a multi-mission observational campaign in the X-ray and ultraviolet (UV) energy ranges during its dwarf nova (DN) outburst in 2015 March-April. Comparing data from quiescence and near outburst, we have found that the maximum plasma temperature decreased from about 26 to 16.2+/-0.4 keV. This is consistent with the previously proposed scenario of increase in mass accretion rate while the inner radius of the magnetically disrupted accretion disc shrinks, thereby lowering the shock temperature. A NuSTAR observation also revealed a high-amplitude WD spin modulation of the very hard X-rays with a single-peaked profile, suggesting an obscuration of the lower accretion pole and an extended shock region on the WD surface. The X-ray spectrum of GK Per measured with the Swift X-Ray Telescope varied on time-scales of days and also showed a gradual increase of the soft X-ray flux below 2 keV, accompanied by a decrease of the hard flux above 2 keV. In the Chandra observation with the High Energy Transmission Gratings, we detected prominent emission lines, especially of Ne, Mg and Si, where the ratios of H-like to He-like transition for each element indicate a much lower temperature than the underlying continuum. We suggest that the X-ray emission in the 0.8-2 keV range originates from the magnetospheric boundary.

We present the first ever X-ray data taken of an intermediate polar, FO Aqr, when in a low accretion state and during the subsequent recovery. The Swift and Chandra X-ray data taken during the low accretion state in July 2016 both show a softer spectrum when compared to archival data taken when FO Aqr was in a high state. The X-ray spectrum in the low state showed a significant increase in the ratio of the soft X-ray flux to the hard X-ray flux due to a change in the partial covering fraction of the white dwarf from $>85\%$ to $70^{+5}_{-8}\%$ and a change in the hydrogen column density within the disc from 19$^{+1.2}_{-0.9}\times 10^{22}$ cm$^{-2}$ to 1.3$^{+0.6}_{-0.3}\times 10^{22}$ cm$^{-2}$. XMM-Newton observations of FO Aqr during the subsequent recovery suggest that the system had not yet returned to its typical high state by November 2016, with the hydrogen column density within the disc found to be 15$^{+3.0}_{-2.0}$ cm$^{-2}$. The partial covering fraction varied in the recovery state between $85\%$ and $95\%$. The spin period of the white dwarf in 2014 and 2015 has also been refined to 1254.3342(8) s. Finally, we find an apparent phase difference between the high state X-ray pulse and recovery X-ray pulse of 0.17, which may be related to a restructuring of the X-ray emitting regions within the system.

Cosmic rays in the energy range $10^{18.0}$ - $10^{18.5}$ eV are thought to have a light, probably protonic, composition. To study their origin one can search for anisotropy in their arrival directions. Extragalactic cosmic rays should be isotropic, but galactic cosmic rays of this type should be seen mostly along the galactic plane, and there should be a shortage of events coming from directions near the galactic anticenter. This is due to the fact that, under the influence of the galactic magnetic field, the transition from ballistic to diffusive behavior is well advanced, and this qualitative picture persists over the whole energy range. Guided by models of the galactic magnetic field that indicate that the enhancement along the galactic plane should have a standard deviation of about 20$^\circ$ in galactic latitude, and the deficit in the galactic anticenter direction should have a standard deviation of about 50$^\circ$ in galactic longitude, we use the data of the Telescope Array surface detector in $10^{18.0}$ to $10^{18.5}$ eV energy range to search for these effects. The data are isotropic. Neither an enhancement along the galactic plane nor a deficit in the galactic anticenter direction is found. Using these data we place an upper limit on the fraction of EeV cosmic rays of galactic origin at 1.3% at 95% confidence level.

Two XMM Newton observations of Nova V4743 Sgr (Nova Sgr 2002) were performed shortly after it returned to quiescence, 2 and 3.5 years after the explosion. The X-ray light curves revealed a modulation with a frequency of ~0.75 mHz, indicating that V4743 Sgr is most probably an intermediate polar (IP). The X-ray spectra have characteristics in common with known IPs, with a hard thermal plasma component that can be fitted only assuming a partially covering absorber. In 2004 the X-ray spectrum had also a supersoft blackbody-like component, whose temperature was close to that of the white dwarf (WD) in the supersoft X-ray phase following the outburst, but with flux by at least two orders of magnitude lower. In quiescent IPs, a soft X-ray flux component originates at times in the polar regions irradiated by an accretion column, but the supersoft component of V4743 Sgr disappeared in 2006, indicating a possible origin different from accretion. We suggest that it may have been due to an atmospheric temperature gradient on the WD surface, or to continuing localized thermonuclear burning at the bottom of the envelope, before complete turn-off. An optical spectrum obtained with SALT 11.5 years after the outburst showed a prominent He II 4686A line and the Bowen blend, which reveal a very hot region, but with peak temperature shifted to the ultraviolet (UV) range. V4743 Sgr is the third post-outburst nova and IP candidate showing a low-luminosity supersoft component in the X-ray flux a few years after the outburst.

Symbiotic star surveys have traditionally relied almost exclusively on low resolution optical spectroscopy. However, we can obtain a more reliable estimate of their total Galactic population by using all available signatures of the symbiotic phenomenon. Here we report the discovery of a hard X-ray source, 4PBC J0642.9+5528, in the Swift hard X-ray all-sky survey, and identify it with a poorly studied red giant, SU Lyn, using pointed Swift observations and ground-based optical spectroscopy. The X-ray spectrum, the optical to UV spectrum, and the rapid UV variability of SU Lyn are all consistent with our interpretation that it is a symbiotic star containing an accreting white dwarf. The symbiotic nature of SU Lyn went unnoticed until now, because it does not exhibit emission lines strong enough to be obvious in low resolution spectra. We argue that symbiotic stars without shell-burning have weak emission lines, and that the current lists of symbiotic stars are biased in favor of shell-burning systems. We conclude that the true population of symbiotic stars has been underestimated, potentially by a large factor.

We describe the X-ray emission as observed with Suzaku from five symbiotic stars that we selected for deep Suzaku observations after their initial detection with ROSAT, ASCA and Swift. We find that the X-ray spectra of all five sources can be adequately fit with absorbed, optically thin thermal plasma models, with either single- or multi-temperature plasmas. These models are compatible with the X-ray emission originating in the boundary layer between an accretion disk and a white dwarf. The high plasma temperatures of kT$~>3$ keV for all five targets were greater than expected for colliding winds. Based on these high temperatures, as well as previous measurements of UV variability and UV luminosity, and the large amplitude of X-ray flickering in 4 Dra, we conclude that all five sources are accretion-powered through predominantly optically thick boundary layers. Our X-ray data allow us to observe a small, optically thin portion of the emission from these boundary layers. Given the time between previous observations and these observations, we find that the intrinsic X-ray flux and the intervening absorbing column can vary by factors of three or more on a time scale of years. However, the location of the absorber and the relationship between changes in accretion rate and absorption are still elusive.

TARA (Telescope Array Radar) is a cosmic ray radar detection experiment colocated with Telescope Array, the conventional surface scintillation detector (SD) and fluorescence telescope detector (FD) near Delta, Utah, U.S.A. The TARA detector combines a 40 kW, 54.1 MHz VHF transmitter and high-gain transmitting antenna which broadcasts the radar carrier over the SD array and within the FD field of view, towards a 250 MS/s DAQ receiver. TARA has been collecting data since 2013 with the primary goal of observing the radar signatures of extensive air showers (EAS). Simulations indicate that echoes are expected to be short in duration (~10 microseconds) and exhibit rapidly changing frequency, with rates on the order of 1 MHz/microsecond. The EAS radar cross-section (RCS) is currently unknown although it is the subject of over 70 years of speculation. A novel signal search technique is described in which the expected radar echo of a particular air shower is used as a matched filter template and compared to waveforms obtained by triggering the radar DAQ using the Telescope Array fluorescence detector. No evidence for the scattering of radio frequency radiation by EAS is obtained to date. We report the first quantitative RCS upper limits using EAS that triggered the Telescope Array Fluorescence Detector.

In this work we are reporting on the measurement of the proton-air inelastic cross section $\sigma^{\rm inel}_{\rm p-air}$ using the Telescope Array (TA) detector. Based on the measurement of the $\sigma^{\rm inel}_{\rm p-air}$ the proton-proton cross section $\sigma_{\rm p-p}$ value is also determined at $\sqrt{s} = 95_{-8}^{+5}$ TeV. Detecting cosmic ray events at ultra high energies with Telescope Array enables us to study this fundamental parameter that we are otherwise unable to access with particle accelerators. The data used in this report is the hybrid events observed by the Middle Drum fluorescence detector together with the surface array detector collected over five years. The value of the $\sigma^{\rm inel}_{\rm p-air}$ is found to be equal to $567.0 \pm 70.5 [{\rm Stat.}] ^{+29}_{-25} [{\rm Sys.}]$ mb. The total proton-proton cross section is subsequently inferred from Glauber Formalism and Block, Halzen and Stanev QCD inspired fit and is found to be equal to $170_{-44}^{+48} [{\rm Stat.}] _{-17}^{+19} [{\rm Sys.}] $mb.

IGR J04571+4527 and Swift J0525.6+2416 are two hard X-ray sources detected in the Swift/BAT and INTEGRAL/IBIS surveys. They were proposed to be magnetic cataclysmic variables of the Intermediate Polar (IP) type, based on optical spectroscopy. IGR J04571+4527 also showed a 1218 s optical periodicity, suggestive of the rotational period of a white dwarf, further pointing towards an IP classification. We here present detailed X-ray (0.3-10 keV) timing and spectral analysis performed with XMM-Newton, complemented with hard X-ray coverage (15-70 keV) from Swift/BAT. These are the first high signal to noise observations in the soft X-ray domain for both sources, allowing us to identify the white dwarf X-ray spin period of Swift J0525.6+2416 (226.28 s), and IGR J04571+4527 (1222.6 s). A model consisting of multi-temperature optically thin emission with complex absorption adequately fits the broad-band spectrum of both sources. We estimate a white dwarf mass of about 1.1 and 1.0 solar masses for IGR J04571+4527 and Swift J0525.6+2416, respectively. The above characteristics allow us to unambiguously classify both sources as IPs, confirming the high incidence of this subclass among hard X-ray emitting Cataclysmic Variables.

We study the X-ray emission from five symbiotic stars observed with Suzaku. These objects were selected for deeper observations with Suzaku after their first detection with ROSAT and Swift. We found that the X-ray spectra can be adequately fit with absorbed optically thin thermal plasma models, either single or multi-temperature. Such a model is compatible with the X-ray emission being originated in the innermost region of the accretion disk, i.e. a boundary layer. Based on the large flickering amplitude (only detected in 4 Dra), the high plasma temperature and previous measurements of UV variability and luminosity, we conclude that all five sources are accretion-powered through predominantly opticall thick boundary layer. Given the time lapse between previous and these observations, we were able to study the long term variability of their X-ray emission and found that the intrinsic X-ray flux and the intervening absorption column can vary by factors of three or more. However, it is still elusive the location of the absorber and how the changes in the accretion rate and absorption are related.

Determining reliable distances to classical novae is a challenging but crucial step in deriving their ejected masses and explosion energetics. Here we combine radio expansion measurements from the Karl G. Jansky Very Large Array with velocities derived from optical spectra to estimate an expansion parallax for nova V959 Mon, the first nova discovered through its gamma-ray emission. We spatially resolve the nova at frequencies of 4.5-36.5 GHz in nine different imaging epochs. The first five epochs cover the expansion of the ejecta from 2012 October to 2013 January, while the final four epochs span 2014 February to 2014 May. These observations correspond to days 126 through 199 and days 615 through 703 after the first detection of the nova. The images clearly show a non-spherical ejecta geometry. Utilizing ejecta velocities derived from 3D modelling of optical spectroscopy, the radio expansion implies a distance between 0.9 +/- 0.2 and 2.2 +/- 0.4 kpc, with a most probable distance of 1.4 +/- 0.4 kpc. This distance implies a gamma-ray luminosity much less than the prototype gamma-ray-detected nova, V407 Cyg, possibly due to the lack of a red giant companion in the V959 Mon system. V959 Mon also has a much lower gamma-ray luminosity than other classical novae detected in gamma-rays to date, indicating a range of at least a factor of 10 in the gamma-ray luminosities for these explosions.

We report on the search for steady point-like sources of neutral particles around 10$^{18}$ eV between 2008 May and 2013 May with the scintillator surface detector of the Telescope Array experiment. We found overall no significant point-like excess above 0.5 EeV in the northern sky. Subsequently, we also searched for coincidence with the Fermi bright Galactic sources. No significant coincidence was found within the statistical uncertainty. Hence, we set an upper limit on the neutron flux that corresponds to an averaged flux of 0.07 km$^{-2}$ yr$^{-1}$ for $E>1$ EeV in the northern sky at the 95% confidence level. This is the most stringent flux upper limit in a northern sky survey assuming point-like sources. The upper limit at the 95% confidence level on the neutron flux from Cygnus X-3 is also set to 0.2 km$^{-2}$ yr$^{-1}$ for $E>0.5$ EeV. This is an order of magnitude lower than previous flux measurements.