The long-term Swift monitoring of ESO 243-49 HLX-1 provides an opportunity to investigate the detailed timing and spectral behaviour of this hyper-luminous X-ray source. Swift has detected 7 outbursts since 2009 mid-August. Using different dynamical timing algorithms, we confirm an increasing trend for the time intervals between outbursts, which is manifest in the delays between the latest outbursts. The X-ray spectra of HLX-1 in quiescence can be described with a single power-law model while the thermal component dominates the X-ray emission during outburst. There is only marginal evidence for photon index (or spectral hardness) changes between quiescent states with about 1{\sigma} deviation. With the updated temporal and spectral features, we re-examine different scenarios to explain the origin of the quasi-periodic modulation of HLX-1. A significantly increasing trend without obvious stochastic fluctuations on the timescale of the detected quasi-period may not fully support an orbital period origin as might be due to mass transfer episodes from a donor star at periastron of an extremely eccentric orbit. The outburst profile seems to be consistent with the effect of tidal-induced-precession of an accretion disc or an oscillating wind scenario in the inner disc. Based on these models, we speculate that the true orbital period is much shorter than the detected quasi-periodicity.

We report a systematic analysis of the spin, orbital, and superorbital modulations of 4U 0114+650, a high-mass X-ray binary consists of one of the slowest spinning neutron stars. Utilizing the dynamic power spectrum, we found that the spin period varied dramatically and anti-correlated with the long-term X-ray flux variation that can be observed using RXTE ASM, Swift BAT, and MAXI. The spin-up rate over the entire dataset is consistent with previously reported values; however, the local spin-up rate is considerably higher. The corresponding local spin-up timescale is comparable to the local spin-up rate of OAO 1657-415, indicating that 4U 0114+650 could also have a transient disk. Moreover, the spin period evolution shows two ~1000-day spin-down/random-walk epochs that appeared together with depressions of the superorbital modulation amplitude. This implies that the superorbital modulation was closely related to the presence of the accretion disk, which is not favored in the spin-down/random-walk epochs because the accretion is dominated by the direct wind accretion.The orbital period is stable during the entire time span; however, the orbital profile changes with time significantly. We found that the depth of the dip near the inferior conjunction of the companion is highly variable, which disfavors the eclipsing scenario. Moreover, the dip was less obvious during the spin-down/random-walk epochs, indicating its correlation with the accretion disk. Further monitoring in both X-ray and optical bands could reveal the establishment of the accretion disk in this system.

A glitch of a pulsar is known as a sudden increase in the spin frequency and spin-down rate (frequency time derivative), and it can be caused by a sudden rel\ ease of the stress built up in the solid crust of the star or pinned vortices in the superfluid interior. PSR J2021+4026 is the first pulsar that shows a significant change in the gamma-ray flux and pulse profile at the glitch that occurred around 2011 October 16. We report the results of timing and spectral analysis of PSR~J2021+4026 using $\sim$ 8~yr Fermi-LAT data. We find that the pulsar stayed at a high spin-down rate ($\sim 4\%$ higher than the pre-glitch value) and a low gamma-ray state ($\sim 18\%$ lower) for about 3~yr after the glitch. Around 2014 December, the spin-down rate and gamma-ray flux gradually returned to pre-glitch values within a time scale of a few months. The phase-resolved spectra and pulse profiles after the relaxation are also consistent with those before the glitch. The observed long-term evolution of the spin-down rate and the gamma-ray flux indicates that the glitch triggered a mode change in the global magnetosphere. We speculate that the glitch changed the local magnetic field structure around the polar cap and/or the inclination angle of the dipole axis, leading to a change in the electric current circulating in the magnetosphere.

We report a possible detection of a 55-day X-ray modulation for the ultraluminous accreting pulsar M82 X-2 from archival Chandra observations. Because M82 X-2 is known to have a 2.5-day orbital period, if the 55-day period is real, it will be the superorbital period of the system. We also investigated variabilities of other three nearby ultraluminous X-ray sources in the central region of M82 with the Chandra data and did not find any evidence of periodicities. Furthermore, we re-examined the previously reported 62-day periodicity near the central region of M82 by performing a systematic timing study with all the archival Rossi X-Ray Timing Explorer and Swift data. Using various dynamic timing analysis methods, we confirmed that the 62-day period is not stable, suggesting that it is not the orbital period of M82 X-1 in agreement with previous work.

We extend the work of Yen et al. (2012) and develop 2nd order formulae to accommodate a nested grid discretization for the direct self-gravitational force calculation for infinitesimally thin gaseous disks. This approach uses a two-dimensional kernel derived for infinitesimally thin disks and is free of artificial boundary conditions. The self-gravitational force calculation is presented in generalized convolution forms for a nested grid configuration. A numerical technique derived from a fast Fourier transform is employed to reduce the computational complexity to be nearly linear. By comparing with analytic potential-density pairs associated with the generalized Maclaurin disks, the extended approach is verified to be of second order accuracy using numerical simulations. The proposed method is accurate, computationally fast and has the potential to be applied to the studies of planetary migration and the gaseous morphology of disk galaxies.

Investigating the evolution of disk galaxies and the dynamics of proto-stellar disks can involve the use of both a hydrodynamical and a Poisson solver. These systems are usually approximated as infinitesimally thin disks using two- dimensional Cartesian or polar coordinates. In Cartesian coordinates, the calcu- lations of the hydrodynamics and self-gravitational forces are relatively straight- forward for attaining second order accuracy. However, in polar coordinates, a second order calculation of self-gravitational forces is required for matching the second order accuracy of hydrodynamical schemes. We present a direct algorithm for calculating self-gravitational forces with second order accuracy without artifi- cial boundary conditions. The Poisson integral in polar coordinates is expressed in a convolution form and the corresponding numerical complexity is nearly lin- ear using a fast Fourier transform. Examples with analytic solutions are used to verify that the truncated error of this algorithm is of second order. The kernel integral around the singularity is applied to modify the particle method. The use of a softening length is avoided and the accuracy of the particle method is significantly improved.

Long-term X-ray modulations on timescales from tens to hundreds of days have been widely studied for X-ray binaries located in the Milky Way and the Magellanic Clouds. For other nearby galaxies, only the most luminous X-ray sources can be monitored with dedicated observations. We here present the first systematic study of long-term X-ray variability of four ultraluminous X-ray sources (ESO 243-49 HLX-1, Holmberg IX X-1, M81 X-6, and NGC 5408 X-1) monitored with Swift. By using various dynamic techniques to analyse their light curves, we find several interesting low-frequency quasi-periodicities. Although the periodic signals may not represent any stable orbital modulations, these detections reveal that such long-term regular patterns may be related to superorbital periods and structure of the accretion discs. In particular, we show that the outburst recurrence time of ESO 243-49 HLX-1 varies over time and suggest that it may not be the orbital period. Instead, it may be due to some kinds of precession, and the true binary period is expected to be much shorter.

We report the X-ray pulsation of ~173.3 ms for the "next Geminga", PSR J1836+5925, with recent XMM-Newton investigations. The X-ray periodicity is consistent wtih the gamma-ray ephemeris at the same epoch. The X-ray folded light curve has a sinusoidal structure which is different from the double-peaked gamma-ray pulse profile. We have also analysed the X-ray phase-averaged spectra which shows the X-ray emission from PSR J1836+5925 is thermal dominant. This suggests the X-ray pulsation mainly originates from the modulated hot spot on the stellar surface.

NGC 1097 is a nearby barred spiral galaxy believed to be interacting with the elliptical galaxy NGC 1097A located to its northwest. It hosts a Seyfert 1 nucleus surrounded by a circumnuclear starburst ring. Two straight dust lanes connected to the ring extend almost continuously out to the bar. The other ends of the dust lanes attach to two main spiral arms. To provide a physical understanding of its structural and kinematical properties, two-dimensional hydrodynamical simulations have been carried out. Numerical calculations reveal that many features of the gas morphology and kinematics can be reproduced provided that the gas flow is governed by a gravitational potential associated with a slowly rotating strong bar. By including the self-gravity of the gas disk in our calculation, we have found the starburst ring to be gravitationally unstable which is consistent with the observation in \citet{hsieh11}. Our simulations show that the gas inflow rate is 0.17 M$_\sun$ yr$^{-1}$ into the region within the starburst ring even after its formation, leading to the coexistence of both a nuclear ring and a circumnuclear disk.

NGC 4945 is a Seyfert 2 galaxy at a distance of 3.82 Mpc. Its relative proximity has permitted a detailed SMA study of the circumnuclear molecular gas in a galaxy exhibiting an AGN. Based on an analysis of the high-resolution velocity field of the central region (20" X 20", 1" = 19 pc), we demonstrate that the S-shaped structure of the isovelocity contours is well reproduced by the numerical results of a two dimensional hydrodynamical simulation. In particular, the velocity structure is represented by the bending produced by a shock along the spiral density waves, which are excited at the outer-inner Lindblad resonance by a fast rotating bar. The simulated density map reveals a pair of tightly wound spirals in the center which pass through most of the ring-like (claimed to be a circumnuclear starburst ring by other authors) high intensity region in the observations as well as intersect several Pa$\alpha$ emission line knots located outside the ring-like region. The calculated mass inflow rate at a scale of 50 pc is about three times the inferred mass accretion rate of the AGN of NGC 4945. We find that self-gravity of the gas is important and should be included in our model for NGC 4945. The model is compared with the gas orbit model discussed in Lim et al. (2009), and it is shown that the hydrodynamic model provides a better match to the observed position-velocity diagram and, hence, provides a more reliable prediction of the outer inner Lindblad resonance position.

We performed a series of hydro-dynamic simulations to investigate the orbital migration of a Jovian planet embedded in a proto-stellar disk. In order to take into account of the effect of the disk's self gravity, we developed and adopted an \textbf{Antares} code which is based on a 2-D Godunov scheme to obtain the exact Reimann solution for isothermal or polytropic gas, with non-reflecting boundary conditions. Our simulations indicate that in the study of the runaway (type III) migration, it is important to carry out a fully self consistent treatment of the gravitational interaction between the disk and the embedded planet. Through a series of convergence tests, we show that adequate numerical resolution, especially within the planet's Roche lobe, critically determines the outcome of the simulations. We consider a variety of initial conditions and show that isolated, non eccentric protoplanet planets do not undergo type III migration. We attribute the difference between our and previous simulations to the contribution of a self consistent representation of the disk's self gravity. Nevertheless, type III migration cannot be completely suppressed and its onset requires finite amplitude perturbations such as that induced by planet-planet interaction. We determine the radial extent of type III migration as a function of the disk's self gravity.